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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,567 records · Page 87

All sources NASA NTRS DOE OSTI

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.↗

NASA NTRS1983

Numerical simulation of separated flows

A new numerical method, based on the Vortex Method, for the simulation of two-dimensional separated flows, was developed and tested on a wide range of gases. The fluid is incompressible and the Reynolds number is high. A rigorous analytical basis for the representation of the Navier-Stokes equation in terms of the vorticity is used. An equation for the control of circulation around each body is included. An inviscid outer flow (computed by the Vortex Method) was coupled with a viscous boundary layer flow (computed by an Eulerian method). This version of the Vortex Method treats bodies of arbitrary shape, and accurately computes the pressure and shear stress at the solid boundary. These two quantities reflect the structure of the boundary layer. Several versions of the method are presented and applied to various problems, most of which have massive separation. Comparison of its results with other results, generally experimental, demonstrates the reliability and the general accuracy of the new method, with little dependence on empirical parameters. Many of the complex features of the flow past a circular cylinder, over a wide range of Reynolds numbers, are correctly reproduced.

Spalart, P. R.↗

NASA NTRS1983

Conceptual design of a data reduction system

A telemetry data processing system was defined of the Data Reduction. Data reduction activities in support of the developmental flights of the Space Shuttle were used as references against which requirements are assessed in general terms. A conceptual system design believed to offer significant throughput for the anticipated types of data reduction activities is presented. The design identifies the use of a large, intermediate data store as a key element in a complex of high speed, single purpose processors, each of which performs predesignated, repetitive operations on either raw or partially processed data. The recommended approach to implement the design concept is to adopt an established interface standard and rely heavily on mature or promising technologies which are considered main stream of the integrated circuit industry. The design system concept, is believed to be implementable without reliance on exotic devices and/or operational procedures. Numerical methods were employed to examine the feasibility of digital discrimination of FDM composite signals, and of eliminating line frequency noises in data measurements.

Source record↗

NASA NTRS1983

Feasibility study of an optically coherent telescope array in space

Numerical methods of image construction which can be used to produce very high angular resolution images at optical wavelengths of astronomical objects from an orbiting array of telescopes are discussed and a concept is presented for a phase-coherent optical telescope array which may be deployed by space shuttle in the 1990's. The system would start as a four-element linear array with a 12 m baseline. The initial module is a minimum redundant array with a photon-counting collecting area three times larger than space telescope and a one dimensional resolution of better than 0.01 arc seconds in the visible range. Later additions to the array would build up facility capability. The advantages of a VLBI observatory in space are considered as well as apertures for the telescopes.

Traub, W. A.↗

NASA NTRS1983

Small Gas Turbine Combustor Primary Zone Study

A development process is described which consists of design, fabrication, and preliminary test evaluations of three approaches to internal aerodynamic primary zone flow patterns: (1) conventional double vortex swirl stabilization; (2) reverse flow swirl stabilization; and (3) large single vortex flow system. Each concept incorporates special design features aimed at extending the performance capability of the small engine combustor. Since inherent geometry of these combustors result in small combustion zone height and high surface area to volume ratio, design features focus on internal aerodynamics, fuel placement, and advanced cooling. The combustors are evaluated on a full scale annular combustor rig. A correlation of the primary zone performance with the overall performance is accomplished using three intrusion type gas sampling probes located at the exit of the primary zone section. Empirical and numerical methods are used for designing and predicting the performance of the three combustor concepts and their subsequent modifications. The calibration of analytical procedures with actual test results permits an updating of the analytical design techniques applicable to small reverse flow annular combustors.

Sullivan, R. E.↗

NASA NTRS1983

A mathematical characterization of vegetation effect on microwave remote sensing from the Earth

In passive microwave remote sensing of the earth, a theoretical model that utilizes the radiative transfer equations was developed to account for the volume scattering effects of the vegetation canopy. Vegetation canopies such as alfalfa, sorghum, and corn are simulated by a layer of ellipsoidal scatterers and cylindrical structures. The ellipsoidal scatterers represent the leaves of vegetation and are randomly positioned and oriented. The orientation of ellipsoids is characterized by a probability density function of Eulerian angles of rotation. The cylindrical structures represent the stalks of vegetation and their radii are assumed to be much smaller than their lengths. The underlying soil is represented by a half-space medium with a homogeneous permittivity and uniform temperature profile. The radiative transfer quations are solved by a numerical method using a Gaussian quadrature formula to compute both the vertical and horizontal polarized brightness temperature as a function of observation angle. The theory was applied to the interpretation of experimental data obtained from sorghum covered fields near College Station, Texas.

Choe, Y.↗

NASA NTRS1982

A time-dependent ice sheet model - Preliminary results

A numerical model of ice sheet flow is developed, and preliminary results are described. This model includes vertical resolution of temperature, stress, and strain rate which represents a considerable improvement over previous vertically averaged ice sheet models. The model follows the flow of ice along a flow line within an ice sheet drainage basin. Longitudinal stresses and basal sliding are included. Basal sliding is dependent on the base shear stress and a specified distribution of basal water pressure. The numerical methods used to solve the coupled set of stress and velocity equations for the static and time-evolutionary cases are discussed. A steady state profile simulating an ice stream is calculated for a particular set of input parameters, and changes in the profile are examined for different choices of parameters. Preliminary studies of response behavior are completed using a simplified ice sheet geometry with a fixed terminus or grounding line. The results of these studies illustrate ice sheet thinning in response to a lowered sea level or to a reduction in the extent of ice rises (or pinning points) within ice shelves.

Bindschadler, R. A.↗

NASA NTRS1982

Noise transmission into semicylindrical enclosures through discretely stiffened curved panels

An analytical study of sound transmission into semicylindrical enclosures through discretely stiffened curved elastic panels is presented. The transmitted sound is estimated by solving the acoustic wave equation for the interior acoustic field, a Galerkin-like method being used. This solution is then coupled to the vibration of the stiffened panels. The response characteristics of these panels are determined by using a modal analysis where the modes are obtained by the finite element-strip method. Numerical results include spectra of the interior sound pressure due to white noise, turbulent boundary layer and propeller noise inputs.

Chang, M. T.↗

NASA NTRS1983

Real gas flow fields about three dimensional configurations

Real gas, inviscid supersonic flow fields over a three-dimensional configuration are determined using a factored implicit algorithm. Air in chemical equilibrium is considered and its local thermodynamic properties are computed by an equilibrium composition method. Numerical solutions are presented for both real and ideal gases at three different Mach numbers and at two different altitudes. Selected results are illustrated by contour plots and are also tabulated for future reference. Results obtained compare well with existing tabulated numerical solutions and hence validate the solution technique.

Balakrishnan, A.↗

NASA NTRS1982

Thermoelastohydrodynamic analysis of an oil pumping ring seal

The selection and evaluation of the numerical methods used in the thermoelastohydrodynamic analysis of an oil pumping ring seal are described, and the results from test problems used to demonstrate the validity of the implementing computer program are presented. A sampling of predicted pumping ring performance values is compared with existing experimental data. Agreement between prediction and test performance is obtained, yielding detailed information concerning fluid and solid parameters.

Smith, P. J.↗