Engineering Papers⌕ Search

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.

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

At least 973 records · Page 54

Electromagnetic diffraction efficiencies for plane reflection diffraction gratings

The theory and computer programs, based on electromagnetic theory, for the analysis and design of echelle gratings were developed. The gratings are designed for instruments that operate in the ultraviolet portion of the spectrum. The theory was developed so that the resulting computer programs will be able to analyze deep (up to 30 wavelengths) gratings by including as many as 100 real or homogeneous diffraction orders. The program calculates the complex amplitude coefficient for each of the diffracted orders. A check on the numerical method used to solve the integral equations is provided by a conservation of energy calculation.

Marathay, A. S.↗

Propagation of solar disturbances in interplanetary space

Time-dependent solutions of a one-fluid model of the interplanetary medium are investigated. This set of unsteady hydrodynamic equations has been written in conservation form in order to apply the Lax-Wendroff method for the solution of this problem. The initial condition is specified by a pulse at 1 solar radius. The equilibrium condition is chosen to be the steady solution of a quiet solar wind. The specified solar disturbances in this calculation are allowed to be both sub- and supersonic by the present theoretical formulation. The results are presented in terms of density, velocity, and temperature profiles of the interplanetary gas flow at heliocentric distances up to about 10 AU at any particular time. The trajectories of disturbances for various initial pulses are shown. Some 1972 solar-flare observational data are compared with these theoretical calculations. From these calculations, the effects on the interplanetary environment, due to the propagation of solar disturbances, can be determined.

Wu, S. T.↗

Nonlinear interaction between two electrostatic harmonics in a plasma

The problem of harmonic and subharmonic generation of electrostatic waves in a general collisionless plasma is treated using coupled-mode theory based on two time scales. A novel feature is that one of the two interacting waves may be a negative-energy wave. Since the model describing the medium need not be specified, only a general linear and nonlinear conductivity or an equivalent description is required. Just by invoking wave-energy conservation, the coupled-mode equations are obtained in such a way that unequivocal conclusions can be drawn. When both waves have positive energy, they exchange part of it in a periodic fashion, provided that both have some energy initially. If initially all the energy is in the fundamental, all of it will eventually end up (irreversibly) in the second harmonic. If all the energy is in the harmonic initially, no generation of the fundamental (or subharmonic) will take place. If one of the two waves is a negative-energy wave, an explosive instability develops, regardless of initial values. For comparable conditions, the instability time depends on whether the negative-energy wave is in the fundamental or in the upper harmonic.

Verheest, F.↗

Stability bounds for the control of large space structures

Balas (1977) has discussed the stability problem of reduced-order regulators and estimators in terms of control and observation 'spillover'. The term 'control spillover' was used to define that part of the feedback control which excites the uncontrolled (or residual) modes, and 'observation spillover' was used to define that part of the measurement which is contaminated by residual modes. In this paper, two sufficient conditions are derived via Lyapunov methods for asymptotic stability of large space structures using a class of reduced-order controllers. These conditions give allowable bounds on the spectral norms of control and observation spillover terms. The sufficient condition given by a specified inequality equation appears to be less conservative, and should be useful as a design tool for the control of large space structures.

Joshi, S. M.↗

Radio emission from nova shells

This paper presents data from multifrequency radio observations of the novae HR Del 1967 and FH Ser 1970 made from June 1970 to December 1978 and of V1500 Cyg 1975 made from September 1975 to December 1978. The data for FH Ser and V1500 Cyg constitute a complete history of the important phases of their radio 'light' curves, while all but the rising portion of the optically thick phase of an expanding nova shell is present for HR Del. A simple solution of the spherically symmetric equations of mass and momentum conservation for an isothermal gas, for which pressure gradients are effectively decoupled from the dynamics, is shown to fit the radio data for all three novae quite well. The main features of this model are strong density and velocity gradients that lead to long periods when the 'observable' shell evolves from an optically thick 'disk' to a completely optically thin shell.

Hjellming, R. M.↗

Electrostatic model of a quiet auroral arc

The physics and equations used in the proposed electrostatic model of a quiet auroral arc are discussed. The procedure for obtaining approximate solutions of Poisson's equation coupled to ionosphere current conservation is given. A comparison with a specific discrete auroral arc as seen from the auroral satellite S3-3 is presented.

Chiu, Y. T.↗

Numerical Modeling of Three-Dimensional Confined Flows

A three dimensional confined flow model is presented. The flow field is computed by calculating velocity and enthalpy along a set of streamlines. The finite difference equations are obtained by applying conservation principles to streamtubes constructed around the chosen streamlines. With appropriate substitutions for the body force terms, the approach computes three dimensional magnetohydrodynamic channel flows. A listing of a computer code, based on this approach is presented in FORTRAN IV language. The code computes three dimensional compressible viscous flow through a rectangular duct, with the duct cross section specified along the axis.

Greywall, M. S.↗

Design of supercritical swept wings

Computational fluid dynamics are used to discuss problems inherent to transonic three-dimensional flow past supercritical swept wings. The formulation for a boundary value problem for the flow past the wing is provided, including consideration of weak shock waves and the use of parabolic coordinates. A swept wing code is developed which requires a mesh of 152 x 10 x 12 points and 200 time cycles. A formula for wave drag is calculated, based on the idea that the conservation form of the momentum equation becomes an entropy inequality measuring the drag, expressible in terms of a small-disturbance equation for a potential function in two dimensions. The entropy inequality has been incorporated in a two-dimensional code for the analysis of transonic flow over airfoils. A method of artificial viscosity is explored for optimum pressure distributions with design, and involves a free boundary problem considering speed over only a portion of the wing.

Garabedian, P.↗

A 3-dimensional Navier-Stokes-Euler code for blunt-body flow computations

The shock-layer flowfield is obtained with or without viscous and heat-conducting dissipations from the conservative laws of fluid dynamics equations using a shock-fitting implicity finite-difference technique. The governing equations are cast in curvilinear-orthogonal coordinates and transformed to the domain between the shock and the body. Another set of equations is used for the singular coordinate axis, which, together with a cone generator away from the stagnation point, encloses the computation domain. A time-dependent alternating direction implicit factorization technique is applied to integrate the equations with local-time increment until a steady solution is reached. The shock location is updated after the flowfield computation, but the wall conditions are implemented into the implicit procedure. Innovative procedures are introduced to define the initial flowfield, to treat both perfect and equilibrium gases, to advance the solution on a coarse-to-fine grid sequence, and to start viscous flow computations from their corresponding inviscid solutions. The results are obtained from a grid no greater than 28 by 18 by 7 and converged within 300 integration steps. They are of sufficient accuracy to start parabolized Navier-Stokes or Euler calculations beyond the nose region, to compare with flight and wind-tunnel data, and to evaluate conceptual designs of reentry spacecraft.

Li, C. P.↗

Two-dimensional internal compressible viscous flows using semi-elliptic analysis

Subsonic viscous flows through planar cascades are examined by means of a semielliptic form of the Navier-Stokes equations, which is obtained by neglecting the streamwise diffusion of momentum and energy in their respective equations and providing upstream interaction through the pressure field. The influence of downstream conditions is propagated upstream in the subsonic flow, via a streamwise implicit determination of the pressure field, in order to permit departure-free solutions for the overall flow problem. The continuity equation is thereby satisfied directly. Governing equations are written in strong conservation law form, in terms of general, nonorthogonal transformed coordinates and Cartesian velocity components.

Ghia, K. N.↗

Application of Runge Kutta time marching scheme for the computation of transonic flows in turbomachines

Numerical solutions of the unsteady Euler equations are obtained using the classical fourth order Runge Kutta time marching scheme. This method is fully explicit and is applied to the governing equations in the finite volume, conservation law form. In order to determine the efficiency of this scheme for solving turbomachinery flows, steady blade-to-blade solutions are obtained for compressor and turbine cascades under subsonic and transonic flow conditions. Computed results are compared with other numerical methods and wind tunnel measurements. The present study also focuses on other important numerical aspects influencing the performance of the algorithm and the solution accuracy such as grid types, boundary conditions, and artificial viscosity. For this purpose, H, O, and C type computational grids as well as characteristic and extrapolation type boundary conditions are included in the solution procedure.

Subramanian, S. V.↗

Interface procedures for overlapping grids

Interpolation at grid boundaries is studied for the purpose of solving partial differential equations using either implicit or conservative explicit finite-difference methods on multi-component overlapping grid systems.

Mastin, C. W.↗

Heat and mass transfer in a dissociated laminar boundary layer of air with consideration of the finite rate of chemical reaction

A numerical investigation of heat and mass transfer in a dissociated laminar boundary layer of air on an isothermal flat plate is carried out for different degrees of cooling of the wall. A finite-difference chemical model is used to study elementary reactions involving NO2 and N2O. The analysis is based on equations of continuity, momentum, energy, conservation and state for the two-dimensional viscous flow of a reacting multicomponent mixtures. Attention is given to the effects of both catalyticity and noncatalyticity of the wall.

Oyegbesan, A. O.↗

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

Landing-Time-Controlled Management Of Air Traffic

Conceptual system controls aircraft with old and new guidance equipment. Report begins with overview of concept, then reviews controller-interactive simulations. Describes fuel-conservative-trajectory algorithm, based on equations of motion for controlling landing time. Finally, presents results of piloted simulations.

Erzberger, Heinz↗

Application of Runge Kutta time marching scheme for the computation of transonic flows in turbomachines

Numerical solutions of the unsteady Euler equations are obtained using the classical fourth order Runge Kutta time marching scheme. This method is fully explicit and is applied to the governing equations in the finite volume, conservation law form. In order to determine the efficiency of this scheme for solving turbomachinery flows, steady blade-to-blade solutions are obtained for compressor and turbine cascades under subsonic and transonic flow conditions. Computed results are compared with other numerical methods and wind tunnel measurements. The study also focuses on other important numerical aspects influencing the performance of the algorithm and the solution accuracy such as grid types, boundary conditions and artificial viscosity. For this purpose, H, O, and C type computational grids as well as characteristic and extrapolation type boundary conditions are included in solution procedures.

Subramanian, S. V.↗

A mathematical model of a lithium/thionyl chloride primary cell

A 1-D mathematical model for the lithium/thionyl chloride primary cell was developed to investigate methods of improving its performance and safety. The model includes many of the components of a typical lithium/thionyl chloride cell such as the porous lithium chloride film which forms on the lithium anode surface. The governing equations are formulated from fundamental conservation laws using porous electrode theory and concentrated solution theory. The model is used to predict 1-D, time dependent profiles of concentration, porosity, current, and potential as well as cell temperature and voltage. When a certain discharge rate is required, the model can be used to determine the design criteria and operating variables which yield high cell capacities. Model predictions can be used to establish operational and design limits within which the thermal runaway problem, inherent in these cells, can be avoided.

Evans, T. I.↗

Flow transition with 2-D roughness elements in a 3-D channel

We develop a new numerical approach to study the spatially evolving instability of the streamwise dominant flow in the presence of roughness elements. The difficulty in handling the flow over the boundary surface with general geometry is removed by using a new conservative form of the governing equations and an analytical mapping. The numerical scheme uses second-order backward Euler in time, fourth-order central differences in all three spatial directions, and boundary-fitted staggered grids. A three-dimensional channel with multiple two-dimensional-type roughness elements is employed as the test case. Fourier analysis is used to decompose different Fourier modes of the disturbance. The results show that surface roughness leads to transition at lower Reynolds number than for smooth channels.

Liu, Zhining↗