The structure of strong shock waves in the krook collision model
Shock wave profile computation using krook collision model equation in iteration scheme starting from navier-stokes solution
KROOK EQUATION↗
SEARCH · Search NASA
Results for “Navier–Stokes equations”
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.
On a class of exact solutions of the equations of motion of a viscous fluid
The general solution is obtained of the equations of motion of a viscous fluid in which the velocity field is inversely proportional to the distance from a certain point. Some particular cases of such motion are investigated.
FLOW, LAMINAR↗
Asymptotic consistency of the WSINDy algorithm in the limit of continuum data
In this work we study the asymptotic consistency of the weak-form sparse identification of nonlinear dynamics algorithm (WSINDy) in the identification of differential equations from noisy samples of solutions. We prove that the WSINDy estimator is unconditionally asymptotically consistent for a wide class of models that includes the Navier–Stokes, Kuramoto–Sivashinsky and Sine–Gordon equations. We thus provide a mathematically rigorous explanation for the observed robustness to noise of weak-form equation learning. Conversely, we also show that, in general, the WSINDy estimator is only conditionally asymptotically consistent, yielding discovery of spurious terms with probability one if the noise level exceeds a critical threshold σ c . We provide explicit bounds on σ c in the case of Gaussian white noise and we explicitly characterize the spurious terms that arise in the case of trigonometric and/or polynomial libraries. Furthermore, we show that, if the data is suitably denoised (a simple moving average filter is sufficient), then asymptotic consistency is recovered for models with locally-Lipschitz, polynomial-growth nonlinearities. Our results reveal important aspects of weak-form equation learning, which may be used to improve future algorithms. We demonstrate our findings numerically using the Lorenz system, the cubic oscillator, a viscous Burgers-growth model and a Kuramoto–Sivashinsky-type high-order PDE.
asymptotic consistency↗
Building a new multiphysics workflow in MOOSE: application to tritium migration, trapping and advection in TMAP8
Fusion devices are anticipated to produce and consume several kilograms of tritium per year. This rare fuel resource is both highly mobile and radioactive, making tracking inventories a priority for operation and safety. The fusion safety program at the Idaho National Laboratory has been developing the Tritium Migration and Analysis Program (TMAP), of which the latest version is a MOOSE-based application. TMAP8 is verified against its predecessors and possesses additional multi-dimensional tritium migration modeling capabilities. As we extend its capabilities towards both whole device (in multiple dimensions) and whole plant (with multiple components) simulations, the syntax of inputs must become compact, descriptive, compatible with quality assurance processes, and as error-proof as achievable. The new Physics system developed MOOSE can set up equations and instantiating them on plant components. The system permits the automatic definition of complex discretization with a consistency between object parameters achieved programmatically. The Physics system can currently instantiate the equations for heat conduction and Navier Stokes weakly compressible flow. In MOOSE-terms, it automates the definition of kernels, boundary conditions, and several core and helper materials and fields. As part of this effort, Physics classes were developed for tritium migration, trapping and advection within either a multi-dimensional Navier Stokes fluid dynamics simulation, or a 1D thermal hydraulics piping system. In this presentation, we will showcase the new syntax, its application to several verification and validation cases which were already studied using the classical TMAP8 syntax, and a demonstration of the new coupling capabilities for the migration of tritium into blanket coolant channels and the subsequent advection into the coolant loop.
70 - PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Hybrid combustion modeling approach for turbulent jet ignition in natural-gas pre-chamber spark-ignition engines at high EGR
Here, this study presented a hybrid modeling approach for simulating turbulent jet ignition and combustion processes in a natural-gas pre-chamber spark-ignition engine operating under exhaust gas recirculation (EGR) diluted conditions. In-depth analyses of experimental data and simulation results from previous work [Chinnathambi et al., ICEF2021-67836; Kim et al., Fuel 409: 137815, 2026] revealed two key findings: (i) the magnitude of pressure difference between the pre-chamber and main chamber ($∆P_{PC-MC}$) was positively correlated with the combustion duration from the moment of $∆P_{PC-MC}=0$ to the point of 5% mass fraction burned, with larger $∆P_{PC-MC}$ associated with longer duration; and (ii) the turbulent combustion regime in the main chamber transitioned from the broken reaction zone to the corrugated flamelet regime, with the Karlovitz number exceeding 100 immediately after turbulent hot jets were ejected from nozzles, coinciding with observed local extinction events. To accurately simulate the entire combustion process, a hybrid approach was developed under Reynolds-Averaged Navier Stokes framework, combining the G-equation model for pre-chamber combustion with the multi-zone well-stirred reactor approach and a turbulence-chemistry interaction (TCI) submodel for main chamber combustion. The TCI submodel accounted for the attenuation of reaction rates due to turbulent strain and modeled local extinction by suppressing reaction rates under certain flow and flame conditions. When applied to three EGR rate conditions toward the dilution limit, the hybrid modeling approach accurately reproduced experimental data in terms of cylinder pressure, apparent heat release rate, and the observed positive correlation, including the delayed onset of main chamber combustion—a feature not captured by existing combustion models.
computational fluid dynamics simulation↗
Potential quantum advantage for simulation of fluid dynamics
Numerical simulation of turbulent fluid dynamics needs to either parametrize turbulence—which introduces large uncertainties—or explicitly resolve the smallest scales—which is prohibitively expensive. Here, we provide evidence through analytic bounds and numerical studies that a potential quantum speedup can be achieved to simulate fluid dynamics using quantum computing. Specifically, we provide a lattice Boltzmann formulation of fluid dynamics for which we give evidence that low-order Carleman linearization is much more accurate than previously believed for these systems. This is achieved via a combination of reformulating the Navier-Stokes nonlinearity (u·$\triangledown$u) to lattice-Boltzmann nonlinearity (u 2 ) and accurately linearizing the dynamical equations, which effectively trades nonlinearity for additional degrees of freedom that add negligible expense in the quantum solver. Based on this, we apply a quantum algorithm for simulating the Carleman-linearized lattice Boltzmann equation and provide evidence that its cost scales logarithmically with system size compared with polynomial scaling in the best known classical algorithms. In this paper, we suggest that a quantum advantage may exist for simulating fluid dynamics, paving the way for simulating nonlinear multiscale transport phenomena in a wide range of disciplines using quantum computing.
42 ENGINEERING↗
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↗
Nonlinear simulation of under-resolved flows with shocks
Here, we consider the numerical simulation of advection-dominated flows whose wide range of physical length scales exceed the memory capacity of finite computers. Simulating flows with shocks and turbulence presented challenges for the earliest computers that were quickly overcome by the development of new numerical methodology. Principal among those new ideas were artificial viscosity and finite volume methods, concepts that remain in common use today. We begin by describing the history of those methods, the innovators and their motivations. We then describe the development of finite scale theory, a reformulation of Navier–Stokes theory that exposes the physical principles on which artificial viscosity is based. We discuss the essential properties of the finite scale equations, the observer, unresolved kinetic energy and inviscid energy dissipation. We briefly consider the implementation of the finite scale equations on the computer from the point of view of Gisin’s conjectures about finite information.
97 MATHEMATICS AND COMPUTING↗
Modeling and Simulation of Fuel Dispersal During the Loss-of-Coolant Accident
This document is the compilation of the milestone portion to a larger end of project NEUP report. The executive summary of the modeling portion is provided below: In the event of cladding rupture during a postulated LOCA in a pressurized water reactor, fuel particles, along with fission gases, can be expelled into the reactor core from the fractured fuel rod, a phenomenon referred to as fuel dispersal. The initial stage of fuel dispersal is strongly influenced by the high-pressure ejection of fuel fragments, the size and geometry of the ruptured cladding, and the depressurization history of the fuel rod during the postulated LOCA transient. Depending on the location of the burst orifice relative to the quench front, the dispersal event represents an intricate three-phase flow and heat transfer phenomenon, where high-temperature fuel particles carried by the fission gases interact with the coolant within the narrow subchannels of the fuel assemblies, inducing localized phase change. Given the unique multiphysics nature of this phenomena, the current study develops a dedicated computational framework to predict the mass distribution and cooling of dispersing fuel particles, facilitating post-accident assessment and management of the fuel assemblies. Considering the scale of nuclear reactor applications, a continuum three-fluid model is proposed for simulating the transport of solids within the reactor core. With high-temperature fuel fragments within the liquid media, nucleation sites inducing phase changes are dispersed within the flow domain. Coupled with the fact that the transient dispersal event occurs on different time scales than other three-phase flow applications, this study derives a time-averaged three-fluid flow model without losing generality. The assumptions regarding the continuum treatment of the solid phase and the modeling of fuel dispersal behavior are incorporated to simplify the governing equations and derive applicable closure relations. The computational validation of the model was conducted using adiabatic experimental results obtained from ongoing research at Oregon State University, focusing on characterizing fuel dispersal behavior during simulated LOCA conditions. Settlement characteristics of the solids, quantified by the probability distribution of equivalent particles, closely matched the probability density functions reported in experimental studies. The transport of fuel particles within a scaled 5 × 5 lattice of a pressurized-water reactor rod bundle geometry was modeled through a two-fluid Eulerian framework. The required boundary conditions were evaluated from the fuel performance code BISON in a postulated large-break LOCA scenario. The modeling framework considered solid fuel particles as granular matter, interacting with the gaseous dry steam phase and fission gases through the governing interfacial momentum exchange between the participating fluids. The simulation results provided the volume fraction of the solids obtained at the bottom surface of the enclosing tank geometry. Postulated LOCA leading to fuel dispersal phenomena involves the strong coupling between fuel thermomechanics, cladding deformation, thermal-hydraulics, and fuel particle transport. Incorporation of such a strong coupling in numerical simulation is performed by coupling the multiphysics solvers. In the case of fuel dispersal, a strong coupled simulation can be performed by coupling the BISON code for fuel performance, the TRACE code for system-level thermal hydraulics, and fuel particle transport in Multiphysics Object-Oriented Simulation Environment (MOOSE). For such intricate infrastructure, the MOOSE Framework eases the data transfer between codes. The recent version of MOOSE has incorporated the Navier-Stokes module for the fluid flow. An exploratory exercise was done to gain familiarity with finite volume capabilities in the MOOSE framework to incorporate the Spalart-Allmaras (SA) turbulence model. New finite-volume and auxiliary kernels were introduced to assemble the SA transport equation, compute turbulent viscosity, and evaluate wall distance and diagnostic turbulence terms, fully integrated with existing Navier-Stokes modules. A turbulent lid-driven cavity at a Reynolds number of approximately 10,000 is used for verification. MOOSE shows the robust solver convergence and produces the turbulent features. But it underpredicts the velocity profile and turbulent quantities, emphasizing the need to develop improved SA near-wall treatments (e.g., low-Re corrections or wall functions) as a key direction for future work.
11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
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.↗
High y + Shear-Stress Turbulence Implementation for High Flux Isotope Reactor Narrow Channel Flows
The research objective of this work was to improve the engineering predictions of the turbulence characteristics of flows in curved narrow channels. Such channel flows are commonly encountered in nuclear research and test reactors, with one of them being the high-flux isotope reactor (HFIR). Research reactors bear high heat fluxes, and the proper computing of turbulence is paramount for safe and reliable reactor operation. The study builds on the results of a previous direct numerical simulation of turbulence to inform a well-known Reynolds-averaged Navier–Stokes shear-stress turbulence model and improves its accuracy in simulating parallel channel flows. A new formulation of the loss term in the dissipation conservation equation is suggested. Combined with high wall distance computational grids, the new implementation provides a fast-running flow solution, suitable for engineering purposes. Model generalization for parallel channel flows, in a broader range of frictional Reynolds numbers, is suggested by introducing a new form of the model constants.
CFD↗
Applying a Compact Porous Media Model to Numerically Derive Resistance Coefficients for Lattice Structures
Additive Manufacturing allows for exploring various geometries to achieve specific engineering criteria. Lattices are one geometry with unique properties, including being periodically repeating structures which allow flow through them to be represented as a porous media according to Darcy-Forchheimer equations. These equation’s coefficients are generally experimentally derived, but this work demonstrates the ability to numerically derive them with CFD. Simulations were performed using three-dimensional stead state Reynolds-averaged Navier-Stokes with a k-ω Shear Stress Transport turbulence model using Ansys Fluent. Three lattice geometries were investigated and drag coefficients were derived. The method was validated against externally published data for similar geometries demonstrating strong agreement, and grid convergence for all simulations was calculated with a Grid Convergence Index method. Wall roughness is demonstrated to have a non-negligible impact on results and roughness values are considered for the primary focus Octahedral geometry where both smooth wall and rough wall coefficients were derived. The porosity coefficients for the Octahedral geometry at 1.0 [m/s] were found to be 2.89×10 6 and 2.90×10 6 [1/(Pa*m*s)] for the permeability coefficients, 6.37×10 1 and 5.44×10 1 [m 2 /kg] for the inertial resistance coefficients, and with a max pressure drop of 5116.7 [Pa] and 4429.5 [Pa] for the smooth walls and rough walls, respectively. The derived numerical method enables rapid exploration and optimization of new lattice designs for diverse engineering applications.
42 ENGINEERING↗
Scalar flux transport models for self-similar turbulent mixing
A common approach to closing turbulent species flux in multicomponent Reynolds-averaged Navier-Stokes models is to use the standard gradient diffusion approximation. While such an approach has been shown to work well when applied to many canonical turbulent mixing configurations, a gradient diffusion approach is fundamentally limited in its ability to capture complex phenomena such as countergradient transport. For this reason, complicated mixing applications may benefit by treating the turbulent diffusivity with a model transport equation in a manner analogous to second-moment momentum closure in Reynolds-stress transport models. Here, the present work explores the development and application of two different scalar flux transport (SFT) models. Self-similarity constraints are derived for these models, and they are evaluated against gradient-diffusion-based models in several one- and two-dimensional problems of turbulent mixing. It is found that the new SFT models out-perform gradient diffusion models in problems involving rapid acceleration reversal and in problems involving anisotropic transport of materials. In addition, it is found that even a hybrid-SFT approach, in which an SFT equation is utilized along with a gradient diffusion closure, provides some measure of improvement over models that transport the mass flux rather than the scalar flux.
Reynolds-averaged Navier Stokes↗
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↗
A robust spectral element implementation of the $k - τ$ RANS model in Nek5000/NekRS
The $k - ω$ Reynolds Averaged Navier Stokes (RANS) model is one of the industry standard approaches for modeling of turbulent flows. It performs better than the $k - ϵ$ model for low Reynolds number flows and is also more suitable for boundary layers with adverse pressure gradients. Major drawback of the model, however, is that the asymptotic value of $ω$ at the walls is singular, necessitating the use of a contrived “sufficiently” large value for $ω$ as the boundary condition for its transport equation. Here, this invariably leads to the solution being sensitive to near wall grid spacing. While an acceptable solution for low order (finite volume) methods, the excessive near wall gradients lead to persistent numerical stability issues in high order codes. To alleviate the problem, specifically in the context of the high order spectral element code Nek5000, a regularized $k - ω$ approach was formulated in our prior work (Tomboulides et al., 2018). The formulation, however, relies on the use of wall distance and its gradients for modeling the closure terms and can pose problems for simulations in complex geometries. This work presents a novel implementation of the $k - τ$ RANS model in Nek5000, where $τ = 1/ω$, eliminating the need for regularization, owing to the asymptotically bounded behavior of the source terms in the $τ$ transport equation, and also eliminating dependence on wall distance. Robustness and stability of the $k - τ$ model is ensured through implicit treatment of the source terms and their careful numerical implementation and demonstrated through several cases aimed at verification and validation. Studies include both canonical and engineering relevant problems, viz., turbulent channel flow, pipe flow, backward facing step, flow over NACA 0012 airfoil and flow in a T-junction. Results from the $k - τ$ model are shown to be consistent with regularized $k - ω$ model and also with the $k - ω$ SST model in OpenFOAM (for select studies). Comparison with experimental data is also shown, where available, to bolster validation efforts for the $k - τ$ model implementation through prediction of key turbulent quantities of interest.
Nek5000↗
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↗
Development and assessment of models for turbulent Rayleigh-Taylor mixing using the macroscopic forcing method
Reynolds-Averaged Navier Stokes (RANS) simulations are a popular method for designing ICF experiments, and accurate mixing models are crucial for these simulations to give good predictions. To this end, the present work seeks to demonstrate the Macroscopic Forcing Method (MFM) as a tool for both improving existing RANS models as well as assessing RANS model forms. First, MFM analysis from Lavacot et al. (Phys. Rev. Fluids, 2025) is used to develop the k–L–F model, an extension of the k–L model of Dimonte and Tipton (Phys. Fluids, 2006) that incorporates nonlocality through addition of a turbulent species flux transport equation. MFM is then applied to the k–L–F model along with the k–L and BHR–4 models to assess their forms and compare the model-implied eddy diffusivity moments to those measured from high-fidelity simulations. Furthermore, the analysis reveals that models incorporating nonlocality (k–L–F and BHR–4) match the high-fidelity simulation data better than purely local models (k–L), both in terms of mean fields and eddy diffusivity moments. However, all of the considered RANS models struggle to match temporal moments at high Atwood numbers, highlighting the importance of temporal nonlocality in these regimes and the need for additional improvement even among models incorporating nonlocality.
general physics↗
Type-II spacetimes and the double copy for fluids metrics
In our previous paper [Keeler , From Navier-Stokes to Maxwell via Einstein, .] we discussed type-D and type-N fluid-dual spacetimes and provided their associated single copies in the context of the Weyl double copy. In this work we extend our analysis to more general fluids thereby requiring the application of the double-copy picture to type-II spacetimes. By combining our previous type-D and type-N fluids via their associated stream functions we construct an example of a viable type-II double copy. We show that the gravity duals of these fluids perturbatively satisfy Einstein’s equations. We use a near-horizon expansion to identify the type-II double copy for the fluid-dual spacetimes. We show a Maxwell spinor ansatz containing a heterogeneous bispinor component is necessary to provide a viable type-II double copy at the lowest order.
Fluid-gravity correspondence↗