SEARCH · Engineering Papers

Results for “conservative”

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 721 records · Page 40

All sources NASA NTRS DOE OSTI

NASA NTRS1989

Grid generation and inviscid flow computation about aircraft geometries

Grid generation and Euler flow about fighter aircraft are described. A fighter aircraft geometry is specified by an area ruled fuselage with an internal duct, cranked delta wing or strake/wing combinations, canard and/or horizontal tail surfaces, and vertical tail surfaces. The initial step before grid generation and flow computation is the determination of a suitable grid topology. The external grid topology that has been applied is called a dual-block topology which is a patched C (exp 1) continuous multiple-block system where inner blocks cover the highly-swept part of a cranked wing or strake, rearward inner-part of the wing, and tail components. Outer-blocks cover the remainder of the fuselage, outer-part of the wing, canards and extend to the far field boundaries. The grid generation is based on transfinite interpolation with Lagrangian blending functions. This procedure has been applied to the Langley experimental fighter configuration and a modified F-18 configuration. Supersonic flow between Mach 1.3 and 2.5 and angles of attack between 0 degrees and 10 degrees have been computed with associated Euler solvers based on the finite-volume approach. When coupling geometric details such as boundary layer diverter regions, duct regions with inlets and outlets, or slots with the general external grid, imposing C (exp 1) continuity can be extremely tedious. The approach taken here is to patch blocks together at common interfaces where there is no grid continuity, but enforce conservation in the finite-volume solution. The key to this technique is how to obtain the information required for a conservative interface. The Ramshaw technique which automates the computation of proportional areas of two overlapping grids on a planar surface and is suitable for coding was used. Researchers generated internal duct grids for the Langley experimental fighter configuration independent of the external grid topology, with a conservative interface at the inlet and outlet.

Smith, Robert E.↗

NASA NTRS1991

The P1-RKDG method for two-dimensional Euler equations of gas dynamics

A class of nonlinearly stable Runge-Kutta local projection discontinuous Galerkin (RKDG) finite element methods for conservation laws is investigated. Two dimensional Euler equations for gas dynamics are solved using P1 elements. The generalization of the local projections, which for scalar nonlinear conservation laws was designed to satisfy a local maximum principle, to systems of conservation laws such as the Euler equations of gas dynamics using local characteristic decompositions is discussed. Numerical examples include the standard regular shock reflection problem, the forward facing step problem, and the double Mach reflection problem. These preliminary numerical examples are chosen to show the capacity of the approach to obtain nonlinearly stable results comparable with the modern nonoscillatory finite difference methods.

Cockburn, Bernardo↗

NASA NTRS1991

Recent improvements in the nonequilibrium VSL scheme for hypersonic blunt-body flows

The nonequilibrium viscous shock-layer (VSL) solution scheme is revisited to improve its solution accuracy in the stagnation point region and also to minimize and control the errors in the conservation of elemental mass. The stagnation-point solution is improved by using a second-order expansion for the normal velocity and the elemental mass conservation is improved by directly imposing the element conservation equations as solution constraints. These modifications are such that the general structure and computational efficiency of the nonequilibrium VSL scheme is not affected. This revised nonequilibrium VSL scheme is used to study the Mach 20 flow over a 7-deg sphere-cone vehicle under zero and 20-deg angle-of-attack conditons. Comparisons are made with the corresponding predictions of Navier-Stokes and Parabolized Navier-Stokes solution schemes. The results of these tests show that the nonequilibrium blunt-body VSL scheme is indeed an accurate, fast and extremely efficient means for generating the blunt-body flowfield over spherical nosetips at small-to-large angles of attack.

Bhutta, Bilal A.↗

NASA NTRS1991

Modal interaction in linear dynamic systems near degenerate modes

In various problems in structural dynamics, the eigenvalues of a linear system depend on a characteristic parameter of the system. Under certain conditions, two eigenvalues of the system approach each other as the characteristic parameter is varied, leading to modal interaction. In a system with conservative coupling, the two eigenvalues eventually repel each other, leading to the curve veering effect. In a system with nonconservative coupling, the eigenvalues continue to attract each other, eventually colliding, leading to eigenvalue degeneracy. Modal interaction is studied in linear systems with conservative and nonconservative coupling using singularity theory, sometimes known as catastrophe theory. The main result is this: eigenvalue degeneracy is a cause of instability; in systems with conservative coupling, it induces only geometric instability, whereas in systems with nonconservative coupling, eigenvalue degeneracy induces both geometric and elastic instability. Illustrative examples of mechanical systems are given.

Afolabi, D.↗

NASA NTRS1992

On recontamination and directional-bias problems in Monte Carlo simulation of PDF turbulence models

Turbulent combustion can not be simulated adequately by conventional moment closure turbulent models. The probability density function (PDF) method offers an attractive alternative: in a PDF model, the chemical source terms are closed and do not require additional models. Because the number of computational operations grows only linearly in the Monte Carlo scheme, it is chosen over finite differencing schemes. A grid dependent Monte Carlo scheme following J.Y. Chen and W. Kollmann has been studied in the present work. It was found that in order to conserve the mass fractions absolutely, one needs to add further restrictions to the scheme, namely alpha(sub j) + gamma(sub j) = alpha(sub j - 1) + gamma(sub j + 1). A new algorithm was devised that satisfied this restriction in the case of pure diffusion or uniform flow problems. Using examples, it is shown that absolute conservation can be achieved. Although for non-uniform flows absolute conservation seems impossible, the present scheme has reduced the error considerably.

Hsu, Andrew T.↗

NASA NTRS1992

Improved nonequilibrium viscous shock-layer scheme for hypersonic blunt-body flowfields

The nonequilibrium viscous shock-layer (VSL) solution scheme is revisited to improve its solution accuracy in the stagnation region and also to minimize and control errors in the conservation of elemental mass. The stagnation-point solution is improved by using a second-order expansion for the normal velocity, and the elemental mass conservation is improved by directly imposing the element conservation equations as solution constraints. These modifications are such that the general structure and computational efficiency of the nonequilibrium VSL scheme is not affected. This revised nonequilibrium VSL scheme is used to study the Mach 20 flow over a 7-deg sphere-cone vehicle under 0- and 20-deg angle-of-attack conditions. Comparisons are made with the corresponding predictions of Navier-Stokes and parabolized Navier-Stokes solution schemes. The results of these tests show that the nonequilibrium blunt-body VSL scheme is indeed an accurate, fast, and extremely efficient means for generating the blunt-body flowfield over spherical nose tips at zero-to-large angles of attack.

Bhutta, Bilal A.↗

NASA NTRS1992

Theory of magnetohydrodynamic waves - The WKB approximation revisited

The present paper revisits the derivation of the WKB approximation for small-amplitude magnetohydrodynamic waves, allowing for possible spatial variation of the background entropy. The equation of variation of wave amplitude is rederived; it is a bilinear equation which, it turns out, can be recast in the action conservation form. It is shown that this action conservation equation is in fact equivalent to the action conservation law obtained from Lagrangian treatments.

Barnes, Aaron↗

NASA NTRS1993

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

The effect of the order-of-accuracy, used for the spatial discretization, on the resolution of the leading edge vortices over sharp-edged delta wings is investigated. The flowfield is computed using a viscous/inviscid zonal approach. The viscous flow in the vicinity of the wing is computed using the conservative formulation of the compressible, thin-layer Navier-Stokes equations. The leeward-side vortical flowfield and the other flow regions away from the surface are computed as inviscid. The time integration is performed with both an explicit fourth-order Runge-Kutta scheme and an implicit, factorized, iterative scheme. High-order-accurate inviscid fluxes are computed using both a conservative and a non-conservative (primitive variable) formulation. The nonlinear, inviscid terms of the primitive variable form of the governing equations are evaluated with a finite-difference numerical scheme based on the sign of the eigenvalues. High-order, upwind-biased, finite difference formulas are used to evaluate the derivatives of the nonlinear convective terms. Computed results are compared with available experimental data, and comparisons of the flowfield in the vicinity of the vortex cores are presented.

Ekaterinaris, J. A.↗

NASA NTRS1993

Flux-Based Finite Volume representations for general thermal problems

Flux-Based Finite Volume (FV) element representations for general thermal problems are given in conjunction with a generalized trapezoidal gamma-T family of algorithms, formulated following the spirit of what we term as the Lax-Wendroff based FV formulations. The new flux-based representations introduced offer an improved physical interpretation of the problem along with computationally convenient and attractive features. The space and time discretization emanate from a conservation form of the governing equation for thermal problems, and in conjunction with the flux-based element representations give rise to a physically improved and locally conservative numerical formulations. The present representations seek to involve improved locally conservative properties, improved physical representations and computational features; these are based on a 2D, bilinear FV element and can be extended for other cases. Time discretization based on a gamma-T family of algorithms in the spirit of a Lax-Wendroff based FV formulations are employed. Numerical examples involving linear/nonlinear steady and transient situations are shown to demonstrate the applicability of the present representations for thermal analysis situations.

Mohan, Ram V.↗

NASA NTRS1954

Analytical investigation of off-design performance of a transonic turbine

The off-design performance and a breakdown of the losses of a transonic turbine were determined by an analytical method that was previously developed for turbines of more conservative design. The analytically obtained performance map is compared with the performance map obtained from an experimental investigation of the turbine. The rotor hub conditions of incidence angle, relative Mach number, and reaction calculated from the analytical results are compared with those calculated from experimental data. The loss breakdown obtained for the transonic turbine did not differ substantially from that previously obtained from a turbine of more conservative design, except that large stator-exit shock losses were predicted for the transonic turbine at low speeds. The trends of the rotor hub incidence angle, relative Mach number, and reaction calculated from the analytical results agreed well with those calculated from the experimental data over the performance range. These trends indicate that, compared with a turbine of more conservative design, the transonic turbine operated over a much smaller range of incidence angle, a much wider range of rotor relative Mach number, and at a considerably lower level of reaction. Good over-all agreement was obtained between the analytically predicted performance and the experimental performance, except at 40-percent design speed, where in the analysis the stator reached limiting loading before the rotor choked. Since this discrepancy resulted from errors in the simplifying assumptions used in the analysis, it is regarded as a limitation in the analytical method as applied to a transonic turbine.

Whitney, Warren J↗

NASA NTRS1993

Space Shuttle Main Engine performance analysis

For a number of years, NASA has relied primarily upon periodically updated versions of Rocketdyne's power balance model (PBM) to provide space shuttle main engine (SSME) steady-state performance prediction. A recent computational study indicated that PBM predictions do not satisfy fundamental energy conservation principles. More recently, SSME test results provided by the Technology Test Bed (TTB) program have indicated significant discrepancies between PBM flow and temperature predictions and TTB observations. Results of these investigations have diminished confidence in the predictions provided by PBM, and motivated the development of new computational tools for supporting SSME performance analysis. A multivariate least squares regression algorithm was developed and implemented during this effort in order to efficiently characterize TTB data. This procedure, called the 'gains model,' was used to approximate the variation of SSME performance parameters such as flow rate, pressure, temperature, speed, and assorted hardware characteristics in terms of six assumed independent influences. These six influences were engine power level, mixture ratio, fuel inlet pressure and temperature, and oxidizer inlet pressure and temperature. A BFGS optimization algorithm provided the base procedure for determining regression coefficients for both linear and full quadratic approximations of parameter variation. Statistical information relative to data deviation from regression derived relations was also computed. A new strategy for integrating test data with theoretical performance prediction was also investigated. The current integration procedure employed by PBM treats test data as pristine and adjusts hardware characteristics in a heuristic manner to achieve engine balance. Within PBM, this integration procedure is called 'data reduction.' By contrast, the new data integration procedure, termed 'reconciliation,' uses mathematical optimization techniques, and requires both measurement and balance uncertainty estimates. The reconciler attempts to select operational parameters that minimize the difference between theoretical prediction and observation. Selected values are further constrained to fall within measurement uncertainty limits and to satisfy fundamental physical relations (mass conservation, energy conservation, pressure drop relations, etc.) within uncertainty estimates for all SSME subsystems. The parameter selection problem described above is a traditional nonlinear programming problem. The reconciler employs a mixed penalty method to determine optimum values of SSME operating parameters associated with this problem formulation.

Santi, L. Michael↗

NASA NTRS1995

Building a better leapfrog

In stellar dynamical computer simulations, as well as other types of simulations using particles, time step size is often held constant in order to guarantee a high degree of energy conservation. In many applications, allowing the time step size to change in time can offer a great saving in computational cost, but variable-size time steps usually imply a substantial degradation in energy conservation. We present a meta-algorithm' for choosing time steps in such a way as to guarantee time symmetry in any integration scheme, thus allowing vastly improved energy conservation for orbital calculations with variable time steps. We apply the algorithm to the familiar leapfrog scheme, and generalize to higher order integration schemes, showing how the stability properties of the fixed-step leapfrog scheme can be extended to higher order, variable-step integrators such as the Hermite method. We illustrate the remarkable properties of these time-symmetric integrators for the case of a highly eccentric elliptical Kepler orbit and discuss applications to more complex problems.

Hut, Piet↗

NASA NTRS1996

Multigrid, Fractional-Step Computation Of Flow

Speed of computer code solving three-dimensional Navier-Stokes equations of flow of incompressible fluid by fractional-step method increased significantly by use of multigrid procedures. In method, equations solved on general nonorthogonal curvilinear coordinate grid, using volume fluxes. At each time step, computations performed in two fractional steps. In first step, equations of conservation of momentum solved by use of gradient of pressure from previous time step via explicit approximate-factorization method, yielding approximate flow field that does not satisfy equation of conservation of mass. In second step, discrete Poisson-like equation with Neumann-type boundary conditions, formed by combining equations of conservation of momentum and mass, solved iteratively.

Kwak, Dochan↗

NASA NTRS1994

A Class of High-Resolution Explicit and Implicit Shock-Capturing Methods

The development of shock-capturing finite difference methods for hyperbolic conservation laws has been a rapidly growing area for the last decade. Many of the fundamental concepts, state-of-the-art developments and applications to fluid dynamics problems can only be found in meeting proceedings, scientific journals and internal reports. This paper attempts to give a unified and generalized formulation of a class of high-resolution, explicit and implicit shock capturing methods, and to illustrate their versatility in various steady and unsteady complex shock waves, perfect gases, equilibrium real gases and nonequilibrium flow computations. These numerical methods are formulated for the purpose of ease and efficient implementation into a practical computer code. The various constructions of high-resolution shock-capturing methods fall nicely into the present framework and a computer code can be implemented with the various methods as separate modules. Included is a systematic overview of the basic design principle of the various related numerical methods. Special emphasis will be on the construction of the basic nonlinear, spatially second and third-order schemes for nonlinear scalar hyperbolic conservation laws and the methods of extending these nonlinear scalar schemes to nonlinear systems via the approximate Riemann solvers and flux-vector splitting approaches. Generalization of these methods to efficiently include real gases and large systems of nonequilibrium flows will be discussed. Some perbolic conservation laws to problems containing stiff source terms and terms and shock waves are also included. The performance of some of these schemes is illustrated by numerical examples for one-, two- and three-dimensional gas-dynamics problems. The use of the Lax-Friedrichs numerical flux to obtain high-resolution shock-capturing schemes is generalized. This method can be extended to nonlinear systems of equations without the use of Riemann solvers or flux-vector splitting approaches and thus provides a large savings for multidimensional, equilibrium real gases and nonequilibrium flow computations.

Yee, H. C.↗

NASA NTRS1995

Bounce-Averaged Hamiltonian for Charged Particles in an Axisymmetric but Nondipolar Model Magnetosphere

In order to facilitate bounce-averaged guiding center simulations of geomagnetically trapped particles, we express the kinetic energy of a particle with magnetic coordinates (L,phi) as an analytic function of the first two adiabatic invariants (M, J) and the L value of the field line. The magnetic field model is axisymmetric, consisting of a dipolar B field plus a uniform southward magnetic field parallel to the dipole moment mu(sub E). This model magnetosphere is surrounded by a circular equatorial neutral line whose radius b is an adjustable parameter. The L value of a field line is (by definition) inversely proportional to the flux enclosed by the corresponding magnetic shell of equatorial radius r(sub 0), and the L value at the neutral line (r(sub 0) = b) is denoted L*. The azimuthal coordinate phi measures magnetic local time. The best functional representation found for the normalized difference (L(exp 3)a(exp 3)/mu(sub E))(B(sub m) - B(sub 0)) between mirror-point field B(sub m) and equatorial field B(sub 0) along any field line is a 5-term expansion in powers (2/3 through 6/3) of the quantity X equivalent to (La/mu(sub E))(exp 1/2)K, where K equivalent to (J(exp 2)/8m(sub 0)M)(exp 1/2) is an adiabatically conserved quantity independent of particle energy, m(sub 0) is the rest mass of the particle, and a is the radius of the Earth. This functional form is motivated by results for limiting cases in which particles mirror very near and very far from the magnetic equator. Expansion coefficients corresponding to various powers of X are obtained from least squares fits to numerically computed results for X as a function of L and B(sub m). These are accurately expressible as fourth-order polynomials in (r(sub 0)/b)(exp 3), hence indirectly as functions of L/L* = 3La/2b. This representation, which leads (except for a manageably small region of parameter space) to better than 1% accuracy in the specification of B(sub m) as a function of K and L, allows bounce-averaged guiding center simulations to be performed without actually tracing the bounce motions of individual particles. Bounce-averaged drifts L' (meridional) and phi' (azimuthal) are proportional to derivatives of the Hamiltonian H (sum of kinetic and potential energies) with respect to phi and L, respectively. Our formulation thus provides a computationally efficient method for tracing the bounce-averaged adiabatic motion (conserving all three invariants) and nonadiabatic transport (violating the third invariant while conserving the first two invariants) of geomagnetically trapped particles in the model magnetosphere.

Schulz, Michael↗

NASA NTRS1996

Towards and FVE-FAC Method for Determining Thermocapillary Effects on Weld Pool Shape

Several practical materials processes, e.g., welding, float-zone purification, and Czochralski crystal growth, involve a pool of molten metal with a free surface, with strong temperature gradients along the surface. In some cases, the resulting thermocapillary flow is vigorous enough to convect heat toward the edges of the pool, increasing the driving force in a sort of positive feedback. In this work we examine this mechanism and its effect on the solid-liquid interface through a model problem: a half space of pure substance with concentrated axisymmetric surface heating, where surface tension is strong enough to keep the liquid free surface flat. The numerical method proposed for this problem utilizes a finite volume element (FVE) discretization in cylindrical coordinates. Because of the axisymmetric nature of the model problem, the control volumes used are torroidal prisms, formed by taking a polygonal cross-section in the (r, z) plane and sweeping it completely around the z-axis. Conservation of energy (in the solid), and conservation of energy, momentum, and mass (in the liquid) are enforced globally by integrating these quantities and enforcing conservation over each control volume. Judicious application of the Divergence Theorem and Stokes' Theorem, combined with a Crank-Nicolson time-stepping scheme leads to an implicit algebraic system to be solved at each time step. It is known that near the boundary of the pool, that is, near the solid-liquid interface, the full conduction-convection solution will require extremely fine length scales to resolve the physical behavior of the system. Furthermore, this boundary moves as a function of time. Accordingly, we develop the foundation of an adaptive refinement scheme based on the principles of Fast Adaptive Composite Grid methods (FAC). Implementation of the method and numerical results will appear in a later report.

Canright, David↗

NASA NTRS1999

Accuracy of the Independent Pixel Approximation at Absorbing Wavelengths

In order to correctly interpret shortwave cloud radiation measured by satellites and ground-based radiometers, or by two aircraft flying above and below clouds, we need to better understand interactions between inhomogeneous clouds and solar radiation. The discrepancies between shortwave absorption inferred from measurements and predicted by models, between cloud optical depths estimated from satellites and ground measurements, between single scattering albedo retrieved from in situ radiation measurements and computed from measured droplet size distribution, among others, are strongly affected by cloud horizontal inhomogeneity. Net horizontal photon transport (i. e., horizontal fluxes) are a direct consequence of the inhomoqeneity in cloud structure. Horizontal fluxes and their effect on the accuracy of the pixel-by-pixel one-dimensional (1 D) radiative transfer calculations has recently undergone close scrutiny for conservative scattering. However, the properties and magnitude of horizontal fluxes in absorbing wavelengths are still poorly understood. As far as we are aware, only Ackerman and Cox and Titov discussed correlations between horizontal fluxes at absorbing wavelengths, though these were far from comprehensive. This paper partly fills this gap. We discuss here of whether the accuracy of the Independent Pixel Approximation (IPA), a 1 D radiative transfer approximation for each pixel, is a better model for multiple scattering at conservative or at absorbing wavelengths. Issues addressed here are: (1) dependence of net horizontal fluxes on single scattering albedo; (2) connection between pixel-by-pixel accuracy of the IPA and horizontal fluxes and (3) radiative smoothing and horizontal fluxes at absorbing wavelengths. In contrast to the traditional understanding of IPA, we study IPA accuracies not only for reflectance but also for transmittance and absorptance at both conservative and absorbing wavelengths. In spite of the apparent similarity between the three processes, dependence of IPA accuracies on single-scattering albedo is completely different. As a result, cloud optical properties retrieved from high resolution satellite images and ground-based measurements using IPA at absorbing channels will have different accuracies.

Marchak, Alexander↗

NASA NTRS1999

Engineered Solutions to Reduce Occupational Noise Exposure at the NASA Glenn Research Center: A Five-Year Progress Summary (1994-1999)

At the NASA John H. Glenn Research Center at Lewis Field (formerly the Lewis Research Center), experimental research in aircraft and space propulsion systems is conducted in more than 100 test cells and laboratories. These facilities are supported by a central process air system that supplies high-volume, high-pressure compressed air and vacuum at various conditions that simulate altitude flight. Nearly 100,000 square feet of metalworking and specialized fabrication shops located on-site produce prototypes, models, and test hardware in support of experimental research operations. These activities, comprising numerous individual noise sources and operational scenarios, result in a varied and complex noise exposure environment, which is the responsibility of the Glenn Research Center Noise Exposure Management Program. Hearing conservation, community noise complaint response and noise control engineering services are included under the umbrella of this Program, which encompasses the Occupational Safety and Health Administration (OSHA) standard on occupational noise exposure, Sec. 29CFR 1910.95, as well as the more stringent NASA Health Standard on Hearing Conservation. Prior to 1994, in the absence of feasible engineering controls, strong emphasis had been placed on personal hearing protection as the primary mechanism for assuring compliance with Sec. 29CFR 1910.95 as well as NASA's more conservative policy, which prohibits unprotected exposure to noise levels above 85 dB(A). Center policy and prudent engineering practice required, however, that these efforts be extended to engineered noise controls in order to bring existing work areas into compliance with Sec. 29CFR 1910.95 and NASA's own policies and to ensure compliance for new installations. Coincident with the establishment in 1995 of a NASA wide multi-year commitment of funding for environmental abatement projects, the Noise Exposure Management Program was established, with its focus on engineering approaches to reducing occupational and community noise exposure. The organization and mission of this Program were documented shortly after its inception, and individual programmatic components have been the subject of subsequent papers. This paper summarizes the status and accomplishments of the engineering aspects of the Program from a five-year retrospective viewpoint and includes a review of retrofit noise control solution strategies, which have not been previously documented.

Cooper, Beth A.↗