Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “FLASH”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 577 records · Page 32
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Evolution and feedback of AGN jets of different cosmic ray composition

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 (P cav –L R ). 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 P cav –L R 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.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

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.

Weis, M. (ORCID:0000000256838860)↗
DOE OSTI2025

2D end-to-end modelling of kilonovae from binary neutron star merger remnants

We investigate the kilonova emission resulting from outflows produced in a 3D general-relativistic magnetohydrodynamic (GRMHD) simulation of a hypermassive neutron star (HMNS) remnant. We map the outflows into the flash hydrodynamics code to model their expansion in axisymmetry, and study the effects of employing different r-process heating rates. Except for the highest heating rate prescription, we find no significant differences with respect to overall ejecta dynamics and morphology compared to the simulation without heating. Once homologous expansion is attained, typically after ~2s for these ejecta, we map the outflows to the sedona radiative transfer code and compute the spectral evolution of the kilonova and broad-band light curves in various Legacy Survey of Space and Time (LSST) bands. The kilonova properties depend on the remnant lifetime, with peak luminosities and peak time-scales increasing for longer lived remnants that produce more massive ejecta. For all models, there is a strong dependence of both the bolometric and broad-band light curves on the viewing angle. While the short-lived (12 ms) remnant produces higher luminosities when viewed from angles closer to the pole, longer lived remnants (240 ms and 2.5 s) are more luminous when viewed from angles closer to the equator. Our results highlight the importance of self-consistent, long-term modelling of merger ejecta, and taking viewing-angle dependence into account when interpreting observed kilonova light curves. We find that magnetized outflows from an HMNS – if it survives long enough – could explain blue kilonovae, such as the blue emission seen in AT2017gfo.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Neutrino heating in 1D, 2D, and 3D core-collapse supernovae: characterizing the explosion of high-compactness stars

Massive stars can end their lives with a successful supernova explosion (leaving behind a neutron star or, more rarely, a black hole), or a failed explosion that leaves behind a black hole. The density structure of the pre-collapse progenitor star already encodes much of the information regarding the outcome and properties of the explosion. However, the complexity of the collapse and subsequent shock expansion phases prevents drawing a straightforward connection between the pre-collapse and post-explosion properties. In order to derive such a connection several explodability studies have been performed in recent years. However, different studies can predict different explosion outcomes. In this article, we show how compactness, which is related to the average density of the star’s core, has an important role in determining the efficiency of neutrino heating, and therefore the outcome of the explosion. Commonly, high-compactness progenitors are assumed to yield failed explosions, due to their large mass accretion rates, preventing the shock from expanding. We show by analysing ~150 2D flash and F ornax simulations and 20 3D F ornax simulations that this is not the case. Instead, due to the rapid increase of neutrino heating with compactness, high-compactness progenitors lead to successful shock revival. We also show that 1D+ simulations that include v-driven convection using a mixing-length theory approach correctly reproduce this trend. Finally, we compare 1D+ models, which we show can reproduce some aspects of multi-D simulations with reasonable accuracy, with other widely used 1D models in the literature.

neutrinos↗
DOE OSTI2024

Development of a Data Overflow Protection System for Super-Kamiokande to Maximize Data from Nearby Supernovae

Neutrinos from very nearby supernovae, such as Betelgeuse, are expected to generate more than ten million events over 10 s in Super-Kamokande (SK). At such large event rates, the buffers of the SK analog-to-digital conversion board (QBEE) will overflow, causing random loss of data that are critical for understanding the dynamics of the supernova explosion mechanism. In order to solve this problem, two new data-acquisition (DAQ) modules were developed to aid in the observation of very nearby supernovae. The first of these, the SN module, is designed to save only the number of hit photomultiplier tubes during a supernova burst and the second, the Veto module, prescales the high-rate neutrino events to prevent the QBEE from overflowing based on information from the SN module. In the event of a very nearby supernova, these modules allow SK to reconstruct the time evolution of the neutrino event rate from beginning to end using both QBEE and SN module data. This paper presents the development and testing of these modules together with an analysis of supernova-like data generated with a flashing laser diode. We demonstrate that the Veto module successfully prevents DAQ overflows for Betelgeuse-like supernovae as well as the long-term stability of the new modules. During normal running the Veto module is found to issue DAQ vetos a few times per month resulting in a total dead-time less than 1 ms, and does not influence ordinary operations. Additionally, using simulation data we find that supernovae closer than 800 pc will trigger the Veto module, resulting in a prescaling of the observed neutrino data.

F20 Instrumentation and technique↗
DOE OSTI2021

Response of extreme precipitation to uniform surface warming in quasi-global aquaplanet simulations at high resolution

Projections of precipitation extremes in simulations with global climate models are very uncertain in the tropics, in part because of the use of parameterizations of deep convection and model deficiencies in simulating convective organization. Here, we analyse precipitation extremes in high-resolution simulations that are run without a convective parameterization on a quasi-global aquaplanet. The frequency distributions of precipitation rates and precipitation cluster sizes in the tropics of a control simulation are similar to the observed distributions. In response to climate warming, 3 h precipitation extremes increase at rates of up to 9 % K - 1 in the tropics because of a combination of positive thermodynamic and dynamic contributions. The dynamic contribution at different latitudes is connected to the vertical structure of warming using a moist static stability. When the precipitation rates are first averaged to a daily timescale and coarse-grained to a typical global climate-model resolution prior to calculating the precipitation extremes, the response of the precipitation extremes to warming becomes more similar to what was found previously in coarse-resolution aquaplanet studies. However, the simulations studied here do not exhibit the high rates of increase of tropical precipitation extremes found in projections with some global climate models. This article is part of a discussion meeting issue ‘Intensification of short-duration rainfall extremes and implications for flash flood risks’.

Science & Technology - Other Topics↗
DOE OSTI2025

Demonstration of a code coupling framework for modeling beam-collimator impacts in the advanced photon source

The high-brightness beams being produced in current and future accelerators present new machine protection concerns with the potential for high-energy-density (HED) conditions ( >100 J/mm 3 ) in beam-intercepting components. Simulating HED conditions in accelerators requires utilizing a suite of physics codes for particle dynamics, particle-matter interactions, and hydrodynamics. This paper describes a method of coupling the codes elegant, fluka, and flash to simulate the effects of a rapid beam loss in the advanced photon source storage ring and the resulting interaction of the beam and collimators. This paper expands previous work [J. Dooling et al., Collimator irradiation studies at the advanced photon source, in Proceedings of the IBIC-2023 (2023), pp. 245–249] by introducing a definition of the evolving geometry of the collimator surface as well as providing methods for simulating the absorption of synchrotron radiation and tracking shower particles produced during beam strikes. We demonstrate this framework by simulating machine conditions of the APS ring before and after its recent upgrade. Simulation results are compared with observed damage to collimators and test samples taken from the APS ring.

Accelerator/storage ring control systems↗
DOE OSTI2025

Ratchet effects in cyclic pattern formation systems with competing interactions

Ratchet effects can appear for particles interacting with an asymmetric potential under ac driving or for a thermal system in which a substrate is periodically flashed. Here, we show that a different type of collective ratchet effect can arise for a pattern-forming system coupled to an asymmetric substrate when the interaction potential between the particles is periodically oscillated in order to cycle the system through different patterns. We consider particles with competing short-range attraction and long-range repulsion subjected to time-dependent oscillations of the ratio between the attractive and repulsive interaction terms, which causes the system to cycle periodically between crystal and bubble states. In the presence of the substrate, this system exhibits both a positive and a reversed ratchet effect, and we show that there is a maximum in the ratchet efficiency as a function of interaction strength, ac drive frequency, and particle density. In conclusion, our results could be realized for a variety of pattern-forming systems on asymmetric substrates where the pattern type or particle interactions can be oscillated.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Experimental evidence of early-time saturation of the ion-Weibel instability in counterstreaming plasmas of CH, Al, and Cu

The collisionless ion-Weibel instability is a leading candidate mechanism for the formation of collisionless shocks in many astrophysical systems, where the typical distance between particle collisions is much larger than the system size. Multiple laboratory experiments aimed at studying this process utilize laser-driven (I≳10 15 W/cm 2 ), counterstreaming plasma flows (V≲ 2000 km/s) to create conditions unstable to Weibel-filamentation and growth. This technique intrinsically produces temporally varying plasma conditions at the midplane of the interaction where Weibel-driven B fields are generated and studied. Experiments discussed herein demonstrate robust formation of Weibel-driven B fields under multiple plasma conditions using CH, Al, and Cu plasmas. Linear theory based on benchmarked radiation-hydrodynamic FLASH calculations is compared with Fourier analyses of proton images taken ~5–6 linear growth times into the evolution. Our new analyses presented here indicate that the low-density, high-velocity plasma-conditions present during the first linear-growth time (~300–500 ps) sets the spectral characteristics of Weibel filaments during the entire evolution. It is shown that the dominant wavelength (~300 μm) at saturation persists well into the nonlinear phase, consistent with theory under these experimental conditions. However, estimates of B-field strength, while difficult to determine accurately due to the path-integrated nature of proton imaging, are shown to be in the ~10–30 T range, an order of magnitude above the expected saturation limit in homogenous plamas but consistent with enhanced B fields in the midplane due to temporally varying plasma conditions in experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Stability of perpendicular magnetohydrodynamic shocks in materials with ideal and nonideal equations of 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 equations of state. For fast magnetosonic shocks in ideal gases, it was found that the magnetic field amplifies the intensity of the perturbations when $\gamma$ > 2 or it weakens them when $\gamma$ < 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}$. Finally, results have been validated with numerical simulations performed with the FLASH code.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Shock wave formation in radiative plasmas

The temporal evolution of weak shocks in radiative media is theoretically investigated in this work. The structure of radiative shocks has traditionally been studied in a stationary framework. Their systematic classification is complex because layers of optically thick and thin regions alternate to form a radiatively-driven precursor and a temperature-relaxation layer, between which the hydrodynamic shock is embedded. In this work, we analyze the formation of weak shocks when two radiative plasmas with different pressures are put in contact. Applying a reductive perturbative method yields a Burgers-type equation that governs the temporal evolution of the perturbed variables including the radiation field. The conditions upon which optically thick and thin solutions exist have been derived and expressed as a function of the shock strength and Boltzmann number. Below a certain Boltzmann number threshold, weak shocks always become optically thick asymptotically in time, while thin solutions appear as transitory structures. The existence of an optically thin regime is related to the presence of an overdense layer in the compressed material. Scaling laws for the characteristic formation time and shock width are provided for each regime. The theoretical analysis is supported by FLASH simulations, and a comprehensive testcase has been designed to benchmark radiative hydrodynamic codes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

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

The effect of initial perturbations on the evolution of the inclined Richtmyer-Meshkov turbulent mixing layer before and after reshock initiated by a shock wave with Mach number 1.55 is investigated through three-dimensional (3D) simulations using the flash code. The 3D simulations aim to reproduce both predominantly single-mode and multimode interfaces between light and heavy gases (N 2 -CO 2 , Atwood number, A≈0.22; amplitude to wavelength ratio of 0.088) which were created in an inclined shock tube facility to analyze the effects of initial conditions on mixing development in the entire flow field. The two-dimensional center slices of 3D simulations are compared with the experimental results to validate the computational code. Mixing width, mixed mass, mixed-mass thickness, and circulation in addition to concentration fields are shown to be in good agreement with the experimental data. The three-dimensional density and vorticity fields are first presented to qualitatively describe the flow behavior before and after reshock. Several measured density/velocity-related quantities indicate that the growth of the mixing material is strongly dependent on initial conditions. Before reshock and at early times after reshock, flow is clearly maintaining the memory of initial perturbations. However, at late time after reshock, although the large wavelength feature still dominates the flow motion, and the morphology of the two different interfaces indicates several differences, by breakdown of large-scale coherent structures to much finer scales, the memory of small scales of the multimode initial perturbation is not as clear as pre-reshock. Regarding three-dimensionality of the flow, before reshock in the multimode case, the baroclinic vorticity production, circulation, turbulent kinetic energy, and turbulent mass flux suggest that the small-scale roll-up features along the large inclined wavelength quickly evolves in all three dimensions. The coherent vortex tubes break down to smaller wormlike vortex structures, and turbulent fluctuations in the out-of-plane dimension are comparable to the spanwise direction. After reshock, this three-dimensionality of mixing growth was observed in the flow for both initial conditions. The results of this work represent a significant extension of previous computational studies performed on this specific topic. A different code with a different numerical method is validated through comparison with the experimental data. The initial perturbations are directly measured from the experimental results. Moreover, the entire three-dimensional experimental shock tube domain is simulated, and more quantities are investigated to understand the mixing mechanism and instability evolution in all three dimensions.

42 ENGINEERING↗
DOE OSTI2024

Supernova Axions Convert to Gamma Rays in Magnetic Fields of Progenitor Stars

It has long been established that axions could have been produced within the nascent proto-neutron star formed following the type II supernova SN1987A, escaped the star due to their weak interactions, and then converted to gamma rays in the Galactic magnetic fields; the nonobservation of a gamma-ray flash coincident with the neutrino burst leads to strong constraints on the axion-photon coupling for axion masses m a ≲10 –10 eV. In this Letter, we use SN1987A to constrain higher mass axions, all the way to m a ~10 –3 eV, by accounting for axion production from the Primakoff process, nucleon bremsstrahlung, and pion conversion along with axion-photon conversion on the still-intact magnetic fields of the progenitor star. Moreover, we show that gamma-ray observations of the next Galactic supernova, leveraging the magnetic fields of the progenitor star, could detect quantum chromodynamics axions for masses above roughly 50 μ⁢eV, depending on the supernova. We propose a new full-sky gamma-ray satellite constellation that we call the GALactic AXion Instrument for Supernova (GALAXIS) to search for such future signals along with related signals from extragalactic neutron star mergers.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Biermann-Battery-Driven Magnetized Collisionless Shock Precursors in Laser-Produced Plasmas

This Letter reports the first complete observation of magnetized collisionless shock precursors formed through the compression of Biermann-battery magnetic fields in laser produced plasmas. At OMEGA, lasers produce a supersonic CH plasma flow which is magnetized with Biermann-battery magnetic fields. The plasma flow collides with an unmagnetized hydrogen gas jet plasma to create a magnetized shock precursor. The situation where the flowing plasma carries the magnetic field is similar to the Venusian bow shock. Imaging 2⁢𝜔 Thomson scattering confirms that the interaction is collisionless and shows density and temperature jumps. Proton radiographs have regions of strong deflections and FLASH magnetohydrodynamic (MHD) simulations show the presence of Biermann fields in the Thomson scattering region. Electrons are accelerated to energies of up to 100 keV in a power-law spectrum. Here, OSIRIS particle-in-cell (PIC) simulations, initialized with measured parameters, show the formation of a magnetized shock precursor and corroborate the experimental observables.

High-energy-density plasmas↗
DOE OSTI2020

Two-stage reflective self-seeding scheme for high-repetition-rate X-ray free-electron lasers

X-ray free-electron lasers (XFELs) open a new era of X-ray based research by generating extremely intense X-ray flashes. To further improve the spectrum brightness, a self-seeding FEL scheme has been developed and demonstrated experimentally. As the next step, new-generation FELs with high repetition rates are being designed, built and commissioned around the world. A high repetition rate would significantly speed up the scientific research; however, alongside this improvement comes new challenges surrounding thermal management of the self-seeding monochromator. Here, a new configuration for self-seeding FELs is proposed, operated under a high repetition rate which can strongly suppress the thermal effects on the monochromator and provides a narrow-bandwidth FEL pulse. Three-dimension time-dependent simulations have been performed to demonstrate this idea. Finally, with this proposed configuration, high-repetition-rate XFEL facilities are able to generate narrow-bandwidth X-ray pulses without obvious thermal concern on the monochromators.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2021

Laser heating system at the Extreme Conditions Beamline, P02.2, $\mathrm{PETRA III}$

A laser heating system for samples confined in diamond anvil cells paired with in situ X-ray diffraction measurements at the Extreme Conditions Beamline of PETRA III is presented. The system features two independent laser configurations (on-axis and off-axis of the X-ray path) allowing for a broad range of experiments using different designs of diamond anvil cells. The power of the continuous laser source can be modulated for use in various pulsed laser heating or flash heating applications. An example of such an application is illustrated here on the melting curve of iron at megabar pressures. The optical path of the spectroradiometry measurements is simulated with ray-tracing methods in order to assess the level of present aberrations in the system and the results are compared with other systems, that are using simpler lens optics. Based on the ray-tracing the choice of the first achromatic lens and other aspects for accurate temperature measurements are evaluated.

47 OTHER INSTRUMENTATION↗
DOE OSTI2023

Process-dependent anisotropic thermal conductivity of laser powder bed fusion AlSi10Mg: impact of microstructure and aluminum-silicon interfaces

Purpose AlSi10Mg alloy is commonly used in laser powder bed fusion due to its printability, relatively high thermal conductivity, low density and good mechanical properties. However, the thermal conductivity of as-built materials as a function of processing (energy density, laser power, laser scanning speed, support structure) and build orientation, are not well explored in the literature. This study aims to elucidate the relationship between processing, microstructure, and thermal conductivity. Design/methodology/approach The thermal conductivity of laser powder bed fusion (L-PBF) AlSi10Mg samples are investigated by the flash diffusivity and frequency domain thermoreflectance (FDTR) techniques. Thermal conductivities are linked to the microstructure of L-PBF AlSi10Mg, which changes with processing conditions. The through-plane exceeded the in-plane thermal conductivity for all energy densities. A co-located thermal conductivity map by frequency domain thermoreflectance (FDTR) and crystallographic grain orientation map by electron backscattered diffraction (EBSD) was used to investigate the effect of microstructure on thermal conductivity. Findings The highest through-plane thermal conductivity (136 ± 2 W/m-K) was achieved at 59 J/mm 3 and exceeded the values reported previously. The in-plane thermal conductivity peaked at 117 ± 2 W/m-K at 50 J/mm 3 . The trend of thermal conductivity reducing with energy density at similar porosity was primarily due to the reduced grain size producing more Al-Si interfaces that pose thermal resistance. At these interfaces, thermal energy must convert from electrons in the aluminum to phonons in the silicon. The co-located thermal conductivity and crystallographic grain orientation maps confirmed that larger colonies of columnar grains have higher thermal conductivity compared to smaller columnar grains. Practical implications The thermal properties of AlSi10Mg are crucial to heat transfer applications including additively manufactured heatsinks, cold plates, vapor chambers, heat pipes, enclosures and heat exchangers. Additionally, thermal-based nondestructive testing methods require these properties for applications such as defect detection and simulation of L-PBF processes. Industrial standards for L-PBF processes and components can use the data for thermal applications. Originality/value To the best of the authors’ knowledge, this paper is the first to make coupled thermal conductivity maps that were matched to microstructure for L-PBF AlSi10Mg aluminum alloy. This was achieved by a unique in-house thermal conductivity mapping setup and relating the data to local SEM EBSD maps. This provides the first conclusive proof that larger grain sizes can achieve higher thermal conductivity for this processing method and material system. This study also shows that control of the solidification can result in higher thermal conductivity. It was also the first to find that the build substrate (with or without support) has a large effect on thermal conductivity.

Engineering↗
DOE OSTI2020

BBOS: Efficient HPC Storage Management via Burst Buffer Over-Subscription

To avoid access to PFS, dedicated BB allocation is preferred despite of severe BB underutilization. Recently, new all-flash HPC storage systems with integrated BB and PFS are proposed, which speed up access to PFS. For this reason, we adopt BB over-subscription allocation method by allowing HPC applications to use BB only for I/O phase for improving BB utilization. Unfortunately, BB over-subscription aggravates I/O interference and demotion overhead from BB to PFS, resulting in degraded performance. To minimize the performance degradation, we develop an I/O scheduler to prevent I/O congestion and a new transparent data management system based on checkpoint/restart characteristics of HPC applications. With the proposed approach, not only the BB utilization can be improved, but also high performance of applications is achieved. In our experiments, we find that BB utilization is improved at least 2.2x, and more stable and higher checkpoint performance is guaranteed compared to other approaches. Besides, we achieve up to 96.4% hit ratio of restart requests on BB and up to 3.1x higher restart performance than others.

Sung, H↗