General time elements for the two body problem

When integrating a perturbed two-body problem, very often the propagation of the numerical error is reduced by using a new time variables defined by dt/ds=|q| ⁿ , (|q| is the radial distance,t the time). This paper introduces a time element for such transformations, i.e., a new variablet _n is defined so that dt _n/ds=1+ (perturbing terms) andt=F _n(t_n), whereF _n is a known function. The time element equation should be useful in reducing the error in the determination of the timet.F _n is given explicitly forn=1, 3/2, 2, 5/2 and 3, and a general expression is given for other values.The work was performed while the author was an NRC Senior Research Associate, Goddard Space Flight Center, Greenbelt, Md., U.S.A.     

Solution of a time-dependent equation of transfer by theF n -method

The time-dependent equation of transfer for a finite, plane-parallel, non-radiating, and isotropically-scattering atmosphere of arbitrary stratification is solved by using theF _n -method developed by Siewert.     

The influence of the great inequality on the secular disturbing function of the planetary system

This paper derives the contributionF ₂ ^* by the great inequality to the secular disturbing function of the principal planets. Andoyer's expansion of the planetary disturbing function and von Zeipel's method of eliminating the periodic terms is employed; thereby, the corrected secular disturbing function for the planetary system is derived. An earlier solution suggested by Hill is based on Leverrier's equations for the variation of elements of Jupiter and Saturn and on the semi-empirical adjustment of the coefficients in the secular disturbing function. Nowadays there are several modern methods of eliminating periodic terms from the Hamiltonian and deriving a purely secular disturbing function. Von Zeipel's method is especially suitable. The conclusion is drawn that the canonicity of the equations for the secular variation of the heliocentric elements can be preserved if there be retained, in the secular disturbing function, terms only of the second and fourth order relative to the eccentricity and inclinations.The Krylov-Bogolubov method is suggested for eliminating periodic terms, if it is desired to include the secular perturbations of the fifth and higher order in the heliocentric elements. The additional part of the secular disturbing functionF ₂ ^* derived in this paper can be included in existing theories of the secular effects of principal planets. A better approach would be to preserve the homogeneity of the theory and rederive all the secular perturbations of principal planets using Andoyer's symbolism, including the part produced by the great inequality.     

Radiation transfer in a diffuse and specular reflecting slab with Rayleigh scattering

A method of analysis is presented for solving the radiative transfer problem in an absorbing, emitting, inhomogeneous, and anisotropically scattering plane-parallel medium with specular and diffuse reflecting boundaries and internal source (problem 1). Exact relations for the radiation heat flux at the boundaries of problem 1 are obtained in terms of the radiation density and albedos of the corresponding source-free medium with specular reflecting boundaries (problem 2). Two coupled integral equations for the radiation density and the second moment of the radiation intensity for problem 2 with Rayleigh phase functions are obtained. The Galerkin method is used to solve these equations. Albedos of problem 2 are compared with theF _n method. Numerical results for radiation heat fluxes at the boundaries of problem 1 are tabulated for different forms of the internal source.     

Dissipation and stability of return currents in solar flares

The dynamics of an electron beam, under the effects of Coulomb collisions and classical Ohmic dissipation of the return current, is analysed for a background plasma with a temperature which is time dependent due to the heating effect of beam dissipation offset by thermal conductive cooling. It is shown that the plasma is heated toward a steady state, in time scales short compared to typical flare beam switch on times, and that in this steady state only two regimes of beam dynamics arise.For moderate values of the ratio of beam flux to plasma density (F ₀/n), beam dynamics is dominated by direct Coulomb collisions as in the usual thick target treatment. With increase of F ₀/n, before return current (classical) Ohmic losses can exceed collisions, the return current becomes unstable. In this latter regime beam dynamics is presumably dominated by wave generation and anomalous Ohmic dissipation of the return current, but no detailed treatment is attempted here.     

Normalisation des systèmes linéaires canoniques et application au problème restreint des trois corps

An algorithm is given for normalizing conservative linear Hamiltonian systems. This one generalizes Siegel's method to the cases where the eigenvalues are multiple. We obtain by a canonical transformation a normal form of two blocks, one of which is the upper Jordan form, and the other, the lower Jordan form. We select real solutions from the solutions of these equations, and we apply the result to the restricted three body problem in the vicinity of the triangular points for Routh's critical mass ratio.

---

---

_n est la matrice unité d'ordren,O _n est la matrice nulle d'ordren.     

Cosmological models with the variable gravitational and cosmological constants

Cosmology with the gravitational and cosmological constants generalized as coupling scalars in Einsteins theory is considered. A general method of solving the field equations is given. Exact solution for Zeldovich fluid satisfying G=G ₀(R/R ₀) ⁿ is given.     

Cosmological models with generalised gravitational and cosmological constants

We consider cosmology with the gravitational and cosmological constants generalized as coupling scalars in Einstein’s theory. A general method of solving the field equations is given. We study here the exact solutions for negative pressure models satisfying G=G ₀(R/R ₀) ⁿ .     

The construction of a third-order secular analytical J-S-U-N theory by the Hori-Lie technique

We expand both parts, the principal and indirect, of the Hamiltonian function up to the third order in the masses for the four major planets Jupiter-Saturn-Uranus-Neptune. Accordingly we write down the secular terms ofF ₁,F ₂,F ₃ and the critical terms ofF ₁,F ₂ in terms of the canonical variables of H. Poincaré neglecting powers higher than the second inH, K, P, Q.     

Numerical solution of the time-dependent heat energy equations for the magnetosphere

A new solution of the magnetospheric heat equations capable of covering the whole region from 300 km along a field line to the equatorial plane has been achieved by adapting the searching procedure of Murphy (1974). It has been found that the protonospheric heat reservoir is sufficient to maintain T_e >T_n down to the height of the F2-peak electron density all through the night at mid-latitudes. Full solution of the equations has also shown that T_i >T_e in the protonosphere at night and the ions constitute a significant source of heat for the electrons.     

Cavitation by Nonlinear Interaction Between Inertial Alfvén Waves and Magnetosonic Waves in Low Beta Plasmas

This paper presents the model equations governing the nonlinear interaction between dispersive Alfvén wave (DAW) and magnetosonic wave in the low-β plasmas (β≪m _e/m _i; known as inertial Alfvén waves (IAWs); here \upbeta = 8pn₀T /B₀²\upbeta = 8\pi n_{0}T /B_{0}^{2} is thermal to magnetic pressure, n ₀ is unperturbed plasma number density, T(=T _e≈T _i) represents the plasma temperature, and m _e(m _i) is the mass of electron (ion)). This nonlinear dynamical system may be considered as the modified Zakharov system of equations (MZSE). These model equations are solved numerically by using a pseudo-spectral method to study the nonlinear evolution of density cavities driven by IAW. We observed the nonlinear evolution of IAW magnetic field structures having chaotic behavior accompanied by density cavities associated with the magnetosonic wave. The relevance of these investigations to low-β plasmas in solar corona and auroral ionospheric plasmas has been pointed out. For the auroral ionosphere, we observed the density fluctuations of ∼ 0.07n ₀, consistent with the FAST observation reported by Chaston et al. (Phys. Scr. T84, 64, 2000). The heating of the solar corona observed by Yohkoh and SOHO may be produced by the coupling of IAW and magnetosonic wave via filamentation process as discussed here.     

MGS electron density profiles: Analysis of the peak magnitudes

We analyze here the behavior of the magnitudes of the F₁ and E peaks of the electron density profiles measured by the Radio Science Subsystem of the Mars Global Surveyor spacecraft, as a function of solar zenith angle χ and solar flux. For each of the 658 days of data in the six occultation seasons in the northern hemisphere, we choose one profile to analyze, which is that for which the F₁ peak is the median value. We assume that the variations of the measured peak densities can be represented as Aa(cosχ) and as Bb(F_10.7), where F_10.7 is the usual solar flux proxy, appropriately shifted to the orbital position of Mars. To minimize the effect of solar activity, we divide the data into 6 F_10.7 bins, fit the data in each bin, and derive the values of the exponent a and the coefficient AF10.7 for each bin. The median values that we derive for the exponent a is 0.46 for the F₁ peak, and 0.395 for the E peak. To minimize the effect of SZA, we divide the data into eight SZA bins, and derive the exponent b and the coefficient Bχ for each SZA bin. We argue that the last three SZA bins should be excluded because the fits were poor, due partly to the small number of data points in each of these bins. If we do so, the median values of b that we derive are 0.27 and 0.40 for the F₁ and E peaks, respectively. Finally we derive a 3-parameter fit to all the data, which expresses the variability of the peak densities as a function of a^′(cosχ) and b^′(F_10.7) simultaneously. The fitted values of the exponents a^′ and b^′ for the F₁ peak are 0.45 and 0.26, respectively; for the E peak, the values are 0.39 and 0.46, respectively. We compare our results to Chapman theory, and to those of other investigators.     

General precession computed by numerical integration of ordinary differential equations

It is suggested that the effects of change in the celestial coordinates due to general precession be computed by solving numerically the differential equations for the instataneous rates of change in right ascension and declination. An approximate solution in closed from is given. Similar differential equations may be solved for the complete rotation matrix, requiring only the knowledge of the variations of the two function m and n with time. This suggests the existence of a relation between the classical precessional elements z, ζ₀ and Θ, and this relation is derived. Some numerical examples are also presented.     

On role of pressure anisotropy for relativistic stars admitting conformal motion

We investigate the spacetime of anisotropic stars admitting conformal motion. The Einstein field equations are solved using different ansatz of the surface tension. In this investigation, we study two cases in details with the anisotropy as: (1) p _t =n p _r , (2) -p_r=(+c_2)p_{t}-p_{r}=\frac{1}{8\pi}(\frac{c_{1}}{r^{2}}+c_{2}) where, n, c ₁ and c ₂ are arbitrary constants. The solutions yield expressions of the physical quantities like pressure gradients and the mass.     

On the Tetrad Formulation of the Equations of Motion in General Relativity

It is shown that the equation of motion Du _j/Ds = (e/mc ²)F _ji u ⁱ , a natural generalization to the curved spacetime of the Heaviside-Lorentz law of ponderomotive force, is equivalent to the metric independent and covariant Van Dantzig's equations of motion dx ^j _[jpi] = 0 or L _v p _i = 0, where p _i is the conjugate momentum 4-vector and v a vector determined by the condition p _i v ⁱ, only with respect to holonomic coordinates. With respect to an anholonomic system, the Heaviside-Lorentz equation is a particular case of the VD equations valid for a privileged class of anholonomic frames, those consisting of orthogonal unit vectors.     

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.     

Variety of Well behaved parametric classes of relativistic charged fluid spheres in general relativity

The paper presents a variety of classes of interior solutions of Einstein–Maxwell field equations of general relativity for a static, spherically symmetric distribution of the charged fluid with well behaved nature. These classes of solutions describe perfect fluid balls with positively finite central pressure, positively finite central density; their ratio is less than one and causality condition is obeyed at the center. The outmarch of pressure, density, pressure–density ratio and the adiabatic speed of sound is monotonically decreasing for these solutions. Keeping in view of well behaved nature of these solutions, two new classes of solutions are being studied extensively. Moreover, these classes of solutions give us wide range of constant K for which the solutions are well behaved hence, suitable for modeling of super dense star. For solution (I₁) the mass of a star is maximized with all degree of suitability and by assuming the surface density ρ _b =2×10¹⁴ g/cm³ corresponding to K=1.19 and X=0.20, the maximum mass of the star comes out to be 2.5M _Θ with linear dimension 25.29 Km and central redshift 0.2802. It has been observed that with the increase of charge parameter K, the mass of the star also increases. For n=4,5,6,7, the charged solutions are well behaved with their neutral counterparts however, for n=1,2,3, the charged solution are well behaved but their neutral counterparts are not well behaved.     

Nonuniform Target Ionization and Fitting Thick Target Electron Injection Spectra to RHESSI Data

Past analyses of flare hard X-ray (HXR) spectra have largely ignored the effect of nonuniform ionization along the electron paths in the thick-target model, though it is very significant for well-resolved spectra. The inverse problem (photon spectrum to electron injection spectrum F ₀(E ₀)) is disturbingly non-unique. However, we show that it is relatively simple to allow for the effect in forward fitting of parametric models of F ₀(E ₀)) and provide an expression to evaluate it for the usual single power-law form of F ₀(E ₀)).The expression involves the column depth N _* of the transition region in the flare loop as one of the parameters so data fitting can enable derivation of N _* (and its evaporative evolution) as part of the fitting procedure. The fit to RHESSI data on four flares for a single power law F ₀(E ₀)) is much improved when ionization structure is included compared to when the usual fully ionized approximation is used. This removes the need, in these events at least, to invoke broken power laws, or other forms, of the acceleration spectrum F ₀(E ₀)) to explain the observed photon spectrum     

TheF N method for solving radiative-transfer problems in plane geometry

TheF _N method of solving problems in radiative transfer for plane-parallel media with anisotropic scattering is established. The method utilizes properties of the exact solution and leads to final equations that are particularly concise and easy to use.