Preprint
Article

Equilibrium Figures for a Rotating Compressible Capillary Two-Layer Liquid

Altmetrics

Downloads

56

Views

30

Comments

0

A peer-reviewed article of this preprint also exists.

Submitted:

09 November 2023

Posted:

10 November 2023

You are already at the latest version

Alerts
Abstract
The paper proves the existence of a family of axisymmetric equilibrium figures as solutions of a stationary problem with unknown boundaries for the Navier–Stokes equations corresponding to the slow rotation of a viscous compressible two-layer liquid mass about some axis. It is assumed that the liquids are barotropic, capillary and have different viscosities, the internal fluid being bounded by a closed surface. This interface does not intersect with the external boundary of the cloud. The proof is based on implicit function theorem and carried out in the Hölder spaces.
Keywords: 
Subject: Computer Science and Mathematics  -   Mathematics

MSC:  Primary 35Q30; 76T06; 76D05; 76D06; Secondary 35R35; 76D03

1. Introduction

Existence of an equilibrium surface for an isolated compressible liquid mass rotating about a fixed axis was first proved in [9]. Our aim is to prove the existence of equilibrium figures for a rotating compressible two-layer fluid.
The problem of the rotation of an isolated incompressible liquid mass about a fixed axis as a rigid body was considered by many famous mathematicians, among them were Newton, Maclaurin, Jacobi, Kovalevskaya, Lyapunov, Poincare and others [1,2,3], who mainly studied the movement without capillarity. The capillary fluids were first investigated by Globa-Mikhailenko [4], Boussinesq and Charrueau in the begining of 20th century. The latter gave a detailed analysis of the problem, calculated the shape of equilibrium figures, including the toroidal case, and considered some aspects of the stability [5,6]. These results were included in a big review on this subject presented in the book of Appell [7]. Stability problem for various ellipsoidal equilibrium figures is analyzed in monograph [8].
Now we state, in a complete setting, the problem on unsteady motion of two compressible barotropic fluids of finite volume separated by a closed unknown interface.
At the initial instant t = 0 , let a fluid with dynamic viscosities μ + , μ 1 + be in a bounded domain Ω 0 + R 3 , and in the domain Ω 0 , surrounding it, there be a fluid with dynamic viscosities μ , μ 1 ;
μ ± > 0 , 2 μ ± + 3 μ 1 ± 0 .
The domain Ω 0 Ω 0 + ¯ Ω 0 is bounded by the free surface Γ 0 and includes the closed interface Γ 0 + Ω + ; Γ 0 ± are given. This two-component cloud rotates about the vertical axis x 3 with an angular velocity ω .
For t > 0 , it is necessary to find the surfaces Γ t , Γ t + , as well as velocity vector field v ( x , t ) = ( v 1 , v 2 , v 3 ) and the density ρ ( x , t ) > 0 of the fluids satisfying diffraction problem for the Navier – Stokes system
Preprints 90135 i001
where
T = p ( ρ ) + μ 1 ± · v I + μ ± S ( v )
is stress tensor, S ( v ) i j = v i x j + v j x i is double strain rate tensor, I is identity matrix; μ ± , μ 1 ± are step functions of dynamic viscosities, equal to μ + , μ 1 + in Ω t + and μ , μ 1 in Ω t respectively; p ( ρ ) is fluid pressure given by a known smooth density function; v 0 and ρ 0 > 0 are initial distributions of velocity and density of the liquids, n is the outward normal vector to the union Γ t ; H ± ( x , t ) are twice the mean curvatures of the surfaces Γ t ± (moreover, H + < 0 at points of convexity Γ t + towards Ω t ); σ , σ + > 0 are surface tension coefficients on Γ t and Γ t + , respectively; V n is the rate of evolution of Γ t in the direction n . We assume that the Cartesian coordinate system { x } is introduced in the space R 3 . The central dot denotes the Cartesian scalar product.
We mean summation over repeated indices from 1 to 3 if they are denoted by Latin letters, and from 1 to 2 if they are Greek. We mark vectors and vector spaces in bold. The notation · T denotes the vector with the components ( · T ) j = T i j x i , j = 1 , 2 , 3 .
The kinematic boundary condition v · n = V n excludes mass transfer across fluid boundaries. It follows from our assumption that the fluid particles do not leave the boundaries Γ t ± during the time.
Local (in time) solvability was proved for problem (1) in the whole space R 3 with a closed interface between the fluids. The result was obtain both in the Sobolev – Slobodetskiǐ classes of functions [10] and in the Hölder ones [11]. One can get similar results for a two-component domain bounded by a free boundary if one takes into account the estimates for a model problem in a half-space [12,13].
As we have mentioned, we suppose the liquids to be barotropic which implies that the pressure p is a known increasing function of the density: p ( ρ ) > 0 . Let, in addition, ρ = ρ ( | x | ) , x = ( x 1 , x 2 , 0 ) .
We assume that equilibrium figures F + , F are nearly globular domains with the radiuses R 0 ± ( R 0 + < R 0 ), and the motion of fluids is close to the state of rest, i.e., the velocity v is small, and the density ρ differs little from a step function ρ ± > 0 . This picture is schematically presented in Figure 1. We denote the balls { x R 3 : | x | R 0 ± } by B R 0 ± .
We are going to prove the existence of G + and G , the boundaries of the figures F + and F , respectively. We follow the plan of paper [9].
Figure 1. Equilibrium Figures for a Two-Layer Compressible Fluid.
Figure 1. Equilibrium Figures for a Two-Layer Compressible Fluid.
Preprints 90135 g001
At rest, the bubble consisted of nested spherical two layers B R 0 + and B R 0 B R 0 + with uniform distributions of densities ρ ± has the piecewise constant pressure:
p ( ρ ) = 2 σ R 0 in B R 0 B R 0 + , p ( ρ + ) = 2 σ + R 0 + + 2 σ R 0 in B R 0 + .
The masses of the layers are
ρ + | B R 0 + | = m + , ρ | B R 0 | | B R 0 + | = m .
Steady motion of a two-layer gaseous body F F + ¯ F uniformly rotating about the axis x 3 with a constant angular velocity ω is governed by the homogeneous stationary Navier–Stokes equations
ρ ( V · ) V · T = 0 , · ( ρ V ) = 0 in F ±
(here the density ρ and velocity V depend only on x) and the boundary conditions
T ( V , P ) n | G = σ H n on G , [ V ] | G + = 0 , [ T ( V , P ) n ] | G + = σ + H + n on G + , V · n = 0 on G = G + G ,
where H , H + are twice the mean curvatures of G , G + , respectively. The last relation follows from the boundary condition v · n = V n . The pressure P depends on ρ .
It is easily seen that velocity vector field
V ( x ) = ω e 3 × x ω ( x 2 , x 1 , 0 )
satisfies (4) together with pressure function gradient
P ( ρ ) = ρ ω 2 x ρ ω 2 2 | x | 2 ,
where e i is the ith basis vector, | x | 2 = x 1 2 + x 2 2 .
First, we consider the simple case when equality (7) coincides with the following one
P ( ρ ) = P ( ρ ) ω 2 2 | x | 2 ,
whence P ( ρ ) = ρ ( x ) and P ± ( ρ ) ρ 2 ( x ) 2 + p ± in F ± with constants p ± , because pressure functions can differ each other in different domains by a constant. These constants can be found from relations (2):
p = 2 σ R 0 ρ 2 2 , p + = 2 σ + R 0 + + 2 σ R 0 ρ + 2 2 .
Let S 1 be the unit sphere in R 3 with the center in zero, ξ = x | x | S 1 . We suppose G ± to be given by functions R ± ( ξ ) on S 1 . In addition, let R ± ( ξ ) be rotationally symmetric, i.e., they depend only on | ξ | = ξ 1 2 + ξ 2 2 and ξ 3 , and be even in ξ 3 .
By substituting V given by (6) and P = P ± into boundary conditions (5), we obtain the equations for the surface G of the domain F and for the interface G + between the fluids:
σ H ( x ) + P ( ρ ) = 0 , x G , σ + H + ( x ) + [ P ( ρ ) ] | G + = 0 , x G + .
Rotationally symmetry implies that R ± ( ξ ) do not depend on arctan ( ξ 2 ξ 1 ) . It is clear that n = ( ξ 1 c 0 , ξ 2 c 0 , n 3 ) , since the first two components of n are proportional to ones of the radial vector of the circles to be horizontal sections of G ± . Therefore V · n = 0 on G ± .
Obviously, the density is given by the formulas
ρ ( x ) = ω 2 2 | x | 2 + c + in F + , ρ ( x ) = ω 2 2 | x | 2 + c in F
with arbitrary positive constants c + and c , equations (9) taking the form
σ H ( x ) + P ω 2 2 | x | 2 + c = 0 , x G , σ + H + ( x ) + P + ω 2 2 | x | 2 + c + P ω 2 2 | x | 2 + c = 0 , x G + .
One can determine the constants c ± by prescribing the masses of fluids to be the same as that of the nested spherical liquid layers (3):
F ± ρ ( x ) d x F ± ω 2 2 | x | 2 + c ± d x = m ± .
We consider the angular momentum to be one more parameter of the problem. It is given:
β F + F ρ ( x ) V · η 3 d x = ω F ± ω 2 2 | x | 2 + c ± | x | 2 d x ,
where η i = e i × x is the ith solid rotation vector, i = 1 , 2 , 3 . Then the angular velocity ω is a function of β .
We denote by C s , s > 0 , s N , the Hölder space of functions f on the sphere S 1 with the norm
| f | C s ( S 1 ) max k = { 1 , . . . , N } | j | < s sup ξ ζ k | D j f ( ξ ) | + | j | = [ s ] sup ξ , ξ ¯ ζ k | ξ ξ ¯ | ( s [ s ] ) | D j f ( ξ ) D j f ( ξ ¯ ) | ,
where D j f is | j | th derivative of f calculated in local coordinates on the subdomain ζ k S 1 , k = 1 N ζ k = S 1 . Under C ˜ s ( S 1 ) , we mean the subspace of C s ( S 1 ) consisting of rotationally symmetric functions that are even with respect to ξ 3 .
Theorem 1.
Let α > 0 , α N , and let the data of problem (4), (5) be such that condition (26) holds. Then for an arbitrary β satisfying the estimate
| β | < ε
with small enough ε, there exists a unique solution ( R ± , ω , c ± ) C ˜ 2 + α ( S 1 ) × C ˜ 2 + α ( S 1 ) × R × R × R to system (11)–(13). It obeys the inequality
± | R ± R 0 ± | C 2 + α ( S 1 ) + | c ± ρ ± | + | ω | < c | β | .

2. Proof of Theorem 1

Proof. 
In order to linearize system (11), we apply the formula for the first variation of a functional δ 0 R [ r ] = d d s R [ s r ] | s = 0 . According to [14,15], the first variation of twice the mean curvature of G ± with respect to the double curvature of the sphere S R 0 ± is
δ 0 H ± [ r ( ξ ) ] + 2 R 0 ± = 1 R 0 ± 2 Δ S 1 r ± + 2 r ± ,
where Δ S 1 is the Laplace – Beltrami operator on S 1 and r ± R ± R 0 ± ; G ± = { x = R 0 ± ξ + r ± ( ξ ) N , ξ S 1 } , N is the outward normal to S 1 . This yields
σ R 0 2 Δ S 1 r + 2 r = k + f , ξ S 1 , σ + R 0 + 2 Δ S 1 r + + 2 r + = k + + f + , ξ S 1 ,
where k = 2 σ R 0 P ( c ) P ( ρ ) P ( c ) , k + = 2 σ + R 0 + P + ( c + ) + P ( c ) P + ( ρ + ) P + ( c + ) + P ( c ) P ( ρ ) in view of (2);
f = σ δ 0 H ( x ) + 2 R 0 H ( x ) + 2 R 0 P ω 2 2 | x | 2 + c + P ( c ) , f + = σ + δ 0 H + ( x ) + 2 R 0 + H + ( x ) + 2 R 0 + P + ω 2 2 | x | 2 + c + + P ω 2 2 | x | 2 + c + P + ( c + ) P ( c ) .
We integrate equations (16) over S 1 by parts, then we get
2 σ ± R 0 ± 2 S 1 r ± d ξ = S 1 ( k ± + f ± ) d ξ .
The integrals S 1 r ± d ξ can be expressed in terms of the differences of the volumes | F | | B R 0 | and | F + | | B R 0 + | , | B R 0 ± | = 4 π 3 R 0 ± 3 . For example,
| F | | B R 0 | = 1 3 S 1 ( R 3 R 0 3 ) d ξ = R 0 2 S 1 r d ξ + R 0 S 1 r 2 d ξ + 1 3 S 1 r 3 d ξ ,
and hence
S 1 r + d ξ = 1 R 0 + 2 | F + | | B R 0 + | + Q [ r + ] , S 1 r d ξ = 1 R 0 2 | F | | B R 0 | + Q [ r ] ,
where Q [ r ± ] 1 R 0 ± S 1 r ± 2 d ξ 1 3 R 0 ± 2 S 1 r ± 3 d ξ .
We rewrite (12) as follows
m + = S 1 d ξ 0 R + ( ξ ) ω 2 2 | ξ | 2 s 4 d s + c + | F + | = ω 2 10 S 1 R + 5 ( ξ ) | ξ | 2 d ξ + c + | F + | , m = S 1 d ξ R + ( ξ ) R ( ξ ) ω 2 2 | ξ | 2 s 4 d s + c | F | = ω 2 10 S 1 ( R 5 R + 5 ) | ξ | 2 d ξ + c | F | ,
where | ξ | 2 = ξ 1 2 + ξ 2 2 . On the other hand, m + = ρ + | B R 0 + | , m = ρ ( | B R 0 | | B R 0 + | ) . That’s why
0 = ω 2 10 S 1 R + 5 ( ξ ) | ξ | 2 d ξ + ρ + | F + | | B R 0 + | + | F + | c + ρ + , 0 = ω 2 10 S 1 ( R 5 R + 5 ) | ξ | 2 d ξ + ρ | F | | B R 0 | + | B R 0 + | + | F | c ρ .
We express | F ± | | B R 0 ± | from (20) and substitute in (19). Then equalities (18) imply that
2 σ + R 0 + 2 1 ρ + R 0 + 2 ω 2 10 S 1 R + 5 ( ξ ) | ξ | 2 d ξ + | F + | c + ρ + + Q [ r + ] = S 1 ( k + + f + ) d ξ
and
2 σ R 0 2 1 R 0 2 | F + | + | F | | B R 0 | ± | B R 0 + | + Q [ r ] = S 1 ( k + f ) d ξ .
Finally, we have
2 σ R 0 2 { 1 ρ R 0 2 { ω 2 10 S 1 ( R 5 R + 5 ) | ξ | 2 d ξ + | F | c ρ + ρ ρ + ω 2 10 S 1 R + 5 ( ξ ) | ξ | 2 d ξ + | F + | c + ρ + } + Q [ r ] } = S 1 ( k + f ) d ξ .
In order to prove the solvability of system (16), (13), (21), (22), we use implicit function theorem. We represent this system as a nonlinear vector equation:
Φ ( φ ) = 0 ,
where Φ ( φ ) = ( Φ 1 ± , Φ 2 , Φ 3 ± ) , φ ( r ± , ω , λ ± ) , λ ± c ± ρ ± ,
Φ 1 ± ( φ ) = σ ± R 0 ± 2 Δ S 1 r ± + 2 r ± k ± f ± , Φ 2 ( φ ) = ω S 1 | ξ | 2 d ξ 0 R 0 + + r + ω 2 2 | ξ | 2 s 2 + c + s 4 d s + ω S 1 | ξ | 2 d ξ R 0 + + r + R 0 + r ω 2 2 | ξ | 2 s 2 + c s 4 d s β , Φ 3 + ( φ ) = 2 σ + R 0 + 2 1 ρ + R 0 + 2 ω 2 10 S 1 R + 5 ( ξ ) | ξ | 2 d ξ + | F + | λ + Q [ r + ] + S 1 ( k + + f + ) d ξ , Φ 3 ( φ ) = 2 σ R 0 2 { 1 ρ R 0 2 { S 1 d ξ R 0 + + r + R 0 + r ω 2 2 | ξ | 2 s 4 d s + | F | λ + ρ ρ + ω 2 10 S 1 R + 5 ( ξ ) | ξ | 2 d ξ + | F + | λ + } Q [ r ] } + S 1 ( k + f ) d ξ .
We linearize (23) at zero and show that the derivative of Φ at this point is not equal to zero. Thus,
Φ ( 0 ) φ + ψ ( φ ) = 0 ,
where Φ ( 0 ) = ( Φ 1 ± , Φ 2 , Φ 3 ± ) ( 0 ) is the Fréchet derivative of Φ , that is,
Φ 1 + ( 0 ) φ = σ + R 0 + 2 Δ S 1 r + + 2 r + + ρ + λ + ρ λ , Φ 1 ( 0 ) φ = σ R 0 2 Δ S 1 r + 2 r + ρ λ , Φ 2 ( 0 ) φ = ω ρ + S 1 | ξ | 2 d ξ 0 R 0 + s 4 d s + ω ρ S 1 | ξ | 2 d ξ R 0 + R 0 s 4 d s = 8 π 15 ω ρ + R 0 + 5 + ρ R 0 5 R 0 + 5 , Φ 3 + ( 0 ) φ = 2 σ + R 0 + 2 | B R 0 + | λ + ρ + R 0 + 2 4 π ρ + λ + + 4 π ρ λ = 4 π P + ( ρ + ) P ( ρ ) 3 ρ + ρ + λ + + ρ λ , Φ 3 ( 0 ) φ = 2 σ R 0 2 1 ρ R 0 2 | B R 0 | | B R 0 + | λ + ρ ρ + | B R 0 + | λ + 4 π ρ λ = 4 π P ( ρ ) R 0 + 3 3 ρ + R 0 3 λ + + R 0 3 R 0 + 3 3 ρ R 0 3 P ( ρ ) ρ λ ,
and ψ ( φ ) = Φ ( φ ) Φ ( 0 ) φ ( ψ 1 ± , ψ 2 , ψ 3 ± ) with
ψ 1 + ( φ ) = P + ( c + ) P + ( ρ + ) P ( c ) + P ( ρ ) ρ + λ + + ρ λ f + , ψ 1 ( φ ) = P ( c ) P ( ρ ) ρ λ f , ψ 2 ( φ ) = ω S 1 | ξ | 2 d ξ 0 R 0 + ω 2 2 | ξ | 2 s 2 + λ + s 4 d s + ω S 1 | ξ | 2 d ξ R 0 + R 0 + + r + ω 2 2 | ξ | 2 s 2 + c + s 4 d s + ω S 1 | ξ | 2 d ξ R 0 + R 0 ω 2 2 | ξ | 2 s 2 + λ s 4 d s + ω S 1 | ξ | 2 d ξ R 0 R 0 + r ω 2 2 | ξ | 2 s 2 + c s 4 d s β , ψ 3 + ( φ ) = 2 σ + R 0 + 2 1 ρ + R 0 + 2 ω 2 10 S 1 R + 5 ( ξ ) | ξ | 2 d ξ + ( | F + | | B R 0 + | ) λ + Q [ r + ] + S 1 ( k + ρ + λ + + ρ λ + f + ) d ξ , ψ 3 ( φ ) = 2 σ R 0 2 { 1 ρ R 0 2 { ω 2 10 S 1 ( R 5 R + 5 ) | ξ | 2 d ξ + ρ ρ + ( ω 2 10 S 1 R + 5 ( ξ ) | ξ | 2 d ξ + ( | F + | | B R 0 + | ) λ + ) + ( | F | | B R 0 | + | B R 0 + | ) λ } Q [ r ] } + S 1 ( k ρ λ + f ) d ξ .
Let X = C ˜ 2 + α × C ˜ 2 + α × R × R × R and Y = C ˜ α × C ˜ α × R × R × R . Using fixed point theorem, we will show that equation (23) is solvable in X .
In [16] (see the corollary to Theorem 2.1), it was proved that for any f C ˜ α ( S 1 ) the equation
Δ S 1 r + 2 r = f
has a unique solution r C ˜ 2 + α ( S 1 ) and
| r | C 2 + α ( S 1 ) c 1 | f | C α ( S 1 ) .
Moreover, due to the assumptions of the theorem, the determinant of the matrix in the left hand sides of two equations Φ 3 ± ( 0 ) φ = ψ 3 ± ( φ ) is not equal to zero:
2 σ + 3 ρ + R 0 + ρ + 2 σ ( R 0 3 R 0 + 3 ) 3 ρ R 0 4 ρ 2 σ ρ R 0 + 3 3 ρ + R 0 4 0 .
(Here we have taken into account relations (8).) Hence, this matrix is invertible. Thus, the whole operator Φ ( 0 ) : X Y is also invertible and
φ X c 2 Φ ( 0 ) φ Y ,
where φ X = | r + | C 2 + α ( S 1 ) + | r | C 2 + α ( S 1 ) + | ω | + | λ + | + | λ | and
Φ ( 0 ) φ Y = | Φ 1 + ( 0 ) φ | C α ( S 1 ) + | Φ 1 ( 0 ) φ | C α ( S 1 ) + | Φ 2 ( 0 ) φ | + | Φ 3 + ( 0 ) φ | + | Φ 3 ( 0 ) φ | .
Therefore, equation (24) can be written in the form:
φ + A ψ ( φ ) = 0 , A = Φ ( 0 ) 1 .
Let us now estimate the non-linear operator ψ . The term f given by (17) can be written as
f = σ 0 1 ( 1 s ) d 2 d s 2 H [ s r ] d s P ( c ¯ ) ω 2 2 R + 2 ( ξ ) | ξ | 2 ,
where H [ s r ] is twice the mean curvature of the surface G s ; G s ± = { x = R 0 ± ξ + s r ± ( ξ ) N , ξ S 1 } , N ξ | ξ | ; c ¯ is some point of the interval [ c , ω 2 2 R + 2 ( ξ ) | ξ | 2 + c ] . One can write a similar formula for f + . The terms f ± have the second order of smallness with respect to φ . After simple calculations for two values of the arguments φ = ( r ± , ω , λ ± ) and φ = ( r ± , ω , λ ± ) X such that
φ X , φ X η , η min { 1 , R 0 + / 2 , ( R 0 R 0 + ) / 2 } ,
it can be shown that the estimates
ψ ( φ ) Y c φ X 1 + α + | β | c 3 η φ X α + | β | , ψ ( φ ) ψ ( φ ) Y c 4 η α φ φ X
hold. This implies that A ψ ( φ ) maps the ball { φ X : φ X η } into itself and it is a contraction operator if
c 2 ( c 3 η 1 + α + | β | ) η , c 2 c 4 η α < 1 .
These inequalities are satisfied for η such that
2 c 2 | β | η min { ( 2 c 2 c 3 ) 1 / α , ( c 2 c 4 ) 1 / α } .
So if
| β | 1 2 c 2 min { ( 2 c 2 c 3 ) 1 / α , ( c 2 c 4 ) 1 / α } ,
then, by fixed point theorem, equation (27) has a unique solution in the ball { φ X η } with η satisfying (29). Thus, equation (23) is also uniquely solvable for β obeying inequality (30), which gives the estimate of ε in condition (14).
We now prove estimate (15). Since
φ X A ψ ( φ ) Y c 2 c 3 η α φ X + | β | ,
then for c 2 c 3 η α 1 2
φ X 2 c 2 | β | ,
which coincides with inequality (15). □
Thus, we have shown that, for certain data of the problem, there are equilibrium figures for a two-layer compressible fluid, pressure function being P ± = ρ 2 ( x ) 2 + p ± and the density ρ being determined by formulas (10).

3. General case

Now, we study the general case of pressure function. We assume only that P ( ρ ) is positive and may have a jump across G + .
Velocity vector field remains the same as above:
V ( x ) = ω ( x 2 , x 1 , 0 ) .
It satisfies (4) together with pressure function gradient P ( ρ ) = ρ ω 2 2 | x | 2 . We introduce Q ( ρ ) such that
Q ( ρ ) P ( ρ ) ρ ρ = ω 2 2 | x | 2 .
The function Q ( ρ ) = ρ 1 ρ P ( s ) s d s , ρ 1 0 . Since Q ( ρ ) = P ( ρ ) ρ > 0 , there is an inverse function Q 1 . And (31) implies
ρ = Q 1 ω 2 2 | x | 2 + C ± in F ±
with arbitrary constants C ± .
Substituting ρ into equations (9), we have
σ H ( x ) + P Q 1 ω 2 2 | x | 2 + C = 0 , x G , σ + H + ( x ) + P + Q 1 ω 2 2 | x | 2 + C + P Q 1 ω 2 2 | x | 2 + C = 0 , x G + .
We prescribe the masses of fluids and apply (32):
m ± = F ± ρ ( x ) d x F ± Q 1 ω 2 2 | x | 2 + C ± d x .
Then one has the equations for determining the constants C ± .
Similar, a given angular momentum β defines the angular velocity ω :
β F + F ρ ( x ) V · η 3 d x = ω F ± Q 1 ω 2 2 | x | 2 + C ± | x | 2 d x .
Let us state the main theorem.
Theorem 2.
Let α > 0 , α N , and let P ± ( ρ ) C α ( R + ) be positive increasing functions such that equalities (2) are satisfied for it. Here R + { x R | x > 0 } . We assume also that the data of problem (4), (5) are subjected to condition (47). Then for an arbitrary β satisfying the estimate
| β | < ε
with ε small enough, there exists a unique solution ( R ± , ω , C ± ) C ˜ 2 + α ( S 1 ) × C ˜ 2 + α ( S 1 ) × R × R × R to system (33)–(35), and the inequality
± | R ± R 0 ± | C 2 + α ( S 1 ) + | Q 1 ( C ± ) ρ ± | + | ω | < c | β |
holds.
Proof. 
After linearisation, system (33) takes the form of (16) with
k = 2 σ R 0 P ( Q 1 ( C ) ) P ( ρ ) P ( Q 1 ( C ) ) , k + = 2 σ + R 0 + P + ( Q 1 ( C + ) ) + P ( Q 1 ( C ) ) P + ( ρ + ) P + ( Q 1 ( C + ) ) + P ( Q 1 ( C ) ) P ( ρ )
due to (2) and
f = σ δ 0 H ( x ) + 2 R 0 H ( x ) + 2 R 0 P Q 1 ω 2 2 | x | 2 + C + P Q 1 ( C ) , f + = σ + δ 0 H + ( x ) + 2 R 0 + H + ( x ) + 2 R 0 + P + Q 1 ω 2 2 | x | 2 + C + + P Q 1 ω 2 2 | x | 2 + C + P + Q 1 ( C + ) P Q 1 ( C ) .
Following Sect. 2, we obtain (18) by integrating equations (16) by parts but now with new k ± (38) and f ± (39). We substitute (19) in equality (18):
2 σ + R 0 + 2 1 R 0 + 2 | F + | | B R 0 + | + Q ˜ [ r + ] = S 1 ( k + + f + ) d ξ , 2 σ R 0 2 1 R 0 2 | F + | + | F | | B R 0 | ± | B R 0 + | + Q ˜ [ r ] = S 1 ( k + f ) d ξ ,
where Q ˜ [ r ± ] 1 R 0 ± S 1 r ± 2 d ξ 1 3 R 0 ± 2 S 1 r ± 3 d ξ .
We can write (34) as follows
m + = S 1 d ξ 0 R + ( ξ ) Q 1 ω 2 2 | ξ | 2 s 2 + C + Q 1 ( C + ) s 2 d s + Q 1 ( C + ) | F + | , m = S 1 d ξ R + ( ξ ) R ( ξ ) Q 1 ω 2 2 | ξ | 2 s 2 + C Q 1 ( C ) s 2 d s + Q 1 ( C ) | F | .
In view of (3), we have
0 = S 1 d ξ 0 R + ( ξ ) Q 1 ω 2 2 | ξ | 2 s 2 + C + Q 1 ( C + ) s 2 d s + ρ + | F + | | B R 0 + | + | F + | Q 1 ( C + ) ρ + , 0 = S 1 d ξ R + ( ξ ) R ( ξ ) Q 1 ω 2 2 | ξ | 2 s 2 + C Q 1 ( C ) s 2 d s + ρ | F | | B R 0 | + | B R 0 + | + | F | Q 1 ( C ) ρ .
By expressing the differences | F + | | B R 0 + | and | F | | B R 0 | + | B R 0 + | from equalities (41) and substituting them in (40), we arrive at
2 σ + R 0 + 2 { 1 ρ + R 0 + 2 ( S 1 d ξ 0 R + ( ξ ) Q 1 ω 2 2 | ξ | 2 s 2 + C + Q 1 ( C + ) s 2 d s + | F + | Q 1 ( C + ) ρ + ) + Q ˜ [ r + ] } = S 1 ( k + + f + ) d ξ
and
2 σ R 0 2 { 1 ρ R 0 2 { S 1 d ξ R + ( ξ ) R ( ξ ) Q 1 ω 2 2 | ξ | 2 s 2 + C Q 1 ( C ) s 2 d s + | F | Q 1 ( C ) ρ + ρ ρ + ( S 1 d ξ 0 R + ( ξ ) Q 1 ω 2 2 | ξ | 2 s 2 + C + Q 1 ( C + ) s 2 d s + | F + | Q 1 ( C + ) ρ + ) } + Q ˜ [ r ] } = S 1 ( k + f ) d ξ .
Next, we represent system (16), (35), (42), (43) in the form of vector equation (23) with Φ ( φ ) = ( Φ 1 ± , Φ 2 , Φ 3 ± ) , φ ( r ± , ω , λ ± ) , λ ± Q 1 ( C ± ) ρ ± ,
Φ 1 ± ( φ ) = σ ± R 0 ± 2 Δ S 1 r ± + 2 r ± k ± f ± , Φ 2 ( φ ) = ω S 1 | ξ | 2 d ξ 0 R 0 + + r + Q 1 ω 2 2 | ξ | 2 s 2 + C + s 4 d s + ω S 1 | ξ | 2 d ξ R 0 + + r + R 0 + r Q 1 ω 2 2 | ξ | 2 s 2 + C s 4 d s β , Φ 3 + ( φ ) = 2 σ + R 0 + 2 { 1 ρ + R 0 + 2 ( S 1 d ξ 0 R 0 + + r + Q 1 ω 2 2 | ξ | 2 s 2 + C + Q 1 ( C + ) s 2 d s + | F + | λ + ) Q ˜ [ r + ] } + S 1 ( k + + f + ) d ξ , Φ 3 ( φ ) = 2 σ R 0 2 { 1 ρ R 0 2 { S 1 d ξ R 0 + + r + R 0 + r Q 1 ω 2 2 | ξ | 2 s 2 + C Q 1 ( C ) s 2 d s + | F | λ + ρ ρ + ( S 1 d ξ 0 R 0 + + r + Q 1 ω 2 2 | ξ | 2 s 2 + C + Q 1 ( C + ) s 2 d s + | F + | λ + ) } Q ˜ [ r ] } + S 1 ( k + f ) d ξ .
Then we apply again implicit function theorem to (23). To this end, we calculate the Fréchet derivative of Φ at zero and linearize (23) at this point. As a result, one has
Φ ( 0 ) φ + ψ ( φ ) = 0
with Φ ( 0 ) = ( Φ 1 ± , Φ 2 , Φ 3 ± ) ( 0 ) ,
Φ 1 + ( 0 ) φ = σ + R 0 + 2 Δ S 1 r + + 2 r + + P + ( ρ + ) λ + P ( ρ ) λ , Φ 1 ( 0 ) φ = σ R 0 2 Δ S 1 r + 2 r + P ( ρ ) λ , Φ 2 ( 0 ) φ = ω ρ + S 1 | ξ | 2 d ξ 0 R 0 + s 4 d s + ω ρ S 1 | ξ | 2 d ξ R 0 + R 0 s 4 d s = 8 π 15 ω ρ + R 0 + 5 + ρ R 0 5 R 0 + 5 , Φ 3 + ( 0 ) φ = 2 σ + R 0 + 2 | B R 0 + | λ + ρ + R 0 + 2 4 π P + ( ρ + ) λ + + 4 π P ( ρ ) λ = 4 π 2 σ + 3 ρ + R 0 + P + ( ρ + ) λ + + P ( ρ ) λ , Φ 3 ( 0 ) φ = 2 σ R 0 2 1 ρ R 0 2 | B R 0 | | B R 0 + | λ + ρ ρ + | B R 0 + | λ + 4 π P ( ρ ) λ = 4 π 2 σ R 0 + 3 3 ρ + R 0 4 λ + + 2 σ ( R 0 3 R 0 + 3 ) 3 ρ R 0 4 P ( ρ ) λ ,
and ψ ( φ ) = Φ ( φ ) Φ ( 0 ) φ ψ 1 ± ( φ ) , ψ 2 ( φ ) , ψ 3 ± ( φ ) , where
ψ 1 + ( φ ) = P + ( Q 1 ( C + ) ) P + ( ρ + ) P ( Q 1 ( C ) ) + P ( ρ ) P + ( ρ + ) λ + + P ( ρ ) λ f + , ψ 1 ( φ ) = P ( Q 1 ( C ) ) P ( ρ ) P ( ρ ) λ f , ψ 2 ( φ ) = ω S 1 | ξ | 2 d ξ 0 R 0 + Q 1 ω 2 2 | ξ | 2 s 2 + Q ( ρ + + λ + ) ρ + s 4 d s + ω S 1 | ξ | 2 d ξ R 0 + R 0 + + r + Q 1 ω 2 2 | ξ | 2 s 2 + Q ( ρ + + λ + ) s 4 d s + ω S 1 | ξ | 2 d ξ R 0 + R 0 Q 1 ω 2 2 | ξ | 2 s 2 + Q ( ρ + λ ) ρ s 4 d s + ω S 1 | ξ | 2 d ξ R 0 R 0 + r Q 1 ω 2 2 | ξ | 2 s 2 + Q ( ρ + λ ) s 4 d s β , ψ 3 + ( φ ) = 2 σ + R 0 + 2 { 1 ρ + R 0 + 2 ( S 1 d ξ 0 R 0 + + r + Q 1 ω 2 2 | ξ | 2 s 2 + C + Q 1 ( C + ) s 2 d s + ( | F + | | B R 0 + | ) λ + ) Q ˜ [ r + ] } + S 1 k + P + ( ρ + ) λ + + P ( ρ ) λ + f + d ξ ,
ψ 3 ( φ ) = 2 σ R 0 2 { 1 ρ R 0 2 { S 1 d ξ R 0 + + r + R 0 + r Q 1 ω 2 2 | ξ | 2 s 2 + C Q 1 ( C ) s 2 d s + ρ ρ + ( S 1 d ξ 0 R 0 + + r + Q 1 ω 2 2 | ξ | 2 s 2 + C + Q 1 ( C + ) s 2 d s + ( | F + | | B R 0 + | ) λ + ) + | F | | B R 0 | + | B R 0 + | λ } Q ˜ [ r ] } + S 1 k P ( ρ ) λ + f d ξ .
We recall the notation: X = C ˜ 2 + α × C ˜ 2 + α × R × R × R and Y = C ˜ α × C ˜ α × R × R × R . By means of fixed point theorem, we prove the solvability of system (23) with (44) in X .
The existence of a unique solution r C ˜ 2 + α ( S 1 ) to two first equations in (25) and the estimate for it are discussed in Sec. 2. The operator Φ 2 ( 0 ) is invertible. In addition, we have assumed that the determinant of two equations Φ 3 ± ( 0 ) φ = ψ 3 ± ( φ ) with relations (46) is not equal to zero:
2 σ + 3 ρ + R 0 + P + ( ρ + ) 2 σ ( R 0 3 R 0 + 3 ) 3 ρ R 0 4 P ( ρ ) 2 σ P ( ρ ) R 0 + 3 3 ρ + R 0 4 0 .
Therefore, the vector value operator Φ ( 0 ) : X Y is invertible too and the solution obeys the inequality
φ X c 2 Φ ( 0 ) φ Y ,
where φ X = | r + | C 2 + α ( S 1 ) + | r | C 2 + α ( S 1 ) + | ω | + | λ + | + | λ | ,
Φ ( 0 ) φ Y = | Φ 1 + ( 0 ) φ | C α ( S 1 ) + | Φ 1 ( 0 ) φ | C α ( S 1 ) + | Φ 2 ( 0 ) φ | + | Φ 3 + ( 0 ) φ | + | Φ 3 ( 0 ) φ | . Hence, one can write equation (45) in the form (27).
We estimate the operator ψ in the same way as in Sec. 2. Using formulas similar to (28) for the new functions f ± , we conclude again that
ψ ( φ ) Y c φ X 1 + α + | β | c 3 η φ X α + | β | , ψ ( φ ) ψ ( φ ) Y c 4 η α φ φ X
for two values of the arguments φ = ( r ± , ω , λ ± ) and φ = ( r ± , ω , λ ± ) X such that
φ X , φ X η , η min { 1 , R 0 + / 2 , ( R 0 R 0 + ) / 2 } .
By repeating the arguments of Sec. 2, we deduce that A ψ ( φ ) is a contraction operator if inequality (30) holds. Hence, fixed point theorem guarantees the existence of a unique solution to equation (27) in the ball { φ X η } with η satisfying (29). Therefore (23) with Φ ( φ ) given by (44) is also uniquely solvable for β such that condition (30) holds. This inequality implies an estimate for ε in (36).
Estimate (37) follows from
φ X A ψ ( φ ) Y c 2 c 3 η α φ X + | β | ,
if c 2 c 3 η α 1 2 . □

Conclusions

Thus, we have shown the existence of axially symmetric equilibrium figures for a two-component compressible fluid when pressure function is given by a smooth growing function of fluid density and the data satisfy some condition. Under a small enough angular momentum these equilibrium figures are close to embedded balls. The next stage of investigation of the problem will consist in proving the existence of a global solution to the nonstationary problem for small initial data and its tendency to the stationary solution ( V , ρ ) , as well as in studying the stability of equilibrium figures obtained. A new paper will deal with this investigation.
In addition, we note that our analysis has been carried out for the case of neglecting the gravity of liquids. This situation is realized in space, and our two-layer gas cloud can be considered, for example, as a gaseous planet or another cosmic rotating body.

References

  1. Lyapunov, A. M.: On stability of ellipsoidal shapes of equilibrium of revolving liquid, Editiion of the Academy of Scien., 1884 (in Russian).
  2. Lyapunov, A. M.: Sur les questions qui appartiennent aux surfaces des figures d’equilibre dérivées des ellopsoïdes, News of the Academy of Scien., (1916) p. 139.
  3. Poincaré, H.: Figures d’équilibre d’une mass fluide, Paris, Gautier–Villars, 1902.
  4. Globa-Mikhailenko, B.: Figures ellipsoïdales d’équilibre d’une masse fluide en rotation quand on tient compte de la pression capillaire, Comptes rendus, 160, 233 (1915).
  5. Charrueau, A. Ètude d’une masse liquide de révolution homogène, sans pesanteur et à tension superficielle, animée d’une rotation uniforme. Ann. de Ècole Normale Supérieure 1926, 43, 129–176. [Google Scholar] [CrossRef]
  6. Charrueau, A. Sur les figures d’équilibre relatif d’une masse liquide en rotation à tension superficielle. C. R. 1927, 184, 1418. [Google Scholar]
  7. Appell, P. Figures d’Équilibre d’une Mass Liquide Homogène en Rotation—Traité de Mécanique Rationnelle, v. IV, Fasc. I, 2nd ed.; Gautier–Villars: Paris, France, 1932. [Google Scholar]
  8. Chandrasekhar, S. Ellipsoidal Figures of Equilibrium, Yale University Press: New Haven and London, 1969.
  9. Solonnikov, V. A.; Tani, A. Equilibrium figures of slowly rotating viscous compressible barotropic capillary liquid, Advances in Math. Sciences and Applications, 1993 2(1), 139–145.
  10. Denisova, I. V. Problem on the motion of two compressible fluids separated by a closed free interface, Zap. Nauch. Semin. POMI 1997, 243, 61–86; English transl. J. Math. Sci. 2000, 99, 837–853. [Google Scholar] [CrossRef]
  11. Denisova, I. V. Solvability in weighted Hölder spaces of a problem governing the evolution of two compressible fluids, Zap. Nauch. Semin. POMI 2003 295, 57–89; English transl. in J. Math. Sci. 2005 127(2), 1849–1868.
  12. Solonnikov, V. A.; Tani, A. Free boundary problem for a viscous compressible flow with surface tension, in: Constantin Carathéodory: An International Tribute, World Scientific (1991), 1270–1303.
  13. Denisova, I.V.; Solonnikov, V.A. Classical solvability of a model problem in the half-space related to the motion of an isolated mass of a compressible liquid. Zap. Nauchn. Semin. Petersburg. Otdel. Mat. Inst. Steklov. (POMI). 2000,271, 92–113 (Russian); English transl. J. Math. Sci. 2003, 115(6), 2753–2765.
  14. Blaschke, W. Vorlesungen über Differentialgeometrie und geometrische Grundlagen von Einsteins Relativitätstheorie. I; Springer: Berlin, Germany, 1924. [Google Scholar]
  15. Solonnikov, V.A. On the problem of non-stationary motion of two viscous incompressible liquids. Probl. Mat. Analiz. 2006, 34, 103–121, Engl. transl. J. Math. Sci. 2007, 142, 1844–1866. [Google Scholar] [CrossRef]
  16. Solonnikov, V. A. Unsteady motions of a finite isolated mass of a self-gravitating fluid, Algebra i Analiz, 1989, 1(1), 207–249; English transl. in Leningrad Math. J. 1990, 1(1), 227–276.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.
Copyright: This open access article is published under a Creative Commons CC BY 4.0 license, which permit the free download, distribution, and reuse, provided that the author and preprint are cited in any reuse.
Prerpints.org logo

Preprints.org is a free preprint server supported by MDPI in Basel, Switzerland.

Subscribe

© 2024 MDPI (Basel, Switzerland) unless otherwise stated