The precession and nutation of deformable bodies,II

In a previous paper of this series (Kopal, 1968a) the Eulerian equations have been set up which govern the precession and nutation of selfgravitating bodies of viscous fluid in inertial coordinates which are at rest in space. In order to facilitate their solution, in the present investigation we shall transform these equations to the rotating body-axes; and shall explicitly evaluate all their coefficients arising as a result of second-harmonic dynamical tides.Following the introductory Section 1 which contains a mathematical statement of the problem, the requisite transformation of coordinates will be outlined in Section 2, and applied to the equations of motion in Section 5. The corresponding moments and products of inertia appropriate for selfgravitating configurations of arbitrary internal structure will be formulated in Section 4; while the deformation terms arising from second-harmonic dynamical tides raised on centrally-condensed configurations will be evaluated in Sections 3 and 6. The concluding Section 7 will then contain a specification of the components of the disturbing force.The next stage of our investigation — namely, a construction of the actual solutions of the equations governing precession and nutation of fluid bodies in different cases of astrophysical interest — has been postponed for a separate paper.     

Dynamical tides in close binary systems

The aim of the present paper will be to investigate the circumstances under which an irreversible dissipation of the kinetic energy into heat is generated by the dynamical tides in close binary systems if (a) their orbit is eccentric; (b) the axial rotation of the components is not synchronized with the revolution; or (c) the equatorial planes are inclined to that of the orbit.In Section 2 the explicit form of the viscous dissipation function will be set up in terms of the velocity-components of spheroidal deformation arising from the tides; in Section 3, the principal partial tides contributing to the dissipation will be detailed; Section 4 will be devoted to a determination of the extent of stellar viscosity — both gas and radiative; while in the concluding Section 5 quantitative estimates will be given of the actual rate at which the kinetic energy of dynamical tides gets dissipated into heat by viscous friction in stellar plasma.The results disclose that the amount of heat produced per unit time by tidal interaction between components of actual close binaries equals only about 10^–10th part of their nuclear energy production; and cannot, therefore, affect the internal structure of evolution of the constituent stars to any appreciable extent. Moreover, it is shown that the kinetic energy of their axial rotation can be influenced by tidal friction only on a nuclear, rather than gravitational (Kelvin) time-scale — as long as plasma or radiative viscosity constitute the sole sources of dissipation. However, the emergence of turbulent viscosity in secondary components of late spectral types, which have evolved away from the Main Sequence, can accelerate the dissipation 10⁵–10⁶ times, and thus give rise to appreciable changes in the elements of the system (particularly, in the orbital periods) over time intervals of the order of 10⁵–10⁶ years. Lastly, it is pointed out that, in close binary systems consisting of a pair of white dwarfs, a dissipation of the kinetic energy through viscous tides in degenerate fermion-gas could produce enough heat to account, by itself, for the observed luminosity of such objects.     

Tidal evolution in close binary systems,II

The aim of the present investigation will be to determine the explicit forms of differential equations which govern secular perturbations of the orbital elements of close binary systems in the plane of the orbit (i.e., of the semi-major axisA, eccentricitye, and longitude of the periastron ), arising from the lag of dynamical tides due to viscosity of stellar material. The results obtained are exact for any value of orbital eccentricity comprised between 0e<1; and include the effects produced by the second, third and fourth-harmonic dynamical tides, as well as by axial rotation with arbitrary inclination of the equator to the orbital plane.In Section 2 following brief introductory remarks the variational equations of the problem of plane motion will be set up in terms of the rectangular componentsR, S, W of disturbing accelerations with respect to a revolving system of coordinates. The explicit form of these coefficients will be established in Section 3 to the degree of accuracy to which squares and higher powers of quantities of the order of superficial distortion can be ignored. Section 4 will be devoted to a derivation of the explicit form of the variational equations for the case of a perturbing function arising from axial rotation; and in Section 5 we shall derive variational equations which govern the perturbation of orbital elements caused by lagging dynamical tides.Numerical integrations of these equations, which govern the tidal evolution of close binary systems prompted by viscous friction at constant mass, are being postponed for subsequent investigations.Prepared at the Lunar Science Institute, Houston, Texas, under the joint support of the Universities Space Research Association, Charlottesville, Virginia, and the National Aeronautics and Space Administration Manned Spacecraft Center, Houston, Texas, under Contract No. NSR 09-051-001. This paper constitutes Lunar Science Institute Contribution no. 100.Normally at the Department of Astronomy, University of Manchester, England.     

Dynamical tides in close binary systems,I

The aim of the present paper will be to develop from the fundamental equations of hydrodynamics a theory of dynamical tides in close binary systems, the components of which are regarded to consist of heterogeneous viscous fluid, and to revolve around their common centre of gravity in eccentric orbits; moreover, the equatorial planes of their axial rotation and the orbital plane need not be co-planar, but all may be inclined to the invariable plane of the system of arbitrary amounts. The changes in the pressure or density invoked by time-dependent deformation will be regarded as adiabatic; but, in the equilibrium state, both the density and viscosity of the material of our components may be arbitrary functions of the radial distance.Following a brief exposition in Section 2 of the fundamental equations linearized to small oscillations — be these free or forced — in Section 3 we shall particularize them to describe spheroidal deformations; with due regard to all terms arising from viscosity. Section 4 will contain a specification of the boundary conditions to be imposed upon such oscillations; and in Section 5 we shall solve the problem of non-radial oscillations of self-gravitating inviscid configurations in terms of hypergeometric series. The remaining Sections 6–8 will be devoted to a discussion of the phenomena arising from viscosity: in particular, we shall solve in a closed form the problem of non-radial oscillations of incompressible viscous globes in the terms of Bessel functions. It will be shown that the effect of viscosity — like those of compressibility — tend to de-stabilize all non-radial oscillations of homogeneous configurations.At the other extreme, a similar treatment of a mass-point model — as well as of one exhibiting high but finite degree of central condensation — is being postponed for a subsequent communication.     

Tidal evolution in close binary systems

The aim of the present paper will be to give a mathematical outline of the theory of tidal evolution in close binary systems of secularly constant total momentum — an evolution activated by viscous friction of dynamical tides raised by the two components on each other. The first section contains a general outline of the problem; and in Section 2 we shall establish the basic expressions for the energy and momenta of close binaries consisting of components of arbitrary internal structure. In Section 3 we shall investigate the maximum and minimum values of the energy (kinetic and potential) which such systems can attain for given amount of total momentum; while in Section 4 we shall compare these results with the actual facts encountered in binaries with components whose internal structure (and, therefore, rotational momenta) are known to us from evidence furnished by the observed rates of apsidal advance.The results show that all such systems — be these of detached or semi-detached type — disclose that more than 99% of their total momenta are stored in the orbital momentum. The sum of the rotational momenta of the constituent components amounts to less than a percent of the total — a situation characteristic of a state close to the minimum energy for given total momentum. This appears, moreover, to be true not only of the systems with both components on the Main Sequence, but also of those possessing evolved components in contact with their Roche limits.Under such conditions, a synchronism between rotation and revolution (characteristic of both extreme states of maximum and minimum energy) is not only possible, but appears to have been actually approached — if not attained — in the majority of cases. In other words, it would appear that — in at least a large majority of known cases — the existing close binaries have already attained orbits of maximum distension consistent with their momenta; and tidal evolution alone can no longer increase the present separations of the components to any appreciable extent.The virtual absence, in the sky, of binary systems intermediate between the stages of maximum and minimum energy for given momentum leads us to conjecture that the process of dynamical evolution activated by viscous tides may enroll on a time-scale which is relatively short in comparison with their total age — even for systems like Y Cygni or AG Persei, whose total age can scarcely exceed 10⁷ yr. A secular increase of the semi-major axes of relative orbits is dynamically coupled with a corresponding variation in the velocity of axial rotation of both components through the tidal lag arising from the viscosity of stellar material. The differential equations of so coupled a system are given in Section 5; but their solution still constitutes a task for the future.The Lunar Science Institute Contribution No. 90. The Lunar Science Institute is operated by the Universities Space Research Association under Contract No. NSR 09-051-001 with the National Aeronautics and Space Administration.     

The effects of viscous friction on the precession and nutation of celestial bodies

The aim of the present study has been to set the system of differential equations which govern the precession and nutation of self-gravitating globes of compressible viscous fluid, due to the attraction exerted on the rotating configuration by its companion; and to construct their approximate solution which are correct to terms of the second order in small dependent variables of the problem. Section 2 contains an explicit formulation of the effects of viscosity arising in this connection, given exactly as far as the viscosity remains a function of radial distancer only; but irrespective of its magnitude. In Section 3 the equations of motion will be linearized for the case of near-circular orbits and small inclinations andi of the equator of the rotating configuration, and of its orbital plane, to the invariable plane of the system; while in Section 4 further simplifications will be introduced which are legitimate for studies of secular (or long-periodic) motions of the nodes and inclinations. The actual solutions of so simplified a system of equations are constructed in Section 5; and these represent a generalization of the results obtained in our previous investigation (Kopal, 1969) of the inviscid case.The physical significance of the new results will be discussed in the concluding Section 6. It is demonstrated that the axes of rotation of deformable components in close binary systems are initially inclined to the orbital plane, viscous dissipation produced by dynamical tides will tend secularly to rectify their positions until perpendicularity to the orbital plane has been established, and the equators as well as orbit made to coincide with the invariable plane of the system-in a similar manner as other effects of tidal friction are bound eventually to synchronize the velocity of axial rotation with that of orbital revolution in the course of time.An application of the results of the present study to the dynamics of the Earth-Moon system discloses that the observed inclination of 1°.5 of the lunar equator to the ecliptic cannot be regarded as being secularly constant, but representing the present deviations from perpendicularity of oscillatory motion of very long period.The Lunar Science Institute is operated by the Universities Space Research Association under Contract No. NSR-09-051-001 with the National Aeronautics and Space Administration. This paper constitutes the Lunar Science Institute Contribution No. 85.     

Dissipation of energy through tidal deformation

The aim of the present paper will be to derive an equation of dissipation of energy for a rotating body of arbitrary viscosity distorted by tides, which arise from the gravitational field of its companion in a close pair of such bodies.By a transformation of the fundamental equation of energy dissipation in terms of velocity of tidal deformation (Section 2), the dissipation function is constructed for a tidally-distorted body (Section 3). From this equation, the rate of dissipation of tidal energy is formulated for a nearly-spherical rotating body distorted by second harmonic longitudinal tides (Section 4); the coefficients of viscosity (or the bulk modulus) are treated as arbitrary functions of spatial coordinates. Finally (Section 5), expressions for the total energy dissipation within the orbital cycle are given for axial rotation of the distorted body, provided its angular velocity is constant (for example, with the Keplerian angular velocity).Research financed in part by the Division of Scientific Research and Development of Ministry of Sciences and Culture of Greece.     

The precession and nutation of deformable bodies,III

In preceding papers of this series (Kopal, 1968; 1969) the Eulerian equations have been set up which govern the precession and nutation of self-gravitating fluid globes of arbitrary structures in inertial coordinates (space-axes) as well as with respect to the rotating body axes; with due account being taken of the effects arising from equilibrium as well as dynamical tides.In Section 1 of the present paper, the explicit form of these equations is recapitulated for subsequent solations. Section 2 contains then a detailed discussion of the coplanar case (in which the equation of the rotating configuration and the plane of its orbit coincide with the invariable plane of the system); and small fluctuations in the angular velocity of axial rotation arising from the tidal breathing in eccentric binary systems are investigated.In Section 3, we consider the angular velocity of rotation about theZ-axis to be constant, but allow for finite inclination of the equator to the orbital plane. The differential equations governing such a problem are set up exactly in terms of the time-dependent Eulerian angles and , and their coefficients averaged over a cycle. In Section 4, these equations are linearized by the assumption that the inclinations of the equator and the orbit to the invariable plane of the system are small enough for their squares to be negligible; and the equations of motion reduced to their canonical form.The solution of these equations — giving the periods of precession and nutation of rotating components of close binary systems, as well as the rate of nodal regression which is synchronised with precession — are expressed in terms of the physical properties of the respective system and of its constituent components; while the concluding Section 6 contains a discussion of the results, in which the differences between the precession and nutation of rigid and fluid bodies are pointed out.     

Distribution of brightness over apparent discs of distorted stars

The aim of the present paper will be to establish the explicit form of the equations of radiative transfer, in plane-parallel atmospheres surrounding the stars which are distorted by axial rotation or tides, in curvilinear coordinates which parallel the distorted surface; with particular attention to the circumstances under which the effects arising from limb- and gravity-darkening are multiplicative and admit of algebraic separation. In Section 2 (which follows a general outline of our problem) the fundamental equations of the radiativetransfer problem will be formulated for the ‘grey’ case; and rewritten in Section 3 in terms of non-orthogonal coordinates in which the potential over a level surface in hydrostatic equilibrium replaces the radial coordinate of spherical polars. In Section 4 we shall proceed to construct an explicit solution of the corresponding transfer problem in a plane-parallel approximation; and to prove that the effects of limb- and gravity-darkening remain factorizable only to terms which are linear in the cosines μ of the angle of foreshortening. Lastly, in Section 5 we shall list additional problems, arising in this connection, which still await appropriate treatment.     

The weak friction approximation and tidal evolution in close binary systems

The aim of this paper is to study the dynamical problem of tidal friction in a binary system consisting of deformable components, with the restriction that the angle of lag or advance of the tidal distortion with respect to the direction of the disturbing companion is small. The fractional distortion of the bodies due to rotation and tidal interaction is also treated as a first-order small quantity, and terms up to the fourth harmonic in the tidal potential are retained. In this linear approximation, the time-dependent tidal potential can be Fourier decomposed into a spectrum of simple harmonic terms, each of which is responsible for raising a partial wave in the body; each such partial wave can then be treated independently of the others. This is the method first employed by Darwin.In Section 2, it is assumed that the phase lag in the response of the body (due to dissipation of kinetic energy of deformation) is proportional to the forcing frequency, which is justified for small amplitude oscillations of a viscous fluid or visco-elastic body. A simple expression is then obtained for the potential function for the distortion in terms of the disturbing potential and the structure of the body.In Section 3, the distortion potential function is employed in deriving the componentsR, S andW of the disturbing force which are then substituted in the Gaussian form of the equations for variation of the elements. In Section 4, the Eulerian equations for motion of deformable bodies are derived, using the so-called mean axes of the body as the rotating axes of reference. In Section 5, it is shown that the dynamical effects of rotational distortion occur on a much shorter time scale than those arising from tidal friction, which allows one to consider the two phenomena as acting independently of one another. The collected set of Gaussian (orbital) and Eulerian (body) equations is re-written in terms of dimensionless variables for the tidal friction case, and the stability of the system is examined on the basis of these equations.In Section 6, the tidal friction equations are integrated numerically for the close binary system AG Persei and for the Earth-Moon system. In the former, the integrations were started from a highly elliptical orbit and the system was found to relax into a circular orbit, with synchronous rotation perpendicular to the orbit. In the latter, the integrations were performed backwards in time from the present day, and it was found that the lunar orbit rapidly becomes highly elliptical at the time of closest approach, thus indicating a probable capture of the Moon by the Earth. This result is in agreement with that obtained by other investigators; however, it is shown that the detailed behaviour of the system at the time of capture, in particular the inclination of the lunar orbit to the ecliptic, depends critically on the chosen rate of dissipation in the Moon's interior. A simple argument is presented which allows an estimation for the mean viscosity of a fluid body from the known age of the system: for the components of AG per, the result is 2×10¹¹ g cm^–1 s^–1, indicating that the stars must have possessed turbulent convective outer regions during some part of their tidal evolution, while for the Earth, the result, is 1.4×10¹² g cm^–1 s^–1. It is shown that the angle of tidal lag in nonsynchronous close binary systems in general is expected to be extremely small, and not observationally detectable.     

月球卫星轨道变化的分析解

由于月球自转缓慢及其引力位的特点,使得讨论月球卫星与人造地球卫星轨道变化的方法有所不同。     

Vibrational stability of the components of close binary systems

The aim of the present paper will be to extend our previous investigation of the vibrational stability of rotating configurations (Kopal, 1981) to a similar investigation of the stability of the components of close binary systems which not only rotate, but also distort each other by tidal action. To this end, differential equations which govern first-order oscillations of arbitrary spherical-harmonic symmetry will be set up in Clairaut coordinates in which the radial coordinate is replaced by the potential which remains constant over level surfaces of equilibrium configurations; introduced by us in an earlier paper (Kopal, 1980), and their form detailed for surface distorted by second-, third-, and fourth-harmonic tides raised by the external mass; and their boundary conditions established. A solution of such differential boundary-value problems arising in connection with the stars of arbitrary structure remains, of course, a task for automatic computers. It may only be added that the tide-generating potential Ψ _T established in this paper should enable us to study, by the same method, not only free, but also forced oscillations of the components of close binary systems, arising from orbital eccentricity of the respective couples, dynamical tides, or other causes likely to be operative in such systems.     

The Roche Coordinates in three dimensions and their application to problems of double-star astronomy

In a preceding paper (Kopal, 1969; in what follows referred to as Paper I) we introduced a new system of curvilinear coordinates-hereafter referred to as Roche Coordinates — in which spheres of constant radius in spherical polars have been replaced by surfaces of constant potential of a rotating gravitational dipole; while the angular coordinates are orthogonal to the equipotentials. In Paper I we established an explicit form of such a transformation, and related the Roche coordinates with polar coordinates (with which they coalesce in the immediate neighbourhood of each one of the two finite mass-points) in the plane case. The aim of the present investigation will be to generalize the definition of the Roche coordinates to three dimensions.The opening Section 1 of this paper will contain a general outline of the proposed three-dimensional transformation; and in Section 2 details of this transformation will be explicitly worked out correctly to quantities of first order in superficial distortion — an approximation which should prove adequate in regions surrounding the two finite masses; while in Section 3 we shall evaluate (to this degree of accuracy) the metric coefficients of the respective transformation, and its direction cosines, in both polar and curvilinear coordinates. Section 4 will then contain a formulation of the fundamental equations of hydrodynamics in terms of the three-dimensional Roche coordinates; and their advantages for a treatment of certain classes of dynamical problems encountered in doublestar astronomy will be illustrated in the concluding Section 5 by an investigation of the vibrational stability of the Roche model. We shall show that this model is capable of performing free radial oscillations which remain barotropic only if its equilibrium form is spherical (i.e., in the absence of any external mass in the neighbourhood); but not if it is distorted to any extent by rotation or tides.     

Tidal evolution in close binary systems,III

In this paper we shall investigate the energy of close binary systems of constant momentum takng into consideration the first-order effects of rotation and tidal attraction of the components of finite size. The equations for the momentum and the energy of the system will be set up in Section 2, making use of terms including the effects of finite size of the components of finite degree of central condensation. In Section 3 perturbation theory is applied to these equations using the results of Kopal (1972b) as our initial values. In Section 4 we shall compare our results with the initial values and then discuss variations in our constants and the application to various real systems.     

An Estimate of the Long-term Tendency of Tidal Evolution of Earth-lunar System

According to the conservation principle of angular momentum, we calculate in this paper the revolution period and the distance between the Earth and the Moon in the equilibrium state of the tidal evolution in the Earth-Moon system. The difference of energy between the current state and the equilibrium state is used to compute the time needed to fulfil the equilibrium state. Then the long-term variations of the Earth-Moon distance and of the Earth rotation rate are further estimated.     

An Estimate on the Lunar Figure

In 1799 Laplace discovered that the three principal moments of the Moon are not in equilibrium with the Moon's current orbital and rotational state. Some authors suggested that the Moon may carry a fossil figure. More than 3 billion years ago, the liquid Moon was closer to the Earth and revolved faster. Then the Moon migrated outwards and its rotation slowed down. During the early stage of this migration, the Moon was continually subjected to tidal and rotational stretching and formed into an ellipsoid. Subsequently the Moon cooled down and solidified quickly. Eventually, the solid Moon's lithosphere was stable and as a result we may see the very early lunar figure.     

Effects of rotation on internal structure of the stars

The aim of the present paper will be to establish the explicit form of the equations which govern the internal structure of stars rotating with constant angular velocity formulated in terms of Clairaut coordinates (cf. Kopal, 1980) in which the radial coordinate is replaced by the total potential, which for equilibrium configurations remains constant over distorted level surfaces. The introductory Section 1 contains an account of previous work on rotating stars, commencing with Milne (1923), von Zeipel (1924) and Chandrasekhar (1933), who all employed orthogonal coordinates for their analysis. In Section 2 we shall apply to this end the curvilinear Clairaut coordinates introduced already in our previous work (cf. Kopal, 1980, 1981); and although these are not orthogonal, this disadvantage is more than offset by the fact that, in their terms, the fundamental equation of our problem will assume the form of ordinary differential equations, subject to very simple boundary conditions. The explicit form of these equations — exact to terms of fourth order in surficial distortion caused by centrifugal force—will be obtained in Section 3; while in the concluding Section 4 these will be particularized (for the sake of comparison with work of previous investigators) to stars of initially polytropic structure. These will prove to be much simpler in Clairaut coordinates than they were in any previously used frame of reference. Lastly, in Appendix A we shall present the explicit forms, in Clairaut coordinates, of the differential operators which were needed to establish the results given in Sections 3–4; while Appendix B will summarize other auxiliary algebraic relations of which use was made to formulate our fourth-order theory developed in Section 3.     

Differential Rotation Driven by Precession

The tidal force in the Earth–Moon system exerted on the Earth's equatorial bulge results in the Earth's precession. It was proposed a long time ago that the strong shear flow driven by the precession of the Earth may power the Earth's dynamo in its liquid core. We present a nonlinear analytical study investigating how the Poincaré force in a rotating, precessing spherical system drives a large-amplitude differential rotation which plays a major role in the modern theory of the geodynamo. The analysis is based on a perturbation approach in terms of the small Poincaré force parameter. It is found that the amplitude of the precession-driven differential rotation is consistent with that estimated from the geomagnetic secular variation.     

On secular stability of tidally distorted stars of arbitrary structure

The aim of the present paper will be to develop a theory which should make it possible to investigate secular stability of close binary systems, consisting of tidally-distorted components of arbitrary internal structure, by a minimization of the potential energy of the system as a whole. In the second section which follows brief introductory remarks, appropriate expressions for the total potential energy of a close binary will be formulated. Section 3 will be concerned mainly with the nature of the tide-generating potential, and its effects on the shape of each star. In Section 4, the amplitudes of partial tides raised by this potential will be specified, for stars of arbitrary structure, correctly to terms of second order in superficial distortion; and in Section 5 we shall investigate the effects of interaction between rotation and tides to the same degree of approximation. The concluding Section 6 will then contain an explicit formulation of different constituents adding up to the total potential energy of the system, which can be used as a basis for its secular stability by the methods outlined already in our previous investigation (Kopal, 1973).     

Rotational dynamics of celestial deformable bodies

In a previous paper of this series (Tokis, 1974b), we have discussed the solution of the Eulerian equation which governs the axial rotation, applied to the effects of viscous friction exhibited in binary systems which consist of a close pair of fluid bodies of arbitrary structure. The aim of the present paper will be to give an application of those results to the Earth-Moon system. It is shown that synchronism between the axial rotation of the Earth and the revolution of the Moon will occur at the value of 650 h, in a time scale which depends strongly on the value of the mean viscosity of the Earth (regarded as spherical or spheroidal). In particular, the variation of rotational angular velocity of the Earth over the next ten centuries commencing from 1900 A.D., depends sensitively on the value of viscosity. On the other hand, the time for synchronism of axial rotation of the Moon is not affected by the viscosity for values between 10²⁴g cm^?1 s^?1 and 10²⁷g cm^?1 s^?1.