Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “numerical method”

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 1,621 records · Page 90

Nonlinear Constitutive Relations for High Temperature Application, 1984

Nonlinear constitutive relations for high temperature applications were discussed. The state of the art in nonlinear constitutive modeling of high temperature materials was reviewed and the need for future research and development efforts in this area was identified. Considerable research efforts are urgently needed in the development of nonlinear constitutive relations for high temperature applications prompted by recent advances in high temperature materials technology and new demands on material and component performance. Topics discussed include: constitutive modeling, numerical methods, material testing, and structural applications.

Source record↗

A survey of unified constitutive theories

The state of the art of time temperature dependent elastic viscoplastic constitutive theories which are based on the unified approach werre assessed. This class of constitutive theories is characterized by the use of kinetic equations and internal variables with appropriate evolutionary equations for treating all aspects of inelastic deformation including plasticity, creep, and stress relaxation. More than 10 such unified theories which are shown to satisfy the uniqueness and stability criteria imposed by Drucker's postulate and Ponter's inequalities are identified. The theories are compared for the types of flow law, kinetic equation, evolutionary equation of the internal variables, and treatment of temperature dependence. The similarities and differences of these theories are outlined in terms of mathematical formulations and illustrated by comparisons of theoretical calculations with experimental results which include monotonic stress-strain curves, cyclic hysteresis loops, creep and stress relaxation rates, and thermomechanical loops. Numerical methods used for integrating these stiff time temperature dependent constitutive equations are reviewed.

Chan, K. S.↗

Determination of cloud parameters from infrared sounder data

The World Climate Research Programme (WCRP) plan is concerned with the need to develop a uniform global cloud climatology as part of a broad research program on climate processes. The International Satellite Cloud Climatology Project (ISCCP) has been approved as the first project of the WCRP. The ISCCP has the basic objective to collect and analyze satellite radiance data to infer the global distribution of cloud radiative properties in order to improve the modeling of cloud effects on climate. Research is conducted to explore an algorithm for retrieving cloud properties by utilizing the available infrared sounder data from polar-orbiting satellites. A numerical method is developed for computing cloud top heights, amount, and emissivity on the basis of a parameterized infrared radiative transfer equation for cloudy atmospheres. Theoretical studies were carried out by considering a synthetic atmosphere.

Yeh, H.-Y. M.↗

Fast numerical calculations of EHD sliding traction forces Application to rolling bearings

Based on analytical calculations assuming isothermal elastohydrodynamic (EHD) lubrication conditions, and on curve-fitting of a thermal correction factor, a new formula is proposed to calculate the sliding traction force developed in a concentrated contact, assuming a nonlinear, viscoelastic lubricant and undirectional sliding. When sliding occurs in the rolling and transverse directions (bidimensional sliding), a simplified numerical method is outlined by which the components of the local shear stress are quickly calculated. The latter method is applied to the calculation of the shear stress and temperature distribution in the ball-raceway contacts of an angular-contact ball bearing. Because of the viscoelastic behavior of the lubricant, a nonzero lateral traction force is obtained.

Houpert, L.↗

Dynamical relations for left ventricular ejection - Flow rate, momentum, force and impulse

An investigation was carried out to quantitatively evaluate left ventricular volume flow rate, momentum, force and impulse derived from application of conservation principles for mass and momentum of blood within the ventricle during the ejection phase. An automated digital image processing system was developed and applied to left ventricular angiograms which are computer processed and analyzed frame by frame to determine the dynamical relations by numerical methods. The initial experience with force and impulse has indicated that neither quantity seemed to be a sensitive indicator of coronary artery disease as evaluated by qualitative angiography for the particular patient group studied. Utilization of the dynamical relations in evaluating human left ventricular performance requires improved means of measurement and interpretation of clinical studies.

Back, L. H.↗

Acoustic shocks in a variable area duct containing near sonic flows

Acoustic shock waves in a variable area duct which contains near sonic flows are considered. The problem is modeled after an aeroengine inlet. Area variation of a duct and high Mach number mean the flow reduces acoustical energy yielding substantial noise reduction. One possible reason for this is acoustic shock. The use of an explicit accurate numerical method which captures shocks is described. Comparison of the results are made with an existing asymptotic theory for Mach numbers close to unity. When shock occurs reduction of sound pressure levels are shown by example.

Hariharan, S. I.↗

The evolution of large-scale structure in a universe dominated by cold dark matter

The results of numerical simulations of nonlinear gravitational clustering in universes dominated by weakly interacting, 'cold' dark matter are presented. The numerical methods used and the way in which initial conditions were generated are described, and the simulations performed are catalogued. The evolution of the fundamental statistical properties of the models is described and their comparability with observation is discussed. Graphical comparisons of these open models with the observed galaxy distribution in a large redshift survey are made. It is concluded that a model with a cosmological density parameter omega equal to one is quite unacceptable if galaxies trace the mass distribution, and that models with omega of roughly two, while better, still do not provide a fully acceptable match with observation. Finally, a situation in which galaxy formation is suppressed except in sufficiently dense regions is modelled which leads to models which can agree with observation quite well even for omega equal to one.

Davis, M.↗

Hypersonic merged-layer flow on a sphere

With the objective to achieve the desired reduction in velocity by aerodynamic forces, an aeroassisted orbital transfer vehicle (AOTV) is expected to fly in the higher region of the atmosphere for a sustained period of time. This aeroassist maneuver will occur in the transitional regime, if the vehicle is designed for a low ballistic coefficient. In this case, a merged-layer (ML) is formed around the space vehicle. The flow characteristics of the vehicle can be studied by making use of the full Navier-Stokes (NS) equations, taking into account slip and temperature jump conditions. The present investigation has the objective to obtain a two-term series solution of the full Navier-Stokes equations with surface slip and temperature jump conditions for the ML flow on a sphere. Attention is given to the mathematical formulation of the problem, the numerical method of integration, and the results.

Jain, A. C.↗

The exact eigenfunctions and eigenvalues of a two-dimensional rigid rotor obtained using Gaussian wave packet dynamics

Exact eigenfunctions for a two-dimensional rigid rotor are obtained using Gaussian wave packet dynamics. The wave functions are obtained by propagating, without approximation, an infinite set of Gaussian wave packets that collectively have the correct periodicity, being coherent states appropriate to this rotational problem. This result leads to a numerical method for the semiclassical calculation of rovibrational, molecular eigenstates. Also, a simple, almost classical, approximation to full wave packet dynamics is shown to give exact results: this leads to an a posteriori justification of the De Leon-Heller spectral quantization method.

Reimers, J. R.↗

Numerical simulation of compressible, turbulent, two-phase flow

A computer program for numerically simulating compressible, turbulent, two-phase flows is described and applied. Special attention is given to flows in which dust is ingested into the turbulent boundary layer behind shock waves moving over the earth's surface. it is assumed that the two phases are interpenetrating continua which are coupled by drag forces and heat transfer. The particle phase is assumed to be dilute, and turbulent effects are modeled by zero- and two-equation eddy viscosity models. An important feature of the turbulence modeling is the treatment of surface boundary conditions which control the ingestion of particles into the boundary layer by turbulent friction and diffusion. The numerical method uses second-order implicit upwind differencing of the inviscid terms of the equations and second-order central differencing of the viscous terms. A diagonal form of the implicit algorithm is used to improve efficiency, and the transformation to a curvilinear coordinate system is accomplished by the finite volume techniques. Applications to a series of representative flows include a two-phase nozzle flow, the steady flow of air over a sand bed, and the air flow behind a normal shock wave in uniform motion over a sand bed. Results of the latter two applications are compared with experimental results.

Coakley, t. J.↗

Higher-order derivative correlations and the alignment of small-scale structures in isotropic numerical turbulence

The classical approach to the investigation of small-scale intermittency in turbulence is based on the higher-order derivative correlations such as skewness and flatness factors. In the study of the small scales, numerial simulations can provide more detail than experiments. In the present paper, a variety of velocity- and scalar-derivative correlations are calculated over a range of Reynolds numbers. Particular attention is given to third- and fourth-order correlations, taking into account also some fifth- and sixth-order correlations to allow comparisons with the phenomenological models. The governing equations are the incompresssible Navier-Stokes equation for the velocity and transport equation for a passive scalar. Two numerical codes are used for the simulations presented. Attention is given to details regarding the numerical method used, forcing, simulation parameters, spectra and skewnesses, and graphics.

Kerr, R. M.↗

Nonlinear Supersonic Full Potential Analysis

Supersonic Implicit Marching Program (SIMP) applies numerical method, based on conservative form of full potential equation, to problem of threedimensional supersonic flows with embedded subsonic regions. Conservative formulation of problem provides ability to capture shocks and to assess accurately impact of sweep, thickness, and lift for conditions where linear theory unsatisfactory. Technique uses characteristic signal-propagation theory to control density biasing for treatment of shocks (including embedded shocks) and mixed elliptic/hyperbolic crossflow. SIMP written in FORTRAN 77.

Shankar, V.↗

Numerical procedure for three-dimensional hypersonic viscous flow over aerobrake configuration

A numerical method, which is simpler and more efficient than others currently in use, is proposed for the computation of the full viscous flow over an aerobrake body in hypersonic stream at high altitude. It treats the shock layer surrounding the blunt forebody and the near wake behind the base simultaneously by formulating the Navier-Stokes equations in conformal and azimuthal-angle coordinates. The computational domain is confined by the body wall, outflow surface and the shock, which is adjusted along the coordinate normal to the wall in the course of iterations. Because of the optimal grid and a well developed alternating direction implicit factorization technique for the governing equations, reasonably accurate results can be obtained with a 28 x 36 x 7 grid and 400 time-marching iterations. Excellent agreement of shock location is found between the present result and the schlieren photograph. Details of the base flow and shear layer impingement on the cylindrical aft body are presented for an adiabatic wall case.

Li, C. P.↗

Theoretical and experimental flow fields for a supersonic cruise fighter forebody

The capability of two numerical methods to predict the flow field about a representative supersonic cruise aircraft was examined. The two codes were a small-disturbance transonic program of Boppe and a conservative form full-potential transonic program of Shankar and Szema. For the former code, comparisons were made with wind-tunnel data at Mach numbers of 0.6, 0.9, and 1.2 and angles of attack of 0 deg, 5 deg, and 10 deg (7.5 deg instead of 10 deg at a Mach number of 1.2). Predictions from the two codes were compared at a Mach number of 1.2 and an angle of attack of 7.5 deg. The comparison criteria were contours of local angle of attack, local angle of sideslip, and local Mach number. The comparisons indicated that both codes may be considered useful for design applications, depending on the degree of accuracy required by the user's solution. Both solutions show an inaccuracy in their predictions, particularly as Mach number and/or angle of attack increases, because of their lack of viscous effects and any mechanisms to predict vortex development.

Yaros, S. F.↗

An LU implicity scheme for high speed inlet analysis

A numerical method is developed to analyze the inviscid flowfield of a high speed inlet by the solution of the Euler equations. The lower-upper implicit scheme in conjunction with adaptive dissipation proves to be an efficient and robust nonoscillatory shock capturing technique for high Mach number flows as well as for transonic flows.

Yoon, S.↗

Simulation of the turbulent Rayleigh-Benard problem using a spectral/finite difference technique

The three-dimensional, incompressible Navier-Stokes and energy equations with the Bousinesq assumption have been directly simulated at a Rayleigh number of 3.8 x 10 to the 5th power and a Prandtl number of 0.76. In the vertical direction, wall boundaries were used and in the horizontal, periodic boundary conditions were used. A spectral/finite difference numerical method was used to simulate the flow. The flow at these conditions is turbulent and a sufficiently fine mesh was used to capture all relevant flow scales. The results of the simulation are compared to experimental data to justify the conclusion that the small scale motion is adequately resolved.

Eidson, T. M.↗

Advances in numerical and applied mathematics

This collection of papers covers some recent developments in numerical analysis and computational fluid dynamics. Some of these studies are of a fundamental nature. They address basic issues such as intermediate boundary conditions for approximate factorization schemes, existence and uniqueness of steady states for time dependent problems, and pitfalls of implicit time stepping. The other studies deal with modern numerical methods such as total variation diminishing schemes, higher order variants of vortex and particle methods, spectral multidomain techniques, and front tracking techniques. There is also a paper on adaptive grids. The fluid dynamics papers treat the classical problems of imcompressible flows in helically coiled pipes, vortex breakdown, and transonic flows.

South, J. C., Jr.↗

A model for the kinetics of a solar-pumped long path laser experiment

A kinetic model for a solar-simulator pumped iodine laser system is developed and compared to an experiment in which the solar simulator output is dispersed over a large active volume (150 cu cm) with low simulator light intensity (approx. 200 solar constants). A trace foreign gas which quenches the upper level is introduced into the model. Furthermore, a constant representing optical absorption of the stimulated emission is introduced, in addition to a constant representing the scattering at each of the mirrors, via the optical cavity time constant. The non-uniform heating of the gas is treated as well as the pressure change as a function of time within the cavity. With these new phenomena introduced into the kinetic model, a best reasonable fit to the experimental data is found by adjusting the reaction rate coefficients within the range of known uncertainty by numerical methods giving a new bound within this range of uncertainty. The experimental parameters modeled are the lasing time, laser pulse energy, and time to laser threshold.

Stock, L. V.↗