Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “numerical techniques”

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 433 records · Page 24

A general perturbation approach for the equations of fluid dynamics

An efficient numerical technique to produce accurate solutions to the equations of fluid dynamics is presented where the governing equations are perturbed about an approximate solution and solved by finite-difference methods on a coarsened grid. The result is a scheme which substantially reduces the number of grid points necessary to accurately resolve the flow. Applications are presented for the two-dimensional Euler equations perturbed about a solution of the transonic full potential equation. However, the concept is applicable to arbitrary equation sets, higher dimensions and for a wide variety of applications.

Chow, L. J.↗

Statistical mechanics of light elements at high pressure. V Three-dimensional Thomas-Fermi-Dirac theory

A numerical technique for solving the Thomas-Fermi-Dirac (TED) equation in three dimensions, for an array of ions obeying periodic boundary conditions, is presented. The technique is then used to calculate deviations from ideal mixing for an alloy of hydrogen and helium at zero temperature and high presures. Results are compared with alternative models which apply perturbation theory to calculation of the electron distribution, based upon the assumption of weak response of the electron gas to the ions. The TFD theory, which permits strong electron response, always predicts smaller deviations from ideal mixing than would be predicted by perturbation theory. The results indicate that predicted phase separation curves for hydrogen-helium alloys under conditions prevailing in the metallic zones of Jupiter and Saturn are very model dependent.

Macfarlane, J. J.↗

Numerical simulations of counterstreaming plasmas and their relevance to interhemispheric flows

The collisionless expansion of counterstreaming plasmas has been studied in order to elucidate the basic physical processes that may be operating during the initial refilling of depleted flux tubes after a magnetic storm. The numerical technique applied is briefly described, and simulation results are presented. The simulation geometry consisted of two high-density H(+)-O(+) electron plasmas separated by a low-density H(+)-electron plasma. The temporal evolution of the expanding plasmas and the electostatic potential in the region between the two sources is described. The main interacting streams are found to be stable with respect to both the ion acoustic and ion cyclotron modes and only the suprathermal forerunner ions are unstable with respect to the ion cyclotron mode. The results also suggest that a localized potential hill can form at the equator and that this potential hill can play an important role in the subsequent trapping and thermalization of the ion streams.

Singh, N.↗

Linear and Nonlinear Aspects of Rotordynamics

Excessive vibrations of the liquid oxygen pump in the Space Shuttle's Main Engine have been recorded during hot firing ground testing. In order to determine mathematical explanations of this possibility, destructive phenomenon differential equations have been examined which describe the rotordynamics of the pump. Modeling the rotor as a random eigenvalue problem was considered. Analytical expressions were derived for the solution in the case of symmetric damping and stiffness. This enables one to determine accuracy estimates when testing numerical techniques to solve both asymmetric and nonlinear problems. Finally, the rotor model has had nonlinear elements incorporated to improve its simulation of the pump and to expand the corresponding mathematical theory.

Day, W. B.↗

Ocean modelling on the CYBER 205 at GFDL

At the Geophysical Fluid Dynamics Laboratory, research is carried out for the purpose of understanding various aspects of climate, such as its variability, predictability, stability and sensitivity. The atmosphere and oceans are modelled mathematically and their phenomenology studied by computer simulation methods. The present state-of-the-art in the computer simulation of large scale oceans on the CYBER 205 is discussed. While atmospheric modelling differs in some aspects, the basic approach used is similar. The equations of the ocean model are presented along with a short description of the numerical techniques used to find their solution. Computational considerations and a typical solution are presented in section 4.

Cox, M.↗

Combustion Fundamentals Research

The various physical processes that occur in the gas turbine combustor and the development of analytical models that accurately describe these processes are discussed. Aspects covered include fuel sprays; fluid mixing; combustion dynamics; radiation and chemistry and numeric techniques which can be applied to highly turbulent, recirculating, reacting flow fields.

Source record↗

Design considerations for a monolithic, GaAs, dual-mode, QPSK/QASK, high-throughput rate transceiver

A monolithic, GaAs, dual mode, quadrature amplitude shift keying and quadrature phase shift keying transceiver with one and two billion bits per second data rate is being considered to achieve a low power, small and ultra high speed communication system for satellite as well as terrestrial purposes. Recent GaAs integrated circuit achievements are surveyed and their constituent device types are evaluated. Design considerations, on an elemental level, of the entire modem are further included for monolithic realization with practical fabrication techniques. Numerous device types, with practical monolithic compatability, are used in the design of functional blocks with sufficient performances for realization of the transceiver.

Kot, R. A.↗

Improvements in the accuracy and stability of algorithms for the small-disturbance and full-potential equations applied to transonic flows

Numerical techniques that improve the accuracy and stability of algorithms for the small disturbance and full potential equations used to calculate transonic flows are described. For the small disturbance equation, the algorithm improvements are: (1) the use of monotone switches in the type dependent finite differencing, and (2) the use of stable and simple second order accurate spatial differencing; these improvements are for steady and unsteady transonic flows. For the steady full potential equation, the improvement is in the use of a monotone switch in the type dependent finite differencing of an approximate factorization (AF2) algorithm. All these improvements are implemented in present computer codes by making minor coding modifications.

Goorjian, P. M.↗

Solar modulation of galactic cosmic rays. V - Time-dependent modulation

A numerical technique is introduced for the solution of the time-dependent equation for the solar modulation of galactic cosmic rays, in which transport in heliographic latitude is ignored. We illustrate this method with a model for the solar cycle variation in cosmic ray intensity. The variation is assumed to result from changes in the number of enhanced cosmic ray scattering regions, which are produced by flare-generated shock waves in the solar wind. The model can account for the observed solar cycle variation in intensity at different energies as well as the observed spatial gradients. This model also provides a natural explanation for the cosmic ray hysteresis effect.

Perko, J. S.↗

Numerical simulation of turbulent flows

Computational models of turbulence in incompressible Newtonian fluids governed by the Navier-Stokes equations are reviewed. The governing equations are presented, and both direct and large-eddy-simulations are examined. Resolution requirements and numerical techniques of spatial representation, definition of initial and boundary conditions, and time advancement are considered. Results of simulations of homogeneous turbulence in uniform shear, the evolution of a turbulent mixing layer, and turbulent channel flow are presented graphically and discussed.

Rogallo, R. S.↗

A new numerical method for the simulation of three-dimensional flow in a pipe

A new numerical technique for simulating three dimensional, unsteady, incompressible pipe flows is presented and its utility and accuracy is shown. Each vector function in the expansion of the velocity field is divergence free and satisfies the boundary conditions for viscous flow. Some of the benefits of the expansion technique are that pressure is eliminated from the dynamics, only two unknowns per mesh point are required, implicit treatment of the viscous terms is provided at no extra computational cost, and no fractional time steps are required. The method uses spectral expansions: Fourier series in the azimuthal and streamwise directions, and Jacobi polynominals in the radial direction. Previously announced in STAR as N82-31644

Leonard, A.↗

Spin dynamics and horizon sensor performance for the Solar Mesosphere Explorer

Numerical techniques for obtaining the trigger altitude of the Solar Mesosphere Explorer (SME) are detailed. The SME has a pair of fixed horizon sensors for attitude determination and initiation of the data acquisition procedures for measuring ozone concentrations at the limb. The spacecraft spin rate is calculated as each of the horizon sensors marks the limb with a pulse telemetered to the ground. A best fit spin dynamics model that takes into account the history of the satellite spin angle and all disturbances to the spin is described. Analyses show that the trigger altitude is calculated to an accuracy of 1.5-2.0 km. A significant indicator of the limb altitude is the stratospheric temperature at the 5 mb pressure level.

Lawrence, G. M.↗

Interactive phenomena in supersonic jet mixing problems. II - Numerical studies

Applications of the numerical techniques discussed in Part I of this article are presented. Fundamental wave/shear-layer interaction calculations are performed exhibiting the waves produced by high-speed turbulent mixing processes and the interactions of an expansion fan and shock wave with a developing shear layer. Calculations for an underexpanded supersonic jet into a supersonic external stream at different velocity ratios are presented which delineate betweeen weakly interactive and strongly interactive processes in the jet near field. Calculations for the multiple-cell shock structure of an underexpanded jet into still air exhibit the damping of wave intensities by the turbulence and compare quite favorably with available data.

Dash, S. M.↗

SUBOPT: A CAD program for suboptimal linear regulators

An interactive software package which provides design solutions for both standard linear quadratic regulator (LQR) and suboptimal linear regulator problems is described. Intended for time-invariant continuous systems, the package is easily modified to include sampled-data systems. LQR designs are obtained by established techniques while the large class of suboptimal problems containing controller and/or performance index options is solved using a robust gradient minimization technique. Numerical examples demonstrate features of the package and recent developments are described.

Fleming, P. J.↗

Geometrically nonlinear analysis of laminated elastic structures

Laminated composite plates and shells that can be used to model automobile bodies, aircraft wings and fuselages, and pressure vessels among many other were analyzed. The finite element method, a numerical technique for engineering analysis of structures, is used to model the geometry and approximate the solution. Various alternative formulations for analyzing laminated plates and shells are developed and their finite element models are tested for accuracy and economy in computation. These include the shear deformation laminate theory and degenerated 3-D elasticity theory for laminates.

Reddy, J. N.↗

Core saturation in a moving medium

A numerical technique for solving the line transfer equation of a two-level atom in static equilibrium is presented. Complete redistribution of emitted photons is assumed, as is saturation at the line core. Emission intensity is calculated either by a generalized Eddington-Barber relation, a first-order differential equation for the specific intensity, or by a formal transfer integral. Sample calculations are performed of the line transfer equation in a semi-infinite atmosphere with a constant Planck function of the collision parameter and for the Mg II resonance line in a model solar atmosphere experiencing shocks. Attention is focused on the line wings in the latter problem. The first order differential equation approach yields the best intensity values and temperature structure.

Kalkofen, W.↗

Line emission from H II blister models

Numerical techniques to calculate the thermal and geometric properties of line emission from H II 'blister' regions are presented. It is assumed that the density distributions of the H II regions are a function of two dimensions, with rotational symmetry specifying the shape in three-dimensions. The thermal and ionization equilibrium equations of the problem are solved by spherical modeling, and a spherical sector approximation is used to simplify the three-dimensional treatment of diffuse ionizing radiation. The global properties of H II 'blister' regions near the edges of a molecular cloud are simulated by means of the geometry/density distribution, and the results are compared with observational data. It is shown that there is a monotonic increase of peak surface brightness from the i = 0 deg (pole-on) observational position to the i = 90 deg (edge-on) position. The enhancement of the line peak intensity from the edge-on to the pole-on positions is found to depend on the density, stratification, ionization, and electron temperature weighting. It is found that as i increases, the position of peak line brightness of the lower excitation species is displaced to the high-density side of the high excitation species.

Rubin, R. H.↗

Numerical study of the response of a compressible, turbulent boundary layer to a short region of surface curvature

Calculations were performed to investigate the supersonic flow of a turbulent boundary layer over short regions of concave surface curvature. Upstream of each curved surface the freestream Mach number was 2.9. Three different constant radii of curvature models were investigated to cover a range of curvatures and turning angles. The numerical technique solved the full, Reynolds-averaged Navier-Stokes equations using two different turbulence models: the algebraic eddy viscosity model due to Baldwin and Lomax (1978), and the one-equation model to Rubesin (1976). The calculations were compared with the experimental data of Taylor and Smits (1984), and the agreement was very satisfactory, especially for the computations using the one-equation model.

Degani, D.↗