Search NASA⌕ Search

SEARCH · Search NASA

Results for “front propagation”

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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Grain boundary widening controls siderite (FeCO3) replacement of limestone (CaCO3)

Abstract The microstructure of minerals and rocks can significantly alter reaction rates. This study focuses on identifying transport paths in low porosity rocks based on the hypothesis that grain boundary widening accelerates reactions in which one mineral is replaced by another (replacement reaction). We conducted a time series of replacement experiments of three limestones (CaCO 3 ) of different microstructures and solid impurity contents using FeCl 2 . Reacted solids were analyzed using chemical imaging, small angle X-ray and neutron scattering and Raman spectroscopy. In high porosity limestones replacement is reaction controlled and complete replacement was observed within 2 days. In low porosity limestones that contain 1–2% dolomite impurities and are dominated by grain boundaries, a reaction rim was observed whose width did not change with reaction time. Siderite (FeCO 3 ) nucleation was observed in all parts of the rock cores indicating the percolation of the solution throughout the complete core. Dolomite impurities were identified to act as nucleation sites leading to growth of crystals that exert force on the CaCO 3 grains. Widening of grain boundaries beyond what is expected based on dissolution and thermal grain expansion was observed in the low porosity marble containing dolomite impurities. This leads to a self-perpetuating cycle of grain boundary widening and reaction acceleration instead of reaction front propagation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Nonlocal effects on thermal transport in hydrodynamic simulations of unmagnetized MagLIF-relevant gaspipes on NIF

We present simulations of heat flow relevant to gaspipe experiments on the National Ignition Facility to investigate kinetic effects on transport phenomena. D 2 and neopentane (C 5 H 12 ) filled targets are used to study the laser preheat stage of a MagLIF scheme where an axial magnetic field is sometimes applied to the target. Simulations were done with the radiation-MHD code HYDRA with a collision-dominated fluid model and the SNB nonlocal electron thermal conduction model. Using the SNB model to evolve the electron temperature increased the heat front propagation of neopentane gas targets compared to a local model by limiting radial heat flow. This increases electron temperature near the axis, which decreases laser absorption. We find that the effect of heat flow models on temperature profiles and laser propagation is modest. Beyond the SNB model, we utilize HYDRA to initialize plasma conditions for the Vlasov–Fokker–Planck K2 code. We run K2 until a quasi-steady state is reached and examine the impact of kinetic effects on heat transport. Although axial heat flow is well predicted by fluid models, the fluid model consistently overpredicts radial heat flow up to 150% in regions with the largest temperature gradient of D 2 filled gaspipes. On the other hand, the SNB nonlocal electron conduction model is found to be adequate for capturing kinetic heat flow in gaspipes.

Lau, Ryan Y. [Univ. of Colorado, Boulder, CO (Unit↗
DOE OSTI2021

Physics of turbulence spreading and explicit nonlocality

In this paper, we systematically derive a model for turbulence spreading from the basic kinetic equation. The model contains explicit nonlocal nonlinear diffusion and nonlocal growth. When the nonlocality scale parameter δ b (banana width) vanishes, this model reduces to the usual turbulence spreading model. We elucidate the mechanisms of nonlinear saturation and nonlocal growth. Results show that nonlocal effects, especially the nonlocal growth, thicken the turbulence spreading front and increase the speed of front propagation. More turbulence intensity penetrates the stable region when δ b increases. The penetration depth Δ p is proportional to $ \delta_b/L_T $, therefore the fraction of turbulence in the unstable region scales as 1 – δ b* . Lastly, the transport scales the same way.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Wire, hybrid, and laser-cut X-pinches as Talbot–Lau backlighters for electron density diagnostics

Talbot–Lau x-ray deflectometry (TXD) enables refraction-based imaging for high-energy-density physics experiments, and thus, it has been studied and developed with the goal of diagnosing plasmas relevant to inertial confinement and magnetic liner inertial fusion. X-pinches, known for reliably generating fast (~1 ns), small (~1 µm) x-ray sources, were driven on the compact current driver generator for ablation structure and implosion studies (~200 kA, 150 ns) as a potential backlighter source for TXD. Considering that different X-pinch configurations have characteristic advantages and drawbacks as x-ray generating loads, three distinct copper X-pinch configurations were studied: the wire X-pinch, the hybrid X-pinch, and the laser-cut X-pinch. Furthermore, the Cu K-shell emission from each configuration was characterized and analyzed regarding the specific backlighter requirements for an 8 keV TXD system: spatial and temporal resolution, number of sources, time of emission, spectrum, and reproducibility. Recommendations for future experimental improvements and applications are presented. The electron density of static objects was retrieved from Moiré images obtained through TXD. This allowed to calculate the mass density of static samples within 4% of the expected value for laser-cut X-pinches, which were found to be the optimal X-pinch configuration for TXD due to their high reproducibility, small source size (≤5 µm), short duration (~1 ns), and up to 106 W peak power near 8 keV photon energy. Plasma loads were imaged through TXD for the first-time using laser-cut X-pinch backlighting. Experimental images were compared with simulations from the x-ray wave-front propagation code, demonstrating that TXD can be a powerful x-ray refraction-based diagnostic for dense Z-pinch loads. Future plans for Talbot–Lau interferometry diagnostics in the pulsed-power environment are described.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Staircase resiliency in a fluctuating cellular array

Inhomogeneous mixing by stationary convective cells set in a fixed array is a particularly simple route to layering. Layered profile structures, or staircases, have been observed in many systems, including drift-wave turbulence in magnetic confinement devices. The simplest type of staircase occurs in passive-scalar advection, due to the existence and interplay of two disparate timescales, the cell turn-over (τ H ), and the cell diffusion (τ D ) time. In this simple system, we study the resiliency of the staircase structure in the presence of global transverse shear and weak vortex scattering. The fixed cellular array is then generalized to a fluctuating vortex array in a series of numerical experiments. The focus is on regimes of low-modest effective Reynolds numbers, as found in magnetic fusion devices. By systematically perturbing the elements of the vortex array, we learn that staircases form and are resilient (although steps become less regular, due to cell mergers) over a broad range of Reynolds numbers. The criteria for resiliency are (a) τ D >>τ H and (b) a sufficiently high profile curvature (κ ≥ 1.5). We learn that scalar concentration travels along regions of shear, thus staircase barriers form first, and scalar concentration "homogenizes" in vortices later. The scattering of vortices induces a lower effective speed of scalar concentration front propagation. The paths are those of the least time. We observe that if background diffusion is kept fixed, the cell geometric properties can be used to derive an approximation for the effective diffusivity of the scalar. Furthermore, the effective diffusivity of the fluctuating vortex array does not deviate significantly from that of the fixed cellular array.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Solid structure of Li 2 BeF 4 (FLiBe) from room temperature to melting studied by neutron and X-ray diffraction

Molten fluoride salts such as Li 2 BeF 4 (FLiBe) are used in molten salt reactors, fluoride-salt-cooled high-temperature reactors and fusion reactors as a fuel solvent, coolant and/or tritium breeding medium. In engineered systems that use molten salt, solid-state material will be present during melting and freezing scenarios, and therefore the temperature-dependent properties of the solid and solid/liquid phase transition merit investigation. To observe the behavior of the solid state of Li 2 BeF 4 from room temperature to melting, this work used neutron and X-ray diffraction to measure the changes in the lattice parameters and volume of the crystalline unit cell and compared the results with prior low-temperature data for solid Li 2 BeF 4 . From neutron diffraction data it is also possible to identify anisotropy: centimetre-scaled crystals align preferentially with the a axes parallel to the direction of freezing front propagation, and the c axes expand 54% more than the a axes. This work provides the lattice constants as a function of temperature, quantifies the thermal expansion, and determines the equation describing the change in density for solid Li 2 BeF 4 from room temperature to 459°C to be ρ solid (kg m −3 ) = 2182 (3) − 0.115 (2) T (°C) and the volume expansion upon melting to be less than 5%. This density changes depending on molecular weight and enrichment.

36 MATERIALS SCIENCE↗
DOE OSTI2021

A Simple Lagrangian Parcel Model for the Initiation of Summer-time Mesoscale Convective Systems over the Central United States

Mesoscale convective systems (MCSs) account for more than 50% of summertime precipitation over the central United States and have a significant impact on local weather and hydrologic cycle. It is hypothesized that the inadequate treatment of MCSs is responsible for the long-standing warm and dry bias over the central United States in coarse-resolution general circulation model (GCM) simulations. In particular, a better understanding of MCS initiation is still lacking. Here a single-column Lagrangian parcel model is first developed to simulate the basic features of a rising parcel. This simple model demonstrates the collective effects of boundary layer moistening and dynamical lifting in triggering convective initiation and reproduces successfully its early afternoon peak with surface equivalent potential temperature as a controlling factor. It also predicts that convection is harder to trigger in the future climate under global warming, consistent with the results from convection-permitting regional climate simulations. Then, a multicolumn model that includes an array of single-column models aligned in the east–west direction and incorporates idealized cold pool interaction mechanisms is developed. The multicolumn model captures readily the cold pool–induced upscale growth feature in MCS genesis from initially scattered convection that is organized into a mesoscale cluster in a few hours. It also highlights the crucial role of lifting effects due to cold pool collision and spreading, subsidence effect, and gust front propagation speed in controlling the final size of mesoscale clusters and cold pool regions. This simple model should be useful for understanding fundamental mechanisms of MCS initiation and providing guidance for improving MCS simulations in GCMs.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Cooling flow regime of a plasma thermal quench

A large class of Laboratory, Space, and Astrophysical plasmas is nearly collisionless. When a localized energy or particle sink, for example, in the form of a radiative cooling spot or a black hole, is introduced into such a plasma, it can trigger a plasma thermal collapse, also known as a thermal quench in tokamak fusion. Here we show that the electron thermal conduction in such a nearly collisionless plasma follows the convective energy transport scaling in itself or in its spatial gradient, due to the constraint of ambipolar transport. As a result, a robust cooling flow aggregates mass toward the cooling spot and the thermal collapse of the surrounding plasma takes the form of four propagating fronts that originate from the radiative cooling spot, along the magnetic field line in a magnetized plasma. The slowest one, which is responsible for deep cooling, is a shock front.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Simulation of emergent compaction banding fronts caused by frictional boundaries

This study examines the role of boundary friction in promoting heterogeneous compaction in soft rock specimens loaded at high confining pressure outside the domain of compaction localisation. An elastoplastic constitutive model characterised by tunable hardening/softening behaviour is used to conduct the analyses. Finite-element simulations suggest that material instability is a non-necessary condition for the emergence of compaction fronts. Such fronts propagated as a result of a severe deviation in the local responses induced by frictional constraints. Overall, these findings suggest that boundary effects can bias the assessment of the extent of the compaction localisation domain. Experimental countermeasures and informed model calibration procedures are therefore necessary to minimise such bias and enable more accurate predictions of soft rock compaction.

42 ENGINEERING↗
DOE OSTI2023

Rapid Simulation of Unconventional Reservoirs by Multidomain Multiresolution Modeling Based on the Diffusive Time of Flight

Summary The fast marching method (FMM)-based rapid flow simulation has been shown to accelerate simulation efficiency by orders of magnitude by transforming 3D simulation to equivalent 1D simulation using the concept of the “diffusive time of flight” (DTOF). However, the 1D transformation does not directly apply to multiwell problems. In this paper, we propose a novel DTOF-based multidomain multiresolution discretization scheme to accelerate multiwell simulation of unconventional reservoirs. Our method formulates multiwell simulation problems based on the DTOF which displays the pressure front propagation in unconventional reservoirs. The DTOF contours are used to partition the reservoir into local and shared domains. A local domain is where the flow is dominated by a single well, and the shared domain is where the fluid flow is influenced by multiple wells. The DTOF contours expand independently in local domains and interfere in the shared domain. After the partitioning, each domain is discretized using a multiresolution scheme whereby the original 3D fine mesh is preserved near the wells to account for detailed physics including gravity, and the rest of the domain is discretized into 1D mesh based on the DTOF contours to alleviate the simulation workload. The power and efficacy of our approach are demonstrated using synthetic and field-scale simulation models with different degrees of geologic and well-completion complexity. The simulation results, number of active cells, and computation time for the proposed discretization scheme are compared with the original high-fidelity 3D model for each case. The results show that the proposed method is suitable for multiwell simulation problems in unconventional reservoirs and can accelerate flow simulations by orders of magnitude with minimal loss of accuracy. The novelty of this work is the creation of DTOF-derived multiresolution discretization with local and shared domains to simplify and accelerate the calculation of subsurface flow problems, especially in unconventional reservoirs. Our workflow can be easily interfaced with commercial simulators, making it suitable for large-scale field applications.

Engineering↗
DOE OSTI2022

Thermal quench of open field plasma intercepting with recycling walls

When a fusion plasma suddenly intercepts a solid surface (wall or pellet), thermal collapse is distinctly kinetic & has novel physics. VPIC simulations & theory revealed the fascinating dynamics of four propagating fronts controlling parallel electron temperature cooling.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Microscopy X-ray imaging enriched with small angle X-ray scattering for few nanometer resolution reveals shock waves and compression in intense short pulse laser irradiation of solids

Understanding how laser pulses compress solids into high-energy-density states requires diagnostics that simultaneously resolve macroscopic geometry and nanometer-scale structure. Here we present a combined X-ray imaging (XRM) and small-angle X-ray scattering (SAXS) approach that bridges this diagnostic gap. Using the Matter in Extreme Conditions end station at LCLS, we irradiated 25 μm copper wires with 45 fs, 0.9 J, 800 nm pulses at 3.5 × 10 19 W/cm 2 while probing with 8.2 keV XFEL pulses. XRM visualizes the evolution of ablation, compression, and inward-propagating fronts with ∼ 200 nm resolution, while SAXS quantifies their nanometer-scale sharpness via the time-resolved evolution of scattering streaks. The joint analysis reveals that an initially smooth compression steepens into a nanometer-sharp shock front after t sh ≈ (18 ± 3) ps, consistent with an analytical steepening model and hydrodynamic simulations. The front reaches a velocity of c sh ≈ 25 k m / s and a lateral width of several tens of microns, demonstrating direct observation of shock formation and decay at solid density for the first time with few-nanometer precision. This integrated XRM–SAXS method establishes a quantitative, multi-scale diagnostic of laser-driven shocks in dense plasmas relevant to inertial confinement fusion, warm dense matter, and planetary physics.

Kluge, Thomas [Helmholtz-Zentrum Dresden-Rossendor↗
DOE OSTI2026

Stochastic Thermo-Hydro Modeling and Neural Network Surrogate Development for Thermal Resource Assessment of the Galleries-to-Calories Geobattery

The Galleries-to-Calories Geobattery concept explores the use of abandoned coal mine workings for large-scale thermal energy transport and storage. The system involves injecting waste heat from a supercomputing facility into flooded mine galleries, where groundwater flow can store and transport thermal energy for potential recovery in downgradient district heating and cooling applications. To evaluate the feasibility and performance of the Geobattery under geological and operational uncertainty, we developed a suite of stochastic thermo-hydrological (TH) simulations using Monte Carlo sampling of key uncertain parameters (e.g., permeability, porosity, thermal conductivity, specific heat capacity) and operating conditions (e.g., injection rate, injection temperature). Results identified injection rate and temperature as the most influential parameters governing thermal front propagation, while the geometry of the room-and-pillar structure played a critical role in directing the extent and orientation of thermal advancement. Optimal combinations of material properties for maximizing heat recovery were also determined. To address the high computational cost of coupled-process stochastic modeling, we trained a neural network surrogate model on 24,000 physics-based realizations, achieving an R² > 0.99 and MAE < 0.1 for temperature predictions at monitoring locations. This surrogate enabled an additional 100,000 realizations for global sensitivity analysis and probabilistic thermal resource assessment. The integrated stochastic physics–surrogate modeling framework offers a computationally efficient tool for quantifying uncertainty, identifying key drivers, and informing early-stage design decisions for Geobattery systems.

15 - GEOTHERMAL ENERGY↗
DOE OSTI2025

Morphogenic Growth 3D Printing

Inspired by nature's morphogenesis, a new 3D printing process –growth printing (GP)– takes advantage of a self‐propagating curing front to produce 3D polymeric parts following a growth‐like development plan. The propagation of the curing front is driven by the exothermic polymerization of dicyclopentadiene (DCPD), which transforms the liquid resin into a stiff polymer as it propagates at 1 mm s −1 . GP is triggered when a heated initiator contacts the uncured liquid resin in an open container. The initiator nucleates the frontal polymerization reaction and the isotropic radial propagation of the growth front. Simultaneously, the initiator is moved up across the free surface of the resin, pulling the cured object out of the uncured resin. The motion trajectory of the initiator with respect to the free resin surface controls the growth morphology of the 3D part. An inverse design algorithm is developed to produce 3D parts by modeling the reaction‐diffusion‐driven solidification process. This process has substantial energy savings and high printing speeds.

3D printing↗
DOE OSTI2020

A Reversible Structural Phase Transition by Electrochemically-Driven Ion Injection into a Conjugated Polymer

We nd that conjugated polymers can undergo reversible structural phase transitions during electrochemical oxidation and ion injection. We study poly[2,5-bis(thiophenyl)-1,4-bis(2-(2-(2- methoxyethoxy)ethoxy)ethoxy)benzene] (PB2T-TEG), a conjugated polymer with glycolated side chains. Using grazing incidence wide angle X-ray scattering (GIWAXS), we show that, in contrast to previously known polymers, this polymer switches between two structurally distinct crystalline phases associated with electrochemical oxidation/reduction in an aqueous electrolyte. Importantly, we show that this unique phase change behavior has important physical consequences for ion transport. Notably, using moving front experi- ments visualized by both optical microscopy and super-resolution photoinduced force microscopy (PiFM), we show that a propagating ion front in PB2T-TEG exhibits non-Fickian transport, retaining a sharp step-edge prole, in stark contrast to the Fickian diusion more commonly observed. This structural phase transition is reminiscent of those accompanying ion uptake in inorganic materials like LiFePO4. We propose that engineering similar properties in future conjugated polymers may enable the realization of new materials with superior performance in electrochemical energy storage or neuromorphic memory applications.

Bischak, Connor G.↗
DOE OSTI2022

Hydrodynamic simulations of electron-capture supernovae: progenitor and dimension dependence

ABSTRACT We present neutrino-transport hydrodynamic simulations of electron-capture supernovae (ECSNe) in flash with new two-dimensional (2D) collapsing progenitor models. These progenitor models feature the 2D modelling of oxygen-flame propagation until the onset of core collapse. We perform axisymmetric simulations with six progenitor models that, at the time of collapse, span a range of propagating flame front radii. For comparison, we also perform a simulation with the same set-up using the canonical, spherically symmetrical progenitor model n8.8. We found that the variations in the progenitor models inherited from simulations of stellar evolution and flame propagation do not significantly alter the global properties of the neutrino-driven ECSN explosion, such as the explosion energy (∼1.36–1.48 × 1050 erg) and the mass (∼0.017–0.018 M⊙) and composition of the ejecta. Due to aspherical perturbations induced by the 2D flame, the ejecta contains a small amount (≲1.8 × 10−3 M⊙) of low-Ye (0.35 < Ye < 0.4) component. The baryonic mass of the protoneutron star is ∼1.34 M⊙ (∼1.357 M⊙) with the new (n8.8) progenitor models when simulations end at ∼400 ms and the discrepancy is due to updated weak-interaction rates in the progenitor evolutionary simulations. Our results reflect the nature of ECSN progenitors containing a strongly degenerate oxygen–neon–magnesium (ONeMg) core and suggest a standardized ECSN explosion initialized by ONeMg core collapse. Moreover, we carry out a rudimentary three-dimensional simulation and find that the explosion properties are fairly compatible with the 2D counterpart. Our paper facilitates a more thorough understanding of ECSN explosions following the ONeMg core collapse, though more three-dimensional simulations are still needed.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Advanced Finite-Volume Numerics and Source Term Assumptions for Kernel and G-Equation Modelling of Propane/Air Flames

Here G-Equation models represent propagating flame fronts with an implicit two-dimensional surface representation (level-set). Level-set methods are fast, as transport source terms for the implicit surface can be solved with finite-volume operators on the finite-volume domain, without having to build the actual surface. However, they include approximations whose practical effects are not properly understood. In this study, we improved the numerics of the FRESCO CFD code’s G-Equation solver and developed a new method to simulate kernel growth using signed distance functions and the analytical sphere-mesh overlap. We analyzed their role for simulating propane/air flames, using three well-established constant-volume configurations: a one-dimensional, freely propagating laminar flame; a disc-shaped, constant-volume swirl combustor; and torch-jet flame development through an orifice from a two-chamber device. We tested the explicit (sub-cycled) vs. implicit formulation for the standard transport operators (advection, diffusion, compressibility). In addition to the accurate flame swept-volume method for chemistry and species source term, we developed a more accurate estimator for the burnt/unburnt split cell composition. Then, we developed a signed-distance-function (SDF) based method which provides a more stable reinitialization of the level-set field at every time-step. We found that simplifying assumptions common to several G-Equation implementations, for straightforward terms such as compressibility and advection, lead to large errors in predicting the propagation of even laminar flames, with deviations up to ~300% in simulated vs. formulated flame speed. Conversely, the enhanced numerics enabled through the SDF field reinitialization and improved chemistry source term improve simulation stability and smooth flame propagation even with significantly larger solver time-steps.

42 ENGINEERING↗