SEARCH · Engineering Papers
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.
Solutions of the Hugoniot shock equations and Meyer expansion equations are plotted in such a manner as to permit the pressure distribution, the local Mach number, and the angles of shock waves on arbitrary sharp-nose airfoils at supersonic speeds to be obtained directly. (author)
Second-year efforts within a three-year study to develop and extend finite element (FE) methodology to efficiently handle the transient/steady state response of rotor-bearing-stator structure associated with gas turbine engines are outlined. The two main areas aim at (1) implanting the squeeze film damper element into a general purpose FE code for testing and evaluation; and (2) determining the numerical characteristics of the FE-generated rotor-bearing-stator simulation scheme. The governing FE field equations are set out and the solution methodology is presented. The choice of ADINA as the general-purpose FE code is explained, and the numerical operational characteristics of the direct integration approach of FE-generated rotor-bearing-stator simulations is determined, including benchmarking, comparison of explicit vs. implicit methodologies of direct integration, and demonstration problems.
The burning of a candle, as typical non-propagating diffusion flame, has been used by a number of researchers to study the effects of electric fields on flame, spontaneous flame oscillation and flickering phenomena, and flame extinction. In normal gravity, the heat released from combustion creates buoyant convection that draws oxygen into the flame. The strength of the buoyant flow depends on the gravitational level and it is expected that the flame shape, size and candle burning rate will vary with gravity. Experimentally, there exist studies of candle burning in enhanced gravity (i.e. higher than normal earth gravity, g(sub e)), and in microgravity in drop towers and space-based facilities. There are, however, no reported experimental data on candle burning in partial gravity (g < g(sub e)). In a previous numerical model of the candle flame, buoyant forces were neglected. The treatment of momentum equation was simplified using a potential flow approximation. Although the predicted flame characteristics agreed well with the experimental results, the model cannot be extended to cases with buoyant flows. In addition, because of the use of potential flow, no-slip boundary condition is not satisfied on the wick surface. So there is some uncertainty on the accuracy of the predicted flow field. In the present modeling effort, the full Navier-Stokes momentum equations with body force term is included. This enables us to study the effect of gravity on candle flames (with zero gravity as the limiting case). In addition, we consider radiation effects in more detail by solving the radiation transfer equation. In the previous study, flame radiation is treated as a simple loss term in the energy equation. Emphasis of the present model is on the gas-phase processes. Therefore, the detailed heat and mass transfer phenomena inside the porous wick are not treated. Instead, it is assumed that a thin layer of liquid fuel coated the entire wick surface during the burning process. This is the limiting case that the mass transfer process in the wick is much faster than the evaporation process at the wick surface.
During the late 1970's and early 1980's, concerns about the potential interference of wind turbine generators with electromagnetic communication signals led to a series of research studies, both in the laboratory and in the field, conducted by the staff of the University of Michigan Radiation Laboratory. These studies were sponsored by organizations such as the U.S. Department of Energy, the Solar Energy Research Institute, and private developers of wind power stations. Research objectives were to identify the mechanisms by which wind turbines might adversely affect communication signals, estimate the severity of these effects for different types of signals (e.g. television, radio, microwave, and navigation), and formulate mathematical models with which to predict the sizes of potential interference zones around wind turbines and wind power plants. This work formed the basis for preliminary standards on assessing electromagnetic interference (EMI) by wind turbines. With the current renewal of interest in wind energy projects, it is appropriate that the many experimental and analytical aspects of this pioneering work be reviewed and correlated. The purpose of this study is to combine test data and theory from previously published and unpublished research reports into a unified and consistent set of equations which are useful for estimating potential levels of television interference from wind turbines. To be comprehensive, these equations will include both horizontal-axis and vertical-axis wind turbines (HAWT's and VAWT's), blade configuration parameters (e.g. number, size, material, twist, and coning), signal frequency and power, and directional characteristics of the receiving antenna. The approach that is followed in this report is as follows. First, some basic equations that describe electromagnetic signals with interference are presented without detailed derivations, since the latter are available in the references. Minor changes in terminology are made for purposes of consistency. Next, the concept of a signal scatter ratio is introduced, which defines the fraction of the signal impinging on a wind turbine that is scattered by its blades onto a nearby receiver. Equations from references are modified for the calculation of experimental scatter ratios (from measured signals containing interference) and idealized scatter ratios (from rotor characteristics and relative locations of the transmitter, the turbine, and the receiver). Experimental and idealized scatter ratios are then calculated and compared for 75 cases from the literature, in which TVI measurements were made around a variety of wind turbines. An empirical equation is then defined for estimating the probability that an actual scatter ratio will differ from an idealized ratio by a given amount. Finally a sample calculation of the size of a potential TV interference zone around a hypothetical wind power station is presented.
An implicit finite-difference scheme is developed for the efficient numerical solution of nonlinear hyperbolic systems in conservation-law form. The algorithm is second-order time-accurate, noniterative, and in a spatially factored form. Second- or fourth-order central and second-order one-sided spatial differencing are accommodated within the solution of a block tridiagonal system of equations. Significant conceptual and computational simplifications are made for systems whose flux vectors are homogeneous functions (of degree one), e.g., the Eulerian gasdynamic equations. Conservative hybrid schemes, which switch from central to one-sided spatial differencing whenever the local characteristic speeds are of the same sign, are constructed to improve the resolution of weak solutions. Numerical solutions are presented for a nonlinear scalar model equation and the two-dimensional Eulerian gasdynamic equations.
This paper describes the performance of a new family of NASA-developed equation solvers used for large-scale (i.e. 551,705 equations) structural analysis. To minimize computer time and memory, the solvers are divided by application and matrix characteristics (sparse/dense, real/complex, symmetric/nonsymmetric, size: in-core/out of core) and exploit the hardware features of current and future computers. In this paper, the equation solvers, which are written in FORTRAN, and are therefore easily transportable, are shown to be faster than specialized computer library routines utilizing assembly code. Twenty NASA structural benchmark models with NASA solver timings reside on World Wide Web with a challenge to beat them.
Thrust chamber performance is evaluated in terms of an analytical model incorporating all the loss processes that occur in a real rocket motor. The important loss processes in the real thrust chamber were identified, and a methodology and recommended procedure for predicting real thrust chamber vacuum specific impulse were developed. Simplified equations for the calculation of vacuum specific impulse are developed to relate the delivered performance (both vacuum specific impulse and characteristic velocity) to the ideal performance as degraded by the losses corresponding to a specified list of loss processes. These simplified equations enable the various performance loss components, and the corresponding efficiencies, to be quantified separately (except that interaction effects are arbitrarily assigned in the process). The loss and efficiency expressions presented can be used to evaluate experimentally measured thrust chamber performance, to direct development effort into the areas most likely to yield improvements in performance, and as a basis to predict performance of related thrust chamber configurations.
An efficient implicit finite-difference algorithm for the gas dynamic equations utilizing matrix reduction techniques is presented. A significant reduction in arithmetic operations is achieved while maintaining the same favorable stability characteristics and generality found in the Beam and Warming approximate factorization algorithm. Steady-state solutions to the conservative Euler equations in generalized coordinates are obtained for transonic flows about a NACA 0012 airfoil. The theoretical extension of the matrix reduction technique to the full Navier-Stokes equations in Cartesian coordinates is presented in detail. Linear stability, using a Fourier stability analysis, is demonstrated and discussed for the one-dimensional Euler equations. It is shown that the method offers advantages over the conventional Beam and Warming scheme and can retrofit existing Beam and Warming codes with minimal effort.
In the first part, the performance of the Absorption Heat Pump (AHP) with water-sulfuric acid and water-magnesium chloride as two new refrigerant-absorbent fluid pairs was investigated. A model was proposed for the analysis of the new working pairs in a heat pump system, subject to different temperature lifts. Computer codes were developed to calculate the Coefficient of Performance (COP) of the system with the thermodynamic properties of the working fluids obtained from the literature. The study shows the potential of water-sulfuric acid as a satisfactory replacement for water-lithium bromide in the targeted temperature range. The performance of the AHP using water-magnesium chloride as refrigerant-absorbent pair does not compare well with those obtained using water-lithium bromide. The second part concentrated on the design and testing of a simple ElectroHydrodynamic (EHD) Pump. A theoretical design model based on continuum electromechanics was analyzed to predict the performance characteristics of the EHD pump to circulate the fluid in the absorption heat pump. A numerical method of solving the governing equations was established to predict the velocity profile, pressure - flow rate relationship and efficiency of the pump. The predicted operational characteristics of the EHD pump is comparable to that of turbomachinery hardware; however, the overall efficiency of the electromagnetic pump is much lower. An experimental investigation to verify the numerical results was conducted. The pressure - flow rate performance characteristics and overall efficiency of the pump obtained experimentally agree well with the theoretical model.
Closed-form equations for the lift, drag, and pitching moment coefficients of two dimensional airfoil sections in steady subsonic flow were obtained from published theoretical and experimental results. A turbulent boundary layer was assumed to exist on the airfoil surfaces. The effects of section angle of attack, Mach number, Reynolds number, and the specific airfoil type were considered. The equations were applicable through an angle of attack range of -180 deg to +180 deg; however, above about + or - 20 deg, the section characteristics were assumed to be functions only of angle of attack. A computer program is presented which evaluates the equations for a range of Mach numbers and angles of attack. Calculated results for the NACA 23012 airfoil section were compared with experimental data.
A mixed finite-difference scheme is presented for the stress and free vibration analysis of simply supported nonhomogeneous and layered orthotropic thick plates. The analytical formulation is based on the linear, three-dimensional theory of orthotropic elasticity and a Fourier approach is used to reduce the governing equations to six first-order ordinary differential equations in the thickness coordinate. The governing equations possess a symmetric coefficient matrix and are free of derivatives of the elastic characteristics of the plate. In the finite difference discretization two interlacing grids are used for the different fundamental unknowns in such a way as to reduce both the local discretization error and the bandwidth of the resulting finite-difference field equations. Numerical studies are presented for the effects of reducing the interior and boundary discretization errors and of mesh refinement on the accuracy and convergence of solutions. It is shown that the proposed scheme, in addition to a number of other advantages, leads to highly accurate results, even when a small number of finite difference intervals is used.
In this paper a distributed parameter model for the estimation of modal characteristics of NASA Mini-Mast truss is proposed. A closed-form solution of the Timoshenko beam equation, for a uniform cantilevered beam with two concentrated masses, is derived so that the procedure and the computational effort for the estimation of modal characteristics are improved. A maximum likelihood estimator for the Timoshenko beam model is also developed. The resulting estimates from test data by using Timoshenko beam model are found to be comparable to those derived from other approaches.
A review of the acoustics of moving media is provided, taking into account a derivation of basic equations, wave-like solutions of acoustic equations, the solutions of acoustic equations by superposition of elementary sources, the source distribution in unbounded regions, the radiation field, energy relations, and moving sound sources. The characteristics of aerodynamic sound are considered along with the effect of solid boundaries, the effect of uniform flow, acoustic field equations, and the effects of nonuniform mean flow on the generation of sound. Attention is given to Lighthill's acoustic analogy, the solution of Lighthill's equation when no solid boundaries are present, the application of Lighthill's theory to turbulent flows, and the physics of jet noise.
Laminar convective heat transfer and boundary layer equations for cones at escape velocity
Characteristics such as storage capacity of the soil, volume of storage in vegetative matter, and volume of storage available in local depressions are expressed in empirical watershed runoff equations as one or more coefficients. Conventional techniques for estimating coefficients representing the spatial distribution of these characteristics over a watershed drainage area are subjective and produce significant errors. Characteristics of the wear surface are described as a single coefficient called the curve number.
A mathematical model taking into account small (and constant) gravitational levels is developed for vaporization of an isolated liquid droplet suspended in a stagnant atmosphere. A goal of the present analysis is to see how small gravitational levels affect droplet gasification characteristics. Attention is focused upon determining the effects on gas-phase phenomena. The conservation equations arc normalized and nondimensionalized, and a small parameter that accounts for the effects of gravity is identified. This parameter is the square of the inverse of a Froude number based on the gravitational acceleration, the droplet radius, and a characteristic gas-phase velocity at the droplet surface. Asymptotic analyses are developed in terms of this parameter. In the analyses, different spatial regions are identified. Near a droplet, gravitational effects are negligible in the first approximation, and the flowfield is spherically symmetric to the leading order. Analysis shows, however, that outer zones exist where gravitational effects cannot be neglected; it is expected that a stagnation point will be present in an outer zone that is not present when gravity is totally absent. The leading order and higher-order differential equations for each zone are derived and solved. The solutions allow the effects of gravity on vaporization rates and temperature, velocity and species fields to be determined.
Solution of the parabolic Navier-Stokes (PNS) equations for supersonic flows is discussed, and compatibility of the PNS method with the triple-deck theory of Stewartson (1974) is demonstrated. Characteristic and stability analyses show that use of an appropriate filter on the pressure term in the x-momentum equation can suppress the departure solutions, giving the usual desired weak interaction solution. An Alternating Direction Explicit procedure, with minimal computer storage requirements compared to the full Navier-Stokes solvers, is proposed to calculate strongly interacting flows using a global iteration procedure for the PNS equations. The PNS equations are used to solve the hypersonic viscous interaction problem, and good agreement is found with experimental results.
The characteristic impedance and other derivative acoustic properties of a material can be derived from impedance tube data using the specific impedance measured from samples with two different thicknesses. In practice, samples are chosen so that their respective thicknesses differ by a factor of 2. This simplifies the solution of the equations relating the properties of the two samples so that the computation of the characteristic impedance is straightforward. This approach has at least two drawbacks. One is that it is often difficult to acquire or produce samples with precisely a factor of 2 difference in thickness. A second drawback is that the phase information contained in the imaginary part of the propagation constant must be unwrapped before subsequent computations are performed. For well-behaved samples, this is not a problem. For ill behaved samples of unknown properties, the phase unwrapping process can be tedious and difficult to automate. Two alternative approaches have been evaluated which remove the factor-of-2 sample thickness requirement and directly compute unwrapped phase angles. One uses a Newton-Raphson approach to solve for the roots of the samples' simultaneous equations. The other produces a wave number space diagram in which the roots are clearly discernable and easily extracted. Results are presented which illustrate the flexibility of analysis provided by the new approaches and how this can be used to better understand the limitations of the impedance tube data.