首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 62 毫秒
1.
Starting with the Einstein-Maxwell field equations in general relativity we construct the general differential equations governing the components of the metric tensor. We do this in a fashion completely analogous to that which we follow in Einstein-Infeld-Hoffmann and Einstein-Infeld theory. These equations allow us to findh if in various orders. An answer to our problem up to the first relativistic corrections is a computational work to find4 h 00, since the other terms2 h 00,2 h and3 h 0 are as in a pure gravitational case. On the other hand, using the defined Einstein-Maxwell tensor, we give the equations of motion of two charged particles in the 0th order; also, the generalization is given in the case ofn particles.  相似文献   

2.
In the present paper n 0 , 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 0 with respect to the modulusX orC, and Section 3 includes the numerical method of the solutions of these differential equations. Theoretical values of 0 0 and 1 0 , with corresponding values ofC, are also added in this section.  相似文献   

3.
The direct problem of dynamics in two dimensions is modeled by a nonlinear second-order partial differential equation, which is therefore difficult to be solved. The task may be made easier by adding some constraints on the unknown function = f y /f x , where f(x, y) = c is the monoparametric family of orbits traced in the xy Cartesian plane by a material point of unit mass, under the action of a given potential V(x, y). If the function is supposed to verify a linear first-order partial differential equation, for potentials V satisfying a differential condition, can be found as a common solution of certain polynomial equations.The various situations which can appear are discussed and are then illustrated by some examples, for which the energy on the members of the family, as well as the region where the motion takes place, are determined. One example is dedicated to a Hénon—Heiles type potential, while another one gives rise to families of isothermal curves (a special case of orthogonal families). The connection between the inverse/direct problem of dynamics and the possibility of detecting integrability of a given potential is briefly discussed.This revised version was published online in October 2005 with corrections to the Cover Date.  相似文献   

4.
A similarity analysis for the free and forced convection hydromagnetic flow over a horizontal semi-infinite flat plate through a non-homogeneous porous medium is presented, taking into account the hydrostatic pressure variation normal to the flat plate. The similarity solution of the problem under consideration is obtained under certain valid simplifying assumptions when, (i) the plate temperature is inversely proportional to the square root of the distance from the leading edge, (ii) the intensity of the applied magnetic field, normal to the plate, changes with the inverse square root of the distance from the leading edge, and (iii) the permeability of the porous medium, occupying a semi-infinite region of the space bounded by the flat plate, is proportional to the distance measured in the direction of the flow. A numerical solution of the resulting system of ordinary differential equations of motion and energy is obtained, depending on the Prandtl number Pr, the magnetic parameterM n ,the bouyancy parameter , and the permeability parameterP m .The variations of the fundamental quantities of the problem are shown graphically followed by a quantitative discussion.  相似文献   

5.
The aim of the present paper will be to detail the explicit form of the equations which govern first-order oscillations of fast-rotating globes of self-gravitating fluids; with due account taken of the effects arising from the centrifugal as well as Coriolis force. As such configurations oscillate in general about distorted figures of equilibrium, the equations governing them can be conveniently expressed in terms of the Clairaut coordinates, associated with distorted spheroidal figures, and introduced in our previous paper (Kopal, 1980) for this purpose.In Section 2 which follows a brief outline of our problem, the equilibrium properties of fast-rotating configurations or arbitrary structure will be formulated. In Section 3 we shall carry out a separation of the variables in the equations of motion, and reduce the partial differential equations of the problem to an equivalent system of ordinary differential equations, by an expansion of expressions for the velocity componentsU, V, W in terms of tesseral harmonicsY n m (, ). The explicit form of such a system, including the effects of all tesseral harmonics of orders up tom=n=4, will be specified in Section 3 for configurations whose equilibrium form is a sphere; while in Section 4 this latter condition will be relaxed to allow for the equilibrium configuration to become a rotational spheroid.In the concluding Section 5 we shall convert the complex form of our equations of motion into real terms, amenable to a solution-analytical or numerical-in terms of real variables; and shall establish the boundary conditions necessary for a specification of the characteristic frequencies of oscillation.  相似文献   

6.
The electromagnetic field produced by a magnetic dipole moment, , which is rotating obliquely surrounded by a corotating plasma sphere, is investigated. This corotating-plasma approximation has the same order of accuracy as the force-free one but has somewhat different physical implications. In the former the effect of non-electromagnetic forces such as the inertial force are included, though in somewhat artificial manner, as a departure from the strict MHD condition and this fact seems to guarantee the existence of physical solutions.Analogous to the relativistic force-free equation, a set of two differential equations (the corotation equation) are derived for the scalar functions associated with the electric and magnetic fields. A self-consistent solution of these equations is given and it is shown that this solution has no singularity, in spite of apparent divergence in the formal solution, on the light cylinder. It is concluded from this solution that, even in the extreme case of the largest possible corotation-radius (i.e.b=r L , wherer L is the light radius), the existence of a corotating plasma does not alter the field structure drastically from the vacuum case. It is also suggested through this treatment that inclusion of the inertial term in generalized Ohm's law might be essential in considering the centrifugal-wind problem.  相似文献   

7.
A recent least squares algorithm, which is designed to adapt implicit models to given sets of data, especially models given by differential equations or dynamical systems, is reviewed and used to fit the Hénon-Heiles differential equations to chaotic data sets.This numerical approach for estimating parameters in differential equation models, called theboundary value problem approach, is based on discretizing the differential equations like a boundary value problem,e.g. by a multiple shooting or collocation method, and solving the resulting constrained least squares problem with a structure exploiting generalized Gauss-Newton-Method (Bock, 1981).Dynamical systems like the Hénon-Heiles system which can have initial values and parameters that lead to positive Lyapunov exponents or phase space filling Poincaré maps give rise to chaotic time series. Various scenarios representing ideal and noisy data generated from the Hénon-Heiles system in the chaotic region are analyzedw.r.t. initial conditions, parameters and Lyapunov exponents. The original initial conditions and parameters are recovered with a given accuracy. The Lyapunov spectrum is then computed directly from the identified differential equations and compared to the spectrum of the true dynamics.presently at IWR, Universität Heidelberg, Im Neuenheimer Feld 368, D-6900 Heidelberg, Germany  相似文献   

8.
Earlier models of compressible, rotating, and homogeneous ellipsoids with gas pressure are generalized to include the presence of radiation pressure. Under the assumptions of a linear velocity field of the fluid and a bounded ellipsoidal surface, the dynamical behaviour of these models can be described by ordinary differential equations. These equations are used to study the finite oscillations of massive radiative models with masses 10M and 30M in which the effects of radiation pressure are expected to be important.Models with two different degrees of equilibrium are chosen: an equilibrium (i.e., dynamically stable) model with an initial asymmetric inward velocity, and a nonequilibrium model with a nonequilibrium central temperature and which falls inwards from rest. For each of these two degrees of equilibrium, two initial configurations are considered: rotating spheroidal and nonrotating spherical models.From the numerical integration of the differential equations for these models, we obtain the time evolution of their principal semi-diametersa 1 anda 3, and of their central temperatures, which are graphically displayed by making plots of the trajectories in the (a 1,a 3) phase space, and of botha 1 and the total central pressureP c against time.It is found that in all the equilibrium radiative models (in which radiation pressure is taken into account), the periods of the oscillations of botha 1 andP c are longer than those of the corresponding nonradiative models, while the reverse is true for the nonequilibrium radiative models. The envelopes of thea 1 oscillations of the equilibrium radiative models also have much longer periods; this result also holds for the nonequilibrium models whenever the envelope is well defined. Further, as compared to the nonradiative models, almost all the radiative models collapse to smaller volumes before rebouncing, with the more massive model undergoing a larger collapse and attaining a correspondingly larger peakP c.When the mass is increased, the dynamical behavior of the radiative model generally becomes more nonperiodic. The ratio of the central radiation pressure to the central gas pressure, which is small for low mass models, increases with mass, and at the center of the more massive model, the radiation pressure can be comparable in magnitude to the gas pressure. In all the radiative models, the average periods as well as the average amplitudes of both thea 1 andP c oscillations also increase with mass.When either rotation or radiation pressure effects or both are included in the equilibrium nonradiative model, the period of the envelope of thea 1 oscillations is increased. The presence of rotation in the equilibrium radiative model, however, decreases this period.Some astrophysical implications of this work are briefly discussed.  相似文献   

9.
Roxburgh  Ian W. 《Solar physics》1974,35(2):481-487
The solution curves of the differential equations determining the behavior of the solar wind are calculated for the case where the heat flux has its maximum value 3/2 nkTv th. All the supersonic solutions are asymptotically adiabatic, T r -4/3.The National Center for Atmospheric Research is sponsored by the National Science Foundation.  相似文献   

10.
In this paper we propose a method for computing the equilibrium structure of differentially rotating polytropic models of the stars. A general law of differential rotation of the type 2=b 0+b 1 s 2+b 2 s 4, which can account for a reasonably large variety of possible differential rotations in the stars has been used. The distortional effects have been incorporated in the structure equations up to second order of smallness in distortion parametersb 0,b 1, andb 2 using Kippenhahn and Thomas' averaging approach in conjunction with Kopal's results on Roche equipotentials in manner similar to the one earlier used by Mohan and Saxena for computing the equilibrium structure of polytropes having solid body rotation. Numerical results have been obtained for various types of differentially rotating polytropic models of stars of polytropic indices 1.5, 3, and 4. Certain differentially rotating models of the Sun which are possible with such a type of law of differential rotation, have also been computed.  相似文献   

11.
The problem on linear waves in a radiating and scattering grey medium is studied using Whitham's method. Analysis of the basic equations distinguishes two limiting cases: the one is theequilibrium case in which the energy exchange between the gas and radiation plays an essential role, and the other is theScattering case in which the effect of energy exchange is negligible. A new type ofradiation acoustic wave with the speed is found in the scattering case. The governing equations for linearized one-dimensional flow are reduced to one equation of radiative acoustics valid to order 1/c, and the criterion for the two limiting cases is derived from studying this equation. The harmonic solution is analytically studied to show that theeffective optical depth corresponding to the wavelength of perturbation gives the measure of the interaction between the gas and radiation. When eff1, the sound speeda g 2 =P g / and the propagating speed of radiative disturbancea f 2 =fc 2 appear as the modified classical and radiation-induced modes respectively, wheref is the Eddington factor. When eff1, the isentropic sound speeda s 2 =(P g +P r / appears in the equilibrium case, and the radiation acoustic speeda A 2 appears in the scattering case. The dispersion relation of the harmonic solution is numerically calculated. The result shows that the wave pattern depends critically on the ratio=P g /(P g +P r ). When , the speeda S anda A arise from the modified classical mode, and when , they originate from the radiation-induced mode.  相似文献   

12.
In this paper the unsteady flow in the Ekman layer of a visco-elastic non-Newtonian fluid near a flat plate is discussed. Laplace transform technique has been employed to show the basic differential equations. Expressions for velocity profile, the skin friction have been calculated. It is shown that the time to attain the steady state increases with the elastic parameter. It is shown that normally the ultimate steady state is reached through a decay of inertial oscillations whose frequency decreases with increase in the elastic parameter. In the present study we examine the following unsteady problem in non-Newtonian fluid. Consider an infinite plate coinciding with the platez=0 and rotating in unison with elasticoviscous liquid occupying the regionz>0 with a uniform angular velocity about thez-axis for timet<-0. At timet>0, the plate starts moving with a uniform velocityU o along thex-axis relative to the rotating frame of reference. The horizontal homogeneity of the problem demands that conditions depend onz andt only. The equation of continuity together with the no slip condition at the plate then shows that thez-component of the velocity vanishes everywhere.  相似文献   

13.
The KS-transformation introduced by P. Kustaanheimo and E. Stiefel into celestial mechanics is derived straight from the Kepler Formulas. There follow the treatment of the inverse Newton problem comprising the derivation of the differential equations of mechanics by J. Hermann and L. Euler and also remarks concerning the fundamental papers by Euler about the planet problem and then-body problem. The conclusion is a simple example given by A. Voss and H. Liebmann, for the differential equations of mechanics with non-holonomic condition, which is of pseudoplanetary quality.This paper originated from several stays at the Eidgenössische Technische Hochschule, Zürich; Seminar of Professor E. Stiefel in 1973/74.  相似文献   

14.
In the ordinary restricted problem of three bodies, the first-order stability of planar periodic orbits may be determined by means of their characteristic exponents, as derived from the condition of a vanishing determinant for the coefficients of an infinite system of homogenous linear equations associated with the exponential series solutionu, v representing any initially small oscillations about the periodic solutionx, y. In the elliptic restricted problem, periodic solutions are possible only for periods which are equal to, or integral multiples of, the periodP of the elliptic motion of the two primary masses. It is shown that the infinite determinant approach to the determination of the characteristic exponents can be extended to the treatment of superposed free oscillations in the elliptic problem, and that in generaltwo exponents appear in any complete solutionu, v for eachone existing in the corresponding ordinary restricted problem. The value of each exponent depends on a series proceeding in even powers of the eccentricitye of the relative orbit of the two primaries, in addition to its basic dependence on the mass ratio . For stable periodic orbits, the oscillation frequenciesn 1 (,e 2),n 2 (,e 2) associated with these two exponents tend, withe0, to certain limiting valuesn 1 (),n 2(), which differ from each other by the amount of the frequencyN=2/P of the orbital motion of the primaries. One of the two frequencies, sayn 1(), is identical with the frequency of the corresponding oscillations in the ordinary restricted problem, while the second one gives rise to oscillations only in the elliptic restricted problem, withe0.The method will be described in more detail, together with its application to two families of small periodie librations about the equilateral points of the elliptic restricted problem (E. Rabe: Two new Classes of Periodic Trojan Librations in the Elliptic Restricted Problem and their Stabilities) in theProceedings of the Symposium on Periodic Orbits, Stability and Resonances, held at the University of São Paulo, Brasil, 4–12 September, 1969.Presented at the Conference on Celestial Mechanics, Oberwolfach, Germany, August 17–23, 1969.  相似文献   

15.
The paper deals with the excitation of the helium singlet level 21 P in the homogeneous and filamentary models of quiescent prominences with following parameters: the optical thickness at the limit of helium Lyman continuum 1c M = 0.1–100, T e = 7000 K, n e = 5 × 1010 cm–3. Assuming a model He atom with seven discrete levels (11 S, 23 S, 21 S, 23 P, 21 P, 33 D, 31 D) and the continuum the steady state equations for the levels 23 S, 21 P and the continuum have been solved together with the radiative transfer equations for the line 584 Å and the continuum 504 Å. The variations with depth of the functions n 2 3 S /n 1 1 S (1 c), n 2 1 P /n 1 1 S (1c ), and n + He n e /n 1 1 S(1c ) as well as the intensities of the triplet (D3, 10830 Å) and singlet (16678, 20581 Å) lines have been calculated. Comparison with observations leads to the following conclusions: (1) The line intensities calculated for filamentary models of prominences agree better with observations than those for homogeneous ones. (2) The helium level 21 P is excited by diffuse field 584 Å being formed by recombinations and spontaneous transitions 21 P – 11 S and escaping from the prominence into the space between the filaments and to the surface. (3) Underpopulation of the singlet level 21 P may be explained by combination of weak excitation mechanism (recombinations and formation of the diffuse field 584 Å) and strong deexcitation mechanism (spontaneous transitions into the level 11 S).  相似文献   

16.
Aimed at the initial value problem of the particular second-order ordinary differential equations,y =f(x, y), the symmetric methods (Quinlan and Tremaine, 1990) and our methods (Xu and Zhang, 1994) have been compared in detail by integrating the artificial earth satellite orbits in this paper. In the end, we point out clearly that the integral accuracy of numerical integration of the satellite orbits by applying our methods is obviously higher than that by applying the same order formula of the symmetric methods when the integration time-interval is not greater than 12000 periods.  相似文献   

17.
An idea is developed that the vacuum in the gravitational field acquires properties of an elastic medium described by a definite tension ik . The vacuum is stated to also participate in the formation of the space-time metric, together with the usual matter. So, the matter, vacuum and metric form a complex unity determined by the solution of the field equations. The vacuum may prove to play an essential role in the extremely strong fields existing in superdense celestial bodies. The tensor ik is not to be identified with the pseudo-tensor of the energy-momentum of the gravitational field the idea of which is preserved.The problem of vacuum is investigated in the case of the central symmetry static field. A number of properties of the tensor ik is found using the symmetry of the field and comparison with the post-Newton limit. The external and internal problems, as well as the procedure of joining the solutions on the surface of a celestial body, have been formulated. The stellar surface is determined in the usual way:P(r) = 0 whereP is the matter pressure. The theory includes three dimensionless parametersa=p/,b=p / (,p, p are the density of the vacuum energy and of its pressures in the radial and transverse directions) and determining the vacuum elastic properties. Generally speaking, they depend on the valueP/c2 in the stellar centre where is the mass density. From general physical considerations it is shown that 0 1 + lim P (l/q). The field equations are solved for the simple version of the theoryb=–a. There are solutions corresponding to superdense celestial bodies with masses considerably exceeding that of the Sun.  相似文献   

18.
This work contains a transformation of Hill-Brown differential equations for the coordinates of the satellite to a type which can be integrated in a literal form using an analytical programming language. The differential equation for the parallax of the satellite is also established. Its use facilitates the computation of Hill's periodic intermediary orbit of the satellite and provides a good check for the expansion of the coordinates and frequencies. The knowledge of the expansion of the parallax facilitates the formation of differential equations for terms with a given characteristic. These differential equations are put into a form which favors the solution by means of iteration on the computer. As in the classical theory we obtain the expansions of the coordinates and of the parallax in the form of trigonometric series in four arguments and in powers of the constants of integration. We expand the differential operators into series in squares of the constants of integration. Only the terms of order zero in these expansions are employed in the integration of the differential equations. The remaining terms are responsible for producing the cross-effects between the perturbations of different order. By applying the averaging operator to the right sides of the differential equations we deduce the expansion of the frequencies in powers of squares of the constants of integration.Basic Notations f the gravitational constant - E the mass of the planet - M the mass of the satellite - t dynamical time - x, y, z planetocentric coordinates of the satellite - u x+y–1 - s x–y–1 - the planetocentric distance of the satellite - w 1/ - 0 the variational part of - w 0 the variational part ofw, - n the mean daily sidereal motion of the satellite - a the mean semi-major axis of the satellite defined by means of the Kepler relation:a 3 n 2=f(E+M) - a the mean semi-major axis defined as the constant factor attached to the variational solution - e the constant of the eccentricity of the satellite - the sine of one half the orbital inclination of the satellite relative to the orbit of the sun - c(n–n) the anomalistic frequency of the satellite - c 0 the part ofc independent frome,e, and - g(n–n) the draconitic frequency of the satellite, - g 0 the part ofg independent frome,e, and - exp (n–n)t–1 - D d/d - e the eccentricity of the solar planetocentric orbit - a the semi-major axis of the solar orbit - n the mean daily motion of the sun in its orbit around the planet - m n/(n–n) - a/a-the parallactic factor - the disturbing function  相似文献   

19.
The stellar equilibrium equations for given surface pressureP * and temperatureT *, and in the absence of convection, are translated into a nonlinear integral equation, in which the radiusR of the star enters as an eigenvalue. We show that under broad mathematical assumptions on the constitutive equations (equation of state, opacity and energy generation) a global existence and uniqueness property can be formulated. If a valueP M is selected, which restricts the allowed pressure and temperature range |P(r)–P *|+E|T(r)–T *P M (E, arbitrary constant of dimensions of a pressure over temperature), thenat least one solutionP(r),T(r) exists in the pressure-temperature range chosen, for anyR<R E . This solution isunique forR<R c .R E andR c are expressed in terms of the constitutive equations, and of the pressure-temperature range adopted. A physical argument in favour of the stability of this solution is presented.  相似文献   

20.
The perturbed motion of a rigid body about its center of mass, is formulated in terms of the six elements:l, the magnitude of the angular momentum vector;h, the total energy; and , two linear functions of the independent variable; and 1 and 1, two Euler angles that orientate the inertial frame with respect to the unperturbed solution. Solutions from the element formulation and the original Euler equations are numerically compared using shuttle-type data. For applied torques smaller than a given magnitude, the element formulation produced the following results: (1) larger step sizes in the numerical integration of the differential equations, resulting in an overall computational time-saving, and (2) more significant figures of accuracy in the computation of the variables describing the state of the rigid body.  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号