SEARCH · Search NASA

Results for “Fluid Analysis”

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 361 records · Page 20

All sources NASA NTRS DOE OSTI

DOE OSTI2025

COMSOL Results for the Nominal Steady-State Operation of the Proposed 95-MW LEU Silicide Core for HFIR Conversion

Engineering design studies are being performed to determine the feasibility of converting the High Flux Isotope Reactor (HFIR) from highly enriched uranium (HEU) to low-enriched uranium (LEU) fuel at Oak Ridge National Laboratory. This activity is sponsored by the Office of Reactor Conversion and Uranium Supply (ORCUS) under the auspices of the US Department of Energy National Nuclear Security Administration’s Office of Material Management and Minimization. HFIR is a very high flux, pressurized, light water–cooled and moderated, flux trap–type research reactor with a core made of involute shaped U 3 O 8 /Al cermet fuel plates and coolant channels. HFIR currently operates at a thermal power of 85 MW and supports key national and international missions in neutron scattering, isotope production, materials/fuels irradiation, neutron activation analysis, gamma irradiation, and neutrino research. Advanced multiphysics computational fluid dynamics models have been developed in the COMSOL Multiphysics software to simulate the steady-state operating conditions for the proposed low-and high-density LEU U 3 Si 2 -Al (uranium silicide dispersion) fuel designs. The COMSOL models for HFIR inner and outer fuel element models incorporate various essential inputs and physics such as spatially dependent nuclear heat deposition, multilayer heat conduction, conjugate heat transfer, turbulent flows (using Reynolds-averaged Navier Stokes turbulence models), structural mechanics (thermal–structural interactions and fuel swelling), and oxide layer build-up. This report presents the best-estimate thermal hydraulics results for the low- and high-density optimized silicide LEU core designs at 95 MW steady-state nominal operation.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

DOE OSTI2020

Fluid characterization using acoustics

A measurement system and a method for determining steam quality (i.e. vapor mass fraction) measurements of multiphase fluid flowing through pipes are described. An acoustic sensor device consists of an acoustic transmitter and an acoustic receiver that are designed to be attached to a pipe. The acoustic transmitter and the acoustic receiver are exposed to an interior space of the pipe through openings in a wall of the pipe. Acoustic waves generated by the transmitter and captured by the receiver traverse the multiphase fluid flowing into the pipe. Swept-frequency acoustic interferometry (SFAI) technique is used to measure ultrasonic acoustic properties of a fluid. Machine-learning techniques based on principal component analysis, support vector machine regression and support vector machine classification are used for determining steam quality.

47 OTHER INSTRUMENTATION↗

DOE OSTI2025

Impacts of Biomass Feedstock Pre-Processing on Heat and Mass Transfer During Pyrolysis Using X-Ray Computed Tomography and Multiscale Modeling

Knowledge of the transport properties of biomass particles such as porosity, tortuosity, and permeability is paramount for high-fidelity modeling of biomass pyrolysis due to the heat and mass transfer limitations imposed by particle microstructure. X-ray computed tomography (XCT) is a non-destructive imaging method that enables full 3D reconstructions of the biomass particle microstructure with high resolution, permitting direct calculation of porosity, tortuosity, and permeability from real particle geometries. In this study, XCT imaging revealed the 3D microstructures of particles and chars from pyrolytic conversion of cylindrically cut or milled/pelletized loblolly pine samples. The porosity, tortuosity, and permeability were calculated directly from the XCT geometries via open-source microstructural analysis tool MATBOX+TauFactor (https://github.com/NREL/MATBOX_Microstructure_analysis_toolbox) and computational fluid dynamics (CFD) simulations using our solver, Mesoflow (https://github.com/NREL/mesoflow). These properties were used in a reactor scale model developed in COMSOL of the single particle reactor at NREL to investigate the impact of feedstock pre-processing on biomass conversion during pyrolysis with rigorous experimental validation.

biomass↗

DOE OSTI2022

Design and performance of a variable gap system for thermal conductivity measurements of high temperature, corrosive, and reactive fluids

We report high-temperature fluids such as molten salts, liquid metals, and gasses are being proposed for many advanced energy systems including thermal energy storage devices, concentrating solar plants, and advanced nuclear reactor designs. However, the chemical behavior and thermophysical properties of many of these fluids have not been well characterized, which hinders the design, modeling, safety analysis, and deployment of these systems. Thermal conductivity is a property that is especially limited by existing measurement capabilities, which are subject to errors caused by convection, material interaction, radiative heat transfer, and instrument degradation. Therefore, there is a lack of standard, systematic measurement techniques for high-temperature, reactive, and corrosive fluids. In this work, the development of a variable gap thermal conductivity measurement system is detailed. The system is designed to measure the thermal conductivity of highly corrosive and reactive fluids, and survive operation between 100 °C and 800 °C. The effects of convection are minimized by limiting the thickness of the specimen to thin sizes (<0.3 mm). Corrections for radiative heat transfer were included in the working equations to consider specimens with varying optical properties. The design, construction, instrumentation, operating principles, and data analysis techniques are discussed in detail. The system was tested up to 500 °C using helium gas and molten KNO 3 -NaNO 3 to verify the measurement technique and determine the sources of error. At 300 and 400 °C KNO 3 -NaNO 3 , results showed maximum relative error of 6% when compared to results in the literature. The helium results were within 13% of those in the literature at 300 and 400 °C. Higher errors were observed at 500 °C for both fluids, and the sources of these errors are discussed.

42 ENGINEERING↗

DOE OSTI2022

Analysis of silica fouling on nonwetting surfaces

Ground water sources used as coolant fluids in a variety of thermal systems such as heat exchangers and power plant condensers contain silica particles that accrete on heat transfer surfaces over time leading to reduction in thermal performance, a problem that is particularly exacerbated with temperature. Nonwetting superhydrophobic, lubricant-infused, and a new class of solid-infused surfaces introduced in this work are candidates for fouling mitigation, by virtue of their water repellency, but little is known about fouling of silica on the surfaces, especially under dynamic flow conditions and as a function of temperature. This article presents, for the first time, a systematic study of dynamic flow fouling of silica on nonwetting surfaces vis-à-vis conventional copper surface over a temperature range 20–50°C. The mechanism of silica aggregate formation and its adherence to the different surfaces is elucidated by scanning electron microscope (SEM) imaging. Sigmoidal growth model is used to describe the time evolution of fouling thermal resistance and an Arrhenius model is presented for the temperature-dependent increase in the asymptotic fouling resistance on nonwetting and conventional surfaces alike. Lubricant-infused and solid-infused surfaces are shown to reduce fouling thermal resistance by up to 25% and 13%, respectively, compared to a conventional surface, whereas superhydrophobic surfaces lose their non-wettability under flow conditions, leading to an adverse increase in the fouling resistance by up to 13%. Considering the possible lubricant depletion in lubricant-infused surfaces over prolonged exposure to a flowing fluid, solid-infused surfaces present a robust alternative.

42 ENGINEERING↗

DOE OSTI2022

Mass Transport in Membrane Systems: Flow Regime Identification by Fourier Analysis

The numerical calculation of local mass distributions in membrane systems by computational fluid dynamics (CFD) offers indispensable benefits. However, the concept to calculate such distributions in response to separate variations of operation conditions (OCs) makes it difficult to address overall, flow-physics-related questions, which require the consideration of the collective interaction of OCs. It is shown that such understanding-related relationships can be obtained by the analytical solution of the advection–diffusion equation considered. A Fourier series model (FSM) is presented, which provides exact solutions of an advection–diffusion equation for a wide range of OCs. On this basis, a new zeroth-order model is developed, which is very simple and as accurate as the complete FSM for all conditions of practical relevance. Advection-dominated blocked and diffusion-dominated unblocked flow regimes are identified (depending on a Péclet number which compares the flow geometry with a length scale imposed by the flow), which implies relevant requirements for the use of lab results for pilot- and full-scale applications. Analyses reveal the equivalence of variations of OCs, which offers a variety of options to accomplish desired flow regime changes.

Heinz, Stefan (ORCID:0000000248712416)↗

DOE OSTI2026

Impact of impurities on peeling–ballooning modes and turbulence in tokamak plasmas

This study investigates the impact of various impurity species on peeling–ballooning (PB) modes and microturbulence in tokamak plasmas through the extension of traditional two-fluid and gyro-landau-fluid (GLF) models. By incorporating finite Larmor radius (FLR) effects, the analysis provides a comprehensive understanding of impurity-driven impact and its interaction with plasma turbulence. Depending on charge state and local plasma conditions, heavy impurities may exhibit gyro-radii larger than those of main ions, which are captured in the extended GLF model presented. Following the presentation of modified two-fluid equations incorporating impurity effects, we systematically analyze the distinctions between impurity and main ion dynamics and their resultant feedback mechanisms on plasma behavior. Derivation of the linear dispersion relation enables quantification of impurity-mediated modifications to: plasma vorticity, diamagnetic drift and gyroviscous effects, electron Hall physics, and FLR effects. BOUT++ – based linear simulations corroborate this formalism, demonstrating systematic stabilization of PB modes upon impurity seeding. And then operational implications for practical impurity control strategies in tokamak devices are proposed. The results underscore the necessity of impurity management to maintain stability and optimize plasma confinement, with specific focus on how FLR effects contribute to transport dynamics. This work paves the way for enhanced modeling and simulation efforts, supporting the development of strategies to control impurity-induced turbulence and improve overall reactor performance.

BOUT++ simulation↗

DOE OSTI2025

Fluoride-Cooled High-Temperature Pebble-Bed Reactor Reference Plant Model Updates

This work presents the latest improvements to, and investigations performed with, the Fluoride-Cooled High-Temperature Pebble-Bed Reactor reference plant models for the United States Nuclear Regulatory Commission. These models, developed with the Comprehensive Reactor Analysis Bundle, or BlueCRAB, serve as the foundation for the future development of detailed design evaluation models based on license applications. BlueCRAB is the code suite proposed for non-light-water reactor systems safety analyses, and it incorporates various simulation tools developed by the Nuclear Energy Advance Modeling and Simulation program, including the Griffin code for reactor physics, the Pronghorn and SAM codes for core thermal fluids, the BISON code for solid conduction and fuel performance, and the SAM code for system analysis. The primary objective of this work is to assess the level of readiness of BlueCRAB for modeling fluoride-cooled high-temperature pebble-bed reactors. We first developed numerical models in BlueCRAB that include the key physics for this technology, ensuring an adequate level of fidelity for modeling the core performance during accident scenarios. This was followed by simulation of transient scenarios, two loss-of-forced-cooling events (one protected and one unprotected), and two control rod withdrawal events (one delayed and one prompt supercritical reactivity insertion). The analysis includes comparisons between the 2-D thermal fluid porous media models in Pronghorn and SAM, comparisons between coupled Pronghorn-Griffin and coupled SAM-Griffin models for two loss-of-forced cooling events and one control rod withdrawal event, and comparisons between SAM single-solve and domain-overlapping approaches for multi-scale thermal fluid coupling. In addition, we performed comparisons between 3-D, 2-D, and 0-D neutronic models for the two control rod withdrawal scenarios with Pronghorn-Griffin. The results show that the BlueCRAB models led to physically intuitive solutions for the scenarios examined. The changes in the various scalar and vector fields such as the neutron flux, power, temperatures, densities, pressures, and velocities are all within the expected ranges, and their distributions can be explained from the system response of the transients and the geometric and material variations. Several comparisons suggest that the porous media models in Pronghorn and SAM can lead to similar solutions, even though they are based on different methodologies. The simulations demonstrate that there are differences between the various levels of fidelity, and it is advisable to have flexible tools that can cover the breadth and depth of needs that may arise in future technical evaluations. We believe that BlueCRAB’s capabilities represent a significant asset for confirmatory analyses aimed at resolving important safety questions.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2024

The Impact of Wettability on the Co-moving Velocity of Two-Fluid Flow in Porous Media

Abstract The impact of wettability on the co-moving velocity of two-fluid flow in porous media is analyzed herein. The co-moving velocity, developed by Roy et al. (Front Phys 8:4, 2022), is a novel representation of the flow behavior of two fluids through porous media. Our study aims to better understand the behavior of the co-moving velocity by analyzing simulation data under various wetting conditions. We analyzed 46 relative permeability curves based on the Lattice–Boltzmann color fluid model and two experimentally determined relative permeability curves. The analysis of the relative permeability data followed the methodology proposed by Roy et al. (Front Phys 8:4, 2022) to reconstruct a constitutive equation for the co-moving velocity. Surprisingly, the coefficients of the constitutive equation were found to be nearly the same for all wetting conditions. On the basis of these results, a simple approach was proposed to reconstruct the relative permeability of the oil phase using only the co-moving velocity relationship and the relative permeability of the water phase. This proposed method provides new information on the interdependence of the relative permeability curves, which has implications for the history matching of production data and the solution of the associated inverse problem. The research findings contribute to a better understanding of the impact of wettability on fluid flow in porous media and provide a practical approach for estimating relative permeability based on the co-moving velocity relationship, which has never been shown before.

Engineering↗

DOE OSTI2022

Transient mixing behavior of a supercritical fluid injected into supercritical and subcritical environments

Transient, near-field mixing of a supercritical fluid injected into supercritical and subcritical environments is investigated experimentally. Fuel (propane) from a 100 μm diameter Bosch diesel injector is injected into a high-pressure test rig with a steady flow of heated, low-velocity nitrogen to control the ambient pressure and temperature. The fuel supply pressure is approximately 190 bar, allowing short injections lasting 3–4 ms to be performed in quick succession to acquire statistically significant datasets to represent the unsteady fuel–air mixing. Two test cases with the same ambient density and injector inlet pressure and temperature were chosen in an attempt to isolate effects of mixture phase behavior. Microscopic rainbow Schlieren deflectometry (RSD) is used for visualization and quantitative data acquisition. In the case of the supercritical ambient, the injector inlet and exit and ambient conditions are all supercritical. The RSD analysis reveals a mixing process similar to the single phase fluid mixing, though condensation is also observed locally, and can be explained by the adiabatic mixing theory presented in the paper. For the subcritical ambient case, the injector inlet is supercritical, but the injector exit and ambient conditions are both subcritical. The results show a rapid expansion in the near-field promoting mixing between propane and nitrogen. At downstream locations, significant condensation is also observed as predicted by the adiabatic mixing theory. As a result, the supercritical ambient case exhibited longer and narrower jet penetration compared to the subcritical ambient case.

02 PETROLEUM↗

DOE OSTI2022

Field Sensitivity Analysis of Turbulence Model Parameters for Flow Over a Wing

Reynolds-averaged-Navier-Stokes (RANS) turbulence models are a critical tool in computational-fluid-dynamics simulations of aerodynamic systems, but simulation results can be highly sensitive to RANS-model parameter choices. Sensitivity analysis can be used to quantify these impacts, and the objective of this study is to demonstrate field sensitivity analysis with respect to ten parameters in the 2003 Menter shear-stress-transport (SST) turbulence model. Here, the analysis is demonstrated for an application relevant to wind energy, namely, flow over a NACA 0015 wing at 12 degree angle of attack and a Reynolds number of 1.5 x 10 6 . We quantify sensitivity using Sobol indices and the mean-squared gradient, which are estimated using polynomial chaos and active subspace models, respectively. Our results indicate that there are substantial spatial variations in parameter sensitivities, with different sets of most-sensitive parameters near the wing, as well as in the downstream wake, consistent with the physical interpretations of the turbulence model inputs. We show that, for this particular turbulence model and flow, simultaneous dimension reduction is possible across all quantities of interest, enabling efficient exploration of model outcomes. Ultimately, this analysis provides new insights into turbulence model parameter sensitivities in incompressible flows, and also demonstrates the implementation of field sensitivity analysis for applications relevant to aerodynamics simulations.

17 WIND ENERGY↗

DOE OSTI2024

EGS Reservoir Modeling for Developing Geothermal District Heating at Cornell University: Preprint

Cornell University is pursuing development of an enhanced geothermal system (EGS) for providing heating to its main campus in Upstate New York. A ~10,000 ft (~3 km) deep vertical observation well ("CUBO") was drilled in 2022 to characterize the subsurface using wellbore logging, borehole imaging, fluid sampling, mini-frac tests, coring and drill cutting analysis. Down-hole temperatures measured at 3 km depth are about 80 degrees C, sufficiently high for direct-use heating. The well drilled through generally low porosity and low permeability Paleozoic sedimentary formations and into metamorphic basement rock, encountered at about 9,400 ft depth. Leveraging subsurface data obtained through CUBO, we investigated technical feasibility and design requirements of a doublet well system with horizontal laterals connected to a fracture network created through hydraulic fracturing. The EGS reservoir is sized to provide a nominal heat output in the range 5 to 10 MWth of continuous heating over a 15 year-lifetime with limited thermal drawdown. We applied the Gringarten Multiple parallel fracture model, the Cornell Discrete Fracture Simulator FOXFEM and the commercial simulator ResFrac to estimate required heat transfer area and design a potential hydraulic stimulation treatment. Reservoir simulations indicate that, depending on fluid flow rate and injection temperature, 2 to 3 km2 of effective fracture heat transfer area is required to supply the target heat output of 5 to 10 MWth over 15 years.

district heating↗

DOE OSTI2025

EGS Reservoir Modeling for Developing Geothermal District Heating at Cornell University

Cornell University is pursuing development of an enhanced geothermal system (EGS) for providing heating to its main campus in Upstate New York. A ~10,000 ft (~3 km) deep vertical observation well ('CUBO') was drilled in 2022 to characterize the subsurface using wellbore logging, borehole imaging, fluid sampling, mini-frac tests, coring and drill cutting analysis. Down-hole temperatures measured at 3 km depth are about 80 degrees C, sufficiently high for direct-use heating. The well drilled through generally low porosity and low permeability Paleozoic sedimentary formations and into metamorphic basement rock, encountered at about 9,400 ft depth. Leveraging subsurface data obtained through CUBO, we investigated technical feasibility and design requirements of a doublet well system with horizontal laterals connected to a fracture network created through hydraulic fracturing. The EGS reservoir is sized to provide a nominal heat output in the range 5 to 10 MWth of continuous heating over a 15 year-lifetime with limited thermal drawdown. We applied the Gringarten Multiple parallel fracture model, the Cornell Discrete Fracture Simulator FOXFEM and the commercial simulator ResFrac to estimate required heat transfer area and design a potential hydraulic stimulation treatment. Reservoir simulations indicate that, depending on fluid flow rate and injection temperature, 2 to 3 km2 of effective fracture heat transfer area is required to supply the target heat output of 5 to 10 MWth over 15 years.

district heating↗

DOE OSTI2020

Experimental Investigation of Polymer Induced Fouling of Heater Tubes in the First-ever Polymer Flood Pilot on Alaska North Slope

Polymer flooding has been validated in the lab and by successful applications in other countries, such as Canada and China, to be an effective EOR technique for heavy oil reservoirs. Currently, polymer flooding is being pilot tested for the first time in the Schrader Bluff viscous oil reservoir at Milne Point field on Alaska North Slope (ANS). One of the major concerns of the operator is the impact of polymer on the oil production system after polymer breakthrough, especially the polymer induced fouling issues in the heat exchanger. This study investigates the propensity of polymer fouling on the heater tubes as a function of different variables, with the ultimate goal of determining safe and efficient operating conditions. A unique experimental set-up was indigenously designed and developed to simulate the fouling process on the heating tube. The influence of heating tube skin temperature, tube material, and polymer concentration on fouling tendency was investigated. Under each test condition, the test was run five times with the same tube, and in each run, the freshly prepared synthetic brine and polymer solution were heated from 77°F to 122°F to mimic field operating conditions. The heating time and fouling amount were recorded for each run. The morphology and composition of the deposit samples were analyzed by environmental scanning electron microscopy (ESEM) and X-ray diffraction (XRD), respectively. It has been found that, in general, the presence of polymer in produced fluids would aggravate the fouling issues on both carbon steel and stainless steel surfaces at all tested skin temperatures (165°F, 250 °F and 350°F), but only higher skin temperatures of 250°F and 350°F could cause polymer induced fouling issues on the copper tube surface and the fouling tendency increased with polymer concentration. At the lower skin temperatures of 165°F, no polymer-induced fouling was identified on the copper tube. A critical temperature that is related to the cloud point of the polymer solution was believed to exist, below which polymer induced fouling would not occur and only mineral scale deposited, but above which the polymer would obviously aggravate the fouling issue. The heating efficiency of the tube would be reduced gradually as more fouling material accumulates on its surface. The ESEM results indicate that if polymer precipitated and deposited on the surface, it would bond to the mineral crystals to form a stronger three- dimensional network structure, and that is why the polymer induced fouling was tougher and more difficult to be removed from the tube surface. The XRD analysis results confirm that the presence of polymer in the fluids can enhance the mineral scale propensity. Overall, the study results have provided practical guidance to the field operator for the ongoing polymer flooding pilot test on ANS. This study may also prove to be valuable for other chemical EOR projects around the world.

Alaska North Slope↗

DOE OSTI2021

Modification of the magneto-hydrodynamic equilibrium by the lower-hybrid wave driven fast electrons on the TST-2 spherical tokamak

Presence of a large population of fast particles may qualitatively modify the tokamak equilibrium from the one that can be described by standard magneto-hydrodynamics (MHD). The kinetic modification of the MHD equilibrium was studied for a lower-hybrid (LH) wave driven plasma on the TST-2 spherical tokamak. The analysis was performed using an equilibrium reconstruction method based on an extended MHD model that considers a two-component plasma of bulk MHD fluid and kinetic fast electrons. Scrape-off-layer current appeared naturally in the extended MHD analysis because of the finite electron orbit excursion from the flux surfaces. This reduced the bulk pressure contribution to the toroidal current substantially. Furthermore, the resulting poloidal flux profile was more consistent with that inferred from the Thomson scattering measurement.

extended magneto-hydrodynamics↗

DOE OSTI2023

Effect of the neutral beam injector operational regime on the Alfven eigenmode saturation phase in DIII-D plasma

The aim of this study is to analyze the effect of the neutral beam injector (NBI) operation regime on the saturation phase of the Alfven Eigenmodes (AEs) in DIII-D plasma. Here, the analysis is done using the linear and nonlinear versions of the gyro-fluid code FAR3d. A set of parametric analyses are performed modifying the nonlinear simulation EP $β$ (NBI injection power), EP energy (NBI voltage) and the radial location of the EP density profile gradient (NBI radial deposition). The analysis indicates a transition from the soft (local plasma relaxation) to the hard MHD (global plasma relaxation) limit if the simulation EP $β$ $\unicode{x2A7E}$ $0.02$, leading to bursting MHD activity caused by radial AEs overlapping. MHD bursts cause an enhancement of the EP transport showing ballistic-like features as avalanche-like events. Simulations in the soft MHD limit show an increment of the EP density gradient as the EP $β$ increases. On the other hand, there is a gradient upper limit in the hard MHD limit, consistent with the critical-gradient behavior. AEs induce shear flows and zonal current leading to the deformation of the flux surfaces and the safety factor profile, respectively, particularly strong for the simulation in the hard MHD limit. Simulations in the hard MHD regime show a decrease of the AE frequency in the saturation phase; this is caused by the destabilization of a transitional mode between a $9/3-10/3$ TAE and a $9/3$ RSAE that may explain the AE frequency down-sweeping observed in some DIII-D discharges. Reducing the EP energy in the nonlinear simulations leads to a weakening of the plasma perturbation. On the other hand, increasing the EP energy causes the opposite effect. Nonlinear simulations of off-axis NBI profiles indicate a lower plasma perturbation as the EP density gradient is located further away from the magnetic axis.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2021

Self-similar Reynolds-averaged mechanical–scalar turbulence models for Rayleigh–Taylor, Richtmyer–Meshkov, and Kelvin–Helmholtz instability-induced mixing in the small Atwood number limit

Analytical self-similar solutions to two-, three-, and four-equation Reynolds-averaged mechanical–scalar turbulence models describing incompressible turbulent Rayleigh–Taylor, Richtmyer–Meshkov, and Kelvin–Helmholtz instability-induced mixing in planar geometry are derived in the small Atwood number (Boussinesq) limit. The models are based on the turbulent kinetic energy K and its dissipation rate ε, together with the scalar (heavy-fluid mass fraction) variance S and its dissipation rate Χ modeled either differentially or algebraically. The models allow for a simultaneous description of mechanical and scalar mixing, i.e., mixing layer growth and molecular mixing, respectively. Mixing layer growth parameters and other physical observables relevant to each instability are obtained explicitly as functions of the model coefficients. The turbulent fields are also expressed in terms of the model coefficients, with their temporal power-law scalings obtained by requiring that the self-similar equations are explicitly time-independent. The model calibration methodology is described and discussed. Expressions for a subset of the various physical observables are used to calibrate each of the two-, three-, and four-equation models, such that the self-similar solutions are consistent with experimental and numerical simulation data corresponding to these values of the observables and to specific canonical Rayleigh–Taylor, Richtmyer–Meshkov, and Kelvin–Helmholtz turbulent flows. A calibrated four-equation model is then used to reconstruct the mean and turbulent fields, and late-time turbulent equation budgets for each instability-induced flow across the mixing layer. The reference solutions derived here can provide systematic calibrations and better understanding of mechanical–scalar turbulence models and their predictions for instability-induced turbulent mixing in the very large Reynolds number limit.

42 ENGINEERING↗

DOE OSTI2022

Automated and efficient local adaptive regression for principal component-based reduced-order modeling of turbulent reacting flows

Principal Component Analysis can be used to reduce the cost of Computational Fluid Dynamics simulations of turbulent reacting flows by reducing the dimensionality of the transported variables through projection of the thermochemical state onto a lower-dimensional manifold. However, because of the nonlinearity of the principal component source terms, nonlinear regression techniques must be utilized for the source terms in terms of the principal components. Unfortunately, widely available and utilized nonlinear regression techniques can have prohibitive computational requirements and/or accuracy that is highly dependent on user experience in ad hoc tuning of model architecture and hyperparameters. Here, in this work, a new nonlinear regression approach is proposed that is both computationally efficient and automated so does not require any user input. The approach is evaluated through a priori prediction of principal component source terms using data from a Direct Numerical Simulation of a turbulent nonpremixed n-heptane/air jet flame. In particular, the proposed framework consists of local regressions whose complexity is adapted according to the local nonlinearity of the data: local linear regression when accurate enough and local Artificial Neural Networks when nonlinear regression is required. The number of local clusters for local regression is determined automatically using the Davies-Bouldin index. In addition, Bayesian optimization is utilized for model training (i.e., to select the best architectures and hyperparameters of the nonlinear regressions in an unsupervised fashion), eliminating ad hoc hand-tuning and/or expensive grid searches. Overall, compared to a single, global neural network, the new local adaptive regression approach is shown to have comparable accuracy but 69% less training time due to the utilization of local linear regression and faster training of local neural networks.

42 ENGINEERING↗