Next Article in Journal
A Review of Combinatorial Optimization Problems in Reverse Logistics and Remanufacturing for End-of-Life Products
Next Article in Special Issue
Optimization Method for Solving Cloaking and Shielding Problems for a 3D Model of Electrostatics
Previous Article in Journal
Topological Indices of Graphs from Vector Spaces
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Multi-Shell Models of Celestial Bodies with an Intermediate Layer of Fluid: Dynamics in the Case of the Large Values of the Ekman Number

by
Vladislav Sidorenko
1,2,* and
Sergey Ramodanov
3
1
Keldysh Institute of Applied Mathematics RAS, Moscow 125047, Russia
2
Sofia Kovalevskaya Northwestern Center for Mathematical Research, Pskov State University, Pskov 180000, Russia
3
Mechanical Engineering Research Institute RAS, Moscow 101000, Russia
*
Author to whom correspondence should be addressed.
Mathematics 2023, 11(2), 296; https://doi.org/10.3390/math11020296
Submission received: 7 December 2022 / Revised: 2 January 2023 / Accepted: 3 January 2023 / Published: 6 January 2023
(This article belongs to the Special Issue Theoretical Research and Computational Applications in Fluid Dynamics)

Abstract

:
We consider a mechanical system that is comprised of three parts: a rigid outer shell with a spherical cavity, a spherical core inside this cavity, and an intermediate layer of liquid between the core and the shell. Such a model provides an adequate description of the behavior of a wide variety of celestial bodies. The centers of the inner and outer liquid’s spherical boundaries are assumed to coincide. Assuming that the viscosity of the liquid is high, we obtained an approximate solution to the Navier–Stokes equations that describes a so called creeping flow of the liquid, which sets on after all transient processes die out. We note that the effect of the liquid on the rotational motion of the system can be modeled as a special torque acting upon the system with “solidified” fluid.
MSC:
74F10; 76D05; 76U05; 76M45

1. Introduction

It is customary to model a celestial object as a three-component structure comprised of an inner core, an intermediate layer of liquid, and an outer solid shell [1,2,3,4,5,6,7,8]. We will deal with the simplest, but fairly common version of such an approach, which does not take into account that in reality the interfaces between liquid and solid components represent a mushy transition zone [9,10]. The mass proportion of each component strongly depends on the object’s type. For example, the core of an Earth-like planet is considerably smaller than the outer shell (mantle) (Figure 1a). At the same time, the icy crust, covering the hypothetical subsurface oceans of some satellites of Jupiter and Saturn, is relatively thin as compared to the satellites’ radii (Figure 1b) [11]. Studying the dynamics of the multi-component object it is convenient to fix a coordinate frame to the most dominant component (either outer shell or inner core) and investigate the behavior of the other components relative to this frame.
In the theory of the rotating mass of liquid, the key role plays Ekman number [12], which reads
E = μ ω r o t ρ f R c 2 .
Here, ω r o t is the characteristic angular velocity of system’s components, R c is the characteristic size (e.g., the radius of the outer boundary of the liquid layer), ρ f is the density of the liquid, and μ is the liquid’s viscosity.
If E 1 , the flow of the liquid within the layer is known to be Stokesian or creeping [13]. In this case, the approximate relations for the liquid’s velocities can be obtained from the Navier–Stokes equations with the relative-acceleration terms neglected. Stokes flow in the gap between two concentric spheres has been considered by many experts. However, in the best part of the studies the motion of the bounding surfaces was assumed to be known, and the effect the liquid invokes on their dynamics was not studied (e.g., see [14,15]).
In [16], F.L. Chernousko showed that when E 1 the analysis of the dynamics of a body containing a cavity filled with liquid can be simplified: the liquid in the cavity is considered “solidified”, and the liquid’s effect on the body is modeled as an additional torque applied to the body. Various aspects of the Chernousko approach are discussed in [17,18].
In the current paper, we show that in the case E 1 , a similar approach can be used for analyzing the rotation of celestial objects containing liquid layers. A useful byproduct of such an approach is an approximate formula for the velocity field of the liquid particles.
As far as we know, there are no multi-component celestial bodies with E 1 . Therefore, at least at present, our results seem to have only a peripheral relationship to the dynamics of actual astrophysical objects. At the same time, our results may be used, for example, for testing software that simulates the dynamics of such objects when it is important to be able to compare numerical results with a known analytical solution, regardless of whether this solution represents a real phenomena [19].
The paper is organized as follows. In Section 2, the governing equations for a three-component system mantle+liquid+core are presented. An approximate formula for the liquid velocity field after Stokes (creeping) flow sets in is obtained in Section 3. A condition of “dynamical agreement” between the motion of the core and the motion of the liquid near the core’s surface is presented in Section 4. This condition helps to establish all of the parameters of the liquid flow. Section 5 discusses some aspects of this flow. In Section 6, the approximate relations describing the regular precession of the mantle relative to the core are given. Some issues associated with the introduction of the additional torque that replicates the action of the liquid (the Chernousko model) are covered in Section 7. Modifications to the formulas for the liquid velocity to suit the case where the thin outer shell moves relative to the massive core are discussed in Section 8. A brief outline of the results obtained is presented in Section 9.

2. Equations of Motion of a “Planet”

We assume that the external forces and torques applied to the system allow for a motion in the course of which the center of mass of the outer shell (referred to further as mantle) coincides with the center of mass of the inner core. Such a regime (at least approximately) occurs when a “planet” orbits a star. Under this assumption, the motion of the core relative to a mantle-fixed frame O x y z is a rotation with angular velocity ω s c r . For simplicity, we assume that the motion of the centers of mass is known (i.e., r O and v O = d r O / d t are known functions of time). Relative to the frame O x y z , the motion of the liquid is governed by the Navier–Stokes equations
1 ν u t + ( u ) u + 2 ω × u + d ω d t × r Δ u + q = 0 ,
q = 1 ν p ρ f + Π 1 2 | ω × r | 2 r , d v O d t .
Here, u ( r , t ) is the liquid velocity, p ( r , t ) is the pressure in the liquid layer, Π is the potential of the gravity attraction forces exerted on the liquid particles, ω is the angular velocity of the frame O x y z , and ν = μ / ρ f is the liquid’s kinematic viscosity. The velocity u ( r , t ) satisfies the standard incompressibility condition
div u = 0
and the kinematic relations
u | r | = R l c = 0 .
u | r | = R s c = ω s c r × r ,
where R l c and R s c stand for the radius of the cavity and the core, respectively.
Within the core, the distribution of the matter’s density is assumed to be spherically symmetric. In the frame O x y z , the evolution of its rotational speed is governed by the equation
I s c d ω s c r d t = M f l u i d I s c d ω d t .
Here, M f l u i d is the friction torque that the liquid exerts on the core. To find M f l u i d , we integrate the tangential stress components over the core’s surface S s c , that is,
M f l u i d = S s c r × T · n d S .
Here, T is the stress tensor and n is the unit outward normal to the surface of the core.
The total angular momentum K p Σ relative to the center of mass can be decomposed as
K p Σ = K p s + K f l u i d r + K s c r .
Here, K p s is the angular momentum of the “solidified” planet in which the liquid and the inner core are fixed relative to the mantle; the internal momentum of the liquid can be written as K f l u i d r = ρ f V l c r × u d V ; the attitude momentum of the core reads K s c r = I s c ω s c r ; and K p s = I p ω , where I p stands for the tensor of inertia of the “solidified” planet. Let M Σ be the net external torque applied to the system; then, the time evolution of the total angular momentum is described by the Euler equation
d d t K p S + K f l u i d r + K s c r + ω × K p S + K f l u i d r + K s c r = M Σ .
The governing equations should be augmented with kinematic relations of the form
d z d t = k ( z , ω ) ,
where the parameters z quantify the orientation of the mantle-fixed frame (e.g., the Euler angles).
Hereafter, the units of mass, length, and time are assumed to be so chosen that ρ f = 1 , R l c R s c 1 and therefore ν E 1 .

3. Stokes Flow in the Liquid Layer

If E 1 , after an initial transient process, the flow in the cavity becomes creeping flow [13]. Omitting the inertial terms in (1) yields the following equation governing the creeping flow:
Δ u = q + d ω d t × r ν ,
where the velocity u satisfies the incompressibility condition (2) and the kinematic relations (3) and (4).
We seek a solution to (8) of the form
u = ν 1 d ω d t × f ,
which upon substitution into (8) provides the following equation in the function f ( r )
Δ f = r .
The velocity field satisfies the boundary condition u | r | = R l c = 0 provided that
f | r | = R l c = 0 .
To meet the boundary condition (3), the rotation of the inner core and the velocity (9) must agree kinematically, that is,
ω s c r = α ν d ω d t ,
f | r | = R s c = α r .
Moreover, Equation (5) implies a certain “dynamical” agreement between the motion of the liquid and the motion of the inner core. In Section 4, we will derive the value of α for which (12) provides an approximate solution to (5), thereby revealing the dynamics of all of the system’s components. For simplicity, put
f | r | = R l c = 0 .
By symmetry, it is quite natural to assume that the solution to the boundary-value problems (10)–(12) and (14) lies among functions that depend solely on r = | r | . For such functions, this problem reduces to the boundary-value problem for the third-order ODE
d d r r 2 d 2 f d r 2 2 d f d r = r 3
of the form
f ( R l c ) = d f d r ( R l c ) = 0 , d f d r ( R s c ) = α R s c .
The solution to the boundary-value problems (15) and (16) looks like
f ( r ) = r 4 40 + β 1 r 2 + β 2 r + β 3 ,
where the constants β 1 , β 2 , β 3 satisfy the linear relations
R l c 4 40 + β 1 R l c 2 + β 2 R l c + β 3 = 0 ,
R l c 2 10 + 2 β 1 β 2 R l c 3 = 0 , R s c 2 10 + 2 β 1 β 2 R s c 3 = 0 .
It is straightforward to obtain
f ( r ) = 1 40 r 2 R l c 2 2 + β 2 1 r 1 R l c + 1 2 R l c 3 r 2 R l c 2
and
α = 1 10 R l c 2 R s c 2 β 2 1 R l c 3 1 R s c 3

4. “Dynamical” Agreement between the Motion of the Inner Core and the Liquid

From Equation (5), it follows that after the creeping flow onset one has
M f l u i d I s c d ω d t
To obtain M f l u i d , one should calculate the integral (6). Following [14], introduce spherical coordinates whose polar axis is aligned with d ω d t . The velocity components of a liquid particle are
u r = u θ = 0 , u φ = ν 1 r 3 10 + 2 β 1 r β 2 r 2 d ω d t sin θ .
The force per unit area of the surface of the inner core can be determined with the help of the τ r φ -component of the stress tensor T
τ r φ = ν u φ r u φ r r = R s c = R s c 2 5 + 3 β 2 R s c 3 d ω d t sin θ .
Integrating over the sphere’s surface yields
M f l u i d = 8 π 3 ρ R s c 3 R s c 2 5 + 3 β 2 R s c 3 d ω d t = I f l u i d 1 + 15 β 2 R s c 5 d ω d t ,
where I f l u i d is the central moment of inertia of the cavity of radius R s c filled with liquid with density ρ f = 1 .
Equating (18) and (19) we obtain
β 2 = R s c 2 15 I s c I f l u i d 1 .
With this expression for β 2 at our disposal, it is easy to find that
α = R s c 2 Ξ s c R s c R l c , I f l u i d I s c
Ξ s c ( ξ , η ) = 1 15 · 1 η 1 1 ξ 3 + 1 10 · 1 ξ 2 1
The graph of Ξ s c is given in Figure 2. It can be easily checked that α > 0 when I f l u i d < I s c , which is quite natural since usually the inner core is more dense than the liquid.

5. Common Properties of the Flow in the Liquid Layer

Using (9), (17) and (20), it is possible to draw the lines of constant velocity magnitude in a plane through the axis of relative rotation of the inner core (Figure 3). Obviously, the magnitude acquires its maximum in the ”equatorial” plane, which is perpendicular to the relative angular velocity of the core and passes through its center.
The most trivial distribution of the velocity magnitude occurs when
I s c = I f l u i d .
Assuming that the core is uniform in density, Equation (21) implies that its density is equal to the density of the ambient liquid. So, if (21) is fulfilled, then in the equatorial plane
u φ = 1 10 ν r 2 R l c 2 r .
Thus, if R s c R l c 3 the maximum of the magnitude is reached on the surface of the inner core, whereas for R s c < R l c 3 the maximum occurs at r = R l c 3 .

6. Motion of the Core Relative to the Mantle

The attitude dynamics of the core is governed by Equation (5). Let us examine the dynamics in greater detail assuming that the mantle is a dynamically symmetric rigid body that is in a precession motion with parameters ω ψ , ω φ and θ ; as usual, ω ψ is the precession angular velocity, ω φ stands for the intrinsic angular velocity, and θ is the angle between the mantle’s axis of symmetry and the precession axis. Suppose that the O z -axis of the mantle-fixed frame is aligned with the symmetry axis; then, in this frame the angular velocity of the mantle and the relative angular velocity of the core (up to an inessential time shift) can be written as
ω = ω ψ sin θ sin ω φ t , ω ψ sin θ cos ω φ t , ω ψ cos θ + ω φ T ,
ω s c r α ν d ω d t = α ν ω ψ ω φ sin θ cos ω φ t , ω ψ ω φ sin θ sin ω φ t , 0 T .
Let three orthogonal unit vectors e x s c , e y s c , e z s c be fixed to the core. The kinematic equations
d e x s c d t = ω s c r × e x s c , d e y s c d t = ω s c r × e y s c , d e z s c d t = ω s c r × e z s c ,
admit a solution
e x s c = 1 + O ν 2 , O ν 2 , α ω ψ ν sin θ cos ω φ t + O ν 2 T
e y s c = O ν 2 , 1 + O ν 2 , α ω ψ ν sin θ sin ω φ t + O ν 2 T
e z s c = α ω ψ ν sin θ cos ω φ t + O ν 2 , α ω ψ ν sin θ sin ω φ t + O ν 2 , 1 + O ν 2 T
This solution can be characterized as the core “swaying” with zero secular drift with respect to the mantle. The other solutions to (23) can be obtained from (24) by the action of the rotation group.
Remark. The secular drift of the core still can occur in the course of the mantle’s moderate despinning or acceleration. For example, suppose α > 0 and the mantle gradually decelerates; then, with respect to the frame O x y z , the inner core rotates with angular velocity ω s c r α ν d ω d t in the direction opposite the vector of angular acceleration. In a sense, this mimics the superrotation of Earth’s inner core [20]. However, the current concept is that quantitatively reliable estimates of the superrotation can only be obtained with due regard for dynamo effects and turbulent convection within the liquid layer.

7. Chernousko’s Approach to Modelling Rotational Motion of the System

The total attitude angular momentum reads
K r = ν 1 V l c r × d ω d t × f d V + I s c ω s c r
The integral evaluates as follows:
V l c r × d ω d t × f d V = V l c ( r , f ) d ω d t f · r , d ω d t d V = V l c ( r , f ) + f d ω d t f · r , d ω d t d V = 2 V l c f d V · d ω d t = κ l c d ω d t .
Here,
κ l c = π 525 · { 15 R l c 7 R s c 7 42 R l c 2 R l c 5 R s c 5 + 30 R l c 4 R l c 3 R s c 3 +
+ 24 R s c 5 I s c I f l u i d 1 5 R l c 2 R s c 2 + 1 R l c 3 R l c 5 R s c 5 + 5 R l c R l c 3 R s c 3 }
From (12) and (26), it follows that
K r = κ Σ ν d ω d t , κ Σ = κ l c + I s c R s c 2 Ξ s c .
The approach proposed by F.L. Chernousko to modeling the dynamics of a rigid body with a cavity entirely filled with a highly viscous liquid is as follows: the liquid is assumed to be frozen or “solidified”; the influence of the liquid is then represented as a special torque acting upon the body with solidified liquid [16]. In the context of this approach, the Euler equations can be rewritten as
d K p r d t + ω × K p r = M Σ + M c ,
where
M c = d K r d t ω × K r κ Σ ν d 2 ω d t 2 + ω × d ω d t .
Suppose that the derivatives d ω / d t and d 2 ω / d t 2 in (28) are found for the “zero approximation” model, that is, when the “planet” is assumed to be a one-piece uniform rigid body; then, Equation (27), together with the kinematic relations (7), forms a finite-dimensional dynamical system that governs our originally infinite-dimensional system. Such models (obtained within Chernousko’s approach) are as accurate as ν 1 during time ν [17,18].

8. Simple Dynamical Model of an “Icy Satellite” with Subsurface Ocean

We now turn to the dynamics of the so-called membrane worlds encountered in the satellite systems of the giant planets [11]. Consider a satellite in the shape of a sphere of radius R o b about which a thin icy shell (crust) is floating on a liquid layer (subsurface ocean). The inner radius of the icy crust will be further denoted as R i s .
By a similar argument as in Section 3 and Section 4, it can be shown that the creeping flow of the ocean is described by formula (9), in which
f ( r ) = 1 40 r 2 R o b 2 2 R i s 5 15 1 + I i s I f l u i d 1 r 1 R o b + 1 2 R o b 3 r 2 R o b 2 .
Here, I i s is the moment of inertia of the icy crust, and I f l u i d = 8 15 π R i s 5 is the moment of inertia of a ball of radius R i s , whose density is equal to the density of the liquid; both moments are about an axis through the center of mass.
Lines of equal velocity magnitude in a plane through the axis of relative rotation of the crust are shown in Figure 4.
The relative angular velocity of the crust is ω i s r = α ν d ω d t , where
α = R i s 2 Ξ i s R o b R i s , I i s I f l u i d
Ξ i s ( ξ , η ) = 1 15 · ( 1 + η ) 1 ξ 3 1 1 10 · 1 ξ 2 .
It is easy to verify that for ξ = 1 Δ , 0 < Δ 1 , the following relation holds
Ξ i s ( 1 Δ , η ) η Δ 5 > 0 .
This means that for a relatively shallow ocean ( R i s R o b R i s ), the coefficient α is positive. Thus, in this case when a satellite moves with angular acceleration (or deceleration), the icy crust will drift relative to the core with angular velocity ω b r = | α | ν d ω d t in the direction opposite the vector d ω d t . The drift of the icy crust of some giant-planet satellites is discussed in, for example, [21,22]. Since for these objects the condition E 1 does not hold, our results can hardly help to obtain new observationally confirmed features of their behavior.

9. Conclusions

We studied the motion of the liquid in the layer between two concentric spherical surfaces. In the case where the motion of one surface is solely determined by the moment of the tangential stress, an approximate expression for the Stokes flow of the liquid is derived. Using this expression, we identified an important class of motions of the three-component system mantle+liquid+core, namely, where the motion of one of the bodies and the liquid are subordinate to the motion imparted by the other body. We showed that such a subordination can be used to obtain a finite dimensional model of the system (the effect of intrinsic motion of liquid particles is modelled as a torque applied to “solidified” system [16]).
Although our study was inspired by real models used to investigate the dynamics of astrophysical objects, our main results were obtained under assumptions that can hardly be met within these models’ framework. Nevertheless, it is worth noting that in the classical problem of describing a creeping flow between two concentric spherical surfaces, we have found a new class of approximate solutions that is valuable in its own right. For example, such solutions can be useful for testing special flow analysis software (as mentioned in the Introduction).
Future research may include the application of perturbation theory for obtaining similar results when the surfaces bounding the flow are not spherical [23,24]. Additionally, our approach may find its application in the study of porous bodies floating within liquid-filled cavities [25,26].

Author Contributions

Conceptualization, V.S. and S.R.; Investigation, V.S. and S.R. All authors have read and agreed to the published version of the manuscript.

Funding

V.V. Sidorenko acknowledges the financial support of his participation in this study by the Scientific and Educational Mathematical Center “Sofia Kovalevskaya Northwestern Center for Mathematical Research” (contract 075-02-2022-891 from 31 January 2022). The work of S.M. Ramodanov is supported by RFBR grant 20-01-00312A “Resonance effects and chaos in the three-body problem”.

Data Availability Statement

Not applicable.

Acknowledgments

The authors thank the anonymous reviewers for helpful comments.

Conflicts of Interest

The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.

References

  1. Busse, F.H. On the free oscillations of the Earth’s inner core. J. Geophys. Res. 1974, 79, 753–757. [Google Scholar] [CrossRef]
  2. Dehant, V.; Van Hoolst, N.; de Veres, O.; Greff-Lefftz, M.; Legros, H.; Defraigne, P. Can a solid inner core of Mars be detected from observations of polar motion and nutations of Mars? J. Geophys. Res. 2003, 108, 1.1–1.12. [Google Scholar] [CrossRef]
  3. Getino, J.; Farto, J.M.; Ferrandiz, J.M. Obtaining the free frequencies of the non-rigid Earth. Celest. Mech. Dyn. Astr. 1999, 71, 95–108. [Google Scholar] [CrossRef]
  4. Grinfeld, P.; Wisdom, J. Motion of the mantle in the translational modes of the Earth and Mercury. Phys. Earth Planet. Inter. 2005, 151, 77–87. [Google Scholar] [CrossRef]
  5. Hussmann, H.; Sohl, F.; Spohn, T. Subsurface oceans and deep interiors of medium-sized outer planet satellites and large transneptunian objects. Icarus 2006, 185, 258–273. [Google Scholar] [CrossRef]
  6. Kitiashvili, I.; Gusev, A. Inner core wobble and free core nutation of pulsar PSR B1828-11. Adv. Space Res. 2003, 42, 1391–1397. [Google Scholar] [CrossRef]
  7. Nimmo, F.; Rappalardo, R.T. Ocean worlds in outer Solar system. J. Geophys. Res. Planets 2016, 121, 1378–1399. [Google Scholar] [CrossRef]
  8. Shematovich, V.I. Ocean worlds in the outer regions of the Solar system (review). Sol. Syst. Res. 2018, 52, 371–381. [Google Scholar] [CrossRef]
  9. Loper, D.E.; Fearn, D.R. A seismic model of a pertial molten inner core. J. Geophys. Res. 1983, 88, 1235–1242. [Google Scholar] [CrossRef]
  10. Peng, Z.R. Effects of a mushy transition zone at the inner core boundary on the Slicher modes. Geophys. J. Int. 1997, 131, 607–617. [Google Scholar] [CrossRef]
  11. Beuthe, M. Tidal Love numbers of membrane worlds: Europa, Titan, and Co. Icarus 2016, 258, 239–266. [Google Scholar] [CrossRef] [Green Version]
  12. Greenspan, H.P. The Theory of Rotating Fluids; Cambridge University Press: Cambridge, UK, 1968. [Google Scholar]
  13. Happel, J.; Brenner, H. Low Reynolds Number Hydrodynamics: With Special Applications to Particulate Media; Prentice-Hall: Englewood Cliffs, NJ, USA, 1965. [Google Scholar]
  14. Landau, L.D.; Lifshitz, E.M. Fluid Mechanics, 2nd ed.; Pergamon Press: Oxford, UK, 1987. [Google Scholar]
  15. Shankar, P. Exact solutions for Stokes flow in and around a sphere and between concentric spheres. J. Fluid Mech. 2009, 631, 363–373. [Google Scholar] [CrossRef]
  16. Chernousko, F.L. Motion of a Rigid Body with Cavities Containing a Viscous Fluid; NASA Technical Translations; Computing Centre of the USSR Academy of Sciences: Moscow, Russia, 1972. [Google Scholar]
  17. Bogatyrev, S.V. Slow motions in problems of the dynamics of a solid with a cavity filled with a viscous liquid. J. Applied Math. Mech. 1994, 58, 849–855. [Google Scholar] [CrossRef]
  18. Ramodanov, S.M.; Sidorenko, V.V. Dynamics of a rigid body with an ellipsoidal cavity filled with viscous fluid. Int. J. Non-Linear Mech. 2017, 95, 42–46. [Google Scholar] [CrossRef]
  19. Child, A.; Hollerbach, R.; Kersale, E. Axisymmetric pulse train solutions in narrow-gap spherical Couette flow. Phys. D Nonlinear Phenom. 2017, 348, 54–59. [Google Scholar] [CrossRef] [Green Version]
  20. Song, X.; Richards, R.G. Seismological evidence for differential rotation of the Earth’s inner core. Nature 1996, 382, 221–224. [Google Scholar] [CrossRef]
  21. Goldreich, P.M.; Mitchell, J.L. Elastic ice shells of synchronous moons. Implications for cracks on Europa and non-synchronous rotation of Titan. Icarus 2010, 209, 631–638. [Google Scholar] [CrossRef] [Green Version]
  22. Patthoff, D.A.; Kattenhorn, S.A.; Cooper, C.M. Implications of nonsynchronous rotation on the deformational history and ice shell properties in the south polar terrain of Enceladus. Icarus 2019, 321, 445–457. [Google Scholar] [CrossRef]
  23. Denisov, G.G.; Novikov, V.V. On fluid flow between rotating surfaces. J. Appl. Mech. Tech. Phys. 2011, 52, 31–36. [Google Scholar] [CrossRef]
  24. Prasad, M.K.; Kaur, M.; Srinivasacharya, D. Slow steady rotation of an approximate sphere in an approximate spherical container with slip surfaces. Int. J. Appl. Comput. Math. 2017, 3, 987–999. [Google Scholar] [CrossRef]
  25. Keh, H.J.; Lu, Y.S. Creeping motion of a porous spherical shell in concentric spherical cavity. J. Fluid. Struct. 2005, 20, 735–747. [Google Scholar] [CrossRef]
  26. Verma, V.K.; Verma, H. Motion of a porous spherical shell in a spherical container. Spec. Top. Rev. Porous Media 2019, 10, 525–537. [Google Scholar] [CrossRef]
Figure 1. Celestial objects modeled as multi-shell structures: (a) “planet” (small inner core), (b) “icy satellite” with thin outer layer of ice, covering the subsurface ocean. In the case of a “planet,” the radii of the inner and outer boundaries of the liquid layer are denoted by R s c and R l c ; in the case of an icy satellite, they are denoted by R o b and R i s , respectively. The visible contours of a body can have any shape; the only thing that matters is that the center of mass of the outer component coincides with the center of the cavity.
Figure 1. Celestial objects modeled as multi-shell structures: (a) “planet” (small inner core), (b) “icy satellite” with thin outer layer of ice, covering the subsurface ocean. In the case of a “planet,” the radii of the inner and outer boundaries of the liquid layer are denoted by R s c and R l c ; in the case of an icy satellite, they are denoted by R o b and R i s , respectively. The visible contours of a body can have any shape; the only thing that matters is that the center of mass of the outer component coincides with the center of the cavity.
Mathematics 11 00296 g001
Figure 2. The graph of the function Ξ s c ( ξ , η ) . The argument ξ characterizes the ratio of the radii of the inner and outer boundaries of the liquid layer, and the argument η represents the ratio of the densities of the liquid and the substance of the inner core (assumed homogeneous).
Figure 2. The graph of the function Ξ s c ( ξ , η ) . The argument ξ characterizes the ratio of the radii of the inner and outer boundaries of the liquid layer, and the argument η represents the ratio of the densities of the liquid and the substance of the inner core (assumed homogeneous).
Mathematics 11 00296 g002
Figure 3. Velocity level lines in a plane through the instantaneous axis of relative rotation of the inner core (which is also the flow’s axis of symmetry) when I s c / I f l u i d = 1 : (a) R s c / R l c = 0.6 > 1 / 3 (the velocity of the liquid reaches its maximum on the surface of the inner core), (b) R s c / R l c = 0.3 < 1 / 3 (the maximal velocity of the flow is attained inside the layer).
Figure 3. Velocity level lines in a plane through the instantaneous axis of relative rotation of the inner core (which is also the flow’s axis of symmetry) when I s c / I f l u i d = 1 : (a) R s c / R l c = 0.6 > 1 / 3 (the velocity of the liquid reaches its maximum on the surface of the inner core), (b) R s c / R l c = 0.3 < 1 / 3 (the maximal velocity of the flow is attained inside the layer).
Mathematics 11 00296 g003
Figure 4. Distribution of velocities in the ocean: lines of equal velocity magnitude in a plane through the instantaneous axis of relative rotation of the icy crust (which is the flow’s axis of symmetry) for R o b / R i s = 0.7 , I i s / I f l u i d = 0.05 , α 0.08308 .
Figure 4. Distribution of velocities in the ocean: lines of equal velocity magnitude in a plane through the instantaneous axis of relative rotation of the icy crust (which is the flow’s axis of symmetry) for R o b / R i s = 0.7 , I i s / I f l u i d = 0.05 , α 0.08308 .
Mathematics 11 00296 g004
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.

Share and Cite

MDPI and ACS Style

Sidorenko, V.; Ramodanov, S. Multi-Shell Models of Celestial Bodies with an Intermediate Layer of Fluid: Dynamics in the Case of the Large Values of the Ekman Number. Mathematics 2023, 11, 296. https://doi.org/10.3390/math11020296

AMA Style

Sidorenko V, Ramodanov S. Multi-Shell Models of Celestial Bodies with an Intermediate Layer of Fluid: Dynamics in the Case of the Large Values of the Ekman Number. Mathematics. 2023; 11(2):296. https://doi.org/10.3390/math11020296

Chicago/Turabian Style

Sidorenko, Vladislav, and Sergey Ramodanov. 2023. "Multi-Shell Models of Celestial Bodies with an Intermediate Layer of Fluid: Dynamics in the Case of the Large Values of the Ekman Number" Mathematics 11, no. 2: 296. https://doi.org/10.3390/math11020296

APA Style

Sidorenko, V., & Ramodanov, S. (2023). Multi-Shell Models of Celestial Bodies with an Intermediate Layer of Fluid: Dynamics in the Case of the Large Values of the Ekman Number. Mathematics, 11(2), 296. https://doi.org/10.3390/math11020296

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop