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,621 records · Page 90

Numerical Simulation of Flow Past a Tapered Cylinder

We have computed the unsteady three-dimensional low Reynolds number flow past a tapered cylinder. The spanwise variation in natural shedding frequency results in interesting three-dimensional flow phenomena. Our computed hot-wire and spectral data are very similar to experimental results. The computation was done on the Connection Machine, a massively parallel computer, we highlight the capabilities of the Connection Machine, for computation and visualization of three-dimensional unsteady flow fields.

Jespersen, Dennis↗

Three-Dimensional Navier-Stokes Method with Two-Equation Turbulence Models for Efficient Numerical Simulation of Hypersonic Flows

A new computational efficient 3-D compressible Reynolds-averaged implicit Navier-Stokes method with advanced two equation turbulence models for high speed flows is presented. All convective terms are modeled using an entropy satisfying higher-order Total Variation Diminishing (TVD) scheme based on implicit upwind flux-difference split approximations and arithmetic averaging procedure of primitive variables. This method combines the best features of data management and computational efficiency of space marching procedures with the generality and stability of time dependent Navier-Stokes procedures to solve flows with mixed supersonic and subsonic zones, including streamwise separated flows. Its robust stability derives from a combination of conservative implicit upwind flux-difference splitting with Roe's property U to provide accurate shock capturing capability that non-conservative schemes do not guarantee, alternating symmetric Gauss-Seidel 'method of planes' relaxation procedure coupled with a three-dimensional two-factor diagonal-dominant approximate factorization scheme, TVD flux limiters of higher-order flux differences satisfying realizability, and well-posed characteristic-based implicit boundary-point a'pproximations consistent with the local characteristics domain of dependence. The efficiency of the method is highly increased with Newton Raphson acceleration which allows convergence in essentially one forward sweep for supersonic flows. The method is verified by comparing with experiment and other Navier-Stokes methods. Here, results of adiabatic and cooled flat plate flows, compression corner flow, and 3-D hypersonic shock-wave/turbulent boundary layer interaction flows are presented. The robust 3-D method achieves a better computational efficiency of at least one order of magnitude over the CNS Navier-Stokes code. It provides cost-effective aerodynamic predictions in agreement with experiment, and the capability of predicting complex flow structures in complex geometries with good accuracy.

Bardina, J. E.↗

Studying Turbulence Using Numerical Simulation Databases: Proceedings of the 1996 Summer Program - Part 6

Topics considered include: New approach to turbulence modeling; Second moment closure analysis of the backstep flow database; Prediction of the backflow and recovery regions in the backward facing step at various Reynolds numbers; Turbulent flame propagation in partially premixed flames; Ensemble averaged dynamic modeling. Also included a study of the turbulence structures of wall-bounded shear flows; Simulation and modeling of the elliptic streamline flow.

Source record↗

Numerical Simulations of High-Speed Chemically Reacting Flow

The essentially nonoscillatory (ENO) shock-capturing scheme for the solution of hyperbolic equations is extended to solve a system of coupled conservation equations governing two-dimensional, time-dependent, compressible chemically reacting flow with full chemistry. The thermodynamic properties of the mixture are modeled accurately, and stiff kinetic terms are separated from the fluid motion by a fractional step algorithm. The methodology is used to study the concept of shock-induced mixing and combustion, a process by which the interaction of a shock wave with a jet of low-density hydrogen fuel enhances mixing through streamwise vorticity generation. Test cases with and without chemical reaction are explored here. Our results indicate that, in the temperature range examined, vorticity generation as well as the distribution of atomic species do not change significantly with the introduction of a chemical reaction and subsequent heat release. The actual diffusion of hydrogen is also relatively unaffected by the reaction process. This suggests that the fluid mechanics of this problem may be successfully decoupled from the combustion processes, and that computation of the mixing problem (without combustion chemistry) can elucidate much of the important physical features of the flow.

Ton, V. T.↗

Numerical Simulation of the Effect of Heating on Supersonic Jet Noise

The axisymmetric linearized Euler equations are used to simulate noise amplification and radiation from a supersonic jet. The effect of heating on the noise field of the jet is studied and compared to experimental results. Special attention was given to boundary treatment, and the resulting solution is stable and nearly free from boundary reflections.

Hixon, R.↗

Numerical Simulation of High-Speed Turbulent Reacting Flows

The purpose of this research is to continue our efforts in advancing the state of knowledge in large eddy simulation (LES) methods for computational analysis of high-speed reacting turbulent flows. We have just completed the first year of Phase 3 of this research.

Jaberi, F. A.↗

Numerical Simulation of the 9-10 June 1972 Black Hills Storm Using CSU RAMS

Strong easterly flow of low-level moist air over the eastern slopes of the Black Hills on 9-10 June 1972 generated a storm system that produced a flash flood, devastating the area. Based on observations from this storm event, and also from the similar Big Thompson 1976 storm event, conceptual models have been developed to explain the unusually high precipitation efficiency. In this study, the Black Hills storm is simulated using the Colorado State University Regional Atmospheric Modeling System. Simulations with homogeneous and inhomogeneous initializations and different grid structures are presented. The conceptual models of storm structure proposed by previous studies are examined in light of the present simulations. Both homogeneous and inhomogeneous initialization results capture the intense nature of the storm, but the inhomogeneous simulation produced a precipitation pattern closer to the observed pattern. The simulations point to stationary tilted updrafts, with precipitation falling out to the rear as the preferred storm structure. Experiments with different grid structures point to the importance of removing the lateral boundaries far from the region of activity. Overall, simulation performance in capturing the observed behavior of the storm system was enhanced by use of inhomogeneous initialization.

Nair, U. S.↗

Numerical Simulation of Fluctuations Leading to Noise in a Flap-Edge Flowfield

In this paper we develop an approximate computational framework for simulation of the fluctuating flowfield associated with the complex vortex system seen at the side edge of a flap in a multi-element high-lift airfoil system. The eventual goal of these simulations is to provide an estimate of the spectral content of these fluctuations, in order that the spectrum of the noise generated by such flowfields may be estimated. Results from simulations utilizing this computational framework are shown.

Streett, C. L.↗

Numerical Simulation of Rotation-Driven Plasma Transport In the Jovian Magnetosphere

A Jupiter version of the Rice Convection Model (RCM-J) was developed with support of an earlier NASA SR&T grant. The conversion from Earth to Jupiter included adding currents driven by centrifugal force, reversing the planetary magnetic field, and rescaling various parameters. A series of informative runs was carried out, all of them solving initial value problems. The simulations followed an initial plasma torus configuration as it fell apart by interchange instability. Some conclusions from the simulations were the following: 1. We confirmed that, for conventional values of the torus density and ionospheric conductance, the torus disintegrates by interchange instability on a time scale of approx. one day, which is 1-2 orders of magnitude shorter than the best estimates of the average residence time of plasma in the torus. 2. In the model, the instability could be slowed to an arbitrary degree by the addition of sufficient impounding energetic particles, as suggested earlier by Siscoe et al (1981). However, the observed energetic particles do not seem sufficient to guarantee impoundment (e.g., Mauk et al., 1996). 3. Whether inhibited by impoundment or not, the interchange was found to proceed by the formation of long fingers, which get thinner as they get longer. This picture differed dramatically from the conventional radial-diffusion picture (e.g., Siscoe and Summers (1981)), more superficially with the outward-moving-blob picture (Pontius and Hill, 1989). The obvious limitation of the original RCM-J was that it could not represent a plasma source. We could represent the decay of a pre-existing torus, but we could not represent the way ionization of material from Io continually replenishes the plasma. We consequently were precluded from studying a whole set of fundamental issues of torus theory, including whether the system can come to a steady state.

Wolf, Richard A.↗

Direct Numerical Simulations of Boundary Layer Transition on a Flat Plate

In recent years the techniques of computational fluid dynamics (CFD) have been used to compute flows associated with geometrically complex configurations. However, success in terms of accuracy and reliability has been limited to cases where the effects of turbulence and transition could be modeled in a straightforward manner. Even in simple flows, the accurate computation of skin friction and heat transfer using existing turbulence models has proved to be a difficult task, one that has required extensive fine-tuning of the turbulence models used. In more complex flows (for example, in turbomachinery flows in which vortices and wakes impinge on airfoil surfaces causing periodic transitions from laminar to turbulent flow) the development of a model that accounts for all scales of turbulence and predicts the onset of transition is an extremely difficult task. Fortunately, current trends in computing suggest that it may be possible to perform direct simulations of turbulence and transition at moderate Reynolds numbers in some complex cases in the near future. This presentation will focus on direct simulations of transition and turbulence using high-order accurate finite-difference methods. The advantage of the finite-difference approach over spectral methods is that complex geometries can be treated in a straightforward manner. Additionally, finite-difference techniques are the prevailing methods in existing application codes. An application of high-order-accurate finite-difference methods to direct simulations of transition and turbulence in a spatially evolving boundary layer subjected to high levels of freestream turbulence will be presented.

Rai, Man Mohan↗

Numerical simulations of infrared emitting regions

The development of CLOUDY, a large scale code designed to compute the spectrum of gas into photoionization or collisional balance, is presented. The major upgrading of CLOUDY are given together with a discussion of the remaining problems and the future plans for improving the atomic data base. The completeness of the line list included in CLOUDY, the results of the Opacity project and the Iron project are discussed.

vanHoof, P. A. M.↗

A Quantitative Investigation of Entrainment and Detrainment in Numerically Simulated Convective Clouds: Simulations of Cumulonimbus Clouds - Pt. 2

Deep cumulonimbus clouds are simulated using a model that makes accurate diagnoses of entrainment and detrainment rates and of the properties of entrained and detrained air. Clouds generated by a variety of initial thermodynamic soundings are compared. In the simulations, updraft entrainment rates are large near and above cloud base, through the entire depth of the conditionally unstable layer. Stronger updrafts in a more unstable environment are better able to entrain relatively undisturbed environmental air, while weaker updrafts can entrain only air that has been modified by the clouds. When the maximum buoyancy is large, the updraft includes parcels with a wide range of buoyancies, while weaker clouds are more horizontally uniform. Strong downdrafts originate from levels at which updrafts detrain, and their mass flux depends on the mass flux of the updraft. The magnitude of mixing between cloud and environment, not the entrainment rate, varies inversely with the cloud radius. How much of the mixed air is entrained depends on the buoyancy.

Cohen, Charles↗

Numerical Simulation of Three-Dimensional Boattail Afterbody Flowfields

The thin shear-layer approximations of the three-dimensional, compressible Navier-Stokes equations are solved for subsonic, transonic, and supersonic now over axisymmetric boattail bodies at moderate angles of attack. The plume is simulated by a solid body configuration identical to those used In experimental tests. An implicit algorithm of second-order accuracy is used to solve the equations on the ILLIAC 4 computer. The turbulence is expressed by an algebraic model applicable to three-dimensional flowfields with moderate separation. The formulation used is attractive in its independence of boundary-layer parameters. Such a simple model, however, is incapable of supporting detailed quantitative descriptions of complex shear flows. Never-the-less, good qualitative comparisons are found with three different sets of experimental date. Quantitative improvement will depend on improved turbulence transport descriptions.

Deiwert, G. S.↗

Acoustic Treatment Design Scaling Methods: Numerical Simulation of the Nonlinear Acoustic Impedance of a Perforated Plate Single-Degree-of-Freedom Resonator Using a Time-Domain Finite Difference Method - Volume 4

Single-degree-of-freedom resonators consisting of honeycomb cells covered by perforated facesheets are widely used as acoustic noise suppression liners in aircraft engine ducts. The acoustic resistance and mass reactance of such liners are known to vary with the intensity of the sound incident upon the panel. Since the pressure drop across a perforated liner facesheet increases quadratically with the flow velocity through the facesheet, this is known as the nonlinear resistance effect. In the past, two different empirical frequency domain models have been used to predict the Sound Pressure Level effect of the incident wave on the perforated liner impedance, one that uses the incident particle velocity in isolated narrowbands, and one that models the particle velocity as the overall velocity. In the absence of grazing flow, neither frequency domain model is entirely accurate in predicting the nonlinear effect that is measured for typical perforated sheets. The time domain model is developed in an attempt to understand and improve the model for the effect of spectral shape and amplitude of multi-frequency incident sound pressure on the liner impedance. A computer code for the time-domain finite difference model is developed and predictions using the models are compared to current frequency-domain models.

Kraft, R. E.↗

Numerical simulation of MHD turbulence in three dimensions

The evolution of Alfvenic turbulence in 3D spherical geometry can now be studied. In simulations, a fast stream is sandwiched between two slower streams. The inflow is both supersonic and superAlfvenic. Alfven waves entering the box are convected into the medium and interact nonlinearly with the velocity shear and with any structures (i.e., flux tubes) that might be present. These initial simulations suggest that velocity shear, even in spherical geometry, is able to drive a turbulent cascade which results in approximately Kolmogoroff-like power spectra.

Goldstein, M. L.↗

Numerical Simulation of High-Speed Turbulent Reacting Flows

The objectives of this research are: (1) to develop and implement a new methodology for large eddy simulation of (LES) of high-speed reacting turbulent flows. (2) To develop algebraic turbulence closures for statistical description of chemically reacting turbulent flows.

Givi, P.↗

The Dependence on Grid Resolution of Numerically Simulated Convective Cloud Systems Using Ice Microphysics

Mesoscale research and forecast models are increasingly being used at horizontal resolutions of 1-8 km to simulate a variety of precipitating systems. When the model is used to simulate convective systems, it is uncertain to what extent the dynamics and microphysics of convective updrafts can be resolved with grids larger than 1 km. In this study, two- and three-dimensional versions of the Goddard Cumulus Ensemble model are used to determine the impact of horizontal grid resolution on the behavior of the simulated storms and on the characteristics of the cloud microphysical fields. It will be shown that as resolution decreases from about 1 km to greater than 3 km, there is a fairly rapid degradation of the storm structure in the form of reduced convective mass fluxes, updraft tilts, and cloud microphysics. A high-resolution simulation of hurricane outer rainbands using the MM5 mesoscale model shows also that there can be a substantial modification of the key microphysical processes that contribute to rainfall as a result of reducing the horizontal resolution.

Braun, Scott A.↗

High Resolution Numerical Simulations and Diagnostic Studies Atmospheric Transport During TRACE-A

The enclosed publications constitute our final report. These are publications completed in referred journals. The completed work includes: Meteorological conditions associated with vertical distributions of aerosols off the west coast of Africa. TRACE-A trajectory intercomparison 2. Isentropic and kinematic methods. Chemical characteristics of tropospheric air over the tropical South Atlantic Ocean: Relationship to trajectory history. Passive Tracer Transport Relevant to the TRACE-A Experiment. The meteorological environment of the tropospheric ozone maximum over the tropical south Atlantic Ocean. : Influence of a middle-latitude cyclone on tropospheric ozone distributions during a period of TRACE-A. All of this work provided the meteorological background for the Transport of Atmospheric Chemistry near the Equator (TRACE-A) project sponsored by Dr,, Joe McNeal of NASA. Our principal findings are that bio-mass burning sources from Africa and South America did contribute to the accumulation of tropospheric ozone over the tropical Atlantic Ocean. Other findings relate to circulations, advections of Ozone and other plausible sources for the TOMS based ozone maximum over the Atlantic Ocean.

Krishnamurti, T. N.↗