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.

SourceResearch connectionFrom yearThrough year

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

At least 631 records · Page 35

All sources NASA NTRS DOE OSTI

NASA NTRS1990

Upward ion flow in ionospheric holes on Venus

The nature of ion flow within two ionospheric holes on the nightside of Venus was investigated using ion composition measurements made by the ion mass spectrometer on the Pioneer Venus Orbiter. A comparison of the altitude profiles of the observed ion densities with those expected under diffusive equilibrium conditions indicates that the major ions O(+), NO(+), and O2(+), as well as the minor ions H(+) and He(+) flow upward and away from Venus along the axes of the holes. This result agrees with a quantitative evaluation of the ion flow speeds appearing in expressions derived from the equations for conservation of mass and momentum of the ions and electrons. The analysis shows that all ion species flow upward in the holes because the upward force produced by the plasma pressure gradient exceeds all downward forces. However, the nature of the ion source required to maintain such flow is not known.

Hartle, R. E.↗

NASA NTRS1990

High-resolution shock-capturing schemes for inviscid and viscous hypersonic flows

Hypersonic computations are presently conducted with an extension of a class of high-resolution implicit TVD algorithms suited to transonic multidimensional Euler and Navier-Stokes equations. These conservative shock-capturing schemes, which are spatially second- and third-order, may be first- and second-order accurate in time and suitable for either steady or unsteady calculations. Attention is given to the enhancement of hypersonic flows' convergence rate and stability; accuracy and efficiency is achieved by these means for very complex two-dimensional hypersonic viscous and inviscid shock interactions.

Yee, H. C.↗

NASA NTRS1991

Numerical investigation of an internal layer in turbulent flow over a curved hill

A numerical investigation of incompressible and compressible turbulent flows over strongly curved surfaces is presented. The turbulent flow equations are solved by a pressure based Navier-Stokes equations solver. In the method, the conservation of mass equation is replaced by a pressure correction equation applicable for both compressible and incompressible flows. The turbulence is described by a multiple time scale turbulence model supplemented with a near-wall turbulence model. The numerical results show that the internal layer is a strong turbulence field which is developed beneath the external boundary layer and is located very close to the wall. The development of the internal layer is attributed to the enormous mean flow strain rate caused by the streamline curvature. The external boundary layer flow responds rather slowly to the streamline curvature. Thus, the turbulence field of the forward corner of the curved hill is characterized by two turbulence fields interacting with each other. The turbulence intensity of the internal layer is much stronger than that of the external boundary layer, so that the development of a new boundary layer in the downstream region of the curved hill depends mostly on the internal layer. These numerical results are in good agreement with the measured data, and show that the turbulence model can resolve the turbulence field subjected to the strong streamline curvature.

Kim, S.-W.↗

NASA NTRS1991

Low thrust viscous nozzle flow field predictions

The low thrust viscous nozzle flows have been studied numerically by using the steady-state Reynolds-averaged Navier-Stokes equations in conservative form. A new Lower-Upper factorization scheme in the axisymmetric coordinate system has been developed and is analyzed in this paper. This new Navier-Stokes solver has shown that the numerical analog is accurate and stable for the convergent and divergent nozzle flowfield calculation. The range of the chamber pressure is (1.85-29.4)psi with the chamber temperature equal to 600 F. The results are compared with the data collected in the NASA/MSFC vacuum chamber, and they are in very reasonable agreement.

Chou, L. C.↗

NASA NTRS1992

Computer Program For Variable-Conductance Heat Pipes

VCHPDA provides accurate mathematical models of transient as well as steady-state performance of variable-conductance heat pipes over wide range of operating conditions. Applies to heat pipes with either cold, wicked or hot, nonwicked gas reservoirs and uses ideal-gas law and "flat-front" (negligible vapor diffusion) gas theory. Calculates length of gas-blocked region and temperature of vapor in active portion of heat pipe by solving set of nonlinear equations for conservation of energy and mass. Written in FORTRAN 77.

Antoniuk, D.↗

NASA NTRS1991

One-way nesting for a primitive equation ocean model

Prognostic numerical models for atmospheric and oceanic circulations require initial fields, boundary conditions, and forcing functions in addition to a consistent set of partial differential equations, including a state relation and equations expressing conservation of mass, momentum, and energy. Depending on the horizontal domain to be modeled, the horizontal boundary conditions are either physically obvious or extremely difficult to specify consistently. If the entire atmosphere is modeled, periodic horizontal boundary conditions are appropriate. On the other hand, the physical horizontal boundaries on the entire ocean are solid walls. Obviously, the normal velocity at a solid wall is zero while the specification of the tangential velocity depends on the mathematical treatment of the horizontal viscous terms. Limitations imposed by computer capacity and cost, as well as research interests, have led to the use of limited area models to study flows in the atmosphere and ocean. The limited area models do not have physical horizontal boundaries, merely numerical ones. Correctly determining these open boundary conditions for limited-area numerical models has both intrigued and frustrated numerical modelers for decades. One common approach is to use the closed or solid wall boundary conditions for a limited-area model. The argument given for this approach is that the boundary conditions affect flow near the walls but that none of these effects are propagated into the interior. Therefore, one chooses a big enough domain that the central region of interest is not corrupted by the boundary flow. Research in progress to model the North Atlantic circulation vividly illustrates the pitfalls of this approach. Two model runs are compared: (1) the southern boundary at 20S between latitudes 0 and 40W is artificially closed; and (2) the same boundary is specified as open with an inward transport of 15 Sv (determined from a global model with the same physics) uniformly spread across the boundary. A comparison of both runs is presented.

Blake, D. W.↗

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 NTRS1994

On the various forms of the energy equation for a dilute, monatomic mixture of nonreacting gases

In the case of gas mixtures, the governing equations become rather formidable and a complete listing of the equations in their various forms and methods to evaluate the transport coefficients is difficult to find. This paper seeks to compile common, as well as less well known, results in a single document. Various relationships between equations describing conservation of energy for a dilute, monatomic, nonreacting gas in local equilibrium are provided. The gas is treated as nonrelativistic, not subject to magnetic or electric fields, or radiative effects.

Kennedy, Christopher A.↗

NASA NTRS1995

Numerical simulations of small amplitude acoustic wave propagation in a converging-diverging nozzle

Numerical simulations of the small amplitude acoustic wave propagation in a converging-diverging nozzle proposed in Category 5 are performed and presented here. The quasi 1-D unsteady flow equations in conservative form are discretized by the DRP scheme with the artificial damping developed by Tam and Webb. Characteristic boundary conditions are used at both the inlet and outlet of the nozzle. The effect of different numerical implementations of the subsonic inflow boundary conditions on the convergence of the solution to the steady state is studied. In the case of a subsonic outflow in which a shock is formed in the nozzle, the interaction between acoustics and the shock is also investigated.

Dong, Thomas Z.↗

NASA NTRS1994

Numerical Simulation of the Instability of an Inviscid Liquid Jet in a Coflowing Compressible Airstream

The nonlinear interfacial instability of a liquid jet in a coflowing compressible airstream is studied numerically. A high-resolution scheme which has second-order accuracy in space and time is coupled with a Lagrangian marker particle algorithm to visualize the large-scale motion of the interfaces in compressible flow. A numerical algorithm based on an approximate equation of state of a compressible liquid is developed to allow this two-fluid system to be governed by the nonlinear unsteady Euler equations in conservative form. The initial growth of small disturbances given by the simulations agrees well with linear theory. The process of jet disruption in compressible flow is demonstrated to consist of the formation of liquid spikes, interweaving of the gas and liquid and stretching and detachment of the liquid main center core.

Li, H. -S.↗

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 NTRS2000

Solution of Third Computational Workshop Internal Propagation Problems using Low Order Schemes

The problems solved are propagation of sound waves through a transonic nozzle and shock-sound interaction. The first problem is solved using a MacCormack scheme. The second problem is solved by perturbation of a conservative Euler equation solution to the steady state. The problems are solved on evenly spaced grids. While solutions were found, the methods selected and not using stretched grids lead to solutions which do not compare well with those found using more accurate schemes. However, the perturbation of the mean flow scheme used to solve the second problem shows promise.

Miles, Jeffrey Hilton↗

NASA NTRS2001

Initial Circulation and Peak Vorticity Behavior of Vortices Shed from Airfoil Vortex Generators

An extensive parametric study of vortices shed from airfoil vortex generators has been conducted to determine the dependence of initial vortex circulation and peak vorticity on elements of the airfoil geometry and impinging flow conditions. These elements include the airfoil angle of attack, chord length, span, aspect ratio, local boundary layer thickness, and free stream Mach number. In addition, the influence of airfoil-to-airfoil spacing on the circulation and peak vorticity has been examined for pairs of co-rotating and counter-rotating vortices. The vortex generators were symmetric airfoils having a NACA-0012 cross-sectional profile. These airfoils were mounted either in isolation, or in pairs, on the surface of a straight pipe. The turbulent boundary layer thickness to pipe radius ratio was about 17 percent. The circulation and peak vorticity data were derived from cross-plane velocity measurements acquired with a seven-hole probe at one chord-length downstream of the airfoil trailing edge location. The circulation is observed to be proportional to the free-stream Mach number, the angle-of-attack, and the span-to-boundary layer thickness ratio. With these parameters held constant, the circulation is observed to fall off in monotonic fashion with increasing airfoil aspect ratio. The peak vorticity is also observed to be proportional to the free-stream Mach number, the airfoil angle-of-attack, and the span-to-boundary layer thickness ratio. Unlike circulation, however, the peak vorticity is observed to increase with increasing aspect ratio, reaching a peak value at an aspect ratio of about 2.0 before falling off again at higher values of aspect ratio. Co-rotating vortices shed from closely spaced pairs of airfoils have values of circulation and peak vorticity under those values found for vortices shed from isolated airfoils of the same geometry. Conversely, counter-rotating vortices show enhanced values of circulation and peak vorticity when compared to values obtained in isolation. The circulation may be accurately modeled with an expression based on Prandtl's relationship between finite airfoil circulation and airfoil geometry. A correlation for the peak vorticity has been derived from a conservation relationship equating the moment at the airfoil tip to the rate of angular momentum production of the shed vortex, modeled as a Lamb (ideal viscous) vortex. This technique provides excellent qualitative agreement to the observed behavior of peak vorticity for low aspect ratio airfoils typically used as vortex generators.

Wendt, Bruce J.↗

NASA NTRS2001

Vorticity and Wave Motion in a Compressible Protoplanetary Disk

The impact of an isolated vortex in a compressible Keplerian disk is examined using higher order numerical solutions of the Euler and entropy-conserving Energy equations. The vortex is stretched by the background shear flow with longer lasting anticyclonic vortices persisting for about 10 vortex revolutions. Simultaneously, the vortex emits transient radial waves consisting mainly of axisymmetrical weak shock waves and a slower, nonaxisymmetric Rossby wave. These waves may contribute to certain transient events in protoplanetary disks. The vortex stretching and waves were found to have little long-term feedback on the baseline 'standard solar nebula' disk structure and confirm the extremely stable structure of non self-gravitating disks.

Davis, Sanford S.↗

NASA NTRS1995

Heliospheric Termination Shock Motion Due to Fluctuations in the Solar Wind Upstream Conditions: Spherically Symmetric Model

Large-scale fluctuations in the solar wind plasma upstream of the heliospheric termination shock (TS) will cause inward and outward motions of the shock. Using numerical techniques, we extend an earlier strictly one-dimensional (planar) analytic gas dynamic model to spherical symmetry to investigate the features of global behavior of shock motion. Our starting point is to establish a steady numerical solution of the gasdynamic equations describing the interaction between the solar wind and the interstellar medium. We then introduce disturbances of the solar wind dynamic pressure at an inner boundary, and follow the subsequent evolution of the system, especially the motion of the termination shock. Our model solves spherically symmetric gasdynamic equations as an initial-boundary value problem. The equations in conservative form are solved using a fully implicit Total Variation Diminishing (TVD) upwind scheme with Roe-type Riemann solver. Boundary conditions are given by the solar wind parameters on an inner spherical boundary, where they are allowed to vary with time for unsteady calculations, and by a constant pressure (roughly simulating the effect of the local interstellar medium) on an outer boundary. We find that immediately after the interaction, the shock moves with speeds given by the earlier analogous analytic models. However, as the termination shock propagates it begins to slow down, seeking a new equilibrium position. In addition, the disturbance transmitted through the TS, either a shock or rarefaction wave, will encounter the heliopause boundary and be reflected back. The reflected signal will encounter the TS, causing it to oscillate. The phenomenon may be repeated for a number of reflections, resulting in a "ringing" of the outer heliosphere.

Ratkiewicz, R.↗

NASA NTRS2003

Modeling Variability of Plasma Conditions in the Io Torus

Telescopic observations an in situ measurements of the Io plasma torus show the density, temperature and composition to vary over time, sometimes up to a factor of 2. While previous models of the physical and chemical processes in the Io plasma torus have reasonably modeled the conditions of the Voyager 1 era, their authors have not addressed the observed variability nor explored the sensitivity of torus conditions to input parameters. In this paper we present a homogeneous torus model parameterized by five variables (transport timescale, neutral source strength, ratio of oxygen sulfur to atoms in the source, fraction of superthermal electrons, temperature of these hot electrons). The model incorporates the most recent data for ionization, recombination, charge exchange and radiative energy losses for the major torus species (S, S(sup +), S(sup ++), S(sup +++), O, O(sup +), O(sup ++). We solve equations of conservation of mass and energy to find equilibrium conditions for a set of input parameters. We compare model plasma conditions with those observed by Voyager 1 Voyager 2, and Cassini. Furthermore, we explore the sensitivity of torus conditions to each parameter. We find that (1) torus conditions are distinctly different for the Voyager 1, Voyager 2 and Cassini eras, (2) unique torus input parameters for any given era are poorly constrained given the wide range of solution space that is consistent with the range of observed torus conditions, (3) ion composition is highly sensitive to the specification of a non-thermal electron distribution, (4) neutral O/S source ratio is highly variable with model values ranging between 1.7 for Cassini to 4.0 for Voyager 1 conditions, (5) transport times range between 23 days for Voyager 2 to 50 days for Voyager 1 and Cassini, (6) neutral source strengths range between 7 to 30 x 10(sup -4) cm (sup -3) s(sup -1) which corresponds to a net production of 0.4 to 1.3 tons/s for a torus volume of 1.4 x 10(sup 31) cm(sup 3), or 38 R(sub j)(sup 3).

Delamere, P. A.↗

NASA NTRS2009

Evaluation of Margins of Safety in Brazed Joints

One of the essential steps in assuring reliable performance of high cost critical brazed structures is the assessment of the Margin of Safety (MS) of the brazed joints. In many cases the experimental determination of the failure loads by destructive testing of the brazed assembly is not practical and cost prohibitive. In such cases the evaluation of the MS is performed analytically by comparing the maximum design loads with the allowable ones and incorporating various safety or knock down factors imposed by the customer. Unfortunately, an industry standard methodology for the design and analysis of brazed joints has not been developed. This paper provides an example of an approach that was used to analyze an AlBeMet 162 (38%Be-62%Al) structure brazed with the AWS BAlSi-4 (Al-12%Si) filler metal. A practical and conservative interaction equation combining shear and tensile allowables was developed and validated to evaluate an acceptable (safe) combination of tensile and shear stresses acting in the brazed joint. These allowables are obtained from testing of standard tensile and lap shear brazed specimens. The proposed equation enables the assessment of the load carrying capability of complex brazed joints subjected to multi-axial loading.

Flom, Yury↗

NASA NTRS2007

Estimating the Inertia Matrix of a Spacecraft

A paper presents a method of utilizing some flight data, aboard a spacecraft that includes reaction wheels for attitude control, to estimate the inertia matrix of the spacecraft. The required data are digitized samples of (1) the spacecraft attitude in an inertial reference frame as measured, for example, by use of a star tracker and (2) speeds of rotation of the reaction wheels, the moments of inertia of which are deemed to be known. Starting from the classical equations for conservation of angular momentum of a rigid body, the inertia-matrix-estimation problem is formulated as a constrained least-squares minimization problem with explicit bounds on the inertia matrix incorporated as linear matrix inequalities. The explicit bounds reflect physical bounds on the inertia matrix and reduce the volume of data that must be processed to obtain a solution. The resulting minimization problem is a semidefinite optimization problem that can be solved efficiently, with guaranteed convergence to the global optimum, by use of readily available algorithms. In a test case involving a model attitude platform rotating on an air bearing, it is shown that, relative to a prior method, the present method produces better estimates from few data.

Acikmese, Behcet↗