Transformations to extend the range of application of power series solutions of differential equations of motion

Let the solution of a differential equation, expanded in powers of the independent variablet, have radius of convergenceT. let , wheret=t(), be a new independent variable, and let the corresponding power series in have radius of convergenceS. Thent(S) will not in general be equal toT. Ift(S)>T, then the series in powers of may have advantages over those in powers oft. Mathematical consequences of this distinction have been appreciated since the time of Poincaré. In this note the practical applications of some transformations are investigated.     

A simplified model solar atmosphere

A simplified representation of the temperature distribution in the solar photosphere is proposed: ( ₀) = ₀ - ₁ log ₀. An expression is derived for the emergent continuous spectrum from the simple model. The limitations and applications of the simple model are discussed.     

A new time element for a general time transformation

The paper introduces a new time element to be used with a general time transformation for satellite equations of motion. The purpose of this time element is to reduce the growth of the numerical errors with respect to the time integration. It is characteristic for the new time element , that the relation between the timet and the element does not depend on the independent variables.     

On the A-transform of the exponential integrals

In this paper, a technique of recursive analysis is developed for the integral transform A of the exponential integral functionsE _n which is denoted as _n (). The main result of this analysis enables us to establish a two-term recurrence formula for _n (0) and a three-term recurrence formula for _n (); 0. A computational algorithm based on these formulae is also constructed and its numerical results forn=2(1)25 are presented to 15-digit accuracy.     

Optically thick accretion onto black holes

Spherically symmetric, steady-state, optically thick accretion onto a nonrotating black hole with the mass of is studied. The gas accreting onto the black hole is assumed to be a fully ionized hydrogen plasma withn ₀=10⁸ cm^–3 andT ₀=10⁴ K far from the black hole, and a new approximate expression for the Eddington factor is introduced. The luminosity is estimated to beL=1.875×10³³ erg s^–1, which primarily arises from the optical surface (1) ofT10⁴ K. The accretion flow is characterized by 1 and (v/c)10. In the optically thin region, the flow remains isothermal, and the increase of temperature occurs at 1. The radiative equilibrium is strictly realized at (v/c)10.     

Gravitational radiation and spiralling time of close binary systems (V)

Forty-six binary systems with their primary component masses between 2M and 3M have been considered for gravitational radiation study. Power output by gravitational radiation (P _B ) and spiral time ₀ for all individual systems have been evaluated. A relation has been given betweenP _B and ₀. The rate of decrease of orbital period (P) has also been given for 10 eccentric orbit systems.     

Theory of the Trojan asteroids

In the previously published Parts I and II of the paper, the author has constructed a formal long-periodic solution for the case of 11 resonance in the restricted problem of three bodies to 0(m ^3/2), wherem is the small mass parameter of the system. The time-dependencet(, ,m), where is the mean synodic longitude and is related to the Jacobi constant, has been expressed by ahyperelliptic integral. It is shown here that with the approximationm=0 in the integrand, the functiont(, , 0) can be expanded in a series involving standardelliptic functions. Then the problem of inversion can be formally solved, yielding the function (t, , 0).Similarly, the normalized period (,m) of the motion can be approximated by theHagihara hyperelliptic integral (, 0), corresponding tom=0. This integral is also expanded into elliptic functions. Asymptotic forms for (, 0) are derived for 0 and for 1, corresponding to the extreme members of thetadpole branch of the family of orbits.     

A novel methodology for radiative transfer in a planetary atmosphere

A novel methodology for evaluating the field of anisotropically scattered radiation within a homogeneous slab atmosphere of arbitrary optical thickness is provided. It departs from the traditional radiative transfer approach in first considering that the atmosphere is illuminated by an isotropic light source. From the solution of this problem, it subsequently proceeds to that for the more conventional case of monodirectional illumination. The azimuthal dependence of the field is separated in the usual manner by an harmonic expansion, leaving a problem in four dimensions (=optical depth, ₀=thickness, , =directions of incidence and scattering) which, as is well known, is numerically extremely inconvenient. Two auxiliary radiative transfer formulations of increasing dimensionality are considered: (i) a transfer equation for the newly introduced functionb ^m(,,₀) with Sobolev's function ^m(,₀) playing the role of a source-function. Because the incident direction does not intervene, ^m is simply expressed as a single integral term involvingb ^m. For bottom illumination, an analogous equation holds for the other new functionh ^m(,,₀). However, simple reciprocity relations link the two functions so that it is only necessary to considerb ^m; (ii) a transfer equation for the other new functiona ^m(,,,₀) with a source-function provided by Sobolev's functionD ^m(,,₀). For bottom illumination, another functionf ^m(,,,₀) is introduced; by a similar argument using reciprocity relations,f ^m is reduced toa ^m rendering necessary only the consideration ofa ^m. However, a fundamental decomposition formula is obtained which shows thata ^m is expressible algebraically in terms of functions of a single angular variable. The functions ^m andD ^m are shown to be the values in the horizontal plane ofb ^m anda ^m, respectively. The other auxiliary functionsX ^m andY ^m are also expressed algebraically in terms ofb ^m. These results enable one to proceed to the final step of evaluating the radiation field for monodirectional illumination. The above reductions toalgebraic relations involving only the functionb ^m appear to be more advantageous than Sobolev's (1972) recent approach; they also circumvent some basic numerical difficulties in it. We believe the present approach may likewise prove to be superior to most (if not all) other methods of solution known heretofore.This paper presents the results of one phase of research carried out at the Jet Propulsion Laboratory under Contract No. NAS-7-100 sponsored by the National Aeronautics and Space Administration.     

Effects of metal abundances on the evolutionary period changes of classical Cepheids

Some peculiarities in the behaviour of a model self-gravitating system described by hydrodynamical equations and isothermal equation of state connected with the presence of thermodynamical fluctuations in real systems were investigated in numerical experiment. The values of density and velocity , , respectively, were computed by numerical code perturbed on each time-step and in each computational cell by random values , for modeling such fluctuations. Perturbed values _i = _i + _i ,v _i = _i + v _i were used to initiate the next step of computations. This procedure is equivalent to an introduction into original hydrodynamical equations of Langevin sources which are random functions. It is shown that these small fluctuations (= v =0,² =v ² = 10^–8) grow many times in marginally-stable state.     

High order symplectic integrators for perturbed Hamiltonian systems

A family of symplectic integrators adapted for the integration of perturbed Hamiltonian systems of the form H=A+B was given in (McLachlan, 1995). We give here a constructive proof that for all integer p, such integrator exists, with only positive steps, and with a remainder of order O(^p + ²²), where is the stepsize of the integrator. Moreover, we compute the analytical expressions of the leading terms of the remainders at all orders. We show also that for a large class of systems, a corrector step can be performed such that the remainder becomes O(^p +⁴²). The performances of these integrators are compared for the simple pendulum and the planetary three-body problem of Sun–Jupiter–Saturn.     

Search for low energyν e,μ,τ andv e,μ,τ in coincidence with BATSE gamma ray bursts

A search for low energy neutrinos of all flavours in correlation with 553 ray bursts detected by BATSE aboard the Compton Observatory has been performed by the LSD (Liquid Scintillator Detector) neutrino telescope. No excess ofe,, orv _e,, candidate has been detected by LSD during the time interval in which BATSE detected the 90% of the photon flux for any of the GRBs analyzed. Upper limits on the neutrino fluxes are given in the paper.     

Dimensional Analysis of the Flare Distribution Problem

Dimensional analysis is used to derive the distribution of solar flare energies,p() = A-^3/2, in accordance with recent observational and numerical results. Several other scalings, notably _fl ² , where _fl is the flare duration, are obtained as well.     

On the solar oblateness: The combined effect of a pole-equator difference in effective temperature and mechanical heating

With the help of a model atmosphere of the Sun we evaluate the pole-equator difference in flux (as measured by Dicke and Goldenberg) assuming the following type of pole-equator temperature difference (T=T _e –T _p ): (a) T 2K for > ₀ (₀ 0.05); (b) T 10K for < ₀.The small T at all optical depths given by (a) could, for example, be due to a pole-equator difference in effective temperatures. At small optical depths a difference in mechanical heating could give rise to the larger temperature difference given by (b). We compare the results of our calculations with Dicke and Goldenberg's observations.The National Center for Atmospheric Research is sponsored by the National Science Foundation.     

Planetary accretions in jet streams

The problem of planetary accretion in a jet stream is studied using the model developed by Alfvén and Arrhenius. We find that there are basically three types of planetary accretion: namely, fast process _c < _i , slow process _c ~ _i and delayed process _c > _i where _c is the characteristic time of the occurrence of catastrophic accretion and _i the time-scale of mass injection to the planetary system (3×10⁸ yr). These different time scales of accretion are found to be closely related to the primordial thermal profiles and equatorial inclinations of the planets. Finally, Venus' retrograde rotational spin is shown to be a possible result of accretion process in a jet stream.     

О природе квазаров и активных ядер галактик

( — , , ). ( , ) . , .

---

---

, 40 48 MAC (, 1976).     

Vertical motions in an intense magnetic flux tube

The recent discovery of localised intense magnetic fields in the solar photosphere is one of the major surprises of the past few years. Here we consider the theoretical nature of small amplitude motions in such an intense magnetic flux tube, within which the field strength may reach 2 kG. We give a systematic derivation of the governing expansion equations for a vertical, slender tube, taking into account the dependence upon height of the buoyancy, compressibility and magnetic forces. Several special cases (e.g., the isothermal atmosphere) are considered as well as a more realistic, non-isothermal, solar atmosphere. The expansion procedure is shown to give good results in the special case of a uniform basic-state (in which gravity is negligible) and for which a more exact treatment is possible.The form of both pressure and velocity perturbations within the tube is discussed. The nature of pressure perturbations depends upon a critical transition frequency, _p , which in turn is dependent upon depth, field strength, pressure and density in the basic (unperturbed) state of the tube. At a given depth in the tube pressure oscillations are possible only for frequencies greater than _p for frequencies below _p exponentially decaying (evanescent) pressure modes occur. In a similar fashion the nature of motions within the flux tube depends upon a transition frequency, _v . At a given depth within the tube vertically propagating waves are possible only for frequencies greater than _v ; for frequencies below _v exponentially decaying (evanscent) motions occur.The dependence of both _v and _p on depth is determined for each of the special cases, and for a realistic solar atmosphere. It is found that the use of an isothermal atmosphere, instead of a more realistic temperature profile, may well give misleading results.For the solar atmosphere it is found that _v is zero at about 12 km above optical depth ₅₀₀₀= 1, thereafter rising to a maximum of 0.04 s^–1 at some 600 km above ₅₀₀₀ = 1. Below ₅₀₀₀ = 1, in the convection zone, _v has a maximum of 0.013 s^–1. The transition frequency, _p , for the pressure perturbations, is peaked at 0.1 s^–1 just below ₅₀₀₀ = 1, falling to a minimum of 0.02 s^–1 at about one scale-height deeper in the tube     

Physical Properties of Dust in the Martian Atmosphere: Analysis of Contradictions and Possible Ways of Their Resolution

Atmospheric aerosols play an important role in forming the Martian climate. However, the basic physical properties of the Martian aerosols are still poorly known; there are many contradictions in their estimates. We present an analytical overview of the published results and potentialities of various methods. We consider mineral dust. Zonally averaged data obtained from mapping IR instruments (TES and IRTM) give the optical thickness of mineral aerosols ₉ = 0.05–0.1 in the 9-m band for quite atmospheric conditions. There is a problem of comparing these estimates with those obtained in the visible spectral range. We suggest that the commonly used ratio _vis/₉ >2 depends on the interpretation and it may actually be smaller. The ratio _vis/₉ 1 is in better agreement with the IRIS data (materials like montmorillonite). If we assume that _vis/₉ = 1 and take into account the nonspherical particle shape, then the interpretation of ground-based integrated polarimetric observations ( < 0.04) can be reconciled with IR measurements from the orbit. However, for thin layers, the sensitivity of both methods to the optical thickness is poorly understood: on the one hand, polarimetry depends on the cloud cover and, on the other hand, the interpretation of IR measurements requires that the atmospheric temperature profile and the surface temperature and emissivity be precisely known. For quite atmospheric conditions, the local optical-thickness estimates obtained by the Bouguer–Lambert–Beer method and from the sky brightness measured from Viking 1 and 2 and Mars Pathfinder landers are much larger: = 0.3–0.6. Estimates of the contrasts in images from theVikingorbiters yield the same values. Thus, there is still a factor of 3 to 10 difference between different groups of optical-thickness estimates for the quiet atmosphere. This difference is probably explained by the contribution of condensation clouds and/or by local/time variations.     

Some elementary geometric aspects in extending the dimension of the space of instants

A local geometric construction is proposed on thepartially ordered set of instantsI. A totally ordered subsetC(I)I is assumed to have 3-dimensional affine coordinate structure,without a specified metric, called the -space ofC(I). Guided by a strong analogy with analytical mechanics theT-configuration space (, ), a real parameter, is constructed where-upon the usual Hamilton-Jacobi theory establishes a simple geometrical construction, viz., the complete figure from the calculus of variations. The duration function, dur:C(I)R is associated withtemporally equidistant hypersurfaces through which pass a congruence of extremal curves to the fundamental integral.     

On Geovenelli's approximate form of j(τ)

The accuracy of the approximate replacement of the frequency meanj() of the total intensityJ () byJ _c+a{J ₀()–J _c} by Geovenelli, whereJ _c is the total intensity at the continuum,J ₀() is that at the line center anda is independent of depth, is studied. Numerical calculations for strong solar lines show very high values of the residual intensities. A slight change in the approximation is then suggested to get a result quite consistent with observational and existing calculated data.     

Околозвездньiе облака поданньiм ультрафиолетовьiх спектров свезл

, ii (2000–3000 Å) i . , i . i (. 2). i i i i + ( 7–10). ii (. 13). ii i i (, 2400 Å) (. 14 15). i i i , iu , i (. 1). i i ii i i . .