Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Navier-Stokes”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

A fourth order sharp immersed method for the incompressible Navier-Stokes equations with stationary and moving boundaries and interfaces

We propose a fourth order Navier-Stokes solver based on the immersed interface method (IIM), for flow problems with stationary and one-way coupled moving boundaries and interfaces. Our algorithm employs a Runge-Kutta-based projection method that maintains high-order temporal accuracy in both velocity and pressure for steady and unsteady velocity boundary conditions. Fourth order spatial accuracy is achieved through a novel fifth order IIM discretization scheme for the advection term, as well as existing high-order interface-corrected finite difference schemes for the other differential operators. Using a set of manufactured flow problems with stationary and moving boundaries, we demonstrate fourth order convergence of velocity and pressure in the infinity norm, both inside the domain and on the immersed boundaries. The solver’s performance is further validated through a range of practical flow simulations, highlighting its efficiency over a second order scheme. Finally, we showcase the ability of our immersed discretization scheme to handle interface-coupled multiphysics problems by solving a conjugate heat transfer problem with multiple immersed solids. Overall, the proposed approach robustly combines the efficiency of high order discretization schemes with the flexibility of immersed discretizations for flow problems with complex, moving boundaries and interfaces.

42 ENGINEERING↗
DOE OSTI2026

Large Eddy Simulation of Low-Reynolds-Number Turbulent Flow of Low-Prandtl-Number Fluid in a Tight Lattice Bundle for Assessment of Reynolds-Averaged Navier-Stokes Turbulence Model

The MARVEL (Microreactor Applications Research Validation and Evaluation) microreactor utilizes natural circulation as core cooling mechanism and liquid metal as a primary coolant. Moreover, the reactor core has a pitch-to-diameter ratio of 1.056, which is considered a tight lattice configuration. Numerous studies have widely reported that Reynolds-Averaged Navier-Stokes (RANS) turbulence models inaccurately predict heat transfer in liquid metals and fail to capture flow pulsations that can occur within tight lattices, leading to further inaccuracies in simulation results. Therefore, evaluating the accuracy of RANS turbulence models in the thermal-hydraulic analysis of the MARVEL microreactor core is crucial for assessing reactor safety. In this study, a Large Eddy Simulation (LES) of the MARVEL microreactor core subchannel was conducted and compared with a RANS simulation to evaluate the accuracies and conservatism of the RANS model. The flow pulsation in a tight lattice predicted by LES enhanced the heat transfer, whereas the RANS model underpredicted it. Consequently, the RANS model predicted the peak cladding temperature higher than the LES model, but the discrepancy between the two approaches was not significant due to the good thermal characteristics of the liquid metal. It can be concluded that the steady-state RANS model is effective for the thermal analysis of liquid-metal-cooled MARVEL microreactor core and can provide conservative predictions from a safety analysis perspective.

21 - SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLAN↗
NASA NTRS1995

A Navier-Stokes boundary element solver

Using global interpolation functions (GIF's) boundary element solutions are obtained for two-dimensional laminar flows. Two schemes are proposed for handling the convective terms. The first treats convection as a forcing function, and converts the flow equations to pseudo-Poisson equations. In the second scheme, some convective effect is incorporated into the fundamental solution used in constructing the pertinent integral equations. The lid-driven cavity flow is selected as the benchmark problem.

Reddy, D. R.↗
DOE OSTI2026

Numerical Simulation of a Natural Convection–Driven Air-Cooled Reactor Cavity Cooling System Experiment

Ensuring the efficient removal of decay heat from the reactor vessel is essential for the safety of advanced reactor technologies. Several Generation-IV concepts incorporate variations in the reactor vessel cooling systems to achieve this objective. High-temperature gas-cooled reactors utilize a reactor cavity cooling system (RCCS), a passive ex-vessel system designed to operate without active components or external power during accident conditions. The RCCS removes decay heat primarily through radiative and convective heat transfer mechanisms. Here, this study presents a comprehensive validation of a computational fluid dynamics Reynolds-averaged Navier-Stokes model for the University of Wisconsin-Madison air-cooled RCCS facility. Validation was conducted for both high- and low-power natural convection cases under a uniform heating profile. Near-wall resolution was found to be critical for accurately modeling natural convection in the RCCS; employing an all-𝑦 + wall treatment resulted in wall temperature discrepancies exceeding 50 °⁢𝐶 compared to a wall-resolved mesh. Thermal-hydraulic behaviors under natural and forced convection conditions were compared within the heated cavity and RCCS. A turbulence model sensitivity analysis indicated that low-Reynolds number k-ɛ, k-ω shear stress transport (SST), and Reynolds stress transport models produce similar wall temperature predictions. A buoyancy modeling sensitivity study revealed that the Boussinesq approximation significantly underpredicted thermal-hydraulic behavior in the RCCS. Based on these findings, modeling recommendations are provided. The validated data set along with identified sensitivities refine the modeling of natural convection in the RCCS. The information produced by this study supports RCCS design, optimization, and safety evaluations, enabling the calibration and verification of reduced-order thermal-hydraulic models.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2025

Analysis of the weighted shifted boundary method for the Poisson and Stokes problems

The Shifted Boundary Method (SBM) belongs to the class of unfitted (or immersed, or embedded) finite element methods, and relies on reformulating the original boundary value problem over a surrogate (approximate) computational domain. Accuracy is maintained by properly shifting the location and values of the boundary conditions. This avoids integration over cut cells and the associated implementation issues. Recently, the Weighted SBM (WSBM) was proposed for the Navier-Stokes equations with free surfaces and the Stokes flow with moving boundaries. The attribute “weighted” in the name WSBM stems from the fact that its variational form is weighted with the elemental volume fraction of active fluid. The motivation for the development of the WSBM was the preservation of the volume of active fluid to a higher degree of accuracy, which in turn resulted in improved stability and robustness characteristics in moving-boundary, time-dependent simulations. In this article, we present the numerical analysis of the WSBM formulations for the Poisson and Stokes problems. We give mathematical conditions under which the bilinear forms defining the discrete variational formulations are uniformly coercive (Poisson problem) or inf-sup stable (Stokes problem). By these results, stability and optimal convergence is proven in the natural norm; L2-error estimates can also be derived.

Approximate domain boundaries↗
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 OSTI2026

Recent progress on coarse graining simulations

We focus on coarse graining simulations based on the primary conservation equations, effectively codesigned physics and algorithms, and low-Mach-number corrected (LMC) hydrodynamics. Simulation methods involve LANL’s x-Radiation-Adaptive-Grid-Eulerian Large-Eddy Simulation, Besnard-Harlow-Rauenzahn (BHR) Reynolds-Averaged Navier-Stokes (RANS) approach, and Dynamic BHR – a paradigm bridging RANS and LES. A relevant question addressed relates to whether 3D RANS and RANS/LES hybrids – the industry standards for aerospace and automotive research, are presently relevant for practical variable-density applications involving shocked and accelerated interface instabilities. Furthermore, recent simulations of the GaTECH inclined mixing-layer shock-tube and NIF ICF-capsule experiments are used to demonstrate issues, challenges, and potential for 3D coarse grained LMC simulation strategies for robustly simulating complex transitional and coupled hydrodynamics-multiphysics with coarser resolution. Present LES readiness to provide accurate predictions at scale is demonstrated – whereas 3D RANS and RANS/LES bridging do not appear impactful in this context.

42 ENGINEERING↗
DOE OSTI2026

Development of MOSCATO: A CFD-Level Electrochemistry and Corrosion Simulator for Molten Salt Systems

For both coolant and fueled variants of molten salt reactors (MSRs), the corrosion of structural materials is a significant challenge. The corrosion stems from chemical and electrochemical reactions initiated by fissile material, fission products, and impurities in the salt. Lower-fidelity models rely on empirical correlations for mass transfer, simplified lumped temperature profiles, and similar assumptions. They do not capture detailed spatial variations in complex geometries, creating the need for high-fidelity modeling to bridge this gap.As we approach the demonstration and possible deployment of MSRs in this decade, the development of a high-fidelity, high-performance simulator becomes imperative. To simulate the complex electrochemical environment and corrosion within molten salt systems, we have developed the Molten Salt Chemistry And TranspOrt (MOSCATO) code. This endeavor is comprised of three essential components. First, mass transfer equations are coupled with the Navier-Stokes equations in order to account for the transport of species in the salt. Second, the diffusion of alloy constituents, such as Cr, Fe, Ni, etc. is simulated within the structural metals. Third, the alloy and salt domains are coupled to account for the heterogeneous chemical and electrochemical reactions that occur at the salt-alloy interface.MOSCATO manages all three components within the framework of the highly scalable, open-source spectral element method computational fluid dynamics code Nek5000/NekRS. This integration enables MOSCATO to harness the immense computational power of modern high-performance computing resources, ensuring both high fidelity and computational speed.In addition to code development, we have initiated a comprehensive verification and validation campaign, utilizing data from diverse sources. First, MOSCATO's electrochemical solver was verified with reference numerical data. Then validation occurred against experiments: one of a thermal galvanic cell and the other for corrosion in flowing molten salt of FLiNaK (LiF-NaF-KF). This campaign verified and validated MOSCATO as a reliable tool for simulating electrochemical environments and corrosion in molten salt systems.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2026

Emergent viscous hydrodynamics from a single quantum particle

We investigate an explicit example of how spatial decoherence can lead to hydrodynamic behavior in the late-time, long-wavelength regime of open quantum systems. We focus on the case of a single nonrelativistic quantum particle linearly coupled to a thermal bath of noninteracting harmonic oscillators at temperature T , a la Caldeira and Leggett. Taking advantage of decoherence in the position representation, we expand the reduced density matrix in powers of the off-diagonal spatial components, so that high-order terms are suppressed at late times. Truncating the resulting power series at second order leads to a set of dissipative transient hydrodynamic equations similar to the nonrelativistic limit of equations widely used in simulations of the quark-gluon plasma formed in ultrarelativistic heavy-ion collisions. Transport coefficients are directly determined by the damping constant γ , which quantifies the influence of the environment. The asymptotic limit of our hydrodynamic equations reduces to the celebrated Navier-Stokes equations for a compressible fluid in the presence of a drag force. Furthermore, our results shed new light on the onset of hydrodynamic behavior in open quantum systems where a system with few degrees of freedom is coupled to a large thermal environment.

Hydrodynamics↗
DOE OSTI2026

Radial oscillations of viscous neutron stars: Zero diffusion case

In this article, the spectrum of radial oscillations of neutron stars is systematically studied within two frameworks of viscous relativistic hydrodynamics: the relativistic Navier-Stokes and Israel-Stewart theories. A correspondence is established between the discrete stellar eigenmodes and the continuous dispersion relation of perturbations around a homogeneous fluid, providing a basis for interpreting our numerical results. We analyze the Newtonian limit and assess the impact of relativistic corrections, such as the gravitational redshifting of microscopic relaxation timescales. We show that bulk viscosity can significantly affect the behavior of both hydrodynamic and nonhydrodynamic fundamental modes, and that, depending on the magnitude of the viscous effects, it is the nonhydrodynamic mode that becomes unstable beyond the turning point in a sequence of equilibrium configurations. These results provide a useful step toward systematic studies of neutron star quasinormal modes in the presence of viscosity.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

HygroThermFEM v1.0

HygroThermFEM is a Finite Element Method-based numerical calculation engine for solving 2-D heat and moisture transfer problems. This numerical engine is used in the THERM software tool, and its primary purpose is for the analysis of building envelopes (e.g., windows, walls, roofs, foundations, etc.). However, the engine can also be used for any heat and moisture transfer problems that require solving fundamental 2-D energy and mass transfer equations. Fluid flow solutions (Navier-Stokes momentum equations) are not included, but the correlations for various convection heat transfer situations are provided, including the translation of complex cavity geometries into those for which correlations are applicable. The calculation engine is written in C++ and includes an API for connecting to third-party tools.

Vidanovic, Dragan [Lawrence Berkeley National Labo↗
DOE OSTI2025

Photochemically Induced Acousto-optics Fluid Simulations

PIAFS is a finite-difference code to solve the compressible Navier-Stokes equations with chemical heating on Cartesian grids. It models chemical reactions of air (oxygen and carbon dioxide) with ozone subject to radiation. It uses a high-order WENO spatial discretization and explicit Runge-Kutta time integration. It is capable of parallel simulations using MPI. The code is written in C/C++.

Oudin, AlbertineN [Lawrence Livermore National Lab↗
DOE OSTI2026

Hydrogen Dispersion Modeling for Development of Smart Distributed Monitoring

Studying hydrogen dispersion is crucial for ensuring the safe and effective deployment of hydrogen as an energy carrier. This study presents a comprehensive CFD modeling framework for simulating hydrogen dispersion at a real-world hydrogen production, storage, and utilization facility. Utilizing the Hydrogen Research Facility under the Advanced Research on Integrated Energy Systems (ARIES) at the National Renewable Energy Laboratory's (NREL) Flatirons campus, controlled hydrogen releases at 27 kg-H2/hr were simulated. The model incorporated site-specific atmospheric conditions, including hourly wind speeds and temperatures recorded between 8 AM and 8 PM from October to December 2023. To reduce computational demands, a statistical reduction technique was applied to condense the dataset to 100 representative scenarios, validated by statistical tests for wind speeds and power law coefficients. Simulations were conducted using the Reynolds-Averaged Navier-Stokes equations. Results demonstrated that wind speed substantially influences hydrogen dispersion, with low wind conditions forming concentrated clouds and higher wind speeds stretching the plume. Additionally, clustering analysis informed optimal sensor placement at various elevations with up to 10 sensor locations on each elevation. This framework offers a robust approach for understanding hydrogen behavior in ambient conditions and informing detection strategies.

08 HYDROGEN↗
NASA NTRS1995

Analysis of developing laminar flows in circular pipes using a higher-order finite-difference technique

A higher-order finite-difference technique is developed to calculate the developing-flow field of steady incompressible laminar flows in the entrance regions of circular pipes. Navier-Stokes equations governing the motion of such a flow field are solved by using this new finite-difference scheme. This new technique can increase the accuracy of the finite-difference approximation, while also providing the option of using unevenly spaced clustered nodes for computation such that relatively fine grids can be adopted for regions with large velocity gradients. The velocity profile at the entrance of the pipe is assumed to be uniform for the computation. The velocity distribution and the surface pressure drop of the developing flow then are calculated and compared to existing experimental measurements reported in the literature. Computational results obtained are found to be in good agreement with existing experimental correlations and therefore, the reliability of the new technique has been successfully tested.

Gladden, Herbert J.↗
NASA NTRS1995

Effect of turbulence models on criticality conditions in swirling flows

The critical state of vortex cores downstream of vortex breakdown has been studied. Base vortical flows were computed using the Reynolds-averaged, axisymmetric Navier-Stokes equations. Standard K - epsilon, RNG and second-order Reynolds stress models were employed. Results indicate that the return to supercriticality is highly dependent on the turbulence model. The K - epsilon model predicted a rapid return of the vortex to supercritical conditions, the location of which showed little sensitivity to changes in the swirl ratio. The Reynolds stress model predicted that the vortex remains subcritical to the end of the domain for each of the swirl ratios employed, and provided results in qualitative agreement with experimental work. The RNG model produced intermediate results, with a downstream movement in the critical location with increasing swirl. Calculations for which area reductions were introduced at the exit in a subcritical flow were also performed using the Reynolds stress model. The structure of the resulting recirculation zone was altered significantly. However, when area reductions were employed within supercritical flows as predicted using the two-equation models, no significant influence on the recirculation zone was noted.

Gatski, Thomas B.↗
NASA NTRS1995

Turbulance boundary conditions for shear flow analysis, using the DTNS flow solver

The effects of different turbulence boundary conditions were examined for two classical flows: a turbulent plane free shear layer and a flat plate turbulent boundary layer with zero pressure gradient. The flow solver used was DTNS, an incompressible Reynolds averaged Navier-Stokes solver with k-epsilon turbulence modeling, developed at the U.S. Navy David Taylor Research Center. Six different combinations of turbulence boundary conditions at the inflow boundary were investigated: In case 1, 'exact' k and epsilon profiles were used; in case 2, the 'exact' k profile was used, and epsilon was extrapolated upstream; in case 3, both k and epsilon were extrapolated; in case 4, the turbulence intensity (I) was 1 percent, and the turbulent viscosity (mu(sub t)) was equal to the laminar viscosity; in case 5, the 'exact' k profile was used and mu(sub t) was equal to the laminar viscosity; in case 6, the I was 1 percent, and epsilon was extrapolated. Comparisons were made with experimental data, direct numerical simulation results, or theoretical predictions as applicable. Results obtained with DTNS showed that turbulence boundary conditions can have significant impacts on the solutions, especially for the free shear layer.

Mizukami, M.↗
NASA NTRSDate not supplied

The NASA Turbulent Heat Flux Experiments: Summary and Lessons Learned

The Turbulent Heat Flux (THX) experiments were conducted at NASA Glenn Research Center (GRC) in order to collect measurements of velocities and temperatures for computational fluid dynamics (CFD) validation of heated flows, with a focus on propulsion system components. The experiments spanned 5 phases; four of which were conducted in the GRC AeroAcoustic Propulsion Laboratory (AAPL) using the Small Hot Jet Flow Rig (SHJAR). In addition to making velocity measurements with Particle Image Velocimetry (PIV), the THX experiments introduced a new Raman-scattering based capability to measure temperatures. Computational studies were also conducted for each of the experimental configurations, in order to provide a baseline of expected CFD results and conduct an assessment of the capability of various CFD approaches for calculating flows where the turbulent transport of heat was important. Two of the collected sets of data were used for American Institute of Aeronautics and Astronautics (AIAA) Propulsion Aerodynamic Workshops (PAWs). The data set from the 5th phase, collected for heated supersonic jets, was also used to construct new validation cases for the NASA Turbulence Model Resource (TMR). This paper provides an overview of the experiments and associated computations for each of the 5 test phases. Key experimental findings are presented. Lessons learned are provided concerning the effect of computational modeling choice on accuracy of predicting turbulent flows where thermal transport is important. Emphasis is placed on comparing Reynolds- averaged Navier-Stokes approaches with large-eddy simulation approaches. The benefits of utilizing a conjugate heat transfer method in conjunction with CFD solver for film cooling is demonstrated.

RANS↗
NASA NTRSDate not supplied

The NASA Turbulent Heat Flux (THX) Experiments: Summary and Lessons Learned

The Turbulent Heat Flux (THX) experiments were conducted at NASA Glenn Research Center (GRC) in order to collect measurements of velocities and temperatures for computational fluid dynamics (CFD) validation of heated flows, with a focus on propulsion system components. The experiments spanned 5 phases; four of which were conducted in the GRC AeroAcoustic Propulsion Laboratory (AAPL) using the Small Hot Jet Flow Rig (SHJAR). In addition to making velocity measurements with Particle Image Velocimetry (PIV), the THX experiments introduced a new Raman-scattering based capability to measure temperatures. Computational studies were also conducted for each of the experimental configurations, in order to provide a baseline of expected CFD results and conduct an assessment of the capability of various CFD approaches for calculating flows where the turbulent transport of heat was important. Two of the collected sets of data were used for American Institute of Aeronautics and Astronautics (AIAA) Propulsion Aerodynamic Workshops (PAWs). The data set from the 5th phase, collected for heated supersonic jets, was also used to construct new validation cases for the NASA Turbulence Model Resource (TMR). This paper provides an overview of the experiments and associated computations for each of the 5 test phases. Key experimental findings are presented. Lessons learned are provided concerning the effect of computational modeling choice on accuracy of predicting turbulent flows where thermal transport is important. Emphasis is placed on comparing Reynolds-averaged Navier-Stokes approaches with large-eddy simulation approaches. The benefits of utilizing a conjugate heat transfer method in conjunction with CFD solver for film cooling is demonstrated.

RANS↗