Results for “Simulation framework” — Page 28

DOE OSTI2025

Modeling Offshore Wind Farm Performance in Coastal Low-Level Jets Using Coupled Mesoscale-Microscale Large Eddy Simulations

Accurately predicting wind farm reliability under complex offshore atmospheric conditions remains a key challenge, particularly during noncanonical meteorological events such as coastal low-level jets (LLJs). LLJs, characterized by strong nonmonotonic vertical shear and directional veer, depart significantly from the simplified inflow assumptions embedded in conventional design standards, low-fidelity engineering models, and microscale large eddy simulations of the atmospheric boundary layer. In this work, we use the virtual wind farm framework—an exascale, graphics processing unit–accelerated large eddy simulation platform coupled with high-fidelity aeroservoelastic turbine models and advanced mesoscale-microscale coupling via the ExaWind software stack—to investigate turbine responses under realistic LLJ forcing. Simulations are performed over the U.S. North Atlantic offshore domain with the use of meteorological inputs from New York State Energy Research and Development Authority buoy data, focusing on a representative LLJ case impacting the International Energy Agency 15 MW reference turbine. Our results show that LLJs can cause up to 50% power deficits in downstream turbine rows and significantly amplify low-speed shaft and tower loads through nonlinear coupling between complex inflow characteristics and turbine structural dynamics. Two primary mechanisms drive these load amplifications: (1) unique LLJ inflow features—including veer and vertical/lateral shear—and (2) the downstream evolution of the flow under stable thermal stratification, which suppresses turbulence mixing and alters wake recovery. These mechanisms produce streamwise variations in turbine loading not captured by standard hub height–based metrics or existing design load case (DLC) definitions. This study highlights the critical role of rotor-scale flow gradients in driving fatigue and system-level aeroelastic responses, challenging current DLC and control strategies. We advocate the integration of full-flow field, environment-aware wind inputs into load modeling and control algorithms. By leveraging exascale computing to resolve mesoscale-microscale coupling, this work lays the groundwork for next-generation offshore wind turbine design and operation in meteorologically complex marine environments.

17 WIND ENERGY↗

DOE OSTI2026

Reactive Transport Modeling with Physics-Informed Machine Learning for Critical Minerals Applications

This study presents a physics-informed neural network (PINN) framework for reactive transport modeling for simulating fast bimolecular reactions in porous media. Accurate characterization of cAhemical interactions and product formation in surface and subsurface environments is essential for advancing critical mineral extraction and related geoscience applications. The proposed methodology sequentially addresses the flow and diffusion–reaction subproblems. The flow field is computed using a mixed formulation, while the diffusion–reaction system is modeled via two uncoupled tensorial diffusion equations reformulated in terms of chemical invariants. PINNs are employed to solve the governing equations, enabling data-efficient, mesh-free prediction of chemical concentration fields. The framework is validated through a series of benchmark problems involving flow in heterogeneous porous media. Initial verification is conducted using patch tests for the flow field, followed by validation of the transport problem with emphasis on preserving non-negativity of concentrations. The complete fast bimolecular reaction scenario is then solved, yielding spatial distributions of reactants and product species. Results demonstrate that the PINNs-based approach effectively captures sharp, mixing-limited reaction fronts and dispersive mixing behavior, offering reliable predictions of reactive plume evolution. These capabilities are crucial for evaluating long-term subsurface behavior in applications such as fluid storage, energy extraction, and efficient extraction of critical minerals.

42 ENGINEERING↗

DOE OSTI2025

Computing chemical potentials with machine-learning-accelerated simulations to accurately predict thermodynamic properties of molten salts

The successful design and deployment of next-generation nuclear technologies heavily rely on thermodynamic data for relevant molten salt systems. However, the lack of accurate force fields and efficient methods has limited the quality of thermodynamic predictions from atomistic simulations. Here we propose an efficient free energy framework for computing chemical potentials, which is the central free energy quantity behind many thermodynamic properties. We accelerate our simulations without sacrificing accuracy by using machine learning interatomic potentials trained on density functional theory (DFT) data. Using lithium chloride as our model system, we compute chemical potentials with DFT-accuracy for solid and liquid phases by transmuting ions into noninteracting particles. Notably, in the liquid phase, we demonstrate consistency whether we transmute one ion pair or the entire system into ideal gas particles. By locating the temperature where the chemical potential of solid and liquid phases cross, we predict a melting point of 880 ± 18 K for lithium chloride, which is remarkably close to the experimental value of 883 K. With this successful demonstration, we lay the foundation for high-throughput thermodynamic predictions of many properties that can be derived from the chemical potentials of the minority and majority components in molten salts.

Gibson, Luke D. [Oak Ridge National Laboratory (OR↗

DOE OSTI2024

Assay-based background projection for the Majorana Demonstrator using Monte Carlo uncertainty propagation

The background index (BI) is an important quantity to project and calculate the half-life sensitivity of neutrinoless double-𝛽 decay (0⁢𝜈⁢𝛽⁢𝛽) experiments. An analysis framework is presented to calculate the BI using the specific activities, masses, and simulated efficiencies of an experiments components as distributions. This Bayesian framework includes a unified approach to combine specific activities from assay. Monte Carlo uncertainty propagation is used to build a BI distribution from the specific activity, mass, and efficiency distributions. This method is applied to the M AJORANA D EMONSTRATOR , which deployed arrays of high-purity Ge detectors enriched in 76 Ge to search for 0⁢𝜈⁢𝛽⁢𝛽. The original assay-based projection is requantified in the new framework, using the as-built geometry of the Demonstrator and additional assay information. While 47% higher than the original projection, the resulting BI of [8.95±0.36]×10 −4 cts/(keVkgyr) from the 232 Th and 238 U decay chains does not account for the higher-than-expected BI observed by the D EMONSTRATOR . Finally, this method enables us to demonstrate the statistical incompatibility between the D EMONSTRATOR 's observed background and the assay results.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2024

A General Framework for Error-controlled Unstructured Scientific Data Compression

Data compression plays a key role in reducing storage and I/O costs. Traditional lossy methods primarily target data on rectilinear grids and cannot leverage the spatial coherence in unstructured mesh data, leading to suboptimal compression ratios. We present a multi-component, error-bounded compression framework designed to enhance the compression of floating-point unstructured mesh data, which is common in scientific applications. Our approach involves interpolating mesh data onto a rectilinear grid and then separately compressing the grid interpolation and the interpolation residuals. This method is general, independent of mesh types and typologies, and can be seamlessly integrated with existing lossy compressors for improved performance. We evaluated our framework across twelve variables from two synthetic datasets and two real-world simulation datasets. The results indicate that the multi-component framework consistently outperforms state-of-the-art lossy compressors on unstructured data, achieving, on average, a 2.3 − 3.5× improvement in compression ratios, with error bounds ranging from 1 × 10 the −6 to 1×10−2. We further investigate impact of hyperparameters, such as grid spacing and error allocation, to deliver optimal compression ratios in diverse datasets.

Gong, Qian↗

DOE OSTI2025

Initial Engineering Framework for Species Tracking in MSRs Involving Fuel Salt and Structures

This report implements a high-fidelity multiphysics modeling framework using the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program tools to track leaching and plating of materials between the fuel salt and structures in Molten Salt Reactors (MSRs). Specifically, this framework is demonstrated in modeling various processes within MSRs including hot and cold leg corrosion, redox potential driven corrosion, colloid precipitation and deposition, and noble metal plating and decay heat modeling. The model integrates neutronics, thermal-hydraulics, depletion, and thermochemistry to simulate the evolving chemical behavior of the fuel salt in interactions with structural materials. Initial findings suggest that hot-cold leg corrosion in MSR systems may become significant if the redox potential of the fuel salt is not sufficiently controlled. Additionally, the impact of noble metal decay heat deposition in large power reactors is explored. The resulting engineering framework for species tracking will be applied in the future toward specific validation work with ongoing experimental efforts focusing on both chloride and fluoride salt loop experiments with in-situ corrosion and redox control instrumentation.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2025

Materials Learning Algorithms (MALA): Scalable machine learning for electronic structure calculations in large-scale atomistic simulations

We present the Materials Learning Algorithms (MALA) package, a scalable machine learning framework designed to accelerate density functional theory (DFT) calculations suitable for large-scale atomistic simulations. Using local descriptors of the atomic environment, MALA models efficiently predict key electronic observables, including local density of states, electronic density, density of states, and total energy. The package integrates data sampling, model training and scalable inference into a unified library, while ensuring compatibility with standard DFT and molecular dynamics codes. We demonstrate MALA's capabilities with examples including boron clusters, aluminum across its solid-liquid phase boundary, and predicting the electronic structure of a stacking fault in a large beryllium slab. Scaling analyses reveal MALA's computational efficiency and identify bottlenecks for future optimization. With its ability to model electronic structures at scales far beyond standard DFT, MALA is well suited for modeling complex material systems, making it a versatile tool for advanced materials research.

Density functional theory↗

DOE OSTI2024

How Well Can Quantum Embedding Method Predict the Reaction Profiles for Hydrogenation of Small Li Clusters?

Quantum computing leverages the principles of quantum mechanics in novel ways to tackle complex chemistry problems that cannot be accurately addressed using traditional quantum chemistry methods. However, the high computational cost and available number of physical qubits with high fidelity limit its application to small chemical systems. This work employed a quantum-classical framework which features a quantum active space-embedding approach to perform simulations of chemical reactions that require up to 14 qubits. This framework was applied to prototypical example metal hydrogenation reactions: the coupling between hydrogen and Li 2 , Li 3 , and Li 4 clusters. Particular attention was paid to the computation of barriers and reaction energies. The predicted reaction profiles compare well with advanced classical quantum chemistry methods, demonstrating the potential of the quantum embedding algorithm to map out reaction profiles of realistic gas-phase chemical reactions to ascertain qualitative energetic trends. Additionally, the predicted potential energy curves provide a benchmark to compare against both current and future quantum embedding approaches.

36 MATERIALS SCIENCE↗

DOE OSTI2024

Simulation of Li-pellet triggered ELMs in EAST with an impurity model implemented under BOUT++ framework

A simple impurity model has been developed under the BOUT++ framework to investigate the Li-pellet triggered edge localized mode (ELM) in EAST configuration. The present impurity model decouples the ion pressure enhancement effect (IPEE) (which affects radial E X B flow shear and gyro-viscosity of deuterium ions and is implicitly included in our previous work) to make the simulation of Li pellet injection more physically reasonable. In addition, the present impurity model also includes the impurity equilibrium effect (IEE) (which induces modifications on vorticity and gyro-viscosity). The simulation results show that without IPEE a turbulent ELM induced by Li pellet is triggered by multiple peeling-ballooning modes (PBMs) rather than by a single dominant mode, and the simultaneous growth of multiple modes is conducive to reducing the time for pedestal entering the energy loss state. When a turbulent ELM occurs, the nonlinearly dominant modes undergo a secondary fast growth during fast crash phase. The finding explains the secondary increase of D α emission observed in the DIII-D Li injection experiments. During the evolution of a turbulent ELM, the continuous pedestal collapse dominates the pedestal energy loss in turbulent transport phase; however, in saturation phase, the stabilizing effects on PBMs (n ≠ 0) by n = 0 mode and radial E X B flow shear within the steepest gradient region are more prominent. It is also found that with IEE and without IPEE, the size threshold of ELM triggering and the magnitude of ELM sizes both show a good agreement with EAST and DIII-D experimental observations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2025

Improvements to MOOSE user workflow through polyhedral elements, automation, and concise physics syntax

The MOOSE framework is a foundational capability used by the NEAMS program to create over 15 different simulation tools for advanced nuclear reactors. Due to MOOSE's broad use, improvements to the framework in support of modeling and simulation goals are critical to the program. Such improvements can take many forms, including optimization, improved user experience, streamlined application programming interfaces (APIs), parallelism, and new capabilities. The work described in this report was conducted in direct support of NEAMS tools and includes: addition of support for polyhedral elements, incorporation of mesh smoothers for mesh repair, integration of the Physics and ActionComponents systems, expansion of the Convergence system, and exploration of automated input file generation. These five areas of development are fundamental capabilities that will be leveraged by many NEAMS applications.

97 - MATHEMATICS AND COMPUTING↗

DOE OSTI2025

ITreeForeCast: An integrated modeling software to simulate tree level growth and forest carbon storage

Healthy trees in forest act as a natural carbon sink, capturing carbon. As they grow, they store carbon in their trunks, leaves and roots. Not all trees store carbon at the same rate, or in the same quantities, as it depends on a variety of biophysical and climatic factors. Furthermore, although carbon estimation in trees can be complex, the precision of estimates is tightly linked to trees growth, both in diameter and height. However, the simulation of carbon uptake by forest and forest growth has each been modeled separately, and independently at differing levels of detail and spatial resolution. In this paper, we introduce ITreeForeCast, a simulation model combining the two types of modeling on a unified platform, enabling the investigation of impacts of management strategies on carbon sequestration and wood products. ITreeForeCast is a user-extendable framework that offers new opportunities to model, simulate, and visualize the dynamics of individual trees in a forest, simulate management strategies over time, and carbon uptake.

09 - BIOMASS FUELS↗

DOE OSTI2026

A bioinspired approach for adaptive solid-solid phase change material coatings with optimized surface features for passive thermal regulation

The necessity to reduce global energy consumption calls for innovative strategies in building thermal management. Passive thermal regulation, particularly through bio-inspired designs, offers a promising avenue by mimicking nature's efficient control of optical properties. This research introduces a novel, climate-responsive coating that integrates optimized bio-inspired surface features with a solid-solid phase change material (SS-PCM) to dynamically manage solar absorptivity without adding additional thickness, enabling both heating and cooling as needed. Drawing on the photonic architectures of the Saharan silver ant and Morpho Didius butterfly, we employed a modeling and multi-objective optimization framework to tailor these surface features. Simulations reveal that surface texture, rather than the intrinsic phase transition of the SS PCM, dominates optical control. Relative to a flat SS PCM coating, optimized isotropic random roughness and broader range features yielded the highest passive heating power increase of about 144 % and 319 % respectively suitable for cold climates. Saharan ant-inspired features enhanced passive cooling for hot climates, achieving a 21.8 % improvement. For moderate climates, Butterfly-wing-inspired surface features provided a balanced enhancement of 19 % for heating and 7 % for cooling. Across all cases, the optimized surface features reduced combined heating and cooling energy demand more effectively than the baseline coating, while preserving material thickness. These findings demonstrate that climate-adaptive, optimized bio-inspired surface features can unlock the full potential of SS PCM coatings, providing a versatile pathway to significant energy savings in buildings and other applications. The methodology establishes a framework for designing next-generation adaptive envelopes that leverage natural photonic principles for high-impact, low-cost thermal regulation.

36 MATERIALS SCIENCE↗

DOE OSTI2026

Microstructure, Transport, and Mechanics of Compacted Clay Simulated at the 0.1 μm Scale (1400 Smectite Clay Particles) Using a Coarse-Grained Model with Explicit Counterions

Clay-rich geomaterials play a critical role in many subsurface systems. The macroscale properties of these materials (low permeability, high ionic conductivity, high swelling pressure, etc.) are sensitive to molecular-level adsorption and hydration interactions at clay−water interfaces. Efforts to develop multiscale simulation approaches to predict these properties reveal a scale gap between atomistic simulations (typically limited to systems smaller than 10 nm) and continuum-scale models (which use computational grid elements with dimensions ≳ 10 μm). In this study, we present a coarse-grained (CG) framework that partly bridges this gap by simulating compacted smectite clay assemblages with dimensions of 0.1 μm containing 1,400 clay particles across a range of dry densities (1,050 to 1,850 kg·m −3 ) and Na/Ca counterion compositions (Na fraction ranging from 0.2 to 1). The simulated systems, along with their reconstructed binary three-dimensional pore networks, are used to evaluate the microstructure, pore size distribution, tortuosity, ion diffusivity, and swelling pressure of compacted smectite clay. Results show that our approach captures important features of the mesoscale heterogeneity of compacted clays, including tactoid formation, hierarchical porosity, and anisotropic pore networks. Results also reveal how compaction and counterion composition govern emergent behaviors, including dominant pore sizes, directional transport, and electrochemical response. This work highlights the potential of CG simulations to bridge molecular and continuum scales and to advance geotechnical and environmental applications involving clay-rich materials as well as related nanoporous media such as geopolymers and calcium-silicate-hydrate. However, the results also suggest that accurate prediction of certain microstructural and mechanical properties (e.g., swelling pressure) may require even larger-scale systems on the order of 1 μm.

02 PETROLEUM↗

DOE OSTI2026

An Integrated Modeling Framework for Sediment Dynamics During Urban Flooding: Application to Hurricane Harvey in Houston

Floodwater can mobilize and redistribute large volumes of sediment from upland to downstream urban areas, threatening infrastructure, water quality, and ecosystem health. However, existing modeling approaches often fail to capture sediment dynamics in urban floodplains due to the lack of integration between upland hydrological processes and riverine sediment transport. This study presents the first integrated modeling framework that couples the Energy Exascale Earth System Model (E3SM) land component, which simulates runoff and hillslope erosion, with TELEMAC-GAIA, a two-dimensional hydrodynamic and sediment transport model. This framework enables the fully distributed, process-based simulation of high-resolution (as fine as 30 m) sediment dynamics from hillslopes to floodplains. Applied to a highly urbanized watershed in Houston during Hurricane Harvey, this framework reproduced observed water levels at 16 USGS gauges (median R 2 = 0.83 and KGE = 0.78), key sediment dynamics such as sediment transport and deposition processes, and reproduced spatial deposition patterns consistent with LiDAR-derived data. Based on the simulation, we estimate 8.0 million m 3 of event-scale sediment deposition, including 5.7 million m 3 trapped in the flood-control reservoirs and 2.3 million m 3 deposited along major channels and floodplains. Using a representative unit removal cost, this corresponds to an estimated dredging cost of $581 million for total deposition. These results provide a first-order, physically based quantification of Harvey-scale sediment impacts. This study provides a valuable tool for the holistic analysis of sediment dynamics triggered by extreme urban flooding, supporting flood-resilience planning. More broadly, it highlights the importance of integrating physically based hydrological processes for urban flooding and sediment research.

Hurricane Harvey↗

DOE OSTI2026

Determining Optimal Magnetometer Configuration on MAGIS-100

Long-baseline atom interferometers such as the Matter-wave Atomic Gradiometer Interferometric Sensor (MAGIS-100) require stringent control and continuous characterization of background magnetic fields and spatial gradients to prevent systemic phase shifts that mimic ultralight dark matter or gravitational wave signatures. Because direct sensor placement within the ultra-high vacuum beam pipe is infeasible, in-situ magnetic field monitoring relies on external sensor arrays situated in the surrounding annular region. This work demonstrates a field reconstruction framework for a 5.3-meter MAGIS-100 modular section using finite-element Opera simulations. Transverse magnetic fields are expanded using a cylindrical multipole framework as informed by Fermilab’s Muon g-2 experiment, with magnetometer array configurations optimized via Fisher information matrix D-optimality. Inverting external sensor readings through a Gauss-Newton scheme recovers interior tube fields across distinct axial positions. In the discontinuity-averse uniform region (slice pair P4), the model achieves sub-noise-floor performance with a cross-validated root-mean-square error (RMSE) of $6.7227 \times 10^{-4}\text{ A/m}$ ($0.845\times$ sensor noise floor) and an interior field coefficient of variation of $1.71\%$. An elbow criterion in the Fisher bounds establishes $n_{\text{max}} = 2$ as the optimal multipole truncation order to prevent noise amplification from over-parameterization, with $n_{\text{max}} = 3$ (sextupole) order chosen for analysis to demonstrate further complexity and cross-pair comparison. Furthermore, analytical differentiation of the fitted multipole coefficients yields dense spatial maps of the transverse Jacobian gradient matrix $\nabla \mathbf{H}$ along with propagated $1\sigma$ uncertainty bounds across the beam region ($r \le 2.75\text{ in}$). This operational framework confirms that external magnetometer arrays can reliably monitor magnetic field uniformity and spatial gradients along the 100-meter flight path given appropriate sampling for any complexity order.

Appleby, Darwin [William Rainey Harper Coll.] (ORC↗

DOE OSTI2026

Determining Optimal Magnetometer Configuration on MAGIS-100

Long-baseline atom interferometers such as the Matter-wave Atomic Gradiometer Interferometric Sensor (MAGIS-100) require stringent control and continuous characterization of background magnetic fields and spatial gradients to prevent systemic phase shifts that mimic ultralight dark matter or gravitational wave signatures. Because direct sensor placement within the ultra-high vacuum beam pipe is infeasible, in-situ magnetic field monitoring relies on external sensor arrays situated in the surrounding annular region. This work demonstrates a field reconstruction framework for a 5.3-meter MAGIS-100 modular section using finite-element Opera simulations. Transverse magnetic fields are expanded using a cylindrical multipole framework as informed by Fermilab’s Muon g-2 experiment, with magnetometer array configurations optimized via Fisher information matrix D-optimality. Inverting external sensor readings through a Gauss-Newton scheme recovers interior tube fields across distinct axial positions. In the discontinuity-averse uniform region (slice pair P4), the model achieves sub-noise-floor performance with a cross-validated root-mean-square error (RMSE) of $6.7227 \times 10^{-4}\text{ A/m}$ ($0.845\times$ sensor noise floor) and an interior field coefficient of variation of $1.71\%$. An elbow criterion in the Fisher bounds establishes $n_{\text{max}} = 2$ as the optimal multipole truncation order to prevent noise amplification from over-parameterization, with $n_{\text{max}} = 3$ (sextupole) order chosen for analysis to demonstrate further complexity and cross-pair comparison. Furthermore, analytical differentiation of the fitted multipole coefficients yields dense spatial maps of the transverse Jacobian gradient matrix $\nabla \mathbf{H}$ along with propagated $1\sigma$ uncertainty bounds across the beam region ($r \le 2.75\text{ in}$). This operational framework confirms that external magnetometer arrays can reliably monitor magnetic field uniformity and spatial gradients along the 100-meter flight path given appropriate sampling for any complexity order.

Appleby, Darwin [William Rainey Harper Coll.] (ORC↗

DOE OSTI2026

In-Situ Magnetic Field Reconstruction in the MAGIS-100 Experiment

Long-baseline atom interferometers such as the Matter-wave Atomic Gradiometer Interferometric Sensor (MAGIS-100) require stringent control and continuous characterization of background magnetic fields and spatial gradients to prevent systemic phase shifts that mimic ultralight dark matter or gravitational wave signatures. Because direct sensor placement within the ultra-high vacuum beam pipe is infeasible, in-situ magnetic field monitoring relies on external sensor arrays situated in the surrounding annular region. This work demonstrates a field reconstruction framework for a 5.3-meter MAGIS-100 modular section using finite-element Opera simulations. Transverse magnetic fields are expanded using a cylindrical multipole framework as informed by Fermilab’s Muon g-2 experiment, with magnetometer array configurations optimized via Fisher information matrix D-optimality. Inverting external sensor readings through a Gauss-Newton scheme recovers interior tube fields across distinct axial positions. In the discontinuity-averse uniform region (slice pair P4), the model achieves sub-noise-floor performance with a cross-validated root-mean-square error (RMSE) of $6.7227 \times 10^{-4}\text{ A/m}$ ($0.845\times$ sensor noise floor) and an interior field coefficient of variation of $1.71\%$. An elbow criterion in the Fisher bounds establishes $n_{\text{max}} = 2$ as the optimal multipole truncation order to prevent noise amplification from over-parameterization, with $n_{\text{max}} = 3$ (sextupole) order chosen for analysis to demonstrate further complexity and cross-pair comparison. Furthermore, analytical differentiation of the fitted multipole coefficients yields dense spatial maps of the transverse Jacobian gradient matrix $\nabla \mathbf{H}$ along with propagated $1\sigma$ uncertainty bounds across the beam region ($r \le 2.75\text{ in}$). This operational framework confirms that external magnetometer arrays can reliably monitor magnetic field uniformity and spatial gradients along the 100-meter flight path given appropriate sampling for any complexity order.

Appleby, Darwin [William Rainey Harper Coll.; Ferm↗

DOE OSTI2026

Time-dependent-bases with local CUR decomposition method for accelerating turbulent combustion simulations

Here, this study presents a novel reduced-order modeling framework, Time-Dependent Bases with Local CUR decomposition (TDB-L-CUR), designed to efficiently and accurately approximate the species transport equations in reacting flow simulations. The method extends the existing TDB-CUR approach for chemically reacting flows (Jung et al. Comput. Methods Appl. Mech. Engrg. 437 (2025) 117758), which leverages matrix decomposition techniques to form a global-in-space, time-dependent low-dimensional manifold. While TDB-CUR performs well in homogeneous systems, it may be less well-suited to spatially heterogeneous systems such as turbulent flames, where higher-rank approximations are typically required. The proposed TDB-L-CUR framework introduces two methodological extensions to the baseline approach. First, it applies unsupervised clustering to partition the physical domain into distinct regions, enabling spatially localized manifold construction, thereby reducing the rank required for the reduced-order representation. Second, it incorporates a computational singular perturbation (CSP)-based scheme for identifying and penalizing fast species, allowing for spatio-temporally adaptive mitigation of chemical stiffness. The proposed framework is validated on a hierarchy of test cases, including a one-dimensional premixed flame, a two-dimensional nonpremixed ignition case with vortex interaction, and a three-dimensional turbulent premixed flame. TDB-L-CUR significantly improves accuracy over TDB-CUR while further reducing computational cost. The fully on-the-fly formulation of TDB-L-CUR (i.e., requiring no offline training or prior knowledge) makes it a robust and scalable tool for reduced-order modeling of reactive flows.

Local manifold↗