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 613 records · Page 34

Statistical mechanics of light elements at high pressure. V Three-dimensional Thomas-Fermi-Dirac theory

A numerical technique for solving the Thomas-Fermi-Dirac (TED) equation in three dimensions, for an array of ions obeying periodic boundary conditions, is presented. The technique is then used to calculate deviations from ideal mixing for an alloy of hydrogen and helium at zero temperature and high presures. Results are compared with alternative models which apply perturbation theory to calculation of the electron distribution, based upon the assumption of weak response of the electron gas to the ions. The TFD theory, which permits strong electron response, always predicts smaller deviations from ideal mixing than would be predicted by perturbation theory. The results indicate that predicted phase separation curves for hydrogen-helium alloys under conditions prevailing in the metallic zones of Jupiter and Saturn are very model dependent.

Macfarlane, J. J.↗

Three-dimensional effects on airfoils

The effects of boundary layer flows along the walls of wind tunnels were studied to validate the transfer of two dimensional calculations to three dimensional transonic flowfield calculations. Results from trials in various wind tunnels were examind to determine the effects of the wall boundary flow on the control surfaces of an airfoil. Models sliding along a groove in the wall of a channel at sub- and transonic speeds were examined, with the finding that with either nonuniformities in the groove, or even if the channel walls are uniform, the lateral boundary layer can cause variations in the central flow region or alter the onset of shock at the transition point. Models for the effects in both turbulence and in the absence of turbulence are formulated, and it is noted that the characteristics of individual wind tunnels must be studied to quantify any existing three dimensional effects.

Chevallier, J. P.↗

Advanced numerics for multi-dimensional fluid flow calculations

In recent years, there has been a growing interest in the development and use of mathematical models for the simulation of fluid flow, heat transfer and combustion processes in engineering equipment. The equations representing the multi-dimensional transport of mass, momenta and species are numerically solved by finite-difference or finite-element techniques. However despite the multiude of differencing schemes and solution algorithms, and the advancement of computing power, the calculation of multi-dimensional flows, especially three-dimensional flows, remains a mammoth task. The following discussion is concerned with the author's recent work on the construction of accurate discretization schemes for the partial derivatives, and the efficient solution of the set of nonlinear algebraic equations resulting after discretization. The present work has been jointly supported by the Ramjet Engine Division of the Wright Patterson Air Force Base, Ohio, and the NASA Lewis Research Center.

Vanka, S. P.↗

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

The Goertler vortex instability mechanism in three-dimensional boundary layers

The two dimensional boundary layer on a concave wall is centrifugally unstable with respect to vortices aligned with the basic flow for sufficiently high values of the Goertler number. However, in most situations of practical interest the basic flow is three dimensional and previous theoretical investigations do not apply. The linear stability of the flow over an infinitely long swept wall of variable curvature is considered. If there is no pressure gradient in the boundary layer the instability problem can always be related to an equivalent two dimensional calculation. However, in general, this is not the case and even for small values of the crossflow velocity field dramatic differences between the two and three dimensional problems emerge. When the size of the crossflow is further increased, the vortices in the neutral location have their axes locally perpendicular to the vortex lines of the basic flow.

Hall, P.↗

Investigation of the three-dimensional flow field within a transonic fan rotor: Experiment and analysis

An experimental investigation of the three dimensional flow field through a low aspect ratio, transonic, axial flow fan rotor has been conducted using an advanced laser anemometer (LA) system. Laser velocimeter measurements of the rotor flow field at the design operating speed and over a range of through flow conditions are compared to analytical solutions. The numerical technique used herein yields the solution to the full, three dimensional, unsteady Euler equations using an explicit time marching, finite volume approach. The numerical analysis, when coupled with a simplified boundary layer calculation, generally yields good agreement with the experimental data. The test rotor has an aspect ratio of 1.56, a design total pressure ratio of 1.629 and a tip relative Mach number of 1.38. The high spatial resolution of the LA data matrix (9 radial by 30 axial by 50 blade to blade) permits details of the transonic flow field such as shock location, turning distribution and blade loading levels to be investigated and compared to analytical results.

Pierzga, M. J.↗

Atmospheric transfer of radiation above an inhomogeneous non-Lambertian reflective ground. I - Theory

A method for solving the three-dimensional equation of transfer for a vertically inhomogeneous atmosphere bounded by a reflecting surface of non-uniform bidirectional reflectance is presented. The technique incorporates a two-dimensional spatial Fourier transform of the transfer equation and solution of the resulting expressions for each Fourier component of the radiation field using the method of Gauss-Seidel iteration. The intensity field in the spatial domain, which is calculated for a variety of altitudes, zenith angles, and azimuths, is reconstructed using the inverse Fourier transform. An empirical surface bidirectional reflectance function is employed, permitting the consideration of non-Lambertian reflective properties.

Diner, D. J.↗

Numerical simulation of controlled flow tunnel for V/STOL testing

A 'smart' wind tunnel concept for V/STOL transition regime studies is presented and performance is projected with numerical simulations. The wind tunnel would provide a flow that faithfully represents free air flow and thereby permits accurate lift calculations for V/STOL aircraft. The design concept is based on free-air tunnels used for planar wings. Modifications are necessary to account for the highly deflected, vortical wake produced in low-speed take-off and landing situations. A three-step numerical procedure covering the potential flow of a high lift system in free-air, controlled flow data with a model in a wind tunnel and addition of the wall conditions to the simulation as a feedback loop to correct the closed-tunnel is described. The initial simulation is for three-dimensional full span jet flapped wings, and is derived from extended two-dimensional calculations. The results indicate that only partial control of the tunnel flow is sufficient to serve as a valid environment for V/STOl studies.

Joppa, R. G.↗

On the Global Distribution of Three-Dimensional Eliassen-Palm Fluxes by Stationary Waves

The global distribution of three dimensional EP fluxes is investigated by using statistics for Dec. 1980 to Feb. 1981 and June to August 1981 calculated from operational analyses. Results at 1000 and 150 mb are given. During the Northern Hemisphere (NH) winter strong upward fluxes can be seen at 1000 mb in the lee of the Himalayas and Rockies, windward of the Canadian Rockies in the east Pacific, to the north of the Atlantic stormtrack and over the west Pacific in a region of strong land sea thermal contrast. Upward fluxes also appear near the west coasts of the Southern Hemisphere (SH) subtropical continents regions of strong land sea thermal contrast. At 150 mb (Fig. 1b) the strongest upward flux now occurs south of Alaska. The horizontal fluxes imply poleward propagation from the convection over Indonesia to 20 N, a region where the assumptions underlying EP fluxes are not well met. Equatorward propagation dominates the NH midlatitude except over the Pacific. The fluxes given offer hints of the sources of stationary waves, but also show several puzzling features and a rather cavalier disregard of regions of easterly wind. The physical meaning and interpretation of three dimensional EP fluxes is not yet clear.

White, G. H.↗

An implicit finite-difference code for a two-equation turbulence model for three-dimensional flows

An implicit finite difference code was developed which solves the transport equations for the turbulence kinetic energy and its dissipation rate in generalized coordinates in three dimensions. The finite difference equations are solved using the Beam-Warming algorithm. The kinetic energy-dissipation code, KEM, provides the closure; i.e., the turbulent viscosity for calculation of either compressible or incompressible flows. Turbulent internal flow over a backward-facing step has been calculated using the present code in conjunction with the Incompressible Navier-Stokes Code, INS3D. The results are in good agreement with experiments and two dimensional computations of other researchers.

Kaul, U. K.↗

Finite difference methods for reducing numerical diffusion in TEACH-type calculations

A methodological evaluation for two-finite differencing schemes for computer-aided gas turbine design is presented. The two computational schemes include; a Bounded Skewed Finite Differencing Scheme (BSUDS); and a Quadratic Upwind Differencing Scheme (QSDS). In the evaluation, the derivations of the schemes were incorporated into two-dimensional and three-dimensional versions of the Teaching Axisymmetric Characteristics Heuristically (TEACH) computer code. Assessments were made according to performance criteria for the solution of problems of turbulent, laminar, and coannular turbulent flow. The specific performance criteria used in the evaluation were simplicity, accuracy, and computational economy. It is found that the BSUDS scheme performed better with respect to the criteria than the QUDS. Some of the reasons for the more successful performance BSUDS are discussed.

Syed, S. A.↗

Potential flow calculations and preliminary wing design in support of an NLF variable sweep transition flight experiment

NASA Langley and NASA Ames-Dryden have defined a variable-sweep transition-flight experiment utilizing the F-14 aircraft to enhance understanding of the interaction of crossflow and Tollmien-Schlichting instabilities on a laminar-boundary-layer transition. The F-14 wing outer panel will be modified to generate favorable pressure gradients on the upper wing surface over a wide range of flight conditions. Extensive computations have been performed using two-dimensional and three-dimensional transonic analysis codes. Flight-test and computational data are compared and shown to validate the applicability of the three-dimensional codes (WBPPW and TAWFIVE). In addition, results from two preliminary glove designs derived from two different approaches to the design problem are presented. Advantages and disadvantages of each approach are identified, and it is concluded that coupling an analysis code with an automated design procedure yields a powerful code with distinct advantages over a 'cut-and-dry' approach.

Waggoner, E. G.↗

Shock-free turbomachinery blade design

A computational method for designing shock-free, quasi-three-dimensional, transonic, turbomachinery blades is described. Shock-free designs are found by implementing Sobieczky's fictitious gas principle in the analysis of a baseline shape, resulting in an elliptic solution that is incorrect in the supersonic domain. Shock-free designs are obtained by combining the subsonic portion of this solution with a characteristic calculation of the correct supersonic flow using the sonic line data from the fictitious elliptic solution. This provides a new, shock-free blade design. Examples presented include the removal of shocks from two blades in quasi-three-dimensional flow and the development of a series of shock-free two-dimensional stators. The new designs all include modifications to the upper surface of an experimental stator blade developed at NASA Lewis Research Center. While the designs presented here are for inviscid flow, the same concepts have been successfully applied to the shock-free design of airfoils and three-dimensional wings with viscous effects. The extension of the present method to viscous flows is straightforward given a suitable analysis algorithm for the flow.

Beauchamp, P. P.↗

Line emission from H II blister models

Numerical techniques to calculate the thermal and geometric properties of line emission from H II 'blister' regions are presented. It is assumed that the density distributions of the H II regions are a function of two dimensions, with rotational symmetry specifying the shape in three-dimensions. The thermal and ionization equilibrium equations of the problem are solved by spherical modeling, and a spherical sector approximation is used to simplify the three-dimensional treatment of diffuse ionizing radiation. The global properties of H II 'blister' regions near the edges of a molecular cloud are simulated by means of the geometry/density distribution, and the results are compared with observational data. It is shown that there is a monotonic increase of peak surface brightness from the i = 0 deg (pole-on) observational position to the i = 90 deg (edge-on) position. The enhancement of the line peak intensity from the edge-on to the pole-on positions is found to depend on the density, stratification, ionization, and electron temperature weighting. It is found that as i increases, the position of peak line brightness of the lower excitation species is displaced to the high-density side of the high excitation species.

Rubin, R. H.↗

Motion analysis report

Human motion analysis is the task of converting actual human movements into computer readable data. Such movement information may be obtained though active or passive sensing methods. Active methods include physical measuring devices such as goniometers on joints of the body, force plates, and manually operated sensors such as a Cybex dynamometer. Passive sensing de-couples the position measuring device from actual human contact. Passive sensors include Selspot scanning systems (since there is no mechanical connection between the subject's attached LEDs and the infrared sensing cameras), sonic (spark-based) three-dimensional digitizers, Polhemus six-dimensional tracking systems, and image processing systems based on multiple views and photogrammetric calculations.

Badler, N. I.↗

Simulation and modeling of homogeneous, compressed turbulence

Low Reynolds number homogeneous turbulence undergoing low Mach number isotropic and one-dimensional compression was simulated by numerically solving the Navier-Stokes equations. The numerical simulations were performed on a CYBER 205 computer using a 64 x 64 x 64 mesh. A spectral method was used for spatial differencing and the second-order Runge-Kutta method for time advancement. A variety of statistical information was extracted from the computed flow fields. These include three-dimensional energy and dissipation spectra, two-point velocity correlations, one-dimensional energy spectra, turbulent kinetic energy and its dissipation rate, integral length scales, Taylor microscales, and Kolmogorov length scale. Results from the simulated flow fields were used to test one-point closure, two-equation models. A new one-point-closure, three-equation turbulence model which accounts for the effect of compression is proposed. The new model accurately calculates four types of flows (isotropic decay, isotropic compression, one-dimensional compression, and axisymmetric expansion flows) for a wide range of strain rates.

Wu, C. T.↗

Assessment of a 3-D boundary layer code to predict heat transfer and flow losses in a turbine

Zonal concepts are utilized to delineate regions of application of three-dimensional boundary layer (DBL) theory. The zonal approach requires three distinct analyses. A modified version of the 3-DBL code named TABLET is used to analyze the boundary layer flow. This modified code solves the finite difference form of the compressible 3-DBL equations in a nonorthogonal surface coordinate system which includes coriolis forces produced by coordinate rotation. These equations are solved using an efficient, implicit, fully coupled finite difference procedure. The nonorthogonal surface coordinate system is calculated using a general analysis based on the transfinite mapping of Gordon which is valid for any arbitrary surface. Experimental data is used to determine the boundary layer edge conditions. The boundary layer edge conditions are determined by integrating the boundary layer edge equations, which are the Euler equations at the edge of the boundary layer, using the known experimental wall pressure distribution. Starting solutions along the inflow boundaries are estimated by solving the appropriate limiting form of the 3-DBL equations.

Anderson, O. L.↗