Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Numerical Methods”

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,081 records · Page 60
All sourcesNASA NTRSDOE OSTI
NASA NTRS1982

Computation of subsonic flow around airfoil systems with multiple separation

A numerical method for computing the subsonic flow around multi-element airfoil systems was developed, allowing for flow separation at one or more elements. Besides multiple rear separation also sort bubbles on the upper surface and cove bubbles can approximately be taken into account. Also, compressibility effects for pure subsonic flow are approximately accounted for. After presentation the method is applied to several examples and improved in some details. Finally, the present limitations and desirable extensions are discussed.

Jacob, K.↗
NASA NTRS1981

Quasi-Lagrangian models of nascent thermals

Slab-symmetric and axisymmetric geometries are considered in the modeling of thermals by considering their boundaries as self-advecting vortex sheets, for the case of a neutral environment without turbulent processes. The initial accelerations for circular cross section thermals are obtained analytically, and a numerical method is developed to simulate the evolution of thermals by time-marching in which the vortex sheet is divided into finite segments whose positions are tracked in a quasi-Lagrangian fashion. Self-advection is considered, and a redistribution procedure is adopted to prevent the segments from bunching unmanageably. The induced field of motion is determined both inside and outside the thermal. It is shown that (1) the axisymmetric thermal rises more quickly than the slab-symmetric thermal, (2) vorticity maxima occur both within the concavities and on the trailing edges, and (3) the leading edges are very smooth.

Rambaldi, S.↗
NASA NTRS1981

Approximate and exact numerical computation of supersonic flow over an airfoil

The computation of inviscid supersonic flow over a two-dimensional airfoil is considered. There are two main nonlinear approaches which lead to approximate solutions. Small-amplitude theory gives valid solutions provided the airfoil thickness is not too great and the Mach number is not too high. The second type of approximation, shock expansion theory, employs the fact that even for flows with strong shocks the effect of the down-running characteristics remains small. This leads to an analytic solution at the airfoil. It is pointed out that the approximate solution is accurate enough for many cases of interest. The numerical method furnishes a rapid correction to the solution in those cases where it is not. The characteristic-streamline coordinate system is useful both for the computation of the approximate solution and the corrections, and is also convenient for displaying and interpreting the results.

Lewis, T. S.↗
NASA NTRS1981

Direct numerical simulations of turbulent shear flows

Numerical simulations of wakes of axisymmetric bodies and of turbulent mixing layers are reported. The flows were assumed to be statistically homogeneous in the mean flow direction, in concert with experimental data and the self-similarity theorem. The nonlinear Navier-Stokes equations were solved by a pseudo-spectral numerical method using a 32 x 32 x 33 point grid and an algorithm for fast Fourier transforms and inverse transforms. Leapfrog time differencing was employed on nonlinear terms and time differencing on viscous terms. Towed wakes and wakes behind a self-propelled body were simulated, showing that the towed wakes exhibited a proper temporal behavior after an initial period of adjustment, including the development of a kurtosis near the wake edge, which is experimentally verifiable. The mixing-layer simulation displayed the laboratory demonstrated presence of large scale features such as vortex cores, while the lateral coherence was weak.

Metcalfe, R. W.↗
NASA NTRS1982

Solar particle propagation and modulation theory

Developments related to solar flare particle propagation theory are considered, taking into account the transport equation currently employed to describe the propagation of solar flare particles. This equation does not take into consideration convection with the solar wind, drifts, and changes in particle energy. The equation was solved numerically to generate typical solar flare particle time profiles at 1 AU. The results, however, were not satisfactory. It was found that the numerical method used was apparently in error. New advances which provide approximate analytical solutions to the equation are reported. Questions of galactic cosmic ray propagation in interplanetary space are also discussed.

Lee, M. A.↗
NASA NTRS1982

Numerical experimentation on spherically symmetric one-dimensional magnetohydrodynamic /MHD/ wave propagation

Radial propagation of one-dimensional magnetohydrodynamic (MHD) waves are analyzed numerically on the basis of the Implicit-Continuous-Fluid-Eulerian (ICE) scheme. Accuracy of the numerical method and other properties are tested through the study of MHD wave propagation. The three different modes of MHD waves (i.e., fast-, slow- and Alfven (transverse) mode) are generated by applying physically consistent boundary perturbations derived from MHD compatibility relations. It is shown that the resulting flow following these waves depend upon the relative configurations of the initial magnetic field and boundary perturbations.

Han, S. M.↗
NASA NTRS1981

Numerical studies of the energy balance in coronal loops

A numerical method is applied to treat the energy balance of quasi-static solar coronal loops, which have been observed to persist for periods much greater than the radiative cooling time. The quasi-static loop model employed takes into account gravity, density-, temperature- or position-dependent energy input, an accurate form of the radiative losses and variable loop cross-sectional area, under assumptions of energy input by coronal heating balanced by radiative and conductive losses, an optically thin plasma, energy conduction along the field lines only and hydrostatic equilibrium. Computations of an emission measure function for various distributions of the energy input and loop geometries are then presented which show that little information on the location of the energy input may be gained from spectral line intensity measurements integrated over a single loop.

Underwood, J. H.↗
NASA NTRS1982

Reynolds averaged Navier-Stokes computations of transonic flows The-state-of-the-art

The state of the art in computer simulations for transonic flowfields requiring solutions for the Navier-Stokes equations is assessed. It is noted that current simulations of transonic flowfields require comparisons with experimental results because the simulations are not free from discretization errors. Cases of turbulence are treated with weighted variables in a time-averaged scheme to yield Reynolds averaged Navier-Stokes equations. The turbulence is modeled in a first-order approach with a Reynolds stress tensor or a second-order approach where the tensor is obtained from the Navier-Stokes equations. The choice of coordinate systems and dependent variables for simulating the flow around airfoils is discussed, with particular attention to curvilinear coordinates. The determinations of boundary conditions is examined, along with numerical methods related to physical phenomena

Mehta, U.↗
NASA NTRS1982

A numerical experiment on the equilibrium and stability of a rotating galactic bar

A self-consistent, three-dimensional numerical experiment is performed on an N-body system whose initial state is a realization of a certain theoretical model of a rotating triaxial galaxy. The model is a stellar-dynamical counterpart of a uniformly rotating polytrope of index equal to 0.5. The aim of the experiment is to study the equilibrium of the system and, in particular, to test its stability. The experimental system behaves in the mean like a realization of the theoretical model for at least seven crossing times. The principal departure of the system from equilibrium is an oscillation which is identified as a radial pulsation. There is no indication in its behavior that the system is unstable with respect to anu mode with an e-folding time shorter than or of the order of two crossing times. Certain changes that occur in the state of the system are interpreted, with the aid of the theoretical model, as secular changes which result from a slight failure of our numerical methods to conserve the mass, energy, and angular momentum of the system; these effects are small enough that they do not vitiate the experiment on a dynamical time scale.

Miller, R. H.↗
NASA NTRS1982

The instability of sheared liquid layers

The critical values of the Reynolds number (Rc) above which unstable travelling waves exist in a steady shear flow are calculated with the inclusion of a term in the normal stress boundary condition omitted by Miles (1960). Results show a value of Rc = 34.2, compared to Miles' value of 203, for the case of a zero-surface tension interface. The critical wavenumber, critical phase speed, and Rc are calculated as functions of a nondimensional surface tension, and the mechanism of the instability is investigated by examining the disturbance energy equation. A long wave instability is found at small values of the Reynolds number when the shear flow has a parabolic velocity profile. In addition, these results are extended to larger values of the wavenumber by numerical methods, and the relation between this long wave instability and the profile curvature is examined.

Smith, M. K.↗
NASA NTRS1982

An analytic approach to two-fixed-impulse transfers between Keplerian orbits

Solutions are obtained for the two-impulse transfer of a vehicle between arbitrary inclined orbits in an inverse-square force field with the restriction that the magnitude of each of the two impulses has a fixed preassigned value. The two magnitudes need not be equal. The equations for the conservation of angular momentum and energy are augmented by the Laplace integral; these equations establish linear relationships between several of the variables. This set of linear equations and one of two quadratic equations constitute the analytically tractable part of the solution. The remaining part, consisting of finding the zeros (if they exist) of a single trigonometric function of one variable, is solved using numerical methods. Explicit lower bounds on each of the magnitudes of the two impulses are obtained by requiring the solution be real. Graphical results are presented to illustrate the solution.

Jezewski, D. J.↗
NASA NTRS1984

Simulating a Three-Dimensional Flow in Pipes

Numerical method simulates in unsteady, incompressible flow in axisymmetric pipe in three dimensions. Method relies on velocity-vector expansion technique in which each vector function in expansion set is divergence-free and satisfies boundary conditions for viscous flow.

Leonard, A.↗
NASA NTRS1981

Higher level simulations of turbulent flows

The fundamentals of large eddy simulation are considered and various approaches to this simulation are compared. The subgrid model required by large eddy simulation is discussed as well as the use of this type of simulation in the development of models for the Reynolds-averaged equations and the application of direct simulation to the testing of both subgrid scale and Reynolds-averaged models. Numerical methods used in large eddy and direct simulation are described with emphasis on special purpose methods. Topics covered include the simulation of homogeneous flow, free shear flows, the mixing layer, wakes, and wall-bounded flows including channel flow and the boundary layer. Applications of large eddy simulation in the laboratory as in meteorological and other environmental flows are examined. Directions in which the work is proceeding and what can be expected from higher levels simulated on are examined.

Ferziger, J. H.↗
NASA NTRS1982

Response studies of rotors and rotor blades with application to aeroelastic tailoring

Various tools for the aeroelastic stability and response analysis of rotor blades in hover and forward flight were developed and incorporated in a comprehensive package capable of performing aeroelastic tailoring of rotor blades in forward flight. The results indicate that substantial vibration reductions, of order 15-40%, in the vibratory hub shears can be achieved by relatively small modifications of the initial design. Furthermore the optimized blade can be up to 20% lighter than the original design. Accomplishments are reported for the following tasks: (1) finite element modeling of rotary-wing aeroelastic problems in hover and forward flight; (2) development of numerical methods for calculating the aeroelastic response and stability of rotor blades in forward fight; (3) formulation of the helicopter air resonance problem in hover with active controls; and (4) optimum design of rotor blades for vibration reduction in forward flight.

Friedmann, P. P.↗
NASA NTRS1982

Geometry and starvation effects in hydrodynamic lubrication

Numerical methods were used to detemine the effects of lubricant starvation on the minimum film thickness under conditions of a hydrodynamic point contact. Starvation was effected by varying the fluid inlet level. The Reynolds boundary conditions were applied at the cavitation boundary and zero pressure was stipulated at the meniscus or inlet boundary. A minimum film thickness equation as a function of both the ratio of dimensionless load to dimensionless speed and inlet supply level was determined. By comparing the film generated under the starved inlet condition with the film generated from the fully flooded inlet, an expression for the film reduction factor was obtained. Based on this factor a starvation threshold was defined as well as a critically starved inlet. The changes in the inlet pressure buildup due to changing the available lubricant supply are presented in the form of three dimensional isometric plots and also in the form of contour plots.

Brewe, D.↗
NASA NTRS1982

Fitting surfaces to scattered data

A variety of numerical methods for fitting a function to data given at a set of points scattered throughout a domain in the plane are surveyed. Four classes of methods are discussed: (1) global interpolation; (2) local interpolation; (3) global approximation; and (4) local approximation. Also, two-stage methods and contouring are discussed. The surfaces constructed include polynomials, spline functions, and rational functions, among others.

Schumaker, L. L.↗
NASA NTRS1983

Prediction of high speed propeller flow fields using a three-dimensional Euler analysis

To overcome the limitations of classical propeller theory, a computer program, NASPROP-E, was developed which solves for the flow field surrounding a multibladed propeller and axisymmetric nacelle combination using a finite difference method. The governing equations are the three dimensional unsteady Euler equations written in a cylindrical coordinate system. They are marched in time until a steady state solution is obtained. The Euler equations require no special treatment to model the blade work vorticity. The equations are solved using an implicit approximate factorization method. Numerical results are presented which have greatly increased the understanding of high speed propeller flow fields. Numerical results for swirl angle downstream of the propeller and propeller power coefficient are higher than experimental results. The radial variation of coefficient are higher than experimental results. The radial variation of swirl angle, however, is in reasonable agreement with the experimental results. The predicted variation of power coefficient with blade angle agrees very well with data.

Bober, L. J.↗