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.

SourceResearch connectionFrom yearThrough year

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

At least 325 records · Page 18

All sources NASA NTRS DOE OSTI

NASA NTRS1976

The oblique impingement of an axisymmetric jet

The mechanics of the oblique impingement of an axisymmetric jet on a plane surface are examined in detail. The stagnation point is discussed. A schematic drawing of the problem and coordinate system used to describe the flow field are given. The kinematic features of the flow above the plate are examined in the context of the conservation of mass, the vorticity of the jet, and the vorticity introduced by the jetplate interaction. The dynamic features of the flow are examined in terms of the surface pressure distribution and the cause-effect relationships which exist between the pressure and velocity/vorticity distributions. Flow calculations performed are given. The investigation is relevant to the flow resulting from the interaction of the propulsion jet with the main airfoil (STOL aircraft), and is appropriate to an over- or under- wing configuration.

Foss, J. F.↗

NASA NTRS1998

Flow Simulation in Secondary Flow Passages of a Rocket Engine Turbopump

This paper presents application of a Generalized Fluid System Simulation Program (GFSSP) to model fluid flow in a very complex secondary flow circuit of a rocket engine turbopump, SIMPLEX. SIMPLEX was a low cost turbopump designed, manufactured and tested to demonstrate the ability to reduce the overall cost and design cycle time of turbomachinery. GFSSP is a general purpose computer program for analyzing flowrate, pressure, temperature and concentration distribution in a complex flow network. The program employs a finite volume formulation of mass, momentum and energy conservation equations in conjunction with thermodynamic equation of state of a real fluid. The secondary flow passages modeled include flow between the impeller shroud and housing, impeller back face and housing, flow through labyrinth seal and bearing, return line flow to inducer and flow adjacent to front and back face of the turbine. Axial load on the bearing are calculated from the predicted pressure on all axial faces of the rotor. The predicted pressure and temperature distributions compared well with the test data.

Schallhorn, Paul↗

NASA NTRS2014

A System of Conservative Regridding for Ice-Atmosphere Coupling in a General Circulation Model (GCM)

The method of elevation classes, in which the ice surface model is run at multiple elevations within each grid cell, has proven to be a useful way for a low-resolution atmosphere inside a general circulation model (GCM) to produce high-resolution downscaled surface mass balance fields for use in one-way studies coupling atmospheres and ice flow models. Past uses of elevation classes have failed to conserve mass and energy because the transformation used to regrid to the atmosphere was inconsistent with the transformation used to downscale to the ice model. This would cause problems for two-way coupling. A strategy that resolves this conservation issue has been designed and is presented here. The approach identifies three grids between which data must be regridded and five transformations between those grids required by a typical coupled atmosphere-ice flow model. This paper develops a theoretical framework for the problem and shows how each of these transformations may be achieved in a consistent, conservative manner. These transformations are implemented in Glint2, a library used to couple atmosphere models with ice models. Source code and documentation are available for download. Confounding real-world issues are discussed, including the use of projections for ice modeling, how to handle dynamically changing ice geometry, and modifications required for finite element ice models.

Ice-atmosphere↗

NASA NTRS1985

On convergence of computation of chemically reacting flows

The computational problems associated with high-temperature flows undergoing finite-rate ionization reactions is investigated. The conservation equations governing chemical species and vibrational and electron energies are solved simultaneously with those for overall mass, momentum, and energy for a one-dimensional subsonic flow, through a constant-area duct, originating behind a normal shock wave, using an implicit time-marching technique. Boundary conditions are imposed in the form of characteristic wave variables accounting for the effects of chemical reactions on the speed of sound. Converging solutions are obtained for cases in which chemical reactions are weak, but difficulty is encountered in other cases. The cause of the difficulty is investigated and shown to be the sharp pressure disturbances produced by such reactions.

Park, C.↗

NASA NTRS1981

Computation of wing-vortex interaction in transonic flow using implicit finite difference algorithm

An implicit delta form finite difference algorithm for Euler equations in conservation law form was used in preliminary calculations of three dimensional wing vortex interaction. Both steady and unsteady transonic flow wing vortex interactions are computed. The computations themselves are meant to guide upcoming wind tunnel experiments of the same flow field. Various modifications to the numerical method that are intended to improve computational efficiency are also described and tested in both two and three dimensions. Combination of these methods can reduce the overall computational time by a factor of 4.

Srinivasan, G.↗

NASA NTRS1998

Nuclear modules for space electric propulsion

Analysis of interplanetary cargo and piloted missions requires calculations of the performances and masses of subsystems to be integrated in a final design. In a preliminary and scoping stage the designer needs to evaluate options iteratively by using fast computer simulations. The Oak Ridge National Laboratory (ORNL) has been involved in the development of models and calculational procedures for the analysis (neutronic and thermal hydraulic) of power sources for nuclear electric propulsion. The nuclear modules will be integrated into the whole simulation of the nuclear electric propulsion system. The vehicles use either a Brayton direct-conversion cycle, using the heated helium from a NERVA-type reactor, or a potassium Rankine cycle, with the working fluid heated on the secondary side of a heat exchanger and lithium on the primary side coming from a fast reactor. Given a set of input conditions, the codes calculate composition. dimensions, volumes, and masses of the core, reflector, control system, pressure vessel, neutron and gamma shields, as well as the thermal hydraulic conditions of the coolant, clad and fuel. Input conditions are power, core life, pressure and temperature of the coolant at the inlet of the core, either the temperature of the coolant at the outlet of the core or the coolant mass flow and the fluences and integrated doses at the cargo area. Using state-of-the-art neutron cross sections and transport codes, a database was created for the neutronic performance of both reactor designs. The free parameters of the models are the moderator/fuel mass ratio for the NERVA reactor and the enrichment and the pitch of the lattice for the fast reactor. Reactivity and energy balance equations are simultaneously solved to find the reactor design. Thermalhydraulic conditions are calculated by solving the one-dimensional versions of the equations of conservation of mass, energy, and momentum with compressible flow.

Difilippo, F. C.↗

NASA NTRS2003

Finite Volume Methods: Foundation and Analysis

Finite volume methods are a class of discretization schemes that have proven highly successful in approximating the solution of a wide variety of conservation law systems. They are extensively used in fluid mechanics, porous media flow, meteorology, electromagnetics, models of biological processes, semi-conductor device simulation and many other engineering areas governed by conservative systems that can be written in integral control volume form. This article reviews elements of the foundation and analysis of modern finite volume methods. The primary advantages of these methods are numerical robustness through the obtention of discrete maximum (minimum) principles, applicability on very general unstructured meshes, and the intrinsic local conservation properties of the resulting schemes. Throughout this article, specific attention is given to scalar nonlinear hyperbolic conservation laws and the development of high order accurate schemes for discretizing them. A key tool in the design and analysis of finite volume schemes suitable for non-oscillatory discontinuity capturing is discrete maximum principle analysis. A number of building blocks used in the development of numerical schemes possessing local discrete maximum principles are reviewed in one and several space dimensions, e.g. monotone fluxes, E-fluxes, TVD discretization, non-oscillatory reconstruction, slope limiters, positive coefficient schemes, etc. When available, theoretical results concerning a priori and a posteriori error estimates are given. Further advanced topics are then considered such as high order time integration, discretization of diffusion terms and the extension to systems of nonlinear conservation laws.

Barth, Timothy↗

NASA NTRS1992

Unsteady transonic Euler solutions using finite elements

A finite element solution of the unsteady Euler equations is presented and demonstrated for 2D airfoil configurations oscillating in transonic flows. Computations are performed by spatially discretizing the conservation equations using the Galerkin weighted residual method and then employing a multistage Runge-Kutta scheme to march forward in time. A mesh deformation scheme has been developed to efficiently move interior points in a smooth fashion as the airfoil undergoes rigid body pitch and plunge motion. Both steady and unsteady results are presented, and a comparison is made with solutions obtained using finite-volume techniques. The effects of using either a lumped or consistent mass matrix are presented; the finite element method provides an accurate solution for unsteady transonic flows about isolated airfoils.

Davis, Gary A.↗

NASA NTRS1993

The combustion program at CTR

Understanding and modeling of turbulent combustion are key problems in the computation of numerous practical systems. Because of the lack of analytical theories in this field and of the difficulty of performing precise experiments, direct numerical simulation (DNS) appears to be one of the most attractive tools to use in addressing this problem. The general objective of DNS of reacting flows is to improve our knowledge of turbulent combustion but also to use this information for turbulent combustion models. For the foreseeable future, numerical simulation of the full three-dimensional governing partial differential equations with variable density and transport properties as well as complex chemistry will remain intractable; thus, various levels of simplification will remain necessary. On one hand, the requirement to simplify is not necessarily a handicap: numerical simulations allow the researcher a degree of control in isolating specific physical phenomena that is inaccessible in experiments. CTR has pursued an intensive research program in the field of DNS for turbulent reacting flows since 1987. DNS of reacting flows is quite different from DNS of non-reacting flows: without reaction, the equations to solve are clearly the five conservation equations of the Navier Stokes system for compressible situations (four for incompressible cases), and the limitation of the approach is the Reynolds number (or in other words the number of points in the computation). For reacting flows, the choice of the equations, the species (each species will require one additional conservation equation), the chemical scheme, and the configuration itself is more complex.

Poinsot, Thierry J.↗

NASA NTRS1994

Modeling jets in cross flow

Various approaches to the modeling of jets in cross flow are reviewed. These are grouped into four classes, namely: empirical models, integral models, perturbation models, and numerical models. Empirical models depend largely on the correlation of experimental data and are mostly useful for first-order estimates of global properties such as jet trajectory and velocity and temperature decay rates. Integral models are based on some ordinary-differential form of the conservation laws, but require substantial empirical calibration. They allow more details of the flow field to be obtained; simpler versions have to assume similarity of velocity and temperature profiles, but more sophisticated ones can actually calculate these profiles. Perturbation models require little empirical input, but the need for small parameters to ensure convergent expansions limits their application to either the near-field or the far-field. Therefore, they are mostly useful for the study of flow physics. Numerical models are based on conservation laws in partial-differential form. They require little empirical input and have the widest range of applicability. They also require the most computational resources. Although many qualitative and quantitative features of jets in cross flow have been predicted with numerical models, many issues affecting accuracy such as grid resolution and turbulence model are not completely resolved.

Demuren, A. O.↗

NASA NTRS1976

A two-region model of the solar wind including azimuthal velocity

The two-region model of the solar wind divides the interplanetary space into two regions: it assumes that the solar wind is one-fluid in an inner region within 0.4 AU and two-fluid in an outer region beyond 0.4 AU. This paper includes the angular motion of the solar wind in the two-region model. The flow in the one-fluid region is governed by the one-fluid magnetohydrodynamic equations. The second and third moment equations of the Vlasov equation together with other conservation equations are used to describe the solar-wind flow in the two-fluid region. The predicted azimuthal velocity at 1 AU is less than 2 km/s. All other macroscopic and microscopic properties from this model are in good agreement with experimental quiet-time observations at 1 AU. The numerical results also confirm that when the azimuthal velocity is included in the analysis, the amount of magnetic-field energy converted into kinetic energy in the solar wind is only a small fraction of the total expansion energy flux and has little effect upon the final radial expansion velocity.

Acuna, M. H.↗

NASA NTRS1985

Three dimensional CSCM method for the compressible Navier-Stokes equations with application to a multi-nozzle exhaust flowfield

The effective design of three-dimensional flow systems such as a clustered multinozzle shrouded exhaust system requires detailed analysis of the internal flowfield. It is pointed out that new numerical techniques on the most modern vector computers are making these analyses practicable. In this paper, the upwind differenced implicit conservative supracharacteristics method (CSCM) for three-dimensional viscous flow is presented. The new method is an implicit 'method of planes' symmetric Gauss-Seidel relaxation scheme. The data is conveniently stored on successive planes along the streamwise coordinate, and the system of equations is solved twice in each successive plane of the streamwise coordinate. The developed new three-dimensional compressible Navier-Stokes algorithm combines the best features of the storage and computationally efficient space marching procedures with the generality of time dependent techniques to solve flows with elliptic and streamwise separated flows.

Bardina, J.↗

NASA NTRS1986

Three-dimensional, conservative, Euler computations using patched grid systems and explicit methods

The method of 'zonal approach' (in which the flow field is partitioned into regions with independent grids) for computation of flow over complex geometries, such as aircraft configuration, requires application of a grid-interfacing procedure. A three-dimensional conservative, boundary scheme for patched grids, applicable in generalized coordinates for arbitrary point distributions on a planar zonal surface is presented. The computation technique is derived within the framework of the Osher (1983) upwind scheme, using Euler equations. The three-dimensional interfacing method is applied to the computation of flow about a wing-canard combination using a two-zone, patched grid.

Hessenius, K. A.↗

NASA NTRS1987

SHARP simulation of discontinuities in highly convective steady flow

For steady multidimesional convection, the Quadratic Upstream Interpolation for Convective Kinematics (QUICK) scheme has several attractive properties. However, for highly convective simulation of step profiles, QUICK produces unphysical overshoots and a few oscillations, and this may cause serious problems in nonlinear flows. Fortunately, it is possible to modify the convective flux by writing the normalized convected control-volume face value as a function of the normalized adjacent upstream node value, developing criteria for monotonic resolution without sacrificing formal accuracy. This results in a nonlinear functional relationship between the normalized variables, whereas standard methods are all linear in this sense. The resulting Simple High Accuracy Resolution Program (SHARP) can be applied to steady multidimensional flows containing thin shear or mixing layers, shock waves, and other frontal phenomena. This represents a significant advance in modeling highly convective flows of engineering and geophysical importance. SHARP is based on an explicit, conservative, control-volume flux formation, equally applicable to one, two, or three dimensional elliptic, parabolic, hyperbolic, or mixed-flow regimes. Results are given for the bench-mark purely convective first-order results and the nonmonotonic predictions of second- and third-order upwinding.

Leonard, B. P.↗

NASA NTRS1988

Simple high-accuracy resolution program for convective modelling of discontinuities

For steady multidimensional convection, the Quadratic Upstream Interpolation for Convective Kinematics (QUICK) scheme has several attractive properties. However, for highly convective simulation of step profiles, QUICK produces unphysical overshoots and a few oscillations, and this may cause serious problems in nonlinear flows. Fortunately, it is possible to modify the convective flux by writing the normalized convected control-volume face value as a function of the normalized adjacent upstream node value, developing criteria for monotonic resolution without sacrificing formal accuracy. This results in a nonlinear functional relationship between the normalized variables, whereas standard methods are all linear in this sense. The resulting Simple High Accuracy Resolution Program (SHARP) can be applied to steady multidimensional flows containing thin shear or mixing layers, shock waves, and other frontal phenomena. This represents a significant advance in modeling highly convective flows of engineering and geophysical importance. SHARP is based on an explicit, conservative, control-volume flux formation, equally applicable to one, two, or three dimensional elliptic, parabolic, hyperbolic, or mixed-flow regimes. Results are given for the bench-mark purely convective first-order results and the nonmonotonic predictions of second- and third-order upwinding.

Leonard, B. P.↗

NASA NTRS1980

Nonlinear axially symmetric circulations in a nearly inviscid atmosphere

The structure of certain axially symmetric circulations in a stably stratified, differentially heated, rotating Boussinesq fluid on a sphere is analyzed. A simple approximate theory (similar to that introduced by Schneider (1977)) is developed for the case in which the fluid is sufficiently inviscid that the poleward flow in the Hadley cell is nearly angular momentum conserving. The theory predicts the width of the Hadley cell, the total poleward heat flux, the latitude of the upper level jet in the zonal wind, and the distribution of surface easterlies and westerlies. Fundamental differences between such nearly inviscid circulations and the more commonly studied viscous axisymmetric flows are emphasized. The theory is checked against numerical solutions to the model equations.

Held, I. M.↗

NASA NTRS1993

Light ion flow in the nightside ionosphere of Venus

The flow characteristics of the light ions H(+) and He(+) have been studied in the midnight region of the ionosphere of Venus. Measurements of ion composition, electron and ion temperatures and magnetic fields by instruments onboard the Pioneer Venus Orbiter have been used in rite electron and ion equations of conservation of mass and momentum to derive the vertical flow velocities of H(+) and He(+). When average height profiles of the measured quantities were used, H(+) was found to flow upward, accelerating to speeds of almost 1 km/s at the ion-exobase. In a similar fashion, He(+) was found to flow downward into the neutral atmosphere where it is readily quenched by charge transfer reactions. The polarization electric field played an important role in forcing H(+) upward, but did not contribute enough to the He(+) force balance to produce upward flow. At the ion-exobase, the outward electric polarization force on H(+) was shown to be five times the gravitational force. Using an analogy with the terrestrial ion-exosphere, H(+) was inferred to flow upward into the ionotail of Venus and accelerate to escape speeds. A planet averaged escape flux of 1.4 x 10 exp 7/sq cm/s was calculated, which is comparable to hydrogen loss rates estimated by other investigators.

Hartle, R. E.↗

NASA NTRS1990

Numerical simulation of steady and unsteady asymmetric vortical flow

The unsteady, compressible, thin-layer, Navier-Stokes equations are solved to simulate steady and unsteady, asymmetric, vortical laminar flow around cones at high incidences and supersonic Mach numbers. The equations are solved by using an implicit, upwind, flux-difference splitting, finite-volume scheme. Locally conical flows are assumed and the solutions are obtained by forcing the conserved components of the flowfield vector to be equal at two axial stations located at 0.95 and 1.0. Computational examples cover steady and unsteady asymmetric flows around a circular cone and its control using side strakes. Unsteady asymmetric flows are also presented for elliptic- and diamond-section cones, which model asymmetric vortex shedding around round and sharp-edged delta wing flows.

Kandil, O. A.↗