Search NASA⌕ Search

SEARCH · Search NASA

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

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

At least 541 records · Page 30
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Revealing complex subsurface dynamics with continuous seismic monitoring: Observations using distributed acoustic sensing and surface orbital vibrators during hydraulic fracturing

Understanding hydraulic fracturing is crucial to improving the stimulation of unconventional reservoirs and increasing fluid production. This study develops a novel seismic monitoring technology using distributed acoustic sensing (DAS) and surface orbital vibrators (SOV) to capture fracture seismic response and mechanical properties at high temporal intervals. We analyze continuous time-lapse vertical seismic profiling (VSP) data acquired every hour during the first nine days of treatment of an unconventional reservoir in the Austin Chalk/Eagle Field Laboratory. The VSP data contain clear seismic signals scattered from the activated fractures. The spatiotemporal changes of the fracture reflectivity revealed by the SOV/DAS data correlate well with the observations of fracture locations inferred from low-frequency DAS data. These results capture the fracture opening and closure processes, as well as highlight potential prestage activations of the fractures due to hydraulic connectivity with preexisting fracture systems. Therefore, analysis of the presented data set provides a unique opportunity to understand fracture initiation and subsequent evolution, not only in the context of unconventional resources but also in enhanced geothermal systems.

Correa, Julia↗
DOE OSTI2021

Design and Optimization of a Gas-Cooled, Airfoil Fin Microchannel Heat Exchanger

High-performance microchannel heat exchangers are needed to supply heat for power conversion for nuclear microreactors. An airfoil fin microchannel design, constructed of Alloy 617 with helium as the working fluid, was analyzed and optimized using a design of experiments with artificial intelligence and machine learning techniques. The use of airfoil fins offers the potential to reduce pressure drop across the heat exchanger, as compared to other types of channel configurations. A framework for topology optimization of airfoil fin PCHEs has been developed that can be readily extended to different fin sizes and shapes, as well as different inlet and operating conditions, materials of construction, and working fluids. An optimization procedure was developed that employs computational fluid dynamics for a set of design points identified using Latin hypercube sampling. STAR-CCM+ was used to analyze a simplified two-channel configuration where five parameters were varied – inlet angle, fin scale, extent of staggering, transverse and longitudinal pitches. Two methods were compared for generating surrogate models – a 5D polynomial and a regression neural network. A response surface approximation was created from the surrogate models and input to a genetic algorithm. The genetic algorithm identified a set of optimal points on the Pareto front. The optimal geometry was found across six channel Reynolds numbers ranging from 1000 to 5000 to analyze how varying inlet conditions affects the optimal design. A set of optimal designs that maximizes heat transfer and minimizes pressure drop was identified, and a thermal stress analysis was performed on the optimal design. This work has developed a digital framework for the expedient topology design and evaluation of PCHE designs for gas-cooled microreactor applications. Correlations for the Nusselt number and Darcy friction factor were developed that can be useful for thermal hydraulic analyses using system codes. A thermal stress analysis was conducted and a brief discussion of the status of code cases of PCHEs for nuclear applications is given. Testing and thermomechanical modeling is needed to facilitate future code compliance of PCHEs for high pressure and high temperature applications.

42 ENGINEERING↗
DOE OSTI2020

Reduced Models for Chemical Kinetics derived from Parallel Ensemble Simulations of Stirred Reactors

Detailed chemical kinetic models are essential to improve the predictive capability of computational fluid dynamics (CFD) simulations and to impact the development of next generation combustion technologies. However, such detailed kinetic models are beyond the reach of CFD codes due to the computational cost of transporting additional chemical degrees of freedom. Reduced order models for combustion kinetics exploit the presence of lower dimensional manifolds in the state space and require far fewer number of scalars to be transported as part of CFD simulations. In this paper, we describe the software infrastructure developed to allow large ensembles of stirred reactor configurations with detailed chemical kinetics. The canonical reactor simulations are incorporated in a workflow with sparse grid and data reduction methodologies aimed at development of reduced order models. We present results from analyzing the results of the stirred reactor calculations for syngas mixture using principal component analysis (PCA) and neural network methodologies. The results show that nonlinear autoencoder model can surpass PCA in representing the original dataset with far fewer degrees of freedom.

Zhang, Pei↗
DOE OSTI2022

Accurate and Efficient Parametric Model-Order Reduction for Turbulent Thermal Transport

This project produced new algorithms and software for low-cost exploration of turbulent thermal-fluids behavior under parametric variation by using reduced-order models (ROMs). The ROMs numerically solve the governing equations for fluid motion by using a small set of basis functions (typically, N=20-200 modes) to represent the solution. The base modes are computed as optimal combinations of solutions from expensive high-fidelity "anchor-point" solutions involving millions of unknowns, which are typically generated by solving the full Navier-Stokes equations on a supercomputer. The ROM solution is itself a combination of the base modes, where the basis coefficients are determined by evolving an NxN system of nonlinear equations. The overarching idea is to use the inexpensive ROM to predict solutions under conditions where the parameters differ from the anchor-point conditions. Several ingredients are required to make ROMs useful for thermal hydraulics analysis. These include: a stable and accurate ROM that is capable of reproducing the large-scale dynamics of turbulent flow, error indicators than can guide the choice of anchor points, and low-cost mechanisms for evaluating nonlinear terms in the reduced equations.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2021

Embedding condition describing shearing, dissipative collapse

Highlights: • Perturbative approach to the time-dependent Karmarkar condition in shearing collapse. • Analysis of the (in)stability of the collapsing fluid using the TOV equation. • Study of the temperature profiles using a causal heat transport equation. In this work we employ the time-dependent Karmarkar condition to model dissipative gravitational collapse in the presence of shear which leads to a final static configuration. We believe that this is a first attempt at using the time-dependent Karmarkar condition describing a radiating star with shear in which the collapsing core dissipates energy in the form of a radial heat flux. The final static configuration which is described by an embedding class I spacetime was obtained by Hansraj and Moodly who showed that this solution satisfies all the physical requirements for a static compact object. As the collapse proceeds and the heat flux tends to zero the dynamical core evolves into a static compact core with the exterior spacetime being described by the vacuum Schwarzschild solution.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Finite Element Modeling of Electrochemical Polishing of Niobium in Hydrofluoric-Sulfuric Acid Electrolyte

Niobium (Nb) used in superconducting radio-frequency cavities requires smooth surface to achieve optimal performance. In this work, a finite element model that coupled electrochemistry, heat transfer, and fluid dynamics was developed to investigate the electrochemical polishing mechanisms of Nb, using experimentally measured polarization results of coupon samples as validations. The current and potential distribution, oxide growth kinetics of Nb in a complex cavity geometry was investigated as a function of temperature and coolant flow. A low temperature coolant with intermediate flow rate was found to reduce surface current and ensure oxide uniformity. These results could shed light on the design of future particle accelerators.

SRF accelerator cavity↗
DOE OSTI2024

Reversible Electron-Beam Patterning of Colloidal Nanoparticles at Fluid Interfaces

The directed self-assembly of colloidal nanoparticles (NPs) using external fields guides the formation of sophisticated hierarchical materials but becomes less effective with decreasing particle size. As an alternative, electron-beam-driven assembly offers a potential avenue for targeted nanoscale manipulation, yet remains poorly controlled due to the variety and complexity of beam interaction mechanisms. Here, we investigate the beam–particle interaction of silica NPs pinned to the fluid–vacuum interface of ionic liquid droplets. In these experiments, scanning electron microscopy of the droplet surface resolves NP trajectories over space and time while simultaneously driving their reorganization. With this platform, we demonstrate the ability to direct particle transport and create transient, reversible colloidal patterns on the droplet surface. By tuning the beam voltage, we achieve precise control over both the strength and sign of the beam–particle interaction, with low voltages repelling particles and high voltages attracting them. This response stems from the formation of well-defined solvent flow fields generated from trace radiolysis of the ionic liquid, as determined through statistical analysis of single-particle trajectories under varying solvent composition. Altogether, electron-beam-guided assembly introduces a versatile strategy for nanoscale colloidal manipulation, offering new possibilities for the design of dynamic, reconfigurable systems with applications in adaptive photonics and catalysis.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Identifying Potential Geochemical and Microbial Impacts of Hydrogen Storage in a Deep Saline Aquifer

Hydrogen is valuable commodity and a promising energy carrier for variable energy production. Storage of hydrogen may occur through injection of hydrogen or a hydrogen/methane gas blend in subsurface reservoirs. However, the geochemical and biological reactions that may impact the stored hydrogen are not yet understood. Therefore, we collected samples from a deep storage aquifer located in the St. Peter Formation in southern Illinois. The reservoir material was primarily quartz with sulphur and iron deposits, while the major constituents of the fluid were chloride and sulphate. 16S rRNA gene amplicon sequencing revealed a low biomass microbial community that contained no obvious hydrogen-consuming bacteria. Next, we enriched a field sample to increase the biomass and completed a metagenomic analysis, finding a low number of genes present that are associated with hydrogen consumption. Then, we completed a series of reactor experiments under reservoir conditions with 15% H2/85% CH4 gas simulating a short-term hydrogen storage, high withdrawal scenario. We found minimal changes in the geochemistry or microbiology for the reactor experiments. This work suggests that short-term storage may be highly successful, although significant additional work needs to be completed in order to accurately evaluate the risks associated with long-term hydrogen storage scenarios. It is essential we continue to expand our understanding of the dynamics present in saline aquifers and provide new insights into how hydrogen storage may impact underground geological storage environments.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Flammability and dispersion of tritium in confined release scenarios

Ignition of a flammable tritium-air mixture is the most probable means to produce the water form (T 2 O or HTO), which is more easily absorbed by living tissue and is hence ~10,000 times more hazardous to human health when uptake occurs compared to the gaseous form (T 2 or HT; per Mishima and Steele, 2002). Tritium-air mixtures with T 2 concentrations below 4 mol% are considered sub-flammable and will not readily convert to the more hazardous water form. It is therefore desirable from a safety perspective to understand the dispersion behavior of tritium under different release conditions, especially since tritium is often stored in quantities and pressures much lower than is typical for normal hydrogen. The formation of a flammable layer at the ceiling is a scenario of particular concern because the rate of dispersion to nonflammable conditions is slowest in this configuration, which maximizes the time window over which the flammable tritium may encounter an ignition source. This report describes the processes of buoyant rise and dispersion of tritium. Accumulation of flammable concentrations of tritium next to the ceiling is a common safety concern for hydrogen, but this situation can only occur if dispersion rates are slow with respect to rates of release and rise. Theory and simulations demonstrate that buoyancy does not cause regions with flammable concentrations to form within buildings from sources that have previously been mixed to sub-flammable concentrations. A simulated series of tritium release events with their associated dispersion behavior are reported herein; these simulations apply computational fluid dynamics to rooms with three different ceiling heights and a variety of tritium release rates. Safety related quantities from these simulations are reported, including the mass and volume of tritium occurring in a flammable mixture, the presence or absence of a flammable layer at the ceiling, and the time required for dispersion to nonflammable conditions after the end of the tritium release event. These safety metrics are influenced by the magnitude and rate of the tritium release with respect to the air volume in the room and also the momentum of the plume or jet with respect to the ceiling height. Several screening criteria are recommended to assess whether a specific tritium release scenario is likely to form a flammable layer at the ceiling. The methods and results in this modeling study have applicability to explosion safety analysis for other buoyant flammable gases, including the lighter isotopes of hydrogen.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2024

Design, Analysis, and Experimental Testing of Hydrogen Lean Direct Injection Nozzles at Elevated Pressure

Abstract There are many challenges of commissioning a hydrogen combustor into future gas turbine engines; especially regarding achieving emissions goals. Previously, Escudero et al. and Tran et al. conducted a study to adapt the liquid fuel Lean Direct Injection (LDI) concept from Jet-A to gaseous natural gas-hydrogen blends and pure hydrogen [1], [2]. Experimental data was collected at atmospheric conditions using a Box Behnken design of experiments. The design of experiments suggested that biasing the air split in favor of the inner air circuit and increasing the swirl strength of this inner air passage resulted in improved NOx emissions, while the inverse was true for stability, which was quantified by studying the lean blowoff point (LBO) [1], [2]. The trends revealed by the original experiment [1], [2] provided a design direction for further iterations of the experimental hardware. The study presented herein describes the further investigation of such LDI injectors through experimental methods and computational fluid dynamic (CFD) simulations at atmospheric conditions, which were used to identify potential flow behaviors driving enhanced emissions performance. Further evaluation of select injectors from both studies was then conducted at elevated pressures up to 6 atmospheres. The results from both experiments are presented in this study, which include flame observations, emissions measurements, and operational challenges. NOx emissions results are reported on a volume basis in ppmvd corrected to 15% O2 and corrected for fuel. A predictive model for relating NOx emissions to test conditions at atmospheric conditions show high significance to adiabatic flame temperature while little to no significance to fuel composition for the best performing configurations. The results illustrate the connection between atmospheric testing and testing elevated pressures. The design direction indicated by the initial tests and CFD results in promising configurations for implementation into a Multi-point LDI array.

08 HYDROGEN↗
DOE OSTI2021

A Transported Livengood–Wu Integral Model for Knock Prediction in Computational Fluid Dynamics Simulation

This work describes the development of a transported Livengood–Wu (L–W) integral model for computational fluid dynamics (CFD) simulation to predict autoignition and engine knock tendency. The currently employed L–W integral model considers both single-stage and two-stage ignition processes, thus can be generally applied to different fuels such as paraffin, olefin, aromatics, and alcohol. The model implementation is first validated in simulations of homogeneous charge compression ignition (HCCI) combustion for three different fuels, showing good accuracy in prediction of autoignition timing for fuels with either single-stage or two-stage ignition characteristics. Then, the L–W integral model is coupled with G-equation model to indicate end-gas autoignition and knock tendency in CFD simulations of a direct-injection spark-ignition engine. This modeling approach is about 10 times more efficient than the ones that based on detailed chemistry calculation and pressure oscillation analysis. Two fuels with same Research Octane Number (RON) but different octane sensitivity are studied, namely, Co-Optima alkylate and Co-Optima E30. Feed-forward neural network model in conjunction with multivariable minimization technique is used to generate fuel surrogates with targets of matched RON, octane sensitivity, and ethanol content. The CFD model is validated against experimental data in terms of pressure traces and heat release rate for both fuels under a wide range of operating conditions. The knock tendency—indicated by the fuel energy contained in the autoignited region—of the two fuels at different load conditions correlates well with the experimental results and the fuel octane sensitivity, implying the current knock modeling approach can capture the octane sensitivity effect and can be applied to further investigation on composition of octane sensitivity.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2021

Computational Fluid Dynamics (CFD) Simulations of Taylor Bubbles in Vertical and Inclined Pipes with Upward and Downward Liquid Flow

Summary Two-phase flow is a common occurrence in pipes of oil and gas developments. Current predictive tools are based on the mechanistic two-fluid model, which requires the use of closure relations to predict integral flow parameters such as liquid holdup (or void fraction) and pressure gradient. However, these closure relations carry the highest uncertainties in the model. In particular, significant discrepancies have been found between experimental data and closure relations for the Taylor bubble velocity in slug flow, which has been determined to strongly affect the mechanistic model predictions (Lizarraga-García 2016). In this work, we study the behavior of Taylor bubbles in vertical and inclined pipes with upward and downward flow using a validated 3D computational fluid dynamics (CFD) approach with level set method implemented in a commercial code. A total of 56 cases are simulated, covering a wide range of fluid properties, pipe diameters, and inclination angles: Eo ∈ [10, 700]; Mo ∈ [1×10–6, 5×103]; ReSL ∈ [–40, 10]; θ ∈ [5°, 90°]. For bubbles in vertical upward flows, the simulated distribution parameter, C0, is successfully compared with an existing model. However, the C0 values of downward and inclined slug flows where the bubble becomes asymmetric are shown to be significantly different from their respective vertical upward flow values, and no current model exists for the fluids simulated here. The main contributions of this work are (1) the relatively large 3D numerical database generated for this type of flow, (2) the study of the asymmetric nature of inclined and some vertical downward slug flows, and (3) the analysis of its impact on the distribution parameter, C0.

Engineering↗
DOE OSTI2024

Shock equation of state experiments in MgO up to 1.5 TPa and the effects of optical depth on temperature determination

Laser-driven shock compression enables an experimental study of phase transitions at unprecedented pressures and temperatures. One example is the shock Hugoniot of magnesium oxide (MgO), which crosses the B1–B2-liquid triple point at 400–600 GPa, 10 000–13 000 K (0.86–1.12 eV). MgO is a major component within the mantles of terrestrial planets and has long been a focus of high-pressure research. Here, we combine time-resolved velocimetry and pyrometry measurements with a decaying shock platform to obtain pressure–temperature data on MgO from 300 to 1500 GPa and 9000 to 50 000 K. Pressure–temperature–density Hugoniot data are reported at 1500 GPa. These data represent the near-instantaneous response of an MgO [100] single crystal to shock compression. We report on a prominent temperature anomaly between 400 and 460 GPa, in general agreement with previous shock studies, and draw comparison with equation-of-state models. We provide a detailed analysis of the decaying shock compression platform, including a treatment of a pressure-dependent optical depth near the shock front. We show that if the optical depth of the shocked material is larger than 1 μm, treating the shock front as an optically thick gray body will lead to a noticeable overestimation of the shock temperature.

36 MATERIALS SCIENCE↗
DOE OSTI2025

CFD Modeling of Concentration Polarization in High-Performance Flat-Sheet Modules

Membrane-based gas separation is a promising, regeneration-free route for industrial CO₂ capture, but its performance is often limited by concentration polarization (CP). Classical 1-D models neglect CP and thus overpredict fluxes, especially for today’s high-permeance, high-selectivity membranes. Therefore, we built a 3-D CFD model that resolves local species transport, boundary-layer resistance, permeate-side effects, and pressure drop in flat-sheet modules. Using the Concentration Polarization Coefficient (CPC) and the overall Effective Membrane CP Index (EAC), the model quantifies how geometry and operating conditions shape CP intensity. Simulations show CP becomes pronounced for advanced membranes, with feed-channel height emerging as the dominant geometric lever: shrinking the gap mitigates CP but raises pressure drop, forcing a trade-off between mass-transfer enhancement and energy cost. These insights emphasize the need for CFD-level analysis in next-generation module design, where 1-D treatments can severely underestimate required membrane area and misguide scale-up.

carbon capture↗
DOE OSTI2024

Investigation of acoustic waves under subsurface conditions to improve the predictions of rock mechanical properties and natural fracture characteristics

Mechanical properties and natural fracture characteristics are critical to investigate for subsurface engineering applications, including carbon storage, well drilling, and stimulation, as they govern rock stability, fluid flow, and mechanical behavior under stress. This dissertation integrates experimental and machine learning approaches to enhance the prediction and understanding of these properties by analyzing acoustic wave behavior under varied subsurface conditions. First, the influence of temperature, pore pressure, and supercritical CO2 (scCO2) saturation on poroelastic properties is examined using Gray Berea sandstone samples. The results show that temperature and pore pressure significantly affect the bulk modulus and Biot’s coefficient, while scCO2 saturation impacts rock compressibility, informing strategies for effective geological carbon storage. The study extends this understanding by experimentally evaluating the impact of reservoir depletion on the dynamic mechanical properties of the emerging Caney shale in South Oklahoma with the employment of unsupervised machine learning to predict static mechanical properties across the Caney shale. Integrating petrophysical data and chemostratigraphy, the workflow—featuring K-means clustering, principal component analysis (PCA), and inverse distance weighting (IDW)—improves stratigraphic characterization and the estimation of static-to-dynamic modulus ratios, which is vital for optimizing drilling and stimulation strategies. Finally, the work explores how natural fracture characteristics in shale influence acoustic waveforms and shear wave splitting (SWS) analysis. Experimental data on fractured samples under different stress and temperature conditions, combined with machine learning models such as K-nearest neighbors (KNN) and extreme gradient boosting (XGBoost), reveal key fracture properties impacting SWS and wave propagation. Together, these studies provide a comprehensive framework for linking acoustic wave behavior with rock properties, advancing the methods for monitoring and predicting geomechanical changes. The insights offered valuable implications for safer, more efficient CO2 injection, hydrocarbon extraction, and subsurface management.

Elkholy, Sherif↗
DOE OSTI2025

High-Temperature Gas-Cooled Reactors Multiphysics Simulation Demonstration and Code Validation

This study presents a comprehensive benchmarking and verification effort of several thermal-hydraulic and multiphysics capabilities for high-temperature gas-cooled reactor applications. The first part of this effort focuses on the running-in verification of Griffin’s multiphysics capabilities, specifically for simulating the evolution of pebble-bed reactor cores from startup to equilibrium. Since Fiscal Year 2024, improvements and enhancements have been implemented in Griffin, including simplifying the process to specify streamlines and developing the online cross-section generation capability. In the absence of validation data, code-to-code comparisons are conducted with kugelpy, showing good agreement for integral quantities like k-eff predictions and predictions for maximum power density. However, accuracy issues are noted for more detailed quantities like the spatial distribution of fission rate densities which will require further work to address. The second part of this report presents an improved System Analysis Module (SAM) core channel model where the effects of cross flow are considered during the pressurized loss of forced cooling transient, resulting in an improved agreement of the predicted pebble temperature with respect to the predictions from the SAM 2D porous media model. Additionally, the wall channeling effect due to variable porosity at the near wall region of the core is also investigated. Furthermore, to demonstrate Griffin’s online cross-section generation capability, a Multiphysics simulation is performed by coupling Griffin to the SAM core channel model. In the third part of the report, as a part of the Organisation for Economic Co-operation and Development/Nuclear Energy Agency (OECD/NEA) thermal-hydraulic code validation benchmark activity for a high-temperature gas-cooled reactor, the High Temperature Test Facility (HTTF) is investigated first using the NekRS computational fluid dynamics (CFD) code to study the flow mixing phenomenon in the lower plenum of the facility. Then, code-to-code and code-to-data comparisons are performed for Test PG27, which is a pressurized conduction cooldown (PCC) test, using five different codes by six organizations from five countries. The different simulations show good agreements in terms of the general trend but there are differences in some results such as the peak temperatures of different regions and heat removal rate.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2024

Macroscopic modeling of gas permeability in hierarchical micro/nanoporous media: A unified characterization of rarefaction using Klinkenberg theory and equivalent diameter

Estimating gas transport through a hierarchical micro/nanoporous system is challenging due to non-equilibrium gas dynamics. The primary difficulty lies in determining the rarefaction level, because identifying a representative flow dimension in a complex porous system with multiple pore scales is not straightforward. Our study performed a pore-level analysis for gas permeability in dual-scale porous media with varying porosity, throat size, and secondary pore size under different rarefaction conditions. We found that secondary porosity negatively affects permeability due to increased friction forces, with this influence growing as the secondary pore size and porosity increase until the secondary pore becomes comparable to the throat. However, rarefaction reduces the effects of secondary pores due to boundary slip. Traditional Knudsen number (Kn) calculations based on Darcy-defined height failed to accurately describe the rarefaction effects on gas permeability. Instead, we introduced an equivalent diameter to calculate the Kn, which provided an accurate normalization of apparent gas permeability independent of pore geometry. Furthermore, the extended Kozeny–Carman–Klinkenberg model developed in our previous study successfully yielded a macroscopic model for apparent gas permeability in hierarchical micro/nanoporous systems as a function of the traditional Darcy height and porosity.

04 OIL SHALES AND TAR SANDS↗