SEARCH · Engineering Papers

Results for “Solving Fluid”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 37 records · Page 2

All sources NASA NTRS DOE OSTI

NASA NTRS2015

Results of Microgravity Fluid Dynamics Captured with the Spheres-Slosh Experiment

This paper provides an overview of the SPHERES-Slosh Experiment (SSE) aboard the International Space Station (ISS) and presents on-orbit results with data analysis. In order to predict the location of the liquid propellant during all times of a spacecraft mission, engineers and mission analysts utilize Computational Fluid Dynamics (CFD). These state-of-the-art computer programs numerically solve the fluid flow equations to predict the location of the fluid at any point in time during different spacecraft maneuvers. The models and equations used by these programs have been extensively validated on the ground, but long duration data has never been acquired in a microgravity environment. The SSE aboard the ISS is designed to acquire this type of data, used by engineers on earth to validate and improve the CFD prediction models, improving the design of the next generation of space vehicles as well as the safety of current missions. The experiment makes use of two Synchronized Position Hold, Engage, Reorient Experimental Satellites (SPHERES) connected by a frame. In the center of the frame there is a plastic, pill shaped tank that is partially filled with green-colored water. A pair of high resolution cameras records the movement of the liquid inside the tank as the experiment maneuvers within the Japanese Experimental Module test volume. Inertial measurement units record the accelerations and rotations of the tank, making the combination of stereo imaging and inertial data the inputs for CFD model validation.

slosh↗

NASA NTRS2015

Result of Microgravity Fluid Dynamics Captured with the SPHERES-Slosh Experiment

slosh↗

NASA NTRS1993

Towards a generalized computational fluid dynamics technique for all Mach numbers

Currently there exists no single unified approach for efficiently and accurately solving computational fluid dynamics (CFD) problems across the Mach number regime, from truly low speed incompressible flows to hypersonic speeds. There are several CFD codes that have evolved into sophisticated prediction tools with a wide variety of features including multiblock capabilities, generalized chemistry and thermodynamics models among other features. However, as these codes evolve, the demand placed on the end user also increases simply because of the myriad of features that are incorporated into these codes. In order for a user to be able to solve a wide range of problems, several codes may be needed requiring the user to be familiar with the intricacies of each code and their rather complicated input files. Moreover, the cost of training users and maintaining several codes becomes prohibitive. The objective of the current work is to extend the compressible, characteristic-based, thermochemical nonequilibrium Navier-Stokes code GASP to very low speed flows and simultaneously improve convergence at all speeds. Before this work began, the practical speed range of GASP was Mach numbers on the order of 0.1 and higher. In addition, a number of new techniques have been developed for more accurate physical and numerical modeling. The primary focus has been on the development of optimal preconditioning techniques for the Euler and the Navier-Stokes equations with general finite-rate chemistry models and both equilibrium and nonequilibrium thermodynamics models. We began with the work of Van Leer, Lee, and Roe for inviscid, one-dimensional perfect gases and extended their approach to include three-dimensional reacting flows. The basic steps required to accomplish this task were a transformation to stream-aligned coordinates, the formulation of the preconditioning matrix, incorporation into both explicit and implicit temporal integration schemes, and modification of the numerical flux formulae. In addition, we improved the convergence rate of the implicit time integration schemes in GASP through the use of inner iteration strategies and the use of the GMRES (General Minimized Resisual) which belongs to the class of algorithms referred to as Krylov subspace iteration. Finally, we significantly improved the practical utility of GASP through the addition of mesh sequencing, a technique in which computations begin on a coarse grid and get interpolated onto successively finer grids. The fluid dynamic problems of interest to the propulsion community involve complex flow physics spanning different velocity regimes and possibly involving chemical reactions. This class of problems results in widely disparate time scales causing numerical stiffness. Even in the absence of chemical reactions, eigenvalue stiffness manifests itself at transonic and very low speed flows which can be quantified by the large condition number of the system and evidenced by slow convergence rates. This results in the need for thorough numerical analysis and subsequent implementation of sophisticated numerical techniques for these difficult yet practical problems. As a result of this work, we have been able to extend the range of applicability of compressible codes to very low speed inviscid flows (M = .001) and reacting flows.

Walters, R. W.↗

NASA NTRS2019

Numerical Modeling of Boiling in a Heated Tube

This paper presents numerical models of boiling in a heated tube using the Generalized Fluid System Simulation Program (GFSSP), a finite-volume-based general-purpose flow network code developed at NASA/Marshall Space Flight Center. The heated tube is discretized into a one-dimensional array of nodes and branches to represent the flow of liquid and vapor in a tube with a prescribed pressure differential. The solid wall is also discretized into solid nodes and conductors to allow for heat transfer between the wall and the fluid. The conservation equations of mass, momentum, and energy of the fluid are solved simultaneously with the energy conservation equation for the solid wall. Two experimental configurations of fluid flowing in a vertical tube have been simulated, one with water and the other with liquid hydrogen. This paper compares experimental data with numerical predictions based on four different published correlations for boiling heat transfer coefficients. Three of these correlations are applicable to the saturated vertical flow conditions of the experiments. One of them is applicable to film boiling and has been used for the liquid hydrogen experiment, which was in film boiling regime. For the case of boiling water, the predictions of wall temperatures using the boiling heat transfer correlations agreed well with the experimental results. However, in the case of boiling hydrogen larger discrepancies were observed between the experimental data and numerical predictions.

Majumdar, Alok↗

NASA NTRS1983

Acoustic responses of coupled fluid-structure system by acoustic-structural analogy

The use of an analogy between structural mechanics and acoustics makes it possible to solve fluid-structural interaction (FSI) problems using an existing structural analysis computer program. This method was implemented in MSC/NASTRAN program and the FSI analysis was performed using two dimensional coupled fluid beam model to assess and evaluate the adequacy of this approach. The coupled modal analysis of 3-D model is also briefly discussed. The normal mode, modal frequency response and transient response analysis of 2-D coupled fluid beam system is presented. The significant reduction of the acoustic pressure response at the fluid structure interface is observed as a result of fluid structure interaction.

Shin, Y. S.↗

NASA NTRS2003

Finite Element Method for Capturing Ultra-relativistic Shocks

While finite element methods are used extensively by researchers solving computational fluid dynamics in fields other than astrophysics, their use in astrophysical fluid simulations has been predominantly overlooked. Current simulations using other methods such as finite difference and finite volume (based on finite difference) have shown remarkable results, but these methods are limited by their fundamental properties in aspects that are important for simulations with complex geometries and widely varying spatial and temporal scale differences. We have explored the use of finite element methods for astrophysical fluids in order to establish the validity of using such methods in astrophysical environments. We present our numerical technique applied to solving ultra-relativistic (Lorentz Factor Gamma >> 1) shocks which are prevalent in astrophysical studies including relativistic jets and gamma-ray burst studies. We show our finite element formulation applied to simulations where the Lorentz factor ranges up to 2236 and demonstrate its stability in solving ultra-relativistic flows. Our numerical method is based on the Flowfield Dependent Variation (FDV) Method, unique in that numerical diffusion is derived from physical parameters rather than traditional artificial viscosity methods. Numerical instabilities account for most of the difficulties when capturing shocks in this regime. Our method results in stable solutions and accurate results as compared with other methods.

Richardson, G. A.↗

NASA NTRS2011

Numerical Modeling of the Chilldown of Cryogenic Transfer Lines Using a Sinda/GFSSP Integrated Solver

An important first step in cryogenic propellant loading is the chilldown of transfer lines. During the chilldown of the transfer line, the flow is two-phase and unsteady, with solid to fluid heat transfer and therefore a coupled thermo-fluid analysis is necessary to model the system. This paper describes a numerical model of pipe chilldown that utilizes the Sinda/GFSSP Conjugate Integrator (SGCI). SGCI is a new analysis tool developed at NASA's Marshall Space Flight Center (MSFC). SGCI facilitates the solution of thermofluid problems in interconnected solid-fluid systems. The solid component of the system is modeled in MSC Patran and translated into an MSC Sinda thermal network model. The fluid component is modeled in GFSSP, the Generalized Fluid System Simulation Program. GFSSP is a general network flow solver developed at NASA/MSFC. GFSSP uses a finite-volume approach to model fluid systems that can include phase change, multiple species, fluid transients, and heat transfer to simple solid networks. SGCI combines the GFSSP Fortran code with the Sinda input file and compiles the integrated model. Sinda solves for the temperatures of the solid network, while GFSSP simultaneously solves the fluid network for pressure, temperature, and flow rate. The two networks are coupled by convection heat transfer from the solid wall to the cryogenic fluid. The model presented here is based on a series of experiments conducted in 1966 by the National Bureau of Standards (NBS). A vacuum-jacketed, 200 ft copper transfer line was chilled by liquid nitrogen and liquid hydrogen. The predictions of transient temperature profiles and chilldown time of the integrated Sinda/GFSSP model will be compared to the experimental measurements.

LeClair, Andre↗

DOE OSTI2024

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↗

NASA NTRS1981

Grid evolution in time asymptotic problems

The selection of the proper coordinate system in solving any fluid flow or heat transfer problem is a very important consideration. A new technique of moving mesh points in physical space is introduced so as to reduce the error in a computed asymptotic solution relative to that obtained using a fixed mesh. The technique has been used to solve the simple viscous Burgers' equation in one and two dimensions. Substantial error reductions are demonstrated. The treatment of boundary points and the effect of using different error criteria in generating grids are discussed.

Rai, M. M.↗

NASA NTRS1993

The NASA Computational Fluid Dynamics (CFD) program - Building technology to solve future challenges

This paper presents the NASA Computational Fluid Dynamics program in terms of a strategic vision and goals as well as NASA's financial commitment and personnel levels. The paper also identifies the CFD program customers and the support to those customers. In addition, the paper discusses technical emphasis and direction of the program and some recent achievements. NASA's Ames, Langley, and Lewis Research Centers are the research hubs of the CFD program while the NASA Headquarters Office of Aeronautics represents and advocates the program.

Richardson, Pamela F.↗

NASA NTRS2001

Modeling of Chill Down in Cryogenic Transfer Lines

A numerical model to predict chill down in cryogenic transfer lines has been developed. Three chill down cases using hydrogen as the working fluid are solved: 1) a simplified model amenable to analytical solution, 2) a realistic model of superheated vapor flow, and 3) a realistic model of initially subcooled liquid flow. The first case compares a numerical model with an analytical solution with very good agreement between the two. Additionally, the analytical solution provides a convenient way to look at parametric effects on the chill down. The second and third cases are numerical models which provide temperature histories of the fluid and solid tube wall during chill down as well as several other quantities of interest such as pressure and mass flow rate. Of great interest is the ability to predict accurate values of chill down time (the time required to achieve steady-state cryogenic flow). The models predict that a 26 in. long, 3/16 in. ID aluminum tube has a shorter chill down time (approx. equal to 100 sec) and uses less hydrogen with superheated vapor flow than with initially subcooled liquid flow (greater than 200 sec for chill down).

Cross, Matthew F.↗

NASA NTRS2004

Rapid Structured Volume Grid Smoothing and Adaption Technique

A rapid, structured volume grid smoothing and adaption technique, based on signal processing methods, was developed and applied to the Shuttle Orbiter at hypervelocity flight conditions in support of the Columbia Accident Investigation. Because of the fast pace of the investigation, computational aerothermodynamicists, applying hypersonic viscous flow solving computational fluid dynamic (CFD) codes, refined and enhanced a grid for an undamaged baseline vehicle to assess a variety of damage scenarios. Of the many methods available to modify a structured grid, most are time-consuming and require significant user interaction. By casting the grid data into different coordinate systems, specifically two computational coordinates with arclength as the third coordinate, signal processing methods are used for filtering the data [Taubin, CG v/29 1995]. Using a reverse transformation, the processed data are used to smooth the Cartesian coordinates of the structured grids. By coupling the signal processing method with existing grid operations within the Volume Grid Manipulator tool, problems related to grid smoothing are solved efficiently and with minimal user interaction. Examples of these smoothing operations are illustrated for reduction in grid stretching and volume grid adaptation. In each of these examples, other techniques existed at the time of the Columbia accident, but the incorporation of signal processing techniques reduced the time to perform the corrections by nearly 60%. This reduction in time to perform the corrections therefore enabled the assessment of approximately twice the number of damage scenarios than previously possible during the allocated investigation time.

Alter, Stephen J.↗

NASA NTRS2006

Rapid Structured Volume Grid Smoothing and Adaption Technique

A rapid, structured volume grid smoothing and adaption technique, based on signal processing methods, was developed and applied to the Shuttle Orbiter at hypervelocity flight conditions in support of the Columbia Accident Investigation. Because of the fast pace of the investigation, computational aerothermodynamicists, applying hypersonic viscous flow solving computational fluid dynamic (CFD) codes, refined and enhanced a grid for an undamaged baseline vehicle to assess a variety of damage scenarios. Of the many methods available to modify a structured grid, most are time-consuming and require significant user interaction. By casting the grid data into different coordinate systems, specifically two computational coordinates with arclength as the third coordinate, signal processing methods are used for filtering the data [Taubin, CG v/29 1995]. Using a reverse transformation, the processed data are used to smooth the Cartesian coordinates of the structured grids. By coupling the signal processing method with existing grid operations within the Volume Grid Manipulator tool, problems related to grid smoothing are solved efficiently and with minimal user interaction. Examples of these smoothing operations are illustrated for reductions in grid stretching and volume grid adaptation. In each of these examples, other techniques existed at the time of the Columbia accident, but the incorporation of signal processing techniques reduced the time to perform the corrections by nearly 60%. This reduction in time to perform the corrections therefore enabled the assessment of approximately twice the number of damage scenarios than previously possible during the allocated investigation time.

Alter, Stephen J.↗

DOE OSTI2026

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↗

DOE OSTI2025

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 twin electrode radio frequency reactor. Our performance studies on three different CPU+GPU architectures indicate approximately 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.

Sitaraman, Hariswaran↗

NASA NTRS1973

Instabilities in fissioning plasmas as applied to the gas-core nuclear rocket-engine

The compressional wave spectrum excited in a fissioning uranium plasma confined in a cavity such as a gas cored nuclear reactor, is studied. Computer results are presented that solve the fluid equations for this problem including the effects of spatial gradients, nonlinearities, and neutron density gradients in the reactor. Typically the asymptotic fluctuation level for the plasma pressure is of order 1 percent.

Source record↗

NASA NTRS1977

Numerical methods for inviscid/viscous flows

A study of numerical schemes for solving viscous fluid flow problems with sizable regions of predominantly inviscid flow is presented. Difficulties associated with the familiar central difference approach for such problems were analyzed and alternative finite difference approaches employing windward concepts are presented. In addition, difference relations based on exponential operators were developed. All such schemes were demonstrated and evaluated through application to the case of Falkner Skan flow with blowing - a problem in which a sizable region of predominantly inviscid flow develops near the injection surface that traditionally causes numerical difficulty.

El-Mistikawy, T. M.↗

NASA NTRS1979

A two-dimensional unsteady Euler-equation solver for flow regions with arbitrary boundaries

A new technique is described for solving supersonic fluid dynamic problems containing multiple regions of continuous flow, each bounded by a permeable or impermeable surface. Region boundaries are, in general, arbitrarily shaped and time dependent. Discretization of such a region for solution by conventional finite difference procedures is accomplished using an elliptic solver which alleviates the dependence on a particular base coordinate system. Multiple regions are coupled together through the boundary conditions. The technique has been applied to a variety of problems including the shock diffraction and pointed wedge with detached bow shock.

Hindman, R. G.↗