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 1,063 records · Page 59

Advances in modeling the pressure correlation terms in the second moment equations

In developing turbulence models, various model constraints were proposed in an attempt to make the model equations more general (or universal). The most recent of these are the realizability principle, the linearity principle, the rapid distortion theory, and the material indifference principle. Several issues are discussed concerning these principles and special attention is payed to the realizability principle. Realizability (defined as the requirement of non-negative energy and Schwarz' inequality between any fluctuating quantities) is the basic physical and mathematical principle that any modeled equation should obey. Hence, it is the most universal, important and also the minimal requirement for a model equation to prevent it from producing unphysical results. The principle of realizability is described in detail, the realizability conditions are derived for various turbulence models, and the model forms are proposed for the pressure correlation terms in the second moment equations. Detailed comparisons of various turbulence models with experiments and direct numerical simulations are presented. As a special case of turbulence, the two dimensional two-component turbulence modeling is also discussed.

Shih, Tsan-Hsing↗

Advances in modeling the pressure correlation terms in the second moment equations

In developing turbulence models, various model constraints were proposed in an attempt to make the model equations more general (or universal). The most recent of these are the realizability principle, the linearity principle, the rapid distortion theory, and the material indifference principle. Several issues are discussed concerning these principles and special attention is payed to the realizability principle. Realizability (defined as the requirement of non-negative energy and Schwarz' inequality between any fluctuating quantities) is the basic physical and mathematical principle that any modeled equation should obey. Hence, it is the most universal, important and also the minimal requirement for a model equation to prevent it from producing unphysical results. The principle of realizability is described in detail, the realizability conditions are derived for various turbulence models, and the model forms are proposed for the pressure correlation terms in the second moment equations. Detailed comparisons of various turbulence models with experiments and direct numerical simulations are presented. As a special case of turbulence, the two dimensional two-component turbulence modeling is also discussed.

Shih, Tsan-Hsing↗

Effects of Gravity on Sheared Turbulence Laden with Bubbles or Droplets

The objective of this numerical/experimental study is to improve the understanding of the effects of gravity on the two-way interaction between dispersed particles (bubbles or liquid droplets) and the carrier turbulent flow. The first phase of the project considers isotropic turbulence. Turbulent homogeneous shear flows laden with droplets/bubbles will be studied in the next phase. The experiments reported here are concerned with the dispersion of liquid droplets by homogeneous turbulence under various gravitational conditions and the effect of these droplets on the evolution of the turbulence of the carrier fluid (air). Direct numerical simulations (DNS) of bubble - laden isotropic decaying turbulence are performed using the two-fluid approach (TF) instead of the Eulerian-Lagrangian approach (EL). The motivation for using the TF formulation is that EL requires considerable computational resources especially for the case of two-way coupling where the instantaneous trajectories of a large number of individual bubbles need to be computed. The TF formulation is developed by spatially averaging the instantaneous equations of the carrier flow and bubble phase over a scale of the order of the Kolmogorov length scale which, in our case, is much larger than the bubble diameter. On that scale, the bubbles are treated as a continuum (without molecular diffusivity) characterized by the bubble phase velocity field and concentration (volume fraction). The bubble concentration, C, is assumed small enough to neglect the bubble-bubble interactions.

Elghobashi, Said↗

DNS of Laminar-Turbulent Transition in Swept-Wing Boundary Layers

Direct numerical simulation (DNS) is performed to examine laminar to turbulent transition due to high-frequency secondary instability of stationary crossflow vortices in a subsonic swept-wing boundary layer for a realistic natural-laminar-flow airfoil configuration. The secondary instability is introduced via inflow forcing and the mode selected for forcing corresponds to the most amplified secondary instability mode that, in this case, derives a majority of its growth from energy production mechanisms associated with the wall-normal shear of the stationary basic state. An inlet boundary condition is carefully designed to allow for accurate injection of instability wave modes and minimize acoustic reflections at numerical boundaries. Nonlinear parabolized stability equation (PSE) predictions compare well with the DNS in terms of modal amplitudes and modal shape during the strongly nonlinear phase of the secondary instability mode. During the transition process, the skin friction coefficient rises rather rapidly and the wall-shear distribution shows a sawtooth pattern that is analogous to the previously documented surface flow visualizations of transition due to stationary crossflow instability. Fully turbulent features are observed in the downstream region of the flow.

Duan, L.↗

Spectral multigrid methods for the solution of homogeneous turbulence problems

New three-dimensional spectral multigrid algorithms are analyzed and implemented to solve the variable coefficient Helmholtz equation. Periodicity is assumed in all three directions which leads to a Fourier collocation representation. Convergence rates are theoretically predicted and confirmed through numerical tests. Residual averaging results in a spectral radius of 0.2 for the variable coefficient Poisson equation. In general, non-stationary Richardson must be used for the Helmholtz equation. The algorithms developed are applied to the large-eddy simulation of incompressible isotropic turbulence.

Erlebacher, G.↗

Spectral multigrid methods for the solution of homogeneous turbulence problems

New three-dimensional spectral multigrid algorithms are analyzed and implemented to solve the variable coefficient Helmholtz equation. Periodicity is assumed in all three directions which leads to a Fourier collocation representation. Convergence rates are theoretically predicted and confirmed through numerical tests. Residual averaging results in a spectral radius of 0.2 for the variable coefficient Poisson equation. In general, non-stationary Richardson must be used for the Helmholtz equation. The algorithms developed are applied to the large-eddy simulation of incompressible isotropic turbulence.

Erlebacher, G.↗

A dynamic subgrid-scale eddy viscosity model

One major drawback of the eddy viscosity subgrid-scale stress models used in large-eddy simulations is their inability to represent correctly with a single universal constant different turbulent field in rotating or sheared flows, near solid walls, or in transitional regimes. In the present work, a new eddy viscosity model is presented which alleviates many of these drawbacks. The model coefficient is computed dynamically as the calculation progresses rather than input a priori. The model is based on an algebraic identity (Germano 1990) between the subgrid-scale stresses at two different filtered levels and the resolved turbulent stresses. The subgrid-scale stresses obtained using the proposed model vanish in laminar flow and at a solid boundary, and have the correct asymptotic behavior in the near-wall region of a turbulent boundary layer. The results of large-eddy simulations of transitional and turbulent channel flow that use the proposed model are in good agreement with the direct simulation data.

Germano, Massimo↗

Chemical reacting flows

Future aerospace propulsion concepts involve the combination of liquid or gaseous fuels in a highly turbulent internal air stream. Accurate predictive computer codes which can simulate the fluid mechanics, chemistry, and turbulence combustion interaction of these chemical reacting flows will be a new tool that is needed in the design of these future propulsion concepts. Experimental and code development research is being performed at Lewis to better understand chemical reacting flows with the long term goal of establishing these reliable computer codes. The approach to understanding chemical reacting flows is to look at separate simple parts of this complex phenomena as well as to study the full turbulent reacting flow process. As a result research on the fluid mechanics associated with chemical reacting flows was initiated. The chemistry of fuel-air combustion is also being studied. Finally, the phenomena of turbulence-combustion interaction is being investigated. This presentation will highlight research, both experimental and analytical, in each of these three major areas.

Mularz, Edward J.↗

Turbulence modeling for high speed compressible flows

Turbulence models are necessary in numerical simulations because of the impracticality of computing all scales of turbulent motion. Approximate methods, or turbulence models, are introduced to simplify and make the computations practical. There are several approaches to turbulence modeling: subgrid scale modeling; modeling of all the scales of turbulent motion; eddy viscosity; and Reynolds stress turbulence closure models. In many instances, computations based on turbulence models are compared with available experimental data. The effort is directed at seeking models which have reasonable accuracy over a limited range of flow conditions. Future needs are briefly discussed. Work to date is outlined.

Chandra, Suresh↗

Wind models for flight simulator certification of landing and approach guidance and control systems

The definition of a model suitable for certification was the main objective of this report. The model was designed to simplify and reduce the wind model parameters to enable evaluation of a large number of aircraft and control system design parameters. Analytical descriptions of wind phenomena were presented. For those parameters defying analytic description, probabilistic descriptions were sought. A brief analysis of the effects of wind on aircraft motion was conducted. The axes transformations required between wind and turbulence components in their inherent axis system and in the airplane's axis system were shown. Techniques of providing a random process on computers for the representation of turbulence were presented. A simulation model was presented that combines all the foregoing components.

Schaeffer, D. R.↗

The effects of particle loading on turbulence structure and modelling

The objective of the present research was to extend the Direct Numerical Simulation (DNS) approach to particle-laden turbulent flows using a simple model of particle/flow interaction. The program addressed the simplest type of flow, homogeneous, isotropic turbulence, and examined interactions between the particles and gas phase turbulence. The specific range of problems examined include those in which the particle is much smaller than the smallest length scales of the turbulence yet heavy enough to slip relative to the flow. The particle mass loading is large enough to have a significant impact on the turbulence, while the volume loading was small enough such that particle-particle interactions could be neglected. Therefore, these simulations are relevant to practical problems involving small, dense particles conveyed by turbulent gas flows at moderate loadings. A sample of the results illustrating modifications of the particle concentration field caused by the turbulence structure is presented and attenuation of turbulence by the particle cloud is also illustrated.

Squires, Kyle D.↗

Analysis of aerodynamic coefficients using gradient data: Spanwise turbulence effects on airplane response

The influence of spanwise turbulence on airplane dynamic behavior is determined. Calculations are based on data collected from test flights with the NASA B-57 research aircraft. The approach is to first compute aerodynamic forces and moments due to a spanwise distribution of angle of attack and airspeed. Secondly, these quantities are incorporated into the equations of motion. Simulation of flights done with the effects of spanwise turbulence included are compared to simulations without any spanwise turbulence. The findings of the study are that the moments developed by turbulence along the span are significant and that more realistic flight simulation can be achieved by including the spanwise turbulence terms.

Ringnes, E. A.↗

Turbulence modeling: Near-wall turbulence and effects of rotation on turbulence

Many Reynolds averaged Navier-Stokes solvers use closure models in conjunction with 'the law of the wall', rather than deal with a thin, viscous sublayer near the wall. This work is motivated by the need for better models to compute near wall turbulent flow. The authors use direct numerical simulation of fully developed channel flow and one of three dimensional turbulent boundary layer flow to develop new models. These direct numerical simulations provide detailed data that experimentalists have not been able to measure directly. Another objective of the work is to examine analytically the effects of rotation on turbulence, using Rapid Distortion Theory (RDT). This work was motivated by the observation that the pressure strain models in all current second order closure models are unable to predict the effects of rotation on turbulence.

Shih, T.-H.↗

Some Results Relevant to Statistical Closures for Compressible Turbulence

For weakly compressible turbulent fluctuations there exists a small parameter, the square of the fluctuating Mach number, that allows an investigation using a perturbative treatment. The consequences of such a perturbative analysis in three different subject areas are described: 1) initial conditions in direct numerical simulations, 2) an explanation for the oscillations seen in the compressible pressure in the direct numerical simulations of homogeneous shear, and 3) for turbulence closures accounting for the compressibility of velocity fluctuations. Initial conditions consistent with small turbulent Mach number asymptotics are constructed. The importance of consistent initial conditions in the direct numerical simulation of compressible turbulence is dramatically illustrated: spurious oscillations associated with inconsistent initial conditions are avoided, and the fluctuating dilatational field is some two orders of magnitude smaller for a compressible isotropic turbulence. For the isotropic decay it is shown that the choice of initial conditions can change the scaling law for the compressible dissipation. A two-time expansion of the Navier-Stokes equations is used to distinguish compressible acoustic and compressible advective modes. A simple conceptual model for weakly compressible turbulence - a forced linear oscillator is described. It is shown that the evolution equations for the compressible portions of turbulence can be understood as a forced wave equation with refraction. Acoustic modes of the flow can be amplified by refraction and are able to manifest themselves in large fluctuations of the compressible pressure.

Ristorcelli, J. R.↗

Simulating transitional flow and heat transfer over the flat plate and circular cylinder using a K-epsilon turbulence model

A boundary layer computer code, called TEXSTAN, has been used to simulate external, transitional, fluid flow and heat transfer over the flat plate and circular cylinder. Turbulence modeling was accomplished by using a low-Reynolds number K-epsilon turbulence model. The study was made for free-stream turbulence intensities less than 15 percent. The comparison of the flat plate heat transfer predictions to experimental data showed good agreement, for the location of transition as well as in the heat transfer coefficient. The typical discrepancy in predicting the heat transfer coefficient was less than 5 percent of the measured value.

Sullivan, T. J.↗

The NASA-Langley Wake Vortex Modelling Effort in Support of an Operational Aircraft Spacing System

Two numerical modelling efforts, one using a large eddy simulation model and the other a numerical weather prediction model, are underway in support of NASA's Terminal Area Productivity program. The large-eddy simulation model (LES) has a meteorological framework and permits the interaction of wake vortices with environments characterized by crosswind shear, stratification, humidity, and atmospheric turbulence. Results from the numerical simulations are being used to assist in the development of algorithms for an operational wake-vortex aircraft spacing system. A mesoscale weather forecast model is being adapted for providing operational forecast of winds, temperature, and turbulence parameters to be used in the terminal area. This paper describes the goals and modelling approach, as well as achievements obtained to date. Simulation results will be presented from the LES model for both two and three dimensions. The 2-D model is found to be generally valid for studying wake vortex transport, while the 3-D approach is necessary for realistic treatment of decay via interaction of wake vortices and atmospheric boundary layer turbulence. Meteorology is shown to have an important affect on vortex transport and decay. Presented are results showing that wake vortex transport is unaffected by uniform fog or rain, but wake vortex transport can be strongly affected by nonlinear vertical change in the ambient crosswind. Both simulation and observations show that atmospheric vortices decay from the outside with minimal expansion of the core. Vortex decay and the onset three-dimensional instabilities are found to be enhanced by the presence of ambient turbulence.

Proctor, Fred H.↗

A New Recycling Method to Generate Turbulent Inflow Profiles

The accuracy of the scale-resolving simulations for practical geometries strongly depends on the inflow boundary conditions. Imposing experimentally observed turbulent inflow profiles for the numerical simulations is a major challenge. Existing methods available in the literature assume self-similar behavior, which is not true for most of the experiments. In the present work, we formulate the turbulent inflow profile generation technique as an optimization problem. An adjoint technique is exploited to evaluate the sensitivities of multiple input parameters for the present problem. The present formulation is then tested to generate a laminar boundary layer profile, turbulent boundary layer profile, and turbulent jet profile.

Garai, Anirban↗

Consistent Initial Conditions for the DNS of Compressible Turbulence

Relationships between diverse thermodynamic quantities appropriate to weakly compressible turbulence are derived. It is shown that for turbulence of a finite turbulent Mach number there is a finite element of compressibility. A methodology for generating initial conditions for the fluctuating pressure, density and dilatational velocity is given which is consistent with finite Mach number effects. Use of these initial conditions gives rise to a smooth development of the flow, in contrast to cases in which these fields are specified arbitrarily or set to zero. Comparisons of the effect of different types of initial conditions are made using direct numerical simulation of decaying isotropic turbulence.

Ristorcelli, J. R.↗