Search NASA⌕ Search

SEARCH · Search NASA

Results for “Computational Fluid Dynamics”

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 289 records · Page 16
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

High-Fidelity CFD Modeling of Cryogenic Hydrogen Isotope Extrusion for Fusion Reactor Pellet Fueling

This study investigates the extrusion processes of deuterium and protium using ANSYS-Polyflow. The geometries and computational fluid dynamics (CFD) settings closely replicate the experimental setups and data acquired from the extruder experiments at Oak Ridge National Laboratory (ORNL) for validation purposes. We explore the impacts of (1) slip versus non-slip boundary conditions and (2) the use of constant, temperature-, and shear rate–dependent viscosities, concluding that the implementation of non-slip wall boundary conditions combined with shear rate–dependent viscosity produced more accurate predictions. The simulations achieved excellent agreement with the experimental data, with relative differences of only 5% for deuterium, and 3% to 6% for protium. This is the first time that experimental extrusion data at ORNL have been accurately predicted through high-fidelity CFD modeling. In conclusion, the advancements offer valuable insights and a foundational modeling tool for optimizing pellet injectors for ITER and other future reactor-scale devices.

ANSYS-Polyflow↗
DOE OSTI2023

In-silico quantification of intersubject variability on aerosol deposition in the oral airway

The extrathoracic oral airway I s not only a major mechanical barrier for pharmaceutical aerosols to reach the lung but also a major source of variability in lung deposition. Using computational fluid dynamics, deposition of 1-30 µm particles was predicted in 11 CT-based models of the oral airways of adults (7 healthy subjects and 4 patients with chronic obstructive pulmonary disease). Simulations were performed for mouth breathing during both inspiration and expiration at two steady-state flow rates representative of resting/light activity (18 L/min) and of that typically used during aerosol inhalation from a dry powder inhaler (45 L/min). For each subject, oral deposition increased with increasing particle size and flow rate. When averaged over all subjects, there was no significant difference in oral deposition efficiency between inspiration and expiration, unlike subregional deposition that show significantly different patterns between the two breathing phases. Consistent with previous in-vitro inspiratory studies, our data showed a large intersubject variability in oral deposition. These results highlight the need for incorporating morphological variation of the upper airway in predictive models of aerosol deposition for accurate predictions of particle dosimetry in the intrathoracic region of the lung.

Computational fluid dynamics (CFD), Inertial impac↗
DOE OSTI2023

Machine Learning-Based Regression Models for Ironmaking Blast Furnace Automation

Computational fluid dynamics (CFD)-based simulation has been the traditional way to model complex industrial systems and processes. One very large and complex industrial system that has benefited from CFD-based simulations is the steel blast furnace system. The problem with the CFD-based simulation approach is that it tends to be very slow for generating data. The CFD-only approach may not be fast enough for use in real-time decisionmaking. To address this issue, in this work, the authors propose the use of machine learning techniques to train and test models based on data generated via CFD simulation. Regression models based on neural networks are compared with tree-boosting models. In particular, several areas (tuyere, raceway, and shaft) of the blast furnace are modeled using these approaches. The results of the model training and testing are presented and discussed. The obtained R 2 metrics are, in general, very high. The results appear promising and may help to improve the efficiency of operator and process engineer decisionmaking when running a blast furnace.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Ocean Energy Systems Wave Energy Modeling Task 10.4: Numerical Modeling of a Fixed Oscillating Water Column

This paper reports on an ongoing international effort to establish guidelines for numerical modeling of wave energy converters, initiated by the International Energy Agency Technology Collaboration Program for Ocean Energy Systems. Initial results for point absorbers were presented in previous work, and here we present results for a breakwater-mounted Oscillating Water Column (OWC) device. The experimental model is at scale 1:4 relative to a full-scale installation in a water depth of 12.8 m. The power-extracting air turbine is modeled by an orifice plate of 1–2% of the internal chamber surface area. Measurements of chamber surface elevation, air flow through the orifice, and pressure difference across the orifice are compared with numerical calculations using both weakly-nonlinear potential flow theory and computational fluid dynamics. Both compressible- and incompressible-flow models are considered, and the effects of air compressibility are found to have a significant influence on the motion of the internal chamber surface. Recommendations are made for reducing uncertainties in future experimental campaigns, which are critical to enable firm conclusions to be drawn about the relative accuracy of the numerical models. It is well-known that boundary element method solutions of the linear potential flow problem (e.g., WAMIT) are singular at infinite frequency when panels are placed directly on the free surface. This is problematic for time-domain solutions where the value of the added mass matrix at infinite frequency is critical, especially for OWC chambers, which are modeled by zero-mass elements on the free surface. A straightforward rational procedure is described to replace ad-hoc solutions to this problem that have been proposed in the literature.

16 TIDAL AND WAVE POWER↗
DOE OSTI2024

Numerical modeling of Joule heated ceramic melter

Computational fluid dynamics and heat transfer models of waste glass melters are being developed to support the vitrification campaign at the Waste Treatment and Immobilization Plant (WTP) as part of the United States nuclear waste stewardship and environmental protection programs. The WTP mission aims to achieve the stable and durable processing of vitrified waste in the safest, most efficient, and economical manner possible. Hanford’s tank waste is far more complex and varied than the waste treated by other vitrification facilities and both startup and long-term operational challenges are likely to be even greater than prior experience has shown. There are thousands of different waste compositions, which makes processing the waste a complex and challenging task. Models of test melters, pilot-scale, and full-scale melters have been developed to provide understanding and insight into the various physical and chemical phenomena occurring during vitrification. So the results of these modeling activities will be used to support the WTP operations to process legacy nuclear waste into a stable form for disposal.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2026

Computational flow modeling of triply periodic minimal surfaces as feed channel spacers in ultra-high pressure reverse osmosis applications

Triply periodic minimal surfaces (TPMS) are a special class of mathematical surfaces characterized by a high surface area-to-volume ratio. They have generated considerable interest in fields such as acoustics, heat transfer, and membrane-based filtration processes. This study evaluates the performance of four different TPMS designs—Schoen Gyroid, Schoen Crossed Layers of Parallels (CLP), Schoen Transverse Crossed Layers of Parallels (tCLP), and Schwarz-Primitive—when used as feed channel spacers under ultra-high pressure reverse osmosis (UHPRO) conditions, at approximately 200 bar. Our experimentally validated computational fluid dynamics model reveal different flow patterns within the feed channels for each of the four TPMS designs, leading to varying hydrodynamic and permeation properties. Under the simulated UHPRO conditions, the Gyroid and tCLP designs yield up to a 23% increase in average permeate velocity and a 14% reduction in average membrane-surface concentration relative to a non-woven spacer of the same porosity. Furthermore, the enhanced performance comes with an increased feed channel pressure drop, although it only constitutes less than 4% of the operating pressure when extrapolated for a meter-long membrane module. Additionally, the study analyzes the effects of varying inlet velocity and spacer porosity on membrane performance. Overall, this research provides valuable insights into the potential use of TPMS spacers in UHPRO applications.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Temperature Rise on Liner Surfaces of Fuel Cell Electric Vehicle Tanks during Fueling Process

Herein, the temperature of the inner tank walls, or plastic liners, of composite pressure tanks in fuel cell electric vehicles during fueling using three‐dimensional computational fluid dynamics (CFD) models is evaluated. The liner materials are limited to 85.0 °C to prevent thermal stress causing material failure that would result in a hydrogen leak. Therefore, the temperature of hydrogen gas must be limited to below 85.0 °C during the fueling process as dictated by the current fueling protocol. However, there are limited experimental or simulation data confirming that the temperature changes do not exceed the threshold. Herein, the liner temperatures with CFD tank models for two sizes of type IV tanks representative of the upper and lower system bounds that are close to the SAE J2601 fueling protocol (36.0 and 244.0 L) are evaluated. First, each model's reliability is validated with experimental data and then analyzed, and the data are used to evaluate the maximum hydrogen and liner temperatures under real‐world fueling conditions. The evaluation shows that the maximum liner surface temperature of each tank model is at least 7 °C lower than that of the hydrogen. Additionally, there is at least 12 °C difference found between the upper limit and actual liner temperatures.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2020

Analysis of combustion acoustic phenomena in compression–ignition engines using large eddy simulation

As computational capabilities continue to grow, exploring the limits of computational fluid dynamics to capture complex and elusive phenomena, which are otherwise difficult to study by experimental techniques, is one of the main targets for the research community. This paper presents a detailed analysis of the physical processes that lead to combustion noise emissions in internal combustion engines. In particular, diesel combustion in a compression-ignition (CI) engine is studied in order to understand the singular behavior of the in-cylinder flow field responsible for the acoustic emissions. The main objective is, therefore, to improve the understanding of the phenomena involved in CI engine noise using large eddy simulations. Several visualization methods are employed to investigate the connection between combustion behavior and its effects on the pressure field. In addition, proper orthogonal decomposition is used to analyze the modal energy distribution among all the acoustic modes. The results show that the acoustic signature is fundamentally conditioned by the intensity of the premixed combustion rather than by the pressure oscillations generated by turbulent fluctuations in the flame surface established during the diffusion stage.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Porous Flow Modeling of Axial Gas Redistribution in Fragmented LWR Fuel Rods using MOOSE

Understanding how gas axially redistributes within fragmented fuel pellets is crucial for predicting the behavior of Light Water Reactor (LWR) fuel rods, particularly during transient and accidental scenarios. The time scale of this phenomenon plays a fundamental role in determining the progression and hazard of a Loss Of Coolant Accident (LOCA), especially when high burn-up fuel in a severe state of fragmentation is involved. Here, this study presents a Computational Fluid Dynamics (CFD) model developed within the Multiphysics Object-Oriented Simulation Environment (MOOSE) to predict the time-scale of plenum depressurization in Light-Water Reactor (LWR) fuel rods driven by axial gas transport through fragmented pellets. The model examines the effects of incorporating non-linearities in the friction term by comparing the results with experimental data. The experiment employed surrogate fuel rods containing pellets subjected to mechanical and/or thermal loadings to simulate various severity of cracking, and aimed at studying the influence of fuel conditions on axial gas redistribution. The results of this analysis indicate that under certain flow regime conditions - determined by the value of an equivalent Reynolds number - accounting for the non-linear friction term in Navier-Stokes equations guarantees better predictions for the time-scale of plenum depressurization. Also, the model enabled the simulation of the pressure decay by assigning distinct permeability values to each pellet instead of a single uniform value. Multiple simulations were run across all possible pellet position combinations, having each pellet assigned with values of permeability extracted from the experimental data. This allows to quantify the impact of the considering various non-uniform distributions of permeability on the dynamics of axial gas redistribution. The present work findings enhance the understanding of axial gas transport, and provide valuable insights for the integration of a model for predicting the axial gas redistribution during a LOCA scenario into the BISON fuel performance code.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2023

A reduced-order model of concentration polarization in reverse osmosis systems with feed spacers

Feed spacers in reverse osmosis systems generate complex fluid flows that limit computational fluid dynamics (CFD) simulations to small length and time scales. That limits our ability to simulate mineral scaling and other membrane fouling phenomena, which occur over longer length and time scales. Thus motivated, we develop a reduced model that replaces the CFD simulation of the velocity field with an analytical model that mimics spacers. This focuses the remaining numerical effort on simulating the advection–diffusion equation governing solute transport. We motivate and validate the model with CFD simulations and bench-scale experiments of spacer filaments in three different arrangements, including cases of unsteady vortex shedding. We show that the model produces a roughly 10,000-fold speedup compared to CFD, and accurately reproduces CFD predictions of not only the average and maximum concentrations, but also the local concentration distribution along the membrane. Here, we also demonstrate the model for simulating a feed channel with a length-to-height ratio of 200. The model provides a simple testbed for exploratory studies of multispecies transport, precipitation, and membrane fouling phenomena for which simulating spacers is often prohibitive.

Computational fluid dynamics↗
DOE OSTI2021

Thermal-striping analysis methodology for sodium-cooled reactor design

Computational fluid dynamics simulations have been performed to study the applicability of engineering CFD methods for thermals striping analysis. Thermal striping is the fluctuating temperature profile in a solid caused by fluctuating fluid temperature, and the resulting fluctuating thermal stresses cause high cycle fatigue and eventual material failure. Here this study presents the methodology for thermal striping analysis that including a transient conjugate heat transfer RANS model and benchmark the results against a sodium triple jet experiment. The multiphysics calculations are implemented to analyze the thermal stresses in the solid domain caused by the coupled heat transfer between the fluid and solid domains. Validation data from a liquid sodium triple jet experiment include time-averaged temperature measurements and power spectra of the temperature signal. The numerical results agree well with these experimental measures, demonstrating key features such as the dominant frequency of temperature fluctuations. The applicability of a low-cost wall treatment method is demonstrated, enabling key computational savings. Finally, the performance of two finite element stress analysis software packages is compared, and the validity of the lower-cost method is confirmed. These results demonstrate the applicability of engineering methods for computational thermal striping calculations, enabling thermal striping estimations in large fluid systems such as the core of a Sodium-cooled nuclear reactor.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2020

Investigation of countercurrent flow profile and liquid holdup in random packed column with local CFD data

Liquid holdup and mass transfer area are critical parameters used in packed column design and CO 2 capture efficiency prediction. Compared to experimental methods, 3D computational fluid dynamics (CFD) can provide high-resolution data for studying countercurrent flow hydrodynamics in a packed column. In this paper, a framework was established for modeling the liquid–gas countercurrent flow in a random packed column with pall rings. Besides the column-averaged information, the radial pall ring distribution, velocity, and liquid holdup profiles are obtained in the simulations. The radial two-phase velocity and liquid holdup profiles show significant fluctuations caused by the non-uniform pall ring distribution. The entrance effect and wall influence are also observed in the packed column. By using local countercurrent flow information, the packing specific area $a_p$ and liquid velocity $u_L$ range are significantly expanded for liquid holdup hL prediction. This enables us to develop the liquid holdup correlation with one packing size and six CFD runs only. The proposed power law of the liquid holdup correlation $h_L~u_L^{0.44}$ indicates the random packed column falls in a viscous to turbulent transition regime with a Reynolds Number range of [6.7–40.2]. The derived liquid holdup correlation is in good agreement with existing correlations that are developed using the column-averaged experimental data.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Machine learning changes the rules for flux limiters

Learning to integrate non-linear equations from highly resolved direct numerical simulations has seen recent interest for reducing the computational load for fluid simulations. Here, we focus on determining a flux-limiter for shock capturing methods. Focusing on flux limiters provides a specific plug-and-play component for existing numerical methods. Since their introduction, an array of flux limiters has been designed. Using the coarse-grained Burgers' equation, we show that flux-limiters may be rank-ordered in terms of their log-error relative to high-resolution data. We then develop a theory to find an optimal flux-limiter and present flux-limiters that outperform others tested for integrating Burgers' equation on lattices with [Formula: see text], and 2x, 3x, 4x, and 8x coarse-grainings. We train a continuous piecewise linear limiter by minimizing the mean-squared misfit to six-grid point segments of high-resolution data, averaged over all segments. While flux limiters are generally designed to have an output of φ(r) = 1 at a flux ratio of r = 1, our limiters are not bound by this rule and yet produce a smaller error than standard limiters. Here we find that our machine learned limiters have distinctive features that may provide new rules-of-thumb for the development of improved limiters. Additionally, we use our theory to learn flux-limiters that outperform standard limiters across a range of values (as opposed to at a specific fixed value) of coarse-graining, number of discretized bins, and diffusion parameter. This demonstrates the ability to produce flux limiters that should be more broadly useful than standard limiters for general applications.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Modeling of Supercritical CO2 Shell-and-Tube Heat Exchangers Under Extreme Conditions: Part II: Heat Exchanger Model

Abstract Heat exchangers play a critical role in supercritical CO2 Brayton cycles by providing necessary waste heat recovery. Supercritical CO2 thermal cycles potentially achieve higher energy density and thermal efficiency operating at elevated temperatures and pressures. Accurate and computationally efficient estimation of heat exchanger performance metrics at these conditions is important for the design and optimization of sCO2 systems and thermal cycles. In this paper (Part II), a computationally efficient and accurate numerical model is developed to predict the performance of shell-and-tube heat exchangers (STHXs). Highly accurate correlations reported in Part I of this study are utilized to improve the accuracy of performance predictions, and the concept of volume averaging is used to abstract the geometry and reduce computation time. The numerical model is validated by comparison with computational fluid dynamics (CFD) simulations and provides high accuracy and significantly lower computation time compared to existing numerical models. A preliminary optimization study is conducted and the advantage of using supercritical CO2 as a working fluid for energy systems is demonstrated.

Engineering↗
DOE OSTI2021

Pore-scale CFD simulations of clay mobilization in natural porous media due to fresh water injection

The present work investigates mechanisms of permeability impairment as a result of low-salinity fluid injection into brine-saturated porous media containing dispersible clays. Here, we present a computational fluid dynamics model at the pore-scale to simulate detachment, migration and straining of fine particles in porous media. The model uses an immersed boundary method to simulate the motion of clay fines in a fluid. In addition to the hydrodynamic forces, we model the Derjaguin-Landau-Verwey-Overbeek forces (DLVO) between clay fines and grains. Our simulations show the impact of the injected fluid's salinity and velocity on the concentration of clay fines retained on the grain surface. We see clay particles dislodging from the grain surface in clusters of up to 12 particles. Our simulation results also demonstrate clogging of narrow pore spaces by the detached particles from the grain surfaces.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Development of a Filtered CFD-DEM Drag Model with Multiscale Markers Using an Artificial Neural Network and Nonlinear Regression

Here, the accuracy of coarse-grained Euler-Lagrangian simulations of fluidized beds heavily depends on the mesoscale drag models to account for the influences of the unresolved sub-grid structures. Traditional filtered drag models are regressed with mesoscale markers such as voidage and slip velocities. In this research, a filtered drag was regressed with both mesoscale and macro-scale markers using fine grid Computational Fluid Dynamics - Discrete Element Method (CFD-DEM) simulations. The traditional non-linear regression method was compared with machine learning regression using an Artificial Neural Network (ANN) implemented in PyTorch and coupled with MFiX. The new drag showed higher accuracy than the Wen-Yu drag and another filtered drag derived from the two-fluid model. The nonlinear regression shows slightly better results than ANN regression in cases with similar R 2 values. The utilization of the gas inlet velocity as an additional macro-scale marker reduced the errors by up to 55.3% in the tested cases.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

High-Fidelity Numerical Investigation on Elucidating Sodium Heat Transfer Characteristics for 37-Pin Wire-Wrapped Fuel Bundle in the PLANDTL Facility

This study involved a Reynolds-averaged Navier-Stokes- (RANS-) based computational fluid dynamics (CFD) analysis of the 37-pin wire-wrapped fuel bundle of the PNC Plant dynamics test loop (PLANDTL) facility. Previously, mainly the hydrodynamic phenomena of the wire-wrapped fuel bundle were analyzed, but the present study additionally included heat transfer analysis through conjugate heat transfer. The main purpose of the study was to benchmark the experimental data of the PLANDTL 37-pin wire-wrapped fuel bundle to investigate the heat transfer phenomena. In addition, the aim was to verify the accuracy of the RANS-based CFD analysis method using the STAR-CCM+ simulation software in comparison with the experimental data. The grid used for verification was an innovative grid system consisting of hexahedra using Fortran-based code. The development of the RANS-based CFD methodology included grid sensitivity analysis, turbulence model sensitivity analysis, and turbulent Prandtl number sensitivity analysis. Information on the temperature, mass flow rate, and area of the CFD results for each subchannel was provided for the top of the heated section and is expected to serve as a reference for future studies aiming to perform the validation and verification of a PLANDTL facility. In addition, the dependence of the peak temperature on the azimuth angle of each pin was analyzed.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Fuel-air mixing in motored CFR engine at research octane number (RON) relevant condition

This paper presents a three-dimensional (3-D) computational fluid dynamics (CFD) study of a motored cooperative fuel research (CFR) engine at research octane number (RON) relevant condition. The boundary conditions for 3-D simulations were generated with a one-dimensional GT-Power model. For the first time in literature, a carburetor was added to a virtual CFR engine model with 3-D CFD. Therefore, the proposed setup can simulate the fuel and thermal stratifications inside the engine cylinder with realistic detail. The transient simulations in this work were performed within the Reynolds-averaged Navier-Stokes (RANS) framework with a Realizable k-ε turbulence model. Major conclusions from the present work are: (1) The in-cylinder flow of the CFR engine is swirl-dominated due to the existence of the intake valve shroud. (2) There is a significant amount of liquid droplets entering the cylinder during the intake stroke. The maximum instantaneous amount of liquid for 50% PRF 87 (containing 87% iso-octane and 13% n-heptane (v/v)) and 50% ethanol mixture is indicated to be around 26% of total injected fuel mass. (3) The heat of vaporization (HoV) of the fuel is responsible for creating both temperature and charge stratification inside the cylinder.

33 ADVANCED PROPULSION SYSTEMS↗