Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “flow conservation”

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 631 records · Page 35

Performance of the Extravehicular Mobility Unit (EMU): Airlock Coolant Loop Recovery (A/L CLR) Hardware - Phase II

An EMU water processing kit (Airlock Coolant Loop Recovery A/L CLR) was developed as a corrective action to Extravehicular Mobility Unit (EMU) coolant flow disruptions experienced on the International Space Station (ISS) in May of 2004 and thereafter. Conservative schedules for A/L CLR use and component life were initially developed and implemented based on prior analysis results and analytical modeling. The examination of postflight samples and EMU hardware in November of 2006 indicated that the A/L CLR kits were functioning well and had excess capacity that would allow a relaxation of the initially conservative schedules of use and component life. A relaxed use schedule and list of component lives was implemented thereafter. Since the adoption of the relaxed A/L CLR schedules of use and component lives, several A/L CLR kit components, transport loop water samples and sensitive EMU transport loop components have been examined to gage the impact of the relaxed requirements. The intent of this paper is to summarize the findings of that evaluation, and to outline updated schedules for A/L CLR use and component life.

Steele, John↗

A High Order Discontinuous Galerkin Method for 2D Incompressible Flows

In this paper we introduce a high order discontinuous Galerkin method for two dimensional incompressible flow in vorticity streamfunction formulation. The momentum equation is treated explicitly, utilizing the efficiency of the discontinuous Galerkin method The streamfunction is obtained by a standard Poisson solver using continuous finite elements. There is a natural matching between these two finite element spaces, since the normal component of the velocity field is continuous across element boundaries. This allows for a correct upwinding gluing in the discontinuous Galerkin framework, while still maintaining total energy conservation with no numerical dissipation and total enstrophy stability The method is suitable for inviscid or high Reynolds number flows. Optimal error estimates are proven and verified by numerical experiments.

Liu, Jia-Guo↗

Two-dimensional hybrid continuum/particle approach for rarefied flows

A hybrid numerical technique previously developed for one-dimensional rarefied gas flows is generalized to two dimensions. The method is based on the fact that the flowfield that develops near a body in a rarefied gas typically contains local regions of continuum, transitional and free-molecular flow. By utilizing the solution technique most appropriate for each region and coupling the techniques in an interface region where both are applicable, more computationally efficient solutions can be obtained, or equivalently, more complex flowfields can be analyzed. The present method combines finite difference solution of the Navier-Stokes equations in the continuum regions, with Direct Simulation Monte Carlo in the more rarefied regions. The two schemes are coupled interactively via a general conservative flux boundary condition. The method is tested by application to the model problem of pressure-driven rarefied flow through a slit. Results show the hybrid scheme offers a speedup of a factor of nearly two for the nominal conditions considered, due to the decrease in the size of the Direct Simulation domain.

Wadsworth, Dean C.↗

A computational method for viscous incompressible flows

An implicit, finite-difference procedure for numerically solving viscous incompressible flows is presented. The pressure-field solution is based on the pseudocompressibility method in which a time-derivative pressure term is introduced into the mass-conservation equation to form a set of hyperbolic equations. The pressure-wave propagation and the spreading of the viscous effect is investigated using simple test problems. Computed results for external and internal flows are presented to verify the present method which has proved to be very robust in simulating incompressible flows.

Kwak, D.↗

Assessment of Pseudoshock Models Against Experiment in a Low-Aspect-Ratio Isolator

A highly confined shock train is investigated in a direct-connect isolator facility with a Mach 2 inflow and a constant-area low-aspect-ratio rectangular test section. High-speed schlieren imaging, wall static pressure measurements, surface oil-flow visualization, and particle image velocimetry from this isolator are synthesized into a three-dimensional schematic of the shock train structure. Against this, the prevailing pseudoshock models in the literature are assessed to evaluate the validity of their underlying assumptions. None of the prevailing pseudoshock models are found to simultaneously model the pressure and Mach number profiles, indicating a gap in the model formation and underlying assumptions when applied to the experimental isolator of interest. The presence of distortion in the isolator flowfield, such as a wall-bounded vortex, is found to skew the structure of the shock train, altering the strength and distribution of the compressive pressure gradient. It is further observed that the separated flow morphology surrounding the shock train is not monolithic, as is typically assumed, adjusting the balance of compressive forces within the shock cells. These findings lead to the conclusion that existing flux-conserved modeling approaches require modification to be effective in distorted and highly confined cases, including closure models that capture the three-dimensional distorted structure of the approach flow and its evolution along the shock train.

Shock Train↗

Navier-Stokes solutions with surface catalysis for Martian atmospheric entry

In this study numerical solutions have been obtained for two-dimensional axisymmetric hypersonic nonequilibrium CO2 flow over a high angle blunt cone with appropriate surface boundary conditions to account for energy and mass conservation at the body surface. The flowfield is described by the Navier-Stokes equations and multicomponent conservation laws which account for both translational and internal vibrational nonequilibrium effects. Complete forebody solutions have been obtained for the peak heating point of the Mars entry trajectory specified in the proposed NASA MESUR (Mars Environmental Survey) project. In these solutions, radiative equilibrium wall temperature and surface heating distributions are determined over the MESUR aeroshell forebody for entry velocity equal to 7 km/sec with varying degrees of surface catalysis. The effects of gas kinetics, surface catalysis, transport properties, and vibrational relaxation times on the surface heating are examined. The results identify some important issues in the prediction of surface heating for flows in thermochemical nonequilibrium and show that the Navier-Stokes code used herein is effective for thermal protection system design and materials selection.

Chen, Y.-K.↗

A hybrid multigrid technique for computing steady-state solutions to supersonic flows

Recently, Li and Sanders have introduced a class of finite difference schemes to approximate generally discontinuous solutions to hyperbolic systems of conservation laws. These equations have the form together with relevant boundary conditions. When modelling hypersonic spacecraft reentry, the differential equations above are frequently given by the compressible Euler equations coupled with a nonequilibrium chemistry model. For these applications, steady state solutions are often sought. Many tens (to hundreds) of super computer hours can be devoted to a single three space dimensional simulation. The primary difficulty is the inability to rapidly and reliably capture the steady state. In these notes, we demonstrate that a particular variant from the schemes presented can be combined with a particular multigrid approach to capture steady state solutions to the compressible Euler equations in one space dimension. We show that the rate of convergence to steady state coming from this multigrid implementation is vastly superior to the traditional approach of artificial time relaxation. Moreover, we demonstrate virtual grid independence. That is, the rate of convergence does not depend on the degree of spatial grid refinement.

Sanders, Richard↗

Numerical Analysis of Flow Evolution in a Helium Jet Injected into Ambient Air

A computational model to study the stability characteristics of an evolving buoyant helium gas jet in ambient air environment is presented. Numerical formulation incorporates a segregated approach to solve for the transport equations of helium mass fraction coupled with the conservation equations of mixture mass and momentum using a staggered grid method. The operating parameters correspond to the Reynolds number varying from 30 to 300 to demarcate the flow dynamics in oscillating and non-oscillating regimes. Computed velocity and concentration fields were used to analyze the flow structure in the evolving jet. For Re=300 case, results showed that an instability mode that sets in during the evolution process in Earth gravity is absent in zero gravity, signifying the importance of buoyancy. Though buoyancy initiates the instability, below a certain jet exit velocity, diffusion dominates the entrainment process to make the jet non-oscillatory as observed for the Re=30 case. Initiation of the instability was found to be dependent on the interaction of buoyancy and momentum forces along the jet shear layer.

Satti, Rajani P.↗

Choked flow effects in the NSI driven pin puller

This paper presents an analysis for pyrotechnic combustion and pin motion in the NASA Standard Initiator (NSI) actuated pin puller. The conservation principles and constitutive relations for a multi-phase system are posed and reduced to a set of eight ordinary differential equations which are solved to predict the system performance. The model tracks the interactions of the unreacted, incompressible solid pyrotechnic, incompressible condensed phase combustion products, and gas phase combustion products. The model accounts for multiple pyrotechnic grains, variable burn surface area, and combustion product mass flow rates through an orifice located within the device. Pressure-time predictions compare favorably with experimental data. Results showing model sensitivity to changes in the cross-sectional area of the orifice are presented.

Gonthier, Keith A.↗

Transonic potential flow and coordinate generation for bodies in a wind tunnel

An accurate method has been developed for computation of transonic potential flow about a 2-D lifting airfoil or an axisymmetric body in a wind tunnel. The computational mesh is nearly orthogonal everywhere and is generated by a sequence of Schwarz-Christoffel transformations and shearings to obtain an '0' grid near the body. A conservative finite-volume scheme for the full-potential equation and exact boundary conditions is used together with 'retarded' density to solve efficiently transonic flow with embedded shocks and large regions of supersonic flow. Supersonic free-stream flows are also solved with captured bow shocks and embedded subsonic regions.

Doria, M. L.↗

Effects on the plasmasphere of irregular electric fields

A conservative convection electric field model developed by Volland (1973) to describe the solar wind induced plasma flow within the inner magnetosphere is modified to include a noisy spatial component. Under steady state conditions such a random component will result in spatial irregularities in the thermal plasma density distribution in the vicinity of the plasmapause. Spatial irregularities in the convection can produce longitudinally restricted perturbations near the plasmapause, some of which are detached from the main body of the plasmasphere. Temporal variations in the midnight to noon flow intensity produce elongated extensions of the plasmasphere, but even short period variations of the overall magnitude of the convection cannot produce longitudinally localized perturbations in the thermal plasma distribution. Convection models based on the 3 hr magnetic index K sub p yield plasmasphere structures which are qualitatively similar to those based on shorter period variations, but the exact location at any given time of the plasmapause is dependent upon the characteristic time scale employed.

Grebowsky, J. M.↗

Lyman-alpha rocket spectra and models of the quiet and active solar chromosphere based on partial redistribution diagnostics

Absolute intensity Lyman-alpha profiles with a spatial resolution of 0.8 min and a spectral resolution of 50 mA were obtained for network and cell regions in the quiet sun, umbral and penumbral areas of a sunspot, two plages, and a quiescent prominence of the limb. Weak limb brightening shown by the Lyman-alpha cores and wings are consistent with predictions derived from partial redistribution line transfer calculations. Through use of a comoving two-level partial redistribution code which conserves mass flux, unequal red and blue Lyman-alpha peak intensities may be interpreted as flow velocities near 20,000 K. Outflows in the plages and downflows in the network are also noted. A model of a mean quiet sun chromosphere consistent with the Lyman-alpha integrated intensity, the Lyman continuum slope, and the millimeter continuum is presented.

Basri, G. S.↗

Conservative treatment of boundary interfaces for overlaid grids and multi-level grid adaptations

Conservative algorithms for boundary interfaces of overlaid grids are presented. The basic method is zeroth order, and is extended to a higher order method using interpolation and subcell decomposition. The present method, strictly based on a conservative constraint, is tested with overlaid grids for various applications of unsteady and steady supersonic inviscid flows with strong shock waves. The algorithm is also applied to a multi-level grid adaptation in which the next level finer grid is overlaid on the coarse base grid with an arbitrary orientation.

Moon, Young J.↗

Computational unsteady aerodynamics for lifting surfaces

Two dimensional problems are solved using numerical techniques. Navier-Stokes equations are studied both in the vorticity-stream function formulation which appears to be the optimal choice for two dimensional problems, using a storage approach, and in the velocity pressure formulation which minimizes the number of unknowns in three dimensional problems. Analysis shows that compact centered conservative second order schemes for the vorticity equation are the most robust for high Reynolds number flows. Serious difficulties remain in the choice of turbulent models, to keep reasonable CPU efficiency.

Edwards, John W.↗

Conservative treatment of boundary interfaces for overlaid grids and multi-level grid adaptations

Conservative algorithms for boundaray interfaces of overlaid grids are presented. The basic method is zeroth order, and is extended to a higher order method using interpolation and subcell decomposition. The present method, strictly based on a conservative constraint, is tested with overlaid grids for various applications of unsteady and steady supersonic inviscid flows with strong shock waves. The algorithm is also applied to a multi-level grid adaptation in which the next level finer grid is overlaid on the coarse base grid with an arbitrary orientation.

Moon, Young J.↗

Comparison of SMAC, PISO, and iterative time-advancing schemes for unsteady flows

Calculations of unsteady flows using a simplified marker and cell (SMAC), a pressure implicit splitting of operators (PSIO), and an iterative time advancing scheme (ITA) are presented. A partial differential equation for incremental pressure is used in each time advancing scheme. Example flows considered are a polar cavity flow starting from rest and self-sustained oscillating flows over a circular and a square cylinder. For a large time step size, the SMAC and ITA are more strongly convergent and yield more accurate results than PSIO. The SMAC is the most efficient computationally. For a small time step size, the three time advancing schemes yield equally accurate Strouhal numbers. The capability of each time advancing scheme to accurately resolve unsteady flows is attributed to the use of new pressure correction algorithm that can strongly enforce the conservation of mass. The numerical results show that the low frequency of the vortex shedding is caused by the growth time of each vortex shed into the wake region.

Kim, Sang-Wook↗

Comparison of the SMAC, PISO, and iterative time-advancing schemes for unsteady flows

Calculations of unsteady flows using a simplified marker and cell (SMAC), a pressure implicit splitting of operators (PISO), and an iterative time advancing scheme (ITA) are presented. A partial differential equation for incremental pressure is used in each time advancing scheme. Example flows considered are a polar cavity flow starting from rest and self-sustained oscillating flows over a circular and a square cylinder. For a large time step size, the SMAC and ITA are more strongly convergent and yield more accurate results than PISO. The SMAC is the most efficient computationally. For a small time step size, the three time advancing schemes yield equally accurate Strouhal numbers. The capability of each time advancing scheme to accurately resolve unsteady flows is attributed to the use of new pressure correction algorithm that can strongly enforce the conservation of mass. The numerical results show that the low frequency of the vortex shedding is caused by the growth time of each vortex shed into the wake region.

Kim, S.-W.↗

Calculating Shocks In Chemically Reacting Flows

Total-variation-diminishing dissipation operator results in crisp shocks. Report discusses method for numerical simulation of shocks in hypersonic, chemically reacting flows. Introduces mathematical model of gas of as many as three species for which species-conservation equations coupled with gas-dynamic equations.

Eberhardt, Scott↗