Search NASASearch

SEARCH · Search NASA

Results for “vortex 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.

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

At least 37 records · Page 2

Multimode turbulent flow measurements using magnetic resonance imaging- and laser-based techniques and computational fluid dynamics simulations

We studied the flow field characteristics of a turbulent flow over a regularized cube array with a perpendicular injection flow through the floor between the second and third cubical elements, representing the complex flow interactions of a 3D jet and the wake flows behind cubical obstacles. Four different experimental measurements were performed: two magnetic resonance imaging-based measurements for three-dimensional three-component velocity (MRV) and concentration (MRC) and two laser-based techniques, particle image velocimetry (PIV) and planar laser-induced fluorescence (PLIF), for two-dimensional two-component velocity and concentration measurement, respectively. The mainstream Reynolds number is Re = 15 000⁠, based on the primary inlet velocity U m and channel height D H ⁠, whereas the injector Reynolds number is Re j = 3400⁠, based on the injector velocity U j and the injector's exit width D j ⁠. Numerical simulations were performed for the studied flow configuration of turbulent flow over a regularized cube array using Reynolds-averaged Navier–Stokes (RANS) and large-eddy simulation (LES) approaches. Results obtained from experimental measurements—including MRV, MRC, PIV, and PLIF—as well as RANS and LES simulations are discussed and compared along several horizontal and vertical planes of the studied configuration. In addition, 3D turbulent flow structures, such as leading-edge vortex, horseshoe vortex, and jet shear-layer vortex, and the isosurfaces of scalar concentration successfully revealed by MRV and MRC techniques were found to be in very good agreement with those 3D features extracted from RANS and LES simulations. In conclusion, the high-resolution experimental and numerical database obtained from this study could be useful for validation and verification of numerical codes.

Computational fluid dynamics

VorLap

SAND2025-10210O VorLap is a vortex-induced vibration overlap prediction tool for static structures, such as parked wind turbines and marine turbines encountering fluid flow (water or air movement around them) that may induce vortex-induced vibration. This tool uses precomputed frequency domain data for specific cross-sectional shapes or a generalized shedding model, along with geometric data and structural natural frequency data, to identify areas and conditions where vortex-induced vibration may occur. Sandia National Laboratories is a multimission laboratory managed and operated by National Technology & Engineering Solutions of Sandia, LLC, a wholly owned subsidiary of Honeywell International Inc., for the U.S. Department of Energy’s National Nuclear Security Administration under contract DE-NA0003525.

Moore, Kevin [Sandia National Lab. (SNL-CA), Liver

Using principal component analysis to distinguish different dynamic phases in superconducting vortex matter

Vortices in type-II superconductors driven over random disorder are known to exhibit a remarkable variety of distinct nonequilibrium dynamical phases that arise owing to the competition between vortex-vortex interactions, the quenched disorder, and the drive. These include pinned states, elastic flows, plastic or disordered flows, and dynamically reordered moving crystal or moving smectic states. The plastic flow phases can be particularly difficult to characterize since the flows are strongly disordered. Here, we perform principal component analysis (PCA) on the positions and velocities of vortex matter moving over random disorder for different disorder strengths and drives. We find that PCA can distinguish the known dynamic phases as well as or better than previous measures based on transport signatures or topological defect densities. In addition, PCA recognizes distinct plastic flow regimes, a slowly changing channel flow and a moving amorphous fluid flow, that do not produce distinct signatures in the standard measurements. In conclusion, our results suggest that this position and velocity-based PCA approach could be used to characterize dynamic phases in a broader class of systems that exhibit depinning and nonequilibrium phase transitions.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND

Staircases of passive and active scalar concentration in cellular flow

This paper develops a unified model for staircase formation in both passive and active scalar systems, building upon prior numerical studies by offering new heuristic and physical insights. While prior studies primarily reported numerical results, they did not explore the underlying unifying physics that governs both types of scalar transport; this work addresses that gap by identifying shared mechanisms across both cases. Results of studies of passive and active scalar staircase formation in cellular flows are presented. Staircase formation in cellular flows occurs due to the interplay of fast mixing within cells and slow transport across the inter-cell boundary. The cell boundary emerges as a de facto transport barrier. Special attention is focused on the effects of cellular fluctuations and noise upon staircase structure. A forced, fluctuating vortex array model is used to drive the underlying flow structure. Cellular Peclet number and staircase profile curvature are identified as figures-of-merit to quantify the resiliency of layering. These are related to simple, multi-scatterer scalar random walk models. Results for Peclet number and curvature scaling with flow excitation are presented. We also study staircases of magnetic potential evolving in two-dimensional magnetohydrodynamics as examples of layering of active scalar concentration. Formation of magnetic potential staircases is indeed observed. Flux expulsion inhibits the intercellular transport of magnetic potential and strengthens staircase barriers. Magnetic staircases can be supported against resistive decay by magnetic potential noise forcing. Implications for staircase formation in magnetic confinement experiments are discussed.

Control theory

Self-Focused Pulse Propagation Is Mediated by Spatiotemporal Optical Vortices

We show that the dynamics of high-intensity laser pulses undergoing self-focused propagation in a nonlinear medium can be understood in terms of the topological constraints imposed by the formation and evolution of spatiotemporal optical vortices (STOVs). STOVs are born from point phase defects on the sides of the pulse nucleated by spatiotemporal phase shear. These defects grow into closed loops of spatiotemporal vorticity that initially exclude the pulse propagation axis, but then reconnect to form a pair of toroidal vortex rings that wrap around it. STOVs constrain the intrapulse flow of electromagnetic energy, controlling the focusing-defocusing cycles and pulse splitting inherent to nonlinear pulse propagation. Here, we illustrate this in two widely studied but very different regimes, relativistic self-focusing in plasma and nonrelativistic self-focusing in gas, demonstrating that STOVs mediate nonlinear propagation irrespective of the detailed physics.

42 ENGINEERING

Numerical Modeling of a Two-Stage Ocean Current Turbine

The Equinox Ocean Turbines (EQOT) current energy converter has a unique design with power generation in two small-diameter turbines attached to the tips of a large-diameter passive rotor. This configuration offers some key advantages for capturing ocean currents. With no centrally placed generator, almost no reaction torque is required at the nacelle of the main large-diameter rotor, and the small-diameter tip turbine generators operate at a higher speed and lower torque. The physics that determine the performance and loads on the turbine are also unique. The interactions of the flow field between the two stages and the general architecture of the system cannot be captured with traditional mid-fidelity modeling tools. For design iterations and large sets of load cases, it is important to have mid-fidelity models that can capture the important phenomenon with enough accuracy to identify global trends. This work uses a limited set of high-fidelity computational fluid dynamics (CFD) simulations to help inform the selection of and construction of a custom mid-fidelity model. Mid-fidelity modeling approaches were verified by comparing key turbine performance quantities to those found with the CFD model. Hydrodynamic interactions of the two-stage rotor were identified through high-fidelity CFD modeling. This highlighted the impact of the main rotor tip vortex and wake on the secondary rotor apparent inflow. This results in a relative flow rotation and sharp deficit, that change the optimal secondary rotor rotation speed and adds unsteadiness to the blade loading respectively. Multiple mid-fidelity approaches were evaluated for their ability to capture these effects. A simple approximation of the combined-stage performance based on single-stage BEM provides a reasonable rough prediction, especially near the peak TSR values, with some larger discrepancy at higher TSRs. Predicting the combined-stage performance based on single-stage CFD data improves this prediction across the TSR range. Although the combined-stage modeling in OLAF was not successful in this stage of the project, it showed promise as a mid-fidelity method, assuming the parameters can be tuned to account for the significant differences in time and length scales between the main and secondary rotors. This may be addressed through code changes in future work. A significant finding from the OLAF work was the agreement between the vortex core radius values found independently via a parameter space search and via CFD. The technique of using single-stage secondary rotor BEM, with a custom inflow taken from single-stage main rotor CFD or OLAF, provides an efficient method to capture one-way coupled flow interactions. This method provided generally good predictions of the impact of the flow rotation on the secondary rotor but struggled to accurately predict the peaks of the unsteady load progression. Future work could include some superposition of a tuned main rotor trailing edge viscous wake into the custom inflow to better predict this interaction.

16 TIDAL AND WAVE POWER

Growth of organized flow coherent motions within a single-stream shear layer: 4D-PTV measurements

Abstract This study investigates the evolution of a single-stream shear layer (SSSL) originating from a wall boundary layer past a backward-facing step. Utilizing a time-resolved 3D-Particle Tracking Velocimetry (4D-PTV) technique, we track the trajectories of fluorescent particles to gain insight into the flow characteristics of the SSSL. A compact water tunnel facility ( $$\textrm{Re}_\tau =1\,240$$ Re τ = 1 240 ) is fabricated to obtain an SSSL with a perpendicular slow entrainment stream past the separation edge. A hybrid interpolation approach that combines ensemble binning and Gaussian weighting is implemented to derive minimally filtered mean and instantaneous lower- and higher-order flow field parameters. Spanwise-dominant coherent motion accompanied by finer flow scales is observed to grow due to flow entrainment through “nibbling” actions of small-scale vortices, “engulfing” by large-scale vortices, and vortex pairing events. Furthermore, the non-zero-speed stream edge grows relatively faster than the zero-speed stream edge, showing a strong asymmetry in mixing composition across a mixing layer. The SSSL reaches self-similarity at a streamwise distance of $$\approx 55\,\theta _{0}$$ ≈ 55 θ 0 , where $$\theta _0$$ θ 0 is the initial momentum thickness from the separation edge, i.e., considerably shorter than reported in previous studies. A literature comparison of growth rate parameters raises intriguing questions regarding a potential inclusive growth scaling unifying the free shear layers. A turbulent kinetic energy (TKE) budget analysis reveals a negative production region immediately downstream of the separation edge attributed to a large positive streamwise gradient of streamwise velocity. In the self-similar region, the phase-averaged flow mapping demonstrates a larger concentration of turbulence production rate around the outer edges of spanwise vortices, specifically at the intersection of braids and vortices. Furthermore, a spatial separation exists in the regions of peak production and dissipation rates within the vortex core region favoring dissipation. The braids exhibit a larger concentration of turbulence diffusion rates, indicating their function as a conduit for exchanging turbulence between neighboring coherent motions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS

Optimizing the Heisenberg Vortex Tube for Hydrogen Cooling (Final Technical Report)

Hydrogen utilization at Plug Power sites ranges between 70-93% with the remainder vented and lost to the atmosphere. The goal of this project was to improve hydrogen utilization at Plug Power fulfillment centers via the patented Heisenberg vortex tube (HVT). The HVT combines conventional counterflow vortex tubes with para-orthohydrogen conversion catalyst to allow the cooling of hydrogen flows with no moving parts or external energy inputs. In year 1 of this project the HVT was evaluated to improve hydrogen utilization via the following concepts: 1) increasing liquid hydrogen pump volumetric efficiency by 20% through vapor separation and subcooling of the liquid, 2) reducing liquid hydrogen storage tank boil-off losses by 20% through thermal vapor shielding (TVS), and 3) increase isentropic efficiency of supercritical hydrogen expansion at 40-50 K, from 31% to over 40%. The TVS concept was selected as the most favorable application for further development. While the conceptual analysis was underway, the cryocatalysis hydrogen experiment facility (CHEF) was retrofitted with a new cryocooler, higher pressure condenser tanks, and in-situ fiber-optic Raman probes for ortho-parahydrogen composition analysis at the inlets and outlets of the HVT. To our knowledge this is the first in-situ implementation of cryogenic Raman probes for ortho-parahydrogen analysis. Subsequent testing of catalyzed and non-catalyzed HVT identified a low flow rate, high-conversion efficacy regime well suited for TVS development. Computational fluid dynamics (CFD), Reduced order modeling (ROM), and quantum Monte Carlo atomistic simulations were applied to optimize the design of the HVT implemented for TVS. The simulations and analysis identified a ruthenium-based catalyst as the most optimal for high conversion with little pressure losses and was matched to experimental measurements. These results indicated the design could achieve the 20% boil-off reduction target. An HVT field-trial was designed and constructed by Plug Power and implemented at a customer site. However, budget limitations reduced the amount of catalyst that could be applied to the HVT. With the reduced amount of catalyst, CFD analysis estimated about 3% reduction in boil-off. A similar amount of boil-off reduction was observed during preliminary measurements during tank commissioning. However, the results have a large margin for error and were operated at relatively low tank temperatures when the HVT has the lowest potential for improvement. This project demonstrated the use of para-orthohydrogen catalysis for reducing liquid hydrogen tank boil-off losses. The field trial system is anticipated to provide long-term experimental measurements on tank operational performance that will allow Plug Power to instrument additional tanks in the future. Recommendations for future work include the development of higher activity, lower cost para-orthohydrogen conversion catalyst for use in the HVT. Tank process optimizations could maximize the efficacy of the HVT and further improve hydrogen utilization.

08 HYDROGEN

Effect of Surface Roughness on Dynamic Stall in Pitching Motion

Dynamic stall plays a critical role in determining the performance and stability of a wide range of fluid-dynamics systems in various engineering applications. This unsteady aerodynamic phenomenon is particularly significant for maneuvering aircraft wings, jet aircraft subjected to gust encounters, helicopter rotor blades, and wind turbine blades [1–3]. The prediction of the dynamic stall vortex (DSV) is challenging due to factors such as unsteady aerodynamics, three-dimensional (3-D) effects, turbulence and flow separation, incoming gust, and surface impact effects [4–6]. Hence, advanced computational techniques and modeling approaches in computational fluid dynamics (CFD) would be required to enhance the accuracy and reliability of DSV predictions in dynamic stall scenarios. In the past, Batther and Lee [7] employed delayed detached eddy simulations (DDES) to understand the flow physics associated with the onset of dynamic stall. Their approach demonstrated that DDES achieves results comparable to those obtained from large-eddy simulations at a reduced computational cost. In another study, Khalifa et al. [8] examined the 3-D aspects of dynamic stall on a NACA 0012 airfoil using DES solvers. In conclusion, the findings underlined the superiority of 3-D simulations over two-dimensional approaches, particularly in predicting the lift coefficient values and capturing dynamic stall stages more precisely.

97 MATHEMATICS AND COMPUTING

Thermovelocimetric characterization of liquid metal convection in a rotating slender cylinder

Rotating turbulent convection occurs ubiquitously in natural convective systems encompassing planetary cores, oceans and atmospheres, as well as in many industrial applications. While the global heat and mass transfer of water-like rotating Rayleigh–Bénard convection is well-documented, the dynamics in low-Prandtl-number liquid metals remain less understood. In this study, we experimentally investigate rotating Rayleigh-Bénard convection in liquid gallium (Prandtl number Pr ≈ 0.027) within a slender cylinder (diameter-to-height aspect ratio Γ = D/H =1/2) using novel thermovelocimetric diagnostic techniques that integrate simultaneous multi-point thermometry and ultrasonic Doppler velocity measurements. Our results reveal the formation of a stable, global-scale azimuthal wavenumber m = 2 quadrupolar vortex at low supercriticality. We propose that enhanced wall modes facilitated by the slender cylinder geometry interact with the bulk flow to create these large-scale axialized vortices. Furthermore, our findings imply a distinct scaling behavior for the wall-mode precession frequency in liquid metals, extending previous results obtained for moderate-Pr fluids. In conclusion, this provides new insights into wall-bulk coupling mechanisms of low-Pr rotating convective turbulence.

Magnetohydrodynamics

Experimental and theoretical evidence of universality in superfluid vortex reconnections

The minimum separation between reconnecting vortices in fluids and superfluids obeys a universal scaling law with respect to time. The prereconnection and the postreconnection prefactors of this scaling law are different, a property related to irreversibility and to energy transfer and dissipation mechanisms. In the present work, we determine the temperature dependence of these prefactors in superfluid helium from experiments and a numeric model which fully accounts for the coupled dynamics of the superfluid vortex lines and the thermal normal fluid component. At all temperatures, we observe a pre- and postreconnection asymmetry similar to that observed in other superfluids and in classical viscous fluids, indicating that vortex reconnections display a universal behavior independent of the small-scale regularizing dynamics. We also numerically show that each vortex reconnection event represents a sudden injection of energy in the normal fluid. Finally we argue that in a turbulent flow, these punctuated energy injections can sustain the normal fluid in a perturbed state, provided that the density of superfluid vortices is large enough.

Reconnections

Three-Dimensional Aerodynamics and Vortex-Shedding Characteristics of Wind Turbine Airfoils over 360-Degree Angles of Attack

In this work, we present the first three-dimensional (3D) computational investigation of wind turbine airfoils over 360° angles of attack to predict unsteady aerodynamic loads and vortex-shedding characteristics. To this end, static–airfoil simulations are performed for the FFA-W3 airfoil family at a Reynolds number of 107 with the Improved Delayed Detached Eddy Simulation turbulence model. Aerodynamic forces reveal that the onset of boundary-layer instabilities and flow separation does not necessarily coincide with the onset of stall. In addition, a comparison with two-dimensional simulation data and flat plate theory extension of airfoil polars, suggest that, in the deep stall regime, 3D effects remain critical for predicting both the unsteady loads and the vortex-shedding dynamics. For all airfoils, the vortex-shedding frequencies are found to be inversely proportional to the wake width. In the case of slender airfoils, the frequencies are nearly independent of the airfoil thickness, and their corresponding Strouhal number St is approximately 0.15. Based on the calculated St, the potential for shedding frequencies to coincide with the natural frequencies of the International Energy Agency 15 MW reference wind turbine blades is investigated. The analysis shows that vortex-induced vibrations occur primarily at angles of attack of around ±90° for all airfoils.

17 WIND ENERGY

Freestream turbulence effects on unsteady wind turbine loads and wakes: An IDDES study

We investigate numerically the effects of freestream turbulence on the unsteady aerodynamics and wakes of the National Renewable Energy Laboratory Phase VI wind turbine rotor for increasing wind speed. Turbulence is modeled using the Improved Delayed Detached-Eddy Simulation (IDDES) method. As a first step, a detailed mesh resolution study is conducted with the decaying freestream turbulence model at turbulence intensity of 0.5%. Our blade-resolved IDDES simulations show that grid-independent average torque and thrust results can be achieved with relatively coarse meshes, whereas dramatically higher mesh resolution is required for grid-independent results for power spectral densities of thrust force, especially in the deep-stall regime. Comparing the loads with the Shear-Stress Transport model demonstrates the superiority of IDDES in predicting massively separated flows. The aerodynamic performance and wake predictions with the decaying freestream turbulence model are compared with the synthetic freestream turbulence model. Both models predict nearly the same loads, spectral energy content, and wake characteristics. The properties of both the near- and far-wake regions are then examined. Furthermore, we show that separated boundary layers accelerate turbulent mixing and entrainment of the external flow, which results in faster wake recovery. The effect of increasing turbulence intensity to 6% is investigated using the synthetic freestream turbulence model. In contrast with the fully attached boundary layer, higher freestream turbulence in deep stall does not significantly affect the loads and vortex-shedding characteristics. However, the turbulent mixing in the wake is enhanced, which further hastens the recovery of the self-similar velocity profile. In general, increasing the wind speed at high turbulence intensity shifts the recovery farther upstream and increases the wake width.

17 WIND ENERGY

Modeling Flow in a Microreactor Core: From High-Fidelity CFD to Subchannel Analysis

The accurate prediction of pressure drop in tightly packed, low–Reynolds number (Re) bare rod bundles is essential for the thermal-hydraulic design of the Microreactor Applications Research Validation and EvaLuation (MARVEL) reactor and other microreactor concepts. However, existing friction factor correlations, particularly the upgraded Cheng-Todreas (UCTD) correlation, have limited validation for the small pitch-to-diameter ratios (P/Ds) and transitional flow conditions characteristic of these systems. In this work, we perform high-fidelity large-eddy simulations (LESs) of both an infinite bare rod bundle and a finite bare 37-pin scalloped rod bundle across the range 1000≤ 𝑅𝑒 ≤5500. The simulations reveal strong gap vortex–driven transitional behavior and indicate that the UCTD may underpredict the friction factor by up to 28% at 𝑃/𝐷 = 1.05. Using the LES-calculated pressure drops, we formulate a new friction factor correlation that follows the Cheng-Todreas functional form but is calibrated for low-Re and tightly packed geometries representative of a MARVEL-like reactor. The correlation is implemented in the MOOSE (Multiphysics Object-Oriented Simulation Environment) subchannel module and compared against both the LES and UCTD predictions. Across all subchannel types, the proposed model reduces the streamwise velocity differences from as high as 44% (UCTD) to below 9%, and decreases the pressure gradient differences from 13% to 25% (UCTD) to 0.7% to 7% relative to the LES results. These results suggest that the new correlation has the potential to improve the pressure drop and flow field predictions for such geometries, highlighting the importance of high-fidelity simulations in supporting microreactor thermal-hydraulic model development and motivating future pressure drop experiments for compact rod bundles to further validate these findings.

21 - SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLAN

Microwave Vortex Motion Characterization of Nb$_{3}$ Sn Coatings for Applications in High Magnetic Fields

In this work, microwave measurements carried out in dielectric-loaded resonators exposed to high magnetic fields are exploited to yield the surface impedance of Nb$_3$ Sn superconducting coatings deposited via two different techniques: vapor tin diffusion, and DC magnetron sputtering. The obtained data lead to qualitative interpretations on both the Nb$_3$ Sn superconducting properties, and vortex-dynamics and pinning, of each coating separately, as well as simple distinctive features when comparing those. When examining the respective surface impedances at varying field, it is expected that the studied films perform at substantially diverse magnitudes of flux-flow resistivity, but also in well-differentiated pinning regimes, yet the obtained surface resistances of both samples are comparable, thus demonstrating that there is room for film optimization at the expense of certain compromise between the parameters involved.

Nb3Sn

Harvesting Subsea Water Motion to Provide Clean Electrical Power

This TEAMER project had three main objectives related to characterizing the performance of a vortex-induced vibration (VIV) marine energy harvesting device called the Whatever Input To Torsion (WITT). The first was to characterize the motion of the WITT device as it vibrates on top of a pipe in a natural tidal flow. The second was to measure how much electricity the WITT device generates when it is deployed on top of a pipe that is secured to the seabed with a lander. The third objective was to measure how much electricity the WITT device generates in the same configuration but with a longer pipe. This TEAMER project succeeded in satisfying the first objective by measuring the motion of the WITT device with an inertial measurement unit (IMU) mounted in the WITT housing. The motion was described by examining changes in the acceleration, pitch, roll, and yaw, and displacement. The frequency of the motion was characterized using spectral analysis. The TEAMER project succeeded in reaching the second objective by measuring the power produced by the WITT. It is reported here as the average and maximum power in hourly intervals. WITT Energy decided to not have us pursue the third objective of testing a longer pipe because the first test was unexpectedly successful at having the system vibrate at the target frequency. A longer pipe would have resulted in a decreased frequency and therefore not have been a useful test. Based on these findings, the decision was made to keep the same pipe length for both tests and improve on other aspects of the test structure that had failed during the first test. The key finding of this TEAMER project is that the WITT is capable of generating electricity in moderate flow speeds when it is in a housing that is attached to a 2.5 m pipe and secured to the seabed by a lander. This report contains the first power values for a VIV instrument that generates electricity from a pendulum swinging at the top of a long pipe. The experiment demonstrated that the vibration frequency of the device varied with the tidal current speed, which changed throughout a given tidal cycle. This meant that the desired peak frequency was not sustained for more than one hour. The success of this TEAMER test provides evidence that continued research in the area of electricity generation from VIV holds promise as a source of marine energy for offshore applications. Future tests should focus on refining the pipe length equations and improving the power take off system of the WITT power electronics. In addition, experiments aimed at expanded testing and evaluation (e.g. longer than one month) are needed to measure the power production at higher tidal current speeds and test the system's durability.

Branch, Ruth A. (ORCID:0000000202356719)

Separating Sensible and Latent Cooling with Electrically Charged Rotating Vortexes and Vapor Capturing Air Handler Technology

The dehumidification of air is a highly energy-intensive process, typically requiring vapor compression systems to cool the air below its dew point temperature. This project investigated air dehumidification achieved by electrostatic and dielectrophoresis forces in electrosprays and centrifugal principles in cyclone separators. As a practical implementation of these processes, Electrospray Vortical Flow eXchanger (EVFX) prototypes were designed, fabricated, and tested. Electrosprays were employed to energize water droplets, which then captured water vapor via dielectrophoresis. A cyclone separator was positioned downstream of the electrosprays. This combination increased the dehumidification efficiency and controlled the flow of the water vapor condensate removed from the air.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI

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