共查询到20条相似文献，搜索用时 31 毫秒
1.
We perform axisymmetric hydrodynamical simulations that describe the nonlinear outcome of the viscous overstability in dense planetary rings. These simulations are particularly relevant for Cassini observations of fine-scale structure in Saturn’s A and B-ring, which take the form of periodic microstructure on the 0.1 km scale, and irregular larger-scale variations on 1-10 km. Nonlinear wavetrains dominate all the simulations, and we associate them with the observed periodic microstructure. The waves can undergo small chaotic fluctuations in their phase and amplitude, and may be punctuated by more formidable ‘wave defects’ distributed on longer scales. It is unclear, however, whether the defects are connected to the irregular larger-scale variations observed by Cassini. The long-term behaviour of the simulations is dominated by the imposed boundary conditions, and more generally by the limitations of the local model we use: the shearing box. When periodic boundary conditions are imposed, the system eventually settles on a uniform travelling wave of a predictable wavelength, while reflecting boundaries, and boundaries with buffer zones, maintain a disordered saturated state. The simulations omit self-gravity, though we examine its influence in future work. 相似文献
2.
Three representative numerical simulations of the growth of the terrestrial planets by accretion of large protoplanets are presented. The mass and relative-velocity distributions of the bodies in these simulations are free to evolve simultaneously in response to close gravitational encounters and occasional collisions between bodies. The collisions between bodies, therefore, arise in a natural way and the assumption of expressions for the relative velocity distribution and the gravitational collision cross section is unnecessary. These simulations indicate that the growth of bodies with final masses approaching those of Venus and the Earth is possible, at least for the case of a two-dimensional system. Simulations assuming an initial uniform distribution of orbital eccentricities on the interval from 0 to emax are found to produce final states containing too many bodies with masses which are too small when emax < 0.10, while simulations with emax > 0.20 result in too many catastrophic collisions between bodies thus preventing rapid accretion of planetary-size bodies. The emax = 0.15 simulation ends with a state surprisingly similar to that of the present terrestrial planets and, therefore, provides a rough estimate of the range of radial sampling to be expected for the terrestrial planets. 相似文献
3.
We present N-body simulations of planetary accretion beginning with 1 km radius planetesimals in orbit about a 1 M⊙ star at 0.4 AU. The initial disk of planetesimals contains too many bodies for any current N-body code to integrate; therefore, we model a sample patch of the disk. Although this greatly reduces the number of bodies, we still track in excess of 105 particles. We consider three initial velocity distributions and monitor the growth of the planetesimals. The masses of some particles increase by more than a factor of 100. Additionally, the escape speed of the largest particle grows considerably faster than the velocity dispersion of the particles, suggesting impending runaway growth, although no particle grows large enough to detach itself from the power law size-frequency distribution. These results are in general agreement with previous statistical and analytical results. We compute rotation rates by assuming conservation of angular momentum around the center of mass at impact and that merged planetesimals relax to spherical shapes. At the end of our simulations, the majority of bodies that have undergone at least one merger are rotating faster than the breakup frequency. This implies that the assumption of completely inelastic collisions (perfect accretion), which is made in most simulations of planetary growth at sizes 1 km and above, is inappropriate. Our simulations reveal that, subsequent to the number of particles in the patch having been decreased by mergers to half its initial value, the presence of larger bodies in neighboring regions of the disk may limit the validity of simulations employing the patch approximation. 相似文献
4.
Three-dimensional numerical simulations of the atmospheric flow on giant planets using the primitive equations show that shallow thermal forcing confined to pressures near the cloud tops can produce deep zonal winds from the tropopause all the way down to the bottom of the atmosphere. These deep winds can attain speeds comparable to the zonal jet speeds within the shallow, forced layer; they are pumped by Coriolis acceleration acting on a deep meridional circulation driven by the shallow-layer eddies. In the forced layer, the flow reaches an approximate steady state where east-west eddy accelerations balance Coriolis accelerations acting on the meridional flow. Under Jupiter-like conditions, our simulations produce 25 to 30 zonal jets, similar to the number of jets observed on Jupiter and Saturn. The simulated jet widths correspond to the Rhines scale; this suggests that, despite the three-dimensional nature of the dynamics, the baroclinic eddies energize a quasi-two-dimensional inverse cascade modified by the β effect (where β is the gradient of the Coriolis parameter). In agreement with Jupiter, the jets can violate the barotropic and Charney-Stern stability criteria, achieving curvatures ∂2u/∂y2 of the zonal wind u with northward distance y up to 2β. The simulations exhibit a tendency toward neutral stability with respect to Arnol'd's second stability theorem in the upper troposphere, as has been suggested for Jupiter, although deviations from neutrality exist. When the temperature varies strongly with latitude near the equator, our simulations can also reproduce the stable equatorial superrotation with wind speeds greater than . Diagnostics show that barotropic eddies at low latitudes drive the equatorial superrotation. The simulations also broadly explain the distribution of jet-pumping eddies observed on Jupiter and Saturn. While idealized, these simulations therefore capture many aspects of the cloud-level flows on Jupiter and Saturn. 相似文献
5.
An iterative approach is used to construct spherically symmetric equilibrium models with an anisotropic velocity distribution. The potentialities of the method have been tested on models with known distribution functions, the Osipkov-Merritt models. It is shown that models that differ significantly from the Osipkov-Merritt models can be constructed. An N-body model of a dark halo with a density distribution that approximates the results of cosmological simulations (the Navarro-Frenk-White model) has been constructed. The anisotropy profile has been taken to be similar to that yielded by cosmological simulations. The constructed models can serve as direct input data for investigating the dynamics and stability of such systems in N-body simulations. 相似文献
6.
Numerical MHD simulations of 3D reconnection events in the solar corona have improved enormously over the last few years, not only in resolution, but also in their complexity, enabling more and more realistic modeling. Various ways to obtain the initial magnetic field, different forms of solar atmospheric models as well as diverse driving speeds and patterns have been employed. This study considers differences between simulations with stratified and non-stratified solar atmospheres, addresses the influence of the driving speed on the plasma flow and energetics, and provides quantitative formulas for mapping electric fields and dissipation levels obtained in numerical simulations to the corresponding solar quantities. The simulations start out from a potential magnetic field containing a null-point, obtained from a Solar and Heliospheric Observatory (SOHO) Michelson Doppler Imager (MDI) magnetogram magnetogram extrapolation approximately 8?hours before a C-class flare was observed. The magnetic field is stressed with a boundary motion pattern similar to?–?although simpler than?–?horizontal motions observed by SOHO during the period preceding the flare. The general behavior is nearly independent of the driving speed, and is also very similar in stratified and non-stratified models, provided only that the boundary motions are slow enough. The boundary motions cause a build-up of current sheets, mainly in the fan-plane of the magnetic null-point, but do not result in a flare-like energy release. The additional free energy required for the flare could have been partly present in non-potential form at the initial state, with subsequent additions from magnetic flux emergence or from components of the boundary motion that were not represented by the idealized driving pattern. 相似文献
7.
The diffusion mechanism of charged particles is one of the very im- portant problems in astrophysics. Recently, Matthaeus et al. (2003) developed a nonlinear guiding center theory (NLGC) of the perpendicular diffusion of en- ergetic charged particles with the parallel diffusion coefficient as an input. This theory agrees with the results of numerical simulations very well. In addition, Qin (2007) developed a nonlinear parallel diffusion theory (NLPA) following the idea of NLGC. Combining the NLGC with the NLPA, a new theory, NLGC-E was developed to solve the parallel and perpendicular diffusion coefficients simultaneously. In the NLGC theory, there is a coefficient a2 which is selected to be 1/3 for the best agreement with the results of numerical simulations. In this work we tried the different values of a2 for the NLGC-E theory to determine the optimal value which ?ts best the results of numerical simulations. 相似文献
8.
Richard H. Durisen Thomas W. Hartquist Megan K. Pickett 《Astrophysics and Space Science》2008,317(1-2):3-8
Numerical hydrodynamics simulations have established that disks which are evolved under the condition of local isothermality will fragment into small dense clumps due to gravitational instabilities when the Toomre stability parameter Q is sufficiently low. Because fragmentation through disk instability has been suggested as a gas giant planet formation mechanism, it is important to understand the physics underlying this process as thoroughly as possible. In this paper, we offer analytic arguments for why, at low Q, fragments are most likely to form first at the corotation radii of growing spiral modes, and we support these arguments with results from 3D hydrodynamics simulations. 相似文献
9.
This article is the third in a series working towards the construction of a realistic, evolving, non-linear force-free coronal-field model for the solar magnetic carpet. Here, we present preliminary results of 3D time-dependent simulations of the small-scale coronal field of the magnetic carpet. Four simulations are considered, each with the same evolving photospheric boundary condition: a 48-hour time series of synthetic magnetograms produced from the model of Meyer et al. (Solar Phys. 272, 29, 2011). Three simulations include a uniform, overlying coronal magnetic field of differing strength, the fourth simulation includes no overlying field. The build-up, storage, and dissipation of magnetic energy within the simulations is studied. In particular, we study their dependence upon the evolution of the photospheric magnetic field and the strength of the overlying coronal field. We also consider where energy is stored and dissipated within the coronal field. The free magnetic energy built up is found to be more than sufficient to power small-scale, transient phenomena such as nanoflares and X-ray bright points, with the bulk of the free energy found to be stored low down, between 0.5?–?0.8 Mm. The energy dissipated is currently found to be too small to account for the heating of the entire quiet-Sun corona. However, the form and location of energy-dissipation regions qualitatively agree with what is observed on small scales on the Sun. Future MHD modelling using the same synthetic magnetograms may lead to a higher energy release. 相似文献
10.
We present results from 44 simulations of late stage planetary accretion, focusing on the delivery of volatiles (primarily water) to the terrestrial planets. Our simulations include both planetary “embryos” (defined as Moon to Mars sized protoplanets) and planetesimals, assuming that the embryos formed via oligarchic growth. We investigate volatile delivery as a function of Jupiter's mass, position and eccentricity, the position of the snow line, and the density (in solids) of the solar nebula. In all simulations, we form 1-4 terrestrial planets inside 2 AU, which vary in mass and volatile content. In 44 simulations we have formed 43 planets between 0.8 and 1.5 AU, including 11 “habitable” planets between 0.9 and 1.1 AU. These planets range from dry worlds to “water worlds” with 100+oceans of water (1 ocean=1.5×1024 g), and vary in mass between 0.23M⊕ and 3.85M⊕. There is a good deal of stochastic noise in these simulations, but the most important parameter is the planetesimal mass we choose, which reflects the surface density in solids past the snow line. A high density in this region results in the formation of a smaller number of terrestrial planets with larger masses and higher water content, as compared with planets which form in systems with lower densities. We find that an eccentric Jupiter produces drier terrestrial planets with higher eccentricities than a circular one. In cases with Jupiter at 7 AU, we form what we call “super embryos,” 1-2M⊕ protoplanets which can serve as the accretion seeds for 2+M⊕ planets with large water contents. 相似文献
11.
Sergey L. Gratiy Andrew C. Walker David B. Goldstein Laurence M. Trafton 《Icarus》2010,207(1):394-408
Conflicting observations regarding the dominance of either sublimation or volcanism as the source of the atmosphere on Io and disparate reports on the extent of its spatial distribution and the absolute column abundance invite the development of detailed computational models capable of improving our understanding of Io’s unique atmospheric structure and origin. Improving upon previous models, Walker et al. (Walker, A.C., Gratiy, S.L., Levin, D.A., Goldstein, D.B., Varghese, P.L., Trafton, L.M., Moore, C.H., Stewart, B. [2009]. Icarus) developed a fully 3-D global rarefied gas dynamics model of Io’s atmosphere including both sublimation and volcanic sources of SO2 gas. The fidelity of the model is tested by simulating remote observations at selected wavelength bands and comparing them to the corresponding astronomical observations of Io’s atmosphere. The simulations are performed with a new 3-D spherical-shell radiative transfer code utilizing a backward Monte Carlo method. We present: (1) simulations of the mid-infrared disk-integrated spectra of Io’s sunlit hemisphere at 19 μm, obtained with TEXES during 2001-2004; (2) simulations of disk-resolved images at Lyman-α obtained with the Hubble Space Telescope (HST), Space Telescope Imaging Spectrograph (STIS) during 1997-2001; and (3) disk-integrated simulations of emission line profiles in the millimeter wavelength range obtained with the IRAM-30 m telescope in October-November 1999. We found that the atmospheric model generally reproduces the longitudinal variation in band depth from the mid-infrared data; however, the best match is obtained when our simulation results are shifted ∼30° toward lower orbital longitudes. The simulations of Lyman-α images do not reproduce the mid-to-high latitude bright patches seen in the observations, suggesting that the model atmosphere sustains columns that are too high at those latitudes. The simulations of emission line profiles in the millimeter spectral region support the hypothesis that the atmospheric dynamics favorably explains the observed line widths, which are too wide to be formed by thermal Doppler broadening alone. 相似文献
12.
The scattering properties of particulate rings with volume filling factors in the interval D=0.001-0.3 are studied, with photometric Monte Carlo ray tracing simulations combining the advantages of direct (photons followed from the source) and indirect methods (brightness as seen from the observing direction). Besides vertically homogeneous models, ranging from monolayers to classical many-particle thick rings, particle distributions obtained from dynamical simulations are studied, possessing a nonuniform vertical profile and a power law distribution of particle sizes. Self-gravity is not included to assure homogeneity in planar directions. Our main goal is to check whether the moderately flattened ring models predicted by dynamical simulations (with central plane D>0.1) are consistent with the basic photometric properties of Saturn's rings seen in ground-based observations, including the brightening near zero phase angle (opposition effect), and the brightening of the B-ring with increasing elevation angle (tilt effect). Our photometric simulations indicate that dense rings are typically brighter in reflected light than those with D→0, due to enhanced single scattering. For a vertically illuminated layer of identical particles this enhancement amounts at intermediate viewing elevations to roughly 1+2D. Increased single scattering is also obtained for low elevation illumination, further augmented at low phase angles α by the opposition brightening when D increases: the simulated opposition effect agrees very well with the Lumme and Bowell (1981, Astron. J. 86, 1694-1704) theoretical formula. For large α the total intensity may also decrease, due to reduced amount of multiple scattering. For the low (α=13°) and high (α=155°) phase angle geometries analyzed in Dones et al. (1993, Icarus 105, 184-215) the brightness change for D=0.1 amounts to 20% and −17%, respectively. In the case of an extended size distribution, dynamical simulations indicate that the smallest particles typically occupy a layer several times thicker than the largest particles. Even if the large particles form a dynamically dense system, a narrow opposition peak can arise due to mutual shadowing among the small particles: for example, a size distribution extending about two decades can account for the observed about 1° wide opposition peak, solely in terms of mutual shadowing. The reduced width of the opposition peak for extended size distribution is in accordance with Hapke's (1986, Icarus 67, 264-280) treatment for semi-infinite layers. Due to vertical profile and particle size distribution, the photometric behavior is sensitive to the viewing elevation: this can account for the tilt-effect of the B-ring, as dense and thus bright central parts of the ring become better visible for larger elevation, whereas in the case of smaller elevation, mainly low volume density upper layers are visible. Since multiple scattering is not involved, the explanation works also for albedo well below unity. Inclusion of nonzero volume density helps also to model some of the Voyager observations. For example, the discrepancy between predicted and observed brightness at large phase angles for much of the A-ring (Dones et al., 1993, Icarus 105, 184-215) is removed when the enhanced low α single scattering and reduced large α multiple scattering is allowed for. Also, a model with vertical thickness increasing with saturnocentric distance offers at least a qualitative explanation for the observed contrast reversal between the inner and outer A-ring in low and high phase Voyager images. Differences in local size distribution and thus on the effective D may also account for the contrast reversal in resonance sites. 相似文献
13.
Asteroid families are groupings of minor planets identified by clustering in their proper orbital elements; these objects have spectral signatures consistent with an origin in the break-up of a common parent body. From the current values of proper semimajor axes a of family members one might hope to estimate the ejection velocities with which the fragments left the putative break-up event (assuming that the pieces were ejected isotropically). However, the ejection velocities so inferred are consistently higher than N-body and hydro-code simulations, as well as laboratory experiments, suggest. To explain this discrepancy between today’s orbital distribution of asteroid family members and their supposed launch velocities, we study whether asteroid family members might have been ejected from the collision at low speeds and then slowly drifted to their current positions, via one or more dynamical processes. Studies show that the proper a of asteroid family members can be altered by two mechanisms: (i) close encounters with massive asteroids, and (ii) the Yarkovsky non-gravitational effect. Because the Yarkovsky effect for kilometer-sized bodies decreases with asteroid diameter D, it is unlikely to have appreciably moved large asteroids (say those with D > 15 km) over the typical family age (1-2 Gyr).For this reason, we numerically studied the mobility of family members produced by close encounters with main-belt, non-family asteroids that were thought massive enough to significantly change their orbits over long timescales. Our goal was to learn the degree to which perturbations might modify the proper a values of all family members, including those too large to be influenced by the Yarkovsky effect. Our initial simulations demonstrated immediately that very few asteroids were massive enough to significantly alter relative orbits among family members. Thus, to maximize gravitational perturbations in our 500-Myr integrations, we investigated the effect of close encounters on two families, Gefion and Adeona, that have high encounter probabilities with 1 Ceres, by far the largest asteroid in the main belt. Our results show that members of these families spreads in a of less than 5% since their formation. Thus gravitational interactions cannot account for the large inferred escape velocities.The effect of close encounters with massive asteroids is, however, not entirely negligible. For about 10% of the simulated bodies, close encounters increased the “inferred” ejection velocities from sub-100 m/s to values greater than 100 m/s, beyond what hydro-code and N-body simulations suggest are the maximum possible initial ejection velocity for members of Adeona and Gefion with D > 15 km. Thus this mechanism of mobility may be responsible for the unusually high inferred ejection speeds of a few of the largest members of these two families.To understand the orbital evolution of the entire family, including smaller members, we also performed simulations to account for the drift of smaller asteroids caused by the Yarkovsky effect. Our two sets of simulations suggest that the two families we investigated are relatively young compared to larger families like Koronis and Themis, which have estimated ages of about 2 Byr. The Adeona and Gefion families seems to be no more than 600 and 850 Myr old, respectively. 相似文献
14.
Stephen J. Kortenkamp 《Icarus》2005,175(2):409-418
Numerical simulations of the gravitational scattering of planetesimals by a protoplanet reveal that a significant fraction of scattered planetesimals can become trapped as so-called quasi-satellites in heliocentric 1:1 co-orbital resonance with the protoplanet. While trapped, these resonant planetesimals can have deep low-velocity encounters with the protoplanet that result in temporary or permanent capture onto highly eccentric prograde or retrograde circumplanetary orbits. The simulations include solar nebula gas drag and use planetesimals with diameters ranging from ∼1 to ∼1000 km. Initial protoplanet eccentricities range from ep=0 to 0.15 and protoplanet masses range from 300 Earth-masses (M⊕) down to 0.1M⊕. This mass range effectively covers the final masses of all planets currently thought to be in possession of captured satellites—Jupiter, Saturn, Neptune, Uranus, and Mars. For protoplanets on moderately eccentric orbits (ep?0.1) most simulations show from 5-20% of all scattered planetesimals becoming temporarily trapped in the quasi-satellite co-orbital resonance. Typically, 20-30% of the temporarily trapped quasi-satellites of all sizes came within half the Hill radius of the protoplanet while trapped in the resonance. The efficiency of the resonance trapping combined with the subsequent low-velocity circumplanetary capture suggests that this trapped-to-captured transition may be important not only for the origin of captured satellites but also for continued growth of protoplanets. 相似文献
15.
Gravitational accretion in the rings of Saturn is studied with local N-body simulations, taking into account the dissipative impacts and gravitational forces between particles. Common estimates of accretion assume that gravitational sticking takes place beyond a certain distance (Roche distance) where the self-gravity between a pair of ring particles exceeds the disrupting tidal force of the central object, the exact value of this distance depending on the ring particles' internal density. However, the actual physical situation in the rings is more complicated, the growth and stability of the particle groups being affected also by the elasticity and friction in particle impacts, both directly via sticking probabilities and indirectly via velocity dispersion, as well as by the shape, rotational state and the internal packing density of the forming particle groups. These factors are most conveniently taken into account via N-body simulations. In our standard simulation case of identical 1 m particles with internal density of solid ice, ρ=900 kg m−3, following the Bridges et al., 1984 elasticity law, we find accretion beyond a=137,000-146,000 km, the smaller value referring to a distance where transient aggregates are first obtained, and the larger value to the distance where stable aggregates eventually form in every experiment lasting 50 orbital periods. Practically the same result is obtained for a constant coefficient of restitution εn=0.5. In terms of rp parameter, the sum of particle radii normalized by their mutual Hill radius, the above limit for perfect accretion corresponds to rp<0.84. Increased dissipation (εn=0.1), or inclusion of friction (tangential force 10% of normal force) shifts the accretion region inward by about 5000 km. Accretion is also more efficient in the case of size distribution: with a q=3 power law extending over a mass range of 1000, accretion shifts inward by almost 10,000 km. The aggregates forming in simulations via gradual accumulation of particles are synchronously rotating. 相似文献
16.
Most stars reside in binary/multiple star systems; however, previous models of planet formation have studied growth of bodies orbiting an isolated single star. Disk material has been observed around both components of some young close binary star systems. Additionally, it has been shown that if planets form at the right places within such disks, they can remain dynamically stable for very long times. Herein, we numerically simulate the late stages of terrestrial planet growth in circumbinary disks around ‘close’ binary star systems with stellar separations 0.05 AU?aB?0.4 AU and binary eccentricities 0?eB?0.8. In each simulation, the sum of the masses of the two stars is 1 M⊙, and giant planets are included. The initial disk of planetary embryos is the same as that used for simulating the late stages of terrestrial planet formation within our Solar System by Chambers [Chambers, J.E., 2001. Icarus 152, 205-224], and around each individual component of the α Centauri AB binary star system by Quintana et al. [Quintana, E.V., Lissauer, J.J., Chambers, J.E., Duncan, M.J., 2002. Astrophys. J. 576, 982-996]. Multiple simulations are performed for each binary star system under study, and our results are statistically compared to a set of planet formation simulations in the Sun-Jupiter-Saturn system that begin with essentially the same initial disk of protoplanets. The planetary systems formed around binaries with apastron distances QB≡aB(1+eB)?0.2 AU are very similar to those around single stars, whereas those with larger maximum separations tend to be sparcer, with fewer planets, especially interior to 1 AU. We also provide formulae that can be used to scale results of planetary accretion simulations to various systems with different total stellar mass, disk sizes, and planetesimal masses and densities. 相似文献
17.
The saturation conditions for bending modes in inhomogeneous thin stellar disks that follow from an analysis of the dispersion relation are compared with those derived from N-body simulations. In the central regions of inhomogeneous disks, the reserve of disk strength against the growth of bending instability is smaller than that for a homogeneous layer. The spheroidal component (a dark halo, a bulge) is shown to have a stabilizing effect. The latter turns out to depend not only on the total mass of the spherical component, but also on the degree of mass concentration toward the center. We conclude that the presence of a compact (not necessarily massive) bulge in spiral galaxies may prove to be enough to suppress the bending perturbations that increase the disk thickness. This conclusion is corroborated by our N-body simulations in which we simulated the evolution of near-equilibrium, but unstable finite-thickness disks in the presence of spheroidal components. The final disk thickness at the same total mass of the spherical component (dark halo + bulge) was found to be much smaller than that in the simulations where a concentrated bulge is present. 相似文献
18.
C. K. Cramphorn 《Astronomy Letters》2001,27(3):135-139
In a recent study of dark mater N-body simulations, a scaling relation between the SZ decrement and the Thomson depth of a cluster of galaxies of the form ΔT r ∝τ T 2 has been found (Diaferio et al. 2000). In this paper, it will be shown that such a scaling relation arises if the intracluster gas is distributed similar to the dark matter density described by the NFW-profile and the finite spatial resolution of the numerical simulation is taken into account. It is furthermore investigated whether the ΔT r ∝τ T 2 relation holds for analytical models of an isothermal gas sphere in the gravitational potential of a dark matter halo distributed according to the NFW-profile, the available experimental data of SZE observations, and recent results from cosmological gas-dynamical simulations of clusters of galaxies. Combining such a relation with temperature estimates from X-ray observations would provide information about a dependence of T e on τT. The Thomson depth might therefore emerge as another important scaling parameter in studies of clusters of galaxies. 相似文献
19.
Ilian T. Iliev Ue-Li Pen Patrick McDonald Paul R. Shapiro Garrelt Mellema Marcelo A. Alvarez 《Astrophysics and Space Science》2009,320(1-3):39-43
Recently the numerical simulations of the process of reionization of the universe at z>6 have made a qualitative leap forward, reaching sufficient sizes and dynamic range to determine the characteristic scales of this process. This allowed making the first realistic predictions for a variety of observational signatures. We discuss recent results from large-scale radiative transfer and structure formation simulations on the observability of high-redshift Ly-α sources. We also briefly discuss the dependence of the characteristic scales and topology of the ionized and neutral patches on the reionization parameters. 相似文献
20.
We have carried out 1.25 pc resolution MHD simulations of the ISM, on a Cartesian grid of 0 ≤ (x, y) ≤ 1 kpc size in the galactic plane and ?10 ≤ z ≤ 10 kpc into the halo, thus being able to fully trace the time-dependent evolution of the galactic fountain. The simulations show that large scale gas streams emerge, driven by SN explosions, which are responsible for the formation and destruction of shocked compressed layers. The shocked gas can have densities as high as 800 cm?3 and lifetimes up to 15 Myr. The cold gas is distributed into filaments which tend to show a preferred orientation due to the anisotropy of the flow induced by the galactic magnetic field. Ram pressure dominates the flow in the unstable branch 102 < T ≤ 103.9 K, whereas for T ≤ 100 K (stable branch) magnetic pressure takes over. Near supernovae thermal and ram pressures determine the dynamics of the flow. Up to 80% of the mass in the disk is concentrated in the thermally unstable regime 102 < T ≤ 103.9 K with ～30% of the disk mass enclosed in the T ≤ 103 K gas. The hot gas in contrast is controlled by the thermal pressure, since magnetic field lines are swept towards the dense compressed walls. 相似文献