Source author record

E. Audit

E. Audit appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

astro-ph.GA astro-ph.SR astro-ph.HE astro-ph astro-ph.CO astro-ph.IM

Catalog footprint

What is connected

14works

6topics

4close collaborators

Actions

Connect this record

Open graph Browse works

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

preprint2020arXiv

A high-performance and portable asymptotic preserving radiation hydrodynamics code with the M1 model

Aims. We present a new radiation hydrodynamics code, called "ARK-RT" which uses a two-moment model with the M1 closure relation for radiative transfer. This code aims at being ready for high-performance computing, on exascale architectures. Methods. The two-moment model is solved using a finite volume scheme. The scheme is asymptotic preserving to capture accurately both optically thick and thin regimes. We also propose a well-balanced discretization of the radiative flux source term able to capture constant flux steady states with discontinuities in opacity. We use the library Trilinos for linear algebra and the package Kokkos allows us to reach high-performance computing and portability across different architectures, such as multi-core, many-core, and GP-GPU. Results. ARK-RT is able to reproduce standard tests in both free-streaming and diffusive limits, including purely radiative tests and radiation hydrodynamics ones. Using a time-implicit solver is profitable as soon as the time step given by the hydrodynamics is 50-100 times larger than the explicit time step for radiative transfer, depending on the preconditioner and the architecture. Albeit more work is needed to ensure stability in all circumstances. Using ARK-RT, we study the propagation of an ionization front in convective dense cores. We show that the ionization front is strongly stable against perturbations even with destabilizing convective motions. As a result, the presence of instabilities should be interpreted with caution. Overall, ARK-RT is well-suited to study many astrophysical problems involving convection and radiative transfer such as the dynamics of H ii regions in massive pre-stellar dense cores and future applications could include planetary atmospheres.

preprint2014arXiv

Ionization compression impact on dense gas distribution and star formation, Probability density functions around H ii regions as seen by Herschel

Ionization feedback should impact the probability distribution function (PDF) of the column density around the ionized gas. We aim to quantify this effect and discuss its potential link to the Core and Initial Mass Function (CMF/IMF). We used in a systematic way Herschel column density maps of several regions observed within the HOBYS key program: M16, the Rosette and Vela C molecular cloud, and the RCW 120 H ii region. We fitted the column density PDFs of all clouds with two lognormal distributions, since they present a double-peak or enlarged shape in the PDF. Our interpretation is that the lowest part of the column density distribution describes the turbulent molecular gas while the second peak corresponds to a compression zone induced by the expansion of the ionized gas into the turbulent molecular cloud. The condensations at the edge of the ionized gas have a steep compressed radial profile, sometimes recognizable in the flattening of the power-law tail. This could lead to an unambiguous criterion able to disentangle triggered from pre-existing star formation. In the context of the gravo-turbulent scenario for the origin of the CMF/IMF, the double peaked/enlarged shape of the PDF may impact the formation of objects at both the low-mass and the high-mass end of the CMF/IMF. In particular a broader PDF is required by the gravo-turbulent scenario to fit properly the IMF with a reasonable initial Mach number for the molecular cloud. Since other physical processes (e.g. the equation of state and the variations among the core properties) have already been suggested to broaden the PDF, the relative importance of the different effects remains an open question.

preprint2013arXiv

Pillars and globules at the edges of H ii regions, Confronting Herschel observations and numerical simulations

Pillars and globules are present in many high-mass star-forming regions, such as the Eagle nebula (M16) and the Rosette molecular cloud, and understanding their origin will help characterize triggered star formation. The formation mechanisms of these structures are still being debated. Recent numerical simulations have shown how pillars can arise from the collapse of the shell in on itself and how globules can be formed from the interplay of the turbulent molecular cloud and the ionization from massive stars. The goal here is to test this scenario through recent observations of two massive star-forming regions, M16 and Rosette. The column density structure of the interface between molecular clouds and H ii regions was characterized using column density maps obtained from far-infrared imaging of the Herschel HOBYS key programme. Then, the DisPerSe algorithm was used on these maps to detect the compressed layers around the ionized gas and pillars in different evolutionary states. Finally, their velocity structure was investigated using CO data, and all observational signatures were tested against some distinct diagnostics established from simulations. The column density profiles have revealed the importance of compression at the edge of the ionized gas. The velocity properties of the structures, i.e. pillars and globules, are very close to what we predict from the numerical simulations. We have identified a good candidate of a nascent pillar in the Rosette molecular cloud that presents the velocity pattern of the shell collapsing on itself, induced by a high local curvature. Globules have a bulk velocity dispersion that indicates the importance of the initial turbulence in their formation, as proposed from numerical simulations. Altogether, this study re-enforces the picture of pillar formation by shell collapse and globule formation by the ionization of highly turbulent clouds.

preprint2013arXiv

Simulations of protostellar collapse using multigroup radiation hydrodynamics. II. The second collapse

Star formation begins with the gravitational collapse of a dense core inside a molecular cloud. As the collapse progresses, the centre of the core begins to heat up as it becomes optically thick. The temperature and density in the centre eventually reach high enough values where fusion reactions can ignite; the protostar is born. This sequence of events entail many physical processes, of which radiative transfer is of paramount importance. Many simulations of protostellar collapse make use of a grey treatment of radiative transfer coupled to the hydrodynamics. However, interstellar gas and dust opacities present large variations as a function of frequency. In this paper, we follow-up on a previous paper on the collapse and formation of Larson's first core using multigroup radiation hydrodynamics (Paper I) by extending the calculations to the second phase of the collapse and the formation of Larson's second core. We have made the use of a non-ideal gas equation of state as well as an extensive set of spectral opacities in a spherically symmetric fully implicit Godunov code to model all the phases of the collapse of a 0.1, 1 and 10 solar mass cloud cores. We find that, for a same central density, there are only small differences between the grey and multigroup simulations. The first core accretion shock remains supercritical while the shock at the second core border is found to be strongly subcritical with all the accreted energy being transfered to the core. The size of the first core was found to vary somewhat in the different simulations (more unstable clouds form smaller first cores) while the size, mass and temperature of the second cores are independent of initial cloud mass, size and temperature. Our simulations support the idea of a standard (universal) initial second core size of 0.003 AU and mass 0.0014 solar masses.

preprint2013arXiv

The influence of frequency-dependent radiative transfer on the structures of radiative shocks

Radiative shocks are shocks in a gas where the radiative energy and flux coming from the very hot post-shock material are non-negligible in the shock's total energy budget, and are often large enough to heat the material ahead of the shock. Many simulations of radiative shocks, both in the contexts of astrophysics and laboratory experiments, use a grey treatment of radiative transfer coupled to the hydrodynamics. However, the opacities of the gas show large variations as a function of frequency and this needs to be taken into account if one wishes to reproduce the relevant physics. We have performed radiation hydrodynamics simulations of radiative shocks in Ar using multigroup (frequency dependent) radiative transfer with the HERACLES code. The opacities were taken from the ODALISC database. We show the influence of the number of frequency groups used on the dynamics and morphologies of subcritical and supercritical radiative shocks in Ar gas, and in particular on the extent of the radiative precursor. We find that simulations with even a low number of groups show significant differences compared to single-group (grey) simulations, and that in order to correctly model such shocks, a minimum number of groups is required. Results appear to eventually converge as the number of groups increases above 50. We were also able to resolve in our simulations of supercritical shocks the adaptation zones which connect the cooling layer to the final post-shock state and the precursor. Inside these adaptation zones, we find that the radiative flux just ahead of the shock in one or several high-opacity groups can heat the gas to a temperature higher than the post-shock temperature. Through the use of Hugoniot curves, we have checked the consistency of our radiation hydrodynamics scheme by showing that conservation of mass, momentum and energy (including radiative flux) holds. ABRIDGED.

preprint2013arXiv

The structure of the thermally bistable and turbulent atomic gas in the local interstellar medium

This paper is a numerical study of the condensation of the warm neutral medium (WNM) into cold neutral medium (CNM) structures under the effect of turbulence and thermal instability. Using low resolution simulations we explored the impact of the WNM initial density and properties of the turbulence (stirring in Fourier with a varying mix of solenoidal and compressive modes) on the cold gas formation. Two sets of initial conditions which match the observations were selected to produce high resolution simulations (1024^3) allowing to study in details the properties of the produced dense structures. For typical values of the density, pressure and velocity dispersion of the WNM in the solar neighborhood, the turbulent motions of the HI can not provoque the phase transition from WNM to CNM, whatever their amplitude and their distribution in solenoidal and compressive modes. On the other hand we show that a quasi-isothermal increase in WNM density of a factor of 2 to 4 is enough to induce the phase transition, leading to the transition of about 40 percent of the gas to the cold phase within 1 Myr. Given the observed properties of the HI in the local ISM, the WNM and individual CNM structures in the local ISM are sub or transsonic and their dynamics are tightly interwoven. The velocity field bears the evidence of subsonic turbulence with a 2D power spectrum following the Kolmogorov law as P(k) \sim k^{-8/3} while the density is highly contrasted with a singificantly shallower power spectrum, reminiscent of what is observed in the cold ISM. Supra-thermal line width observed for CNM might be the result of relative velocity between cold structures. Finally, the cold structures denser than 5 cm^{-3} reproduce well the laws M \sim L^{2.25-2.28} and sigma(v) \sim 0.5-0.8L^{1/3} generally observed in molecular clouds.

preprint2012arXiv

3D simulations of globules and pillars formation around HII regions: turbulence and shock curvature

We investigate the interplay between the ionization radiation from massive stars and the turbulence inside the surrounding molecular gas thanks to 3D numerical simulations. We used the 3D hydrodynamical code HERACLES to model an initial turbulent medium that is ionized and heated by an ionizing source. Three different simulations are performed with different mean Mach numbers (1, 2 and 4). A non-equilibrium model for the ionization and the associated thermal processes was used. This revealed to be crucial when turbulent ram pressure is of the same order as the ionized-gas pressure. The density structures initiated by the turbulence cause local curvatures of the dense shell formed by the ionization compression. When the curvature of the shell is sufficient, the shell collapse on itself to form a pillar while a smaller curvature leads to the formation of dense clumps that are accelerated with the shell and therefore remain in the shell during the simulation. When the turbulent ram pressure of the cold gas is sufficient to balance the ionized-gas pressure, some dense-gas bubbles have enough kinetic energy to penetrate inside the ionized medium, forming cometary globules. This suggests a direct relation in the observations between the presence of globules and the relative importance of the turbulence compared to the ionized-gas pressure. The probability density functions present a double peak structure when the turbulence is low relative to the ionized-gas pressure. This could be used in observations as an indication of the turbulence inside molecular clouds.

preprint2012arXiv

Globules and pillars seen in the [CII] 158 micron line with SOFIA

Molecular globules and pillars are spectacular features, found only in the interface region between a molecular cloud and an HII-region. Impacting Far-ultraviolet (FUV) radiation creates photon dominated regions (PDRs) on their surfaces that can be traced by typical cooling lines. With the GREAT receiver onboard SOFIA we mapped and spectrally resolved the [CII] 158 micron atomic fine-structure line and the highly excited 12CO J=11-10 molecular line from three objects in Cygnus X (a pillar, a globule, and a strong IRAS source). We focus here on the globule and compare our data with existing Spitzer data and recent Herschel Open-Time PACS data. Extended [CII] emission and more compact CO-emission was found in the globule. We ascribe this emission mainly to an internal PDR, created by a possibly embedded star-cluster with at least one early B-star. However, external PDR emission caused by the excitation by the Cyg OB2 association cannot be fully excluded. The velocity-resolved [CII] emission traces the emission of PDR surfaces, possible rotation of the globule, and high-velocity outflowing gas. The globule shows a velocity shift of ~2 km/s with respect to the expanding HII-region, which can be understood as the residual turbulence of the molecular cloud from which the globule arose. This scenario is compatible with recent numerical simulations that emphazise the effect of turbulence. It is remarkable that an isolated globule shows these strong dynamical features traced by the [CII]-line, but it demands more observational studies to verify if there is indeed an embedded cluster of B-stars.

preprint2011arXiv

3D simulations of pillars formation around HII regions: the importance of shock curvature

Radiative feedback from massive stars is a key process to understand how HII regions may enhance or inhibit star formation in pillars and globules at the interface with molecular clouds. We aim to contribute to model the interactions between ionization and gas clouds to better understand the processes at work. We study in detail the impact of modulations on the cloud-HII region interface and density modulations inside the cloud. We run three-dimensional hydrodynamical simulations based on Euler equations coupled with gravity using the HERACLES code. We implement a method to solve ionization/recombination equations and we take into account typical heating and cooling processes at work in the interstellar medium and due to ionization/recombination physics. UV radiation creates a dense shell compressed between an ionization front and a shock ahead. Interface modulations produce a curved shock that collapses on itself leading to stable growing pillar-like structures. The narrower the initial interface modulation, the longer the resulting pillar. We interpret pillars resulting from density modulations in terms of the ability of these density modula- tions to curve the shock ahead the ionization front. The shock curvature is a key process to understand the formation of structures at the edge of HII regions. Interface and density modulations at the edge of the cloud have a direct impact on the morphology of the dense shell during its formation. Deeper in the cloud, structures have less influence due to the high densities reached by the shell during its expansion.

preprint2011arXiv

A numerical model for multigroup radiation hydrodynamics

We present in this paper a multigroup model for radiation hydrodynamics to account for variations of the gas opacity as a function of frequency. The entropy closure model (M1) is applied to multigroup radiation transfer in a radiation hydrodynamics code. In difference from the previous grey model, we are able to reproduce the crucial effects of frequency-variable gas opacities, a situation omnipresent in physics and astrophysics. We also account for the energy exchange between neighbouring groups which is important in flows with strong velocity divergence. These terms were computed using a finite volume method in the frequency domain. The radiative transfer aspect of the method was first tested separately for global consistency (reversion to grey model) and against a well established kinetic model through Marshak wave tests with frequency dependent opacities. Very good agreement between the multigroup M1 and kinetic models was observed in all tests. The successful coupling of the multigroup radiative transfer to the hydrodynamics was then confirmed through a second series of tests. Finally, the model was linked to a database of opacities for a Xe gas in order to simulate realistic multigroup radiative shocks in Xe. The differences with the previous grey models are discussed.

preprint2009arXiv

On the Structure of the Turbulent Interstellar Clouds

{It is well established that the atomic interstellar hydrogen is filling the galaxies and constitutes the building blocks of molecular clouds.} {To understand the formation and the evolution of molecular clouds, it is necessary to investigate the dynamics of turbulent and thermally bistable as well as barotropic flows.} {We perform high resolution 3-dimensional hydrodynamical simulations of 2-phase, isothermal and polytropic flows.} {We compare the density PDF and Mach number density relation in the various simulations and conclude that 2-phase flows behave rather differently than polytropic flows. We also extract the clumps and study their statistical properties such as mass spectrum, mass-size relation and internal velocity dispersion. In each case, it is found that the behaviour is well represented by a simple powerlaw. While the various exponents inferred are very similar for the 2-phase, isothermal and polytropic flows, we nevertheless find significant differences, such as for example the internal velocity dispersion which is smaller for 2-phase flows.} {The structures statistics are very similar to what has been inferred from observations, in particular the mass spectrum, the mass-size relation and the velocity dispersion-size relation are all powerlaws whose indices well agree with the observed values. Our results suggest that in spite of various statistics being similar for 2-phase and polytropic flows, they nevertheless present significant differences, stressing the necessity to consider the proper thermal structure of the interstellar atomic hydrogen for computing its dynamics as well as the formation of molecular clouds.}

preprint2009arXiv

Protostellar collapse: radiative and magnetic feedbacks on small scale fragmentation

It is established that both radiative transfer and magnetic field have a strong impact on the collapse and the fragmentation of prestellar dense cores, but no consistent calculation exists yet at such scales. We present original AMR calculations including magnetic field (in the ideal MHD limit) and radiative transfer, within the Flux Limited Diffusion approximation, of the collapse of a 1 solar mass dense core. We compare the results with calculations performed with a barotropic EOS. We show that radiative transfer has an important impact on the collapse and the fragmentation, through the cooling or heating of the gas, and is complementary of the magnetic field. A larger field yields a stronger magnetic braking, increasing the accretion rate and thus the effect of the radiative feedback. Even for a strongly magnetized core, where the dynamics of the collapse is dominated by the magnetic field, radiative transfer is crucial to determine the temperature and optical depth distributions, two potentially accessible observational diagnostics. A barotropic EOS cannot account for realistic fragmentation. The diffusivity of the numerical scheme, however, is found to strongly affect the output of the collapse, leading eventually to spurious fragmentation. Both radiative transfer and magnetic field must be included in numerical calculations of star formation to obtain realistic collapse configurations and observable signatures. Nevertheless, the numerical resolution and the robustness of the solver are of prime importance to obtain reliable results. When using an accurate solver, the fragmentation is found to always remain inhibited by the magnetic field, at least in the ideal MHD limit, even when radiative transfer is included.

preprint2009arXiv

The dusty, albeit ultraviolet bright infancy of galaxies

The largest galaxies acquire their mass early on, when the Universe is still youthful. Cold streams violently feed these young galaxies a vast amount of fresh gas, resulting in very efficient star formation. Using a well resolved hydrodynamical simulation of galaxy formation, we demonstrate that these mammoth galaxies are already in place a couple of billion years after the Big Bang. Contrary to local starforming galaxies, where dust re-emits a large part of the stellar ultraviolet (UV) light at infrared and sub-millimetre wavelengths, our self-consistent modelling of dust extinction predicts that a substantial fraction of UV photons should escape from primordial galaxies. Such a model allows us to compute reliably the number of high redshift objects as a function of luminosity, and yields galaxies whose UV luminosities closely match those measured in the deepest observational surveys available. This agreement is remarkably good considering our admittedly still simple modelling of the interstellar medium (ISM) physics. The luminosity functions (LF) of virtual UV luminous galaxies coincide with the existing data over the whole redshift range from 4 to 7, provided cosmological parameters are set to their currently favoured values. Despite their considerable emission at short wavelengths, we anticipate that the counterparts of the brightest UV galaxies will be detected by future sub-millimetre facilities like ALMA

preprint2006arXiv

On the structure of the turbulent interstellar atomic hydrogen. I- Physical characteristics

{We study in some details the statistical properties of the turbulent 2-phase interstellar atomic gas.{We present high resolution bidimensional numerical simulations of the interstellar atomic hydrogen which describe it over 3 to 4 orders of magnitude in spatial scales.}{The simulations produce naturally small scale structures having either large or small column density. It is tempting to propose that the former are connected to the tiny small scale structures observed in the ISM. We compute the mass spectrum of CNM structures and find that ${\cal N}(M) dM \propto M ^{-1.7} dM$, which is remarkably similar to the mass spectrum inferred for the CO clumps. We propose a theoretical explanation based on a formalism inspired from the Press & Schecter (1974) approach and used the fact that the turbulence within WNM is subsonic. This theory predicts ${\cal N}(M) \propto M ^{-5/3}$ in 2D and ${\cal N}(M) \propto M ^{-16/9}$ in 3D. We compute the velocity and the density power-spectra and conclude that, although the latter is rather flat, as observed in supersonic isothermal simulations, the former follows the Kolmogorov prediction and is dominated by its solenoidal component. This is due to the bistable nature of the flow which produces large density fluctuations even when the rms Mach number (of WNM) is not large. We also find that, whereas the energy at large scales is mainly in the WNM, at smaller scales, it is dominated by the kinetic energy of the CNM fragments.}

E. Audit

What is connected

Connect this record

See the researcher in context

Building this map preview

14 published item(s)

A high-performance and portable asymptotic preserving radiation hydrodynamics code with the M1 model

Ionization compression impact on dense gas distribution and star formation, Probability density functions around H ii regions as seen by Herschel

Pillars and globules at the edges of H ii regions, Confronting Herschel observations and numerical simulations

Simulations of protostellar collapse using multigroup radiation hydrodynamics. II. The second collapse

The influence of frequency-dependent radiative transfer on the structures of radiative shocks

The structure of the thermally bistable and turbulent atomic gas in the local interstellar medium

3D simulations of globules and pillars formation around HII regions: turbulence and shock curvature

Globules and pillars seen in the [CII] 158 micron line with SOFIA

3D simulations of pillars formation around HII regions: the importance of shock curvature

A numerical model for multigroup radiation hydrodynamics

On the Structure of the Turbulent Interstellar Clouds

Protostellar collapse: radiative and magnetic feedbacks on small scale fragmentation

The dusty, albeit ultraviolet bright infancy of galaxies

On the structure of the turbulent interstellar atomic hydrogen. I- Physical characteristics