A simple model of a supernova

In this paper we formulate the problem of the collapse of a spherically-symmetric, radiating body in general relativity. The requirement that the metric and its normal derivative be continuous across the boundary imposes conditions upon the evolution of the star and allows identification of physical phenomena measured by a distant observer. A solution to Einstein's field equations for the exterior of a spherically-symmetric radiating body is that derived originally by Vaidya in 1951. By requiring the continuity described above we identify the mass, luminosity, velocity, and time increment measured by a distant observer in terms of the metric parameters evaluated in a frame comoving with the outer boundary. We also assume that the interior metric is a sum of products of functions of the radius and time. The continuity requirements allow the evolution of two of the three functions of time to be determined. The evolution of the third function, describing the motion of the core, is determined by the imposition of an equation of state at the center. The adiabatic index derived from the Baym-Bethe-Pethick equation of state was used to provide this last equation. A major result is obtaining an analytic solution to Einstein's field equations describing the core of a collapsing star. As a consequence of this solution we found that for the relatively small values of the adiabatic index (_max1.6), the star smoothly made the transition to a final collapsed state. Neither bounce nor shock wave was obtained. Also, there is a readily understood connection between the adiabatic index, and such parameters of the edge of the core as the velocity and acceleration. Finally, the analytic solutions provide the time-scales for the collapse which are significantly different from that of free-fall. The retarding effects of pressure upon the collapse are apparent. It is hoped that such analytic solutions will provide insight into more complicated dynamic systems in general relativity.     

A continuum model for the orbit evolution of self-propelled ‘smart dust’ swarms

A continuity equation is developed to model the evolution of a swarm of self-propelled ‘smart dust’ devices in heliocentric orbit driven by solar radiation pressure. These devices are assumed to be MEMs-scale (micro-electromechanical systems) with a large area-to-mass ratio. For large numbers of devices it will be assumed that a continuum approximation can be used to model their orbit evolution. The families of closed-form solutions to the resulting swarm continuity equation then represent the evolution of the number density of devices as a function of both position and time from a set of initial data. Forcing terms are also considered which model swarm sources and sinks (device deposition and device failure). The closed-form solutions presented for the swarm number density provide insights into the behaviour of swarms of self-propelled ‘smart dust’ devices an can form the basis of more complex mission design methodologies.     

Exact and numerical solutions of the kompaneets equation: Evolution of the spectrum and average frequencies

The evolution of the spectrum of isotropic uniform radiation in an infinite space filled with a homogeneous, nonrelativistic electron gas is calculated by solving the Kompaneets equation. For an infinitely narrow initial spectrum, the time dependence of the average frequency and frequency dispersion is determined in a linear approximation of the equation. Characteristic times corresponding to changes in the character of this dependence are introduced. Two schemes are proposed for the numerical solution of the nonlinear equation: a nonconservative scheme with a grid that is uniform in frequency and a conservative scheme with automatic selection of an adaptive grid in frequency and time. For the linear equation the method yields results consistent with calculations of its solutions in terms of an eigenfunction expansion of the Kompaneets operator calculated in [D. I. Nagirner and V. M. Loskutov, Astrofizika, 40, 97 (1977)]. The influence of nonlinearity on the evolution of the spectrum of initially monochromatic radiation of various intensities is traced as an example of the application of the method.     

On the 'coarse-grained' evolution of collisionless stellar systems

We have suggested in a previous article that the coarse-grained evolution of a collisionless stellar system could be viewed as a diffusion process in velocity space compensated by an appropriate friction. Using a quasi-linear theory, we calculate the diffusion coefficient associated with this evolution. This provides a new self-consistent relaxation equation for f , the locally averaged distribution function. This equation bears some resemblance to the conventional Fokker–Planck equation of collisional systems but the friction term is non-linear in f (accounting for degeneracy effects) and the relaxation time is much smaller (in agreement with the concept of 'violent relaxation'). Under the condition that the diffusion current vanishes identically at equilibrium, we recover Lynden-Bell's distribution function; but if we allow stars to escape from the system at a constant rate, we can derive a truncated model which coincides with Lynden-Bell's solution in the core but provides a depletion of high-energy stars in the halo. This distribution function has a finite mass and is the generalization of the Michie–King model to the case of (possibly degenerate) collisionless stellar systems.     

Behaviour of viscous fluid in string cosmological models in the framework of Lyra geometry

In this communication, we studied the aspects of bulk viscous fluid cosmological model with quadratic equation of state in the presence of strings loaded with particles in a higher dimensional (5- dimensional) Bianchi type-III geometry in Lyra’s Manifold (Lyra, 1951). Using physically plausible circumstances, an exact model of the universe is presented by obtaining the solutions of the Einstein’s field equations. Important geometrical and dynamical parameters of the model universes are premeditated and physical significance regarding their prospect in modern cosmology are discussed in details. Interestingly it is seen that both bulk viscosity and quadratic equation of state are acting crucial jobs throughout the evolution of the model which is expanding with acceleration so it represents dark energy model universe. Hence our model can be thought as a realistic universe.     

Evolution of the cometary cloud under the action of stellar perturbations

The probability of variation of the integrals of the orbit as a result of an encounter was found for a two dimensional system. A method of solution of the Kolmogorov-Feller's equation is obtained using this probability function as a kernel, and it allows us to obtain the distribution of the integrals of the orbit as a function of time. The method is applied to the investigation of the evolution of orbits in the outer cometary cloud under the action of galactic stars. We consider the variations of orbits as a purely discontinuous random process, so we take into account not only distant but also close interactions.     

Evolution of the cometary cloud under the action of stellar perturbations

The probability of variation of the integrals of the orbit as a result of an encounter was found for a two dimensional system. A method of solution of the Kolmogorov-Feller's equation is obtained using this probability function as a kernel, and it allows us to obtain the distribution of the integrals of the orbit as a function of time. The method is applied to the investigation of the evolution of orbits in the outer cometary cloud under the action of galactic stars. We consider the variations of orbits as a purely discontinuous random process, so we take into account not only distant but also close interactions.     

Effects of electromagnetic field on the stability of locally isotropic gravastars

We have studied a new solution of charged gravastars with isotropic matter configuration in the framework of f(R, T) theory of gravity. For this purpose, we have assumed the electric charge as a constant. This stellar structure divided into three different regions: The preliminary part shows the interior charged region in which pressure equals to the negative density, second is the intermediate charged shell which is assumed to be very thin and filled with ultrarelativistic stiff fluid and the last corresponds to the electrovacuum region which is defined by an exterior Reissner-Nordström solution. Under these assumptions, we have found some physical aspects like length, energy, entropy and equation of state for charged spherical gravastar distribution. Moreover, we present an exact solution that free from event horizon and non-singular for this our new model.     

Virial oscillations of celestial bodies

A general approach to the solution of the perturbed oscillation problem for celestial bodies is considered. The solution sought describes unperturbed virial oscillations (zero approximation) affected by external perturbing effects. In the general case, these perturbations can be expressed by an arbitrary given function of time, Jacobi's function and its first derivative. Standard methods and modes of perturbation theory are used for solution of the problem.It is shown that while studying the evolution of a celestial body as a dissipative system in the framework of perturbed virial oscillations, the analytical expression for perturbing function can be derived, assuming the celestial body to be an oscillating electrical dipole emitting electromagnetic energy.The general covariant form of Jacobi's equation is derived and its spur is examined. It is shown that the scalar form of Jacobi's equation appears to be more universal than Newton's laws of motion from which it is derived.     

Exact solution of a key equation in a finite stellar atmosphere by the method of Laplace transform and linear singular operators

The key equation which commonly appears for radiative transfer in a finite stellar atmosphere having ground reflection according to Lambert's law is considered in this paper. The exact solution of this equation is obtained for surface quantities in terms of theX-Y equations of Chandrasekhar by the method of Laplace transform and linear singular operators. This exact method is widely applicable for obtaining the solution for surface quantities in a finite atmosphere.     

Soft X-Ray Absorption by High-Redshift Intergalactic Helium

Radio galaxies are uniquely useful as probes of large-scale structure since their uniform identification with giant elliptical galaxies out to high redshift means that the evolution of their bias factor can be predicted. As the initial stage in a project to study large-scale structure with radio galaxies, we have performed a small redshift survey, selecting 29 radio galaxies in the range of 0.19相似文献   

Constraining cold dark matter halo merger rates using the coagulation equations

We place additional constraints on the three parameters of the dark matter halo merger rate function recently proposed by Parkinson, Cole & Helly by utilizing Smoluchowski's coagulation equation, which must be obeyed by any binary merging process which conserves mass. We find that the constraints from Smoluchowski's equation are degenerate, limiting to a thin plane in the three-dimensional parameter space. This constraint is consistent with those obtained from fitting to N -body measures of progenitor mass functions, and provides a better match to the evolution of the overall dark matter halo mass function, particularly for the most massive haloes. We demonstrate that the proposed merger rate function does not permit an exact solution of Smoluchowski's equation and, therefore, the choice of parameters must reflect a compromise between fitting various parts of the mass function. The techniques described herein are applicable to more general merger rate functions, which may permit a more accurate solution of Smoluchowski's equation. The current merger rate solutions are most probably sufficiently accurate for the vast majority of applications.     

On the holographic dark energy in chameleon scalar-tensor cosmology

We study the holographic dark energy (HDE) model in generalized Brans-Dicke scenario with a non-minimal coupling between the scalar field and matter Lagrangian namely Chameleon Brans Dicke (CBD) mechanism. In this study we consider the interacting and non-interacting cases for two different cutoffs. The physical quantities of the model such as, equation of state (EoS) parameter, deceleration parameter and the evolution equation of dimensionless parameter of dark energy are obtained. We shall show that this model can describe the dynamical evolution of fraction parameter of dark energy in all epochs. Also we find the EoS parameter can cross the phantom divide line by suitable choices of parameters without any mines kinetic energy term.     

The evolution of adiabatic density perturbations

The evolution of adiabatic perturbations has been studied in Friedmann-Lemaître universes containing a mixture of non-interacting matter and radiation. The differential equation governing the growth of infinitely long adiabatic fluctuations has been derived. In the case of the Einstein-de Sitter universe, an exact solution is given to this differential equation.     

Energetics and thermal structure of the middle atmosphere

A study of a large number of temperature measurements in the middle atmosphere shows a much more complex thermal structure of this region than described in the U.S. Standard Atmosphere, 1976. The mesopause height which is generally assumed to be at 80 km varies between 70–100 km, often with two minima in temperature at about 70 and 100 km and a maximum between 80–85 km. By solving the energy balance equation and the equations of continuity, the physical significance of the observed thermal structure is discussed in terms of the energetics of the various regions of the middle atmosphere. It is shown that the solar u.v. radiation plays a major role only in the energy budget of the stratosphere and the lower thermosphere. The energetics of the mesosphere is primarily influenced by the dissipation of eddy energy. The temperature in this region is a good indicator of the eddy diffusivity and can be used in deriving the eddy diffusion coefficient.     

Status of physical observables of a friedmann universe in the classical and quantum hamiltonian formalisms

The paper is devoted to an investigation of the relationships between the classical Friedmann cosmology and the Dirac Hamiltonian approach to quantization of the universe, based on the simple but important example of a homogeneous universe filled with excitations of a scalar field. The method of gaugeless reduction is used to completely separate the sector of physical variables from the purely gauge sector, making it possible to find the relationship between cosmological observables in the Friedmann — Einstein sense and observables of the Dirac Hamiltonian formalism in the Narlikar conformai reference frame. Gaugeless reduction enabled us to establish that in the process of reduction, one of the variables of the nonphysical sector is converted into an invariant time parameter and cannot be treated as a dynamical variable in either the functional or the operator approach to quantization. It is shown that in this conversion of a variable into a time parameter, the Hartle-Hawking functional integral is the reason why the wave function of the Wheeler—De Witt (WDW) equation cannot be normalized and why an infinite gauge factor arises. The gaugeless reduction provides a certain recipe for mathematical and physical interpretation of the WDW equation and wave functions, the use of which makes their relationship to observational cosmology clear and transparent. It is shown, in particular, how the WDW wave function describes the Friedmann evolution with respect to proper time. Translated from Astrofizika, Vol. 40, No. 2, pp. 303–321, April–June, 1997.     

The vector approach to the problem of physical libration of the Moon. II. The nonlinear problem

By the new vector method in a nonlinear setting, a physical libration of the Moon is studied. Using the decomposition method on small parameters we derive the closed system of nine differential equations with terms of the first and second order of smallness. The conclusion is drawn that in the nonlinear case a connection between the librations in a longitude and latitude, though feeble, nevertheless exists; therefore, the physical libration already is impossible to subdivide into independent from each other forms of oscillations, as usually can be done. In the linear approach, ten characteristic frequencies and two special invariants of the problem are found. It is proved that, taking into account nonlinear terms, the invariants are periodic functions of time. Therefore, the stationary solution with zero frequency, formally supposing in the linear theory a resonance, in the nonlinear approach gains only small (proportional to e) periodic oscillations. Near to zero frequency of a resonance there is no, and solution of the nonlinear equations of physical libration is stable. The given nonlinear solution slightly modifies the previously unknown conical precession of the Moon’s spin axis. The character of nonlinear solutions near the basic forcing frequency Ω₁, where in the linear approach there are beats, is carefully studied. The average method on fast variables is obtained by the linear system of differential equations with almost periodic coefficients, which describe the evolution of these coefficients in a nonlinear problem. From this follows that the nonlinear components only slightly modify the specified beats; the interior period T ≈ 16.53 days appears 411 times less than the exterior one T ≈ 18.61 Julian years. In particular, with such a period the angle between ecliptic plane and Moon orbit plane also varies. Resonances, on which other researches earlier insisted, are not discovered. As a whole, the nonlinear analysis essentially improves and supplements a linear picture of the physical libration.     

Compressible Fluid Model for the Seismic Waves Generated by A Sunquake

The solar convection zone is modeled as a horizontally stratified atmosphere with a constant gravitational field and an adiabatic temperature gradient (a neutrally stratified polytrope). At equilibrium, the gas pressure and density decreases to zero at the solar surface so that the solar surface is treated as a free surface which is bounded by vacuum. The evolution of small amplitude perturbations about the equilibrium state is described by the linearized Euler equations for an inviscid compressible fluid with an adiabatic equation of state. A sunquake is initiated at time zero by means of an initial perturbation with a Gaussian velocity profile and the exact solution of the initial value problem is obtained in terms of a Fourier integral. Comparisons between theory and observations indicate that this highly simplified model is able to predict the propagation of sunquake waves across the solar surface with an error of approximately 10% or 20%.     

The solution for seasonal variations in the Earth's rate of rotation

New physical principles for an explanation of seasonal variations in the Earth's rate of rotation are proposed. It is thought that the variations are caused by a variation of the total energy of the Earth's atmosphere in the course of the planet's revolution about the Sun in elliptic orbit. Jacobi's virial equation for the Earth's atmosphere is derived from the Eulerian equations. The virial theorem is obtained. The existence of the relationship between Jacobi's function and potential energy of the atmosphere is confirmed. In the framework of this relationship, Jacobi's equation is reduced to the equation of unperturbed virial oscillations. The solution of the above-mentioned equation expresses the periodic virial oscillations of Jacobi's function (moment of inertia) of the Earth's atmosphere with time. The solution of the perturbed virial oscillation problem of the atmosphere-solid Earth system is obtained. The perturbation term in Jacobi's virial equation regards, in explicit form, the energy changes occurring in the atmosphere in the course of the planet's revolution about the Sun in elliptic orbit. The annual and semi-annual periodic variations in the Earth's rate of rotation can be considered as an astrometrical result following from the obtained solution. A satisfactory accord of the theoretical results with experimental data is shown.     

An efficient computational method for the approximate solution of nonlinear Lane-Emden type equations arising in astrophysics

The key purpose of this article is to introduce an efficient computational method for the approximate solution of the homogeneous as well as non-homogeneous nonlinear Lane-Emden type equations. Using proposed computational method given nonlinear equation is converted into a set of nonlinear algebraic equations whose solution gives the approximate solution to the Lane-Emden type equation. Various nonlinear cases of Lane-Emden type equations like standard Lane-Emden equation, the isothermal gas spheres equation and white-dwarf equation are discussed. Results are compared with some well-known numerical methods and it is observed that our results are more accurate.