首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 31 毫秒
1.
We report on an application of gas-kinetic BGK scheme to the computation of turbulent compressible convection in the stellar interior. After incorporating the Sub-grid Scale (SGS) turbulence model into the BGK scheme, we tested the effects of numerical parameters on the quantitative relationships among the thermodynamic variables, their fluctuations and correlations in a very deep, initially gravity-stratified stellar atmosphere. Comparison indicates that the thermal properties and dynamic properties are dominated by different aspects of numerical models separately. An adjustable Deardorff constant in the SGS model cμ = 0.25 and an amplitude of artificial viscosity in the gas-kinetic BGK scheme C2 = 0 are appropriate for the current study. We also calculated the density-weighted auto-and cross-correlation functions in Xiong's turbulent stellar convection theory based on which the gradient type of models of the non-local transport and the anisotropy of the turbulence were preliminarily studied. No universal relations or con-stant parameters were found for these models.  相似文献   

2.
We have analyzed the formation, structure, and dynamical evolution of the population of stars that escaped from open clusters by numerical simulations using S. Aarseth’s modified NBODY6 code. In the Galactic tidal field, the population of stars that escaped from a cluster is shown to be elongated along the orbit of the cluster symmetrically about its core in the form of stellar tails of increasing sizes. We analyze the parameters of stellar tails as a function of such initial simulation conditions as the number of stars, the cluster density, the eccentricity of the Galactic cluster orbit in the plane of the Galactic disk, and the z velocity component. As a result, we constructed a grid of model stellar tails of open clusters. The grid includes such time-dependent parameters of the stellar tails as the length, the cross section, the number of stars, the velocity distribution, etc. Our simulations allow us to clarify the origin of moving clusters and stellar streams and to assess the role of star clusters in forming the stellar velocity field in the solar neighborhood.  相似文献   

3.
We review current understanding of the interaction of magnetic fields with convective motions in stellar convection zones. Among the most exciting recent results is the discovery that magnetic fields need not primarily be confined to the stable layer below the convection zone; numerical simulations have shown that surprisingly, strong magnetic fields can be maintained in the interior of the convection zone.  相似文献   

4.
Flux expulsion is an important consequence of the interaction of magnetic fields with fluid convection and has been well studied for particular cases of steady, single-cell flows. Here we examine a related phenomenon in inhomogeneous turbulence using direct numerical simulations. To understand our numerical results, we analyse average properties of our model, and obtain mean transport coefficients which can be used to describe the approach of the system to its final state. For the kinematic problem these transport coefficients give an excellent prediction of the expulsion process; however, the enhanced transport is suppressed by dynamical back-reaction of the Lorentz force. Finally, we discuss the astrophysical implications for magnetic fields in stellar convection zones. Segregation of magnetic fields from turbulent motion not only allows strong toroidal fields to accumulate in regions of convective overshoot but also permits significant poloidal fields to be maintained by dynamo action in stars like the Sun.  相似文献   

5.
The evolution of a stellar, initially dipole type magnetosphere interacting with an accretion disk is investigated using numerical ideal MHD simulations. The simulations follow several 1000 Keplerian periods of the inner disk (for animated movies see http://www.aip.de～cfendt).Our model prescribes a Keplerian disk around a rotating star as a fixed boundary condition. The initial magnetic field distribution remains frozen into the star and the disk. The mass flow rate into the corona is fixed for both components. The initial dipole type magnetic field develops into a spherically radial outflow pattern with two main components – a disk wind and a stellar wind – both evolving into a quasi-stationary final state. A neutral field line divides both components, along which small plasmoids are ejected in irregular time intervals. The half opening angle of the stellar wind cone varies from 30° to55° depending on the ratio of the mass flow rates of disk wind and stellar wind. The maximum speed of the outflow is about the Keplerian speed at the inner disk radius. An axial jet forms during the first decades of rotations. However, this feature does not survive on the very long time scale and a pressure driven low velocity flow along the axis evolves. Within a cone of 15° along the axis the formation of knots may be observed if the stellar wind is weak. With the chosen mass flow rates and field strength we see almost no indication for a flow self-collimation. This is due to the weak net poloidal electric current in the magnetosphere which is in difference to typical jet models.  相似文献   

6.
Dynamical friction arises from the interaction of a perturber and the gravitational wake it excites in the ambient medium. We study the effects of the presence of a boundary on dynamical friction by studying analytically the interaction of perturber with uniform rectilinear motion in a uniform homogeneous medium with a reflecting planar boundary. Wake reflection at a medium’s boundary may occur at the edges of truncated disks perturbed by planetary or stellar companions as well as in numerical simulations of planet-disk interaction with no-outflow boundary conditions. In this paper, we show that the presence of the boundary modifies the behaviour of dynamical friction significantly. We find that perturbers are invariably pushed away from the boundary and reach a terminal subsonic velocity near Mach 0.37 regardless of initial velocity. Dynamical friction may even be reversed for Mach numbers less than 0.37 thereby accelerating instead of decelerating the perturber. Perturbers moving parallel to the boundary feel additional friction orthogonal to the direction of motion that is much stronger than the standard friction along the direction of motion. These results indicate that the common use of the standard Chandrasekhar formula as a short hand estimate of dynamical friction may be inadequate as observed in various numerical simulations.  相似文献   

7.
  相似文献   

8.
We have developed a new stellar evolution and oscillation code YNEV,which calculates the structures and evolutions of stars,taking into account hydrogen and helium burning.A nonlocal turbulent convection theory and an updated overshoot mixing model are optional in this code.The YNEV code can evolve low-and intermediate-mass stars from the pre-main sequence to a thermally pulsing asymptotic branch giant or white dwarf.The YNEV oscillation code calculates the eigenfrequencies and eigenfunctions of the adiabatic oscillations for a given stellar structure.The input physics and numerical scheme adopted in the code are introduced.Examples of solar models,stellar evolutionary tracks of low-and intermediate-mass stars with different convection theories(i.e.mixing-length theory and nonlocal turbulent convection theory),and stellar oscillations are shown.  相似文献   

9.
We investigate simulated turbulent flow within thermally driven stellar convection zones. Different driving sources are studied, including cooling at the top of the convectively unstable region, as occurs in surface convection zones; and heating at the base by nuclear burning. The transport of enthalpy and kinetic energy, and the distribution of turbulent kinetic energy dissipation are studied. We emphasize the importance of global constraints on shaping the quasi-steady flow characteristics, and present an analysis of turbulent convection which is posed as a boundary value problem that can be easily incorporated into standard stellar evolution codes for deep, efficient convection. Direct comparison is made between the theoretical analysis and the simulated flow and very good agreement is found. Some common assumptions traditionally used to treat quasi-steady turbulent flow in stellar models are briefly discussed. The importance and proper treatment of convective boundaries are indicated.  相似文献   

10.
We performed high-resolution simulations of two stellar collisions relevant for stars in globular clusters. We considered one head-on collision and one off-axis collision between two 0.6-M⊙ main-sequence stars. We show that a resolution of about 100 000 particles is sufficient for most studies of the structure and evolution of blue stragglers. We demonstrate conclusively that collision products between main-sequence stars in globular clusters do not have surface convection zones larger than 0.004 M⊙ after the collision, nor do they develop convection zones during the 'pre-main-sequence' thermal relaxation phase of their post-collision evolution. Therefore, any mechanism which requires a surface convection zone (i.e. chemical mixing or angular momentum loss via a magnetic wind) cannot operate in these stars. We show that no disc of material surrounding the collision product is produced in off-axis collisions. The lack of both a convection zone and a disc proves a continuing problem for the angular momentum evolution of blue stragglers in globular clusters.  相似文献   

11.
In our previous papers, we showed that at the final phases of the dynamical evolution of an open cluster, an extended population of stars elongated along its Galactic orbit, the stellar tail of the cluster, is formed. The tail stars that escaped from the cluster at different times move in a common orbit with low relative velocities. Experiencing a weak interaction with Galactic field stars, these objects, the relics of open clusters, can exist for a fairly long time. In this paper, we investigate the structures of such stellar tails in the nearest open clusters: Hyades, Pleiades, Praesepe, Alpha Persei, Coma, IC 2391, and IC 2602. To this end, we performed several numerical simulations of the dynamical evolution of these clusters in the tidal field of the Galaxy. Our computations of the dynamical evolution were based on known cluster age estimates and real Galactic orbits. The initial conditions were chosen in such a way that the parameters of the simulated clusters corresponded to their observed parameters. As a result, we obtained models of the stellar tails for the nearest open clusters and estimated such parameters of the tails as their sizes, densities, locations relative to the solar neighborhood, and others.  相似文献   

12.
We present a new method for constructing equilibrium phase models for stellar systems, which we call the iterative method. It relies on constrained, or guided evolution, so that the equilibrium solution has a number of desired parameters and/or constraints. This method is very powerful, to a large extent due to its simplicity. It can be used for mass distributions with an arbitrary geometry and a large variety of kinematical constraints. We present several examples illustrating it. Applications of this method include the creation of initial conditions for N -body simulations and the modelling of galaxies from their photometric and kinematic observations.  相似文献   

13.
E. Papini  L. Gizon  A. C. Birch 《Solar physics》2014,289(6):1919-1929
Linear time-domain simulations of acoustic oscillations are unstable in the stellar convection zone. To overcome this problem it is customary to compute the oscillations of a stabilized background stellar model. The stabilization affects the result, however. Here we propose to use a perturbative approach (running the simulation twice) to approximately recover the acoustic wave field while preserving seismic reciprocity. To test the method we considered a 1D standard solar model. We found that the mode frequencies of the (unstable) standard solar model are well approximated by the perturbative approach within 1 μHz for low-degree modes with frequencies near 3 mHz. We also show that the perturbative approach is appropriate for correcting rotational-frequency kernels. Finally, we comment that the method can be generalized to wave propagation in 3D magnetized stellar interiors because the magnetic fields have stabilizing effects on convection.  相似文献   

14.
We combine a convectively driven dynamo in a spherical shell with a nearly isothermal density-stratified cooling layer that mimics some aspects of a stellar corona to study the emergence and ejections of magnetic field structures. This approach is an extension of earlier models, where forced turbulence simulations were employed to generate magnetic fields. A spherical wedge is used which consists of a convection zone and an extended coronal region to ???1.5 times the radius of the sphere. The wedge contains a quarter of the azimuthal extent of the sphere and 150° in latitude. The magnetic field is self-consistently generated by the turbulent motions due to convection beneath the surface. Magnetic fields are found to emerge at the surface and are ejected to the coronal part of the domain. These ejections occur at irregular intervals and are weaker than in earlier work. We tentatively associate these events with coronal mass ejections on the Sun, even though our model of the solar atmosphere is rather simplistic.  相似文献   

15.
We report on a very large set of simulations of collisions between two main-sequence (MS) stars. These computations were carried out with the smoothed particle hydrodynamics method. Realistic stellar structure models for evolved MS stars were used. In order to sample an extended domain of initial parameters space (masses of the stars, relative velocity and impact parameter), more than 14 000 simulations were carried out. We considered stellar masses ranging between 0.1 and  75 M⊙  and relative velocities up to a few thousand km s−1. To limit the computational burden, a resolution of 1000–32 000 particles per star was used. The primary goal of this study was to build a complete data base from which the result of any collision can be interpolated. This allows us to incorporate the effects of stellar collisions with an unprecedented level of realism into dynamical simulations of galactic nuclei and other dense stellar clusters. We make the data describing the initial condition and outcome (mass and energy loss, angle of deflection) of all our simulations available on the Internet. We find that the outcome of collisions depends sensitively on the stellar structure and that, in most cases, using polytropic models is inappropriate. Published fitting formulae for the collision outcomes, established from a limited set of collisions, prove of limited use because they do not allow robust extrapolation to other stellar structures or relative velocities.  相似文献   

16.
Combination frequencies are observed in the Fourier spectra of pulsating DA and DB white dwarfs, along with frequencies that are associated with stellar gravity modes. They appear at the sum and difference frequencies of the stellar modes. Brickhill proposed that the combination frequencies result from mixing of the eigenmode signals by a depth-varying surface convection zone when undergoing pulsation. The depth changes cause time-dependent thermal impedance.
Following Brickhill's proposal, we developed analytical expressions for the amplitudes and phases of these combination frequencies. The parameters that appear in these expressions are the depth of the stellar convection zone when at rest, the sensitivity of this depth towards changes in the stellar effective temperature, the inclination angle of the stellar pulsation axis with respect to the line of sight, and lastly the spherical degrees of the eigenmodes involved in the mixing. Adopting credible values for these parameters, we apply our expressions to DA and DB variable white dwarfs. We find reasonable agreement between theory and observation, although some discrepancies remain unexplained. It is possible to identify the spherical degrees of the pulsation modes using the combination frequencies.  相似文献   

17.
We study the evolution of ionization fronts around the first protogalaxies by using high-resolution numerical cosmological (Λ+ cold dark matter, CDM, model) simulations and Monte Carlo radiative transfer methods. We present the numerical scheme in detail and show the results of test runs from which we conclude that the scheme is both fast and accurate. As an example of interesting cosmological application, we study the reionization produced by a stellar source of total mass M =2×108 M⊙ turning on at z ≈12, located at a node of the cosmic web. The study includes a spectral energy distribution of a zero-metallicity stellar population, and two initial mass functions (IMFs; Salpeter/Larson). The expansion of the ionization front (I-front) is followed as it breaks out from the galaxy and is channelled by the filaments into the voids, assuming (in a 2D representation) a characteristic butterfly shape. The ionization evolution is very well tracked by our scheme, as realized by the correct treatment of the channelling and shadowing effects resulting from overdensities. We confirm previous claims that both the shape of the IMF and the ionizing power metallicity dependence are important to correctly determine the reionization of the Universe.  相似文献   

18.
The bar formation is still an open problem in modern astrophysics. In this paper we present numerical simulations performed with the aim of analyzing the growth of the bar instability inside stellar-gaseous disks, where the star formation is triggered, and a central black hole is present. The aim of this paper is to point out the impact of such a central massive black hole on the growth of the bar. We use N-body-SPH simulations of the same isolated disk-to-halo mass systems harboring black holes with different initial masses and a different energy feedback on the surrounding gas. We compare the results of these simulations with the one of the same disk without a black hole in its center. We make the same comparison (disk with and without black hole) for a stellar disk in a fully cosmological scenario. A stellar bar, lasting 10 Gyrs, is present in all our simulations.  相似文献   

19.
We analyze the relationship between the mass of a spherical component and the minimum possible thickness of stable stellar disks. This relationship for real galaxies allows the lower limit on the dark halo mass to be estimated (the thinner the stable stellar disk is, the more massive the dark halo must be). In our analysis, we use both theoretical relations and numerical N-body simulations of the dynamical evolution of thin disks in the presence of spherical components with different density profiles and different masses. We conclude that the theoretical relationship between the thickness of disk galaxies and the mass of their spherical components is a lower envelope for the model data points. We recommend using this theoretical relationship to estimate the lower limit for the dark halo mass in galaxies. The estimate obtained turns out to be weak. Even for the thinnest galaxies, the dark halo mass within four exponential disk scale lengths must be more than one stellar disk mass.  相似文献   

20.
The results of numerical simulations of a gaseous disk in the potential of a stellar spiral density wave are presented. The conditions under which straightened spiral arm segments (rows) form in the gas component are studied. These features of the spiral structure were identified in a series of works by A.D. Chernin with coauthors. Gas-dynamic simulations have been performed for a wide range of model parameters: the pitch angle of the spiral pattern, the amplitude of the stellar spiral density wave, the disk rotation speed, and the temperature of the gas component. The results of 2D- and 3D-disk simulations are compared. The rows in the numerical simulations are shown to be an essentially nonstationary phenomenon. A statistical analysis of the distribution of geometric parameters for spiral patterns with rows in the observed galaxies and the constructed hydrodynamic models shows good agreement. In particular, the numerical simulations and observations of galaxies give 〈α〉 }~ 120° for the average angles between straight segments.  相似文献   

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

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