Source author record

S. A. E. G. Falle

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

20works
7topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

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.

Published work

20 published item(s)

preprint2022arXiv

Shocking interactions of supernova remnants with atomic and molecular clouds -- the interplay between shocks, thermal instability and gravity in the large cloud regime

Using the adaptive mesh refinement code MG, we perform 3D hydrodynamic simulations of a supernova-cloud interaction in the "large cloud regime". The cloud is initially atomic and evolving due to the thermal instability (TI) and gravity. We study interactions in a "pre-TI" and "post-TI" stage when cold and dense clumps are present, and compare these results to idealised shock-cloud scenarios in the "small cloud regime", and a scenario without shocks. On aggregate, the supernova disruption is significantly weaker than that from an idealised shock due to the supernova impact being instantaneous, and not continuous. In both supernova-cloud interactions, we observe two shocks impact the cloud, followed by the development of a weak 10 km s$^{-1}$ upstream flow on the cloud interface, and a global ambient pressure drop. When the cloud is still atomic, it expands due to this drop. Additionally, the TI is triggered at the front of the cloud, causing the formation of a cap-like structure with clumps embedded inside. The upstream flow converges in this region, resulting in a lobe-like cloud morphology. When the cloud is molecular, the transmitted shock disrupts the inter-clump material and causes the clumps' outer envelopes to expand slightly and form tail-like morphologies. These effects are less pronounced than those in our shock-cloud scenarios, and more pronounced that those in our un-shocked scenario. After 3.5 Myrs, the effects from the supernova decay and the cloud returns to an almost indistinguishable state from an un-shocked cloud, in spite of the global ambient pressure drop. In neither supernova-cloud scenario do we see any local gravitational collapse.

preprint2020arXiv

Interactions of a shock with a molecular cloud at various stages of its evolution due to thermal instability and gravity

Using the adaptive mesh refinement code MG, we perform hydrodynamic simulations of the interaction of a shock with a molecular cloud evolving due to thermal instability and gravity. To explore the relative importance of these processes, three case studies are presented. The first follows the formation of a molecular cloud out of an initially quiescent atomic medium due to the effects of thermal instability and gravity. The second case introduces a shock whilst the cloud is still in the warm atomic phase, and the third scenario introduces a shock once the molecular cloud has formed. The shocks accelerate the global collapse of the clouds with both experiencing local gravitational collapse prior to this. When the cloud is still atomic, the evolution is shock dominated and structures form due to dynamical instabilities within a radiatively cooled shell. While the transmitted shock can potentially trigger the thermal instability, this is prevented as material is shocked multiple times on the order of a cloud crushing time-scale. When the cloud is molecular, the post-shock flow is directed via the pre-existing structure through low-density regions in the inter-clump medium. The clumps are accelerated and deformed as the flow induces clump-clump collisions and mergers that collapse under gravity. For a limited period, both shocked cases show a mixture of Kolmogorov and Burgers turbulence-like velocity and logarithmic density power spectra, and strongly varying density spectra. The clouds presented in this work provide realistic conditions that will be used in future feedback studies.

preprint2016arXiv

Ambipolar diffusion regulated collapse of filaments threaded by perpendicular magnetic fields

We numerically reproduce the density profiles for filaments that are in magnetohydrostatic and pressure equilibrium with their surroundings obtained in Tomisaka (2014) and show that these equilibria are dynamically stable. If the effect of ambipolar diffusion is considered, these filaments lose magnetic support initiating cloud collapse. The filaments do not lose magnetic flux. Rather the magnetic flux is redistributed within the filament from the centre towards the envelope. The rate of the collapse is inversely proportional to the fractional ionisation and two gravitationally-driven ambipolar diffusion regimes for the collapse are observed as predicted in Mouschovias & Morton (1991). For high values of the ionisation coefficient, that is $X \geq 10^{-7}$, the gas is strongly coupled to the magnetic field and the Jeans length is larger than the ambipolar diffusion length scale. Then the collapse is governed by magnetically-regulated ambipolar diffusion. For $X \lesssim 10^{-8}$, the gas is weakly coupled to the magnetic field and the magnetic support is removed by gravitationally-dominated ambipolar diffusion. Here, neutrals and ions only collide sporadically, that is the ambipolar diffusion length scale is larger than the Jeans length. When decaying turbulence is included, additional support is provided to the filament. This slows down the collapse of the filament even in the absence of a magnetic field. When a magnetic field is present, the collapse rate increases by a ratio smaller than for the non-magnetic case. This is because of a speed-up of the ambipolar diffusion due to larger magnetic field gradients generated by the turbulence and because the ambipolar diffusion aids the dissipation of turbulence below the ambipolar diffusion length scale. The highest increase in the rate is observed for the lowest ionisation coefficient and the highest turbulent intensity.

preprint2016arXiv

MHD simulation of the formation of clumps and filaments in quiescent diffuse medium by thermal instability

We have used the AMR hydrodynamic code, MG, to perform idealised 3D MHD simulations of the formation of clumpy and filamentary structure in a thermally unstable medium without turbulence. A stationary thermally unstable spherical diffuse atomic cloud with uniform density in pressure equilibrium with low density surroundings was seeded with random density variations and allowed to evolve. A range of magnetic field strengths threading the cloud have been explored, from beta=0.1 to beta=1.0 to the zero magnetic field case (beta=infinity), where beta is the ratio of thermal pressure to magnetic pressure. Once the density inhomogeneities had developed to the point where gravity started to become important, self-gravity was introduced to the simulation. With no magnetic field, clouds and clumps form within the cloud with aspect ratios of around unity, whereas in the presence of a relatively strong field (beta=0.1) these become filaments, then evolve into interconnected corrugated sheets that are predominantly perpendicular to the magnetic field. With magnetic and thermal pressure equality (beta=1.0), filaments, clouds and clumps are formed. At any particular instant, the projection of the 3D structure onto a plane parallel to the magnetic field, i.e. a line of sight perpendicular to the magnetic field, resembles the appearance of filamentary molecular clouds. The filament densities, widths, velocity dispersions and temperatures resemble those observed in molecular clouds. In contrast, in the strong field case beta=0.1, projection of the 3D structure along a line of sight parallel to the magnetic field reveals a remarkably uniform structure.

preprint2016arXiv

The interaction of hydrodynamic shocks with self-gravitating clouds

We describe the results of 3D simulations of the interaction of hydrodynamic shocks with Bonnor-Ebert spheres performed with an Adaptive Mesh Refinement code. The calculations are isothermal and the clouds are embedded in a medium in which the sound speed is either four or ten times that in the cloud. The strengths of the shocks are such that they induce gravitational collapse in some cases and not in others and we derive a simple estimate for the shock strength required for this to occur. These results are relevant to dense cores and Bok globules in star forming regions subjected to shocks produced by stellar feedback.

preprint2014arXiv

Numerical Simulations of a Shock Interacting with Multiple Magnetized Clouds

We present 2D adiabatic magnetohydrodynamic (MHD) simulations of a shock interacting with groups of two or three cylindrical clouds. We study how the presence of a nearby cloud influences the dynamics of this interaction, and explore the resulting differences and similarities in the evolution of each cloud. The understanding gained from this small-scale study will help to interpret the behaviour of systems with many 10's or 100's of clouds. We observe a wide variety of behaviour in the interactions studied, which is dependent on the initial positions of the clouds and the orientation and strength of the magnetic field. We find: i) some clouds are stretched along their field-lines, whereas others are confined by their field-lines; ii) upstream clouds may accelerate past downstream clouds (though magnetic tension can prevent this); iii) clouds may also change their relative positions transverse to the direction of shock propagation as they "slingshot" past each other; iv) downstream clouds may be offered some protection from the oncoming flow as a result of being in the lee of an upstream cloud; v) the cycle of cloud compression and re-expansion is generally weaker when there are nearby neighbouring clouds; vi) the plasma $β$ in cloud material can vary rapidly as clouds collide with one another, but low values of $β$ are always transitory. This work is relevant to studies of multi-phase regions, where fast, low-density gas interacts with dense clouds, such as in circumstellar bubbles, supernova remnants, superbubbles and galactic winds.

preprint2013arXiv

Convergence of AMR and SPH simulations - I. Hydrodynamical resolution and convergence tests

We compare the results for a set of hydrodynamical tests performed with the AMR finite volume code, MG and the SPH code, SEREN. The test suite includes shock tube tests, with and without cooling, the non-linear thin-shell instability and the Kelvin-Helmholtz instability. The main conclusions are : (i) the two methods converge in the limit of high resolution and accuracy in most cases. All tests show good agreement when numerical effects (e.g. discontinuities in SPH) are properly treated. (ii) Both methods can capture adiabatic shocks and well-resolved cooling shocks perfectly well with standard prescriptions. However, they both have problems when dealing with under-resolved cooling shocks, or strictly isothermal shocks, at high Mach numbers. The finite volume code only works well at 1st order and even then requires some additional artificial viscosity. SPH requires either a larger value of the artificial viscosity parameter, alpha_AV, or a modified form of the standard artificial viscosity term using the harmonic mean of the density, rather than the arithmetic mean. (iii) Some SPH simulations require larger kernels to increase neighbour number and reduce particle noise in order to achieve agreement with finite volume simulations. However, this is partly due to the need to reduce noise that can corrupt the growth of small-scale perturbations. In contrast, instabilities seeded from large-scale perturbations do not require more neighbours and hence work well with standard SPH formulations and converge with the finite volume simulations. (iv) For purely hydrodynamical problems, SPH simulations take an order of magnitude longer to run than finite volume simulations when running at equivalent resolutions, i.e. when they both resolve the underlying physics to the same degree. This requires about 2-3 times as many particles as the number of cells.

preprint2013arXiv

Ram pressure stripping of the hot gaseous haloes of galaxies using the k-epsilon sub-grid turbulence model

We perform three dimensional hydrodynamic simulations of the ram pressure stripping of the hot extended gaseous halo of a massive galaxy using the k-epsilon sub-grid turbulence model at Mach numbers 0.9, 1.1 and 1.9. The k-epsilon model is used to simulate high Reynolds number flows by increasing the transport coefficients in regions of high turbulence. We find that the initial, instantaneous stripping is the same whether or not the k-epsilon model is implemented and is in agreement with the results of other studies. However the use of the k-epsilon model leads to five times less gas remaining after stripping by a supersonic flow has proceeded for 10 Gyr, which is more consistent with what simple analytic calculations indicate. Hence the continual Kelvin-Helmholtz stripping plays a significant role in the ram pressure stripping of the haloes of massive galaxies. To properly account for this, simulations of galaxy clusters will require the use of sub-grid turbulence models

preprint2013arXiv

The Responses of Magnetically Sub-Critical Cores to Shocks

An ideal magnetohydrodynamics (MHD) code with adaptive mesh refinement (AMR) was used to investigate the interactions of fast-mode shocks with self-gravitating, isothermal cores with mass-to-flux ratios that are somewhat below the minimum value required for gravitational collapse. We find that shock focussing produces colliding flows along the field lines that generate very high densities, even for relatively weak shocks. Self-gravity plays only a minor role in determining the highest density that is reached, but it does play a role in the subsequent evolution. The densities at comparable times differ by a factor of a few for shocks initially propagating perpendicularly or obliquely to the magnetic field in the ambient medium.

preprint2012arXiv

Magnetic Fields and Star Formation

Research performed in the 1950s and 1960s by Leon Mestel on the roles of magnetic fields in star formation established the framework within which he and other key figures have conducted subsequent investigations on the subject. This short tribute to Leon contains a brief summary of some, but not all, of his ground breaking contributions in the area. It also mentions of some of the relevant problems that have received attention in the last few years. The coverage is not comprehensive, and the authors have drawn on their own results more and touched more briefly on those of others than they would in a normal review. Theirs is a personal contribution to the issue honouring Leon, one of the truly great gentlemen, wits, and most insightful of astrophysicists.

preprint2012arXiv

Numerical simulations of shocks encountering clumpy regions

We present numerical simulations of the adiabatic interaction of a shock with a clumpy region containing many individual clouds. Our work incorporates a sub-grid turbulence model which for the first time makes this investigation feasible. We vary the Mach number of the shock, the density contrast of the clouds, and the ratio of total cloud mass to inter-cloud mass within the clumpy region. Cloud material becomes incorporated into the flow. This "mass-loading" reduces the Mach number of the shock, and leads to the formation of a dense shell. In cases in which the mass-loading is sufficient the flow slows enough that the shock degenerates into a wave. The interaction evolves through up to four stages: initially the shock decelerates; then its speed is nearly constant; next the shock accelerates as it leaves the clumpy region; finally it moves at a constant speed close to its initial speed. Turbulence is generated in the post-shock flow as the shock sweeps through the clumpy region. Clouds exposed to turbulence can be destroyed more rapidly than a similar cloud in an "isolated" environment. The lifetime of a downstream cloud decreases with increasing cloud-to-intercloud mass ratio. We briefly discuss the significance of these results for starburst superwinds and galaxy evolution.

preprint2012arXiv

Sputtering in oblique C-type shocks

We present the first results for the sputtering of grain mantles and cores obtained with self-consistent multifluid hydromagnetic models of C-type shocks propagating through dusty media. The threshold shock speed for mantle sputtering is about 10 km/s and is independent of density. The mantles are completely vapourised in shocks with speeds of 20-25 km/s. At such shock speeds core sputtering commences and gas-phase SiO forms. Core destruction is not total in any C-type shock because grains are not completely destroyed in shocks with speeds near the minimum speeds at which J-type shocks appear. Due to the density-dependence of the critical shock speed for this transition, higher gas-phase SiO fractional abundances are produced behind shocks propagating in lower density gas. For shock speeds near the threshold speeds for both core and mantle sputtering, sputtering is much greater for shock velocities at smaller angles relative to the upstream magnetic field. At higher shock speeds, the angular variation is still present but less pronounced.

preprint2011arXiv

Cosmic Rays and Radiative Instabilities

In the absence of magnetic fields and cosmic rays, radiative cooling laws with a range of dependences on temperature affect the stability of interstellar gas. For about four and a half decades, astrophysicists have recognised the importance of the thermal instablity for the formation of clouds in the interstellar medium. Even in the past several years, many papers have concerned the role of the thermal instability in the production of molecular clouds. About three and a half decades ago, astrophysicists investigating radiative shocks noticed that for many cooling laws such shocks are unstable. Attempts to address the effects of cosmic rays on the stablity of radiative media that are initially uniform or that have just passed through shocks have been made. The simplest approach to such studies involves the assumption that the cosmic rays behave as a fluid. Work based on such an approach is described. Cosmic rays have no effect on the stability of initially uniform, static media with respect to isobaric perturbations, though they do affect the stability of such media with respect to isentropic perturbations. The effect of cosmic rays on the stability of radiative shocked media depends greatly on the efficiency of the conversion of energy in accelerated cosmic rays into thermal energy in the thermalized fluid. If that efficiency is low, radiative cooling makes weak shocks propagating into upstream media with low cosmic-ray pressures more likely to be cosmic-ray dominated than adiabatic shocks of comparable strength. The cosmic-ray dominated shocks do not display radiative overstability. Highly efficient conversion of cosmic-ray energy into thermal energy leads shocked media to behave as they do when cosmic rays are absent.

preprint2010arXiv

Dusty Magnetohydrodynamics in Star Forming Regions

Star formation occurs in dark molecular regions where the number density of hydrogen nuclei, nH, exceeds 10^4 cc and the fractional ionization is 10^-7 or less. Dust grains with sizes ranging up to tenths of microns and perhaps down to tens of nanometers contain just under one percent of the mass. Recombination on grains is important for the removal of gas phase ions, which are produced by cosmic rays penetrating the dark regions. Collisions of neutrals with charged grains contribute significantly to the coupling of the magnetic field to the neutral gas. Consequently, the dynamics of the grains must be included in the magnetohydrodynamic models of large scale collapse, the evolution of waves and the structures of shocks important in star formation.

preprint2010arXiv

Shock-triggered formation of magnetically-dominated clouds. II. Weak shock-cloud interaction in three dimensions

To understand the formation of a magnetically dominated molecular cloud from an atomic cloud, we study the interaction of a weak, radiative shock with a magnetised cloud. The thermally stable warm atomic cloud is initially in static equilibrium with the surrounding hot ionised gas. A shock propagating through the hot medium then interacts with the cloud. We follow the dynamical evolution of the shocked cloud with a time-dependent ideal magnetohydrodynamic code. By performing the simulations in 3D, we investigate the effect of different magnetic field orientations including parallel, perpendicular and oblique to the shock normal. We find that the angle between the shock normal and the magnetic field must be small to produce clouds with properties similar to observed molecular clouds.

preprint2010arXiv

The Turbulent Destruction of Clouds - II. Mach Number Dependence, Mass-loss Rates, and Tail Formation

The turbulent destruction of a cloud subject to the passage of an adiabatic shock is studied. We find large discrepancies between the lifetime of the cloud and the analytical result of Hartquist et al. (1986). These differences appear to be due to the assumption in Hartquist et al. that mass-loss occurs largely as a result of lower pressure regions on the surface of the cloud away from the stagnation point, whereas in reality Kelvin-Helmholtz (KH) instabilities play a dominant role in the cloud destruction. We find that the true lifetime of the cloud (defined as when all of the material from the core of the cloud is well mixed with the intercloud material in the hydrodynamic cells) is about 6 times t_KHD, where t_KHD is the growth timescale for the most disruptive, long-wavelength, KH instabilities. These findings have wide implications for diffuse sources where there is transfer of material between hot and cool phases. The properties of the interaction as a function of Mach number and cloud density contrast are also studied. The interaction is milder at lower Mach numbers with the most marked differences occuring at low shock Mach numbers when the postshock gas is subsonic with respect to the cloud (i.e. M < 2.76). Material stripped off the cloud only forms a long "tail-like" feature if the density contrast of the cloud to the ambient medium, chi > 1e3.

preprint2009arXiv

Transient evolution of C-type shocks in dusty regions of varying density

Outflows of young stars drive shocks into dusty, molecular regions. Most models of such shocks assume that they are steady and propagating perpendicular to the magnetic field. Real shocks often violate both of these assumptions and the media through which they propagate are inhomogeneous. We use the code employed previously to produce the first time-dependent simulations of fast-mode, oblique C-type shocks interacting with density perturbations. We include a self-consistent calculation of the thermal and ionisation balances and a fluid treatment of grains. We identify features that develop when a multifluid shock encounters a density inhomogeneity to investigate whether any part of the precursor region ever behaves in a quasi-steady fashion. If it does the shock may be modelled approximately without solving the time-dependent hydromagnetic equations. Simulations were made for initially steady oblique C-type shocks encountering density inhomogeneities. For a semi-finite inhomogeneity with a density larger than the surrounding medium, a transmitted shock evolves from being J-type to a steady C-type shock on a timescale comparable to the ion-flow time through it. A sufficiently upstream part of the precursor of an evolving J-type shock is quasi-steady. The ion-flow timescale is also relevant for the evolution of a shock moving into a region of decreasing density. The models for shocks propagating into regions in which the density increases and then decreases to its initial value cannot be entirely described in terms of the results obtained for monotonically increasing and decreasing densities. For the latter model, the long-term evolution to a C-type shock cannot be approximated by quasi-steady models.

preprint2008arXiv

A cosmic-ray precursor model for a Balmer-dominated shock in Tycho's supernova remnant

We present a time-dependent cosmic-ray modified shock model for which the calculated H-alpha emissivity profile agrees well with the H-alpha flux increase ahead of the Balmer-dominated shock at knot g in Tycho's supernova remnant, observed by Lee et al (2007). The backreaction of the cosmic ray component on the thermal component is treated in the two-fluid approximation, and we include thermal particle injection and energy transfer due to the acoustic instability in the precursor. The transient state of our model that describes the current state of the shock at knot g, occurs during the evolution from a thermal gas dominated shock to a smooth cosmic-ray dominated shock. Assuming a distance of 2.3 kpc to Tycho's remnant we obtain values for the cosmic ray diffusion coefficient, the injection parameter, and the time scale for the energy transfer of 10^{24} cm^{2} s^{-1}, 4.2x10^{-3}, and 426 y, respectively. We have also studied the parameter space for fast (300 km s^{-1} - 3000 km s^{-1}), time-asymptotically steady shocks and have identified a branch of solutions, for which the temperature in the cosmic ray precursor typically reaches 2-6x10^{4} K and the bulk acceleration of the flow through the precursor is less than 10 km s^{-1}. These solutions fall into the low cosmic ray acceleration efficiency regime and are relatively insensitive to shock parameters. This low cosmic ray acceleration efficiency branch of solutions may provide a natural explanation for the line broadening of the H-alpha narrow component observed in non-radiative shocks in many supernova remnants.

preprint2008arXiv

Dust in Interstellar Clouds, Evolved Stars and Supernovae

Outflows of pre-main-sequence stars drive shocks into molecular material within 0.01 - 1 pc of the young stars. The shock-heated gas emits infrared, millimeter and submillimeter lines of many species including. Dust grains are important charge carriers and play a large role in coupling the magnetic field and flow of neutral gas. Some effects of the dust on the dynamics of oblique shocks began to emerge in the 1990s. However, detailed models of these shocks are required for the calculation of the grain sputtering contribution to gas phase abundances of species producing observed emissions. We are developing such models. Some of the molecular species introduced into the gas phase by sputtering in shocks or by thermally driven desorption in hot cores form on grain surfaces. Recently laboratory studies have begun to contribute to the understanding of surface reactions and thermally driven desorption important for the chemistry of star forming clouds. Dusty plasmas are prevalent in many evolved stars just as well as in star forming regions. Radiation pressure on dust plays a significant role in mass loss from some post-main-sequence stars. The mechanisms leading to the formation of carbonaceous dust in the stellar outflows are similar to those important for soot formation in flames. However, nucleation in oxygen-rich outflows is less well understood and remains a challenging research area. Dust is observed in supernova ejecta that have not passed through the reverse shocks that develop in the interaction of ejecta with ambient media. Dust is detected in high redshift galaxies that are sufficiently young that the only stars that could have produced the dust were so massive that they became supernovae. Consequently, the issue of the survival of dust in strong supernova shocks is of considerable interest.

preprint2003arXiv

The Evolution of Mass Loaded Supernova Remnants. II. Temperature Dependent Mass Injection Rates

We investigate the evolution of spherically symmetric supernova remnants in which mass loading takes place due to conductively driven evaporation of embedded clouds. Numerical simulations reveal significant differences between the evolution of conductively mass loaded and the ablatively mass loaded remnants studied in Paper I. A main difference is the way in which conductive mass loading is extinguished at fairly early times, once the interior temperature of the remnant falls below ~10 million K. Thus, at late times remnants that ablatively mass load are dominated by loaded mass and thermal energy, while those that conductively mass load are dominated by swept-up mass and kinetic energy. Simple approximations to the remnant evolution, complementary to those in Paper I, are given.