SEARCH · Engineering Papers

Results for “THREE-DIMENSIONAL CALCULATIONS”

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 451 records · Page 25

All sources NASA NTRS DOE OSTI

NASA NTRS1980

On the stability of three-dimensional compressible nonparallel boundary layers

A compressible linear stability theory is presented for three-dimensional nonuniform boundary layers. The amplitude, phase, and wavenumber equations which govern the motion of the disturbance are obtained by using the method of multiple scales. Group velocity trajectories are used to identify the disturbance growth direction. The spatial stability theory is applied to the flow on a laminar flow control supercritical sweptback wing of infinite span. Three different methods are used to calculate the absolute maximum logarithmic amplitude ratio N. In the front crossflow instability region, the method of maximum spatial growth rate predicts large difference in N compared with the method of fixed wavelength and the method of fixed spanwise component of wavelength. This difference decreases in the middle streamwise instability region, and almost vanishes in the rear crossflow instability region. Compressibility of the medium reduces N by about 15% in both the front and rear regions, and by about 40% in the middle region of the wing. Nonuniformity of the medium has large effects specially in the rear region.

El-Hady, N. M.↗

NASA NTRS1980

Illustration of airfoil shape effect on forward-swept wing divergence

A static aeroelastic analysis is presented of the divergence of untapered wings with conventional and supercritical airfoil sections at sweep angles of zero and -15 deg. One bending and one torsion mode were employed for a uniform rectangular cantilevered beam with the elastic axis at midchord, and calculations were based on a two-dimensional differential equations formulation in the structural coordinate system and in simple strip theory. A minimum divergence speed in the transonic range is obtained which is associated with the rearward shift of the aerodynamic center, and a 17% difference in minimum divergence dynamic pressure is found between a supercritical and a conventional wing. It is noted that although the strip method employed allows the assessment of the sensitivity of airfoil shapes to divergence, three-dimensional transonic aerodynamic methods should be used to predict wing divergence characteristics.

Bland, S. R.↗

NASA NTRS1981

Numerical simulation of flows in curved diffusers with cross-sectional transitioning using a three-dimensional viscous analysis

A three dimensional analysis for fully viscous, subsonic, compressible flow is evaluated. An approximate form of the Navier Stokes equations is solved by an implicit spatial marching technique. Calculations were made for flow in a circular S duct and in the F 16 inlet duct. The computed total pressure contours and secondary flow velocity vectors are presented. Qualitative comparisons with experiment are shown for both ducts. The analysis is used to show how the cross section transitioning in the F 16 inlet suppresses the development of a secondary flow vortex.

Towne, C. E.↗

NASA NTRS1980

Full-coverage film cooling. I - Three-dimensional measurements of turbulence structure. II - Prediction of the recovery-region hydrodynamics

Hydrodynamic measurements of turbulence structure were performed with a triaxial hot wire in the full coverage and the recovery regions following an array of injection holes under isothermal conditions at ambient temperature and pressure for blowing ratios of 0.9 and 0.4. High levels of turbulence kinetic energy (TKE) were determined for low blowing, and low TKE levels were found for the high blowing levels; in the recovery region, the flow can be represented by a model with an outer boundary layer and a 2-dimensional inner boundary layer. Recovery region hydrodynamics can be modelled by considering that a new boundary layer started to grow immediately after the end of blowing; the Prandtl mixing length distributions calculated from the values of mean velocity and turbulent shear stresses were consistent with the presence of a dual boundary layer structure in the recovery region. The program used here contains a one-equation model of turbulence, using turbulence kinetic energy with an algebraic mixing length; this 2-dimensional, finite difference program can predict the mean velocity and turbulence kinetic energy profiles based on initial values, boundary conditions, and a closure condition.

Yavuzkurt, S.↗

NASA NTRS1981

Development of simplified airborne computations for fuel conservative descents in a time-based metered air traffic environment

The NASA has developed and flight-tested a simple flight management descent algorithm designed to improve the accuracy of delivering an airplane in a fuel-conservative manner to a metering fix at a time designated by air traffic control. This algorithm provides a three-dimensional path with terminal area time constraints (four-dimensional) for an airplane to make an idle-thrust, clean-configured (landing gear up, flaps zero, and speed brakes retracted) descent to arrive at the metering fix at a predetermined time, altitude, and airspeed. The descent path is calculated for a constant Mach/airspeed schedule from linear approximations of airplane performance with considerations given for gross weight, wind, and nonstandard pressure and temperature effects. Applications of the four-dimensional and descent planning capabilities of the algorithm to conventional airplanes is being investigated. This report describes the flight management descent algorithm and presents the results of the flight tests flown with the Terminal Configured Vehicle airplane.

Knox, C. E.↗

NASA NTRS1982

Three-dimensional relativistic field-electron interaction in a multicavity high-power klystron. 1: Basic theory

A theoretical investigation of three dimensional relativistic klystron action is described. The relativistic axisymmetric equations of motion are derived from the time-dependent Lagrangian function for a charged particle in electromagnetic fields. An analytical expression of the fringing RF electric and magnetic fields within and in the vicinity of the interaction gap and the space-charge forces between axially and radially elastic deformable rings of charges are both included in the formulation. This makes an accurate computation of electron motion through the tunnel of the cavities and the drift tube spaces possible. Method of analysis is based on Lagrangian formulation. Bunching is computed using a disk model of electron stream in which the electron stream is divided into axisymmetric disks of equal charge and each disk is assumed to consist of a number of concentric rings of equal charges. The Individual representative groups of electrons are followed through the interaction gaps and drift tube spaces. Induced currents and voltages in interacting cavities are calculated by invoking the Shockley-Ramo theorem.

Kosmahl, H. G.↗

NASA NTRS1981

Numerical Modeling of Three-Dimensional Confined Flows

A three dimensional confined flow model is presented. The flow field is computed by calculating velocity and enthalpy along a set of streamlines. The finite difference equations are obtained by applying conservation principles to streamtubes constructed around the chosen streamlines. With appropriate substitutions for the body force terms, the approach computes three dimensional magnetohydrodynamic channel flows. A listing of a computer code, based on this approach is presented in FORTRAN IV language. The code computes three dimensional compressible viscous flow through a rectangular duct, with the duct cross section specified along the axis.

Greywall, M. S.↗

NASA NTRS1982

Computation of incompressible, three-dimensional turbulent boundary layers and comparison with experiment

Incompressible three dimensional, turbulent boundary layer (3DTBL) experiments were simulated numerically by integrating the boundary layer equations together with an algebraic eddy viscosity turbulence model. For the flow treated, the downstream portion, where the crossflow was large, was not predicted with the present computational method; the flow was significantly influenced by elliptic flow field effects. Departures from the boundary layer concept are indicated. Calculations agreed reasonably well with the mean flow development up to separation. In one experiment the normal pressure gradients were found to be neligible in regions with large skewing and allowed testing turbulence models using the boundary layer equations. The simulation of this flow compared favorably with the experimental data throughout the flow field and suggested the applicability of algebraic eddy viscosity models for 3DTBLs.

Mueller, U. R.↗

NASA NTRS1981

Edge delamination in angle-ply composite laminates, part 5

A theoretical method was developed for describing the edge delamination stress intensity characteristics in angle-ply composite laminates. The method is based on the theory of anisotropic elasticity. The edge delamination problem is formulated using Lekhnitskii's complex-variable stress potentials and an especially developed eigenfunction expansion method. The method predicts exact orders of the three-dimensional stress singularity in a delamination crack tip region. With the aid of boundary collocation, the method predicts the complete stress and displacement fields in a finite-dimensional, delaminated composite. Fracture mechanics parameters such as the mixed-mode stress intensity factors and associated energy release rates for edge delamination can be calculated explicity. Solutions are obtained for edge delaminated (theta/-theta theta/-theta) angle-ply composites under uniform axial extension. Effects of delamination lengths, fiber orientations, lamination and geometric variables are studied.

Wang, S. S.↗

NASA NTRS1982

Fast generation of three-dimensional computational boundary-conforming periodic grids of C-type

A fast computer program, GRID3C, was developed to generate multilevel three dimensional, C type, periodic, boundary conforming grids for the calculation of realistic turbomachinery and propeller flow fields. The technique is based on two analytic functions that conformally map a cascade of semi-infinite slits to a cascade of doubly infinite strips on different Riemann sheets. Up to four consecutively refined three dimensional grids are automatically generated and permanently stored on four different computer tapes. Grid nonorthogonality is introduced by a separate coordinate shearing and stretching performed in each of three coordinate directions. The grids are easily clustered closer to the blade surface, the trailing and leading edges and the hub or shroud regions by changing appropriate input parameters. Hub and duct (or outer free boundary) have different axisymmetric shapes. A vortex sheet of arbitrary thickness emanating smoothly from the blade trailing edge is generated automatically by GRID3C. Blade cross sectional shape, chord length, twist angle, sweep angle, and dihedral angle can vary in an arbitrary smooth fashion in the spanwise direction.

Dulikravich, D. S.↗

NASA NTRS1982

An implicit solution of the three-dimensional Navier-Stokes equations for an airfoil spanning a wind tunnel

An implicit finite-difference algorithm is developed for the numerical solution of the incompressible three dimensional Navier-Stokes equations in the non-conservative primitive-variable formulation. The flow field about an airfoil spanning a wind-tunnel is computed. The coordinate system is generated by an extension of the two dimensional body-fitted coordinate generation techniques of Thompson, as well as that of Sorenson, into three dimensions. Two dimensional grids are stacked along a spanwise coordinate defined by a simple analytical function. A Poisson pressure equation for advancing the pressure in time is arrived at by performing a divergence operation on the momentum equations. The pressure at each time-step is calculated on the assumption that continuity be unconditionally satisfied. An eddy viscosity coefficient, computed according to the algebraic turbulence formulation of Baldwin and Lomax, simulates the effects of turbulence.

Moitra, A.↗

NASA NTRS1982

Ion current in a magnetic neutral region - Generation of an incipient magnetopause

The current contributed by ions trapped in the vicinity of a magnetic X line is calculated. The three dimensional configuration of the neutral region is found to be critical in determining the current in that the escape mechanism and trapping times depend on the three-dimensional aspects. A trapping criterion is defined. In the neutral region the ions can gain substantial kinetic energy, and the current will change the X line configuration in such a way that there will be a positive feedback effect, rapidly forming an extended magnetopauselike structure for even very small incident plasma densities.

Whipple, E. C.↗

NASA NTRS1982

Plasma and magnetic field characteristics of magnetic flux transfer events

Plasma and magnetic field data from ISEE 1 and 2 are examined for 5 passes of the magnetopause region at 20 and 40 deg northern latitudes, and are presented in terms of moments of the distribution function, calculated from two-dimensional or three-dimensional data. Flux transfer events are characterized by a mixture of magnetosheath and magnetospheric particles, which supports the hypothesis that flux transfer events represent encounters of reconnected flux tubes. An excess pressure appears to be balanced by the tension of the ambient magnetic field lines as they are draped around the reconnected flux tube, and the different observed magnetic field signatures are consistent with expectations for encounters of the flux tubes at different relative locations. It is suggested that increased flow speeds are caused by continued reconnection at the low-latitude boundaries of the flux tubes.

Paschmann, G.↗

NASA NTRS1982

A conservative finite difference algorithm for the unsteady transonic potential equation in generalized coordinates

An implicit, approximate-factorization, finite-difference algorithm has been developed for the computation of unsteady, inviscid transonic flows in two and three dimensions. The computer program solves the full-potential equation in generalized coordinates in conservation-law form in order to properly capture shock-wave position and speed. A body-fitted coordinate system is employed for the simple and accurate treatment of boundary conditions on the body surface. The time-accurate algorithm is modified to a conventional ADI relaxation scheme for steady-state computations. Results from two- and three-dimensional steady and two-dimensional unsteady calculations are compared with existing methods.

Bridgeman, J. O.↗

NASA NTRS1982

Comparison of two- and three-dimensional flow computations with laser anemometer measurements in a transonic compressor rotor

Two and three dimensional inviscid solutions for the flow in a transonic axial compressor rotor at design speed are compared with probe and laser anemometers measurements at near-stall and maximum-flow operating points. Experimental details of the laser anemometer system and computational details of the two dimensional axisymmetric code and three dimensional Euler code are described. Comparisons are made between relative Mach number and flow angle contours, shock location, and shock strength. A procedure for using an efficient axisymmetric code to generate downstream pressure input for computationally expensive Euler codes is discussed. A film supplement shows the calculations of the two operating points with the time-marching Euler code.

Chima, R. V.↗

NASA NTRS1982

Three-dimensional stability of vortex arrays

The stability to three-dimensional disturbances of three classical steady vortex configurations in an incompressible inviscid fluid is studied in the limit of small vortex cross-sectional area and long axial disturbance wavelength. The configurations examined are the single infinite vortex row, the Karman vortex street of staggered vortices and the symmetric vortex street. It is shown that the single row is most unstable to a two-dimensional disturbance, while the Karman vortex street is most unstable to a three-dimensional disturbance over a significant range of street spacing ratios. The symmetric vortex street is found to be most unstable to three-dimensional or two-dimensional symmetric disturbances depending on the spacing ratio of the street. Short remarks are made concerning the relevance of the calculations to the observed instabilities in free shear layer, wake and boundary-layer type flows.

Robinson, A. C.↗

NASA NTRS1981

Viscous-inviscid interaction method including wake effects for three-dimensional wing-body configurations

A viscous-inviscid interaction method has been developed by using a three-dimensional integral boundary-layer method which produces results in good agreement with a finite-difference method in a fraction of the computer time. The integral method is stable and robust and incorporates a model for computation in a small region of streamwise separation. A locally two-dimensional wake model, accounting for thickness and curvature effects, is also included in the interaction procedure. Computation time spent in converging an interacted result is, many times, only slightly greater than that required to converge an inviscid calculation. Results are shown from the interaction method, run at experimental angle of attack, Reynolds number, and Mach number, on a wing-body test case for which viscous effects are large. Agreement with experiment is good; in particular, the present wake model improves prediction of the spanwise lift distribution and lower surface cove pressure.

Streett, C. L.↗

NASA NTRS1982

COSAL: A black-box compressible stability analysis code for transition prediction in three-dimensional boundary layers

A fast computer code COSAL for transition prediction in three dimensional boundary layers using compressible stability analysis is described. The compressible stability eigenvalue problem is solved using a finite difference method, and the code is a black box in the sense that no guess of the eigenvalue is required from the user. Several optimization procedures were incorporated into COSAL to calculate integrated growth rates (N factor) for transition correlation for swept and tapered laminar flow control wings using the well known e to the Nth power method. A user's guide to the program is provided.

Malik, M. R.↗