SEARCH · Engineering Papers

Results for “numerics”

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,171 records · Page 65

All sources NASA NTRS DOE OSTI

NASA NTRS1993

Euler solutions for blunt bodies using triangular meshes - Artificial viscosity forms and numerical boundary conditions

A finite volume method is used to calculate compressible inviscid flows over blunt bodies using, in general, unstructured grids. Artificial viscosity forms are derived based on a simplified least squares procedure. The extra second order terms are consistent with the governing equations, hence a systematic treatment of the numerical boundary conditions can be easily implemented. A special treatment of blunt bodies may be required. The discrete equations are linearized and the resulting system is solved by a relaxation method. Preliminary results indicate that the effect of the numerical dissipation is minimal. For subsonic flows over smooth bodies, the solution is practically vorticity-free and the total pressure loss is of the same order as the truncation error. Finally, some extensions of the present method are briefly discussed.

Winterstein, R.↗

NASA NTRS1992

Direct numerical simulations of a reacting turbulent mixing layer by a pseudospectral-spectral element method

An account is given of the implementation of the spectral-element technique for simulating a chemically reacting, spatially developing turbulent mixing layer. Attention is given to experimental and numerical studies that have investigated the development, evolution, and mixing characteristics of shear flows. A mathematical formulation is presented of the physical configuration of the spatially developing reacting mixing layer, in conjunction with a detailed representation of the spectral-element method's application to the numerical simulation of mixing layers. Results from 2D and 3D calculations of chemically reacting mixing layers are given.

Mcmurtry, Patrick A.↗

NASA NTRS1993

Numerical simulation of laminar-turbulent transition in a spatially-developing flat plate wake

Laminar-turbulent transition of an incompressible flat-plate wake is investigated by direct numerical integration of the Navier-Stokes equations. For the numerical integration, a combination of finite-difference and spectral methods along with an ADI/Crank-Nicolson/Adams-Bashforth time integration scheme is employed. Subject to 2D forcing, the wake exhibited a rapidly-growing fundamental disturbance that quickly saturated. This saturation was due partly to the stabilizing effect of the mean flow distortion. Downstream of the saturation point, disturbance energy was concentrated in the fundamental disturbance, the second harmonic, and the mean flow distortion component. At large amplitude levels, a Karman vortex street formed. Variations in the 2D forcing level did not alter the qualitative behavior of the disturbances. Simulations of 3D breakdown indicates that the presence of large-amplitude, 2D disturbances tends to initially suppress small-amplitude 3D disturbance growth. Following this initial suppression, a resumption of 3D growth is observed that may have been due to a secondary instability mechanism. For high levels of 3D disturbance energy, lambda-shaped vortical structures formed between adjacent Karman vortices.

Dratler, D. I.↗

NASA NTRS1993

A numerical study of the effect of geometry variation, turbulence models, and dissipation on the flow past control surfaces

The steady 3D viscous flow past the ONERA M6 wing and a slender delta wing-body with trailing edge control surfaces has been computed. A cell-centered finite-volume Navier-Stokes patched zonal method has been used for the numerical simulation. Both diagonalized and LUSGS schemes have been implemented. Besides the standard nonplanar zonal interfacing techniques, a new virtual zone capability has been employed. For code validation, the transonic flow past the ONERA M5 wing is calculated for angles-of-attack of 3.06 deg and 5.06 deg and compared with the available experiments. The wing-body computational results are compared with experimental data for both trailing-edge flaps deflected. The experimental flow conditions are M subinfinity = 0.4, a turbulent Reynolds number of 5.41 million based on a mean aerodynamic chord of 25.959 inches, adiabatic wall, and angles-of-attack varying from 0 deg to 23.85 deg. The computational results are presented for the 23.85 deg angle-of-attack case. The effects of the base flow due to a model sting, the varying second and fourth order numerical dissipation, and the turbulence model are all considered.

Chaussee, Denny S.↗

NASA NTRS1993

Experimental and numerical investigation of supersonic turbulent flow in an annular duct

Experimental and numerical results are presented for developing supersonic turbulent flow in an annular duct formed by a circular centerbody and outer shroud. The experimental results are based on data taken in a new flow facility that was designed to generate a shock-free, supersonic annular flow. Numerical computations were performed using the Baldwin-Lomax turbulence model for comparison with experimentally measured profiles. The results demonstrate that computed and measured profiles are in excellent agreement, so that studies can now be conducted of shock wave/boundary layer interaction phenomena within the duct, such as those induced by changes in downstream duct geometry or by the placement of struts between the duct walls.

Williams, K. E.↗

NASA NTRS1991

Two-dimensional numerical modeling and solution of convection heat transfer in turbulent He II

Numerical schemes are employed to investigate heat transfer in the turbulent flow of He II. FEM is used to solve a set of equations governing the heat transfer and hydrodynamics of He II in the turbulent regime. Numerical results are compared with available experimental data and interpreted in terms of conventional heat transfer parameters such as the Prandtl number, the Peclet number, and the Nusselt number. Within the prescribed Reynolds number domain, the Gorter-Mellink thermal counterflow mechanism becomes less significant, and He II acts like an ordinary fluid. The convection heat transfer characteristics of He II in the highly turbulent regime can be successfully described by using the conventional turbulence and heat transfer theories.

Zhang, Burt X.↗

NASA NTRS1991

Numerical simulation of the helium gas spin-up channel performance of the relativity gyroscope

The dependence of the spin-up system efficiency on each geometrical parameter of the spin-up channel and the exhaust passage of the Gravity Probe-B (GPB) is individually investigated. The spin-up model is coded into a computer program which simulates the spin-up process. Numerical results reveal optimal combinations of the geometrical parameters for the ultimate spin-up performance. Comparisons are also made between the numerical results and experimental data. The experimental leakage rate can only be reached when the gap between the channel lip and the rotor surface increases beyond physical limit. The computed rotating frequency is roughly twice as high as the measured ones although the spin-up torques fairly match.

Karr, Gerald R.↗

NASA NTRS1993

Numerical modeling of preburner flowfield

This work is intended to numerically predict the flowfields inside the preburner of the Space Shuttle Main Engine. The computer code (FDNS) based on pressure correction method is modified and adapted with an elliptic grid generator. The original configuration of the preburner in conjunction with downstream gas turbines has been simplified geometrically and numerically modeled at its full power in this work. The computational results are presented and qualitatively discussed with test data collected in NASA/MSFC.

Chow, A. S.↗

NASA NTRS1993

Numerical simulation and physical aspects of supersonic vortex breakdown

Existing numerical simulations and physical aspects of subsonic and supersonic vortex-breakdown modes are reviewed. The solution to the problem of supersonic vortex breakdown is emphasized in this paper and carried out with the full Navier-Stokes equations for compressible flows. Numerical simulations of vortex-breakdown modes are presented in bounded and unbounded domains. The effects of different types of downstream-exit boundary conditions are studied and discussed.

Liu, C. H.↗

NASA NTRS1939

Some notes on the numerical solution of shear-lag and mathematically related problems

The analysis of box beams with shear deformation of the flanges can be reduced to the solution of a differential equation. The same equation is met in other problems of stress analysis. No analytical solutions of this equation can be given for practical cases, and numerical methods of evaluation must be used. Available methods are briefly discussed. Two numerical examples show the application of the step-by-step method of integration to shear-lag problems.

Kuhn, Paul↗

NASA NTRS1952

An iterative transformation procedure for numerical solution of flutter and similar characteristics-value problems

An iterative transformation procedure suggested by H. Wielandt for numerical solution of flutter and similar characteristic-value problems is presented. Application of this procedure to ordinary natural-vibration problems and to flutter problems is shown by numerical examples. Comparisons of computed results with experimental values and with results obtained by other methods of analysis are made.

Gossard, Myron L↗

NASA NTRS1993

Numerical simulation of the non-Newtonian mixing layer

This work is a continuing effort to advance our understanding of the effects of polymer additives on the structures of the mixing layer. In anticipation of full nonlinear simulations of the non-Newtonian mixing layer, we examined in a first stage the linear stability of the non-Newtonian mixing layer. The results of this study show that, for a fluid described by the Oldroyd-B model, viscoelasticity reduces the instability of the inviscid mixing layer in a special limit where the ratio (We/Re) is of order 1 where We is the Weissenberg number, a measure of the elasticity of the flow, and Re is the Reynolds number. In the present study, we pursue this project with numerical simulations of the non-Newtonian mixing layer. Our primary objective is to determine the effects of viscoelasticity on the roll-up structure. We also examine the origin of the numerical instabilities usually encountered in the simulations of non-Newtonian fluids.

Azaiez, Jalel↗

NASA NTRS1993

Surrogates for numerical simulations; optimization of eddy-promoter heat exchangers

Although the advent of fast and inexpensive parallel computers has rendered numerous previously intractable calculations feasible, many numerical simulations remain too resource-intensive to be directly inserted in engineering optimization efforts. An attractive alternative to direct insertion considers models for computational systems: the expensive simulation is evoked only to construct and validate a simplified, input-output model; this simplified input-output model then serves as a simulation surrogate in subsequent engineering optimization studies. A simple 'Bayesian-validated' statistical framework for the construction, validation, and purposive application of static computer simulation surrogates is presented. As an example, dissipation-transport optimization of laminar-flow eddy-promoter heat exchangers are considered: parallel spectral element Navier-Stokes calculations serve to construct and validate surrogates for the flowrate and Nusselt number; these surrogates then represent the originating Navier-Stokes equations in the ensuing design process.

Patera, Anthony T.↗

NASA NTRS1993

Computer code for scattering from impedance bodies of revolution. Part 3: Surface impedance with s and phi variation. Analytical and numerical results

The third phase of the development of the computer codes for scattering by coated bodies that has been part of an ongoing effort in the Electromagnetics Laboratory of the Electrical Engineering and Computer Science Department at the University of Illinois at Chicago is described. The work reported discusses the analytical and numerical results for the scattering of an obliquely incident plane wave by impedance bodies of revolution with phi variation of the surface impedance. Integral equation formulation of the problem is considered. All three types of integral equations, electric field, magnetic field, and combined field, are considered. These equations are solved numerically via the method of moments with parametric elements. Both TE and TM polarization of the incident plane wave are considered. The surface impedance is allowed to vary along both the profile of the scatterer and in the phi direction. Computer code developed for this purpose determines the electric surface current as well as the bistatic radar cross section. The results obtained with this code were validated by comparing the results with available results for specific scatterers such as the perfectly conducting sphere. Results for the cone-sphere and cone-cylinder-sphere for the case of an axially incident plane were validated by comparing the results with the results with those obtained in the first phase of this project. Results for body of revolution scatterers with an abrupt change in the surface impedance along both the profile of the scatterer and the phi direction are presented.

Uslenghi, Piergiorgio L. E.↗

NASA NTRS1993

Numerical Studies of Collisionless Current Layers

The purpose of this proposal was to investigate collisionless current layers using a variety of analytic and numerical tools. The first year of the contract was dedicated to analytical studies, to the porting and adaption of codes being used in this study, and to the numerical simulation of collisionless current layers. The second year entailed the development of multi-dimensional hybrid algorithms as well as the re-examination of the problem of integro-differential equations that occur in the linear stage of plasma instabilities.

Quest, Kevin B.↗

NASA NTRS1993

Numerical simulation of three-dimensional self-gravitating flow

The three-dimensional flow of a self-gravitating fluid is numerically simulated using a Fourier pseudospectral method with a logarithmic variable formulation. Two cases with zero total angular momentum are studied in detail, a 323 simulation (Run B). Other than the grid size, the primary difference between the two cases are that Run A modeled atomic hydrogen and had considerably more compressible motion initially than Run B, which modeled molecular hydrogen. The numerical results indicate that gravitational collapse can proceed in a variety of ways. In the Run A, collapse led to an elongated tube-like structure, while in the Run B, collapse led to a flatter, disklike structure.

Shebalin, John V.↗

NASA NTRS1993

Multigrid direct numerical simulation of the whole process of flow transition in 3-D boundary layers

A new technology was developed in this study which provides a successful numerical simulation of the whole process of flow transition in 3-D boundary layers, including linear growth, secondary instability, breakdown, and transition at relatively low CPU cost. Most other spatial numerical simulations require high CPU cost and blow up at the stage of flow breakdown. A fourth-order finite difference scheme on stretched and staggered grids, a fully implicit time marching technique, a semi-coarsening multigrid based on the so-called approximate line-box relaxation, and a buffer domain for the outflow boundary conditions were all used for high-order accuracy, good stability, and fast convergence. A new fine-coarse-fine grid mapping technique was developed to keep the code running after the laminar flow breaks down. The computational results are in good agreement with linear stability theory, secondary instability theory, and some experiments. The cost for a typical case with 162 x 34 x 34 grid is around 2 CRAY-YMP CPU hours for 10 T-S periods.

Liu, Chaoqun↗

NASA NTRS1993

Numerical investigation of complex, transitional, and chemically reacting flows

This document is divided into two parts. Part A, on turbulent transport of chemical species in compressible reacting flows, includes the following topics: mixing of chemical species in reacting flows; chemically reacting turbulent flows; turbulent mixing of chemical species in reacting flows, which includes large eddy mixing, turbulent mixing, and molecular mixing; multiple time scale turbulence equations; species conservation equations for reacting flows; and a finite volume numerical method. Part B, on unsteady transitional flows over forced oscillatory surfaces, discusses unsteady transitional flow equations with moving boundaries and the application of the finite volume numerical solution method.

Kim, S.-W.↗