SEARCH · Search NASA

Results for “ONE-DIMENSIONAL FLOW”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 91 records · Page 5

All sources NASA NTRS DOE OSTI

DOE OSTI2025

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↗

NASA NTRS1965

Attenuation of Sinusoidal Perturbations Superimposed on Laminar Flow of a Liquid in a Long Line

The attenuation constant for sinusoidal pressure and flow perturbations superimposed on the laminar flow of a viscous liquid was measured in a system consisting of a long, straight, cylindrical hydraulic line. The upstream and downstream ends of the line were securely fastened t o the ground. A sinusoidal perturbation was imposed on the mean flow at the upstream end by means of a s m a l l oscillation of a throttle valve abmt a partly open mean position. The downstream end was terminated in a restricting orifice. Pressure perturbations were measured at three locations along the line for frequencies from 15 t o 100 cps. These pressure measurements were reduced by use of a pair of complex damped acoustic one-dimensional wave equations to obtain the attenuation constant along with the phase constant and the dimensionless downstream admittance. For the range of frequencies investigated, the experimental values of the attenuation constant are in good agreement with classical theory.

PERTURBATION↗

DOE OSTI2026

Simulations of a hypersonic turbulent boundary layer over wavy surfaces

Here, we conduct large-eddy simulations of a Mach 5.84 cold wall turbulent boundary layer over one-dimensional wavy walls with varying amplitudes and wavelengths. Across all wall topologies, a series of alternating shock and expansion waves is shown to influence the entire boundary layer, and generate repeating wave patterns in the turbulent stresses, dispersive stresses, and turbulent kinetic energy budget. The series of alternating shocks and expansions imposes repeating adverse and favourable pressure gradients across the wavy wall, and at sufficient wall amplitude, triggers flow separation in the trough of the wave. Flow separation is demonstrated to influence the behaviour of wall pressure fluctuations over the wavy wall. In attached flows, the prominent frequencies are consistent with integral-scale boundary layer turbulence, whereas in separated flows, a two-decade frequency range is present, akin to two-dimensional shock–boundary layer interactions. Counter-rotating streamwise-oriented structures are observed on the windward side of the wave, which diminish over the wave crest. A conditional analysis demonstrates that these structures are present in the upstream boundary layer, and are amplified with increasing wall amplitude. An examination of the Görtler number and probability density function (PDF) of the fluctuating lateral wall shear stress demonstrates the strong correlation between a large Görtler number and growth of the PDF tail density, suggesting that the amplification of the counter-rotating streamwise-oriented structures are linked to centrifugal instabilities in regions of streamline concavity.

boundary layers↗

NASA NTRS1954

Thrust characteristics of a series of convergent-divergent exhaust nozzles at subsonic and supersonic flight speeds

An experimental investigation of a series of four convergent-divergent exhaust nozzles was conducted in the Lewis 8-by-6 foot supersonic wind tunnel at Mach numbers of 0.1, 0.6, 1.6, and 2.0 over a range of nozzle pressure ratios. The thrust characteristics of these nozzles were determined by a pressure-integration technique. From a thrust standpoint, a nozzle designed to give uniform parallel flow at the exit had no advantage over the simple geometric design with conical convergent and divergent sections. The rapid-divergent nozzles might be competitive with the more gradual-divergent nozzles since the relatively short length of these nozzles would be advantageous from a weight standpoint and might result in smaller thrust losses due to friction. The thrusts, with friction losses neglected, were predicted satisfactorily by one-dimensional theory for the nozzles with relatively gradual divergence. The thrusts of the rapid-divergent designs were several percentages below the theoretical values at the design pressure ratio or above, while at low pressure ratios there was a considerable effect of free-stream Mach number, with thrusts considerably above theoretical values at subsonic speeds and somewhat above theoretical values at supersonic speeds. This Mach numb effect appeared to be related to the variation of the model base pressure with free-stream Mach number.

Fradenburgh, Evan A↗

DOE OSTI2025

Origin of Pulsed Radio Emission from Magnetars

Extended periods of radio pulsations have been observed for six magnetars, displaying characteristics different from those of ordinary pulsars. In this Letter, we argue that radio emission is generated in a closed, twisted magnetic flux bundle originating near the magnetic pole and extending beyond 100 km from the magnetar. The electron–positron flow in the twisted bundle has to carry electric current and, at the same time, experiences a strong drag from the radiation field of the magnetar. This combination forces the plasma into a “radiatively locked” state with a sustained two-stream instability, generating radio emission. We demonstrate this mechanism using novel first-principles simulations that follow the plasma behavior by solving the relativistic Vlasov equation with the discontinuous Galerkin method. First, using one-dimensional simulations, we demonstrate how radiative drag induces the two-stream instability, sustaining turbulent electric fields. When extended to two dimensions, the system produces electromagnetic waves, including superluminal modes capable of escaping the magnetosphere. We measure their frequency and emitted power and incorporate the local simulation results into a global magnetospheric model. The model explains key features of the observed radio emission from magnetars: its appearance after an X-ray outburst, wide pulse profiles, luminosities ∼10 30 erg s −1 , and a broad range of frequencies extending up to ∼100 GHz.

Astronomical simulations↗

DOE OSTI2025

Vacuum-assisted carbon molecular sieve membrane reactor for non-oxidative ethane dehydrogenation

Non-oxidative ethane dehydrogenation (EDH) is equilibrium-limited and endothermic. Selective hydrogen removal using a gas-permeable membrane within the EDH reaction zone can overcome the thermodynamic equilibrium, enabling higher ethane conversions. Employing vacuum as the permeation driving force, rather than a sweep gas, enhances the industrial viability of membrane reactors by eliminating additional post-reaction separation units. This study presents a membrane reactor that integrates H 2 -permeable carbon molecular sieve (CMS) hollow fiber membranes embedded in a fixed bed of cobalt in a dealuminated beta zeolite (Co@DeAl-BEA) catalyst, utilizing a vacuum to remove hydrogen efficiently. The CMS membrane exhibits high hydrogen permeance and an excellent H 2 /C 2 H 6 separation factor. The membrane reactor significantly enhanced the ethane conversion under reaction conditions comparable to those reported in the literature. A Langmuir-Hinshelwood kinetic rate expression was developed and incorporated into a one-dimensional steady-state reactor model. The experimentally validated model indicates that increasing the number of hollow fibers improves ethane conversion, although ethane loss to the permeate limits the benefit. The contact area between the catalyst and the membrane limits the reactor performance more than the catalytic throughput. Furthermore, we find that the location of the catalyst packing relative to the hollow fiber membranes influences ethane loss and conversion. Higher reactor pressures and inlet ethane flow rates improve space-time yield at the expense of lower ethane conversion. Increasing reactor temperature or packing length promotes both performance metrics. The EDH membrane reactor demonstrated durability over 200 h of continuous operation, maintaining record-low deactivation rates and high ethylene selectivity. Protocols for catalyst regeneration were developed.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

DOE OSTI2026

Measurement of Local Impedance Characteristics of In-Plane PEM Electrolyzer Component Features via Segmented Cell

Segmented cells are a class of advanced diagnostic devices that enable the measurement of current distribution within the active area of an electrochemical cell, allowing for the evaluation of localized impacts from operating conditions, flow fields and component inhomogeneities. Segmented cell devices have been applied to support R&D efforts for the advancement of polymer electrolyte membrane fuel cell technologies and more recently proton exchange membrane water electrolysis (PEMWE). Beyond current density distributions, segmented cells have been successfully deployed to collect the distribution of other relevant parameters such as high frequency resistance and electrode potentials. Very recently, a promising impedance method has been reported for water electrolysis where a one-dimensional segmented cell is coupled to a multichannel potentiostat, to enable the measurement of local cell impedances, further expanding the diagnostic capabilities of the segmented cell. In this study, we will present recent efforts to accurately measure local impedance characteristics for a PEMWE cell using a two-dimensional segmented cell device coupled with a multichannel potentiostat. The distribution of local resistances, i.e. high frequency, charge transfer, mass transport and catalyst layer, and iR-free overpotentials were successfully measured after eliminating sources of random and systematic error. We will discuss the application of this diagnostic on a cell containing a 1 cm feature simulating an ionomer skin. The results from these tests give insights if such an irregularity needs to be classified as a defect, and highlight the usefulness of coupling impedance spectroscopy with segmented cells for holistic spatial diagnostics of electrochemical devices.

08 HYDROGEN↗

DOE OSTI2026

Insights on carbon dioxide adsorption in a flexible 1-D coordination polymer from in situ X-ray scattering and density functional theory

A combination of in situ small-angle X-ray scattering (SAXS) microstructure characterization over scales ranging from 1 nm to 10 µm and powder X-ray diffraction (XRD) structure characterization under various gas/pressure/temperature conditions with density functional theory (DFT) calculations provides new insights for the CO2 sorption behavior of a one-dimensional porous coordination polymer: catena-bis(dibenzoylmethanato)(4,4′-bipyridyl)nickel(II), denoted NiDBM-Bpy. The NiDBM-Bpy chains are held together by van der Waals forces, but the structure of guest-free NiDBM-Bpy is unsolved due to a lack of suitable crystals and high-quality powder XRD patterns. Nevertheless, SAXS and powder XRD can follow microstructural and structural changes as a function of gas pressure, composition and temperature. Both mixed-gas flow and static supercritical CO2 regimes are explored experimentally. DFT calculations are used to model the structural variation associated with XRD changes and hysteresis in the sorption isotherms. XRD and DFT calculations suggest that an orthorhombic Fddd structure with two CO2 per Ni emerges following a transition from the structure with lower molar volume and symmetry that exists without CO2 present.

Allen, Andrew J.↗

DOE OSTI2025

Modeling aerosol bolus inhalations in the human lung with the multiple path particle deposition model: Comparison with experimental data

Existing one-dimensional (1D) models of aerosol dosimetry often ignore mixing mechanisms of inhaled aerosols during their transport in the lung. This mixing or aerosol dispersion results from different physical mechanisms in different regions of the lung. It is a higher order effect, which cannot be directly captured in 1D modeling approaches, and thus is sometimes modeled as a diffusive process. Here, in this study, we improved our recently developed alveolar mixing module incorporated in the multiple path particle dosimetry model (MPPD) to account for flow irreversibility and particle trapping in the alveolar spaces, as well as mixing occurring in the tracheobronchial region. This new version of MPPD was coupled with CFPD-based predictions of aerosol bolus dispersion in the oral airway. The model was used to predict the deposition, dispersion, and mode shift of aerosol bolus inhaled at different penetration depths within the lung for breathing patterns and particle size matching those used in a previous experimental study (Darquenne et al., 2016). Even though a quite simplified approach was used, the computations appear to describe subject-specific and test-specific experimental data reasonably well. The proposed combined dispersion-deposition model can be a useful tool for targeted drug delivery and also for exposure health risk assessment.

MPPD↗

DOE OSTI2026

Uncertainty-Aware Machine Learning for Small-Angle X-ray Scattering Analysis in Autonomous Experimentation

Small-angle X-ray scattering (SAXS) is a powerful high-throughput characterization tool for probing nanoscale structure in native sample environments, providing real-time morphological information such as nanoparticle size and shape during synthesis. However, automated SAXS data analysis for extracting meaningful structural parameters is non-trivial and remains a bottleneck in closed-loop experimentation towards autonomous materials discovery, which demands fast, reliable, and uncertainty-aware data analysis. Here, we develop a machine-learning approach for automated SAXS analysis tailored to closed-loop nanoparticle synthesis. A Random Forest (RF) regression model is trained on 100,000 synthetic SAXS curves generated from polydisperse spherical nanoparticles with realistic background contributions. Using normalized one-dimensional SAXS intensity profiles as input, the RF model directly predicts nanoparticle radius, size polydispersity, and background parameters, while the ensemble standard deviation across trees provides built-in uncertainty quantification (UQ). On synthetic data, we show that combining fit-quality metrics (R 2 , MAE) with thresholds on prediction uncertainty reliably identifies accurate parameter estimates without access to ground truth. We then apply the trained model to 365 experimental SAXS profiles of citrate-reduced gold nanoparticles synthesized using an automated droplet-flow microreactor with in situ SAXS at a synchrotron beamline, classifying the results into high- and low-confidence subsets based on UQ metrics. Finally, we integrate RF-based SAXS analysis into a simulated closed-loop optimization campaign using Gaussian process Bayesian optimization to minimize nanoparticle polydispersity, benchmarking against conventional automated Levenberg–Marquardt fitting. The RF-guided campaign exhibits substantially faster convergence and lower relative opportunity cost (∼0.07 vs ∼0.3), demonstrating that uncertainty-aware machine-learning SAXS analysis significantly enhances the efficiency and robustness of autonomous nanomaterials synthesis workflows.

Bayesian optimization↗

DOE OSTI2026

A Data-Driven Method for Modeling Creep-Fatigue Stress- Strain Behavior Using Neural ODEs

In this paper, we introduce a data-driven machine learning approach for modeling one-dimensional stress–strain behavior under cyclic loading, utilizing experimental data from the nickel-based Alloy 617. The study employs uniaxial creep–fatigue test data acquired under various loading histories and compares two distinct neural network-based ODE models. The first model, known as the black-box model, comprehensively describes the strain–stress relationship using a Neural ODE equation. To interpret this black-box model, we apply the Sparse Identification of Nonlinear Dynamical Systems (SINDy) technique, transforming the black-box model into an equation-based model using symbolic regression. The second model, the Neural flow rule model, incorporates Hooke’s Law for the linear elastic component, with the nonlinear part characterized by a Neural ODE. Both models are trained with experimental data to accurately reflect the observed stress–strain behavior. We conduct a detailed comparison with the standard Chaboche model, which includes three back stresses. Our results demonstrate that the neural network-based ODE models precisely capture the experimental creep–fatigue mechanical behavior, exceeding the standard Chaboche model’s accuracy. Furthermore, an interpretable model derived from the black-box neural ODE model through symbolic regression achieves accuracy comparable to the Chaboche model, enhancing its interpretability. The results highlight the potential of neural network-based ODE models to depict complex creep–fatigue behavior, eliminating the necessity for experts to define a specific, material-focused model form.

creep-fatigue↗

DOE OSTI2026

Coupled heat and mass transfer in internally cooled silica gel-coated aluminum foam heat exchangers for dehumidification

Solid desiccant dehumidifiers offer significant potential to improve humidity control in buildings due to their ability to independently regulate latent (humidity) and sensible (temperature) loads. However, adsorption of water vapor by desiccants is an exothermic process; the resulting heat of adsorption increases both the desiccant bed temperature and the temperature of the air flowing over the bed, thereby limiting the performance of desiccant dehumidifiers. To address this problem, desiccant-coated heat exchanger prototypes were developed by coating silica gel onto aluminum foam substrates and integrated with chilled water circulation for effective removal of the heat of adsorption. Prototype heat exchangers with two different pore densities (10 and 20 pores per inch (PPI)) and two different desiccant loadings (∼10.36% and ∼17.80%) were fabricated. Experiments were conducted by varying the inlet air face velocity, chilled water inlet temperature, and chilled water flow rate to investigate the coupled heat and mass transfer characteristics of the prototype heat exchangers over a wide range of operating conditions. The underlying transport mechanisms were analyzed using transient desiccant bed temperature, outlet air temperature, and outlet humidity profiles. The results demonstrate that internal cooling significantly enhances dehumidification performance by maintaining more favorable thermodynamic conditions for water vapor adsorption. Under identical operating conditions, a 20 PPI heat exchanger coated with 25.05 g of silica gel achieved 98.98% higher cumulative water uptake and 97.75% higher moisture removal capacity than the corresponding non-cooled case after 5 min of adsorption. This substantial enhancement was further supported by a validated one-dimensional coupled heat and mass transfer model developed in the present study. To date, this is one of the first experimental studies to systematically investigate the combined effects of internal cooling, foam pore density, and desiccant loading on the performance of silica-gel-coated metal foam heat exchangers under HVAC-relevant operating conditions.

42 ENGINEERING↗

DOE OSTI2025

Understanding the Effects of Inhomogeneities at the Back Interface of CdTe‐Based Solar Cells Using 2D Modeling

One-dimensional modeling cannot capture lateral inhomogeneities in CdTe-based devices. Here, we use 2D modeling to investigate the role of varying energetics at the back interface. We consider improvements in the back interface layer (BIL) through either reducing back surface recombination velocity (BSRV) or decreasing the downward band bending near the back interface. We show that when the BSRV is reduced, but strong downward band bending remains, there is no change in the device performance until the BSRV of 90% of the back interface is improved by the BIL. On the other hand, any coverage with a BIL that improves band bending results in device improvements. We use band bending, back interface recombination current densities, and voltage dependent current flow through the device to understand these improvements. The modeling shows that lateral flow of carriers greatly affects device performance, which is not captured in parallel diode modeling, and demonstrates improved understanding with 2D modeling.

2D modelling↗

DOE OSTI2025

Development, Verification, and Validation of an OpenFOAM-Based Solver for Modeling Inertial Fusion Energy Chambers

Our work seeks to introduce a computational tool tailored to the physics of inertial fusion energy chambers, in particular, those concepts based on thick liquid walls. In this approach, the structural materials are protected by several neutron mean-free-paths of renewable liquid and thus will be able to survive much longer than un-shielded walls, with virtually all structures lasting for the life of the plant and enabling the use of commercially available and qualified materials. The OpenFOAM-based solver named rhoCentralFoam has been used as a starting point. rhoCentralFoam belongs to the standard OpenFOAM solver toolset. It is a high-speed, explicit compressible flow solver with shock-capturing capability. While the main features have been retained, the solver had to be restructured to make use of tabular data for equations of states, a necessary addition to model the complex thermo-physical properties of ionized gasses. This entailed the need to change the independent state variables used by the solver, resulting in a new thermodynamic library and slightly different solution algorithm. Moreover, a radiation heat transfer model based on the P-1 approximation was added to the solver. The solver is verified against an analytical solution from the Sedov-Taylor-Neumann test problem to showcase the ability of the hydrodynamic solvers to handle strong shocks, whereas the P-1 model was verified using a simple one-dimensional problem with an analytical solution. Additionally, a validation case involving shock-wave propagation through jet array is presented, and the results are compared with experimental data from the open literature. Lastly, in order to showcase the utility of the solver for practical cases, we applied the refined solver to two representative scenarios: gas venting within the HYLIFE-II chamber and the compression of the gas following the partial ablation of the liquid wall.

Chamber dynamics↗

DOE OSTI2025

A real-time distributed solid oxide electrolysis cell (SOEC) model for cyber-physical simulation

System integration and dynamic operability between SOEC and balance-of-plant (BoP) components are major technical challenges before realizing rapid load following of SOEC systems. Cyber-physical simulation (CPS) is a leading-edge digital engineering approach and is regarded as the next step beyond Digital Twins. CPS approach can be used to research SOEC system integration and develop dynamic controls prior to actual pilot testing without using a real SOEC. To seamlessly couple with BoP hardware and access non-observable operational parameters (e.g., local temperature gradient) during transients, a distributed one-dimensional (1D) real-time SOEC model was developed. Its real-time execution was demonstrated for 20 to 640 nodes at the fixed time step of 5 ms. A higher excess air ratio enabled smaller local temperature gradients on SOEC solid materials and faster transients upon current density step change from 0.15 to 0.55 A cm -2 . During the transients, the magnitude of the peak temperature gradient nearly doubled in 10 s from -3.5 to -5.9 °C cm -1 . This represents a significant operating risk that can impact the dynamic operability of SOEC systems. In addition, the local temperature gradient was found to change directions on all nodes in SOEC solid materials, with the greatest impact on the upstream nodes. The SOEC model was also tested at the thermal neutral voltage using actual process air flow parameters as variable model inputs. Variable process air temperatures were found to induce alternating local temperature gradients on SOEC solid materials. These are new operational mechanisms for SOEC degradation relevant for load following operational modes yet distinct from previous reports. To mitigate these unfavorable features, the SOEC can be operated at voltages that are slightly (±20 mV) deviated from the thermal neutral voltage. Here, the corresponding net thermal energy change was less than 1.6% of the electric power consumption. This 1D real-time SOEC model established the basis of cyber-physical simulation of SOEC hybrid systems.

24 POWER TRANSMISSION AND DISTRIBUTION↗

DOE OSTI2025

Generative AI models for learning flow maps of stochastic dynamical systems in bounded domains

Simulating stochastic differential equations (SDEs) in bounded domains, presents significant computational challenges due to particle exit phenomena, which requires accurate modeling of interior stochastic dynamics and boundary interactions. Despite the success of machine learning-based methods in learning SDEs, existing learning methods are not applicable to SDEs in bounded domains because they cannot accurately capture the particle exit dynamics. We present a unified hybrid data-driven approach that combines a conditional diffusion model with an exit prediction neural network to capture both interior stochastic dynamics and boundary exit phenomena. Our ML model consists of two major components: a neural network that learns exit probabilities using binary cross-entropy loss with rigorous convergence guarantees, and a training-free diffusion model that generates state transitions for non-exiting particles using closed-form score functions. The two components are integrated through a probabilistic sampling algorithm that determines particle exit at each time step and generates appropriate state transitions. Here, the performance of the proposed approach is demonstrated via three test cases: a one-dimensional simplified problem for theoretical verification, a two-dimensional advection-diffusion problem in a bounded domain, and a three-dimensional problem of interest to magnetically confined fusion plasmas.

Bounded domains↗

DOE OSTI2025

Unraveling membrane electrode assembly design for electrochemical conversion of carbon dioxide to formate/formic acid

This work presents a one-dimensional continuum modeling approach to investigate various cell architectures used for electrochemical conversion of CO 2 to formate/formic acid. Ion transport is simulated by a system of generalized modified Poisson–Nernst–Planck (GMPNP) equations that reflect the reactive transport phenomena including steric effects as the electrolyte solutions become concentrated. In the cathode catalyst layer, ionic current contributions from both the supporting electrolyte and solid-state ionomer are considered. Voltage and CO 2 utilization breakdowns are utilized to deconvolute the impacts of the cell architecture. The origins of (bi)carbonate formation in the cathode are explored, as the subsequent decrease in CO 2 availability is a key reason for low faradaic efficiencies to formate/formic acid. In addition, the role of a supporting electrolyte (KOH) is investigated to understand its tradeoffs: while the K + ions can improve both conductivity and electrochemically active surface area in the cathode, the presence of OH − ions raises the pH and leads to deleterious formation of (bi)carbonates. To this end, we also present parametric studies on the concentration and flow rate of supplied KOH to the cell, to establish a path towards eliminating the need for a supporting electrolyte.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

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↗