FLASH Code Research Articles

Active Galactic Nucleus Jet-inflated Bubbles as Possible Origin of Odd Radio Circles

Odd radio circles (ORCs) are newly discovered extragalactic radio objects with an unknown origin. In this work, we carry out three-dimensional cosmic-ray (CR) magnetohydrodynamic simulations using the FLASH code and predict the radio morphology of end-on active galactic nucleus (AGN) jet-inflated bubbles considering hadronic emission. We consider CR proton (CRp)-dominated jets as they tend to inflate oblate bubbles, promising to reproduce the large inferred sizes of the ORCs when viewed end-on. We find that powerful and long-duration CRp-dominated jets can create bubbles with similar sizes (∼300–600 kpc) and radio morphology (circular and edge-brightened) to the observed ORCs in low-mass (M vir ∼ 8 × 1012 − 8 × 1013 M ⊙) halos. Given the same amount of input jet energy, longer-duration (thus lower-power) jets tend to create larger bubbles since high-power jets generate strong shocks that carry away a significant portion of the jet energy. The edge-brightened feature of the observed ORCs is naturally reproduced due to efficient hadronic collisions at the interface between the bubbles and the ambient medium. We further discuss the radio luminosity, X-ray detectability, and the possible origin of such strong AGN jets in the context of galaxy evolution. We conclude that end-on CR-dominated AGN bubbles could be a plausible scenario for the formation of ORCs.

AHKASH: a new Hybrid particle-in-cell code for simulations of astrophysical collisionless plasma.

We introduce Astrophysical Hybrid-Kinetic simulations with the flash code ([Formula: see text]) - a new Hybrid particle-in-cell (PIC) code developed within the framework of the multiphysics code flash. The new code uses a second-order accurate Boris integrator and a predictor-predictor-corrector algorithm for advancing the Hybrid-kinetic equations, using the constraint transport method to ensure that magnetic fields are divergence-free. The code supports various interpolation schemes between the particles and grid cells, with post-interpolation smoothing to reduce finite particle noise. We further implement a [Formula: see text] method to study instabilities in weakly collisional plasmas. The new code is tested on standard physical problems such as the motion of charged particles in uniform and spatially varying magnetic fields, the propagation of Alfvén and whistler waves, and Landau damping of ion acoustic waves. We test different interpolation kernels and demonstrate the necessity of performing post-interpolation smoothing. We couple the turbgen turbulence driving module to the new Hybrid PIC code, allowing us to test the code on the highly complex physical problem of the turbulent dynamo. To investigate steady-state turbulence with a fixed sonic Mach number, it is important to maintain isothermal plasma conditions. Therefore, we introduce a novel cooling method for Hybrid PIC codes and provide tests and calibrations of this method to keep the plasma isothermal. We describe and test the 'hybrid precision' method, which significantly reduces (by a factor [Formula: see text]) the computational cost, without compromising the accuracy of the numerical solutions. Finally, we test the parallel scalability of the new code, showing excellent scaling up to 10,000 cores.

MHD modeling of the molecular filament evolution

We perform numerical magnetohydrodynamic (MHD) simulations of the gravitational collapse and fragmentation of a cylindrical molecular cloud with the help of the FLASH code. The cloud collapses rapidly along it’s radius without any signs of fragmentation in the simulations without magnetic field. The radial collapse of the cloud is stopped by the magnetic pressure gradient in the simulations with parallel magnetic field. Cores with high density form at the cloud’s edges during further evolution. The core densities are n ≈ 1.7×108 and 2×10-7 cm–3 in the cases with initial magnetic field strengths B = 1.9×10-4 and 6×10-4 G, respectively. The cores move toward the cloud’s centre with supersonic speeds |vz| = 3.6 and 5.3 km/s. The sizes of the cores along the cloud’s radius and cloud’s main axis are dr = 0.0075 pc and dz = 0.025 pc, dr = 0.03 pc and dz = 0.025 pc, respectively. The masses of the cores increase during the filament evolution and lie in range of ≈(10-20)Me. According to our results, the cores observed at the edges of molecular filaments can be a result of the filament evolution with parallel magnetic field.

Properties of molecular clumps and cores in colliding magnetized flows

ABSTRACT We simulate the formation of molecular clouds in colliding flows of warm neutral medium with the adaptive mesh refinement code flash in eight simulations with varying initial magnetic field strength, between 0.01–5 μG. We include a chemical network to treat heating and cooling and to follow the formation of molecular gas. The initial magnetic field strength influences the fragmentation of the forming cloud because it prohibits motions perpendicular to the field direction and hence impacts the formation of large-scale filamentary structures. Molecular clump and core formation occurs anyhow. We identify 3D clumps and 3D cores, which are defined as connected, CO-rich regions. Additionally, 3D cores are heavily shielded. While we do not claim those 3D objects to be directly comparable to observations, this enables us to analyse their full virial state. With increasing field strength, we find more fragments with a smaller average mass; yet the dynamics of the forming clumps and cores only weakly depends on the initial magnetic field strength. The molecular clumps are mostly unbound, probably transient objects, which are weakly confined by ram pressure or thermal pressure, indicating that they are swept up by the turbulent flow. They experience significant fluctuations in the mass flux through their surface, such that the Eulerian reference frame shows a dominant time-dependent term due to their indistinct nature. We define the cores to encompass highly shielded molecular gas. Most cores are in gravitational-kinetic equipartition and are well described by the common virial parameter $\alpha _\mathrm{vir}$, while some undergo minor dispersion by kinetic surface effects.

Magnetohydrodynamics simulation of magnetic reconnection process based on the laser-driven Helmholtz capacitor-coil targets

Magnetic reconnection is an important rapid energy release mechanism in astrophysics. Magnetic energy can be effectively converted into plasma kinetic energy, thermal energy, and radiation energy. This study is based on the magnetohydrodynamics simulation method and utilizes the FLASH code to investigate the laser-driven magnetic reconnection physical process of the Helmholtz capacitor-coil target. The simulation model incorporates the laser driving effect, and the external magnetic field consistent with the Helmholtz capacitor-coil target is written in. This approach achieves a magnetic reconnection process that is more consistent with the experiment. By changing the resistivity, subtle differences in energy conversion during the evolution of magnetic reconnection are observed. Under conditions of low resistivity, there is a more pronounced increase in the thermal energy of ions compared to other energy components. In simulations with high resistivity, the increase in electrons thermal energy is more prominent. The simulation gives the evolution trajectory of magnetic reconnection, which is in good agreement with the experimental results. This has important reference value for experimental research on the low-β magnetic reconnection.

Stability of perpendicular magnetohydrodynamic shocks in materials with ideal and nonideal equationsof state.

Magnetized target fusion approach to inertial confinement fusion involves the formation of strong shocks that travel along a magnetized plasma. Shocks, which play a dominant role in thermalizing the upstream kinetic energy generated in the implosion stage, are seldom free from perturbations, and they wrinkle in response to upstream or downstream disturbances. In Z-pinch experiments, significant plasma instability mitigation was observed with pre-embedded axial magnetic fields. To isolate effects, in this work we theoretically study the impact of perpendicular magnetic fields on the planar shock dynamics for different equationsof state. For fast magnetosonic shocks in ideal gases, it was found that the magnetic field amplifies the intensity of the perturbations when γ>2 or it weakens them when γ<2. Weak shocks have been found to be stable regardless of the magnetic plasma intensity and gas compressibility; however, for sufficiently strong shocks the magnetic fields can promote a neutral stability/SAE at the shock if the adiabatic index is higher than 1+sqrt[2]. Results have been validated with numerical simulations performed with the FLASH code.

Energy transfer and scale dynamics in 2D and 3D laser-driven jets

We demonstrate a methodology for diagnosing the multiscale dynamics and energy transfer in complex HED flows with realistic driving and boundary conditions. The approach separates incompressible, compressible, and baropycnal contributions to energy scale-transfer and quantifies the direction of these transfers in (generalized) wavenumber space. We use this to compare the kinetic energy (KE) transfer across scales in simulations of 2D axisymmetric vs fully 3D laser-driven plasma jets. Using the FLASH code, we model a turbulent jet ablated from an aluminum cone target in the configuration outlined by Liao et al. [Phys. Plasmas, 26 032306 (2019)]. We show that, in addition to its well known bias for underestimating hydrodynamic instability growth, 2D modeling suffers from significant spurious energization of the bulk flow by a turbulent upscale cascade. In 2D, this arises as vorticity and strain from instabilities near the jet's leading edge transfer KE upscale, sustaining a coherent circulation that helps propel the axisymmetric jet farther (≈25% by 3.5 ns) and helps keep it collimated. In 3D, the coherent circulation and upscale KE transfer are absent. The methodology presented here may also help with inter-model comparison and validation, including future modeling efforts to alleviate some of the 2D hydrodynamic artifacts highlighted in this study.

A tool and a methodology to use macros for abstracting variations in code for different computational demands

Scientific software used on high-performance computing platforms is in a phase of transformation because of the combined increase in the heterogeneity and complexity of models and hardware platforms. Having separate implementations for different platforms can easily lead to combinatorial explosions; therefore, the computational science community has been looking for mechanisms to express code through abstractions that can be specialized for different platforms. Most existing approaches use template meta-programming in C++, and are, therefore language specific. We have developed a tool that uses customized expansion of macros to mimic some of C++ behaviour in other languages. It enables unification of any code variants that may be necessary to run efficiently on different target architectures and different computational environments through use of macros with multiple alternative definitions and ability to arbitrate on definition selection for expansion. Combined with two other tools, a custom runtime, and a user specified recipe translator, our custom macroprocessor becomes a part of an overall performance portability solution that does not depend on any specific programming language. We also use macros as code-shorthand that lets code snippets become building blocks that allow variations in control flow to explore performance options. We demonstrate use of macros in Flash-X, a multiphysics multicomponent code with many Fortran legacy components derived from an earlier community code FLASH.

How do supernova remnants cool? – I. Morphology, optical emission lines, and shocks

ABSTRACT Supernovae (SNe) inject ∼1051 erg in the interstellar medium, thereby shocking and heating the gas. A substantial fraction of this energy is later lost via radiative cooling. We present a post-processing module for the flash code to calculate the cooling radiation from shock-heated gas using collisional excitation data from mappings v. When applying this tool to a simulated SN remnant (SNR), we find that most energy is emitted in the EUV. However, optical emission lines ([O iii], [N ii], [S ii], H α, H β) are usually best observable. Our shock detection scheme shows that [S ii] and [N ii] emissions arise from the thin shell surrounding the SNR, while [O iii], H $\rm \alpha$, and H $\rm \beta$ originate from the volume-filling hot gas inside the SNR bubble. We find that the optical emission lines are affected by the SNR’s complex structure and its projection on to the plane of the sky because the escaping line luminosity can be reduced by 10–80 per cent due to absorption along the line of sight. Additionally, the subtraction of contaminating background radiation is required for the correct classification of an SNR on the oxygen or sulphur BPT diagrams. The electron temperature and density obtained from our synthetic observations match well with the simulation but are very sensitive to the assumed metallicity.

Tree-based solvers for adaptive mesh refinement code FLASH – III: a novel scheme for radiation pressure on dust and gas and radiative transfer from diffuse sources

ABSTRACT Radiation is an important contributor to the energetics of the interstellar medium, yet its transport is difficult to solve numerically. We present a novel approach towards solving radiative transfer of diffuse sources via backwards ray tracing. Here, we focus on the radiative transfer of infrared radiation and the radiation pressure on dust. The new module, TreeRay/RadPressure, is an extension to the novel radiative transfer method TreeRay implemented in the grid-based Magneto-Hydrodynamics code Flash. In TreeRay/RadPressure, every cell and every star particle is a source of infrared radiation. We also describe how gas, dust, and radiation are coupled via a chemical network. This allows us to compute the local dust temperature in thermal equilibrium, leading to a significantly improvement over the classical grey approximation. In several tests, we demonstrate that the scheme produces the correct radiative intensities as well as the correct momentum input by radiation pressure. Subsequently, we apply our new scheme to model massive star formation from a collapsing, turbulent core of 150 M⊙. We include the effects of both, ionizing and infrared radiation on the dynamics of the core. We find that the newborn massive star prevents fragmentation in its proximity due to radiative heating. Over time, dust and radiation temperature equalize, while the gas temperature can be either warmer due to shock heating or colder due to insufficient dust–gas coupling. Compared to gravity, the effects of radiation pressure are insignificant for the stellar mass on the simulated time-scale in this work.

Intricate structure of the plasma Rayleigh–Taylor instability in shock tubes

Spikes and bubbles grow on unstable interfaces that are accelerated in high-energy-density conditions. If a shock propagates ahead of the interface, the plasma can be heated to extreme conditions where conduction and radiation fluxes influence the hydrodynamics. For example, a National Ignition Facility experiment found reduced single-mode nonlinear mixed-width growth in conditions scaled from a supernova explosion [Kuranz et al., Nat. Commun. 9, 1564 (2018)]. We present high-resolution two-dimensional radiation hydrodynamic simulations with the Flash code that quantitatively reproduce the experiment. Radiative fluxes are primarily responsible for ablating the spike and removing the mushroom caps. The ablated plasma increases the mixed mass and forms a low-density halo with spikes forming in both directions. This is considerably more complex than classical instability. The halo is sensitive to ablative physics, so radiographing it may aid in the verification of energy transport modeling. Although ablation changes the spike shape, it has little effect on the overall mixed width for these parameters. This is because ablation enhances the bubble velocity but it has the opposite effect on the spike. The radiation transport instead suppresses the growth via increasing the shocked foam density, thus decreasing the Atwood number. A terminal velocity model including the rarefaction expansion agrees with the experimental mixed-width growth.

Evolution and feedback of AGN jets of different cosmic ray composition

ABSTRACT Jet feedback from active galactic nuclei (AGN) is one of the most promising mechanisms for suppressing cooling flows in cool-core clusters. However, the composition of AGN jets and bubbles remains uncertain; they could be thermally dominated, or dominated by cosmic ray proton (CRp), cosmic ray electron (CRe), or magnetic energy. In this work, we investigate the evolution and feedback effects of CRp and CRe dominated jets by conducting 3D magnetohydrodynamic simulations of AGN jet-inflated bubbles in the intracluster medium using the FLASH code. We present the evolution of their energies, dynamics, and heating, and model their expected cavity-power versus radio-luminosity relation (Pcav–LR). We find that bubbles inflated by CRe dominated jets follow a very similar dynamical evolution to CRp dominated bubbles even though CRe within bubbles suffer significantly stronger synchrotron and inverse-Compton cooling. This is because, as CRe lose their energy, the jet-inflated bubbles quickly become thermally dominated within ∼30 Myr. Their total energy stops decreasing with CR energy and evolves similarly to CRp dominated bubbles. The ability of CRe and CRp dominated bubbles to heat the intracluster medium is also comparable; the cold gas formed via local thermal instabilities is well suppressed in both cases. The CRp and CRe bubbles follow different evolutionary trajectories on the Pcav–LR plane, but the values are broadly consistent with observed ranges for FR-I sources. We also discuss observational techniques that have potential for constraining the composition of AGN jets and bubbles.

Impact of rotation on the multimessenger signatures of a hadron-quark phase transition in core-collapse supernovae

We study the impact of rotation on the multimessenger signals of core-collapse supernovae (CCSNe) with the occurrence of a first-order hadron-quark phase transition (HQPT). We simulate CCSNe with the \texttt{FLASH} code starting from a 20~$M_\odot$ progenitor with different rotation rates, and using the RDF equation of state from \textit{Bastian} 2021 that prescribes the HQPT. Rotation is found to delay the onset of the HQPT and the resulting dynamical collapse of the protocompact star (PCS) due to the centrifugal support. All models with the HQPT experience a second bounce shock which leads to a successful explosion. The oblate PCS as deformed by rotation gives rise to strong gravitational-wave (GW) emission around the second bounce with a peak amplitude larger by a factor of $\sim10$ than that around the first bounce. The breakout of the second bounce shock at the neutrinosphere produces a $\bar{\nu}_e$-rich neutrino burst with a luminosity of serveral 10$^{53}$~erg~s$^{-1}$. In rapidly rotating models the PCS pulsation following the second bounce generates oscillations in the neutrino signal after the burst. In the fastest rotating model with the HQPT, a clear correlation is found between the oscillations in the GW and neutrino signals immediately after the second bounce. In addition, the HQPT-induced collapse leads to a jump in the ratio of rotational kinetic energy to gravitational energy ($\beta$) of the PCS, for which persistent GW emission may arise due to secular nonaxisymmetric instabilities.

The effect of saturated thermal conduction on clouds in a hot plasma

ABSTRACT We numerically investigate the internal evolution of multiphase clouds, which are at rest with respect to an ambient, highly ionized medium (HIM) representing the hot component of the circumgalactic medium. Time-dependent saturated thermal conduction and its implications like condensation rates and mixing efficiency are assessed in multiphase clouds. Our simulations are carried out by using the adaptive mesh refinement code Flash. . The model clouds are initially in both hydrostatic and thermal equilibrium and are in pressure balance with the HIM. Thus, they have steep gradients in both temperature and density at the interface to HIM leading to non-negligible thermal conduction. Several physical processes are considered numerically or semi-analytically: thermal conduction, radiative cooling and external heating of gas, self-gravity, mass diffusion, and dissociation of molecules and ionization of atoms. It turns out that saturated thermal conduction triggers a continuous condensation irrespective of cloud mass. Dynamical interactions with ambient HIM all relate to the radial density gradient in the clouds: (1) mass flux due to condensation is the higher the more homogeneous the clouds are; (2) mixing of condensed gas with cloud gas is easier in low-mass clouds, because of their shallower radial density gradient; and thus (3) accreted gas is distributed more efficiently. A distinct and sub-structured transition zone forms at the interface between cloud and HIM, which starts at smaller radii and is much narrower as deduced from analytical theory.

Three-dimensional simulations of reshocked inclined Richtmyer-Meshkov instability: Effects of initial perturbations

A Mach 1.55 shock and subsequent reshock interacting with predominantly single-mode and multi-mode interfaces between nitrogen and carbon dioxide is studied through three-dimensional adaptive mesh refinement simulations using the FLASH code. The simulations are directly compared with the PLIF/PIV experimental measurements to validate the computational code. Development of instability after first shock and chaotic turbulent mixing after reshock are discussed for the entire shock tube domain.

Mitigating laser imprint with a foam overcoating

In direct-drive inertial confinement fusion, laser imprint can cause areal density perturbations on the target shell that seed the Rayleigh–Taylor instability and further degrade the implosion. To mitigate the effect of laser imprint, a foam overcoating layer outside the target shell has been suggested to increase the thermal smoothing of the conduction region (between the ablation front and the critical density surface) and mass ablation of the ablation front. In this paper, we use a two-dimensional radiation hydrodynamic code FLASH to investigate the laser imprint mitigation performance and find other physical mechanisms of foam overcoatings. First, radiation ablation dynamically modulates density distribution not only to increase the frequency of the perturbed ablation front oscillation but also to decrease the amplitude of the oscillation. Second, a larger length of the shocked compression region reduces the amplitude of the perturbed shock front oscillation. The areal density perturbations decrease with the decrease in the perturbations of the ablation front and shock front. Based on the abovementioned physical mechanisms, we propose the optimal ranges of foam parameters to mitigate laser imprint with the aid of dimensional analysis: the foam thickness is about two to three times that of the perturbation wavelength, and the foam density is about 1/2–3/2 times that of the critical density.

Formation of the Asymmetric Accretion Disk from Stellar Wind Accretion in an S-type Symbiotic Star

The accretion process in a typical S-type symbiotic star, targeting AG Draconis, is investigated through 3D hydrodynamical simulations using the FLASH code. Regardless of the wind velocity of the giant star, an accretion disk surrounding the white dwarf is always formed. In models where the wind is faster than the orbital velocity of the white dwarf, the disk size and accretion rate are consistent with the predictions under Bondi–Hoyle–Lyttleton (BHL) conditions. In slower-wind models, unlike the BHL predictions, the disk size does not grow, and the accretion rate increases to a considerably higher level, up to >20% of the mass-loss rate of the giant star. The accretion disk in our fiducial model is characterized by a flared disk with a radius of 0.16 au and a scale height of 0.03 au. The disk mass of ∼5 × 10−8 M ⊙ is asymmetrically distributed, with the density peak toward the giant star being about 50% higher than the density minimum in the disk. Two inflowing spiral features are clearly identified, and their relevance to the azimuthal asymmetry of the disk is pointed out. The flow in the accretion disk is found to be sub-Keplerian, at about 90% of the Keplerian speed, which indicates a caveat of overestimating the O vi emission region from the spectroscopy of Raman-scattered O vi features at 6825 and 7082 Å.

Particle-in-cell modeling of plasma jet merging in the large-Hall-parameter regime

The merging process of magnetized plasma jets with parameters relevant to the plasma-jet-driven magneto-inertial fusion (PJMIF) design and the plasma liner experiment (PLX) is modeled by fully kinetic particle-in-cell (PIC) simulations in one and two spatial dimensions. The modified two-stream instability is identified to be the main mechanism responsible for stopping the plasma jets and preventing species interpenetration. The electron and ion Hall parameters of the merged plasma are greater than unity, and the plasma β is close to unity, which is the desired characteristic of planned experiments at PLX. Our 2D PIC simulations validate the results of the radiation magneto-hydrodynamics code FLASH, which will be the primary tool for modeling various stages of future PJMIF experiments.

Code-to-code comparison and validation of the radiation-hydrodynamics capabilities of the FLASH code using a laboratory astrophysical jet

The potential for laser-produced plasmas to yield fundamental insights into high energy density physics (HEDP) and deliver other useful applications can sometimes be frustrated by uncertainties in modeling the properties and behavior of these plasmas using radiation-hydrodynamics codes. In an effort to overcome this and to corroborate the accuracy of the HEDP capabilities in the publicly available FLASH radiation-hydrodynamics code, we present detailed code-to-code comparisons between FLASH and the HYDRA code developed at Lawrence Livermore National Laboratory using previously published HYDRA simulations from Grava et al. [Phys. Rev. E 78, 016403 (2008)]. That study describes a laser experiment that produced a jet-like feature that the authors compare to astrophysical jets. Importantly, the Grava et al. [Phys. Rev. E 78, 016403 (2008)] experiment included detailed x-ray interferometric measurements of electron number densities and a time-integrated measurement of the soft x-ray spectrum. Despite markedly different methods for treating the computational mesh, and different equations of state and opacity models, the FLASH results resemble the results from HYDRA and, most importantly, the experimental measurements of electron density. Having validated the FLASH code in this way, we use the code to further investigate and understand the formation of the jet seen in the Grava et al. [Phys. Rev. E 78, 016403 (2008)] experiment and discuss its relation to the Wan et al. [Phys. Rev. E 55, 6293 (1997)] experiment at the NOVA laser.

Numerical investigation of laser-driven shock interaction with a deformable particle

A laser-driven shock propagating through an isolated particle embedded in a plastic (CH) target was studied using the radiation-hydrodynamic code FLASH. Preliminary simulations using IONMIX equations of state (EOS) showed significant differences in the shock Hugoniot of aluminum compared to experimental data in the low-pressure regime [O(10) GPa], resulting in higher streamwise compression and deformation of an aluminum particle. Hence, a simple modification to the ideal gas EOS was developed and employed to describe the target materials and examine the particle dynamics. The evolution of the pressure field demonstrated a complex wave interaction, resulting in a highly unsteady particle drag which featured two drag minima due to shock focusing at the rear end of the particle and rarefaction stretching due to laser shut-off. Although ∼30% lateral expansion and ∼25% streamwise compression were observed, the aluminum particle maintained considerable integrity without significant distortion. Additional simulations examined the particle response for a range of particle densities, sizes, and acoustic impedances. The results revealed that lighter particles such as aluminum gained significant momentum, reaching up to ∼96% of the shocked CH's speed, compared to ∼29% for the heavier tungsten particles. Despite the differences seen in the early stage of shock interaction, particles with varying acoustic impedances ultimately reached the same peak velocity. This identified particle-to-host density ratio is an important factor in determining the inviscid terminal velocity of the particle. In addition, the modified EOS model presented in this study could be used to approximate solid materials in hydrocodes that lack material strength models.