Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Direct 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 811 records · Page 45

A Stereo Imaging Velocimetry Technique for Analyzing Structure of Flame Balls at Low Lewis-Number (SOFBALL) Data

Stereo Imaging Velocimetry (SIV) is a NASA Glenn Research Center (GRC) developed fluid physics technique for measuring threedimensional (3-D) velocities in any optically transparent fluid that can be seeded with tracer particles. SIV provides a means to measure 3-D fluid velocities quantitatively and qualitatively at many points. This technique provides full-field 3-D analysis of any optically clear fluid or gas experiment using standard off-the-shelf CCD cameras to provide accurate and reproducible 3-D velocity profiles for experiments that require 3-D analysis. A flame ball is a steady flame in a premixed combustible atmosphere which, due to the transport properties (low Lewis-number) of the mixture, does not propagate but is instead supplied by diffusive transport of the reactants, forming a premixed flame. This flame geometry presents a unique environment for testing combustion theory. We present our analysis of flame ball phenomena utilizing SIV technology in order to accurately calculate the 3-D position of a flame ball(s) during an experiment, which can be used as a direct comparison of numerical simulations.

McDowell, Mark↗

Numerical simulation of rarefied flow through a slit. I - Direct simulation Monte Carlo results

The pressure-driven flow of a rarefied monatomic gas through a two-dimensional slit is simulated using the direct simulation Monte Carlo technique. Of particular interest is the change in flow field structure as pressure ratio and Knudsen number are varied. Comparisons are made to quantify the limits of validity of free-molecular theory and approximate, nearly free-molecular iterative methods. Also addressed is the sensitivity of the numerical solutions to grid structure and boundary conditions. The free-molecular theory is found to predict quantitative flow field properties (e.g., centerline velocities or downstream flux profiles) reasonably well for large finite Knudsen number with the error dependent on the pressure ratio. The nearly free-molecular corrections are shown to have limited range of applicability. A previously derived parameter is found to correlate total mass flux well as a function of pressure ratio and Knudsen number over a large portion of the transitional regime.

Wadsworth, D. C.↗

Investigation of turbulent shear flows by numerical simulation

A data base (Kim et al., 1986) generated by direct numerical integration of the Navier-Stokes equations for unsteady three-dimensional turbulent channel flow at low Reynolds numbers is applied to investigate the physical processes acting in wall-bounded shear flows. The techniques used to construct the data base are outlined; the modeling of turbulence, streaks and streamwise vortices, horseshoe and inverted-horseshoe vortices, and vortex rings is described; the conditional sampling techniques employed are discussed; and typical results are presented graphically. An organized vorticity-field structure is revealed and found to be associated with the bursting process and a local peak in the wall pressure, indicating a possible approach to reducing viscous drag.

Kim, J.↗

Spacecraft outgas ambient flow interaction

An investigation has been made of the extent to which outgassed or emitted molecules return to a spacecraft as a result of intermolecular collisions with ambient freestream molecules. The governing parameters are discussed and heuristic arguments are used to predict the combination of parameters that is most likely to describe the return flux ratio. This is then tested by an extensive set of numerical calculations using the direct simulation Monte Carlo method. Computations have been made for the flow past spheres and past circular cylinders with their axes normal in the stream. The numerical results lead to empirical expressions that can be used for engineering estimates in spacecraft contamination studies.

Bird, G. A.↗

Parallel implementation of a particle simulation for modeling rarefied gas dynamic flow

When the conditions of flow are rarefied and hypersonic, a more suitable alternative to the use of the Navier-Stokes equations for developing a numerical solution is the Direct Simulation Monte Carlo method (DSMC), a method of simulation which employs a large number of particles in modeling a rarefied gas. The performance of a parallel DSMC code developed for the Intel iPSC/860 Touchstone Gamma prototype computer is studied and the scaleup is found to be very nearly over the range of 16-128 processors.

Fallavollita, M. A.↗

Axisymmetric confined turbulent jet directed towards the liquid surface from below

A numerical simulation is presented of an axisymmetric turbulent jet discharging axially from below into a cylindrical tank and directed towards the liquid vapor interface. The liquid vapor interface is assumed to be flat and shear free. The k-epsilon turbulence model is used to calculate the eddy viscosity. The turbulence intensity distribution and the length scale associated with the k-epsilon model are calculated as functions of jet flow rates and systems parameters. Numerical results are compared with appropriate experimental data. The problems associated with the free surface boundary conditions for the turbulent quantities are discussed.

Hasan, Mohammad M.↗

Axisymmetric confined turbulent jet directed towards the liquid surface from below

A numerical simulation is presented of an axisymmetric turbulent jet discharging axially from below into a cylindrical tank and directed towards the liquid vapor interface. The liquid vapor interface is assumed to be flat and shear free. The k-epsilon turbulence model is used to calculate the eddy viscosity. The turbulence intensity distribution and the length scale associated with the k-epsilon model are calculated as functions of jet flow rates and systems parameters. Numerical results are compared with appropriate experimental data. The problems associated with the free surface boundary conditions for the turbulent quantities are discussed.

Hasan, Mohammad M.↗

Advances in turbulent transport modeling based on direct simulations of turbulence

Recent work using numerical simulations to study the modeling of scalar transport is outlined. A model developed by reference to simulations for homogeneous shear flow is presented. This model predicts cross-gradient transports that can give rise to the illusion of counter-gradient transport. The model is found to work remarkably well when compared to simulations for channel flow, but fails to predict important features of the scalar fluxes in mixing layers. Reasons for this failure are discussed and suggestions are made for future work.

Reynolds, W. C.↗

Numerical simulation of nonlinear interactions in channel and boundary-layer transition

The basic numerical techniques used for the direct simulation of transition to turbulence are described. Applications are presented for the computation of finite amplitude, two-dimensional steady states, the simulation of the secondary instability of channels and boundary-layers, and the simulation of the breakdown of detached shear layers.

Zang, T. A.↗

Measurement and analysis of a small nozzle plume in vacuum

Measurements of Pitot pressure and flow angle were made in the plume of a nozzle flowing nitrogen and exhausting to a vacuum. The measurements were compared to results from a numerical simulation of the flow that was based on kinetic theory and used the direct-simulation Monte Carlo (DSMC) method. Numerical results were compared with measurements made in the plume at various axial and radial stations. Total pressure measurements were made with Pitot tubes sized for specific regions of the plume. Flow angle measurements were made with a conical probe. The measurement area for flow angle extended to 160 mm (5 exit diameters) downstream of the nozzle exit plane and radially to 60 mm (1.9 exit diameters) from the plume axis. The total pressure measurements extended 480 mm (16 exit diameters) downstream and radially to 60 mm. Comparisons of computed results from the DSMC method with measurements of flow angle displayed improved agreement with increasing distance from the exit plane. Pitot pressures computed from the DSMC method were in reasonably good agreement with experimental results over the entire measurement area.

Penko, Paul F.↗

DSMC simulations of Mach 20 nitrogen flows about a 70 degree blunted cone and its wake

Numerical results obtained with the direct simulation Monte Carlo (DSMC) method are presented for Mach 20 nitrogen flow about a 70-deg blunted cone. The flow conditions simulated are those that can be obtained in existing low-density hypersonic wind tunnels. Three sets of flow conditions are simulated with freestream Knudsen numbers ranging from 0.03 to 0.001. The focus is to characterize the wake flow under rarefied conditions. This is accomplished by calculating the influence of rarefaction on wake structure along with the impact that an afterbody has on flow features. This data report presents extensive information concerning flowfield features and surface quantities.

Moss, James N.↗

Transient nucleation in condensed systems

Using classical nucleation theory we consider transient nucleation occurring in a one-component, condensed system under isothermal conditions. We obtain an exact closed-form expression for the time dependent cluster populations. In addition, a more versatile approach is developed: a numerical simulation technique which models directly the reactions by which clusters are produced. This simulation demonstrates the evolution of cluster populations and nucleation rate in the transient regime. Results from the simulation are verified by comparison with exact analytical solutions for the steady state. Experimental methods for measuring transient nucleation are assessed, and it is demonstrated that the observed behavior depends on the method used. The effect of preexisting cluster distributions is studied. Previous analytical and numerical treatments of transient nucleation are compared to the solutions obtained from the simulation. The simple expressions of Kashchiev are shown to give good descriptions of the nucleation behavior.

Kelton, K. F.↗

Computational fluid dynamics - Current capabilities and directions for the future

Computational fluid dynamics (CFD) has made great strides in the detailed simulation of complex fluid flows, including some of those not before understood. It is now being routinely applied to some rather complicated problems and starting to affect the design cycle of aerospace flight vehicles and their components. It is being used to complement, and is being complemented by, experimental studies. Several examples are presented in the paper to illustrate the current state of the art. Included is a discussion of the barriers to accomplishing the basic objective of numerical simulation. In addition, the directions for the future in the discipline of computational fluid dynamics are addressed.

Kutler, Paul↗

On the Generation of Turbulent Inflow Conditions for Boundary Layer Simulations

Turbulent flows that exhibit inhomogeneities in the streamwise direction pose a particular challenge to numerical simulation approaches due to the need to prescribe time-dependent turbulent inflow conditions. In most cases the flow downstream is more or less 'driven' by the conditions at the inlet, making it necessary to specify realistic turbulent fluctuations that are in equilibrium with the assumed mean flow. This requirement often dictates that the inflow data should satisfy the Navier-Stokes equations, which in turn implies that an independent simulation be used to generate them. Detailed simulations for the purpose of creating inflow conditions can be costly and thus certain levels of approximation are desirable. In this paper we shall focus on an approximate yet accurate method for generating inflow conditions for spatially-developing boundary layer simulations. The proposed method is essentially a simplification of the method of Spalart and Leonard (1985), who devised an ingenious transformation that allows for the calculation of spatially evolving boundary layers in conjunction with periodic boundary conditions applied in the streamwise direction. While this method is elegant and highly accurate, it is more complicated than is necessary for the purpose of generating inflow data. A few key approximations are used in this work to arrive at a 'modified Spalart method' that is very easy to implement and efficient to use. The new method is shown to yield results that compare well with the computations of Spalart (1988). When used as a means of generating inflow data, the modified Spalart method is shown to be superior to existing approaches.

Lund, T. S.↗

SIM_EXPLORE: Software for Directed Exploration of Complex Systems

Physics-based numerical simulation codes are widely used in science and engineering to model complex systems that would be infeasible to study otherwise. While such codes may provide the highest- fidelity representation of system behavior, they are often so slow to run that insight into the system is limited. Trying to understand the effects of inputs on outputs by conducting an exhaustive grid-based sweep over the input parameter space is simply too time-consuming. An alternative approach called "directed exploration" has been developed to harvest information from numerical simulators more efficiently. The basic idea is to employ active learning and supervised machine learning to choose cleverly at each step which simulation trials to run next based on the results of previous trials. SIM_EXPLORE is a new computer program that uses directed exploration to explore efficiently complex systems represented by numerical simulations. The software sequentially identifies and runs simulation trials that it believes will be most informative given the results of previous trials. The results of new trials are incorporated into the software's model of the system behavior. The updated model is then used to pick the next round of new trials. This process, implemented as a closed-loop system wrapped around existing simulation code, provides a means to improve the speed and efficiency with which a set of simulations can yield scientifically useful results. The software focuses on the case in which the feedback from the simulation trials is binary-valued, i.e., the learner is only informed of the success or failure of the simulation trial to produce a desired output. The software offers a number of choices for the supervised learning algorithm (the method used to model the system behavior given the results so far) and a number of choices for the active learning strategy (the method used to choose which new simulation trials to run given the current behavior model). The software also makes use of the LEGION distributed computing framework to leverage the power of a set of compute nodes. The approach has been demonstrated on a planetary science application in which numerical simulations are used to study the formation of asteroid families.

Burl, Michael↗

Nonlinear transient finite element analysis of rotor-bearing-stator systems

This paper extends the finite element scheme to handle the highly nonlinear interfacial fields generated in the fluid filled annulli of squeeze film and journal bearings so as to model the transient response of rotor-bearing-stator systems. Since such simulations are highly nonlinear, direct numerical integration schemes are employed to generate the overall response. In this context, the paper gives consideration to such items as (1) numerical efficiency/stability, (2) comparison of implicit and explicit schemes, (3) determines extent of response nonlinearity as well as (4) extensively benchmarks the overall concept/methodologies.

Padovan, J.↗

Nonequilibrium effects for hypersonic transitional flows

Presented are the results of numerical simulations of hypersonic flow about blunt cones and hemispherical nose configurations for reentry velocities of 7.5 and 10 km/s. Cone half angles 0, 5, and 10 deg are considered at zero angle of incidence; however, the focus is for the 5 deg cone. The body size and altitude ranges considered (70 to 110 km) are such that the flow is in the transitional regime. Translational, thermodynamic, and chemical nonequilibrium effects are considered in the numerical simulation by utilizing the direct simulation Monte Carlo (DSMC) method of Bird. The DSMC results are compared with those obtained with viscous shock-layer and Navier-Stokes methods. Comparisons between the DSMC and continuum calculations show the altitude range where differences in flowfield structure and surface quantities become significant. The current calculations show that the binary scaling similitude provides a means of correlating the blunt body surface quantities in the hypersonic, transitional regime. Furthermore, for the higher velocity entry conditions, the results highlight some of the concerns in the application of multitemperature continuum formulations, particularly the use of some proposed functional relations for the chemical rate constants under thermodynamic nonequilibrium conditions.

Moss, James N.↗

Rarefied flow past a flat plate at incidence

Results of a numerical study using the direct simulation Monte Carlo (DSMC) method are presented for the transitional flow about a flat plate at 40 deg incidence. The plate has zero thickness and a length of 1.0 m. The flow conditions simulated are those experienced by the Shuttle Orbiter during reentry at 7.5 km/s. The range of freestream conditions are such that the freestream Knudsen number values are between 0.02 and 8.4, i.e., conditions that encompass most of the transitional flow regime. The DSMC simulations show that transitional effects are evident when compared with free molecule results for all cases considered. The calculated results demonstrate clearly the necessity of having a means of identifying the effects of transitional flow when making aerodynamic flight measurements as are currently being made with the Space Shuttle Orbiter vehicles. Previous flight data analyses have relied exclusively on adjustments in the gas-surface interaction models without accounting for the transitional effect which can be comparable in magnitude. The present calculations show that the transitional effect at 175 km would increase the Space Shuttle Orbiter lift-drag ratio by 90 percent over the free molecule value.

Dogra, Virendra K.↗