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 883 records · Page 49

Quasi 1-D Analysis of a Circular, Compressible, Turbulent Jet Laden with Water Droplets

Recent experimental studies indicate that presence of small amount of liquid droplets reduces the Overall Sound Pressure Level (OASPL) of a jet. Present study is aimed at numerically investigating the effect of liquid particles on the overall flow quantities of a heated, compressible round jet. The jet is assumed perfectly expanded. A quasi-1D model was developed for this purpose which uses area-averaged quantities that satisfy integral conservation equations. Special attention is given to represent the early development region since it is acoustically important. Approximate velocity and temperature profiles were assumed in this region to evaluate entrainment rate. Experimental correlations were used to obtain spreading rate of shear layer. The base flow thus obtained is then laden with water droplets at the exit of the nozzle. Mass, momentum and energy coupling between the two phases is represented using empirical relations. Droplet size and mass loading are varied to observe their effect on flow variables.

Source record↗

A longitudinally-patched grid approach with applications to high speed flows

The combined use of an implicit, upwind finite-volume scheme with an 'equivalent' gamma formulation (for real gas effects) and with a three-dimensional conservative patched grid formulation is discussed. Results are presented over a wide Mach number range on both single and patched grids.

Walters, R. W.↗

Pulsed electromagnetic gas acceleration

Experimental data were combined with one-dimensional conservation relations to yield information on the energy deposition ratio in a parallel-plate accelerator, where the downstream flow was confined to a constant area channel. Approximately 70% of the total input power was detected in the exhaust flow, of which only about 20% appeared as directed kinetic energy, thus implying that a downstream expansion to convert chamber enthalpy into kinetic energy must be an important aspect of conventional high power MPD arcs. Spectroscopic experiments on a quasi-steady MPD argon accelerator verified the presence of A(III) and the absence of A(I), and indicated an azimuthal structure in the jet related to the mass injection locations. Measurements of pressure in the arc chamber and impact pressure in the exhaust jet using a piezocrystal backed by a Plexiglas rod were in good agreement with the electromagnetic thrust model.

Jahn, R. G.↗

Blood flow and oxygen uptake during exercise

A model is developed for predicting oxygen uptake, muscle blood flow, and blood chemistry changes under exercise conditions. In this model, the working muscle mass system is analyzed. The conservation of matter principle is applied to the oxygen in a unit mass of working muscle under transient exercise conditions. This principle is used to relate the inflow of oxygen carried with the blood to the outflow carried with blood, the rate of change of oxygen stored in the muscle myoglobin, and the uptake by the muscle. Standard blood chemistry relations are incorporated to evaluate venous levels of oxygen, pH, and carbon dioxide.

Mitchell, J. W.↗

Numerical computation of transonic flows with shock waves

Some recent developments in numerical methods for calculating solutions to the transonic potential flow equation are reviewed, including (1) the construction of the stable coordinate independent difference schemes; (2) the use of conservation form to insure proper shock jump conditions; (3) analysis of the relaxation method by the time dependent analogy; (4) accelerated iterative schemes.

Jameson, A.↗

An alternative approach to the numerical simulation of steady inviscid flow

A numerical procedure for the efficient simulation of steady inviscid flow is described and its utility demonstrated. Application of the surrogate equation technique allows the formulation of stable, fully conservative, type dependent finite difference equations for use in obtaining numerical solutions to systems of first order partial differential equations, such as the steady state Euler equations or their various approximations. Computational results are presented for the full Euler equations and for the transonic disturbance equations. For the latter case, a computational efficiency greater than that obtained by means of the standard perturbation potential approach is indicated.

Johnson, G. M.↗

Solution of the two-dimensional Euler equations with generalized coordinate transformation using flux vector splitting

An implicit finite difference code using flux vector splitting has been developed for solving the two-dimensional inviscid gas dynamics equations. The method is spatially second-order acurate, fully conservative, and uses body-conforming generalized coordinates for treating complex geometries. Numerical results have been obtained for transonic flow over a circular cylinder and airfoils. Steady results for a half cylinder (top and bottom symmetry-imposed) range from critical flow to a strong shock case with rotationally induced flow separation. Full cylinder solutions at freestream Mach number values of 0.5, however, show unsteady oscillation. A perturbation form of the method has also been developed and used to compute both fore and aft inviscid flow separation about a cylinder for a nonuniform incoming stream.

Buning, P. G.↗

Transonic viscous and inviscid solutions using a composite velocity procedure

Solutions for transonic viscous and inviscid flows using a composite velocity procedure are presented. The velocity components of the compressible Navier-Stokes equations are written in terms of a multiplicative composite consisting of a viscous or rotational velocity and an inviscid, irrotational, potential like function. For viscous flows, the governing equations are solved using a coupled strongly implicit procedure and the outer inviscid flow corresponds to an irrotational, potential flow. Further investigation of the inviscid flow solution is made so that the correct rotational flow solution may be obtained for the transonic case. The Euler equations are solved for the inviscid flow calculations. Solutions to the full potential equation are recovered unless the axial momentum equation is written in full conservation form.

Gordnier, R. E.↗

Stability of stationary barotropic modons by Lyapunov's direct method

A new Liapunov stability condition is formulated for the shallow-water equations, using a gage-variable formalism. This sufficient condition is derived for the class of perturbations that conserve the total mass. It is weaker than existing stability criteria, i.e., it applies to a wider class of flows. Formal stability to infinitesimally small perturbations of arbitrary shape is obtained for two classes of large-scale geophysical flows: pseudo-eastward flow with constant shear, and localized coherent structures of modon type.

Sakuma, H.↗

A diagonal implicit scheme for computing flows with finite-rate chemistry

A new algorithm for solving steady, finite-rate chemistry, flow problems is presented. The new scheme eliminates the expense of inverting large block matrices that arise when species conservation equations are introduced. The source Jacobian matrix is replaced by a diagonal matrix which is tailored to account for the fastest reactions in the chemical system. A point-implicit procedure is discussed and then the algorithm is included into the LU-SGS scheme. Solutions are presented for hypervelocity reentry and Hydrogen-Oxygen combustion. For the LU-SGS scheme a CFL number in excess of 10,000 has been achieved.

Eberhardt, Scott↗

Validation of the space shuttle main engine steady state performance model

The primary objective was to present methods for validating predictions of Rocketdyne's most current version of the Space Shuttle Main Engine (SSME) Power Balance Model (PBM) with respect to physical relations governing flow systems. This required the development and implementation of postprocessors to check results of PBM computations for satisfaction of conservation relations. A cursory uncertainty analysis of PBM predictions with respect to mass and energy balances was performed. In addition, an effort to identify the empirical relations and physical assumptions within PBM which impact the ability of the model to attain rigorous balance was begun.

Santi, L. Michael↗

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.↗

Modeling of InP metalorganic chemical vapor deposition

The growth of InP by metalorganic chemical vapor deposition (MOCVD) in a horizontal reactor is being modeled with a commercially available computational fluid dynamics modeling code. The mathematical treatment of the MOCVD process has four primary areas of concern: 1) transport phenomena, 2) chemistry, 3) boundary conditions, and 4) numerical solution methods. The transport processes involved in CVD are described by conservation of total mass, momentum, energy, and atomic species. Momentum conservation is described by a generalized form of the Navier-Stokes equation for a Newtonian fluid and laminar flow. The effect of Soret diffusion on the transport of particular chemical species and on the predicted deposition rate is examined. Both gas-phase and surface chemical reactions are employed in the model. Boundary conditions are specified at the inlet and walls of the reactor for temperature, fluid flow and chemical species. The coupled set of equations described above is solved by a finite difference method over a nonuniform rectilinear grid in both two and three dimensions. The results of the 2-D computational model is presented for gravity levels of zero- and one-g. The predicted growth rates at one-g are compared to measured growth rates on fused silica substrates.

Black, Linda R.↗

Design optimization of a brush turbine with a cleaner/water based solution

Recently, a fluid turbine which has a brush attached to it has been designed and tested with water as fluid. The purpose of the turbine-brush is to clean up fouling in a tube. The Montreal Protocol prohibits the use of CFC products from refrigeration industry or from industry in general as a cleanser in 1996. Alternatives for the cleansers, devices or a combination of alternative devices with a cleanser should be found. One of the methods is to develop a device which cleans fouling with a cleaning medium. In this paper, we describe a turbine connected with a brush. However, the turbine with the brush should be simple and easy to install. This device is a combined small liquid turbine with a brush. The turbine is activated by the liquid flowing through the tube. Then the turbine turns the brush cleaning fouling along the tube. Based on the energy conservation and the Bernoulli equation along with an empirical relationship of drag force obtained from an experimental apparatus, a relationship of the rotational speed, the number of blades, and geometric variables of the turbine-brush was obtained. The predicted rotational speeds were compared with the experimental observations. Further work was recommended for improvements.

Kim, Rhyn H.↗

Effects of spatial order of accuracy on the computation of vortical flowfields

The objective of the present investigation is to assess the effect of the spatial order of accuracy used for the evaluation of the inviscid fluxes on the resolution of higher order quantitites, such as velocity gradients. The viscous terms are computed as second-order accurate with central difference formulas, even though for the explicit part of the algorithm higher order approximations may be used. A viscous/inviscid method is used, and the outer part of the flowfield is computed with the inviscid flow equations. The viscous boundary-layer type flow region close to the body surface is computed with an algebraic eddy viscosity model. Results obtained with the conservative and nonconservative formulations and the viscous/inviscid approach are compared with available experimental data. The effect of grid refinement on the accuracy of the solution is also presented.

Ekaterinaris, J. A.↗

Efficient implementation of weighted ENO schemes

In this paper, we further analyze, test, modify and improve the high order WENO (weighted essentially non-oscillatory) finite difference schemes of Liu, Osher and Chan. It was shown by Liu et al. that WENO schemes constructed from the r-th order (in L1 norm) ENO schemes are (r+1)-th order accurate. We propose a new way of measuring the smoothness of a numerical solution, emulating the idea of minimizing the total variation of the approximation, which results in a 5-th order WENO scheme for the case r = 3, instead of the 4-th order with the original smoothness measurement by Liu et al. This 5-th order WENO scheme is as fast as the 4-th order WENO scheme of Liu et al., and both schemes are about twice as fast as the 4-th order ENO schemes on vector supercomputers and as fast on serial and parallel computers. For Euler systems of gas dynamics, we suggest computing the weights from pressure and entropy instead of the characteristic values to simplify the costly characteristic procedure. The resulting WENO schemes are about twice as fast as the WENO schemes using the characteristic decompositions to compute weights, and work well for problems which do not contain strong shocks or strong reflected waves. We also prove that, for conservation laws with smooth solutions, all WENO schemes are convergent. Many numerical tests, including the 1D steady state nozzle flow problem and 2D shock entropy wave interaction problem, are presented to demonstrate the remarkable capability of the WENO schemes, especially the WENO scheme using the new smoothness measurement, in resolving complicated shock and flow structures. We have also applied Yang's artificial compression method to the WENO schemes to sharpen contact discontinuities.

Jiang, Guang-Shan↗

Large-eddy simulation of flow past a circular cylinder

Some of the most challenging applications of large-eddy simulation are those in complex geometries where spectral methods are of limited use. For such applications more conventional methods such as finite difference or finite element have to be used. However, it has become clear in recent years that dissipative numerical schemes which are routinely used in viscous flow simulations are not good candidates for use in LES of turbulent flows. Except in cases where the flow is extremely well resolved, it has been found that upwind schemes tend to damp out a significant portion of the small scales that can be resolved on the grid. Furthermore, it has been found that even specially designed higher-order upwind schemes that have been used successfully in the direct numerical simulation of turbulent flows produce too much dissipation when used in conjunction with large-eddy simulation. The objective of the current study is to perform a LES of incompressible flow past a circular cylinder at a Reynolds number of 3900 using a solver which employs an energy-conservative second-order central difference scheme for spatial discretization and compare the results obtained with those of Beaudan & Moin (1994) and with the experiments in order to assess the performance of the central scheme for this relatively complex geometry.

Mittal, R.↗

Computing Axisymmetric Jet Screech Tones using Unstructured Grids

The purpose of this paper is to show that computations with an aeroacoustic feedback loop, the jet screech noise, can be obtained using truly unstructured grid technology. Numerical results are presented for a nozzle with two different lip thicknesses which will be referred to in this paper as a thin and a thick lip nozzle respectively. The space-time conservation element and solution element (CE/SE) method is used to solve the conservation laws of the compressible axisymmetric Navier-Stokes equations. The equations are time marched to predict the unsteady flow and the near-field screech tone noise issuing from an underexpanded circular jet. The CE/SE method uses an unstructured grid based data structure. The unstructured grids for these calculations are generated based on the method of Delaunay triangulation. Comparisons of numerical results with available experimental data are shown for flows corresponding to several different jet Mach numbers. Generally good agreement is obtained in terms of flow physics, screech tone frequency, and sound pressure level.

Jorgenson, Philip C. E.↗