Electrohydrodynamic charge relaxation and interfacial perpendicular-field instability.
Small amplitude motions of plane interface between fluids stressed by initially perpendicular electric field, modeling fluids as ohmic conductors
Melcher, J. R.↗
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.
Uncertainty Quantification for Multiphase Computational Fluid Dynamics Closure Relations with a Physics-Informed Bayesian Approach
Multiphase Computational Fluid Dynamics (MCFD) based on the two-fluid model is considered a promising tool to model complex two-phase flow systems. MCFD simulation can predict local flow features without resolving interfacial information. As a result, the MCFD solver relies on closure relations to describe the interaction between the two phases. Those empirical or semi-mechanistic closure relations constitute a major source of uncertainty for MCFD predictions. In this paper, we leverage a physics-informed uncertainty quantification (UQ) approach to inversely quantify the closure relations’ model form uncertainty in a physically consistent manner. This proposed approach considers the model form uncertainty terms as stochastic fields that are additive to the closure relation outputs. Combining dimensionality reduction and Gaussian processes, the posterior distribution of the stochastic fields can be effectively quantified within the Bayesian framework with the support of experimental measurements. As this UQ approach is fully integrated into the MCFD solving process, the physical constraints of the system can be naturally preserved in the UQ results. Here, in a case study of adiabatic bubbly flow, we demonstrate that this UQ approach can quantify the model form uncertainty of the MCFD interfacial force closure relations, thus effectively improving the simulation results with relatively sparse data support.
42 ENGINEERING↗
An improved ten-moment closure for reconnection and instabilities
The integration of kinetic effects in fluid models is important for global simulations of Earth's magnetosphere. The use of the two-fluid ten-moment model, which includes the pressure tensor for both species, has had some success in simulating Ganymede and Mercury with a simple closure model. Here, we discuss a heat flux closure which accounts for some limitations of the earlier work while remaining computationally tractable. Comparisons with kinetic simulations for magnetic reconnection and lower-hybrid drift instabilities show good agreement with kinetic results and improvements on previous closure models.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Model of a Hollow Cathode Insert Plasma
A 2-D axisymmetric fluid model of the plasma in the insert region of a hollow cathode is presented. The level of sophistication included in the model is motivated in part by the need to determine quantitatively plasma fluxes to the emitter surface. The ultimate goal is to assess whether plasma effects can degrade the life of impregnated inserts beyond those documented throughout the 30-50 year history of vacuum cathode technologies. Results from simulations of a 1.2-cm diameter cathode operating at a discharge current of 25 A, and a gas flow rate of 5 sccm, suggest that approximately 10 A of electron current, and 3.5 A of ion current return to the emitter surface. The total emitted electron current computed by the model is about 35 A. Comparisons with plasma measurements suggest that anomalous heating of the plasma due to two-stream instabilities is possible near the orifice region. Solution to the heavy species energy equation, with classical transport and no viscous effects, predicts heavy species temperatures as high as 2640 K.
hollow cathodes↗
Vidyut3d: A GPU accelerated fluid solver for non-equilibrium plasmas on adaptive grids
We present the numerical methods, programming methodology, verification, and performance assessment of a non-equilibrium plasma fluid solver that can effectively utilize current and upcoming central processing and graphics processing unit (CPU+GPU) architectures, in this work. Our plasma fluid model solves the coupled conservation equations for species transport, electrostatic Poisson and electron temperature on adaptive Cartesian grids. Our solver is written using performance portable adaptive-grid/particle management library, AMReX, and is portable over widely available vendor specific GPU architectures. We present verification of our solver using method of manufactured solutions that indicate formal second order accuracy with central diffusion and fifth-order weighted-essentially-non-oscillatory (WENO) advection scheme. We also verify our solver with published literature on capacitive discharges and atmospheric pressure streamer propagation. We demonstrate the use of our solver on two 3D simulation cases: an atmospheric streamer propagation in Ar-H2 mixtures and a low pressure three-electrode radio frequency reactor. Our performance studies on three different CPU+GPU architectures indicate ~ 150-400X speed-up using AMD and NVIDIA GPUs per time step compared to a single CPU core for a 4 million cell simulation with 15 species.
71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
Development of a Performance Portable Non-Equilibrium Plasma Fluid Solver on Adaptive Grids
This presentation will describe the numerical techniques, programming paradigms, verification, and performance of a non-equilibrium plasma fluid solver that can effectively utilize current and upcoming central processing and graphics processing unit (CPU+GPU) architectures. Our plasma fluid model solves the conservation equations for self-consistent electrostatic Poisson, electron and heavy species transport, and electron temperature on adaptive Cartesian grids. Our solver is written using performance portable adaptive mesh management library, AMReX (Zhang et al., JOSS, 4 (37) 1370, 2019), and can be built and run on widely available vendor specific GPU architectures (NVIDIA/AMD/Intel). We utilize a non-subcycled second order semi-implicit time-stepping method where all adaptive mesh refinement (AMR) levels are advanced with the same time step. The composite multi-level multigrid solver from within AMReX is used for each of the governing equations that are cast into a Helmholtz equation form. We have also developed a python based chemical mechanism parser framework that uses a similar format as CANTERA (Goodwin et al., Zenodo, 2018) yaml files as input. Our custom parser reads the yaml file and provides C++ files with transport and production rate functions that can be executed on both host (CPU) and device (GPU). We present verification of our solver using method of manufactured solutions that indicate formal second order accuracy with central diffusion and fifth order weighted-essentially-non-oscillatory (WENO) advection scheme. We also verify our solver with published literature on low-pressure capacitive and high-pressure streamer discharges. Our initial performance studies indicate 10X speed-up using 20 NVIDIA GPUs versus 200 CPUs for an atmospheric streamer discharge problem solved on a 512 x 1024 x 512 grid.
graphics processing units↗
Development of an Efficient CFD Model for Nuclear Thermal Thrust Chamber Assembly Design
The objective of this effort is to develop an efficient and accurate computational methodology to predict both detailed thermo-fluid environments and global characteristics of the internal ballistics for a hypothetical solid-core nuclear thermal thrust chamber assembly (NTTCA). Several numerical and multi-physics thermo-fluid models, such as real fluid, chemically reacting, turbulence, conjugate heat transfer, porosity, and power generation, were incorporated into an unstructured-grid, pressure-based computational fluid dynamics solver as the underlying computational methodology. The numerical simulations of detailed thermo-fluid environment of a single flow element provide a mechanism to estimate the thermal stress and possible occurrence of the mid-section corrosion of the solid core. In addition, the numerical results of the detailed simulation were employed to fine tune the porosity model mimic the pressure drop and thermal load of the coolant flow through a single flow element. The use of the tuned porosity model enables an efficient simulation of the entire NTTCA system, and evaluating its performance during the design cycle.
Cheng, Gary↗
Helium Pressurization Modeling of the EDU Submerged Diffuser Using Thermal Desktop
NASA is currently working on developing its predictive modeling capabilities to support future missions that require long term storage of cryogenic propellants. Through the Evolvable Cryogenics project, Development and Validation of Analysis Tools (DVAT) task, both CFD and multi-node fluid modeling tools are being validated against to increase the following cryogenic fluid management technologies: Self-Pressurization, Pressure Control, Helium Pressurization, de-stratification, jet mixing, Tank and transfer line chill-down, tank filling and transfer. NASA is using previous and future micro-gravity experiments to validate both types of models. Recent activities include: Validation of the Tank Pressure Control Experiment (TPCE), Hydrogen and Helium pressurization, and Completion of the Zero Boil Off Tank (ZBOT) experiment on ISS. Future activities include validation of the ZBOT experiment with both CFD and multi-node analysis tools, validation of the Cryogenic Demonstration System of the Robotic Refueling Mission (RRM3).
low-gravity fluid modeling↗
Fluids and Combustion Facility: Fluids Integrated Rack Modal Model Correlation
The Fluids Integrated Rack (FIR) is one of two racks in the Fluids and Combustion Facility on the International Space Station. The FIR is dedicated to the scientific investigation of space system fluids management supporting NASA s Exploration of Space Initiative. The FIR hardware was modal tested and FIR finite element model updated to satisfy the International Space Station model correlation criteria. The final cross-orthogonality results between the correlated model and test mode shapes was greater than 90 percent for all primary target modes.
McNelis, Mark E.↗
A parameter study of the two-fluid solar wind
A two-fluid model of the solar wind was introduced by Sturrock and Hartle (1966) and Hartle and Sturrock (1968). In these studies the proton energy equation was integrated neglecting the heat conductive term. Later several authors solved the equations for the two-fluid solar wind model keeping the proton heat conductive term. Methods where the equations are integrated simultaneously outward and inward from the critical point were used. The equations were also integrated inward from a large heliocentric distance. These methods have been applied to cases with low coronal base electron densities and high base temperatures. In this paper we present a method of integrating the two-fluid solar wind equations using an iteration procedure where the equations are integrated separately and the proton flux is kept constant during the integrations. The technique is applicable for a wide range of coronal base densities and temperatures. The method is used to carry out a parameter study of the two-fluid solar wind.
Sandbaek, Ornulf↗
Drift-kinetic effects of tungsten on plasma response to RMP in ITER
Here, effects of high- Z ( Z is the particle charge number) tungsten impurity ions on the plasma response to the resonant magnetic perturbation (RMP) field are numerically investigated for the ITER 15 MA baseline scenario, where the tungsten contribution to the plasma response is computed with a drift-kinetic model while the bulk thermal particle contributions follow the fluid approximation. The study yields three highlights: (i) the drift-kinetic contribution of the tungsten impurity exerts minor influence on the plasma response compared to that computed by the pure fluid model without tungsten; (ii) a new figure of merit, based on the resonant spectrum perturbation at the plasma boundary surface, results in different optimal coil phasing compared to that previously obtained by maximizing the edge-peeling plasma response; (iii) the optimal $n = 3$ RMP (for edge localized mode (ELM) control, $n$ is the toroidal mode number) is found to induce a large tungsten particle influx near the plasma edge associated with the neoclassical toroidal viscosity. The study thus provides useful data on the compatibility of the full tungsten wall with RMP ELM control in ITER.
ITER↗
A statistical formulation of one-dimensional electron fluid turbulence
A one-dimensional electron fluid model is investigated using the mathematical methods of modern fluid turbulence theory. Non-dissipative equilibrium canonical distributions are determined in a phase space whose co-ordinates are the real and imaginary parts of the Fourier coefficients for the field variables. Spectral densities are calculated, yielding a wavenumber electric field energy spectrum proportional to k to the negative second power for large wavenumbers. The equations of motion are numerically integrated and the resulting spectra are found to compare well with the theoretical predictions.
Fyfe, D.↗
Statistical formulation of one-dimensional electron fluid turbulence
A one-dimensional electron fluid model is investigated using the mathematical methods of modern fluid turbulence theory. Nondissipative equilibrium canonical distributions are determined in a phase space whose coordinates are the real and imaginary parts of the Fourier coefficients for the field variables. Spectral densities are calculated, yielding a wavenumber electric field energy spectrum proportional to 1/k squared for large wavenumbers. The equations of motion are numerically integrated and the resulting spectra are found to compare well with the theoretical predictions.
Fyfe, D.↗
The structure of mass-loading shocks. 1: Comets
A new multifluid model to describe the solar wind-cometary ion plasma in the outer cometary coma is derived. This model is distinguished from previous multifluid models in that wave-particle effects are included explicitly. By considering hydrodynamic timescales, self-consistent dissipative terms (which correspond to the spatial diffusion of cometary ions) are derived, thereby allowing the structure of a cometary bow shock to be resolved. In the hypersonic limit, our diffusive multifluid model reveals an attractive connection to the one-fluid model formulated by Biermann et al. (1967), a connection which serves as an important guide when investigating our more complicated nonhypersonic diffusive model. The detailed model consists of cometary ions diffusing in a warm solar wind which is heated by the dissipation of Alfvenic turbulence (in situ and generated by the initial ring beam distribution). The quasi-parallel cometary shock is found to be at least an order of magnitude thicker than a quasi-perpendicular cometary shock. It is furthermore found that for undisturbed solar wind (sonic) Mach numbers in excess of approximately 6, the cometary bow shock (for both perpendicular and parallel cases) is smoothed completely by the cometary ions and no 'proton' or thermal subshock is necessary. Finally, a detailed comparison of the model with the plasma observations made at Halley by Giotto is undertaken. It is found that the observed shock structure, thickness, location, and plasma parameters compare well with the predictions of the theoretical model, particularly in the case of the quasi-parallel shock.
Zank, G. P.↗
LTF/Fluid and Cryogenic Systems Branch
Functions of the Branch: Cross-Cutting Fluids Engineering; Research and Technology Development of Cryogenic Fluid Management (CFM) Systems; Component and System level; Design, Development Test and Evaluation (DDT&E) for Fluid, CFM, Propellant and Gaseous Systems; Fluid Component Specification and Selection; Code Compliance of Pressurized Systems; Fluid/Cryogenic System Performance Analysis and Trade Studies; Hazards Analysis and Risk Assessment of Fluidic Systems; Fluids Modeling and Simulation Capabilities
Fluid↗
Nodal Modeling of Liquid Propellant Feed and Pressurization System
Nodal or network flow modeling plays an important role in the design and operation of the feed and pressurization system of a liquid rocket engine. Model development and execution time is relatively short for nodal codes in comparison to the Navier-Stokes based CFD codes. Nodal models also allow the inclusion of several components into one model to predict the behavior of a larger system. Unlike CFD models, the nodal models do not need very fine discretization of the flow field because they use empirical correlations to model fluid friction and heat transfer. This paper presents several applications of nodal modeling of liquid propellant feed and pressurization systems using the Generalized Fluid System Simulation Program (GFSSP), a nodal code developed at NASA/ Marshall Space Flight Center. GFSSP discretizes the flow field into nodes which are connected by branches. The mass and energy conservation equations and the equation of state are solved to calculate pressure, temperature, and resident mass at the nodes. The momentum equations are solved at the branches to calculate flow rates. Applications include a) tank pressurization by inert gas as well as autogenous pressurization by gaseous propellant, b) submerged pressurization by helium, c) self-pressurization due to boil-off of cryogenic propellant, d) chilldown of a transfer line of a cryogenic propellant feed system, and e) chilldown and filling of a cryogenic tank. Each of the above-mentioned models was verified and validated by comparing with test data.
Nodal Model↗
Nodal Modeling of Liquid Propellant Feed and Pressurization System
Nodal or network flow modeling plays an important role in the design and operation of the feed and pressurization system of a liquid rocket engine. Model development and execution time is relatively short for nodal codes in comparison to the Navier-Stokes based CFD codes. Nodal models also allow the inclusion of several components into one model to predict the behavior of a larger system. Unlike CFD models, the nodal models do not need very fine discretization of the flow field because they use empirical correlations to model fluid friction and heat transfer. This paper presents several applications of nodal modeling of liquid propellant feed and pressurization systems using the Generalized Fluid System Simulation Program (GFSSP), a nodal code developed at NASA/ Marshall Space Flight Center. GFSSP discretizes the flow field into nodes which are connected by branches. The mass and energy conservation equations and the equation of state are solved to calculate pressure, temperature, and resident mass at the nodes. The momentum equations are solved at the branches to calculate flow rates. Applications include a) tank pressurization by inert gas as well as autogenous pressurization by gaseous propellant, b) submerged pressurization by helium, c) self-pressurization due to boil-off of cryogenic propellant, d) chilldown of a transfer line of a cryogenic propellant feed system, and e) chilldown and filling of a cryogenic tank. Each of the above-mentioned models was verified and validated by comparing with test data.
Nodal Model↗
Fluid dynamic and thermodynamic analysis of a model pertaining to cryogenic fluid management in low gravity environments for a system with dynamically induced settling
The settling behavior of the liquid and gaseous phases of a fluid in a propellant and in a zero-g environment, when such settling is induced through the use of a dynamic device, in this particular case, a helical screw was studied. Particular emphasis was given to: (1) the description of a fluid mechanics model which seems applicable to the system under consideration, (2) a First Law of Thermodynamics analysis of the system, and (3) a discussion of applicable scaling rules.
Rios, J.↗