Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “local oscillator”

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 217 records · Page 12

The reverse quantum limit and its implications for unconventional quantum oscillations in YbB 12

The quantum limit in a Fermi liquid, realized when a single Landau level is occupied in strong magnetic fields, gives rise to unconventional states, including the fractional quantum Hall effect and excitonic insulators. Stronger interactions in metals with nearly localized $f$-electron degrees of freedom increase the likelihood of these unconventional states. However, access to the quantum limit is typically impeded by the tendency of $f$-electrons to polarize in a strong magnetic field, consequently weakening the interactions. In this study, we propose that the quantum limit in such systems must be approached in reverse, starting from an insulating state at zero magnetic field. In this scenario, Landau levels fill in the reverse order compared to regular metals and are closely linked to a field-induced insulator-to-metal transition. We identify YbB 12 as a prime candidate for observing this effect and propose the presence of an excitonic insulator state near this transition.

36 MATERIALS SCIENCE↗

A Reversed Impedance-Based Stability Criterion for IBR Grids

The existing impedance-based stability criterion is effective for analyzing local control interactions; however, it is difficult to scale the existing criterion to analyze wide-area control interactions among numerous IBRs through a complex power system network. The scaled version of the existing criterion requires the impedance response of each IBR in the system as well as of the network looking from all the IBRs. It is quite challenging to obtain all these impedance responses because of the computational effort and the requirement of separately scanning the impedance of the network and the IBRs. We propose a reversed criterion for the impedance-based stability analysis to address these problems. In contrast to the existing criterion, the reversed criterion analyzes the stability of a power system when an IBR is disconnected from the system. The reversed criterion estimates the impact of an IBR on the frequency and damping of power system oscillation modes using the impedance scans of only the IBR and the grid at its terminal. It can be sequentially applied at different IBRs to evaluate their impact on the power system stability. In addition to scalability, the reversed criterion gives flexibility to focus only on a few selected IBRs, depending on their rating, the magnitude of oscillations observed at their terminals, and the vendor support available for implementing stabilizing control system updates. The reversed criterion is demonstrated on a 14-bus power system with 100% IBRs.

control interactions↗

Positive Neutrino Masses with DESI DR2 via Matter Conversion to Dark Energy

The Dark Energy Spectroscopic Instrument (DESI) is a massively parallel spectroscopic survey on the Mayall telescope at Kitt Peak, which has released measurements of baryon acoustic oscillations determined from over 14 million extragalactic targets. We combine DESI Data Release 2 with CMB datasets to search for evidence of matter conversion to dark energy (DE), focusing on a scenario mediated by stellar collapse to cosmologically coupled black holes (CCBHs). In this physical model, which has the same number of free parameters as Λ⁢CDM, DE production is determined by the cosmic star formation rate density (SFRD), allowing for distinct early- and late-time cosmologies. Using two SFRDs to bracket current observations, we find that the CCBH model: accurately recovers the cosmological expansion history, agrees with early-time baryon abundance measured by BBN, reduces tension with the local distance ladder, and relaxes constraints on the summed neutrino mass ∑𝑚 𝜈 . For these SFRDs, we find a peaked positive ∑𝑚 𝜈 < 0.149 eV (95% confidence) and ∑𝑚 𝜈 = 0.106$^{+0.050}_{−0.069}$ eV, respectively, in good agreement with lower limits from neutrino oscillation experiments. A peak in ∑𝑚 𝜈 > 0 results from late-time baryon consumption in the CCBH scenario and is expected to be a general feature of any model that converts sufficient matter to dark energy during and after reionization.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

A Multiagent Deep Reinforcement Learning-Enabled Dual-Branch Damping Controller for Multimode Oscillation

Here, this study develops a multiagent deep reinforcement learning (MADRL)-enabled framework for the decentralized cooperative control of a novel dual-branch (DB) damping controller for both low-frequency oscillation (LFO) and ultralow-frequency oscillation (ULFO). It has two branches, each of which consists of a proportional resonance (PR) and a second-order polynomial that is designed to handle target oscillation modes. To improve the robustness of the controller to system uncertainties, MADRL is developed, where multiagents are centrally trained to obtain the coordinated adaptive control policy while being executed in a decentralized manner to provide the optimal parameter setting for each controller with only local states. Comparisons with the IEEE 10-machine 39-bus system demonstrate that the proposed method achieves better robustness to uncertainties, lower communication delay, and single-point failure, as well as damping control performances for both LFO and ULFO.

97 MATHEMATICS AND COMPUTING↗

Medium-Range Ordering in the Ionic Glass Electrolytes LiPON and LiSiPON

In this work, we provide an in-depth structural characterization of the amorphous ionic glasses LiPON and LiSiPON with high Li content. Based on ab initio molecular dynamics simulations, the structure of these materials is an inverted structure with either isolated polyanion tetrahedra or polyanion dimers in a Li+ matrix. Based on neutron scattering data, this type of inverted structure leads to a significant amount of medium-range ordering in the structure, as demonstrated by two sharp diffraction peaks and a periodic structural oscillation in the density function G(r). While this medium-range ordering is commonly observed in liquids and metallic glasses, it has not previously been observed in oxides. On a local scale, adding N and Si increases the number of anion bridges and polyanion dimer structures, leading to higher ionic conductivity. In the medium-range ordering, the addition of Si leads to more disorder in the polyanion substructure but a significant increase in the ordering of the O substructure. Finally, we demonstrate that this inverted structure with medium-range ordering results in a glassy material that is both mechanically stiff and ductile on the nanoscale.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Reversible Phasonic Control of a Quantum Phase Transition in a Quasicrystal

Periodic driving can tune the quasistatic properties of quantum matter. A well-known example is the dynamical modification of tunneling by an oscillating electric field. Here we show experimentally that driving the phasonic degree of freedom of a cold-atom quasicrystal can continuously tune the effective quasidisorder strength, reversibly toggling a localization-delocalization quantum phase transition. Measurements agree with fit-parameter-free theoretical predictions, and illuminate a fundamental connection between Aubry-André localization in one dimension and dynamic localization in the associated two-dimensional Harper-Hofstadter model. These results open up new experimental possibilities for dynamical coherent control of quantum phase transitions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Grid-Connected Modular Soft-Switching Solid State Transformers (M-S4T)

The objective of this project is to develop and verify the concept of a flexible and modular soft-switching solid-state transformer (M-S4T) for direct grid-connected applications. The ability to directly connect power electronics converters to the medium voltage grid (4 kV – 13 kV), and to potentially replace the passive and bulky, but ubiquitous 60 hertz service transformer in the 25 kVA to 100 kVA range, with a more flexible and controllable device, has been regarded as the ‘holy grail’ in grid control. However, this has proven to be extremely difficult. This project has developed the solutions to several key challenges of the direct grid-connected power electronics and realized a 7.2 kV M-S4T prototype. First, a protection method to protect the M-S4T from the high voltages (110 kV for the 13 kV system) that occur on the grid due to transients and lightning strikes have been developed and experimentally verified. Second, the realization and the operation of the M-S4T based on high-voltage SiC devices (>3.3 kV) and a medium-frequency medium-voltage low-leakage transformer in a single-stage solid-state transformer with zero-voltage switching, low dv/dt, and low electromagnetic interference has been successfully demonstrated up to 7.5 kV peak. Third, an oil-cooling system and stable communication and distributed control system for converter module voltage sharing have been developed and experimentally verified. The developed M-S4T has realized a modular universal high-performance power conversion system. This conversion system is scalable to different voltage and power levels and adaptable to four-quadrant bidirectional operation. Moreover, the use of passive cooling techniques meets the equipment life requirements, and the lightning protection scheme fulfills the basic insulation level specifications for direct grid connection. Such power conversion system opens up near-term opportunities, including energy storage, solar PV, or electric vehicle charging with significant cost and footprint savings. In the longer term, the possibility of replacing the utility distribution transformer with an M-S4T will be transformative for future distribution grids with a compact footprint and full controllability to enable high renewable energy and storage penetration. In addition to the main project, this report expands on the Plus-Up projected including as part of the main award. This project developed and demonstrated the technology for autonomous collaborative inverters that can be connected in an ad hoc manner to the grid. The aim of the project was to: (1) evaluate the existing techniques for grid-connected inverters and find their limitations; (2) develop detailed requirements for grid-connected inverters in the modern grid with millions of active nodes; (3) design a unified control strategy that brings more autonomy and intelligence to grid-connected inverters, and addresses parts of the issues with the existing techniques. The proposed technique, called UniCon, enables inverters to 1) connect/disconnect to/from the grid in an ad hoc manner; (2) work based on local sensing. Slow communication could be used for a more optimized behavior; (3) work automatically in both grid-forming/grid-following mode; (4) handle large disturbances, e.g., big load step and fault, in an oscillation-free manner; (5) work collaboratively with other inverters in steady-state and during transients. UniCon can be implemented in the middle-level control; hence it is agnostic to the vendor and to the implementation of the inner voltage/current and protection loops. Furthermore, a new synchronization scheme, based on deep learning, was developed that can extract the grid voltage phase and amplitude in a stable manner. The method is cheap to implement can improve the dynamic performance of the grid-connected inverters during fast transients, e.g., fault. The proposed control scheme was validated by (1) MATLAB/Simulink; (2) hardware-in-the-loop results, and; (3) experimental results using three inverters that form a microgrid in a down-scaled feeder. Lastly, both the M-S4T and UniCon have achieved promising tangible paths to markets. In the case of the M-S4T, the underlying technology — the Soft Switching Solid State Transformer (S4T) developed at the Georgia Tech Center for Distributed Energy (GT-CDE) has been licensed by GridBlock from the Georgia Tech Research Corporation, and GridBlock has been working with manufacturing partner Jabil (one of the largest US-based contract manufacturers) and system integrator Power Secure (largest deployer of microgrids in the US with 4.7 GW under management), to meet the strong initial demand. Similarly, GridBlock has an exclusive license to the UniCon technology, developed under this award by GT-CDE. The UniCon provides an intermediate control layer that enables the implementation of the higher-level ‘transactive’ control commands for the system. The architecture of the system - slow communications with the cloud for system optimization and setpoints, and the use of locally measured quantities for real-time control, provide a very robust and secure way of implementing a real-time must-run grid that is also secure and stable. This is a brand-new functionality that is critical for the future grid and key to GridBlock’s business model.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Exact and locally implicit source term solvers for multifluid-Maxwell systems

Recently, a family of models that couple multifluid systems to the full Maxwell equations have been used in laboratory, space, and astrophysical plasma modeling. These models are more complete descriptions of the plasma than reduced models like magnetohydrodynamic (MHD) since they are derived more closely from the full kinetic Vlasov-Maxwell system, without assumptions like quasi-neutrality, negligible electron mass, etc. Thus these models naturally retain non-ideal MHD effects like electron inertia, Hall term, pressure anisotropy/nongyrotropy, displacement current, among others. One obstacle to broader application of these model is that an explicit treatment of their source terms leads to the need to resolve rapid processes like plasma oscillation and electron cyclotron motion, even when these are not important. In this paper, we suggest two ways to address this issue. First, we derive the analytic solutions to the source update equations, which can be implemented as a practical, but less generic solver. We then develop a time-centered, locally implicit algorithm to update the source terms, allowing stepping over the fast kinetic time-scales. For a plasma with S species, the locally implicit algorithm involves inverting a local (3 S + 3) × (3 S + 3) matrix only, thus is very efficient. The performance can be further increased by using the direct update formulas to skip null calculations. In this paper, we present benchmarks illustrating the exact energy-conservation of the locally implicit solver, as well as its efficiency and robustness for both small-scale, idealized problems and largescale, complex systems. The locally implicit algorithm can be also easily extended to include other local sources, like collisions and ionization, which are difficult to solve analytically.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Real-time tracking and analysis of gas bubble dynamics in laser powder bed fusion using in-situ X-ray characterization and machine learning

Porosity defects remain a significant challenge in the laser powder bed fusion (LPBF) process, adversely affecting the mechanical properties and reliability of additively manufactured components. Here, this study investigates the real-time formation and trajectory of gas bubbles during LPBF of Al6061 alloy using advanced in-situ X-ray characterization and machine learning. The unsupervised Gaussian mixture model and particle tracking algorithm developed are able to precisely track and quantify the properties of gas bubbles and keyhole pores. Our analysis identified five distinct types of gas bubble formation and movement patterns, emphasizing the diverse origins and behaviors of these defects. It enables precise quantification of trajectories, velocities, and morphological changes of gas bubbles, offering a granular view of the subsurface dynamics within the melt pool. Additionally, we explored keyhole-induced pore dynamics, revealing the critical role of keyhole oscillation and collapse for the formation of both large and small gas pores. It defines four different regions of gas bubble movement within the melt pool, providing a clearer understanding of how local fluid dynamics affect pore behavior. The results underscore the importance of integrating in-situ experimental observation and automated machine learning to develop a more robust predictive model for defect formation in LPBF.

In-situ X-ray imaging↗

Optimizing Magnetic Force Microscopy Resolution and Sensitivity to Visualize Nanoscale Magnetic Domains

Magnetic force microscopy (MFM) enables mapping local magnetic fields across a sample surface with nanoscale resolution. To perform MFM, an atomic force microscopy (AFM) probe whose tip has been magnetized vertically (i.e., perpendicular to the probe cantilever) is oscillated at a fixed height above the sample surface. The resultant shifts in the oscillation phase or frequency, which are proportional to the magnitude and sign of the vertical magnetic force gradient at each pixel location, are then tracked and mapped. Although the spatial resolution and sensitivity of the technique increases with decreasing lift height above the surface, this seemingly straightforward path to improved MFM images is complicated by considerations such as minimizing topographical artifacts due to shorter range van der Waals forces, increasing the oscillation amplitude to further improve sensitivity, and the presence of surface contaminants (in particular water due to humidity under ambient conditions). In addition, due to the orientation of the probe's magnetic dipole moment, MFM is intrinsically more sensitive to samples with an out-of-plane magnetization vector. Here, high-resolution topographical and magnetic phase images of single and bicomponent nanomagnet artificial spin-ice (ASI) arrays obtained in an inert (argon) atmosphere glovebox with <0.1 ppm O 2 and H 2 O are reported. Further, optimization of lift height and drive amplitude for high resolution and sensitivity while simultaneously avoiding the introduction of topographical artifacts is discussed, and detection of the stray magnetic fields emanating from either end of the nanoscale bar magnets (~250 nm long and <100 nm wide) aligned in the plane of the ASI sample surface is shown. Likewise, using the example of a Ni-Mn-Ga magnetic shape memory alloy (MSMA), MFM is demonstrated in an inert atmosphere with magnetic phase sensitivity capable of resolving a series of adjacent magnetic domains each ~200 nm wide.

47 OTHER INSTRUMENTATION↗

Convergence analysis for a nonlocal gradient descent method via directional Gaussian smoothing

We analyze the convergence of a nonlocal gradient descent method for minimizing a class of high-dimensional non-convex functions, where a directional Gaussian smoothing (DGS) is proposed to define the nonlocal gradient (also referred to as the DGS gradient). The method was first proposed in [Zhang et al., Enabling long-range exploration in minimization of multimodal functions, UAI 2021], in which multiple numerical experiments showed that replacing the traditional local gradient with the DGS gradient can help the optimizers escape local minima more easily and significantly improve their performance. However, a rigorous theory for the efficiency of the method on nonconvex landscape is lacking. In this work, we investigate the scenario where the objective function is composed of a convex function, perturbed by deterministic oscillating noise. We provide a convergence theory under which the iterates exponentially converge to a tightened neighborhood of the solution, whose size is characterized by the noise wavelength. Here, we also establish a correlation between the optimal values of the Gaussian smoothing radius and the noise wavelength, thus justifying the advantage of using moderate or large smoothing radii with the method. Furthermore, if the noise level decays to zero when approaching the global minimum, we prove that DGS-based optimization converges to the exact global minimum with linear rates, similarly to standard gradient-based methods in optimizing convex functions. Several numerical experiments are provided to confirm our theory and illustrate the superiority of the approach over those based on the local gradient.

Tran, Hoang [Oak Ridge National Laboratory (ORNL),↗

Terahertz-Driven Local Dipolar Correlation in a Quantum Paraelectric

Light-induced ferroelectricity in quantum paraelectrics is a new avenue of achieving dynamic stabilization of hidden orders in quantum materials. In this Letter, we explore the possibility of driving a transient ferroelectric phase in the quantum paraelectric KTaO 3 via intense terahertz excitation of the soft mode. We observe a long-lived relaxation in the terahertz-driven second harmonic generation (SHG) signal that lasts up to 20 ps at 10 K, which may be attributed to light-induced ferroelectricity. Through analyzing the terahertz-induced coherent soft-mode oscillation and finding its hardening with fluence well described by a single-well potential, we demonstrate that intense terahertz pulses up to 500 kV=cm cannot drive a global ferroelectric phase in KTaO 3 . Instead, we find the unusual long-lived relaxation of the SHG signal comes from a terahertz-driven moderate dipolar correlation between the defect-induced local polar structures. We discuss the impact of our findings on current investigations of the terahertz-induced ferroelectric phase in quantum paraelectrics.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Charge-neutral fermions and magnetic field-driven instability in insulating YbIr 3 Si 7

Kondo lattice materials, where localized magnetic moments couple to itinerant electrons, provide a very rich backdrop for strong electron correlations. They are known to realize many exotic phenomena, with a dramatic example being recent observations of quantum oscillations and metallic thermal conduction in insulators, implying the emergence of enigmatic charge-neutral fermions. Here, we show that thermal conductivity and specific heat measurements in insulating YbIr 3 Si 7 reveal emergent neutral excitations, whose properties are sensitively changed by a field-driven transition between two antiferromagnetic phases. In the low-field phase, a significant violation of the Wiedemann-Franz law demonstrates that YbIr 3 Si 7 is a charge insulator but a thermal metal. In the high-field phase, thermal conductivity exhibits a sharp drop below 300 mK, indicating a transition from a thermal metal into an insulator/semimetal driven by the magnetic transition. These results suggest that spin degrees of freedom directly couple to the neutral fermions, whose emergent Fermi surface undergoes a field-driven instability at low temperatures.

36 MATERIALS SCIENCE↗

RKKY coupled local-moment magnetism in NaFe 1 – x Cu x As

A central question in a large class of strongly correlated electron systems, including heavy fermion compounds and iron pnictides, is the identification of different phases and their origins. It has been shown that the antiferromagnetic (AFM) phase in some heavy fermion compounds is induced by Ruderman-Kittel-Kasuya-Yosida (RKKY) interaction between localized moments, and that the competition between this interaction and Kondo effect is responsible for quantum criticality. However, conclusive experimental evidence of the RKKY interaction in pnictides is lacking. Here, using high resolution 23 Na NMR measurements on lightly Cu-doped metallic single crystals of NaFe 1–x Cu x As (x ≈ 0.01) and numerical simulation, we show direct evidence of the RKKY interaction in this pnictide system. Aided by computer simulation, we identify the 23 Na NMR satellite resonances with the RKKY oscillations of spin polarization at Fe sites. Our 23 Na spin-lattice and spin-spin relaxation data exhibit a signature of an itinerant and inhomogeneous AFM phase in this system, accompanied by a simulation of Cu-induced perturbation to the ordered moments on the Fe sites. Our NMR results indicate coexistence of local and itinerant magnetism in lightly Cu-doped NaFe 1–x Cux As.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Understanding Predictability of Daily Southeast U.S. Precipitation Using Explainable Machine Learning

Abstract We investigate the predictability of the sign of daily southeastern U.S. (SEUS) precipitation anomalies associated with simultaneous predictors of large-scale climate variability using machine learning models. Models using index-based climate predictors and gridded fields of large-scale circulation as predictors are utilized. Logistic regression (LR) and fully connected neural networks using indices of climate phenomena as predictors produce neither accurate nor reliable predictions, indicating that the indices themselves are not good predictors. Using gridded fields as predictors, an LR and convolutional neural network (CNN) are more accurate than the index-based models. However, only the CNN can produce reliable predictions that can be used to identify forecasts of opportunity. Using explainable machine learning we identify which variables and grid points of the input fields are most relevant for confident and correct predictions in the CNN. Our results show that the local circulation is most important as represented by maximum relevance of 850-hPa geopotential heights and zonal winds to making skillful, high-probability predictions. Corresponding composite anomalies identify connections with El Niño–Southern Oscillation during winter and the Atlantic multidecadal oscillation and North Atlantic subtropical high during summer.

Pegion, Kathy↗

Estimating the local dark matter density in a non-axisymmetric wobbling disc

The density of dark matter near the Sun, ρDM, ⊙, is important for experiments hunting for dark matter particles in the laboratory, and for constraining the local shape of the Milky Way’s dark matter halo. Estimates to date have typically assumed that the Milky Way’s stellar disc is axisymmetric and in a steady-state. Yet the Milky Way disc is neither, exhibiting prominent spiral arms and a bar, and vertical and radial oscillations. Here, we assess the impact of these assumptions on determinations of ρDM, ⊙ by applying a free-form, steady-state, Jeans method to two different N-body simulations of Milky Way-like galaxies. In one, the galaxy has experienced an ancient major merger, similar to the hypothesized Gaia–Sausage–Enceladus; in the other, the galaxy is perturbed more recently by the repeated passage and slow merger of a Sagittarius-like dwarf galaxy. We assess the impact of each of the terms in the Jeans–Poisson equations on our ability to correctly extract ρDM, ⊙ from the simulated data. We find that common approximations employed in the literature – axisymmetry and a locally flat rotation curve – can lead to significant systematic errors of up to a factor ~1.5 in the recovered surface mass density ~2 kpc above the disc plane, implying a fractional error on ρDM, ⊙ of the order of unity. However, once we add in the tilt term and the rotation curve term in our models, we obtain an unbiased estimate of ρDM, ⊙, consistent with the true value within our 95 percent confidence intervals for realistic 20 percent uncertainties on the baryonic surface density of the disc. Other terms – the axial tilt, 2nd Poisson and time-dependent terms – contribute less than 10 percent to ρDM, ⊙ (given current data) and can be safely neglected for now. In the future, as more data become available, these terms will need to be included in the analysis.

79 ASTRONOMY AND ASTROPHYSICS↗

Interplay of Kondo physics with incommensurate charge density waves in CeTe 3

CeTe 3 is a 2-dimensional (2D) Van der Waals (VdW) material with incommensurate charge density waves (CDW), extremely high transition temperature (T CDW ), and a large momentum-dependent CDW gap that leaves a significant portion of the Fermi surface intact. It is also considered to be a weak Kondo system, a property unexpected for a material with incommensurate CDW, where each atomic site is slightly different. Here, we study the properties of the CDW state in several RTe3 (R is rare earth) materials and examine the hybridization of itinerant states with the localized Ce 4f multiplet in CeTe 3 by using angle resolved photoemission spectroscopy. We find that the renormalization of the itinerant states originating from the hybridization with the deeper localized 4f states at −260 meV is k−dependent and extends to the Fermi level. As these localized states are far from the Fermi level, the observed hybridization affects the effective masses only marginally and does not lead to heavy fermions. However, since the same renormalizing mechanism normally leads to the heavy fermion physics when the localized 4f states are near the Fermi level, our observation of its strong k−dependence suggests that this could be the reason for discrepancy between the heavy masses in specific heat and light ones in Shubnikov de Haas oscillations, often observed in heavy fermions.

36 MATERIALS SCIENCE↗

Adaptive clipping‐and‐redistribution algorithms for bounded and conservative high‐order interpolations applied to discontinuous and reactive flows

Abstract A new adaptive clipping‐and‐redistribution method is presented which provides bounds‐preservation for multidimensional interpolation in the context of high‐order finite‐volume discretizations with adaptive mesh refinement (AMR). The underlying finite‐volume method (FVM) for the computational fluid dynamics applications is fourth‐order accurate for smooth solutions and utilizes AMR for computational efficiency in solving multiscale problems involving turbulence and combustion. High‐order interpolation between different AMR levels is required. However, this operation often leads to numerical issues because combustion species must have physical bounds preserved. The present study overcomes two major challenges in the development of the high‐order interpolation method. First, the method needs to be bound‐preserving near extrema or discontinuities to prevent the emergence of unphysical oscillations while maintaining fourth‐order accuracy in smooth flows. Second, the method needs to satisfy the conservation requirement in multiple dimensions, particularly in the context of curvilinear coordinate transformations. Additionally, the method is designed to be localized and computationally inexpensive. The new interpolation scheme is demonstrated by solving reacting flows, which are extremely sensitive to unphysical overshoots in conserved quantities. The test problems are shock‐induced ‐ combustion and a ‐air flame in a practical bluff‐body combustor. Results show the method prevents new extrema near discontinuities while maintaining high‐order accuracy in smooth regions. In particular, the method is extremely beneficial for combustion with stiff chemistry. With the proposed new method, even if flame fronts cross AMR interfaces or new grids are created in the vicinity of the flame, solution stability is retained.

97 MATHEMATICS AND COMPUTING↗