Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Conservation 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.

At least 487 records · Page 27
All sourcesNASA NTRSDOE OSTI
NASA NTRS2018

One-Dimensional Modeling Methodology for Shock Tubes: Application to the EAST Facility

In this work, a one-dimensional methodology for simulating shock tubes is developed. The model accounts for the viscous interactions of the shock with the shock tube wall by adding an area change source term in the 1-D conservation equations corresponding to the boundary layer growth. This source term corresponds to the mass and energy going into the boundary layer. The boundary layer growth is computed using a simple model with a scaling factor. This scale factor is used to tailor a solution to match the deceleration profile of a shock tube test. In doing so, not only will the source term take into account boundary layer losses, it will also cover any effect due to radiative cooling loses from the gas. For this study, the Electric Arc Shock Tube(EAST) facility at NASA Ames Research Center is modeled for Earth reentry conditions. The purpose of this paper is to investigate if anomalies identified for certain conditions in the EAST data are due to shock deceleration. These anomalies include measuring electron number density above equilibrium predictions and observing that radiance profiles can continually increase behind the shock, never reach steady state, for certain shots (typically those less than 10 km/s). An eleven species air mixture is chosen to study the chemistry of the flow. Comparisons of the simulations to the experimental results are presented. Good agreement with the shock deceleration profiles was achieved by tuning in the boundary layer scale factor. The temperature as well as electron number density increases behind the shock, as has also been observed in the experiments. Finally, radiance comparisons between results from NEQAIR and experiments also show good agreement for some shots, but significant discrepancies are still observed for others.

Sharma, Maitreyee↗
NASA NTRSDate not supplied

Field Reconstruction from PIV Measurements Employing Bernstein Polynomial Derived Operators

A fluid-dynamic reconstruction algorithm is presented that generates a least-squares best-fit, two-dimensional density field from a prespecified two-dimensional velocity field. This method recasts the mass-conservation equation as a modified Sylvester equation employing high-order operators derived from modified Bernstein polynomial expansions. To demonstrate its practical utility, this analytic methodology is applied to two canonical cases and a Particle Image Velocimetry dataset obtained from a Mach-2, mechanically back-pressured, isolator experiment. This methodology is envisioned to be used in conjunction with hypersonic-diagnostic techniques to aid in the quantification of isolator flow fields. However, also note that this reconstruction technique is well suited to other applications relevant to fluid dynamics, such as obtaining three-dimensional flow field reconstructions.

Bernstein Polynomials↗
NASA NTRSDate not supplied

Constrained Field Construction Using Bernstein Polynomial Derived Operators

This method recasts the mass-conservation equation as a Sylvester equation employing high-order derivative operators derived from modified Bernstein polynomial expansions. Given prescribed velocity fields, the algorithm yields a constrained, least-squared solution for the associated density field. To demonstrate the practical utility of this methodology, it is applied to a computationally-derived two-dimensional isolator dataset. This reconstruction method is envisioned to be used in conjunction with diagnostic techniques to aid in the quantification of isolator flow fields, since obtaining highly characterized datasets within this engine component is exceedingly difficult.

Isolator↗
NASA NTRS2024

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↗
NASA NTRS2024

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↗
DOE OSTI2020

Tracking Radioactive Isotopes in HVAC and Application for Hot Cell Analyses - 20155

In recent years, the surge in the number of isotope production facilities under design has increased the need to analyze isotope migration in facility Heating, Ventilation, and Air Conditioning (HVAC) during normal operations and accident scenarios. These facilities are used to produce isotopes for medical, security, and industrial applications and as such are subject to license and regulatory requirements 10CFR20, 10CFR30, 10CFR50, and 10CFR70. GOTHIC, a general-purpose thermal-hydraulics software package, includes the ability to model isotopic tracers, radioactive decay and isotope migration as well as HEPA and charcoal filters for isotope retention. A GOTHIC model was developed by Zachry Nuclear Engineering (ZNE) to examine the effect of negative room pressure, HEPA filtration, and HVAC fluctuations on radiation areas, hot cells and gloveboxes. Radioactive tracers were used to simulate the concentration of spills within contaminated areas, track the migration of isotopes of interest, and determine the isotopic retention and buildup on facility HEPA filters. A variety of isotopes with concern to dose (e.g., Kr-85, Sr-90, I-131, etc.), including their decay and progeny, are included in the analysis. Negative pressures are maintained in the regions of interest by a representative central HVAC system equipped with a volumetric fan that exhausts to the environment after a series of isolation valves and HEPA filters. GOTHIC is an industry trusted tool for providing engineering solutions for a variety of applications, including fission product tracking, aerosol and particulate transport and ventilation assessments. The software provides an integrated analysis environment that includes a graphical user interface (GUI) for constructing analysis models, a numerical solver that includes parallel processing capabilities and a post-processor for evaluating simulation results. It solves the conservation equations for mass, momentum and energy for multicomponent, multi-phase flow in lumped parameter and multi-dimensional geometries (1, 2, or full 3D), including the effects of turbulence, diffusion and buoyancy. It has been developed and maintained under a Quality Assurance program in compliance with the requirements of 10CFR50 Appendix B and applicable portions of ASME NQA-1 since 1995. GOTHIC has been used for assessing both forced and natural convection conditions for a wide range of applications, including: - Tracking concentration of hazardous gases and chemicals for habitability and safety assessments - Determining ventilation and filtration requirements and optimizing location and arrangement of these systems - Room heat-up, including diverse and Flexible coping strategies for Extended Loss of AC Power (FLEX/ELAP) - Equipment Qualification (EQ) A distinctive feature of GOTHIC is the ability to track many different fields/substances in a simulation, including user defined tracer elements, in the liquid, vapor and droplet fields as well as surfaces and filters. This capability allows GOTHIC to model fission product transport and release or the removal of particulates or harmful toxins from exhaust gases using a spray scrubber or other types of filtration systems. GOTHIC also includes models for engineered equipment, such as fans, filters, charcoal filters, dryers/demisters, dampers, etc. The aerosols and other filtered material are removed or accumulated in these components. The range of aerosol and radiological applications that GOTHIC has been used for includes: - Source Term: Primary Coolant (Equilibrium) Activity; Non-Water Coolant Source/Leakage. - Conditions for Iodine Re-evolution: Sump/Suppression Pool Conditions and pH; - RWST Conditions and pH. - Isotope Removal Mechanisms: Containment Sprayed and Unsprayed Region Mixing; Charcoal Filter Heating due to Iodine decay. - Radionuclide Transport and Decay: Post-LOCA Release in containment; Transport between connected Compartments and vent systems; Groundwater transport of radionuclides. - Non-Newtonian Fluid modeling for sludge, waste tanks, etc. ADAMS ML071581053 (titled 'Best Practice Guidelines for the use of CFD in Nuclear Reactor Safety Applications') poses guidelines for applying single phase CFD codes in nuclear reactor safety problems and GOTHIC is listed as a 'tool for 3D flows' and 'dispersal and deposition of radionuclides.' The Nuclear Quality Assurance (NQA) pedigree of GOTHIC is an important aspect for applications in the nuclear industry. The fundamental tracer models (convective transport, molecular and turbulent diffusion, removal mechanisms, etc.) have been verified using analytical solutions and validated against applicable separate effects tests. Also, GOTHIC has been benchmarked to many integrated effects tests, including Phebus FP (Fission Product). GOTHIC gives good agreement for the buildup and decay of fission products in Phebus Test 3. The model developed by ZNE demonstrates GOTHIC's applicability and acceptability for use in analyzing the migration and retention of radioactive isotopes and their progeny in normal operation and accident scenario analyses for isotope production facilities, hot cells, and gloveboxes. The tracer activities calculated by GOTHIC can then be used in downstream radiation transport and shielding codes like RADTRAD-NAI{sup C}, MCNP{sup R}, and MicroShield{sup R} to determine on-site and of-site doses. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
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 NTRS1988

Universal limiter for transient interpolation modeling of the advective transport equations: The ULTIMATE conservative difference scheme

A fresh approach is taken to the embarrassingly difficult problem of adequately modeling simple pure advection. An explicit conservative control-volume formation makes use of a universal limiter for transient interpolation modeling of the advective transport equations. This ULTIMATE conservative difference scheme is applied to unsteady, one-dimensional scalar pure advection at constant velocity, using three critical test profiles: an isolated sine-squared wave, a discontinuous step, and a semi-ellipse. The goal, of course, is to devise a single robust scheme which achieves sharp monotonic resolution of the step without corrupting the other profiles. The semi-ellipse is particularly challenging because of its combination of sudden and gradual changes in gradient. The ULTIMATE strategy can be applied to explicit conservation schemes of any order of accuracy. Second-order schemes are unsatisfactory, showing steepening and clipping typical of currently popular so-called high resolution shock-capturing of TVD schemes. The ULTIMATE third-order upwind scheme is highly satisfactory for most flows of practical importance. Higher order methods give predictably better step resolution, although even-order schemes generate a (monotonic) waviness in the difficult semi-ellipse simulation. Little is to be gained above ULTIMATE fifth-order upwinding which gives results close to the ultimate for which one might hope.

Leonard, B. P.↗
NASA NTRS1991

Conservative multizonal interface algorithm for the 3-D Navier-Stokes equations

A conservative zonal interface algorithm using features of both structured and unstructured mesh CFD technology is presented. The flow solver within each of the zones is based on structured mesh CFD technology. The interface algorithm was implemented into two three-dimensional Navier-Stokes finite volume codes and was found to yield good results.

Klopfer, G. H.↗
NASA NTRS1979

Fast, conservative schemes for the full potential equation applied to transonic flows

Implicit approximate factorization techniques (AF) are investigated for the solution of matrix equations resulting from finite-difference approximations to the full potential equation in conservation form. For transonic flows, an artificial viscosity, required to maintain stability in supersonic regions, is introduced by an upwind bias of the density. Two implicit AF procedures are presented, and their convergence performance is compared with that of the standard transonic solution procedure: successive line overrelaxation (SLOR). Subcritical and supercritical test cases are considered. Results indicate a substantial improvement in convergence rate for AF schemes relative to SLOR.

Holst, T. L.↗
NASA NTRS1978

Conservative implicit schemes for the full potential equation applied to transonic flows

Implicit approximate factorization techniques (AF) were investigated for the solution of matrix equations resulting from finite difference approximations to the full potential equation in conservation form. For transonic flows, an artificial viscosity, required to maintain stability in supersonic regions, was introduced by an upwind bias of the density. Two implicit AF procedures are presented and their convergence performance is compared with that of the standard transonic solution procedure, successive line overrelaxation (SLOR). Subcritical and supercritical test cases are considered. The results indicate that the AF schemes are substantially faster than SLOR.

Holst, T. L.↗
NASA NTRS1985

Analysis of two-dimensional incompressible flow past airfoils using unsteady Navier-Stokes equations

The conservative form of the unsteady Navier-Stokes equations in terms of vorticity and stream function in generalized curvilinear coordinates are used to analyze the flow structure of steady separation and unsteady flow with massive separation. The numerical method solves the discretized equations using an ADI-BGE method. The method is applied to a symmetric 12 percent thick Joukowski airfoil. A conformal clustered grid is generated; several 1-D stretching transformations are used to obtain a grid that attempts to resolve many of the multiple scales of the unsteady flow with massive separation, while maintaining the transformation metrics to be smooth and continuous in the entire flow field. Detailed numerical results are obtained for three flow configurations (1) Re = 1000, alpha = 5 deg, (2) Re = 1000, alpha = 15 deg, (3) Re = 10,000, alpha = 5 deg. No artificial dissipation was added; however, lack of a fine grid in the normal direction has presently led to results which are considered qualitative, especially for case (3).

Ghia, K. N.↗
NASA NTRS1986

Analysis of two-dimensional incompressible flow past airfoils using unsteady Navier-Stokes equations

The conservative form of the unsteady Navier-Stokes equations in terms of vorticity and stream function in generalized curvilinear coordinates are used to analyze the flow structure of steady separation and unsteady flow with massive separation. The numerical method solves the discretized equations using an ADI-BGE method. The method is applied to a symmetric 12 percent thick Joukowski airfoil. A conformal clustered grid is generated; several 1-D stretching transformations are used to obtain a grid that attempts to resolve many of the multiple scales of the unsteady flow with massive separation, while maintaining the transformation metrics to be smooth and continuous in the entire flow field. Detailed numerical results are obtained for three flow configurations (1) Re = 1000, alpha = 5 deg., (2) Re =1000, alpha = 15 deg., (3) Re = 10,000, alpha = 5 deg. No artificial dissipation was added; however, lack of a fine grid in the normal direction has presently led to results which are considered qualitative, especially for case (3).

Ghia, K. N.↗
NASA NTRS1979

Computations of unsteady transonic flow governed by the conservative full potential equation using an alternating direction implicit algorithm

A development was the time linearization of the density function. This linearization reduces the solution process from solving just a single equation. Two sample cases were computed. First, a one dimensional traveling shock wave was computed and compared with the analytic solution. Second, a two dimensional case was calculated for a flow field which resulted from a thickening and subsequently, thinning airfoil. The resulting flow field, which included a traveling shock wave, was compared to the flow field obtained from the low frequency, small disturbance, transonic equation.

Goorjian, P. M.↗
NASA NTRS1988

Unsteady vortex-dominated flows around maneuvering wings over a wide range of Mach numbers

The problem of unsteady flow around maneuvering wings is solved using the unsteady Euler equations. The unsteady conservative Euler equations are derived for the flow relative motion with respect to a moving (translating and rotating) frame of reference. The resulting equations can handle the most general case for unsteady three-dimensional flow around maneuvering wings or wing-body configurations undergoing six degrees of freedom motion; three translations and three rotations. The equations are solved using two computational schemes; an explicit multistage finite-volume scheme and an implicit approximately factored finite-volume scheme. The computational applications cover two cases. The first case is for a locally conical supersonic flow of rolling oscillation of a sharp-edged delta wing at zero angle of attack. The second case is for a pitching oscillation around a mean angle of attack of a NACA 0012 airfoil in transonic flow.

Kandil, Osama A.↗
NASA NTRS1988

Generalized energy and potential enstrophy conserving finite difference schemes for the shallow water equations

The conditions under which finite difference schemes for the shallow water equations can conserve both total energy and potential enstrophy are considered. A method of deriving such schemes using operator formalism is developed. Several such schemes are derived for the A-, B- and C-grids. The derived schemes include second-order schemes and pseudo-fourth-order schemes. The simplest B-grid pseudo-fourth-order schemes are presented.

Abramopoulos, Frank↗