Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Turbulence 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 793 records · Page 44

The NASA Ames Hypersonic Combustor-Model Inlet CFD Simulations and Experimental Comparisons

Computations have been performed on a three-dimensional inlet associated with the NASA Ames combustor model for the hypersonic propulsion experiment in the 16-inch shock tunnel. The 3-dimensional inlet was designed to have the combustor inlet flow nearly two-dimensional and of sufficient mass flow necessary for combustion. The 16-inch shock tunnel experiment is a short duration test with test time of the order of milliseconds. The flow through the inlet is in chemical non-equilibrium. Two test entries have been completed and limited experimental results for the inlet region of the combustor-model are available. A number of CFD simulations, with various levels of simplifications such as 2-D simulations, 3-D simulations with and without chemical reactions, simulations with and without turbulent conditions, etc., have been performed. These simulations have helped determine the model inlet flow characteristics and the important factors that affect the combustor inlet flow and the sensitivity of the flow field to these simplifications. In the proposed paper, CFD modeling of the hypersonic inlet, results from the simulations and comparison with available experimental results will be presented.

Venkatapathy, E.↗

On the dynamics of small-scale vorticity in isotropic turbulence

It was previously shown that the strong vorticity in isotropic turbulence is organized into tubular vortices ('worms') whose properties were characterized through the use of full numerical simulations at several Reynolds numbers. At the time most of the observations were kinematic, and several scaling laws were discovered for which there was no theoretical explanation. In the meantime, further analysis of the same fields yielded new information on the generation of the vortices, and it was realized that even if they had to be formed by stretching, they were at any given moment actually compressed at many points of their axes. This apparent contradiction was partially explained by postulating axial inertial waves induced by the nonuniformity of the vortex cores, which helped to 'spread' the axial strain and allowed the vortices to remain compact even if not uniformly stretched. The existence of such solutions was recently proved numerically. The present report discusses a set of new numerical simulations of isotropic turbulence, and a reanalysis of the old ones, in an effort to prove or disprove the presence of these waves in actual turbulent flows and to understand the dynamics, as opposed to the kinematics, of the vortices.

Jimenez, Javier↗

The Effect of Boundary Layer Character on Stochastic Rotor Blade Vortex Shedding Noise

This work illustrates the effect of a rotor blade’s boundary layer on the broadband LBL-VS self-noise emitted from an optimum hovering rotor through experimental and multi-fidelity computational studies. Blade surface roughness effects associated with different manufacturing techniques and the effect of adding a boundary layer trip were shown to decrease LBL-VS noise by upwards of 30 dB at the frequency of maximum emission with a slight penalty in aerodynamic performance when compared with smooth rotor blades. Low-fidelity 2-D viscous flow analysis verified the presence of laminar separation bubbles on the rotor blades which are responsible for LBL-VS noise. Three high-fidelity lattice-Boltzmann simulations were conducted with different wall-functions to predict the boundary layer character correspondent to their experimental counterpart and the resultant presence or absence of LBL-VS noise. Excellent aerodynamic and aeroacoustic agreement was seen between the lattice-Boltzmann simulations and the experimental data for the cases with surface roughness and the boundary layer trip, with the broadband noise results from the simulation with fully turbulent wall-functions diverging from the experimental results above 5 kHz. The transitional wall-function simulation which emulated the smooth experimental blades underpredicted thrust by 14% and broadband noise by a minimum of 10 dB with an accurately predicted broadband noise trend.

Christopher S Thurman↗

Physics-Based Spectra of Accretion Disks around Black Holes

The purpose of this grant was to begin the process of deriving the light output of accretion disks around black holes directly from the actual processes that inject heat into the accreting matter, rather than from guessed dependences of heating rate on physical parameters. At JHU, the effort has focussed so far on models of accretion onto "intermediate mass black holes", a possible class of black holes, examples of which may have recently been discovered in nearby galaxies. There, Krolik and his student (Yawei Hui) have computed stellar atmospheres for uniformly-heated disks around this class of black holes. Their models serve two purposes: they are the very first serious attempts to compute the spectrum from accreting black holes in this mass range; and a library of such models can be used later in this program as contrasts for those computed on the basis of real disk dynamics. The output from these local disk calculations has also been successfully coupled to a program that applies the appropriate relativistic transformations and computes photon trajectories in order to predict the spectrum received by observers located at different polar angles. The principal new result of these calculations is the discovery of potentially observable ionization edges of H-like C and O at frequencies near the peak in flux from these objects. Most of the grant money at UCSB was spent on supporting graduate student Shane Davis. In addition. some money was spent on supporting two other students: Ari Socrates (now a Hubble Fellow at Princeton), and Laura Melling. Davis spent the year constructing stellar atmosphere models of accretion disks appropriate for the high/soft (thermal) state of black hole X-ray binaries. As with AGN models published previously by our collaboration with NASA support. our models include a complete general relativistic treatment of both the disk structure and the propagation of photons from the disk to a distant observer. They also include all important continuum opacity sources, including Compton scattering and bound-free opacity from abundant metal species. The principal new result is that bound-free opacity is very significant in altering the continuum spectral shape, resulting for example in quite different "color correction factors" compared to those predicted previously. In addition, the models predict a relationship between luminosity and inner disk temperature that is, for the first time, in accord with that observed. The primary purpose of the grant was to incorporate more realistic accretion disk physics, learned largely from simulations, into such spectral models. The Davis et al. paper includes consideration of a vertical dissipation profile computed from radiation magneto-hydrodynamic simulations of MRI turbulence by N. J. Turner (2004). So long as the disk is effectively thick, such dissipation profiles do not affect the predicted spectrum significantly. (More work needs to be done on these simulations, however.) A potentially more serious issue is that MRI turbulence produces substantial inhomogeneities, as do photon bubble instabilities. These inhomogeneities can affect the spectra by enhancing the effects of absorption opacity over scattering opacity. We have done some preliminary Monte Carlo calculations to explore these effects.

Krolik, Julian H.↗

One-equation near-wall turbulence modeling with the aid of direct simulation data

The length scales appearing in the relations for the eddy viscosity and dissipation rate in one-equation models were evaluated from direct numerical simulation data for developed channel and boundary-layer flow at two Reynolds numbers each. To prepare the ground for the evaluation, the distribution of the most relevant mean-flow and turbulence quantities is presented and discussed with respect to Reynolds-number influence and to differences between channel and boundary-layer flow. An alternative model is also examined in which bar-(v'(exp 2))(exp 1/2) is used as velocity scale instead of k(exp 1/2). With this velocity scale, the length scales now appearing in the model follow very closely a linear relationship near the wall so that no damping is necessary. For the determination of bar-v'(exp 2) in the context of a one-equation model, a correlation is provided between bar-(v'(exp 2))/k and bar-(u'v')/k.

Rodi, W.↗

One-equation near-wall turbulence modeling with the aid of direct simulation data

The length scales appearing in the relations for the eddy viscosity and dissipation rate in one-equation models were evaluated from direct numerical (DNS) simulation data for developed channel and boundary-layer flow at two Reynolds numbers each. To prepare the ground for the evaluation, the distribution of the most relevant mean-flow and turbulence quantities is presented and discussed, also with respect to Reynolds-number influence and to differences between channel and boundary-layer flow. An alternative model is tested as near wall component of a two-layer model by application to developed-channel, boundary-layer and backward-facing-step flows.

Rodi, W.↗

Simulator Investigations of the Problems of Flying a Swept-Wing Transport Aircraft in Heavy Turbulence

An investigation of several factors which may contribute to the problem of piloting jet transport aircraft in heavy turbulence was conducted by using a piloted simulator that included the most significant airplane response and cockpit vibrations induced by rough air. Results indicated that the primary fuselage structural frequency contributed significantly to a distracting cockpit environment, and there was obtained evidence of severely reduced instrument flight proficiency during simulated maneuvering flight in heavy turbulence. It is concluded that the addition of similar rough-air response capabilities to training simulators would be of value in pilot indoctrination in turbulent-flight procedures.

Bray, Richard S.↗

A spatial model of wind shear and turbulence

The purpose of the spatial model considered in the present investigation is to generate the wind environment for use by others for flight simulation. Winds and gusts are provided over any finite area (e.g., aircraft body) from which aircraft loads and moments may be calculated. Three-dimensional autospectral information and correlation are contained in the data. It is pointed out that the three-dimensionality as contained in the spatial model affords much greater realism than widely used one-dimensional models. The resulting simulated wind is a nonlinear, non-Gaussian combination of real atmospheric winds and Gaussian, three-dimensional turbulence modulated by gust intensities which may vary freely as desired over space. The turbulence as represented by a product of a varying gust intensity and simulated turbulence is nonlinear and non-Gaussian.

Campbell, C. W.↗

On the active control of wall-bounded turbulent flows

The direct numerical simulation technique is used to explore concepts for manipulation of turbulent boundary layers. Significant drag reduction was achieved when at each instant the normal component of velocity at the wall was prescribed to be 180 deg out of phase with the normal velocity slightly above the wall. The drag reduction is 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 are examined.

Moin, P.↗

Numerical computation of transition to turbulence

Current applications of numerical simulations, including studies of mildly nonlinear and strongly nonlinear effects, and incipient turbulence are reviewed, and it is observed that one mildly nonlinear application for which direct simulations need to take a leading part is in establishing confidence in the dominant instability of a flow, while direct simulations are the primary tool when nonlinear effects are stronger than can be handled by existing asymptotic methods. The characteristics of existing algorithms for transition simulations are discussed, along with algorithmic improvements which would extend the range of problems amenable to numerical simulations. Current and future issues and problems are outlined.

Zang, Thomas A.↗

Uncertainty Propagation for Turbulent, Compressible Flow in a Quasi-1D Nozzle Using Stochastic Methods

This paper describes a fully spectral, Polynomial Chaos method for the propagation of uncertainty in numerical simulations of compressible, turbulent flow, as well as a novel stochastic collocation algorithm for the same application. The stochastic collocation method is key to the efficient use of stochastic methods on problems with complex nonlinearities, such as those associated with the turbulence model equations in compressible flow and for CFD schemes requiring solution of a Riemann problem. Both methods are applied to compressible flow in a quasi-one-dimensional nozzle. The stochastic collocation method is roughly an order of magnitude faster than the fully Galerkin Polynomial Chaos method on the inviscid problem.

Zang, Thomas A.↗

Wall-modeled LES of Flow Over a Gaussian Bump

We perform wall-modeled large eddy simulations (WMLES) of turbulent flow over a Gaussian-shaped bump geometry, to assess its performance in the acceleration and separation regions of the flow. An unstructured finite-volume solver is used along with an equilibrium wall model. The Mach number of the oncoming flow selected for the simulation is 0.176, and two Reynolds numbers (Re) are simulated that are about 10000 and 36000 based on boundary layer thickness upstream of the bump. Spanwise periodic results for the lower Re case are compared with the available DNS data, while the full three-dimensional simulation results for the higher Re case are compared to available experimental data. Sensitivity of WMLES results are assessed for a number of factors including grid resolution, wall model exchange location, type of wall model, and unstructured grid topology. The WMLES results for the lower Re case agree well with available DNS data in terms of the wall pressure variation, velocity and turbulent stress profile comparisons. The skin-friction predictions show a reduced tendency to separate. This appears to be accentuated for the higher Re case, for the medium resolution grid used in this study, where WMLES does not predict any flow separation in the mean, which is in contrast to the large separation bubble observed in the experiments.

turbulence↗

Base Heating Sensitivity Study for a 4-Cluster Rocket Motor Configuration in Supersonic Freestream

In support of launch vehicle base heating and pressure prediction efforts using the Loci-CHEM Navier-Stokes computational fluid dynamics solver, 35 numerical simulations of the NASA TND-1093 wind tunnel test have been modeled and analyzed. This test article is composed of four JP-4/LOX 500 lbf rocket motors exhausting into a Mach 2 - 3.5 wind tunnel at various ambient pressure conditions. These water-cooled motors are attached to a base plate of a standard missile forebody. We explore the base heating profiles for fully coupled finite-rate chemistry simulations, one-way coupled RAMP (Reacting And Multiphase Program using Method of Characteristics)-BLIMPJ (Boundary Layer Integral Matrix Program - Jet Version) derived solutions and variable and constant specific heat ratio frozen flow simulations. Variations in turbulence models, temperature boundary conditions and thermodynamic properties of the plume have been investigated at two ambient pressure conditions: 255 lb/sq ft (simulated low altitude) and 35 lb/sq ft (simulated high altitude). It is observed that the convective base heat flux and base temperature are most sensitive to the nozzle inner wall thermal boundary layer profile which is dependent on the wall temperature, boundary layer s specific energy and chemical reactions. Recovery shock dynamics and afterburning significantly influences convective base heating. Turbulence models and external nozzle wall thermal boundary layer profiles show less sensitivity to base heating characteristics. Base heating rates are validated for the highest fidelity solutions which show an agreement within +/-10% with respect to test data.

Mehta, Manish↗

Effects of turbulence mixing, variable properties, and vaporization on spray droplet combustion

Combustion of liquid fuels in the form of spray droplets is simulated numerically. Various vaporization models are examined as to their performance in finite element calculations involving a turbulent flow field. The Eulerian coordinate for the gas and Lagrangian coordinate for the liquid spray droplets are coupled through source terms being updated in the equations of continuity, momentum, and energy. The k-epsilon and modified eddy breakup models are used for simulating turbulent spray combustion flow field. Numerical results for the droplet trajectories, droplet heating, recirculation characteristics, and effects of evaporation models are evaluated. It is also shown that the finite element method is advantageous in dealing with complex geometries, complex boundary conditions, adaptive unstructured grids.

Kim, Y. M.↗

Velocity and pressure fields associated with near-wall turbulence structures

Computer generated databases containing velocity and pressure fields in three-dimensional space at a sequence of time-steps, were used for the investigation of near-wall turbulence structures, their space-time evolution, and their associated pressure fields. The main body of the results were obtained from simulation data for turbulent channel flow at a Reynolds number of 180 (based on half-channel height and friction velocity) with a grid of 128 x 129 x and 128 points. The flow was followed over a total time of 141 viscous time units. Spanwise centering of the detected structures was found to be essential in order to obtain a correct magnitude of the associated Reynolds stress contribution. A positive wall-pressure peak is found immediately beneath the center of the structure. The maximum amplitude of the pressure pattern was, however, found in the buffer region at the center of the shear-layer. It was also found that these flow structures often reach a maximum strength in connection with an asymmetric spanwise motion, which motivated the construction of a conditional sampling scheme that preserved this asymmetry.

Johansson, Arne V.↗

Conditional statistics in a turbulent premixed flame derived from direct numerical simulation

The objective of this paper is to briefly introduce conditional moment closure (CMC) methods for premixed systems and to derive the transport equation for the conditional species mass fraction conditioned on the progress variable based on the enthalpy. Our statistical analysis will be based on the 3-D DNS database of Trouve and Poinsot available at the Center for Turbulence Research. The initial conditions and characteristics (turbulence, thermo-diffusive properties) as well as the numerical method utilized in the DNS of Trouve and Poinsot are presented, and some details concerning our statistical analysis are also given. From the analysis of DNS results, the effects of the position in the flame brush, of the Damkoehler and Lewis numbers on the conditional mean scalar dissipation, and conditional mean velocity are presented and discussed. Information concerning unconditional turbulent fluxes are also presented. The anomaly found in previous studies of counter-gradient diffusion for the turbulent flux of the progress variable is investigated.

Mantel, Thierry↗

TPSAS-NF1676L-10767-DND

Version 2.1ae of the NASA CFD code OVERFLOW has been used to simulate the flow past tandem cylinders. The multiblock grid system consisted of 15 overset meshes with a total grid count of 15.7 million points. The grid extended 3 cylinder diameters in the spanwise direction with periodic boundary conditions at the ends. Riemann conditions were used at the outer boundaries, and a no-slip condition on the cylinder surfaces. The 3rd order HLLC upwind spatial discretization scheme was used with a dual-time stepping algorithm. A modified version of the Delayed Detached Eddy Simulation hybrid RANS/LES turbulence model proposed by Spalart was used in a fully turbulent mode. The code was run on 53 Core 2 Duo nodes using MPI across nodes and OPENMP within each node. Simulation results will be compared with the available experimental data.

David P. Lockard↗

Numerical Simulation of High Drag Reduction in a Turbulent Channel Flow with Polymer Additives

The addition of small amounts of long chain polymer molecules to wall-bounded flows can lead to dramatic drag reduction. Although this phenomenon has been known for about fifty years, the action of the polymers and its effect on turbulent structures are still unclear. Detailed experiments have characterized two distinct regimes (Warholic et al. 1999), which are referred to as low drag reduction (LDR) and high drag reduction (HDR). The first regime exhibits similar statistical trends as Newtonian flow: the log-law region of the mean velocity profile remains parallel to that of the Newtonian ow but its lower bound moves away from the wall and the upward shift of the log-region is a function of drag reduction, DR. Although streamwise fluctuations are increased and transverse ones are reduced, the shape of the rms velocity profiles is not qualitatively modified. At higher drag reductions, of the order of 40-50%, the ow enters the HDR regime for which the slope of the log-law is dramatically augmented and the Reynolds shear stress is small (Warholic et al. 1999; Ptasinski et al. 2001). The drag reduction is eventually bounded by a maximum drag reduction (MDR) (Virk & Mickley 1970) which is a function of the Reynolds number. While several experiments report mean velocity profiles very close to the empirical profile of Virk & Mickley (1970) for MDR conditions, the observations regarding the structure of turbulence can differ significantly. For instance, Warholic et al. (1999) measured a near-zero Reynolds shear stress, whereas a recent experiment (Ptasinski et al. 2001) shows evidence of non-negligible Reynolds stress in their MDR flow. To the knowledge of the authors, only the LDR regime has been documented in numerical simulations (Sureshkumar et al. 1997; Dimitropoulos et al. 1998; Min et al. 2001; Dubief & Lele 2001; Sibilla & Baron 2002). This paper discusses the simulation of polymer drag reduced channel ow at HDR using the FENE-P (Finite Elastic non-linear extensibility-Peterlin) model which was used for the first LDR simulation by Sureshkumar et al. (1997). Flow and polymer parameters are close to realistic polymer drag reducing conditions. High drag reductions are achieved by using finite differences and a robust time stepping technique. A minimal channel flow is also used as a numerical experiment to investigate the effect of the outer region turbulent structures on the overall drag at HDR. The drag reducing action of the model is finally studied through the structure of energy transfers from the polymers to the velocity components. This investigation sheds some light on the details of polymer drag reduction.

Dubief, Yves↗