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

A computer program for wing subsonic aerodynamic performance estimates including attainable thrust and vortex lift effects

Numerical methods incorporated into a computer program to provide estimates of the subsonic aerodynamic performance of twisted and cambered wings of arbitrary planform with attainable thrust and vortex lift considerations are described. The computational system is based on a linearized theory lifting surface solution which provides a spanwise distribution of theoretical leading edge thrust in addition to the surface distribution of perturbation velocities. The approach used relies on a solution by iteration. The method also features a superposition of independent solutions for a cambered and twisted wing and a flat wing of the same planform to provide, at little additional expense, results for a large number of angles of attack or lift coefficients. A previously developed method is employed to assess the portion of the theoretical thrust actually attainable and the portion that is felt as a vortex normal force.

Carlson, H. W.↗
NASA NTRS1982

Geometry and starvation effects in hydrodynamic lubrication

Numerical methods were used to determine 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. E.↗
NASA NTRS1981

Computation of wing-vortex interaction in transonic flow using implicit finite difference algorithm

An implicit delta form finite difference algorithm for Euler equations in conservation law form was used in preliminary calculations of three dimensional wing vortex interaction. Both steady and unsteady transonic flow wing vortex interactions are computed. The computations themselves are meant to guide upcoming wind tunnel experiments of the same flow field. Various modifications to the numerical method that are intended to improve computational efficiency are also described and tested in both two and three dimensions. Combination of these methods can reduce the overall computational time by a factor of 4.

Srinivasan, G.↗
NASA NTRS1982

Numerical simulation of one-dimensional heat transfer in composite bodies with phase change

A numerical simulation was developed to investigate the one dimensional heat transfer occurring in a system composed of a layered aircraft blade having an ice deposit on its surface. The finite difference representation of the heat conduction equations was done using the Crank-Nicolson implicit finite difference formulation. The simulation considers uniform or time dependent heat sources, from heaters which can be either point sources or of finite thickness. For the ice water phase change, a numerical method which approximates the latent heat effect by a large heat capacity over a small temperature interval was applied. The simulation describes the temperature profiles within the various layers of the de-icer pad, as well as the movement of the ice water interface. The simulation could also be used to predict the one dimensional temperature profiles in any composite slab having different boundary conditions.

Dewitt, K. J.↗
NASA NTRS1982

Computer Science Research at Langley

A workshop was held at Langley Research Center, November 2-5, 1981, to highlight ongoing computer science research at Langley and to identify additional areas of research based upon the computer user requirements. A panel discussion was held in each of nine application areas, and these are summarized in the proceedings. Slides presented by the invited speakers are also included. A survey of scientific, business, data reduction, and microprocessor computer users helped identify areas of focus for the workshop. Several areas of computer science which are of most concern to the Langley computer users were identified during the workshop discussions. These include graphics, distributed processing, programmer support systems and tools, database management, and numerical methods.

Voigt, S. J.↗
NASA NTRS1982

Fluid flow in crystallization processes

Investigations of the interactions of heat, mass, and momentum transport in crystal growth from the melt are described. Studies of the mall-scale floating zone process and on a prototype of the vertical Bridgman growth system were emphasized. In both systems detailed numerical calculations are used to dissect the interplay between fluid convection and dopant segregation. These calculations are based on finite element techniques that make feasible the complete solution of solidification problems which include convection. Other numerical methods were developed for solving thermal models of crystal growth processes with melt/gas menisci and for simulating the microscale instabilities in solidification interfaces.

Brown, R. A.↗
NASA NTRS1982

Three dimensional flow field inside compressor rotor, including blade boundary layers

The Reynolds stress equation, pressure strain correlation, and dissipative terms and diffusion are discussed in relation to turbulence modelling using the Reynolds stress model. Algebraic modeling of Reynolds stresses and calculation of the boundary layer over an axial cylinder are examined with regards to the kinetic energy model for turbulence modelling. The numerical analysis of blade and hub wall boundary layers, and an experimental study of rotor blade boundary layer in an axial flow compressor rotor are discussed. The Patankar-Spalding numerical method for two dimensional boundary layers is included.

Galmes, J. M.↗
NASA NTRS1982

Numerical simulation of wall-bounded turbulent shear flows

Developments in three dimensional, time dependent numerical simulation of turbulent flows bounded by a wall are reviewed. Both direct and large eddy simulation techniques are considered within the same computational framework. The computational spatial grid requirements as dictated by the known structure of turbulent boundary layers are presented. The numerical methods currently in use are reviewed and some of the features of these algorithms, including spatial differencing and accuracy, time advancement, and data management are discussed. A selection of the results of the recent calculations of turbulent channel flow, including the effects of system rotation and transpiration on the flow are included.

Moin, P.↗
NASA NTRS1982

Theoretical analyses of Baroclinic flows

A stability analysis of a thin horizontal rotating fluid layer which is subjected to arbitrary horizontal and vertical temperature gradients is presented. The basic state is a nonlinear Hadley cell which contains both Ekman and thermal boundary layers; it is given in closed form. The stability analysis is based on the linearized Navier-Stokes equations, and zonally symmetric perturbations in the form of waves propagating in the meridional direction are considered. Numerical methods were used for the stability problem. It was found that the instability sets in when the Richardson number is close to unity and that the critical Richardson number is a non-monotonic function of the Prandtl number. Further, it was found that the critical Richardson number decreases with increasing Ekman number until a critical value of the Ekman number is reached beyond which the fluid is stable.

Antar, B.↗
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.↗