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 271 records · Page 15
All sourcesNASA NTRSDOE OSTI
NASA NTRS1984

Numerical simulations of double layers and auroral electric fields

Recent one-dimensional and two-dimensional numerical simulations of double layers (DLs) in the electric fields of the auroral plasma are reviewed, with reference to observational data. It is found that two-dimensional DLs driven by current sheets of finite thickness have different characteristics, depending on whether the layer thickness is less than or much greater than the ion gyroradius: When thickness is less than ion gyroradius, V-shaped DLs form with nearly equal parallel and perpendicular potential drops; when layer thickness is much greater than ion gyroradius the major parallel potential drop occurs outside the current sheet and the perpendicular electric fields are localized at the edges of the current sheet. It is shown that some features of the simulated fields, such as the amplitudes and scale lengths, are qualitatively similar to those observed in space.

Singh, N.↗
NASA NTRS1989

A review of quasi-coherent structures in a numerically simulated turbulent boundary layer

Preliminary results of a comprehensive study of the structural aspects of a numerically simulated number turbulent boundary layer are presented. A direct Navier-Stokes simulation of a flat-plate, zero pressure gradient boundary layer at Re0 = 670 was used. Most of the known nonrandom, coherent features of turbulent boundary layers are confirmed in the simulation, and several new aspects of their spatial character are reported. The spatial relationships between many of the various structures are described, forming the basis for a more complete kinematical picture of boundary layer physics than has been previously known. In particular, the importance of vortex structures of various forms to the generation of Reynolds shear stress is investigated.

Robinson, S. K.↗
NASA NTRS1990

Active turbulence control in wallbounded flow using direct numerical simulations

An exploratory study of concepts for active control of turbulent boundary layers using the direct numerical simulation technique was performed. Significant drag reduction was achieved when the surface boundary condition was modified such that it could suppress the large-scale structures present in the wall region. This was achieved by prescribing the normal component of velocity at the wall to be 180 deg out of phase with the normal velocity slightly above the wall at each instant. The drag reduction was accompanied with significant reduction in the intensity of the wall-layer structures and reductions in the magnitude of Reynolds stresses throughout the flow. Suitability of wall-pressure and shear-stress fluctuations for detection of flow structures above the wall was examined. A preliminary result obtained by applying the present control strategy to a transitional flow is also briefly described, from which one can infer a possible linkage between the control strategy and flow stability.

Kim, J.↗
NASA NTRS1992

The three-point function in an ensemble of numerical simulations

The three-point function in Fourier space is evaluated for an ensemble of 2D numerical simulations with nine initial power spectra distinguished by spectral index and cutoff. To remove the main dependences on scale and time, results are presented as the reduced amplitude Q in the hierarchical model. To lowest nonvanishing order in perturbation theory, Q is a constant, independent of length scale, of time, and of initial spectrum. The results show that in the nonlinear regime of evolution, normalizing to Q does remove the main variations, but systematic dependences on spectral index and cutoff remain at late times that depend on properties of the initial spectrum. Models with initial power cut off on a scale that has been evolved far into the nonlinear regime are indistinguishable from models with no cutoff.

Fry, J. N.↗
NASA NTRS1993

Numerical simulation of the Space Shuttle Launch Vehicle flowfield with real gas solid rocket plume effects

A numerical simulation of the external ascent flowfield of the Space Shuttle Launch Vehicle (SSLV) was carried out at the freestream Mach number 1.25, the angle of attack -5.1 deg, and the flight Reynolds number 3.25 x 10 exp 6/ft. The simulation is an extension of the solution by Kandula et al. (1991) and incorporates variable gamma effects with a high fidelity grid appropriate for a flight Reynolds number flow solution. Three-dimensional SSLV steady-state simulations with either perfect-gas or variable-gamma-gas Redesigned-Solid-Rocket-Motor (RSRM) plumes are computed on a 14-grid overlapping body-conforming grid system, and the influence of the RSRM exhaust plumes on the SSLV component pressure distributions and integrated loads is examined.

Slotnick, J. P.↗
NASA NTRS2001

Numerical Simulations of the Slingatron

The slingatron mass accelerator is described for several track configurations (shapes), and numerical simulations of this accelerating mass traversing a given track configuration are presented. The sled is modeled as a point mass that interacts with the slingatron track using both a conventional and a new empirical velocity dependent friction law. The closed loop circular slingatron was found to produce high maximum sled velocities provided the gyration angular speed is always increasing. In contrast several spiral shaped slingatron tracks reveal that high maximum sled velocities are obtainable with the gyration speed held constant. In fact, a slingatron constructed out of semi-circles is shown capable of generating high velocity sleds in such a way that no initial sled injection is necessary. Choosing the proper initial gyration phase with an empirically determined friction model allows the mass sled to gain ever-increasing velocities when placed in a semi-circle slingatron. The sled bearing pressure and its total acceleration are examined and presented.

Cooper, Gene R.↗
NASA NTRS2012

Three-Dimensional Numerical Simulations of Equatorial Spread F: Results and Observations in the Pacific Sector

A three-dimensional numerical simulation of plasma density irregularities in the postsunset equatorial F region ionosphere leading to equatorial spread F (ESF) is described. The simulation evolves under realistic background conditions including bottomside plasma shear flow and vertical current. It also incorporates C/NOFS satellite data which partially specify the forcing. A combination of generalized Rayleigh-Taylor instability (GRT) and collisional shear instability (CSI) produces growing waveforms with key features that agree with C/NOFS satellite and ALTAIR radar observations in the Pacific sector, including features such as gross morphology and rates of development. The transient response of CSI is consistent with the observation of bottomside waves with wavelengths close to 30 km, whereas the steady state behavior of the combined instability can account for the 100+ km wavelength waves that predominate in the F region.

ALTAIR↗
NASA NTRS1983

Numerical simulation of two-dimensional heat transfer in composite bodies with application to de-icing of aircraft components

Transient, numerical simulations of the de-icing of composite aircraft components by electrothermal heating were performed for a two dimensional rectangular geometry. The implicit Crank-Nicolson formulation was used to insure stability of the finite-difference heat conduction equations and the phase change in the ice layer was simulated using the Enthalpy method. The Gauss-Seidel point iterative method was used to solve the system of difference equations. Numerical solutions illustrating de-icer performance for various composite aircraft structures and environmental conditions are presented. Comparisons are made with previous studies. The simulation can also be used to solve a variety of other heat conduction problems involving composite bodies.

Chao, D. F. K.↗
NASA NTRS1987

Numerical simulation of transonic propeller flow using a 3-dimensional small disturbance code employing novel helical coordinates

The numerical simulation of three-dimensional transonic flow about propeller blades is discussed. The equations for the unsteady potential flow about propellers is given for an arbitrary coordinate system. From this the small disturbance form of the equation is derived for a new helical coordinate system. The new coordinate system is suited to propeller flow and allows cascade boundary conditions to be applied straightforward. A numerical scheme is employed which solves the steady flow as an asymptotic limit of unsteady flow. Solutions are presented for subsonic and transonic flow about a 5 percent thick bicircular arc blade of an eight bladed cascade. Both high and low advance ratio cases are given which include a lifting case as well as nonlifting cases. The nonlifting cases are compared to solutions from a Euler code.

Snyder, Aaron↗
NASA NTRS1989

Numerical simulation of hyperbolic heat conduction with convection boundary conditions and pulse heating effects

The paper describes the numerical simulation of hyperbolic heat conduction with convection boundary conditions. The effects of a step heat loading, a sudden pulse heat loading, and an internal heat source are considered in conjunction with convection boundary conditions. Two methods of solution are presened for predicting the transient behavior of the propagating thermal disturbances. In the first method, MacCormack's predictor-corrector method is employed for integrating the hyperbolic system of equations. Next, the transfinite element method, which employs specially tailored elements, is used for accurately representing the transient response of the propagating thermal wave fronts. The agreement between the results of various numerical test cases validate the representative behavior of the thermal wave fronts. Both methods represent hyperbolic heat conduction behavior by effectively modeling the sharp discontinuities of the propagating thermal disturbances.

Glass, David E.↗
NASA NTRS1985

Numerical simulation and comparison with experiment for self-excited oscillations in a diffuser flow

This paper describes numerical simulations of self-excited oscillations in a two-dimensional transonic diffuser flow obtained by solving the Navier-Stokes equations with a two-equation turbulence model. Comparisons were made between the computational results and experimental data. For the mean flowfields, the agreement between computation and experiment is good for the wall pressures, shock location, and the separation and reattachment points. However, the thickness of the computed recirculation zone is about 50 percent of the measured thickness. For the fluctuating flowfields, a great deal of qualitative similarity exists between the computation and experiment; however, the predicted oscillation frequency is about 50 percent higher than the measured value. The formation of a succession of downstream-traveling counter-rotating vortices, as seen experimentally, is also vividly displayed in the numerical results.

Hsieh, T.↗
NASA NTRS1985

Direct numerical simulations of a reacting mixing layer with chemical heat release

In order to study the coupling between chemical heat release and fluid dynamics, direct numerical simulations of a chemically reacting mixing layer with heat release are performed. The fully compressible equations as well as an approximate set of equations that is asymptotically valid for low-Mach-number flows are treated. These latter equations have the computational advantage that high-frequency acoustic waves have been filtered out, allowing much larger time steps to be taken in the numerical solution procedure. A detailed derivation of these equations along with an outline of the numerical solution technique is given. Simulation results indicate that the rate of chemical product formed, the thickness of the mixing layer, and the amount of mass entrained into the layer all decrease with increasing rates of heat release.

Mcmurtry, P. A.↗
NASA NTRS1987

Numerical simulation of channel flow transition, resolution requirements and structure of the hairpin vortex

Three-dimensional, nonlinear numerical simulations are presented for the K-type and H-type transitions for channel flow. There are two objectives. The first is to establish firmly the resolution requirements for the various stages in the transition process. Comparisons between calculations on various grids suggest a set of guidelines for maintaining a physically meaningful calculation. The second objective is to map out the structure of the hairpin vortices which arise in K-type and H-type transitions in channel flow, to the latest stage currently feasible. Flow field details are presented for both a subcritical Reynolds number of 1500 and a supercritical Reynolds number of 8000. The diagnostics include illustrations of the vertical shear, streamwise and spanwise vorticity, helicity, vortex stretching, and vortex diffusion fields.

Krist, Steven E.↗
NASA NTRS1991

The thermal stability of coronal loops - Numerical simulations

The radiative stability of thermally isolated coronal loops with free-flow boundary conditions by nonlinear numerical simulation are studied. A chromosphere-to-corona loop equilibrium (including the option of a deep chromosphere) is first established by following the nonlinear evolution from an initial isothermal state with rigid boundaries. The end conditions are then changed, to allow free flow and to fix the temperature, and investigate the response to nonisobaric perturbations. Within a family of loops of the same pressure, long hot loops to be stable and short cool loops to be unstable to the thermal chromospheric-expansion mode are found. The stable cases remain so, even when long chromospheric ends and/or gravity are added. In those cases which are unstable, the subsequent nonlinear evolution which exhibits swelling of the chromosphere until the entire loop becomes cool and dense are followed.

Mok, Y.↗
NASA NTRS1992

The program FANS-3D (finite analytic numerical simulation 3-dimensional) and its applications

In this study, the program named FANS-3D (Finite Analytic Numerical Simulation-3 Dimensional) is presented. FANS-3D was designed to solve problems of incompressible fluid flow and combined modes of heat transfer. It solves problems with conduction and convection modes of heat transfer in laminar flow, with provisions for radiation and turbulent flows. It can solve singular or conjugate modes of heat transfer. It also solves problems in natural convection, using the Boussinesq approximation. FANS-3D was designed to solve heat transfer problems inside one, two and three dimensional geometries that can be represented by orthogonal planes in a Cartesian coordinate system. It can solve internal and external flows using appropriate boundary conditions such as symmetric, periodic and user specified.

Bravo, Ramiro H.↗
NASA NTRS1994

Scalability study of parallel spatial direct numerical simulation code on IBM SP1 parallel supercomputer

The implementation and the performance of a parallel spatial direct numerical simulation (PSDNS) code are reported for the IBM SP1 supercomputer. The spatially evolving disturbances that are associated with laminar-to-turbulent in three-dimensional boundary-layer flows are computed with the PS-DNS code. By remapping the distributed data structure during the course of the calculation, optimized serial library routines can be utilized that substantially increase the computational performance. Although the remapping incurs a high communication penalty, the parallel efficiency of the code remains above 40% for all performed calculations. By using appropriate compile options and optimized library routines, the serial code achieves 52-56 Mflops on a single node of the SP1 (45% of theoretical peak performance). The actual performance of the PSDNS code on the SP1 is evaluated with a 'real world' simulation that consists of 1.7 million grid points. One time step of this simulation is calculated on eight nodes of the SP1 in the same time as required by a Cray Y/MP for the same simulation. The scalability information provides estimated computational costs that match the actual costs relative to changes in the number of grid points.

Hanebutte, Ulf R.↗
NASA NTRS1996

Computational aeroacoustics and numerical simulation of supersonic jets

The research project has been a computational study of computational aeroacoustics algorithms and numerical simulations of the flow and noise of supersonic jets. During this study a new method for the implementation of solid wall boundary conditions for complex geometries in three dimensions has been developed. In addition, a detailed study of the simulation of the flow in and noise from supersonic circular and rectangular jets has been conducted. Extensive comparisons have been made with experimental measurements. A summary of the results of the research program are attached as the main body of this report in the form of two publications. Also, the report lists the names of the students who were supported by this grant, their degrees, and the titles of their dissertations. In addition, a list of presentations and publications made by the Principal Investigators and the research students is also included.

Morris, Philip J.↗