Boundary curves for families of planar orbits

For monoparametric familiesf(x,y)=c of planar orbits, created by a planar potentialV(x,y), we introduce the notion of the family boundary curves (FBC). All members of the familyf(x,y)=c are traced in an allowable region of thexy plane, defined by the corresponding FBC, with total energyE=E(c) varying along the family. Family boundary curves are also found for two-parametric familiesf(x,y,b)=c. The relation of equilibrium points and asymptotic orbits, possibly possessed by the potentialV(x,y), to be FBC is studied.     

An equation for the characteristic curve of a family of symmetric periodic orbits

Szebehely's equation for the inverse problem of Dynamics is used to obtain the equation of the characteristic curve of a familyf(x,y)=c of planar periodic orbits (crossing perpendicularly thex-axis) created by a certain potentialV(x,y). Analytic expressions for the characteristic curves are found both in sideral and synodic systems. Examples are offered for both cases. It is shown also that from a given characteristic curve, associated with a given potential, one can obtain an analytic expression for the slope of the orbit at any point.     

Generalization of Szebehely's equation

Szebehely's partial differential equation for the force functionU=U(x,y) which gives rise to a given family of planar orbitsf(x,y)=Constant is generalized to account for velocity-dependent potentials V^*=V^*(x,y, ). The new partial differential equation is quasi-linear and of the first order. An example is given and a comparison is made of the two equations.     

Families of planar orbits generated by homogeneous potentials

The second order partial differential equation which relates the potentialV(x,y) to a family of planar orbitsf(x,y)=c generated by this potential is applied for the case of homogeneousV(x,y) of any degreem. It is shown that, if the functionf(x,y) is also homogeneous, there exists, for eachm, a monoparametric set of homogeneous potentials which are the solutions of an ordinary, linear differential equation of the second order. Iff(x,y) is not homogeneous, in general, there is not a homogeneous potential which can create the given family; only if =f _y /f _x satisfies two conditions, a homogeneous potential does exist and can be determined uniquely, apart from a multiplicative constant. Examples are offered for all cases.     

Determination of autonomous three-dimensional force fields from a two-parameter family of orbits

Any two of the componentsX, Y, andZ of an autonomous force field which gives rise to the space orbitsF(x, y, z)=c ₁,G(x, y, z)=c ₂ are related by a partial differential equation with coefficients depending on the functionsF andG. This is a generalization of the corresponding equation for planar orbits (Bozis, 1983). The above partial differential equation is accompanied by the algebraic linear equation inX, Y, andZ expressing the fact that the force vector is lying in the osculating plane at each point of the orbit. The two equations constitute a generalization of the corresponding Szebehely's equations in the three dimensional space (Érdi, 1982). The generalization is meant in the sense that the dynamical system is not necessarily assumed to be conservative.     

On the adelphic potentials compatible with a set of planar orbits

Given a planar potentialB=B(x, y), compatible with a monoparametric family of planar orbitsf(x, y)=c, we face the problem of producing potentialsA=A(x, y), adelphic toB(x, y), i.e. nontrivial potentials which have in common withB(x, y) the given set of orbits. We establish a linear, second order partial differential equation for a functionP(x, y) and we prove that, to any definite positive solution of this equation, there corresponds a potentialA(x, y) adelphic toB(x, y).     

On numerical evaluation of theH-functions of transport problems by kernel approximation for the albedo 0<ω≤1

Das Gupta represented theH-functions of transport problems for the albedo [0, 1] in the formH(z)=R(z)–S(z) (see Das Gupta, 1977) whereR(z) is a rational function ofz andS(z) is regular on [–1, 0] ^c . In this paper we have representedS(z) through a Fredholm integral equation of the second kind with a symmetric real kernelL(y, z) as . The problem is then solved as an eigenvalue problem. The kernel is converted into a degenerate kernel through finite Taylor's expansion and the integral equation forS(z) takes the form: (which is solved by the usual procedure) where _r 's are the discrete eigenvalues andF _r 's the corresponding eigenfunctions of the real symmetric kernelL(y, z).     

OnH-functions of radiative transfer

Chandrasekhar'sH-functionH(z) corresponding to the dispersion functionT(z)=| _rs –frs(z)|, where [f _rs (z)] is of rank 1, is obtained in terms of a Cauchy integral whose density functionQ(x, ₁, ₂,...) can be approximated by approximating polynomials (uniformly converging toQ(x)) having their coefficients expressed as known functions of the parameters _r 's. A closed form approximation ofH(z) to a sufficiently high degree of accuracy is then readily available by term by term integration.     

Inverse problem with two-parametric families of planar orbits

As a possible extension of recent work we study the following version of the inverse problem in dynamics: Given a two-parametric familyf(x, y, b)=c of plane curves, find an autonomous dynamical system for which these curves are orbits.We derive a new linear partial differential equation of the first order for the force componentsX(x, y) andY(x, y) corresponding to the given family. With the aid of this equation we find that, depending on the given functionf, the problem may or may not have a solution. Based on given criteria, we present a full classification of the various cases which may arise.     

Transient Flow of a Magnetised Plasma due to Thermal perturbation

At time t<0, a steady stationary condition of a fully ionised gas exists in space, such that the velocity components, the induced magnetic field and a uniform temperature are given by (0,0,0),(0, 0,H_o) and Θ, respectively, in the rectangular coordinate (x,y,z). At t =0, a sudden increase of temperature is applied at z = 0. If the Eckert terms are neglected, an integral transform technique is employed to solve for the z–axis related transient flow (U,V,0),(H_o), and temperature Θ (z,t). For large applied magnetic field, H_o the flow is observed to exhibit disturbance modes some of which are oscilatory, and some streaming modes are seen accompanied by the expected decays. For a thermally perturbed plasma the flow is seen to be largely governed by conduction parameters if viscous terms are not neglected under the MHD approximations. This revised version was published online in July 2006 with corrections to the Cover Date.     

Velocity structure of the shallow lunar crust

The data from the Apollo-14 and Apollo-16 Active Seismic Experiments have been reanalyzed and show that a power-law velocity variation with depth,v(z)110z ^1/6 m s^–1 (0<z<10 m), is consistent with both the travel times and amplitudes of the first arrivals for source-to-geophone separations up to 32 m. The data were improved by removing spurious glitches, by filtering and stacking. While this improved the signal-to-noise ratios, it was not possible to measure the arrival times or amplitudes of the first arrivals beyond 32 m. The data quality precludes a definitive distinction between the power-law velocity variation and the layered-velocity model proposed previously. However, the physical evidence that the shallow lunar regolith is made up of fine particles adds weight to the 1/6-power velocity model because this is the variation predicted theoretically for self-compacting spheres.The 1/6-power law predicts the travel time,t(x), varies with separation,x, ast(x)=t ₀(x/x ₀)^5/6 and, using a first-order theory, the amplitude,A(x), varies asA(x)=A ₀(x/x ₀)^{–(13–m)/12},m>1; the layervelocity model predictst(x)=t ₀(x/x ₀) andA(x)=A ₀(x/x ₀)^–2, respectively. The measured exponents for the arrival times were between 0.63 and 0.84 while those for the amplitudes were between –1.5 and –2.2. The large variability in the amplitude exponent is due, in part, to the coarseness with which the amplitudes are measured (only five bits are used per amplitude measurement) and the variability in geophone sensitivity and thumper-shot strengths.A least-squares analysis was devised which uses redundancy in the amplitude data to extract the geophone sensitivities, shot strengths and amplitude exponent. The method was used on the Apollo-16 ASE data and it indicates there may be as much as 30 to 40% variation in geophone sensitivities (due to siting and coupling effects) and 15 to 20% variability in the thumper-shot strengths. However, because of the low signal-to-noise ratios in the data, there is not sufficient accuracy or redundancy in the data to allow high confidence in these results.     

Instability of an area-preserving polynomial mapping

Fixed points and eigencurves have been studied for the Hénon-Heiles mapping:x′=x+a (y?y ³),y′=y(x′?x′ ³). Eigencurves of order 21 proceed rapidly to infinity fora=1.78, but as ‘a’ decreases, they spiral around the origin repeatedly before escaping to infinity. Fixed pointsx _f on thex-axis have been located for the range 1≤a≤2.4, for ordersn up to 100. Their locations vary continuously witha, as do the eigencurves, and hyperbolic points remain hyperbolic. Forn=3 and 2.4≥a≥2.37, a very detailed study has been made of how escape occurs, with segments of an eigencurve mapping to infinity through various escape channels. Further calculations with ‘a’ decreasing to 2.275 show that this instability is preserved and that the eigencurve will spiral many times around the origin before reaching an escape channel, there being more than 34 turns fora=2.28. The rapid increase of this number is associated with the rapid decrease of the intersection angle between forward and backward eigencurves (at the middle homoclinic point), with decreasing ‘a’, this angle governing the outward motion. By a semi-topological argument, it is shown that escape must occur if the above intersection angle is nonzero. In the absence of a theoretical expression for this angle, one is forced to rely on the numerical evidence. If the angle should attain zero for a valuea=a _c>a_m,wherea _m .is the minimum value for which the fixed points exist, then no escape would be possible fora _c However, on the basis of calculations by Jenkins and Bartlett (1972) forn=6, and the results of the present article forn=3, it appears highly probable thata _c=a_m,and that escape from the neighborhood of a hyperbolic point is always possible. If there is escape from the hyperbolic fixed point forn=4,a=1.6, located atx _f=0.268, then the eigencurve must cross the apparently closed invariant curve of Hénon-Heiles which intersects thex-axis atx?±0.4, so that this curve cannot in fact be closed.     

The proton-proton reaction and related reactions in an intense magnetic field

The reaction rates for the proton-proton reaction and the related electron capture reaction in a strongly magnetized relativistic electron gas of arbitrary degree of degeneracy are computed. The proton-proton reaction rates are unaffected by the presence of the magnetic field for field strengths up to the critical valueH _q =m ² c ³/e=4.414×10¹³G. For fields greater thanH _q , the proton-proton reaction rates are enhanced linearly with (=H/H _q ).The PeP reaction is investigated in detail for a wide range of temperatures, densities and magnetic field strengths that are of interest. The main results are as follows: In the non-degenerate regime the reaction rates are significantly reduced for high temperatures (T ₉5) and low fields (1). For instance, _pep(H)=0.04 _pep(O) at =10^–3,T ₉=10. For relatively high fields (>1) and low temperatures (T ₉2), the reaction rates are enhanced approximately linearly with . In the complete degenerate regime the reaction rates are reduced up to one-third of the field-free value for moderate densities (₆/_e10). At high densities (₆/_e10) the reaction rates are unaffected by the magnetic field.     

Analysis of red and infrared wide-band photometry of Algol

The basic physical picture of the Algol system is reviewed, and, using collected red and infrared observations, photometric curve fits are investigated by applying numerical quadratures to determine theoretical light curves appropriate to Roche model stars. The contact nature of Algol B appears to be confirmed, and effective temperatures of the three components areT _A=11 400°,T _B=5300°,T _C=7600°. In terms of a Lambert's law approach to the reflection effect, the effective heat-albedo is required to be reduced from unity to one half; and it is also found that the averaged gravity-darkening coefficient is close to a value appropriate for a diffusion type of heat-transfer mechanism operating in sub-photospheric layers.     

The interstellar band at λ 4430 Å and the abundance of interstellar iron,titanium, and molecular hydrogen

For the direction to a number of stars the depletion of interstellar gaseous iron and titanium as well as the relative abundance of molecular hydrogen and the strength of the interstellar band at 4430 Å were determined by different authors and can be found in the literature. In this paper it is shown that the difference (A _c-A_co) is a more reliable measurement of the column density of 4430 Å absorbers than the usually used central depth,A _c, because the positive valueA _co ofA _c forE(B-V)=0, i.e. the intercept with theA _c axis of a least-squares fit to the observedA _c vsE(B-V) data, is with a high probability not caused by an interstellar effect.There was no correlation found between the interstellar depletion of iron and titanium, respectively, and (A _c-A_co), whereas a tendency exists that with increasing relative molecular hydrogen abundance the number of 4430 Å absorbers per hydrogen atom decreases. If the carriers of the 4430 Å absorbers are interstellar grains, then these grains must be altered during the same process in which molecular hydrogen is built. The found correlation is also compatible with the assumption that the 4430 Å absorber is related to an interstellar gaseous species.     

Two‐dimension potentials which generate spatial families of orbits

We consider the following case of the 3D inverse problem of dynamics: Given a spatial two‐parametric family of curves f (x, y, z) = c₁, g (x, y, z) = c₂, find possibly existing two‐dimension potentials under whose action the curves of the family are trajectories for a unit mass particle. First we establish the conditions which must be fulfilled by the family so that potentials of the form w (y, z) give rise to the curves of the family, and we present some applications. Then we examine briefly the existence of potentials depending on (x, z), respectively (x, y), which are compatible with the given family (© 2009 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

A new integration theory for conservative two degree-of-freedom systems

A two degree-of-freedom, conservative system is reduced to a single degree-of-freedom, kinematic system with Hamiltonian integral under the change of independent variable: $$dt = \zeta dt (\zeta = \upsilon _x - \upsilon _y )$$ where ζ is the curl (or vorticity) of the velocity field with cartesian inertial componentsu(x, y, t) andv(x, y, t). In the autonomous case whenu _t=v _t=0, orbits are globally represented by the level curves of an autonomous Hamiltonian functionH(x,y) satisfying a second-order quasilinear partial differential equation (Szebehely's Equation): $$2(H + U)\left( {H_{xx} H_y^2 - 2H_{xy} H_x H_y + H_{yy} H_x^2 } \right) + (H_x U_x + H_y U_y )\left( {H_x^2 + H_y^2 } \right) = 0$$ whereU(x, y) is the autonomous potential function. An inversion of dependent and independent variables reduces this equation to a second-order, ordinary differential equation for a function specifying the orbital curve. The true time variable is recovered by evaluating a quadrature. Fundamental differences exist between this approach and Hamilton-Jacobi theory.     

The gravitational field of a disk

This note gives the gravitational potential of the disk {(x, y, z):x ² +y ² p ² , z=0} and the gravitational field at the point (x, y, z). Formulas for a ring can be obtained as the difference of our results for two different values ofp. Results are obtained in terms of elliptic integrals and we indicate how these functions can be computed efficiently. Formulas necessary for the computation of partial derivatives are also given.This paper presents the results of one phase of research carried out at the Jet Propulsion Laboratory, California Institute of Technology, under Contract NAS7-100, sponsored by the National Aeronautics and Space Administration.     

The lunar mean moment of inertia and the size of the Moon's core

In this paper the relation between the uncertainty of the Moon's mean moment of inertia (I/Ma ²) and that of the core density _c is discussed with a two-layer model of the Moon - a mantle obeying Roche's law of the density distribution and a homogeneous core (Fe-core or Fe-FeS-core). When the uncertainty of I/Ma ² is 0.0023 (that is the accuracy in present observation), a core with radius of 450 km will be appropriate to the limitation of _c about 1 g cm^–3. Considering the accuracy obtained in space explorations, and the compressibility and the quasi-homogeneity of the Moon, we suggest that the parameters C ₂₀, , , a, and GM of the Moon should define as primary constants, but C ₂₂ and C/Ma ² as derived constants. Therefore, the ratio of mass of Moon to that of Earth in the IAU (1976) system of astronomical constants will become a deducible constant.     

Planar electron-acoustic solitary waves and double layers in a two-electron-temperature plasma with nonthermal ions

A rigorous theoretical investigation of nonlinear electron-acoustic (EA) waves in a plasma system (containing cold electrons, hot electrons obeying a Boltzmann distribution, and hot ions obeying a nonthermal distribution) is studied by the reductive perturbation method. The modified Gardner (MG) equation is derived and numerically solved. It has been found that the basic characteristics of the EA Gardner solitons (GSs), which are shown to exist for α around its critical value α _c [where α is the nonthermal parameter, α _c is the value of α corresponding to the vanishing of the nonlinear coefficient of the Korteweg-de Vries (K-dV) equation, e.g. α _c ≃0.31 for μ=n _h0/n _i0=0.5, σ=T _h /T _i =10, n _h0, n _i0 are, respectively, hot electron and nonthermal ion number densities at equilibrium, T _h (T _i ) is the hot electron (ion) temperature], are different from those of the K-dV solitons, which do not exist for α around α _c , and mixed K-dV solitons, which are valid around α∼α _c , but do not have any corresponding double layers (DLs) solution. The parametric regimes for the existence of the DLs, which are found to be associated with positive potential, are obtained. The present investigations can be observed in various space plasma environments (viz. the geomagnetic tail, the auroral regions, the cusp of the terrestrial magnetosphere, etc.).