Exact analytic solution for the rotation of a rigid body having spherical ellipsoid of inertia and subjected to a constant torque

The exact analytic solution is introduced for the rotational motion of a rigid body having three equal principal moments of inertia and subjected to an external torque vector which is constant for an observer fixed with the body, and to arbitrary initial angular velocity. In the paper a parametrization of the rotation by three complex numbers is used. In particular, the rows of the rotation matrix are seen as elements of the unit sphere and projected, by stereographic projection, onto points on the complex plane. In this representation, the kinematic differential equation reduces to an equation of Riccati type, which is solved through appropriate choices of substitutions, thereby yielding an analytic solution in terms of confluent hypergeometric functions. The rotation matrix is recovered from the three complex rotation variables by inverse stereographic map. The results of a numerical experiment confirming the exactness of the analytic solution are reported. The newly found analytic solution is valid for any motion time length and rotation amplitude. The present paper adds a further element to the small set of special cases for which an exact solution of the rotational motion of a rigid body exists.     

A fourth-order solution of the ideal resonance problem

The second-order solution of the Ideal Resonance Problem, obtained by Henrard and Wauthier (1988), is developed further to fourth order applying the same method. The solutions for the critical argument and the momentum are expressed in terms of elementary functions depending on the time variable of the pendulum as independent variable. This variable is related to the original time variable through a Kepler-equation. An explicit solution is given for this equation in terms of elliptic integrals and functions. The fourth-order formal solution is compared with numerical solutions obtained from direct numerical integrations of the equations of motion for two specific Hamiltonians.     

Analysis of solar radiation pressure induced coupled liberations of a gravity stabilized axisymmetric satellite

A planar, fixed-orbit model of the rotation of the planet Mercury is analyzed. The model includes only the solar torques on the planet's permanent asymmetry and its solar tidal bulge. For this model, it is shown that the zero of the averaged tidal torque corresponds to an asymptotically stable periodic solution of the second kind which, for two tidal torque representations, is close to the asymptotically stable equilibrium point corresponding to an exact 32 spin-orbit resonance. A conjecture that the current rotation state of Mercury is due to transfer from capture by the zero of the averaged tidal torque to 32 resonance capture with changes in the eccentricity of the planet's orbit is discussed briefly.     

Short-axis-mode rotation of a free rigid body by perturbation series

A simple rearrangement of the torque free motion Hamiltonian shapes it as a perturbation problem for bodies rotating close to the principal axis of maximum inertia, independently of their triaxiality. The complete reduction of the main part of this Hamiltonian via the Hamilton–Jacobi equation provides the action-angle variables that ease the construction of a perturbation solution by Lie transforms. The lowest orders of the transformation equations of the perturbation solution are checked to agree with Kinoshita’s corresponding expansions for the exact solution of the free rigid body problem. For approximately axisymmetric bodies rotating close to the principal axis of maximum inertia, the common case of major solar system bodies, the new approach is advantageous over classical expansions based on a small triaxiality parameter.     

Comparison of the classical and the global solutions of the Ideal Resonance Problem

The Ideal Resonance Problem is defined by the Hamiltonian $$F = B(y) + 2\varepsilon A(y) \sin ^2 x,\varepsilon \ll 1.$$ The classical solution of the Problem, expanded in powers of ε, carries the derivativeB′ as a divisor and is, therefore, singular at the zero ofB′, associated with resonance. With α denoting theresonance parameter, defined by $$\alpha \equiv - B'/|4AB''|^{1/2} \mu ,\mu = \varepsilon ^{1/2} ,$$ it is shown here that the classical solution is valid only for $$\alpha ^2 \geqslant 0(1/\mu ).$$ In contrast, the global solution (Garfinkelet al., 1971), expanded in powers ofμ=ε^1/2, removes the classical singularity atB′=0, and is valid for all α. It is also shown here that the classical solution is an asymptotic approximation, for largeα ², of the global solution expanded in powers ofα ^?2. This result leads to simplified expressions for resonancewidth and resonantamplification. The two solutions are compared with regard to their general behavior and their accuracy. It is noted that the global solution represents a perturbed simple pendulum, while the classical solution is the limiting case of a pendulum in a state offast circulation.     

The origin of scaling in the galaxy distribution

I present an analytic model for non‐linear clustering of the luminous (baryonic) material in a universe in which the gravitational field is dominated by dark matter. The model is based on a two-component generalization of the adhesion approximation in which the gravitational potential of the dark component is determined by the standard Zel'dovich approximation or one of its variants, or by an N ‐body simulation. The baryonic matter flow is dissipative and is driven by this dark matter gravitational potential. The velocity potential of the matter is described by a generalization of the Burgers equation: the random heat equation ('RH equation') with a spatially correlated Gaussian driving potential.
The properties of the RH equation are well understood: it is closely related to the equation for the Anderson model and to Brownian motion in a random potential: the solution can be expressed in terms of path integrals. Using this it is possible to derive the scaling properties of the solution and, in particular, those of the resultant velocity field. Even though the flow is non‐linear, the velocity field remains Gaussian and inherits its scaling properties from the gravitational potential. This provides an underlying dynamical reason why the density field in the baryonic component is lognormally distributed and manifests multifractal scaling.
By explicitly putting dark and luminous matter on different footings, the model provides an improved framework for considering the growth of large‐scale cosmic structure. It provides a solution for the velocity potential of the baryonic component in closed form (albeit a path integral) from which the statistical properties of the baryonic flow can be derived.     

A generic dynamical model of gamma-ray burst remnants

The conventional generic model is considered to explain the dynamics of gamma-ray burst remnants very well, no matter whether they are adiabatic or highly radiative. However, we find that, for adiabatic expansion, the model cannot reproduce the Sedov solution in the non-relativistic phase, and thus it needs to be revised. In this paper a new differential equation is derived. The generic model based on this equation is shown to be correct for both radiative and adiabatic fireballs, and in both ultrarelativistic and non-relativistic phases.     

Rotational evolution of solar-type stars

The aim of the present investigation has been to consider rotational evolution of solar-type stars simulated by a polytropic model that possesses differential rotation of Clement's type. A properly defined reduction factor moderates the effects of such a rotation. The present treatment is based upon the general Eulerian equation, governing nonuniform (i.e., nonrigid-body) rotation, which has been set up in a previous investigation. Nonconservative terms, arising when stellar wind torque is under consideration, are taken into account. Data available for the viscosity of the Sun are used to construct a plausible viscosity model. Certain assumptions are made that remove the mathematical difficulties and simplify the physical ground. The obtained results are compared to corresponding estimates of recent observations.     

Reduction the secular solution to periodic solution in the generalized restricted three-body problem

The aim of the present work is to find the secular solution around the triangular equilibrium points and reduce it to the periodic solution in the frame work of the generalized restricted thee-body problem. This model is generalized in sense that both the primaries are oblate and radiating as well as the gravitational potential from a belt. We show that the linearized equation of motion of the infinitesimal body around the triangular equilibrium points has a secular solution when the value of mass ratio equals the critical mass value. Moreover, we reduce this solution to periodic solution, as well as some numerical and graphical investigations for the effects of the perturbed forces are introduced. This model can be used to examine the existence of a dust particle near the triangular points of an oblate and radiating binary stars system surrounded by a belt.     

First-order theory of satellite attitude motion application to HIPPARCOS

This work aims at finding an analytic solution corresponding to the attitude evolution in space of a satellite submitted to disturbing torques. This paper presents a basic frame applicable to any perturbed rotation satellite, and a method of resolution leading to a formal solution which is given here to the first order. Thus, the main problem is the slow rotation of a body with three unequal axes of inertia, essentially submitted to a dominant solar radiation pressure torque, with the axis pointing far away from a position of equilibrium. The comparison of the results with a numerical integration based upon a HIPPARCOS model is convincing.     

Magnetized, mass-loaded, rotating accretion flows

We present a semi-analytical investigation of a simple one-dimensional, steady-state model for a mass-loaded, rotating, magnetized, hydrodynamical flow. Our approach is analogous to one used in early studies of magnetized winds. The model represents the infall towards a central point mass of the gas generated in a cluster of stars surrounding it, as is likely to occur in some active nuclei and starburst galaxies. We describe the properties of the different classes of infall solutions. We find that the flow becomes faster than the fast-mode speed, and hence decoupled from the centre, only for a limited range of parameter values, and when magnetic stresses are ineffective. Such flow is slowed as it approaches a centrifugal barrier, implying the existence of an accretion disc. When the flow does not become super-fast and the magnetic torque is insufficient, no steady solution extending inward to the centre exists. Finally, with a larger magnetic torque, solutions representing steady sub-Alfvénic flows are found, which can resemble spherical hydrodynamical infall. Such solutions, if applicable, would imply that rotation is not important and that any accretion disc formed would be of very limited size.     

Out-of-plane librations of spinning tethered satellite systems

We analyze the out-of-plane librations of a tethered satellite system that is nominally rotating in the orbit plane. To isolate the librational dynamics, the system is modeled as two point masses connected by a rigid rod with the system mass center constrained to an unperturbed circular orbit. For small out-of-plane librations, the in-plane motion is unaffected by the out-of-plane librations and a solution for the in-plane motion is determined in terms of Jacobi elliptic functions. This solution is used in the linearized equation for the out-of-plane librations, resulting in a Hill’s equation. Floquet theory is used to analyze the Hill’s equation, and we show that the out-of-plane librations are unstable for certain ranges of in-plane spin rate. For relatively high in-plane spin rates, the out-of-plane librations are stable, and the Hill’s equation can be approximated by a Mathieu’s equation. Approximate solutions to the Mathieu’s equation are determined, and we analyze the dominant characteristics of the out-of-plane librations for high in-plane spin rates. The results obtained from the analysis of the linearized equations of motion are compared to numerical simulations of the nonlinear equations of motion, as well as numerical simulations of a more realistic system model that accounts for tether flexibility. The instabilities discovered from the linear analysis are present in both the nonlinear system and the more realistic system model. The approximate solutions for the out-of-plane librations compare well to the nonlinear system for relatively high in-plane rotation rates, and also capture the significant qualitative behavior of the flexible system.     

A viable dark fluid model

We consider a cosmological model based on a generalization of the equation of state proposed by Nojiri and Odintsov (2004) and ?tefan?i? (2005, 2006). We argue that this model works as a dark fluid model which can interpolate between dust equation of state and the dark energy equation of state. We show how the asymptotic behavior of the equation of state constrained the parameters of the model. The causality condition for the model is also studied to constrain the parameters and the fixed points are tested to determine different solution classes. Observations of Hubble diagram of SNe Ia supernovae are used to further constrain the model. We present an exact solution of the model and calculate the luminosity distance and the energy density evolution. We also calculate the deceleration parameter to test the state of the universe expansion.     

Applications of Szebehely's equation

It is shown, that the potential obtained from Joukovsky's formula, corresponding to a given family of orbits is a general solution of Szebehely's equation. Then it is shown how a general solution of Szebehely's equation can be obtained from its particular solution. This method is applied to several examples. Potentials generating families of concentric elliptic orbits and families of orbits of conic sections are determined. Finally, the inverse Keplerian problem is solved using Szebehely's equation in polar coordinates.     

On spin-up/spin-down torque reversals in disc accreting pulsars

The recent BATSE observations of the spin-up and spin-down of accreting pulsars have shown that the standard formulation for the accretion torque as proposed by Ghosh &38; Lamb may need to be revised. The observations indicate alternate spin-up and spin-down phases driven by torques of similar magnitude and typically larger than the mean torque. The variations of the torque in systems such as Cen X-3 are difficult to explain in terms of changes of the mass accretion rate. The implication is that the torque does not depend on the accretion rate as in the GL model. In this paper we argue that the observed changes in the spin rate can result from stochastic transitions between two magnetospheric states. In particular, we show that intermediate magnetospheric systems are not admissible, because of a disc-induced magnetospheric instability which exists in a star–disc magnetic interaction system. This explains why torque reversal occurs in disc accreting pulsars with similar magnitudes.     

Magnetic torques between accretion discs and stars

I show in this paper that two types of magnetic torques can appear in the interaction between an accretion disc and a magnetic accretor. There is the well-known torque resulting from the difference in angular velocity between the accretion disc and the star, but in addition there is a torque coming from the interaction between the stellar magnetic field and the disc's own magnetic field. The latter form of magnetic torque decreases in strength more slowly with increasing radius, and will therefore dominate at large radii. The direction of the disc field is not determined by the difference in angular velocity between the star and the disc as in the Ghosh &38; Lamb model, but rather is a free parameter. The magnetic torque may therefore either spin up or spin down the star, and the torque changes sign if the magnetic field in the disc reverses. I suggest that this mechanism can explain the torque reversals that have been observed in some disc-fed X-ray pulsars.     

Control Of Chaos In Hamiltonian Systems

We present a technique to control chaos in Hamiltonian systems which are close to integrable. By adding a small and simple control term to the perturbation, the system becomes more regular than the original one. We apply this technique to a forced pendulum model and show numerically that the control is able to drastically reduce chaos.     

Monoenergetic-source solutions of the steady-state cosmic-ray equation of transport

A non spherically-symmetric monoenergetic-point-source solution of the steady-state equation of transport for cosmic-rays in the interplanetary region, in which monoenergetic particles are released isotropically and continuously from a fixed heliocentric position is derived by a Laplace transform method. The solution is for a spherically-symmetric model of the propagating region incorporating anisotropic diffusion, with a diffusion tensor symmetric about the radial direction, and the solar wind velocity is radial and of constant speedV. The spherically-symmetric monoenergeticsource solution of Webb and Gleeson (1973) and of Toptygin (1973) is regained from the spherically-symmetric component of the point-source solution.     

Wave disturbances near the temperature jump in the transition region of the solar atmosphere

We theoretically analyzed the properties of surface waves near the temperature jump in the solar atmosphere whose dispersion relation is identical in form to the equation for waves on deep water. We found an exact solution to the model equation for the vertical velocity of the medium in such a wave. Based on the derived space-time dependence of the vertical velocity of the medium, we quantitatively explained one of the events recorded in the SOI/MDI experiment onboard the SOHO spacecraft that accompanied coronal mass ejection.     

An efficient computational method for the approximate solution of nonlinear Lane-Emden type equations arising in astrophysics

The key purpose of this article is to introduce an efficient computational method for the approximate solution of the homogeneous as well as non-homogeneous nonlinear Lane-Emden type equations. Using proposed computational method given nonlinear equation is converted into a set of nonlinear algebraic equations whose solution gives the approximate solution to the Lane-Emden type equation. Various nonlinear cases of Lane-Emden type equations like standard Lane-Emden equation, the isothermal gas spheres equation and white-dwarf equation are discussed. Results are compared with some well-known numerical methods and it is observed that our results are more accurate.