Search NASA⌕ Search

SEARCH · Search NASA

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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
NASA NTRS1952

Generalization of Boundary-Layer Momentum-Integral Equations to Three-Dimensional Flows Including Those of Rotating System

The Navier-Stokes equations of motion and the equation of continuity are transformed so as to apply to an orthogonal curvilinear coordinate system rotating with a uniform angular velocity about an arbitrary axis in space. A usual simplification of these equations as consistent with the accepted boundary-layer theory and an integration of these equations through the boundary layer result in boundary-layer momentum-integral equations for three-dimensional flows that are applicable to either rotating or nonrotating fluid boundaries. These equations are simplified and an approximate solution in closed integral form is obtained for a generalized boundary-layer momentum-loss thickness and flow deflection at the wall in the turbulent case. A numerical evaluation of this solution carried out for data obtained in a curving nonrotating duct shows a fair quantitative agreement with the measures values. The form in which the equations are presented is readily adaptable to cases of steady, three-dimensional, incompressible boundary-layer flow like that over curved ducts or yawed wings; and it also may be used to describe the boundary-layer flow over various rotating surfaces, thus applying to turbomachinery, propellers, and helicopter blades.

Mager, Arthur↗
NASA NTRS1947

Effects of sweepback on boundary layer and separation

Following a law of stress adopted in the Navier-Stokes equations, the configuration of the viscous flow in planes at right angles to the axis of an infinite cylinder is found to be independent of the axial motion of the cylinder. In the limiting case of a yawed or swept wing of very high aspect ratio, certain boundary-layer and separation phenomena are thus determined independently by the crosswise component of velocity. It follows that the effect of sweepback is to increase the area of stable laminar flow and to decrease the lift coefficient at which flow separation occurs.

Jones, Robert T↗
NASA NTRS1951

Sideslip in a viscous compressible gas

By an analysis of the Navier-Stokes equations it is shown that the aerodynamic coefficients of an infinite rectangular swept wing in an isothermal or adiabatic flow of a compressible gas can be determined from the aerodynamic coefficients of the unswept wing. When the flow is neither isothermal nor adiabatic, a three-dimensional boundary layer theory is developed and applied to the special case of a swept flat plate.

THEORIES - WINGS, SWEEPBACK & SWEEPFORWARD↗
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 NTRS1961

Analysis of Multipoint-Multitime Correlations and Diffusion in Decaying Homogeneous Turbulence

Two-point, two-time correlation equations are obtained by considering the Navier-Stokes equations for two points in a fluid at two time. By neglecting the triple correlations in the equations, a solution is obtained for the final period of decay. The analysis is extended to earlier times by considering three points at three different times. The set of equations is made determinate by neglecting the quadruple correlations in comparison with the triple correlations. The diffusion of particles from a source in a decaying turbulent field is calculated approximately by assuming that the velocity fluctuations are small.

Deissler, Robert G.↗
NASA NTRS1961

Applicability of Mixing Length Theory to a Turbulent Vortex System

The ability of mixing length theory to correlate vortex data is evaluated. Expressions are derived for eddy diffusivity by applying the techniques of von Karman and Prandtl which have been established for pipe flow. Total and static pressures were measured from the outer radius to the exhaust-nozzle radius of a vortex generator for a range of mass flows. These data are combined with Navier-Stokes solutions for this region of a compressible vortex to determine turbulent Reynolds numbers. The Reynolds number is related to Prandtl and Karman functions for various assumed boundary conditions, and the experimental data are used to determine the usefulness of these expressions. The following conclusions were reached: (1) Mixing length functions developed by applying von Karman's similarity hypothesis to vortex motion correlate the data better than do Prandtl functions obtained with the assumption that mixing length is proportional to radius. (2) Some of the expressions developed do not adequately represent the experimental data. (3) The data are correlated with acceptable scatter by evaluating the fluid radial inertia at the outer boundary and the shear stress at the inner boundary. The universal constant K was found to be 0.04 to 0.08, rather than the value of 0.4 which is accepted for rectilinear flow. (4) The data are best correlated by a modified Karman expression which includes an effect of radial inertia, as well as shear stress, on eddy diffusivity.

Ragsdale, Robert G.↗
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↗
NASA NTRSDate not supplied

Improvements to RANS Modeling for Aeroheating Predictions on Blunt Bodies

Accurate predictions of aeroheating are critical for designing thermal protection systems for planetary entry vehicles. For larger vehicles, turbulence in the boundary layer can substantially increase convective heating. This turbulence must be accurately modeled to ensure the thermal protection system is sufficient. The majority of hypersonic turbulence model development and validation focuses on boundary layers developing over flat-plates or sharp cones; these cases are substantially different than the boundary layer that develops over the heatshield of a blunt body traveling at hypersonic speeds. Planetary missions often use blunt body geometries, such as the 70-degree sphere-cone favored by Mars missions or the 45-degree sphere-cone planned for the upcoming DAVINCI mission. Due to smaller vehicle size and the lower velocities in the stagnation region, planetary entry vehicles have relatively low Reynolds numbers. Surface curvature and high enthalpy gradients create additional challenges. These difficulties must be addressed to obtain high accuracy needed for the ambitious planetary missions in the upcoming decade. This work focuses on both assessing and improving Reynolds-averaged Navier-Stokes (RANS) turbulence models for blunt-body geometries typical of planetary entry vehicles, with a focus on one-equation and two-equation formulations.

Mars2020↗
DOE OSTI2024

Implementation of Surface Tension on a Reacting Flow Solver, PeleLM: Preprint

In liquid rocket engines, the fuel is supplied to the combustion chamber in the liquid state though injectors. Such fuel undergoes atomization, vaporization, and combustion processes. To design reliable and efficient injectors, it is required to understand the full processes. This research is part of an effort to develop a full atomization-vaporization-combustion solver from first principles. As an initial step to tackle the atomization process, a multiphase flow solver is under development. For the development, a library of the volume of fluid scheme for multiphase, IRL is coupled with a reacting Navier-Stokes equation solver, PeleLM. Furthermore, as the surface tension has considerable effects on spray breakup. surface tension is implemented in the momentum equation using the continuum surface force model and the improved height function technique.

height function↗
DOE OSTI2025

Computational Modeling of Molten Salt Infiltration and Oxidation in Nuclear Graphite

Graphite is utilized as a moderator and reflector in advanced nuclear reactor designs due to its high thermal conductivity, neutron moderation properties, and resistance to radiation damage. However, its longterm performance and reliability are challenged by degradation mechanisms such as molten salt infiltration in molten salt reactors (MSRs) and oxidation in gas-cooled reactors (GCRs). These mechanisms can compromise the structural integrity and operational lifetime of graphite components, necessitating a more detailed assessment of their physical behavior. This report focuses on the development of computational models for molten salt infiltration and oxidation of graphite to aid the design and performance analysis of graphite components. For molten salt infiltration, a computational framework is developed that couples incompressible Navier-Stokes and phase-field model to simulate the penetration of molten salt into graphite?s interconnected pore structure. Initial model verification is performed using two-phase flows in two dimensions, demonstrating the models ability to capture fundamental physical behavior and agree with analytical solution. This framework is then applied to a realistic IG110 nuclear graphite , where a computed tomography extracted pore geometry is used to analyse the infiltration behavior of FLiNaK molten salt. This model provides insights into how the microstructure and other relevant parameters influence the transport pathways of molten salt into graphite, potentially offering a means to rapidly evaluate a graphite grade?s resistance to infiltration. For oxidation, the report details pore-scale mass and heat transport models, describing the diffusion of gases, reaction kinetics, and thermal effects. Additionally, this report highlights inconsistencies in the existing volume-averaged macroscopic model, particularly in upscaling of reaction kinetics and flux terms, and surface to volume transformations. These inconsistencies suggest that current formulations may not accurately capture the experimentally observed graphite oxidation process, highlighting the need for improved model development. This work advances the development of physics-based computational models for graphite degradation, contributing to improved predictive models for next-generation nuclear reactor designs. Future efforts will focus on refining the infiltration model to address non-physical behaviors and enhance its robustness. Additionally, for oxidation, further studies will employ the principles of volume averaging to rigorously derive the upscaled equations, potentially in collaboration with subject matter experts.

Computational Modeling of Molten Salt Infiltration↗
DOE OSTI2025

Summary Report Of The FY25 Computational Fluid Dynamics Verification And Validation Exercises In The Advanced Reactor Technologies - Gas-cooled Reactor Program

Verification and Validation (V&V) of numerical tools is critical for ensuring reasonable predictions during design, safety analysis, and licensing. Recent work in the Advanced Reactor Technologies - Gas-cooled Reactor (ART-GCR) program has focused on V&V of common Computational Fluid Dynamics (CFD) tools that are used within the Untied States. This report presents an update on these CFD V&V activities. These Generation IV Forum (GIF) Very High Temperature Reactor (VHTR) Computational Methods, Validation, and Benchmarking (CMVB) is an international organization dedicated to the verification and validation of High Temperature Gas-Cooled Reactor (HTGR) simulation tools. Participation in the CMVB provides additional value to the V&V activities, as it allows for access to a wider range of data, and provides valuable benchmarking exercises. Three HTGR phenomena are targeted: Reactor Cavity Cooling System (RCCS) performance, core bypass flow, and lower plenum mixing. Simulations of the University of Wisconsin-Madison (UW-Madison) RCCS facilities are performed with Reynolds Averaged Navier-Stokes (RANS) in StarCCM+. Results are compared for both forced and natural convection conditions, with both exhibiting good agreement with experimental measurements. The Idaho National Laboratory (INL) matched index of refraction (MIR) and Korean Atomic Energy Research Institute (KAERI) bypass flow expeirments are used to validation CFD predictions of bypass flow. Simulations are performed with RANS in StarCCM+ and with Large Eddy Simulation (LES) in NekRS. Finally, preliminary simulations of the Institute of Nuclear and New Energy Technology (INET) lower plenum mixing facilities are presented. Initial work has developed models with LES, RANS, and porous media models. These preliminary models are presented and compared to each other to gauge differences in predictions with each of the three methods.

and Benchmarking (CMVB)↗
DOE OSTI2025

Towards stabilization of pressurized oxy-coal diffusion flames

Designing an effective burner is vital for the development of pressurized oxy-fuel combustion technologies. Turbulent jet diffusion burners are adopted for the pressurized oxy-combustor in this research and a bluff-body is employed to create a recirculation zone to stabilize the flame in pressurized oxy-combustor. The objective of this numerical work is to perform a systematic analysis of the characteristics of this pressurized non-premixed flame. Specifically, a 15-bar POFC combustor of power 100 kWth is modeled employing the Ansys FLUENT commercial platform, using Reynolds-averaged Navier-Stokes (RANS) modeling. The present computational work focuses on identifying the aerothermodynamic features of the isothermal and reacting flow with a disk-shape bluff-body. Results show that fuel-oxidizer stream momentum ratio has great impacts on the temperature profile of the down-fired, co-axial, pressurized oxy-coal, diffusion flame. Parametric study of the width of the burner tip identifies an optimal range for this pressurized burner.

Li, Lei↗
DOE OSTI2025

Recent Improvements in Pronghorn for Advanced Reactor Modeling

Pronghorn is a thermal-hydraulics computational tool developed using the Idaho National Laboratory's Multiphysics Object-Oriented Simulation Environment (MOOSE). It is designed to support Computational Fluid Dynamics (CFD) modeling, ranging from subchannel and porous media analysis to Reynolds Averaged Navier-Stokes (RANS) turbulence modeling. As an integral part of the MOOSE-based suite of tools, Pronghorn seamlessly couples with other MOOSE-based applications to simulate a variety of physical phenomena. This article highlights recent significant enhancements to Pronghorn's CFD modeling capabilities and demonstrates their application to advanced nuclear reactor designs. The recent improvements in Pronghorn primarily focus on modifications to its turbulence modeling capabilities, near-wall corrections and numerical schemes. In terms of turbulence modeling, the two-equation $k-\epsilon$ and $k-\omega$ SST models have been implemented and validated with both equilibrium and non-equilibrium wall treatments. Additionally, corrections for wall roughness, and curvature, and wall-channeling in pebble beds have been introduced in the near-wall modeling. These developments enable more accurate simulations of advanced nuclear reactors. Two case studies are presented in this work: a pool-type Molten Chloride Reactor and a salt-cooled Pebble-Bed High Temperature Reactor. In both cases, the previous models in Pronghorn are compared with the new implementations, demonstrating the improved accuracy achieved with the updated models.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2025

Amplification of turbulence through multiple planar shocks

We study the amplification of isotropic, incompressible turbulence through multiple planar, collisional shocks, using analytical linear theory. There are two limiting cases we explore. The first assumes shocks occur rapidly in time such that the turbulence does not evolve between shocks. Whereas the second case allows enough time for turbulence to isotropize between each shock. For the latter case, through a quasi-equation-of-state, we show that the weak multishock limit is agnostic to the distinction between thermal and vortical turbulent pressures, like an isotropic volumetric compression. When turbulence does not return to isotropy between shocks, the generated anisotropy—itself a function of shock strength—can feedback on amplification by further shocks, altering choices for maximal or minimal amplification. In addition for this case, we find that amplification is sensitive to the shock ordering. As a result, we map how choices of shock strength can impact these amplification differences due to ordering, finding, for example, shock pairs which lead to identical mean postshock fields (density, temperature, pressure) but maximally distinct turbulent amplification.

Astrophysical fluid dynamics↗
DOE OSTI2025

Thermal relaxation and the complete set of second-order transport coefficients for the unitary Fermi gas from kinetic theory

We compute the complete set of second-order transport coefficients of the unitary Fermi gas, a dilute gas of spin-1/2 particles interacting via an 𝑠 -wave interaction tuned to infinite scattering length. The calculation is based on kinetic theory and the Chapman-Enskog method at second order in the Knudsen expansion. We take into account the exact two-body collision integral. We extend previous results on second-order coefficients related to shear stress by including terms related to heat flow and gradients of the fugacity. We confirm that the thermal relaxation time is given by the simple estimate 𝜏 𝜅 = 𝜅⁢𝑚/(𝑐 𝑃 ⁢𝑇) even if the full collision kernel is taken into account. Furthermore, 𝜅 is the thermal conductivity, 𝑚 is the mass of the particles, 𝑐𝑃 is the specific heat at constant pressure, and 𝑇 is the temperature.

Kinetic theory↗