On a new form for the differential equations of relative motion of the three-body problem

A new method for the development of the disturbing function of the three-body problem is outlined in this paper. A special process is devised to get the distance between two planetsP ₁ andP ₂ in terms of their heliocentric distances. It is then shown that the differential equations of relative motion of this problem can be brought in an homogeneous set of differential equations.     

Analytical solution to the perturbedJ 2 motion of Artificial Satellite in terms of Euler parameters

In this paper an analytical solution for the differential equations governing the motion of Artificial Satellite under the oblateness of the earth will be developed. We start with the differential equations in terms of Euler parameters. To compact algebra, we introduce the Cayley-Klein type complex variables. Comparison between numerical and analytical final states withJ ₂ will be given for a test case.     

Poisson equations of rotational motion for a rigid triaxial body with application to a tumbling artificial satellite

Differential equations are derived for studying the effects of either conservative or nonconservative torques on the attitude motion of a tumbling triaxial rigid satellite. These equations, which are analogous to the Lagrange planetary equations for osculating elements, are then used to study the attitude motions of a rapidly spinning, triaxial, rigid satellite about its center of mass, which, in turn, is constrained to move in an elliptic orbit about an attracting point mass. The only torques considered are the gravity-gradient torques associated with an inverse-square field. The effects of oblateness of the central body on the orbit are included, in that, the apsidal line of the orbit is permitted to rotate at a constant rate while the orbital plane is permitted to precess (either posigrade or retrograde) at a constant rate with constant inclination.A method of averaging is used to obtain an intermediate set of averaged differential equations for the nonresonant, secular behavior of the osculating elements which describe the complete rotational motions of the body about its center of mass. The averaged differential equations are then integrated to obtain long-term secular solutions for the osculating elements. These solutions may be used to predict both the orientation of the body with respect to a nonrotating coordinate system and the motion of the rotational angular momentum about the center of mass. The complete development is valid to first order in (n/w ₀)², wheren is the satellite's orbital mean motion andw ₀ its initial rotational angular speed.     

Motion of artificial satellites in the set of Eulerian redundant parameters (III)

In this paper, the classical and generalized Sundman time transformations are used to establish new generating set of differential equations of motion in terms of the Eulerian redundant parameters. The implementation of this set on digital computers for the commonly used independent variables is developed once and for all. Motion prediction algorithms based on these equations are developed in a recursive manner for the motions in the Earth's gravitational field with axial symmetry whatever the number of the zonal harmonic terms may be. Applications for the two types of short and long term predictions are considered for the perturbed motion in the Earth's gravitational field with axial symmetry with zonal harmonic terms up to J ₃₆ . Numerical results proved the very high efficiency and flexibility of the developed equations.     

Notes on Hori Method for Canonical Systems

In this paper a slightly different approach is proposed for the process of determining the functions S _m and H _m ^* of the algorithm of the canonical version of Hori method. This process will be referred to as integration theory of the mth order equation of the method. It will be shown that the ordinary differential equation with an auxiliary parameter t ^* as independent variable, introduced through Hori auxiliary system, can be replaced by a partial differential equation in the time t. In this way, the mth order equation of the algorithm assumes a form very similar to the one of other perturbation methods. In virtue of this new approach of the integration theory for Hori method, Lagrange's variational equations introduced by Sessin are revised. As an example, the Duffing equation is solved through this new approach.     

Equilibrium Structures of Differentially Rotating Primary Components of Binary Stars

In this paper a method is proposed for computing the equilibrium structures and various other observable physical parameters of the primary components of stars in binary systems assuming that the primary is more massive than the secondary and is rotating differentially about its axis. Kippenhahn and Thomas averaging approach (1970) is used in a manner earlier used by Mohan, Saxena and Agarwal (1990) to incorporate the rotational and tidal effects in the equations of stellar structure. Explicit expressions for the distortional terms appearing in the stellar structure equations have been obtained by assuming a general law of differential rotation of the typeω² = b ₀+b ₁ s ²+b ₂ s ⁴, where ω is the angular velocity of rotation of a fluid element in the star at a distance s from the axis of rotation, and b ₀, b ₁, b ₂ are suitably chosen numerical constants. The expressions incorporate the effects of differential rotation and tidal distortions up to second order terms. The use of the proposed method has been illustrated by applying it to obtain the structures and observable parameters of certain differentially rotating primary components of the binary stars assuming the primary components to have polytropic structures. This revised version was published online in July 2006 with corrections to the Cover Date.     

Viscous ΛCDM universe models

We explore flat ΛCDM models with bulk viscosity, and study the role of the bulk viscosity in the evolution of these universe models. The dynamical equations for these models are obtained and solved for some cases of bulk viscosity. We obtain differential equations for the Hubble parameter H and the energy density of dark matter ρ _m , for which we give analytical solutions for some cases and for the general case we give a numerical solution. Also we calculate the statefinder parameters for these models and display them in the s–r-plane.     

On the Einstein-Maxwell field equations

Starting with the Einstein-Maxwell field equations in general relativity we construct the general differential equations governing the components of the metric tensor. We do this in a fashion completely analogous to that which we follow in Einstein-Infeld-Hoffmann and Einstein-Infeld theory. These equations allow us to findh _if in various orders. An answer to our problem up to the first relativistic corrections is a computational work to find₄ h ₀₀, since the other terms₂ h ₀₀,₂ h and₃ h ₀ are as in a pure gravitational case. On the other hand, using the defined Einstein-Maxwell tensor, we give the equations of motion of two charged particles in the 0th order; also, the generalization is given in the case ofn particles.     

Equilibrium structure of differentially rotating polytropic models of the stars

In this paper we propose a method for computing the equilibrium structure of differentially rotating polytropic models of the stars. A general law of differential rotation of the type ²=b ₀+b ₁ s ²+b ₂ s ⁴, which can account for a reasonably large variety of possible differential rotations in the stars has been used. The distortional effects have been incorporated in the structure equations up to second order of smallness in distortion parametersb ₀,b ₁, andb ₂ using Kippenhahn and Thomas' averaging approach in conjunction with Kopal's results on Roche equipotentials in manner similar to the one earlier used by Mohan and Saxena for computing the equilibrium structure of polytropes having solid body rotation. Numerical results have been obtained for various types of differentially rotating polytropic models of stars of polytropic indices 1.5, 3, and 4. Certain differentially rotating models of the Sun which are possible with such a type of law of differential rotation, have also been computed.     

The color-redshift relation for giant elliptical galaxies

The colors of giant elliptical (gE) galaxies in clusters out to redshiftz=0.2, observed by Oke and Sandage (1968), are studied for systematic color-redshift effects. To reveal any intrinsic changes, theK term is subtracted from each color, after correction for galactic reddening and the resulting (B-V) _c -K _B-V versus z relation analysed. If the Oke and SandageK terms (relevant to nuclear colors) are used, the best fitting linear relation shows negligible change withz. But if the Whitford (1970)K terms (relevant to integrated colors) are used, there is a trend to bluer colors, by 0.07 mag.atz=0.2 if higher weight is given the better observed clusters. An upper limit, of ±0.08 mag. atz=0.2, is set to possible systematic aperture effects, by the total change between nuclear and integrated B-V and U-B of nearbygE galaxies.The color-redshift trends, interpreted as evolutionary changes, are related to evolution in the magnitude-redshift relation by means of models of stellar evolution in agE galaxy. If one uses the linear fit to the color-redshift relation obtained with WhitfordK terms and neglect of any aperture effects (which is the appropriate case if Oke and Sandage used large enough apertures), and if one adopts plausible limits to the ratio between color and magnitude evolution, the effect of evolution in the magnitude-redshift relation results in a negative value ofq ₀. There are still substantial observational and theoretical uncertainties affecting this conclusion.     

Fourier analysis of the light curves of eclipsing variables,XX

In the present paper _n ⁰ , for occulation and transit eclipses of partial phases, are evaluated numerically by means of the Runge-Kutta methods. Section 2 contains the required differential equations of _n ⁰ with respect to the modulusX orC, and Section 3 includes the numerical method of the solutions of these differential equations. Theoretical values of ₀ ⁰ and ₁ ⁰ , with corresponding values ofC, are also added in this section.     

Vibrational stability of rotating stars

The aim of the present paper will be to detail the explicit form of the equations which govern first-order oscillations of fast-rotating globes of self-gravitating fluids; with due account taken of the effects arising from the centrifugal as well as Coriolis force. As such configurations oscillate in general about distorted figures of equilibrium, the equations governing them can be conveniently expressed in terms of the Clairaut coordinates, associated with distorted spheroidal figures, and introduced in our previous paper (Kopal, 1980) for this purpose.In Section 2 which follows a brief outline of our problem, the equilibrium properties of fast-rotating configurations or arbitrary structure will be formulated. In Section 3 we shall carry out a separation of the variables in the equations of motion, and reduce the partial differential equations of the problem to an equivalent system of ordinary differential equations, by an expansion of expressions for the velocity componentsU, V, W in terms of tesseral harmonicsY _n ^m (, ). The explicit form of such a system, including the effects of all tesseral harmonics of orders up tom=n=4, will be specified in Section 3 for configurations whose equilibrium form is a sphere; while in Section 4 this latter condition will be relaxed to allow for the equilibrium configuration to become a rotational spheroid.In the concluding Section 5 we shall convert the complex form of our equations of motion into real terms, amenable to a solution-analytical or numerical-in terms of real variables; and shall establish the boundary conditions necessary for a specification of the characteristic frequencies of oscillation.     

Special Families of Orbits in the Direct Problem of Dynamics

The direct problem of dynamics in two dimensions is modeled by a nonlinear second-order partial differential equation, which is therefore difficult to be solved. The task may be made easier by adding some constraints on the unknown function = f _y /f _x , where f(x, y) = c is the monoparametric family of orbits traced in the xy Cartesian plane by a material point of unit mass, under the action of a given potential V(x, y). If the function is supposed to verify a linear first-order partial differential equation, for potentials V satisfying a differential condition, can be found as a common solution of certain polynomial equations.The various situations which can appear are discussed and are then illustrated by some examples, for which the energy on the members of the family, as well as the region where the motion takes place, are determined. One example is dedicated to a Hénon—Heiles type potential, while another one gives rise to families of isothermal curves (a special case of orthogonal families). The connection between the inverse/direct problem of dynamics and the possibility of detecting integrability of a given potential is briefly discussed.This revised version was published online in October 2005 with corrections to the Cover Date.     

A note on the solution of the saturation flux limited solar wind equations

The solution curves of the differential equations determining the behavior of the solar wind are calculated for the case where the heat flux has its maximum value 3/2 nkTv _th. All the supersonic solutions are asymptotically adiabatic, T r ^-4/3.The National Center for Atmospheric Research is sponsored by the National Science Foundation.     

Anharmonicity effects in bounce motion of particles

Equations for parallel motion for a particle trapped in a magnetic field have been considered and improved solutions of differential equations have been derived. The expressions for the change in energy of the particle (Δw) and diffusion coefficient (D_ww) have been presented in a simple form using the improved solution.     

Theory of the motion of an artificial Earth satellite

An improved analytical solution is obtained for the motion of an artificial Earth satellite under the combined influences of gravity and atmospheric drag. The gravitational model includes zonal harmonics throughJ ₄, and the atmospheric model assumes a nonrotating spherical power density function. The differential equations are developed through second order under the assumption that the second zonal harmonic and the drag coefficient are both first-order terms, while the remaining zonal harmonics are of second order.Canonical transformations and the method of averaging are used to obtain transformations of variables which significantly simplify the transformed differential equations. A solution for these transformed equations is found; and this solution, in conjunction with the transformations cited above, gives equations for computing the six osculating orbital elements which describe the orbital motion of the satellite. The solution is valid for all eccentricities greater than 0 and less than 0.1 and all inclinations not near 0^o or the critical inclination. Approximately ninety percent of the satellites currently in orbit satisfy all these restrictions.     

Hydrodynamic equations of differentially rotating stellar systems and the velocity ellipsoid

Anisotropic hydrodynamic equations for differentially rotating collisionless stellar systems are derived. These equations can describe the evolution of the systems in a time span longer than their rotation periods.As a by-product of derivation of hydrodynamic equations, the well-known relation that the ratio of the principal axes of the velocity ellipse in a differentially rotating stellar disk is [B/(B-A)]^1/2 is re-found if the system is in a purely circular rotation, whereA andB are the Oort's constants. In addition, we find a systematic mean motion superposed on a purely circular differential rotation makes the directions of axes of the velocity ellipse deviate from the radial and the transverse direction. The observed deviation of directions of axes in our neighbourhood in the Galaxy can be explained if in the mean motion superposed on a purely circular differential rotatin the gas of stars near us is compressed in the radial direction or rarefied in the transverse directions, with irregularities of the order of 5 km/sec in amplitude of velocity and 1 kpc in size. These magnitudes of irregularities agree with those actually observed or with those anticipated from other theoretical considerations.     

Self‐similar evolutionary solutions for an accreting magneto‐fluid around a compact object with finite electrical conductivity

In this paper, we investigate the time evolution of an accreting magneto‐fluid with finite conductivity. For the case of a thin disk, the fluid equations along with Maxwell's equations are derived in a simplified, one‐dimensional model that neglects the latitudinal dependence of the flow. The finite electrical conductivity of the plasma is taken into account by Ohm's law; however, the shear viscous stress is neglected, as well as the self‐gravity of the disk. In order to solve the integrated equations that govern the dynamical behaviour of the magneto‐fluid, we have used a self‐similar solution. We introduce two dimensionless variables, S₀ and ε_ϱ, which represent the size of the electrical conductivity and the density behaviour with time, respectively. The effect of each of these on the structure of the disk is studied. While the pressure is obtained simply by solving an ordinary differential equation, the density, the magnetic field, the radial velocity, and the rotational velocity are presented analytically. The solutions show that the S0 and ε_ϱ parameters affect the radial thickness of the disk. Also, radial velocity and gas pressure are more sensitive to the electrical conductivity in the inner regions of disk. Moreover, the parameter ε_ϱ has a more significant effect on the physical quantities for small radii. (© 2015 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Stabilized Kepler motion connected with analytic step adaptation

In the publication Baumgarte and Stiefel (1974a) a stabilization of the Keplerian motion was offered by making use of a manipulation of the Hamiltonian. By this stabilization technique the given HamiltonianH(p _i,q _i) is replaced by a new HamiltonianH ^*, which leads to Lyapunov-stable differential equations of motion.Whereas, in the quoted publication, the physical timet was used as the independent variable we now develop a generalization which allows to combine the stabilization with the introduction of a new independent variables. Such a fictitious times is popular for achieving an analytic step-size adaptation (Baumgarte and Stiefel, 1974c). Perturbations of Kepler motion are discussed.     

A class of semi-analytical solutions for a rotating toroidal plasma

In this paper, the problem of stationary MHD flow for a rotating toroidal plasma is investigated by assuming that the entropy is a surface quantity. Then, the system of ideal MHD equations is reduced to a single second-order elliptic partial differential equation known as the modified Grad-Shafranov (or Maschke-Perrin) equation. Under the assumption that both the function,P _s andf ² are quadratic polynomials of the flux function, a class of semi-analytical solutions is obtained for a plasma contained in a perfectly conducting toroidal boundary with a rectangular cross section. The flux function, poloidal current and the generalized pressure are obtained and discussed for relevant values of the parameters.