Supersonic combustion in premixed hydrogen-air flows.
Ignition and propagation of turbulent flame front in jet flow investigated for combustible hydrogen- air mixture
FLAME PROPAGATION↗
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.
Air-coupled tsunamis generated from impacts and airbursts: Our understanding before Hunga-Tonga Hunga-Ha'apai
The effort to prevent or mitigate the effects of an impact on Earth is known as planetary defense. A significant component of planetary defense research involves risk assessment. Much of our understanding of the risk from near-Earth objects comes from the geologic record in the form of impact craters, but not all asteroid impacts are crater-forming events. Small asteroids explode before reaching the surface, generating an airburst, and most impacts into the ocean do not penetrate the water to form a crater in the sea floor. The risk from these non-crater-forming ocean impacts and airbursts is difficult to quantify and represents a significant uncertainty in our assessment of the overall threat. We are currently working to better understand impact scenarios that can generate dangerous tsunamis. One of the suggested mechanisms for the production of asteroid–generated tsunamis is by direct coupling of the pressure wave to the water, analogous to the means by which a moving weather front can generate a meteotsunami. To test this hypothesis, we ran a series of airburst simulations and provided time-resolved pressure and wind profiles to use as source functions for tsunami models. We used the CTH hydrocode to model the various airburst scenarios to compare to the results of other simulations and provide time dependent boundary conditions as input to shallow-water wave propagation codes. The strongest and most destructive meteotsunamis are generated by atmospheric pressure oscillations with amplitudes of only a few hPa1 (mbar), corresponding to changes in sea level of a few cm. The resulting wave is strongest when there is a resonance between the ocean and the atmospheric forcing. A Proudman resonance takes place when the atmospheric disturbance’s translational speed (U) equals the longwave phase speed $\sqrt{gh}$ of shallow water wave. Coupling is strongest when the Froude number (Fr=U/c) is unity. A weather front propagates much slower than the speed of sound, so meteotsunamis are most common and dangerous in shallow bodies of water such as the Mediterranean Sea or Lake Michigan. By contrast, the blast wave from an airburst or crater-forming impact propagates at a speed faster than a tsunami in the deepest ocean, and a Proudman resonance cannot be achieved even though the overpressures are orders of magnitude greater. However, blast wave profiles are N-waves in which a sharp shock wave leading to overpressure is followed by a more gradual rarefaction to a much longer-duration underpressure phase. Even though the blast outruns the water wave it is forcing, the tsunami should continue to be driven by the out-of-resonance gradient associated with the suction phase, which may depend strongly on the details of the airburst or impact scenario. The open question is whether there are any conditions under which such an airburst-driven tsunami can be dangerous enough to contribute to the overall impact risk. We have also identified other potential mechanisms for airburst-generated tsunamis: 1) reaction force at the surface from the plume ejected into space, which carries significant momentum, 2) expanding toroidal vortices at the surface, which travel more slowly than the shock wave and can generate a Proudman resonance in relatively shallow ocean (such as continental shelf), and 3) steam explosion from seawater ablation by a “Type II” (Libyan Desert Glass-type) airburst in which the hot vapor jet descends to the surface. On January 15, 2022, the Hunga-Tonga Hunga-Ha’apai volcano, located approximately 60 km north of Tongatapu, the main island of Tonga, violently erupted with a powerful explosion, culminating the period of volcanic activity that started in December of 2021. This event and resulting tsunamis provided an existence proof for the air pressure wave coupling mechanism we proposed. It also suggests that it can be stronger and more significant over much greater distances than we contemplated, leading to global tsunamis associated with impact events on land as well as in the water. Large atmospheric explosions generate global Lamb waves with larger amplitudes, longer periods, and slower speeds than the local and regional blast waves we modeled prior to that event. This paper reviews our analysis and modeling of airburst-driven tsunamis prior to the 2022 Hunga-Tonga Hunga-Ha’apai tsunami, which was the subject of two presentations at the 2023 Planetary Defense Conference and is the subject of another paper currently in preparation.
54 ENVIRONMENTAL SCIENCES↗
An experimental and analytical investigation of ignition and axisymmetric turbulent flame propagation in premixed hydrogen-air flows with application to supersonic combustion.
Ignition and propagation of turbulent flame front in jet flow investigated for combustible hydrogen- air mixture
GAS MIXTURE↗
The delayed-detonation model of a type Ia supernovae. 1: The deflagration phase
The nature of the 'delayed detonation' mechanism of Khokhlov for the explosion of Type Ia supernovae is investigated by using two-dimensional numerical hydrodynamics simulations. A new algorithm is used to treat the deflagration front. Assuming that it propagates locally at the laminar flame speed, the deflagration is insufficient to unbind the star. Expansion shuts of the flame; much of this small production of iron group nuclei occurs at lower densities, which reduces the electron-capture problem. The burning front does become wrinkled, but the wavelength of the instability is much larger than the computational grid size and is resolved; this is consistent with previous analysis. Because the degenerate star has an adiabatic exponent only slightly above 4/3, the energy released by deflagration drives a pulsation of large amplitude. During the first expansion phase, adiabatic cooling shuts off the burning, and a Rayleigh-Taylor instability then gives mixing of high-entropy ashes with low-entropy fuel. During the first contraction phase, compressional heating reignites the material. This paper deals with the deflagration phase, from the onset of burning, through expansion and quenching of the flame, to the first contraction.
Arnett, David↗
“Tailing” Phenomenon During Polymer Propagation at The Milne Point Polymer Flood
Chemical flooding simulators traditionally assume that polymer retention follows the Langmuir isotherm, while fractional-flow calculations assume concentration-independent retention. Both models predict a delay in propagation of the polymer front (in proportion to the retention value), followed by a rapid rise in produced polymer concentration to the injected level. In contrast, laboratory retention studies for the Milne Point polymer flood (North Slope of Alaska) consistently show virtually no delay in polymer propagation, a rapid rise in produced polymer to 70-90% of the injected concentration, followed by produced concentration gradually approaching the injected value over many pore volumes. From a practical viewpoint, this behavior means that retention causes no significant delay in propagation of the polymer bank (and therefore the oil bank), but the effective viscosity and displacement efficiency are less than originally planned. This paper explores why this polymer “tailing” phenomenon occurs. During core floods, effluent was analyzed for brine-tracer concentration, viscosity, total organic carbon, and total nitrogen. The latter three items allowed three independent measures of polymer concentration and polymer retention. Experiments were performed with Milne Point core material, oil, brine, and conditions. Mechanistic experiments were also performed with packs of glass beads with controlled levels of various minerals. In Milne Point cores, the presence of the tailing phenomenon was not sensitive to flow rate, polymer concentration, core heterogeneity, or whether the core was preserved, cleaned of oil, or cleaned and resaturated with oil, or cleaned, re-saturated and aged with oil. The three independent measures of polymer concentration confirmed that the tailing phenomenon was not an artifact associated with polymer degradation or detection interferences. Tailing was also observed during mechanistic floods using glassbead packs when sufficient levels of illite or kaolinite were present. This observation was consistent with high levels of illite noted in Milne Point cores. The tailing was not noted during mechanistic floods that contained beads with montmorillonite, chlorite, calcium carbonate, dolomite, siderite, pyrite, or calcium sulfate. A model is proposed to account for the tailing phenomenon. This work suggests that mineralogy analysis (especially for illite and kaolinite) may reveal whether tailing should be accounted for during simulations of polymer propagation/retention in a given field application.
Seright, Randy↗
Modeling Delamination in Postbuckled Composite Structures Under Static and Fatigue Loads
The ability of the Abaqus progressive Virtual Crack Closure Technique (VCCT) to model delamination in composite structures was investigated for static, postbuckling, and fatigue loads. Preliminary evaluations were performed using simple Double Cantilever Beam (DCB) and Mixed-Mode Bending (MMB) specimens. The nodal release sequences that describe the propagation of the delamination front were investigated. The effect of using a sudden or a gradual nodal release was evaluated by considering meshes aligned with the crack front as well as misaligned meshes. Fatigue simulations were then performed using the Direct Cyclic Fatigue (DCF) algorithm. It was found that in specimens such as the DCB, which are characterized by a nearly linear response and a pure fracture mode, the algorithm correctly predicts the Paris Law rate of propagation. However, the Abaqus DCF algorithm does not consider different fatigue propagation laws in different fracture modes. Finally, skin/stiffener debonding was studied in an aircraft fuselage subcomponent in which debonding occurs deep into post-buckling deformation. VCCT was shown to be a robust tool for estimating the onset propagation. However, difficulties were found with the ability of the current implementation of the Abaqus progressive VCCT to predict delamination propagation within structures subjected to postbuckling deformations or fatigue loads.
Bisagni, Chiara↗
Modal propagation angles in a cylindrical duct with flow and their relation to sound radiation
The angles of propagation for the wave front making up a duct mode are presented with the Mach number in the duct. Approximate equations are derived to provide simple utilitarian expressions. These expressions are valid only near the outer wall which is the most important region since the bulk of the acoustic intensity is located there. Exact solutions using Hankel functions are given in an appendix. These data corroborate the approximate solution accuracy near the outer wall. The axial propagation angle is used to infer information about the far-field radiation pattern. The resultant axial angle of propagation in the duct is shown to agree exactly with the peak of the principal lobe to far-field radiation obtained from formal radiation calculations when the Mach number is uniform everywhere. The obtained solution is extended to cover the case of different Mach numbers inside and outside the duct for which exact calculations have not been available for engine inlet configurations.
Rice, E. J.↗
Negative streamer development in FEP teflon
A computational model is developed which describes the evolution and propagation of an ionizing front (negative streamer) in solid materials. The ionization front consists of drifting avalanching electrons moving self-consistently under the influence of their own space-charge field together with an applied external field. The required input information for the model consists of the functional dependence of the macroscopic transport coefficients on the local electric field, the initial conditions for beginning the calculation, and the strength of the applied field. A computational approach for specifying the transport coefficients and initional conditions is also described. The approach has been implemented by constructing three computer codes which sequentially interface, beginning with single electron scattering, and ending with streamer development. Computational results are presented for model calculations in Teflon. The overall model is perceived to provide a picture of the initiation phase of a propagating discharge in electron-irradiated dielectrics.
Beers, B. L.↗
Possible observation of a disconnected magnetic structure in a coronal transient
Solar Maximum Mission coronagraph/polarimeter observations of the coronal transient of March 15-16, 1980, which show a bright front convex toward the sun propagating out through the corona at 175 km/sec, are reported. A fan-shaped bright region connecting this front to the inner corona collapses into a narrow, bright ray over the next few hours. These observations are interpreted as resulting from the disconnection of magnetic loops from the underlying prominencelike loop structures seen rising earlier in the event. Other similar inverted arches exist in the C/P data, but are usually seen with poorer temporal resolution.
Illing, R. M. E.↗
Characterization of mechanical discontinuities based on data-driven classification of compressional-wave travel times
Wave propagation and diffusive transport phenomena are influenced by the mechanical discontinuities in material. Here this study shows that certain bulk properties of the network of low-velocity mechanical discontinuities (e.g. air-filled cracks) in a material can be characterized by processing compressional-wave travel times using traditional data-driven classification techniques. To that end, we perform three tasks in chronological order: (1) use the discrete fracture network (DFN) method to create two-dimensional (2D) numerical models of crack-bearing material embedded with various types of low-velocity mechanical discontinuities, (2) use the fast marching method (FMM) to simulate the propagation of the wave/diffusion front from a single source through the 2D crack-bearing material to multiple receivers placed on the boundary of the material, and (3) train 9 data-driven classifiers to characterize the crack-bearing materials (i.e. bulk properties of the network of mechanical discontinuities in the crack-bearing material) by learning from the simulations of travel times detected by multiple receivers placed around the crack-bearing material. The classifiers identified the orientation, spatial distribution, and dispersion of the low-velocity mechanical discontinuities. Voting classifier performs the best among the 9 classifiers. For the characterization of bulk dispersion and distribution of discontinuities, the sensors located on the adjacent boundaries are more important; whereas for the characterization of bulk orientation of discontinuities, the sensors located on the opposite side are more important.
42 ENGINEERING↗
In-situ re-crystallization of heavily-irradiated Gd 2 Ti 2 O 7
Gadolinium titanate (Gd 2 Ti 2 O 7 , or GTO) and other lanthanide pyrochlores are interesting candidates for actinide waste disposal and fast ion conduction because the relevant material properties are intimately dependent on local cation structure. Therefore, a deep understanding of the kinetics associated with cation ordering and disordering is required if such material properties are to be tuned for specific device designs. To this end, single crystals of GTO were irradiated with 190 keV helium ions to a total fluence of 1 × 10 17 ions/cm 2 , amorphizing the sample surface to a depth of ~1 μm and resulting in significant He bubble accumulation. FIB lamellae lifted out from the irradiated sections were examined during heat treatment in the (scanning) transmission electron microscope. Two distinct stages of the re-crystallization of the amorphized material were observed. The material near the end of the ions’ range transformed first and with the same orientation as the pristine material. This was due to the close proximity of the pristine material and the presence of small defect fluorite seeds, but the propagation of this growth front was frustrated by the large pores (formerly bubbles) in the He accumulation layer. This was followed by heterogeneous nucleation of new crystallites at random orientations at the top of the He accumulation layer, which is attributed to the high surface area associated with the many small He bubbles in that region. Finally, it is inferred that the kinetics of grain growth in this material are significantly faster than the kinetics of grain nucleation.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Enhancing photoionization rate calculations in low-temperature plasmas using spectral methods
Photoionization plays a central role in the development of streamer discharges and other non-equilibrium plasma phenomena. It creates seed electrons, which are essential for positive streamer propagation, allowing the ionization front to move forward. Because of this, accurate modeling of photoionization is very important for predicting streamer behavior and plasma evolution. The photoionization process in air (N 2 – O 2 mixture) is often described by the Zheleznyak model (1982). This model is usually solved through Helmholtz-type equations that approximate the Zheleznyak photoionization model (Zheleznyak et al. 1982) as Partial Differential Equations (PDEs). Conventional numerical methods, such as the Finite Difference Method (FDM) or Finite Volume Method (FVM), are widely used to solve these equations. Although they are prevalent, the computational cost of these methods due to their need for matrix operations and iterative solver is demanding. To address this challenge, this work develops a spectral solver based on the Fast Fourier Transform (FFT) combined with Discrete Cosine Transform (DCT) and Discrete Sine Transform (DST) to calculate the photoionization rate efficiently in an axisymmetric cylindrical domain. This method naturally satisfies the boundary conditions used in the model and converts the PDE into algebraic ones in spectral space. Thus, avoids the need for iterative matrix solvers. When compared with FDM results, it is demonstrated that the new solver not only maintains accuracy, but also reduces the computational cost, showing a performance increase of approximately 100 compared to FDM over a wide range of problem sizes. The method is parallelized using Message Passing Interface (MPI) and has been integrated into a fluid plasma model for streamer simulation. Here, this FFT-based approach provides a fast and reliable alternative for calculating photoionization in fluid models, helping large-scale plasma simulations run faster and efficiently, and allows higher-resolution simulation without extra computational cost.
Axisymmetric system↗
A fast reduced model for a shell-and-tube based latent heat thermal energy storage heat exchanger and its application for cost optimal design by nonlinear programming
Numerical simulation of latent heat thermal energy storage (LHTES) systems plays a fundamental role in studying the physical process and guiding the engineering design. Discretization of the PDEs describing the nonlinear solidification/melting process of phase change materials (PCMs) leads to a large-scale complex dynamics system, where the system behavior depends on a set of parameters. In a design setting, repeated model evaluations are required over the set of parameters results in significant computational burden. In this paper, an explicit analytic solution was built for the propagation of the solidification front in a cylindrical coordinate. The analytic solution approach is further employed to develop a low computational reduced model (RM) as a module for a shell-and-tube based LHTES heat exchanger. The levelized Cost of Energy (LCOE) is used as a design metric and the RM model is used to apply system-level constraints in the nonlinear programming formulation that facilitates efficient global optimal design of the PCM properties, flow conditions and tube geometries. The use of LCOE as the design metric prevents over design of the heat transfer rate and also establishes a fair ground for evaluation of different thermal storage technologies and their integrated applications with other systems. Optimal results showed that a higher effectiveness results in a higher LCOE; the velocity of the HTF and the length of the channel are highly correlated with each other; both larger PCM latent energy and conductivity result in lower LCOE.
25 ENERGY STORAGE↗
High-Fidelity Simulation of the Light-to-Dense Stratification Transient in the HiRJET Facility
Density stratification in a large enclosure is a crucial phenomenon to heat transfer and sustainable passive heat removal of a sodium fast reactor during reactivity transients. However, engineering turbulence models were identified to have unsatisfactory performance in predicting propagation of a stratified front. Yet, the scarcity of high-resolution data for stratification hampers the development of models. To explor e applications of leveraging direct numerical simulation (DNS) data to support turbulence model development, this work conducted DNS using NekRS to study a long stratification transient in the High-Resolution Jet (HiRJET) experimental facility. Here, this work considers an experiment run where light fluid is injected into a tank containing a denser fluid with a relative density difference of 1.5%. Formation of the stratified layer is identified as impingement of the buoyant jet promoting mixing of the two fluids. Based on the transient statistics, transport of the concentration can be characterized by regions with dominating effects of turbulent mixing, buoyant dissipation, and molecular diffusion, respectively, as moving away from the elevation of jet impingement. Concentration near the stratified front also exhibits oscillation at Brunt-Väisälä frequency. Preliminary validation of the simulation showed encouraging agreement of the concentration distribution with the reference experiment.
buoyant mixing↗
Unified wavefront singularity characterization of three-dimensional elastodynamic time-domain half-space Green's function under impulsive boundary and internal loads
Founded on a novel analytical formulation that led to a rigorous yet compact path-integral representation of the time-domain elastodynamic half-space Green's function, a unified analysis of the possible occurrence of different singular wavefront behaviour in the response under arbitrary impulsive internal or surface point loads at arbitrary source-receiver locations is presented. With the decomposition of the general solution into distinct initiating and reflected wave group integrals that share a common factored format and simple contour definitions, the mathematical framework is shown to allow a straightforward identification of the specific conditions and the particular wave groups that are responsible for the singular wavefront phenomena without resorting to advanced analytic function theories or asymptotic methods. Analytic characterizations of the nature, strength and direction of all intrinsic singular wavefront behaviours of the three-dimensional Green's function in three canonical cases of source-receiver configurations are given in a dual integral-closed form format to facilitate their theoretical understanding as well as computational applications. Graphical illustrations of their variation with the source-receiver configuration and the medium's Poisson's ratio together with relevant comparison and clarifications of some classical treatments are included.
Green’s functions↗
Single-shot spatiotemporal plasma density measurements with a chirped probe pulse
In this work, we present the development and demonstration of a diagnostic for the measurement of the spatial and temporal evolution of plasma density in a single shot. Single-shot Advanced Plasma Probe HolographIc REconstruction (SAPPHIRE) utilizes a chirped probe pulse, a diffractive optical element, a self-referenced interferometer, and an interference bandpass filter to achieve high-fidelity electron density measurements suitable for underdense plasmas that exhibit cylindrical symmetry. The method overcomes limitations in conventional diagnostics, such as reliance on shot-to-shot reproducibility, while capturing plasma dynamics on picosecond timescales with micron-level spatial resolution. The capabilities of SAPPHIRE are demonstrated through measurements of laser-driven plasma channels in helium–nitrogen gas jets. SAPPHIRE demonstrates the formation and expansion of plasma channels in a single shot and the propagation of supersonic ionization fronts while revealing shot-to-shot variations in the plasma profiles. Experimental results are validated against theoretical models and scaling laws, underscoring the robustness and accuracy of this technique. By enabling ultrafast, high-resolution plasma diagnostics in a single exposure, SAPPHIRE represents a transformative advancement in plasma measurement technology.
Grace, Elizabeth S. [Lawrence Livermore National L↗
Three-Phase Modeling in Sockeye
The heat pipe simulation tool Sockeye is given the capability to allow the presence of the solid phase of the heat pipe working fluid, which is necessary for modeling transients to and from room temperature, since microreactor designs feature high-temperature heat pipes employing metal working fluids. The liquid phase equations are converted to equations for the mixture of liquid and solid phases, with a simple latent heat relation determining the fractions of these phases. This capability is introduced while preserving equivalent solutions for the liquid-vapor regime. Theory is presented, along with numerical examples to illustrate the propagation of a melt front and startup from the solid state.
97 MATHEMATICS AND COMPUTING↗
Atomic Mixing of Initially Separated Mixtures
We consider the atomic mixing of two materials, initially separated in space, at fixed temperature– pressure conditions [T, p], and zero mass-mean velocity. Each material is composed of a specified set of atomic species whose density corresponds to the specieswise real Equation of State (EOS) $\tilde{ρ}$ i [T, p]. Because there are no spatial gradients in the mean [T, p] state, the only mechanism for mixing is a specieswise nonequilibrium thermal velocity distribution. The recipe for nonequilibrium species mass fluxes is taken directly from Hirschfelder, Curtiss, and Bird, who express the diffusion velocities in terms of gradients in the molar concentrations and pairwise coefficients of diffusion. EOS data from the SESAME library is used for all species. Data for the specieswise self diffusion coefficients is adopted from the Orbital Free Molecular Dynamics (OFMD) simulations of Ticknor et al.; the so-called Darken approximation is used to express the mutual diffusion coefficients of a mixture in the manner used by White et al. A one-dimensional spatial domain is used; with the left half occupied by one multispecies material and the right half occupied by the other, and with zero–gradient boundary conditions at the ends of the domain. Hence the total mass in the domain is a constant. At time t = 0 the materials are separated. Accordingly, spatial gradients at the separation point are infinite and the species fluxes are also infinite. As t → $\textbf{∞}$ the spatial gradients approach zero and a uniform set of species concentrations is obtained; at this point the system is said to be at “thermodynamic equilibrium.” Analytic solution of the problem is not possible in general due to the nonlinear concentration dependence of the species fluxes. Hence we utilize the numerical solution apparatus in Mathematica for obtaining accurate time-dependent solutions for the mixing of various initial species masses for conditions of interest in laser-driven Inertial Confinement Fusion. Concentration profiles and the speed of propagation of the mixing fronts are discussed in light of the nonlinearity of the problem.
73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗