共查询到20条相似文献，搜索用时 31 毫秒
1.
Rabindra Nath Das 《Astrophysics and Space Science》1980,71(1):25-35
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). 相似文献
2.
Byron D. Tapley 《Celestial Mechanics and Dynamical Astronomy》1970,2(3):319-333
A completely regular form for the differential equations governing the three-dimensional motion of a continuously thrusting space vehicle is obtained by using the Kustaanheimo-Stiefel regularization. The differential equations for the thrusting rocket are transformed using the K-S transformation and an optimal trajectory problem is posed in the transformed space. The canonical equations for the optimal motion in the transformed space are regularized by a suitable change of the independent variable. The transformed equations are regular in the sense that the differential equations do not possess terms with zero divisors when the motion encounters a gravitational force center. The resulting equations possess symmetry in form and the coefficients of the dependent variables are slowly varying quantities for a low-thrust space vehicle.Presented at the Conference on Celestial Mechanics, Oberwolfach, Germany, August 17–23, 1969. 相似文献
3.
It is shown that it is possible to make a change of variables in a Lagrangian in such a way that the number of variables is increased. The Euler-Lagrange equations in the redundant variables are obtained in the standard way (without the use of Lagrange Multipliers!). These equations are not independent but they are all valid and consistent. In some cases they are simpler than if the minimum number of variables are used. The redundant variables are supposed to be related to each other by several constraints (not necessarily holonomic), but these constraints are not used in the derivation of the equations of motion. The method is illustrated with the well known Kustaanheimo-Stiefel Regularization. Some interesting applications to perturbation theory are also described.The present research was carried out partially at the University of California and partially at the Jet Propulsion Laboratory under contract NAS7-100 with NASA.Presently Visiting Associate Professor at the University of Texas. 相似文献
4.
G. H. Born E. J. Christensen L. K. Seversike 《Celestial Mechanics and Dynamical Astronomy》1974,9(1):41-53
The concept of employing osculating reference position and velocity vectors in the numerical integration of the equations of motion of a satellite is examined. The choice of the reference point is shown to have a significant effect upon numerical efficiency and the class of trajectories described by the differential equations of motion. For example, when the position and velocity vectors on the osculating orbit at a fixed reference time are chosen, a universal formulation is yielded. For elliptical orbits, however, this formulation is unattractive for numerical integration purposes due to Poisson terms (mixed secular) appearing in the equations of motion. Other choices for the reference point eliminate this problem but usually at the expense of universality. A number of these formulations, including a universal one, are considered here. Comparisons of the numerical characteristics of these techniques with those of the Encke method are presented. 相似文献
5.
J. M. A. Danby 《Celestial Mechanics and Dynamical Astronomy》1970,2(3):311-318
Matrix methods for computing perturbations of non-linear perturbed systems, as formulated by Alexeev, involve an expression for the full solution of the first variational equations of the system evaluated about a reference orbit. These cannot be immediately applied to a regularized system of equations where perturbations about Keplerian motion are considered since the solution of the variational equations of regularized Keplerian motion does not in general correspond to the solution of the variational equations of the unregularized equations. But, as Kustaanheimo and Stiefel have pointed out, the regularized equations of Keplerian motion should be excellent for the initiation of a perturbation theory since they are linear in form. This paper describes a method for applying Alexeev's theorem to a regularized system where full advantage is taken of the basic linear form of the unperturbed equations.Presented at the Conference on Celestial Mechanics, Oberwolfach, Germany, August 17–23, 1969. 相似文献
6.
The possible existence of strong magnetic fields in stars is discussed and a method of constructing highly distorted models of magnetic, rotating stars developed. For stars with both poloidal and toroidal fields at the surface a force-free outer boundary condition is necessary. Non-linear solutions of the force-free equations must be used. The force-free equations and the structure equations for a white dwarf are solved simultaneously by a finite difference method.The National Center for Atmospheric Research is sponsored by the National Science Foundation. 相似文献
7.
Angelo J. Skalafuris 《Astrophysics and Space Science》1969,3(2):234-257
We examine the problem of a shock wave propagating in a gravitational field in the presence of pressure and density gradients by attacking the non-linear equations of fluid flow. Our approach is analytical rather than numerical, and we analyze the characteristic equations of a fluid in the presence of gravity with radiative dissipation. Because the radiation field enters the fluid equations in the form of an integral, radiative dissipation may be considered an inhomogeneity which does not affect the characteristic directions. The fluid equations remain hyperbolic and thus are amenable to solution by the standard techniques of gas analysis.We give an equation of path for a shock wave and we enumerate the physical conditions which lead to stability or instability. We find that shock waves are generally unstable in most stellar atmospheres unless they are very weak. The form of the instability is that of a spicule deformation similar to that observed in the upper solar chromosphere.This work was carried out at the Smithsonian-Harvard Astrophysical Observatory and was presented in a thesis to Brandeis University, May 1963. 相似文献
8.
Joseph J. F. Liu Philip M. Fitzpatrick 《Celestial Mechanics and Dynamical Astronomy》1975,12(4):463-487
Differential equations are derived for studying the effects of either conservative or nonconservative torques on the attitude motion of a tumbling triaxial rigid satellite. These equations, which are analogous to the Lagrange planetary equations for osculating elements, are then used to study the attitude motions of a rapidly spinning, triaxial, rigid satellite about its center of mass, which, in turn, is constrained to move in an elliptic orbit about an attracting point mass. The only torques considered are the gravity-gradient torques associated with an inverse-square field. The effects of oblateness of the central body on the orbit are included, in that, the apsidal line of the orbit is permitted to rotate at a constant rate while the orbital plane is permitted to precess (either posigrade or retrograde) at a constant rate with constant inclination.A method of averaging is used to obtain an intermediate set of averaged differential equations for the nonresonant, secular behavior of the osculating elements which describe the complete rotational motions of the body about its center of mass. The averaged differential equations are then integrated to obtain long-term secular solutions for the osculating elements. These solutions may be used to predict both the orientation of the body with respect to a nonrotating coordinate system and the motion of the rotational angular momentum about the center of mass. The complete development is valid to first order in (n/w
0)2, wheren is the satellite's orbital mean motion andw
0 its initial rotational angular speed. 相似文献
9.
We present a new numerical method for solving the system of partial differential equations describing the structure and evolution of a spherically symmetric star. As usual, we employ the transversal method of lines in order to split the equations into a coupled spatial and temporal part. The novel features of the algorithm are the following: (a) Instead of using the Lagrangian picture we formulate the system of partial differential equations in the Eulerian picture. (b) We reformulate the equations of stellar structure as a multipoint boundary-value problem. By means of this reformulation the rather clumsy iterative matching procedure of stellar atmosphere and interior is avoided. (c) The multipoint boundary-value problem is solved by the multiple shooting method. This approach not only ensures a high accuracy of the stellar models calculated at each time step but also allows the free boundaries inside the star due to different energy transport mechanisms to be located exactly. (d) The time derivatives involved in the stellar-structure equations are discretized implicitly to second order accuracy. Moreover, at each time step, the chemical abundances are determined by using a sophisticated update procedure. In this way, a high accuracy is achieved with respect to the integration in time. The algorithm has turned out to be exceedingly reliable and numerically accurate. This is shown by the evolution of a 1 M⊙ star up to the hydrogen-shell burning phase. In this example, the virial theorem, the law of mass conservation, and the law of energy conservation is fulfilled to a hitherto unattainable degree of accuracy. Since the multiple shooting method, which is at the heart of our approach, is a perfect example of a parallel algorithm, the computational speed of the algorithm might be substantially improved provided easy-to-program, high-performance parallel computers with sufficiently many processors become available in the near future. 相似文献
10.
M. S. Abdel Krim 《Astrophysics and Space Science》1990,164(1):69-77
Radiative transfer equation in a plane-parallel medium with isotropic boundary conditions for linearly anisotropic scattering phase function is considered. Two coupled integral equations for total density of radiation and total radiation flux are obtained. The Galerkin method is used to solve these equations. Numerical results for the radiative fluxes at the boundaries show that the Galerkin method yields accurate results compared well with other exact methods. 相似文献
11.
We have compared solutions obtained from the bi-Maxwellian based 16-moment transport equations with those obtained from the Maxwellian based 13-moment transport equations for conditions leading to the steady state, subsonic flow of a fully-ionized electron-proton plasma along geomagnetic field lines in the vicinity of the plasmapause. The bi-Maxwellian based equations can account for large temperature anisotropies and the flow of both parallel and perpendicular thermal energy, while the Maxwellian based equations account for small temperature anisotropies and only the total heat flow. Our comparison indicates that for Stable Auroral Red arc (SAR-arc) conditions leading to strong field-aligned heat flows (temperatures of 8000 K and temperature gradients of4K. km−1 at 1500 km), the bi-Maxwellian based equations predict a different thermal structure in the topside ionosphere than the less rigorous Maxwellian based equations. In particular, the bi-Maxwellian based equations predict proton and electron temperature anisotropies with T > T, while the Maxwellian based equations predict the opposite behavior for the same boundary conditions. This difference is related to the way in which the temperature anisotropies and heat flows are treated in the two formulations. For the bi-Maxwellian based equations, the inclusion of separate heat flows for parallel and perpendicular thermal energy allows for the development of a pronounced tail in both the electron and proton distribution functions, which leads to temperature anisotropies with T > T. For the Maxwellian based equations, on the other hand, the tail development is restricted because only the total heat flow is considered. Consequently, as the heat flows down, the presence of an increasing magnetic field acts to produce an anisotropy with T > T, and this process dominates tail formation for the Maxwellian based equations. 相似文献
12.
13.
Relativistic flows both in stellar interior and stellar atmospheres are obviously in action, hence, their treatment is of interest. The present communication aims at deriving by a different method the general equations that govern such flows.The derivation is based merely on the first two principles of thermodynamics and the inertia-energy equivalence principles; while the end results are the formulation of the energy tensor and the energy-flux and momentum equations. The assumptions of ideal flow and adiabatic behaviour of the medium adopted here are restrictive. 相似文献
14.
A detailed derivation of the effect of solar radiation pressure on the orbit of a body about a primary orbiting the Sun is
given. The result is a set of secular equations that can be used for long-term predictions of changes in the orbit. Solar
radiation pressure is modeled as a Fourier series in the body’s rotation state, where the coefficients are based on the shape
and radiation properties of the body as parameters. In this work, the assumption is made that the body is in a synchronous
orbit about the primary and rotates at a constant rate. This model is used to write explicit variational equations of the
energy, eccentricity vector, and angular momentum vector for an orbiting body. Given that the effect of the solar radiation
pressure and the orbit are periodic functions, they are readily averaged over an orbit. Furthermore, the equations can be
averaged again over the orbit of the primary about the Sun to give secular equations for long-term prediction. This methodology
is applied to both circular and elliptical orbits, and the full equations for secular changes to the orbit in both cases are
presented. These results can be applied to natural systems, such as the binary asteroid system 1999 KW4, to predict their
evolution due to the Binary YORP effect, or to artificial Earth orbiting, nadir-pointing satellites to enable more precise
models for their orbital evolution. 相似文献
15.
There are two parts to this paper. In the first we calculate the hydrodynamic response of the solar atmosphere to the injection of an intense beam of electrons in a numerical simulation of a solar flare. In the second we predict the spectroscopic consequences of the hydrodynamic behaviour calculated in the first part. The hydrodynamics is predicted by solving the equations of conservation of mass, momentum, and energy. The latter is expressed as two temperature equations; one for the electrons and the other for the neutral atoms and positive ions of hydrogen. The equations are solved in one dimension and the geometric form is of a semi-circular loop having its ends in the photosphere. The results show how the loop is filled at supersonic speed with plasma at temperatures characteristic of flares. At the same time a compression wave is predicted to propagate down towards the photosphere. After the heating pulse stops, the plasma that has risen into the loop, starts to decay and return to the condition it was in before the pulse started. In predicting the spectrum that would be emitted by such a plasma calcium was chosen for illustration. The first and main part of this calculation was setting up and solving the time-dependent equations of ionization/recombination. In order to provide a standard for comparison the same ionization and recombination rate coefficients are used to predict the steady-state distribution of populations of ionization stages. This is then compared with the distribution found from the time-dependent solution and shows that there is a negligibly small time lag predicted by the time-dependent result. However the more significant comparisons to make are between the temperatures of the peak abundances of the various ions under the assumptions of steady-state and time-dependent ionization. For the particular circumstances chosen here the temperature differences are predicted to be in the neighbourhood of 10% or less and in view of the overall accuracy of the atomic data are not significant. It would appear therefore that the much simpler assumption of steady-state ionization balance leads to results of acceptable accuracy for the particular case considered. 相似文献
16.
The general-relativistic equations of hydrodynamics for a nonideal fluid are derived in the modified Lagrangian form. Together with the zeroth and first moment equations of radiative transfer derived by Morita and Kaneko (1986), the equations provide a complete set of the modified Lagrangian equations of radiation hydrodynamics. The equations of hydrodynamics are specialized for a thermally conducting, Newtonian viscous fluid in the modified Lagrangian form, which are the generalization of the specialrelativistic equations of hydrodynamics derived by Greenberg (1975). 相似文献
17.
《Chinese Astronomy and Astrophysics》1982,6(1):60-63
We discuss the excitation of polar motion by earthquake displacement field. Instead of the usual static equilibrium equations in the literature, we use an improved set as given in /1/, which guarantee continuity at the core-mantle boundary. We take the parameter values of three earthquakes from /2/.To obviate the singularity at r = 0, we use asymptotic solutions by power series within a small sphere around the centre. Outisde this sphere, the equations are numerically integrated by the Runge-Kutta algorithm. Our equations /1/ gave polar shifts some 3 times larger than Dahlen's equations /2/. 相似文献
18.
Non-linear equations governing the temporal evolution of the vector of instantaneous rotation are developed for an Earth with a homogeneous mantle having a viscoelastic Maxwell rheology and with a homogeneous inviscid fluid core.This general theory is investigated using the angular momentum theorem applied to the coupled core-mantle system. It allows to study the influence upon the planetary rotation of a quasi-rigid rotational motion in the liquid core. It also enables to investigate the consequences of excitation sources (e.g. pressure), located at the core-mantle interface. Especially, the influence of viscoelastic variations in the inertia tensors resulting from the rotation itself or from various excitation sources are detailed with the help of a Love number formalism. The equations of the linear theory for an elastic Earth with a liquid core, and the non-linear theory for a viscous planet with a quasi-fluid behavior are shown to be particular cases of our generalized system of equations. Some planetological applications may be derived from the quasi-fluid approximation. 相似文献
19.
Ryan E. Sherrill Andrew J. Sinclair S. C. Sinha T. Alan Lovell 《Celestial Mechanics and Dynamical Astronomy》2014,119(1):55-73
The relative motion of chief and deputy satellites in close proximity with orbits of arbitrary eccentricity can be approximated by linearized time-periodic equations of motion. The linear time-invariant Hill–Clohessy–Wiltshire equations are typically derived from these equations by assuming the chief satellite is in a circular orbit. Two Lyapunov–Floquet transformations and an integral-preserving transformation are here presented which relate the linearized time-varying equations of relative motion to the Hill–Clohessy–Wiltshire equations in a one-to-one manner through time-varying coordinate transformations. These transformations allow the Hill–Clohessy–Wiltshire equations to describe the linearized relative motion for elliptic chief satellites. 相似文献
20.
Victor R. Bond 《Celestial Mechanics and Dynamical Astronomy》1974,10(3):303-318
The Newtonian differential equations of motion for the two-body problem can be transformed into four, linear, harmonic oscillator equations by simultaneously applying the regularizing time transformation dt/ds=r and the Kustaanheimo-Stiefel (KS) coordinate transformation. The time transformation changes the independent variable from time to a new variables, and the KS transformation transforms the position and velocity vectors from Cartesian space into a four-dimensional space. This paper presents the derivation of uniform, regular equations for the perturbed twobody problem in the four-dimensional space. The variation of parameters technique is used to develop expressions for the derivatives of ten elements (which are constants in the unperturbed motion) for the general case that includes both perturbations which can arise from a potential and perturbations which cannot be derived from a potential. These element differential equations are slightly modified by introducing two additional elements for the time to further improve long term stability of numerical integration.Originally presented at the AAS/AIAA Astrodynamics Specialists Conference, Vail, Colorado, July 1973 相似文献