Modeling and testing of a highly sensitive surface acoustic wave pressure sensor for liquid depth measurements
Not Available
Sreejith, Vattaparambil Sreedharan↗
SEARCH · Search NASA
Results for “PRESSURE WAVE”
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.
Numerical simulation of compressible fluid-dynamics in the chamber of inertial fusion energy systems
Here, this paper aims to establish new and innovative modeling capabilities for analyzing chambers in Inertial Fusion Energy (IFE) systems. IFE is emerging as a promising method to achieve fusion power production, but several challenges must be overcome to develop an IFE pilot plant or deploy commercial IFE systems. These challenges are both theoretical and technical, encompassing a deeper understanding of the underlying physical phenomena and the development of new technologies and materials. One of the needs is to develop mathematical models to describe IFE systems and numerical tools to simulate them. This paper contributes to this endeavor by presenting a new OpenFOAM solver for IFE systems, focusing on gas dynamics in their chambers. The analysis and development of chamber designs will play a significant role in the transition from single-shot experiments to high-repetition rates, as there is a need to protect the chamber walls from the intense radiation fields produced by fusion reactions. A promising design option, normally referred to as thick wall chamber design, consists in using lithium or molten salt jet arrays within the chamber. A critical phenomenon is the venting of high-pressure gases from the center to the external part of the chamber, passing through the blanket array. This process involves the propagation and attenuation of strong pressure waves, requiring suitable modeling approaches for compressible fluid-dynamics. The solver proposed in this work implements a multi-material hydrodynamics model tailored to accurately describe the non-linear propagation of pressure waves while avoiding numerical oscillation issues typical of high-velocity compressible simulation. This solver is verified against numerical test cases, validated against experimental data, and applied to the analysis of the High-Yield Lithium-Injection Fusion-Energy (HYLIFE-I) concept. The relevance of this paper is threefold. Firstly, it contributes to developing and testing modeling approaches for compressible fluid-dynamics phenomena, with specific focus on the new and unexplored topic of IFE thick-liquid-wall blanket modeling. Secondly, it marks one of the first applications of the OpenFOAM library in the research field of IFE systems. Finally, the investigated problem is of practical interest for IFE developers, as it provides useful indications about relevant phenomena in pressure wave propagation in the chamber of these systems.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Predicting multiphase flow and tracer transport for an underground chemical explosive test
Detecting radionuclide gas seepage from clandestine underground nuclear tests is central to nonproliferation explosion monitoring research. Yet, early-time (<6 day) gas transport driven by the explosive pressure wave remains poorly constrained due to scarcity of field data. We simulate multi-phase gas transport in the vadose zone using pre-shot data from a recent chemical explosion in P-Tunnel at the Nevada National Security Site, USA. Despite using a simplified 2D-radial model, predictions of tracer arrival matched observations within one order-of-magnitude. Our results show how transient blast forcing rapidly mobilizes gases from the cavity into surrounding rock – critical for optimizing sensor placement and test planning. This unique integration of field data and modeling represents a significant improvement in our ability to predict gas migration from underground explosions. More broadly, it offers insights into the coupled dynamics of pressure waves and contaminant transport in the vadose zone, with implications for monitoring and hazard assessment.
54 ENVIRONMENTAL SCIENCES↗
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↗
Nondestructive Evaluation of Concrete: Elastic Property Imaging Through Full Waveform Inversion
Concrete is a major construction material worldwide and plays a crucial role in the nuclear industry. The elastic properties of concrete are prone to change and degrade while in service, as it is often subjected to extreme operational and environmental conditions. An accurate evaluation of concrete's elastic properties is thus essential to ensure structural integrity and safety. This is especially true for concrete in nuclear power plants, where irradiation effects significantly impact concrete mechanical properties. There are various methods to assess these properties, with ultrasound-based techniques showing high potential due to their nondestructive nature, cost-effectiveness, and safety. While several nondestructive evaluation methods exist, most rely on idealizations such as assuming homogeneous material and plane wavefronts. In this work, we address these issues by introducing an ultrasound-based nondestructive method aimed at reconstructing spatially varying images of concrete mechanical properties. By accurately modeling wave physics, including scattering and reflection, we overcome several of the aforementioned idealizations and aim to utilize the full waveform for imaging material properties through depth, resulting in more reliable images. Full waveform inversion (FWI) was first introduced by geophysicists to reconstruct subsurface elastic property images. The goal is to minimize the difference between simulated and recorded wavefield signals, often through gradient-based optimization algorithms. While FWI is primarily conducted using the acoustic approximation of the wave equation, few works focus on elastic FWI, where the goal is to reconstruct images of not only the pressure wave speed but also the shear wave speed and density (or their equivalents). This work explores the potential of using elastic FWI to predict concrete mechanical properties as an initial effort for a more accurate monitoring of concrete conditions in service. Reconstructing images of different elastic parameters enables more specificity and accurate condition diagnosis. This paper will detail this approach and provide examples demonstrating the effectiveness of elastic FWI in reconstructing comprehensive maps of concrete mechanical properties.
42 - ENGINEERING↗
Shock-induced bubble jets: a dual perspective of bubble collapse and interfacial instability theory
Interactions between shock waves and gas bubbles in a liquid can lead to bubble collapse and high-speed liquid jet formation, relevant to biomedical applications such as shock wave lithotripsy and targeted drug delivery. This study reveals a complex interplay between acceleration-induced instabilities that drive jet formation and radial accelerations causing overall bubble collapse under shock wave pressure. Using high-speed synchrotron X-ray phase contrast imaging, the dynamics of micrometre-sized air bubbles interacting with laser-induced underwater shock waves are visualised. These images offer full optical access to phase discontinuities along the X-ray path, including jet formation, its propagation inside the bubble, and penetration through the distal side. Jet formation from laser-induced shock waves is suggested to be an acceleration-driven process. A model predicting jet speed based on the perturbation growth rate of a single-mode Richtmyer–Meshkov instability shows good agreement with experimental data, despite uncertainties in the jet-driving mechanisms. The jet initially follows a linear growth phase, transitioning into a nonlinear regime as it evolves. To capture this transition, a heuristic model bridging the linear and nonlinear growth phases is introduced, also approximating jet shape as a single-mode instability, again matching experimental observations. Upon piercing the distal bubble surface, jets can entrain gas and form a toroidal secondary bubble. Linear scaling laws are identified for the pinch-off time and volume of the ejected bubble relative to the jet’s Weber number, characterising the balance of inertia and surface tension. At low speeds, jets destabilise due to capillary effects, resulting in ligament pinch-off.
Drops and Bubbles: Bubble dynamics↗
Experimental shock Hugoniot and melt boundary temperatures of aluminum
Aluminum is a ubiquitous component in dynamic compression, pulsed power, and other high energy density physics studies. Its high-pressure behavior and phase diagram are extensively studied standards in shockwave physics. While theoretical calculations and multiphase equations of state have been benchmarked to velocity measurements of loading and unloading waves, pressure and density under shock, and other mechanical data, experimental temperature data under these conditions have not been reported. We conducted a series of experiments shocking and releasing aluminum 6061 and 1100 samples into lithium fluoride windows. We measured temperature at the sample–window interface under steady compression and subsequent isentropic release. These results allow us to constrain the temperature of the solid Hugoniot and the boundary between the liquid and face-centered cubic solid phases.
Hartsfield, Thomas Murray [Sandia National Laborat↗
Nondestructive Evaluation of Concrete: Elastic Property Imaging Through Full-Waveform Inversion
Concrete is a vital material in construction—especially in the nuclear industry, where it is used in critical structures such as containment vessels. Over time, concrete can degrade due to harsh operational and environmental conditions, necessitating that its elastic properties be accurately evaluated to ensure structural integrity and safety. Traditional nondestructive evaluation methods such as ultrasound-based techniques often rely on simplifying assumptions that may not hold true for concrete. This paper presents an advanced ultrasound-based method that uses elastic full-waveform inversion (EFWI) to create detailed images of concrete’s mechanical properties. By accurately modeling wave behaviors such as scattering and reflection, we aim to overcome the limitations of conventional ultrasonic-based methods. In this work, the imaging problem involved reconstructing the various elastic properties of a heterogenous concrete block with three steel rebars embedded in it. The ultrasonic measurements were synthetically generated from multiple sources and receivers, and the reconstruction process was performed using a gradient-based optimization algorithm. Our approach leveraged EFWI to reconstruct high-resolution images of the pressure wave speed, shear wave speed, and density. Multiple misfit functions—including L2-norm, cross-correlation (CC), and L1-norm—combined with total variation (TV) regularization and parameter constraints using a Sigmoid function—were explored for the reconstruction. The results demonstrated that using the L1-norm misfit function in conjunction with TV regularization and Sigmoid constraints significantly improved the reconstruction quality in comparison to traditional methods. This approach provided clearer images with fewer artifacts and better captured background heterogeneity. Our findings highlight that, when properly designed, EFWI carries great potential for providing comprehensive, more accurate, and more reliable assessments of concrete conditions, as is crucial for the maintenance and safety of nuclear power plant structures.
97 - MATHEMATICS AND COMPUTING↗
Characterization of blast waves induced by femtosecond laser irradiation in solid targets
Blast waves have been produced in solid target by irradiation with short-pulse high-intensity lasers. The mechanism of production relies on energy deposition from the hot electrons produced by laser–matter interaction, producing a steep temperature gradient inside the target. Hot electrons also produce preheating of the material ahead of the blast wave and expansion of the target rear side, which results in a complex blast wave propagation dynamic. Several diagnostics have been used to characterize the hot electron source, the induced preheating and the velocity of the blast wave. Results are compared to numerical simulations. These show how blast wave pressure is initially very large (more than 100 Mbar), but it decreases very rapidly during propagation.
46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
Multi-energy ion irradiation effects and distinctive features of color center generation in yttria-stabilized zirconia driven by structural discrepancy
Regarding a requirement for inert matrix fuel (IMF) to burn minor actinides and plutonium and high-level waste immobilization, Yttria-stabilized zirconia (ZrO 2 -6.5 wt% Y 2 O 3 ) (YSZ) crystals with cubic phase are selected to comprehensively evaluate the irradiation resistance by irradiation-induced structural discrepancies across multiple energy conditions. Electronic excitation and nuclear collisions reveal distinct color centers, influenced by varying component ratios that reflect their formation mechanisms. Key factors driving ultra-fast structural transitions include collision cascades, pressure waves, and energy dissipation from electronic excitation. In this work, we find that in the electronic energy loss (E ele )-dominant region, internal latent tracks with unique surface nanostructures emerge, distinct from the defect clusters seen in the nuclear energy loss (E nuc )-dominant region, while the threshold for latent track formation and melting are identified. In conclusion, structural discrepancies driven by different energy-loss mechanisms promote the generation of singly ionized and nearest-neighbor doubly ionized oxygen vacancies, which respectively dominate the formation of F⁺ and T-centers, resulting in bandgap narrowing (4.80 eV → 4.73 eV), and furthermore triggering visible light emission (2.17 eV → 2.20 eV) in irradiated YSZ crystals, thus, enabling the design of novel irradiation-tailored material functionalities.
36 MATERIALS SCIENCE↗
How Did Westward Volcaniclastic Deposits Accumulate in the Deep Sea Following the January 2022 Eruption of Hunga Volcano?
Most volcanic eruptions on Earth take place below the ocean surface and remain largely unobserved. Reconstruction of past submerged eruptions has thus primarily been based on the study of seafloor deposits. Rarely before the 15 January 2022 eruption of Hunga volcano (Kingdom of Tonga) have we been able to categorically link deep-sea deposits to a specific volcanic source. This eruption was the largest in the modern satellite era, producing a 58-km-tall plume, a 20-m high tsunami, and a pressure wave that propagated around the world. The eruption induced the fastest submarine density currents ever measured, which destroyed submarine telecommunication cables and traveled at least 85 km to the west to the neighboring Lau Basin. Here we report findings from a series of remotely operated vehicle dives conducted 4 months after the eruption along the Eastern Lau Spreading Center-Valu Fa Ridge. Hunga-sourced volcaniclastic deposits 7–150 cm in thickness were found at nine sites, and collected. Study of the internal structure, grain size, componentry, glass chemistry, and microfossil assemblages of the cores show that these deposits are the distal portions of at least two ~100-km-runout submarine density currents. We identify distinct physical characteristics of entrained microfossils that demonstrate the dynamics and pathways of the density currents. Microfossil evidence suggests that even the distal parts of the currents were erosive, remobilizing microfossil-concentrated sediments across the Lau Basin. Remobilization by volcaniclastic submarine density currents may thus play a greater role in carbon transport into deep sea basins than previously thought.
58 GEOSCIENCES↗
Effects of spatiotemporal plasma power distribution on the modeling of ignition kernel evolution in quiescent and turbulent methane/air mixtures
Abstract The present work improves a phenomenological plasma-assisted combustion model by integrating the spatiotemporal distribution of plasma power density, thereby considering the evolution of plasma streamers in the modeling, and subsequently, better predicting the ignition kernel evolution. The improved phenomenological model is validated against experiments representing the plasma discharge and post-discharge ignition kernel evolution. Specifically, the new model demonstrates a more accurate prediction of ultrafast gas heating and O 2 dissociation during the plasma discharge, compared to the original model. In addition, the new model is found to closely match the experimental pressure wave and heated channel profiles post-discharge without the need for tuning the energy deposition (unlike the original model), highlighting its accuracy of post-discharge ignition kernel dynamics. The improved phenomenological model is then employed to investigate ignition kernel evolution for a stoichiometric methane-air discharge across various discharge gap configurations. Simulations reveal a non-uniform temperature and streamer distribution progressing from the electrode tips toward the center, contrasting uniform cylindrical discharges previously described in the original model. Streamer propagation is observed to be faster for larger gaps when maintained at the same average electric field for different discharge gaps. The tendency of smaller gaps to produce detached toroidal ignition kernels is observed, while larger gaps promote cylindrical and attached ignition kernels. Interactions between successive ignition kernels from consecutive discharges varied significantly, with the smallest gap (1 mm) promoting the quenching of the preceding ignition kernel due to the initial kernel–kernel separation. The intermediate gap (2 mm) promotes detached kernel growth. In contrast, in the largest gap (4 mm), kernels consistently combine and expand attached to electrodes. The impact of homogeneous isotropic turbulence is also explored, showing the persistence of ignition kernels early on but eventually quenching due to enhanced radical and heat losses with pronounced turbulence intensity.
Johnson, Praise Noah↗
Evolving Outer Heliosphere: Tracking Solar Wind Transients from 1 au to the VLISM with IBEX and Voyager 1
Abstract Interstellar Boundary Explorer (IBEX) observations of energetic neutral atom (ENA) fluxes from the heliosphere have greatly enriched our understanding of the interaction of the solar wind (SW) with the local interstellar medium (LISM). However, there has been recent controversy surrounding the inability of most ENA models to produce as high an intensity of ∼0.5–6 keV ENAs as IBEX observes at 1 au, especially as a function of time. In our previous study (E. J. Zirnstein et al.), we introduced a new model that utilizes a data-driven magnetohydrodynamic simulation of the SW–LISM interaction to propagate pickup ions through the heliosheath (HS) after they are nonadiabatically heated at the heliospheric termination shock. E. J. Zirnstein et al. only simulated and analyzed IBEX observations from the direction of Voyager 2. In this study, we expand our model to include fluxes from the direction of Voyager 1, as well as in the low-latitude part (middle) of the ribbon (10° below the ecliptic plane). We show that the model results at Voyager 1 are consistent with E. J. Zirnstein et al.’s results at Voyager 2 in terms of a secondary ENA source contribution of ≲20% from both directions. Our results in the middle of the ribbon also reproduce the data, when including a time-dependent secondary ENA source. Finally, we demonstrate with our simulation that three large pressure waves likely merged in the VLISM and were observed by Voyager 1 as “pf2,” while at least one of the wave’s effects in the HS was observed by IBEX as a brief enhancement in ENA flux in early 2016.
79 ASTRONOMY AND ASTROPHYSICS↗
Gigacycle Fatigue Strength Evaluation of Welded 316L Stainless Steels for Mercury Target Vessel
At the Materials and Life Science Experimental Facility (MLF) in J-PARC, liquid mercury target for the pulsed spallation neutron source is in operation. An enclosure vessel for the liquid mercury target made of type 316L stainless steel (SS316L) suffers two kinds of cyclic stress during operation. One is the thermal stress due to the internal heating and swings by proton beam trip. The other is the impulsive stress by the pressure waves generated by the proton beam injection. The total number of loading cycles for the former is ∼104, and the latter is ∼4 × 108 for a year operation in the J-PARC mercury target vessel. The target vessel is assembled by an electron beam welding (EBW) and a gas tungsten arc welding (GTAW). However, fatigue data of welded SS316L up to gigacycle is limited. Ultrasonic fatigue testing, applying load cycles by utilizing ultrasonic resonance, for the welded SS316L was performed to investigate the effect of welding on fatigue behavior up to a gigacycle. The result showed that the fatigue strength degradation by EBW and EBW with GTAW were not recognized up to 109 cycles. Crack initiation in welded specimens nucliated at off-center areas of the specimen whereas the cracks in base metal specimen originated at the specimen center.
Naoe, Takashi [Japan Atomic Energy Agency (JAEA)]↗
Pressure suppresses the density wave order in kagome metal LuNb 6 Sn 6
The density waves that develop in kagome metals ScV 6 Sn 6 and LuNb 6 Sn 6 at low temperature appear to arise from underfilled atomic columns within a V-Sn or Nb-Sn scaffolding. Compressing this network with applied pressure in ScV 6 Sn 6 suppressed the structural transition temperature by constraining atomic rattling and inhibiting the shifts that define the structural modulation. We predicted that the density wave transition in LuNb 6 Sn 6 at 68 K would be suppressed by pressure as well. In this Letter, we examine the pressure dependence of the density wave transition by measuring resistance vs temperature up to 2.26 GPa. We found the transition temperature is smoothly depressed and disappears around 1.9 GPa. In conclusion, this result not only addresses our prediction, but strengthens the rattling chains origin of structural instabilities in the HfFe 6 Ge 6 -type kagome metals.
Charge density waves↗
Adaptive Stiffness Power Take-Off System for a Hydraulic Hose-Pump Point Absorber Wave Energy Converter
Two crucial challenges faced in the design of ocean wave energy converters (WECs) are maintaining high capacity factors and remaining robust in everchanging, harsh ocean environments. This study introduces a novel mechanism for achieving reactive WEC control with a soft hydraulic power take-off, thereby addressing both challenges. The device of interest involves a hydraulic pump, constructed from a fluidic flexible matrix composite (F2MC), which is cyclically stretched by driving waves and pressurizes working fluid through a turbine-generator. This class of fiber-reinforced stretch hoses is established for ocean use as robust moorings and shock absorbers. The dynamic properties of the F2MC pump depend on internal pressure, which can be controlled by a continuously variable transmission between the turbine and generator, allowing the system to be adapted for better performance over a wide range of wave conditions. This work presents a dynamic model for operation and adaptive behavior of this power take-off system, experimentally validates components of the model, and performs dynamic simulations of the system that show increases in power production and capacity factor compared to non-adaptive operation for a variety of driving wave conditions. This work develops the basis for future large-scale experimental testing and controller design of this system.
Wave energy conversion, fluidic flexible matrix co↗
Charge-density-wave quantum critical point under pressure in 2 H -TaSe 2
The presence of a quantum-critical point (QCP) at which a nearby ordered phase is suppressed to zero temperature is often invoked to explain emergent quantum phases, e.g. superconductivity. Yet, identifying a QCP and establishing its correlation with superconductivity remains challenging. Materials featuring charge-density-wave (CDW) order and superconductivity offer a clear scenario as both states can be associated with electron-phonon coupling. Here, we uncover a CDW-QCP and demonstrate its interrelation with superconductivity in the prototypical transition-metal dichalcogenide 2 H -TaSe 2 . We determine the evolution of the CDW state up to and beyond its suppression at the critical pressure p c = 19.9(1) GPa by means of X-ray diffraction and inelastic X-ray scattering measurements providing a full crystallographic refinement of the commensurate CDW superstructure. The pressure-induced CDW-QCP in close vicinity to the maximum superconducting transition temperature. Ab-initio lattice dynamical calculations corroborate that 2 H -TaSe 2 features order-parameter fluctuation enhanced superconductivity and can serve as a paradigm to investigate superconductivity near a CDW-QCP.
critical phenomena↗
Chemical timescale effects on detonation convergence
Numerical simulations of detonation-containing flows have emerged as crucial tools for designing next- generation power and propulsion devices. As these tools mature, it is important for the combustion community to properly understand and isolate grid resolution effects when simulating detonations. To this end, the objective of this work is to provide a comprehensive analysis of the numerical convergence of unsteady detonation simulations, with focus on isolating the impacts of chemical timescale modifications on convergence characteristics in the context of operator splitting. With the aid of an AMReX-based adaptive mesh refinement flow solver-which enables resolutions up to ($\mathcal{O}$ (1000) cells-per-induction length-the convergence analysis is conducted using two kinetics configurations: (1) the simplified three-step Arrhenius-based model mechanism of Short and Quirk (1997), where chemical timescales in the detonation are modified by adjusting activation energies in the initiation and branching reactions, and (2) a detailed hydrogen- air mechanism, where the chemical timescales are adjusted by varying the ambient pressure. The convergence of unsteady self-sustained detonations in one-dimensional channels is then analyzed with reference to steady-state theoretical baseline solutions using these mechanisms. The goal of the analysis is to provide a detailed comparison of the effects of grid resolution on both macroscopic (peak pressures and wave speeds) and microscopic (wave structure) quantities of interest, drawing connections between the deviations from steady-state baselines and minimum chemical timescales. In particular, chemical timescale reductions were found to have minimal impact on the convergence of macroscopic properties. Furthermore, analyses of microscopic convergence trends, particularly in the reaction front location, revealed a key insight: maintaining the induction time while eliminating prohibitive chemical timescales through mechanism simplifications and combustion modeling can significantly enhance detonation convergence properties. Ultimately, this work uncovers resolution-dependent unsteady detonation convergence regimes and highlights the important role played by not only the chemical timescales, but also the ratio between the chemical timescale and induction time on the numerical convergence of the detonation wave structure.
Adaptive mesh refinement↗