Analytic continuation of quasi-periodic orbits from non-linear periodic modes

A new theory is formulated for the analytic continuation of quasi-periodic orbits from non-linear, periodic modes of a two degree-of-freedom system in rotating coordinates. Formal solutions of the equations of motion are developed in series expansions, with the non-linear modes as reference solutions. Conditions are determined on the existence and boundedness of the time-dependent coefficients of the velocity field. Boundedness properties of the orbits are specified by an orbit function derived from the Jacobi integral.     

On magnetostatic equilibrium in a stratified atmosphere

This is a study of the relationship between a magnetic field and its embedding plasma in static equilibrium in a uniform gravity. The ideal gas law is assumed. A system invariant in a given direction is treated first. We show that an exact integral of the equation for force balance across field lines can be derived in a closed form. Using this integral, exact solutions can be generated freely by integrating directly for the distributions of pressure, density and temperature necessary to keep a given magnetic field in equilibrium. Particular solutions are presented for illustration with the solar atmosphere in mind. Extending the treatment to the general system depending on all three spatial coordinates, we arrive at the general form of a theorem of Parker that a magnetic field in static equilibrium must possess certain symmetries. We derive an equation involving the Euler potentials of the magnetic field stipulating these necessary symmetries. Only those magnetic fields satisfying this equation can be in static equilibrium and for these fields, the endowed symmetries make the construction of exact solutions an essentially two dimensional problem as exemplified by the special case of invariance in a given direction.     

Equilibria in the secular, non-co-planar two-planet problem

We investigate the secular dynamics of a planetary system composed of the parent star and two massive planets in mutually inclined orbits. The dynamics are investigated in wide ranges of semimajor axes ratios (0.1–0.667) and planetary masses ratios (0.25–2), as well as in the whole permitted ranges of the energy and total angular momentum. The secular model is constructed by semi-analytic averaging of the three-body system. We focus on equilibria of the secular Hamiltonian (periodic solutions of the full system) and we analyze their stability. We attempt to classify families of these solutions in terms of the angular momentum integral. We identified new equilibria, yet unknown in the literature. Our results are general and may be applied to a wide class of three-body systems, including configurations with a star and brown dwarfs and substellar objects. We also describe some technical aspects of the seminumerical averaging. The HD 12661 planetary system is investigated as an example configuration.     

The second integral of the restricted problem in regularized elements

A perturbation series integral for the restricted problem of three bodies is derived by use of a new set of canonical elements for the regularized two-body problem. These elements are similar to theKS elements of Stiefel and Scheifele, but they contain small parameters other than the semimajor axis. The variable analogous to the longitude of perihelion not only remains well defined as the orbit approaches a circle, but also it can be used as a second small parameter. Regularized elements permit canonical use of the eccentric anomaly as independent variable, but most of the major benefits of regularization in the two-body problem do not carry over to perturbation theory.     

A new technique for exact and unique solution of transfer equations in finite media

In this paper we develop a new method, combined with Laplace transformation and Wiener-Hopf technique, to obtain unique solutions of transport equations in finite media. For this purpose we consider the simple transfer equation for diffuse reflection by a plane-parallel finite atmosphere scattering radiation with moderate anisotropy. It is transformed, by Laplace transformation, into two coupled linear integral equations which are then reduced to two uncoupled Fredholm integral equations admitting of unique solutions by the method of iteration for values of the breadth of the atmosphere greater than that specified, depending on the scattering process.     

Circular restricted three-body problem when both the primaries are heterogeneous spheroid of three layers and infinitesimal body varies its mass

The circular restricted three-body problem, where two primaries are taken as heterogeneous oblate spheroid with three layers of different densities and infinitesimal body varies its mass according to the Jeans law, has been studied. The system of equations of motion have been evaluated by using the Jeans law and hence the Jacobi integral has been determined. With the help of system of equations of motion, we have plotted the equilibrium points in different planes (in-plane and out-of planes), zero velocity curves, regions of possible motion, surfaces (zero-velocity surfaces with projections and Poincaré surfaces of section) and the basins of convergence with the variation of mass parameter. Finally, we have examined the stability of the equilibrium points with the help of Meshcherskii space–time inverse transformation of the above said model and revealed that all the equilibrium points are unstable.     

Coherent backscattering and opposition effects observed in some atmosphereless bodies of the solar system

The results of photometric and polarimetric observations carried out for some bright atmosphere-less bodies of the Solar system near the zero phase angle reveal the simultaneous existence of two spectacular optical phenomena, the so-called brightness and polarization opposition effects. In a number of studies, these phenomena were explained by the influence of coherent backscattering. However, in general, the interference concept of coherent backscattering can be used only in the case where the particles are in the far-field zones of each other, i.e., when the scattering medium is rather rarefied. Because of this, it is important to prove rigorously and to demonstrate that the coherent backscattering effect may also exist in densely packed scattering media like regolith surface layers of celestial bodies. From the results of the computer modeling performed with the use of numerically exact solutions of the macroscopic Maxwell equations for discrete random media with different packing densities of particles, we studied the origin and evolution of all the opposition phenomena predicted by the coherent backscattering theory for low-packing-density media. It has been shown that the predictions of this theory remain valid for rather high packing densities of particles that are typical, in particular, of regolith surfaces of the Solar system bodies. The results allow us to conclude that both opposition effects observed simultaneously in some high-albedo atmosphereless bodies of the Solar system are caused precisely by coherent backscattering of solar light in the regolith layers composed of microscopic particles.     

Equilibrium points and zero velocity surfaces in the restricted four-body problem with solar wind drag

We have analyzed the motion of an infinitesimal mass in the restricted four-body problem with solar wind drag. It is assumed that the forces which govern the motion are mutual gravitational attractions of the primaries, radiation pressure force and solar wind drag. We have derived the equations of motion and found the Jacobi integral, zero velocity surfaces, and particular solutions of the system. It is found that three collinear points are real when the radiation factor 0<β<0.1 whereas only one real point is obtained when 0.125<β<0.2. The stability property of the system is examined with the help of Poincaré surface of section (PSS) and Lyapunov characteristic exponents (LCEs). It is found that in presence of drag forces LCE is negative for a specific initial condition, hence the corresponding trajectory is regular whereas regular islands in the PSS are expanded.     

On a Bernoulli's integral pertaining to gas flow in close binary systems

A Bernoulli's integral supplemented with the equation of continuity provides a solution for the motion of gas surrounding a binary system.There exist two velocity modes whose streamlines are confined within appropriate equipotential surfaces.     

The use of integrals in numerical integrations of theN-body problem

The numerical integration of systems of differential equations that possess integrals is often approached by using the integrals to reduce the number of degrees of freedom or by using the integrals as a partial check on the resulting solution, retaining the original number of degrees of freedom.Another use of the integrals is presented here. If the integrals have not been used to reduce the system, the solution of a numerical integration may be constrained to remain on the integral surfaces by a method that applies corrections to the solution at each integration step. The corrections are determined by using linearized forms of the integrals in a least-squares procedure.The results of an application of the method to numerical integrations of a gravitational system of 25-bodies are given. It is shown that by using the method to satisfy exactly the integrals of energy, angular momentum, and center of mass, a solution is obtained that is more accurate while using less time of calculation than if the integrals are not satisfied exactly. The relative accuracy is ascertained by forward and backward integrations of both the corrected and uncorrected solutions and by comparison with more accurate integrations using reduced step-sizes.     

Exact solution of a basic equation in finite atmosphere by the method of Laplace transform and linear singular operators

A finite atmosphere having distribution of intensity at both surfaces with definite form of scattering function and source function is considered here. The basic integro-differential equation for the intensity distribution at any optical depth is subjected to the finite Laplace transform to have linear integral equations for the surface quantities under interest. These linear integral equations are transformed into linear singular integral equations by use of the Plemelj's formulae. The solution of these linear singular integral equations are obtained in terms of theX-Y equations of Chandrasekhar by use of the theory of linear singular operators which is applied in Das (1978a).     

Solution of Riemann–Hilbert problems for determination of new decoupled expressions of Chandrasekhar’s <Emphasis Type="Italic">X</Emphasis>- and <Emphasis Type="Italic">Y</Emphasis>-functions for slab geometry in radiative transfer

In radiative transfer, the intensities of radiation from the bounding faces of a scattering atmosphere of finite optical thickness can be expressed in terms of Chandrasekhar’s X- and Y-functions. The nonlinear nonhomogeneous coupled integral equations which the X- and Y-functions satisfy in the real plane are meromorphically extended to the complex plane to frame linear nonhomogeneous coupled singular integral equations. These singular integral equations are then transformed into nonhomogeneous Riemann–Hilbert problems using Plemelj’s formulae. Solutions of those Riemann–Hilbert problems are obtained using the theory of linear singular integral equations. New forms of linear nonhomogeneous decoupled expressions are derived for X- and Y-functions in the complex plane and real plane. Solutions of these two expressions are obtained in terms of one known N-function and two new unknown functions N ₁- and N ₂- in the complex plane for both nonconservative and conservative cases. The N ₁- and N ₂-functions are expressed in terms of the known N-function using the theory of contour integration. The unknown constants are derived from the solutions of Fredholm integral equations of the second kind uniquely using the new linear decoupled constraints. The expressions for the H-function for a semi-infinite atmosphere are obtained as a limiting case.     

The transfer of polarized radiation in spectral lines: Formalism and solutions in simple cases

A general treatment of the transfer of polarized radiation in spectral lines assuming a Rayleigh phase function and a general law of frequency redistribution is derived. It is shown how nine families of coupled integral equations for the moments of the radiation field arise which are necessary to fully describe the state of polarization of the emergent radiation from a plane-parallel, semi-infinite atmosphere. The special case of angle independent redistribution functions is derived from the general formalism, and it is shown how the nine families of integral equations reduce to the six linearly independent integral equations derived by Collins (1972). To serve as a test of the formulation, solutions for isothermal atmospheres are given.     

A new class of relativistic stellar models

Einstein field equations for a static and spherically symmetric perfect fluid are considered. A formulation given by Patiño and Rago is used to obtain a class of nine solutions, two of them are Tolman solutions I, IV and the remaining seven are new. The solutions are the correct ones corresponding to expressions derived by Patiño and Rago which have been shown by Knutsen to be incorrect. Similar to Tolman solution IV each of the new solutions satisfies energy conditions inside a sphere in some range of two independent parameters. Besides, each solution could be matched to the exterior Schwarzschild solution at a boundary where the pressure vanishes and thus the solutions constitute a class of new physically reasonable stellar models.     

The Parker Problem and the Theory of Coronal Heating

To illustrate his theory of coronal heating, Parker initially considers the problem of disturbing a homogeneous vertical magnetic field that is line-tied across two infinite horizontal surfaces. It is argued that, in the absence of resistive effects, any perturbed equilibrium must be independent of z. As a result random footpoint perturbations give rise to magnetic singularities, which generate strong Ohmic heating in the case of resistive plasmas. More recently these ideas have been formalized in terms of a magneto-static theorem but no formal proof has been provided. In this paper we investigate the Parker hypothesis by formulating the problem in terms of the fluid displacement. We find that, contrary to Parker's assertion, well-defined solutions for arbitrary compressibility can be constructed which possess non-trivial z-dependence. In particular, an analytic treatment shows that small-amplitude Fourier disturbances violate the symmetry ∂_z = 0 for both compact and non-compact regions of the (x, y) plane. Magnetic relaxation experiments at various levels of gas pressure confirm the existence and stability of the Fourier mode solutions. More general footpoint displacements that include appreciable shear and twist are also shown to relax to well-defined non-singular equilibria. The implications for Parker's theory of coronal heating are discussed.     

Chemically reacting flow of a compressible thermally radiating two-component plasma

The paper studies the compressible flow of a hot two-component plasma in the presence of gravitation and chemical reaction in a vertical channel. For the optically thick gas approximation, closed form analytical solutions are possible. Asymptotic solutions are also obtained for the general differential approximation when the temperatures of the two bounding walls are the same. In the general case the problem is reduced to the solution of standard nonlinear integral equations which can be tackled by iterative precedure. The results are discussed quantitatively. The problem may be applicable to the understanding of explosive hydrogen-burning model of solar flares.     

Partial Reduction in the N-Body Planetary Problem using the Angular Momentum Integral

We present a new set of variables for the reduction of the planetary n-body problem, associated to the angular momentum integral, which can be of any use for perturbation theory. The construction of these variables is performed in two steps. A first reduction, called partial is based only on the fixed direction of the angular momentum. The reduction can then be completed using the norm of the angular momentum. In fact, the partial reduction presents many advantages. In particular, we keep some symmetries in the equations of motion (d'Alembert relations). Moreover, in the reduced secular system, we can construct a Birkhoff normal form at any order. Finally, the topology of this problem remains the same as for the non-reduced system, contrarily to Jacobi's reduction where a singularity is present for zero inclinations. For three bodies, these reductions can be done in a very simple way in Poincaré's rectangular variables. In the general n-body case, the reduction can be performed up to a fixed degree in eccentricities and inclinations, using computer algebra expansions. As an example, we provide the truncated expressions for the change of variable in the 4-body case, obtained using the computer algebra system TRIP.     

Exact solution of the problem of diffuse reflection and transmission by the method of Laplace transform and linear singular operators

The basic integro-differential equation is subjected to a one-sided finite Laplace transform to obtain linear integral equations of angular distribution of bounding faces. These linear integral equations have been transformed into linear singular integral equations which have been solved exactly to get the emergent distributions from the bounding faces by the theory of linear singular operators. Some solutions of linear singular integral equations have also been derived for future use in radiative transfer problems.     

The theory of asynchronous relative motion I: time transformations and nonlinear corrections

Using alternative independent variables in lieu of time has important advantages when propagating the partial derivatives of the trajectory. This paper focuses on spacecraft relative motion, but the concepts presented here can be extended to any problem involving the variational equations of orbital motion. A usual approach for modeling the relative dynamics is to evaluate how the reference orbit changes when modifying the initial conditions slightly. But when the time is a mere dependent variable, changes in the initial conditions will result in changes in time as well: a time delay between the reference and the neighbor solution will appear. The theory of asynchronous relative motion shows how the time delay can be corrected to recover the physical sense of the solution and, more importantly, how this correction can be used to improve significantly the accuracy of the linear solutions to relative motion found in the literature. As an example, an improved version of the Clohessy-Wiltshire (CW) solution is presented explicitly. The correcting terms are extremely compact, and the solution proves more accurate than the second and even third order CW equations for long propagations. The application to the elliptic case is also discussed. The theory is not restricted to Keplerian orbits, as it holds under any perturbation. To prove this statement, two examples of realistic trajectories are presented: a pair of spacecraft orbiting the Earth and perturbed by a realistic force model; and two probes describing a quasi-periodic orbit in the Jupiter-Europa system subject to third-body perturbations. The numerical examples show that the new theory yields reductions in the propagation error of several orders of magnitude, both in position and velocity, when compared to the linear approach.     

The structure of MHD shock waves in a thermally-radiating gas at high mach numbers

The structure of a strong MHD shock wave which radiates thermally downstream of the shock is studied by asymptotic expansion. The exact integral equation for radiation is adopted for the study. Hence, the optically thick (and thin), the general differential approximate and the exact integral equation solutions may now be compared.