Search NASA⌕ Search

SEARCH · Search NASA

Results for “fluid model”

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 OSTI2025

Fluid modeling of low-temperature plasmas

Fluid models are essential for understanding and predicting low-temperature plasma (LTP) behavior in various scientific and industrial settings. This paper provides an introductory tutorial on fluid modeling of LTPs, covering model formulation, implementation, and computational simulations. The tutorial focuses on five main components of the formulation of LTP fluid models: fluid flow, energy, chemistry, electromagnetism, and material properties, as well as in essential aspects of model implementations, including multiscale phenomena, multiphysics coupling, and numerical convergence. Designed for students and early-career researchers, this work offers a practical foundation for developing and using fluid models, from in-house computational codes to commercial software, bridging fundamental theory with real-world applications.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Ionization waves in low-current dc discharges in noble gases obtained with a hybrid kinetic-fluid model

A hybrid kinetic-fluid model is used to study ionization waves (striations) in a low-current plasma column of dc discharges in noble gases. Coupled solutions of a kinetic equation for electrons, a drift-diffusion equation for ions, and a Poisson equation for the electric field are obtained to clarify the nature of plasma stratification in the positive column. A simplified two-level excitation-ionization model is used for the conditions when the nonlinear effects due to stepwise ionization, gas heating, and Coulomb interactions among electrons are negligible. It is confirmed that the nonlocal effects are responsible for the formation of moving striations in dc discharges at low plasma densities and low values of pR (the product of gas pressure and tube radius). Here, the calculated properties of self-excited waves of S–, P–, and R types in neon and S type in argon agree with available experimental data. The reason for helium plasma stability to stratification is clarified. It is shown that sustaining stratified plasma is more efficient than striation-free plasma when the ionization rate is a nonlinear function of the electric field. However, the nonlinear dependence of the ionization rate on the electric field is not required for plasma stratification. Striations of S–, P–, and R types in neon exist with minimal or no ionization enhancement. Effects of the column length and plasma density on the wave properties are demonstrated.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

A thermodynamically consistent discretization of 1D thermal-fluid models using their metriplectic 4-bracket structure

Thermodynamically consistent models in continuum physics, i.e. models which satisfy the first and second laws of thermodynamics, may be expressed using the metriplectic formalism. In this work, we leverage the structures underlying this modeling formalism to preserve thermodynamic consistency in discretizations of a fluid model. The procedure relies (1) on ensuring that the spatial semi-discretization retains certain symmetries and degeneracies of the Poisson and metriplectic 4-brackets, and (2) on the use of an appropriate energy conserving time-stepping method. Here, the minimally simple yet nontrivial example of a one-dimensional thermal-fluid model is treated. It is found that preservation of the requisite symmetries and degeneracies of the 4-bracket is relatively simple to ensure in Galerkin spatial discretizations, suggesting a path forward for thermodynamically consistent discretizations of more complex fluid models using more specialized Galerkin methods.

Hamiltonian structure↗
DOE OSTI2026

An implicit-explicit time splitting strategy for the far SOL plasma fluid model with DG-FEM discretization

We consider a far scrape-off layer (SOL) plasma fluid model of ions that is governed by a Braginskiitype model: a one-dimensional, nonlinear system of advection-diffusion equations coupled with a diffusion equation for neutral particles. Our motivation for studying this system arises from the coupling between the edge plasma and radio-frequency (RF) heating, where solving a far SOL plasma fluid model provides critical insights into edge plasma dynamics. Numerical simulations of plasma fluid models require advanced computational techniques to achieve both efficiency and accuracy, especially when resolving the boundary layer in magnetically confined plasmas. In this work, we propose an implicit-explicit time operator splitting strategy that allows for an efficient solution algorithm, where the diffusive terms are treated semi-implicitly requiring only a linear solve, while the advection part is handled explicitly using a strong-stability-preserving Runge-Kutta (SSP-RK3) scheme. This leads to a fully decoupled system in which the diffusion and advection sub-problems can be solved separately, simplifying the overall solution procedure and allowing for efficient parallelization, which is particularly relevant for exploring the impact of RF heating on the SOL plasma. The main challenge of the discretization is due to the strong coupling between diffusion and advection, particularly through the boundary conditions. This makes implementation of such a scheme in an accurate and stable manner nontrivial. We discuss in detail how to split the equations and manage boundary conditions to maintain stability and well-posedness for each subsystem. We also describe a spatial discretization approach, based on the discontinuous Galerkin finite element method (DG-FEM) and present numerical results for a one-dimensional system.

Burkovska, Olena [ORNL] (ORCID:0000000163101130)↗
DOE OSTI2022

Electrostatic Relativistic Fluid Models of Electron Emission in a Warm Diode

A semi-analytic fluid model has been developed for characterizing relativistic electron emission across a warm diode gap. Here we demonstrate the use of this model in (i) verifying multi-fluid codes in modeling compressible relativistic electron flows (the EMPIRE-Fluid code is used as an example; see also Ref. 1), (ii) elucidating key physics mechanisms characterizing the influence of compressibility and relativistic injection speed of the electron flow, and (iii) characterizing the regimes over which a fluid model recovers physically reasonable solutions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Fluid model for the piezothermal effect

When a gas in an externally imposed potential field is compressed, temperature gradients appear. This has been called the piezothermal effect. It is possible to analytically calculate the time-dependent behavior of the piezothermal effect using a linearized fluid model. Quantitative differences between the fluid-model results and previous numerical calculations can be explained by the effects of viscosity and heat conductivity. Here, the fluid model casts the piezothermal effect as a spectrum of buoyancy oscillations, which yields new physical insights into the effect.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Neural network representability of fully ionized plasma fluid model closures

The closure problem in fluid modeling is a well-known challenge to modelers aiming to accurately describe their systems of interest. Over many years, analytic formulations in a wide range of regimes have been presented but a practical, generalized fluid closure for magnetized plasmas remains an elusive goal. Here, as a first step toward constructing a novel data-based approach to this problem, we apply ever-maturing machine learning methods to assess the capability of neural network architectures to reproduce crucial physics inherent in known magnetized plasma closures. We find encouraging results, indicating the applicability of neural networks to closure physics, but also arrive at recommendations on how one should choose appropriate network architectures for the given locality properties dictated by the underlying physics of the plasma.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Theoretical Review of the MFIX Fluid and Two-Fluid Models

MFIX (Multiphase Flow with Interphase eXchanges) is a general purpose code that can be used for describing the hydrodynamics, chemical reactions, and heat transfer of dense or dilute multiphase fluid flows. MFIX calculations give detailed information on pressure, temperature, composition, and velocity distributions in the defined system. This report provides an overview of MFIX’s expressions for single phase flow and continuous dispersed multiphase flow (i.e., the two-fluid model). In addition to the two fluid model (MFIX-TFM), MFIX contains a discrete element model (MFIX-DEM) and a multi-phase particle in cell model (MFIX-PIC) that are not reviewed within this document.

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

Thermo-Fluid Modeling Framework for Supercomputer Digital Twins: Part 1, Demonstration at Exascale

A thermo-fluid modeling framework is being developed for ExaDigiT---an open-source framework for developing comprehensive digital twins of liquid-cooled supercomputers. The work is being conducted in two parts, and discussion is divided into two companion papers. The work documented in this paper focuses on the development of a cooling system library in Dymola for the Frontier supercomputer at Oak Ridge National Laboratory. The second part, outlined in a companion paper, focuses on a templating structure called Auto-CSM for easily creating model-agnostic, physics-based thermo-fluid cooling system models for liquid-cooled supercomputers using a text-based schema. The cooling model is being developed using primarily the open-source Transient Simulation Framework of Reconfigurable Models (TRANSFORM) library. The library follows the templating architecture developed within the TRANSFORM library for modeling subsystems. A full-system validation was performed to validate a very simple model that is integrated with the system controls, and the results are presented herein.

Kumar, Vineet↗
DOE OSTI2021

Final report of activities for the LDRD-express project #223796 titled: “Fluid models of charged species transport: numerical methods with mathematically guaranteed properties”, PI: Ignacio Tomas, Co-PI: John Shadid

This report summarizes the findings and outcomes of the LDRD-express project with title “Fluid models of charged species transport: numerical methods with mathematically guaranteed properties”. The primary motivation of this project was the computational/mathematical exploration of the ideas advanced aiming to improve the state-of-the-art on numerical methods for the one-fluid Euler-Poisson models and gain some understanding on the Euler-Maxwell model. Euler-Poisson and Euler-Maxwell, by themselves are not the most technically relevant PDE plasma-models. However, both of them are elementary building blocks of PDE-models used in actual technical applications and include most (if not all) of their mathematical difficulties. Outside the classical ideal MHD models, rigorous mathematical and numerical understanding of one-fluid models is still a quite undeveloped research area, and the treatment/understanding of boundary conditions is minimal (borderline non-existent) at this point in time. This report focuses primarily on bulk-behaviour of Euler-Poisson’s model, touching boundary conditions only tangentially.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Modeling fluid phase equilibria of carbon dioxide-methanol binary system

Accurate fluid phase equilibria modeling of carbon dioxide-methanol binary system is essential for numerous industrial applications. Prior modeling studies were limited in scope in terms of temperature and pressure ranges. In this study, the phase equilibria behavior of the binary system is modeled with cubic equations-of-state (EoS) including Peng-Robinson (PR) and Soave-Redlich-Kwong with various mixing rules, Predictive Soave-Redlich-Kwong, Cubic Plus Association, and Perturbed-Chain Statistical Associating Fluid Theory. Among them, the classical PR EoS and its variants yield the overall best results in representing the phase behavior at temperature above 330 K. Furthermore, PR incorrectly predicts two liquid phases when the CO 2 mole fraction exceeds ~ 0.4 at temperature lower than 330 K. Raoult’s law with the classical nonrandom two-liquid excess Gibbs energy model is recommended for the low temperature conditions.

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

Benchmark of interfacial area concentration approaches for the two-fluid model in gas-dispersed condensing flow

Accurate predictions of interfacial area concentration (IAC) are required for closure of the condensation phase change term. The interfacial area transport equation (IATE) can dynamically predict interfacial area concentration and is increasingly being used in the benchmarking of phase change flows over static based correlations which require a flow regime map. For higher flow regimes, the two-group two-fluid IATE model can separately model the transport of spherical bubbles and cap/slug bubbles. Recently, a modified formulation of the one-group and two-group two-fluid model with interfacial area transport equation was proposed for gas-dispersed condensing flows. Three datasets have been used for validation consisting of condensing two-phase flow data in the unheated section of a vertical annulus with conditions spanning a range of pressures, inlet subcoolings and mass fluxes for benchmarking the various IAC approaches. For bubbly flows, the coupled void transport - IATE simulation results in good agreement in the predictions of the void fraction and the interfacial area concentration under nearly all sixty conditions. Considering the proposed two-group two-fluid IATE model formulation, the group-2 Nusselt number correlation is investigated to account for the appropriate heat transfer length scale for group-2 bubbles. Validation of the proposed two-group two-fluid model against the new dataset shows good agreement with experimental data with significant improvement over the one-group IATE model for conditions with substantial group-2 void fraction throughout the simulation domain.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2025

Ten-moment fluid modeling of the Weibel instability

We investigate the one-dimensional non-relativistic Weibel instability through the capture of anisotropic pressure tensor dynamics using an implicit 10-moment fluid model that employs the electromagnetic Darwin approximation. The results obtained from the 10-moment model are compared with an implicit particle-in-cell simulation. The linear growth rates obtained from the numerical simulations are in good agreement with the theoretical fluid and kinetic dispersion relations. The fluid dispersion relations are derived using Maxwell’s equations and the Darwin approximation. We also show that the magnetohydrodynamic approximation can be used to model the Weibel instability if one accounts for an anisotropic pressure tensor and unsteady terms in the generalised Ohm’s law. In addition, we develop a preliminary theory for the saturation magnetic field strength of the Weibel instability, showing good agreement with the numerical results.

Kuldinow, D. A. (ORCID:0000000319730196)↗
DOE OSTI2021

nn-PINNs: Non-Newtonian physics-informed neural networks for complex fluid modeling

Time- and rate-dependent material functions in non-Newtonian fluids in response to different deformation fields pose a challenge in integrating different constitutive models into conventional computational fluid dynamic platforms. Considering their relevance in many industrial and natural settings alike, robust data-driven frameworks that enable accurate modeling of these complex fluids are of great interest. The main goal is to solve the coupled Partial Differential Equations (PDEs) consisting of the constitutive equations that relate the shear stress to the deformation and fully capture the behavior of the fluid under various flow protocols with different boundary conditions. Here, in this work, we present non-Newtonian physics-informed neural networks (nn-PINNs) for solving systems of coupled PDEs adopted for complex fluid flow modeling. The proposed nn-PINN method is employed to solve the constitutive models in conjunction with conservation of mass and momentum by benefiting from Automatic Differentiation (AD) in neural networks, hence avoiding the mesh generation step. nn-PINNs are tested for a number of different complex fluids with different constitutive models and for several flow protocols. These include a range of Generalized Newtonian Fluid (GNF) empirical constitutive models, as well as some phenomenological models with memory effects and thixotropic timescales. nn-PINNs are found to obtain the correct solution of complex fluids in spatiotemporal domains with good accuracy compared to the ground truth solution. We also present applications of nn-PINNs for complex fluid modeling problems with unknown boundary conditions on the surface, and show that our approach can successfully recover the velocity and stress fields across the domain, including the boundaries, given some sparse velocity measurements.

42 ENGINEERING↗
DOE OSTI2024

Establishing criteria for the transition from kinetic to fluid modeling in hollow cathode analysis

Hollow cathodes for plasma switch applications are investigated via 2D3V particle-in-cell simulations of the channel and plume region. The kinetic nature of the plasma within the channel is dependent on the thermalization rate of electrons, emitted from the insert. When Coulomb collisions occur at a much greater rate than ionization or excitation collisions, the electron energy distribution function rapidly relaxes to a Maxwellian and the plasma within the channel can be described accurately via a fluid model. In contrast, if inelastic processes are much faster than Coulomb collisions, then the electron energy distribution function in the channel exhibits a notable high-energy tail, and a kinetic treatment is required. This criterion is applied to hollow cathodes from the literature, revealing that a fluid approach is suitable for most electric propulsion applications, whereas a kinetic treatment can be more critical to accurate modeling of plasma switches.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Two-fluid model of rf current condensation in magnetic islands

The stabilization of tearing modes with rf waves is subject to a nonlinear effect, termed rf current condensation, that has the potential to greatly enhance and localize current driven within magnetic islands. Here we extend previous investigations of this effect with a two fluid model that captures the balance of diffusive and thermal equilibration processes within the island. We show that the effective power and resulting strength of the condensation effect can be greatly enhanced by avoiding collisional heat loss to the ions. The relative impact of collisions on the overall power balance within the island depends on the ratio of the characteristic diffusion timescale and the electron–ion equilibration time, rather than the latter alone. Although relative heat loss to ions increases with island size, the heating efficiency does as well. In particular, we show that the latter safely dominates for large deposition profiles, as is typically the case for lower hybrid current drive. This supports the possibility of passive stabilization of neoclassical tearing modes without the precise aiming of the rf waves required for electron cyclotron current drive stabilization.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Data-driven, multi-moment fluid modeling of Landau damping

Deriving governing equations of complex physical systems based on first principles can be quite challenging when there are certain unknown terms and hidden physical mechanisms in the systems. In this work, we apply a deep learning architecture to learn fluid partial differential equations (PDEs) of a plasma system based on the data acquired from a fully kinetic model. Here, the learned multi-moment fluid PDEs are demonstrated to incorporate kinetic effect such as Landau damping. Based on the learned fluid closure, the data-driven, multi-moment fluid modeling can well reproduce all the physical quantities derived from the fully kinetic model. The calculated damping rate of Landau damping is consistent with both the fully kinetic simulation and the linear theory. The data-driven fluid modeling of PDEs for complex physical systems may be applied to improve the fluid closure and reduce the computational cost of multi-scale modeling of global systems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗