Engineering Papers⌕ Search

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.

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

At least 469 records · Page 26

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

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

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

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

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

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

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

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

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

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

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

Validation of a three-dimensional viscous analysis of axisymmetric supersonic inlet flow fields

A three-dimensional viscous marching analysis for supersonic inlets was developed. To verify this analysis several benchmark axisymmetric test configurations were studied and are compared to experimental data. Detailed two-dimensional results for shock-boundary layer interactions are presented for flows with and without boundary layer bleed. Three dimensional calculations of a cone at angle of attack and a full inlet at attack are also discussed and evaluated. Results of the calculations demonstrate the code's ability to predict complex flow fields and establish guidelines for future calculations using similar codes.

Benson, T. J.↗

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

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

A time-split finite-volume algorithm for three-dimensional flow-field simulation

A general finite-volume algorithm is developed for solving three-dimensional, time-dependent, compressible Navier-Stokes equations for high Reynolds number flows over an arbitrary geometry. This algorithm adapts MacCormack's (1982) explicit-implicit scheme to a time-split, three-dimensional finite-volume concept in a general coordinate system. It is shown that the thin-layer approximation in all three spatial directions significantly reduces the evaluation of viscous terms and allows the algorithm to solve more complicated geometries with all boundaries in two or all three directions. The calculated results using this method are found to be in good agreement with the experimental measurements of a blunt-fin induced shock wave and boundary-layer interaction problems. Observations of the existence of peak pressure, primary horseshoe and secondary vortices, and reversed supersonic zones show that computational fluid dynamics can effectively supplement the wind tunnel tests for aerodynamic design as well as for understanding basic fluid dynamics.

Hung, C. M.↗

An aeroacoustic model for high-speed, unsteady blade-vortex interaction

A three-dimensional aeroacoustic model is developed to predict the sound pulse radiated by the passage of a helicopter blade over a potential vortex. The linearized analysis assumes that either the blade-vortex separation is small or that the blade-tip Mach number is close to 1, or both, so that an acoustically noncompact situation exists. The three-dimensional blade loading due to blade-vortex interaction is constructed through a spanwise superposition of two-dimensional solutions with strength linearly increasing from hub to tip. Such a loading overestimates somewhat the strength of tip region dipoles in the acoustic calculation that follows. The final expression for the predicted far-field signature is obtained in closed form, and thus permits a relatively inexpensive calculation of the directivity of peak acoustic pressures in three dimensions.

Martinez, R.↗

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

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