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
Here,
is the characteristic angular velocity of system’s components,
is the characteristic size (e.g., the radius of the outer boundary of the liquid layer),
is the density of the liquid, and
is the liquid’s viscosity.
If
, 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
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 , 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
. 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
is a rotation with angular velocity
. For simplicity, we assume that the motion of the centers of mass is known (i.e.,
and
are known functions of time). Relative to the frame
, the motion of the liquid is governed by the Navier–Stokes equations
Here,
is the liquid velocity,
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
, and
is the liquid’s kinematic viscosity. The velocity
satisfies the standard incompressibility condition
and the kinematic relations
where
and
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
, the evolution of its rotational speed is governed by the equation
Here,
is the friction torque that the liquid exerts on the core. To find
, we integrate the tangential stress components over the core’s surface
, that is,
Here,
T is the stress tensor and
is the unit outward normal to the surface of the core.
The total angular momentum
relative to the center of mass can be decomposed as
Here,
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
; the attitude momentum of the core reads
; and
, where
stands for the tensor of inertia of the “solidified” planet. Let
be the net external torque applied to the system; then, the time evolution of the total angular momentum is described by the Euler equation
The governing equations should be augmented with kinematic relations of the form
where the parameters
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 and therefore
3. Stokes Flow in the Liquid Layer
If
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:
where the velocity
satisfies the incompressibility condition (2) and the kinematic relations (3) and (4).
We seek a solution to (8) of the form
which upon substitution into (8) provides the following equation in the function
The velocity field satisfies the boundary condition
provided that
To meet the boundary condition (3), the rotation of the inner core and the velocity (9) must agree kinematically, that is,
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
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
. For such functions, this problem reduces to the boundary-value problem for the third-order ODE
of the form
The solution to the boundary-value problems (15) and (16) looks like
where the constants
satisfy the linear relations
It is straightforward to obtain
and
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
To obtain
, one should calculate the integral (6). Following [
14], introduce spherical coordinates whose polar axis is aligned with
. The velocity components of a liquid particle are
The force per unit area of the surface of the inner core can be determined with the help of the
-component of the stress tensor
TIntegrating over the sphere’s surface yields
where
is the central moment of inertia of the cavity of radius
filled with liquid with density
.
Equating (18) and (19) we obtain
With this expression for
at our disposal, it is easy to find that
The graph of
is given in
Figure 2. It can be easily checked that
when
, 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
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
Thus, if the maximum of the magnitude is reached on the surface of the inner core, whereas for the maximum occurs at .
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
-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
Let three orthogonal unit vectors
be fixed to the core. The kinematic equations
admit a solution
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
and the mantle gradually decelerates; then, with respect to the frame
, the inner core rotates with angular velocity
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
The integral evaluates as follows:
Here,
From (12) and (26), it follows that
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
where
Suppose that the derivatives
and
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
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
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
.
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
Here,
is the moment of inertia of the icy crust, and
is the moment of inertia of a ball of radius
, 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
, where
It is easy to verify that for
, the following relation holds
This means that for a relatively shallow ocean (
), 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
in the direction opposite the vector
The drift of the icy crust of some giant-planet satellites is discussed in, for example, [
21,
22]. Since for these objects the condition
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].