Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Numerical Simulation”

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 1,477 records · Page 82
All sourcesNASA NTRSDOE OSTI
NASA NTRS1987

Numerical simulations of trailing vortex bursting

Solutions of the steady-state Navier-Stokes equations for the axisymmetric bursting of a laminar trailing vortex are computed with Newton's method and the pseudo-arc length continuation method for wide ranges of vortex strength and Reynolds number. The results indicate that a trailing vortex can undergo a transition from a state in which the core slowly diffuses to a state marked by large amplitude, spatial oscillations of core radius and core axial velocity. At the transition point the core grows rapidly in size. This event is interpreted as vortex bursting. The results also suggest that when the maximum core swirl velocity is sufficiently large the centerline axial flow downstream of transition will be reversed.

Beran, Philip S.↗
NASA NTRS1987

A numerical simulation of finite-length Taylor-Couette flow

The processes leading to laminar-turbulent transition in finite-channel-length Taylor-Couette flow are investigated analytically, solving the unsteady incompressible Navier-Stokes equations by spectral-collocation methods. A time-split algorithm, implementable in both axisymmetric and fully three-dimensional time-accurate versions, and an algorithm based on the staggered-mesh discretization of Bernardi and Maday (1986) are described in detail, and results obtained by applying the axisymmetric version of the first algorithm and a steady-state version of the second are presented graphically and compared with published experimental data. The feasibility of full three-dimensional simulations of the progression through chaotic states to turbulence under the constraints of Taylor-Couette flow is demonstrated.

Streett, C. L.↗
NASA NTRS1987

Numerical simulation of flow path in the oxidizer side hot gas manifold of the Space Shuttle main engine

The purpose of this study is to examine in detail incompressible laminar and turbulent flows inside the oxidizer side Hot Gas Manifold of the Space Shuttle Main Engine. To perform this study, an implicit finite difference code cast in general curvilinear coordinates is further developed. The code is based on the method of pseudo-compressibility and utilize ADI or implicit approximate factorization algorithm to achieve computational efficiency. A multiple-zone method is developed to overcome the complexity of the geometry. In the present study, the laminar and turbulent flows in the oxidizer side Hot Gas Manifold have been computed. The study reveals that: (1) there exists large recirculation zones inside the bowl if no vanes are present; (2) strong secondary flows are observed in the transfer tube; and (3) properly shaped and positioned guide vanes are effective in eliminating flow separation.

Lin, S. J.↗
NASA NTRS1987

Numerical simulation of the flow field and fuel sprays in an IC engine

A two-dimensional model for axisymmetric piston-cylinder configurations is developed to study the flow field in two-stroke direct-injection Diesel engines under motored conditions. The model accounts for turbulence by a two-equation model for the turbulence kinetic energy and its rate of dissipation. A discrete droplet model is used to simulate the fuel spray, and the effects of the gas phase turbulence on the droplets is considered. It is shown that a fluctuating velocity can be added to the mean droplet velocity every time step if the step is small enough. Good agreement with experimental data is found for a range of ambient pressures in Diesel engine-type microenvironments. The effects of the intake swirl angle in the spray penetration, vaporization, and mixing in a uniflow-scavenged two-stroke Diesel engine are analyzed. It is found that the swirl increases the gas phase turbulence levels and the rates of vaporization.

Nguyen, H. L.↗
NASA NTRS1987

Numerical simulations of unsteady, viscous, transonic flow over isolated and cascaded airfoils using a deforming grid

A compressible, unsteady, full Navier-Stokes, finite difference code was developed for modeling transonic flow through two-dimensional, oscillating cascades. The procedure introduces a deforming grid technique to capture the motion of the airfoils. Results using a deforming grid are presented for both isolated and cascaded airfoils. The load histories and unsteady pressure distributions are predicted for the NASA 64A010 isolated airfoil and compared with existing experimental data. Results show that the deforming grid technique can be used to successfully predict the unsteady flow properties around an oscillating airfoil. The deforming grid technique was extended for modeling unsteady flow in a cascade. The use of a deforming grid simplifies the specification of boundary conditions. Unsteady flow solutions similar to the isolated airfoil predictions are found for a NACA 0012 cascade with zero interblade phase angle and zero stagger. Experimental data for these cases are not available for code validation, but computational results are presented to show sample predictions from the code. Applications of the code to typical turbomachinery flow conditions will be presented in future work.

Huff, Dennis L.↗
NASA NTRS1987

Numerical simulation of the generation of turbulence from cometary ion pick-up

Observations of magnetic field fluctuations near Comet Halley have revealed a rapid development of a Kolmogoroff-like turbulence spectrum extending from below 0.01 Hz to above 0.1 Hz. Spectra obtained far from the comet have a strong peak in power near the Doppler-shifted ion-cyclotron frequency of singly ionized water. Closer to the comet, the spectrum at higher frequencies is enhanced in power level over the background solar wind spectrum by approximately an order of magnitude. The equations of incompressible MHD are solved using a two-dimensional 256 x 256 mode spectral method code to simulate this spectral evolution as an inertial range turbulent cascade. The initial conditions contained a constant magnetic field and a single coherent wave mode at a low wave number. The solar wind turbulence was modeled by a background noise spectrum having a Kolmogoroff spectral index. The coherent mode decayed into an inertial range spectrum with Kolmogoroff slope within a few eddy-turnover times. Both the time scale and the increase in power level of the turbulence seen in the simulation are in accord with the Giotto observations.

Goldstein, M. L.↗
NASA NTRS1986

Theory of electrophoretic separations. II - Construction of a numerical simulation scheme and its applications

The implementation of a previously derived model describing electrophoretic transport processes is reformulated, and its usefulness is demonstrated by simulating moving boundary electrophoresis and isoelectric focusing with immobilized species. The resulting solution is stable, evolves smoothly when suitable small step sizes are used, and has a loosely defined convergence as the step sizes are made smaller. The major computational problem involves a need to use small spatial step sizes at high current levels.

Palusinski, O. A.↗
NASA NTRS1987

Collisional interchange instability. II - A comparison of the numerical simulations with the in situ experimental data

A two-dimensional model applicable to F region ionosphere plasma instabilities has been developed and described in a companion paper. Here the model is applied to equatorial F region irregularities and is tested against rocket and satellite data. As a diagnostic, simulated data sets are created which are similar to the one-dimensional measurements of plasma density performed by space probes and the Fourier transform of these data is taken in the same manner and used by the space experimentors. Unlike previous simulations of this phenomenon, an inherent anisotropy is found in the instability development which is mirrored in the in situ data. Evidence that the shallow spectral slopes which often characterize spread F rocket spectra near the F peak may be due to a change in the angle between the rocket velocity vector and the characteristic directions in the medium is presented.

Kelley, M. C.↗
NASA NTRS1986

Numerical simulation of scramjet inlet flow fields

A computer program was developed to analyze supersonic combustion ramjet (scramjet) inlet flow fields. The program solves the three-dimensional Euler or Reynolds averaged Navier-Stokes equations in full conservation form by either the fully explicit or explicit-implicit, predictor-corrector method of MacCormack. Turbulence is modeled by an algebraic eddy-viscosity model. The analysis allows inclusion of end effects which can significantly affect the inlet flow field. Detailed laminar and turbulent flow results are presented for a symmetric-wedge corner, and comparisons are made with the available experimental results to allow assessment of the program. Results are then presented for two inlet configurations for which experimental results exist at the NASA Langley Research Center.

Kumar, Ajay↗
NASA NTRS1987

Numerical simulation of self-sustained and forced oscillations in jet shear layers

Excitation of axisymmetric jet shear layers has been investigated by solving the time-dependent compressible Navier-Stokes equations. MacCormack's explicit finite difference scheme is used in performing the computations on a CRAY X-MP computer. This work examines the influence of various means of excitation of the shear layer, to produce large scale vortex structures in such a manner that the mixing between two flows will be enhanced or diminished. While successful pressure excitation has been achieved in previous studies only limited success with temperature excitation at the inflow boundary has been accomplished. The current study focuses on the use of temperature excitation through the boundary layer immediately inside the jet lip. In order to achieve the required resolution significant grid refinement is necessary. The results of the computation with and without temperature excitation are analyzed and compared with results obtained using pressure excitation. The influence of amplitude and frequency on specific flow characteristics are also examined. The excitation frequency for both the temperature and pressure corresponds to a Strouhal number of .45 and the inflow jet Mach number is approximately .3.

Scott, J. N.↗
NASA NTRS1987

Advances in magnetospheric physics by numerical simulations - A critical analysis by an outsider

This is a selective, somewhat editorialized assessment, based on a literature survey, of computer simulation as it exists today in magnetospheric physics. Both large (MHD) and small scale dynamical simulations are described and considered from the perspective of what they are trying to do and with what success. Several specific problem areas where simulations are being carried out are called for commentary: global magnetospheric structure, magnetotail and magnetopause reconnection, Kelvin-Helmholtz instability, hydrodynamic expansion of ionospheric ions, electric double layers, heavy ion heating, and auroral kilometric radiation. It is concluded that simulation is a necessary tool for understanding magnetospheric physics and that significant progress has been made in simulation development. However, results should be evaluated knowing that many factors, some real and physical, others structural, may contribute to such output. A plea is made for greater intercalibration among different simulators working in parallel areas, so that facts can be distinguished from artifacts.

Birmingham, Thomas J.↗
NASA NTRS1988

Numerical simulation of wing lift augmentation with spanwise tip blowing

Navier-Stokes simulations are used to study the phenomena of wing lift augmentation with spanwise tip blowing. Flowfields calculated for varying blowing conditions are analyzed to determine the mechanisms responsible for providing lift augmentation. It is found that both viscous and inviscid effects are important.

Caruso, S. C.↗
NASA NTRS1988

Numerical simulation of aircraft rotary aerodynamics

A potential-flow based panel method was used to compute the aerodynamic loads along a three-dimensional flight path. Comparisons were made with rotary rig tests of an airplane performing a coning motion and found to compare well with the main components of the aerodynamic forces. The effect of the sting-type wind tunnel mounting of the model on the measured loads was briefly investigated, and the possible benefits of using this method combined with wind tunnel experiments were highlighted.

Katz, Joseph↗
NASA NTRS1988

The numerical simulation of transonic separated flow about the complete F-16A

The thin-layer, Reynolds-averaged, Navier-Stokes equations are used to simulate the transonic viscous flow about the complete F-16A fighter aircraft. These computations demonstrate how computational fluid dynamics (CFD) can be used to simulate turbulent viscous flow about realistic aircraft geometries. A zonal grid approach is used to provide adequate viscous grid clustering on all aircraft surfaces. Zonal grids extend inside the F-16A inlet and up to the compressor face while power on conditions are modeled by employing a zonal grid extending from the exhaust nozzle to the far field. A simple solution adaptive grid procedure is used on the wing surface and good agreement with experimental data is obtained. Computations for the F-16A in side slip are also presented.

Flores, Jolen↗
NASA NTRS1988

Effects of stratospheric lapse rate on thunderstorm cloud-top structure in a three-dimensional numerical simulation. I - Some basic results of comparative experiments

The effects of stratospheric temperature lapse rate on cloud top height/temperature structure for strongly sheared, mature, isolated midlatitude thunderstorms are investigated by performing three different experiments with an anelastic, three-dimensional model: (1) with an assumed stratospheric lapse rate of 0 K/km (i.e., the isothermal case), (2) with 3 K/km, and (3) with -3 K/km (i.e., the case of inversion). Kinematic storm structure is very similar in all three cases, especially in the troposphere; a strong quasi-steady updraft evolves and splits into a dominant cyclonic overshooting right-mover and a weaker, anticyclonic left-mover that does not reach the tropopause.

Schlesinger, Robert E.↗