Search NASA⌕ Search

Engineering topics

Ferziger, J. H.

Publications and source records attributed to Ferziger, J. H..

At least 19 records
All sourcesNASA NTRSDOE OSTI
NASA NTRS2003

Testing of RANS Turbulence Models for Stratified Flows Based on DNS Data

In most geophysical flows, turbulence occurs at the smallest scales and one of the two most important additional physical phenomena to account for is strati cation (the other being rotation). In this paper, the main objective is to investigate proposed changes to RANS turbulence models which include the effects of stratifi- cation more explicitly. These proposed changes were developed using a DNS database on strati ed and sheared homogenous turbulence developed by Shih et al. (2000) and are described more fully in Ferziger et al. (2003). The data generated by Shih, et al. (2000) (hereinafter referred to as SKFR) are used to study the parameters in the k- model as a function of the turbulent Froude number, Frk. A modified version of the standard k- model based on the local turbulent Froude number is proposed. The proposed model is applied to a stratified open channel flow, a test case that differs significantly from the flows from which the modified parameters were derived. The turbulence modeling and results are discussed in the next two sections followed by suggestions for future work.

Venayagamoorthy, S. K.↗
NASA NTRS2003

Numerical Simulation of Regional Circulation in the Monterey Bay Region

The objective of this study is to produce a high-resolution numerical model of Mon- terey Bay area in which the dynamics are determined by the complex geometry of the coastline, steep bathymetry, and the in uence of the water masses that constitute the CCS. Our goal is to simulate the regional-scale ocean response with realistic dynamics (annual cycle), forcing, and domain. In particular, we focus on non-hydrostatic e ects (by comparing the results of hydrostatic and non-hydrostatic models) and the role of complex geometry, i.e. the bay and submarine canyon, on the nearshore circulation. To the best of our knowledge, the current study is the rst to simulate the regional circulation in the vicinity of Monterey Bay using a non-hydrostatic model. Section 2 introduces the high resolution Monterey Bay area regional model (MBARM). Section 3 provides the results and veri cation with mooring and satellite data. Section 4 compares the results of hydrostatic and non-hydrostatic models.

Tseng, Y. H.↗
NASA NTRS2002

Waves in Turbulent Stably-Stratified Shear Flow

Two approaches for the identification of internal gravity waves in sheared and unsheared homogeneous stratified turbulence are investigated. First, the phase angle between the vertical velocity and density fluctuations is considered. It is found, however, that a continuous distribution of the phase angle is present in weakly and strongly stratified flow. Second, a projection onto the solution of the linearized inviscid equations of motion of unsheared stratified flow is investigated. It is found that a solution of the fully nonlinear viscous Navier-Stokes equations can be represented by the linearized inviscid solution. The projection yields a decomposition into vertical wave modes and horizontal vortical modes.

Jacobitz, F. G.↗
NASA NTRS2002

Waves in Turbulent Stably Stratified Shear Flow

Two approaches for the identification of internal gravity waves in sheared and unsheared homogeneous stratified turbulence are investigated. First, the phase angle between the vertical velocity and density fluctuations is considered. It was found, however, that a continuous distribution of the phase angle is present in weakly and strongly stratified flow. Second, a projection onto the solution of the linearized inviscid equations of motion of unsheared stratified flow is investigated. It was found that a solution of the fully nonlinear viscous Navier-Stokes equations can be represented by the linearized inviscid solution. The projection yields a decomposition into vertical wave modes and horizontal vortical modes.

Jacobitz, F. G.↗
NASA NTRS1996

Dynamic Models for LES of Turbulent Front Propagation With a Spectral Method

Direct numerical simulation of turbulent reacting flows places extreme demands on computational resources. At the present time, simulations can be performed only for greatly simplified reaction systems and for very low Reynolds numbers. Direct simulation of more realistic cases occurring at higher Reynolds number and including multiple species and numerous chemical reactions will exceed available computational resources far into the future. Because of this, there is a clear need to develop the technique of large eddy simulation for reacting flows. Unfortunately this task is complicated by the fact that combustion arises from chemical reactions that occur at the smallest scales of the flow. Capturing the large-scale behavior without resolving the small-scale details is extremely difficult in combustion problems. Thus LES modeling for turbulent combustion encounters difficulties not present in modeling momentum transport, in which the main effect of the small scales is to provide dissipation. The difficulty is more pronounced in premixed combustion, where detailed chemistry plays an essential role in determining the flame speed (or overall burning rate); in nonpremixed combustion infinite rate chemistry can be assumed, eliminating the small scale features to a first approximation.

Im, H. G.↗
NASA NTRS1996

A dynamic subgrid-scale model for LES of the G-equation

Turbulent combustion is a difficult subject as it must deal with all of the issues found in both turbulence and combustion. (We consider only premixed flames in this paper, but some of the ideas can be applied to the non-premixed case.) As in many other fields, there are two limiting cases that are easier to deal with than the general case. These are the situations in which the chemical time scale is either much shorter or much longer than the time scale associated with the turbulence. We deal with the former case. In this limit, the flame is thin compared to the turbulence length scales and can be idealized as an infinitely thin sheet. This is commonly called the flamelet regime; it has been the subject of many papers and the basis for many models (see, e.g., Linan & Williams 1993). In the flamelet model, the local flame structure is assumed to be identical to the laminar flame structure; thus the flame propagates normal to itself at the laminar flame speed, S(sub L). This allows the use of simple approximations. For example, one expects the rate of consumption of fuel to be proportional to the area of the flame surface. This idea allowed Damkohler (1940) to propose that the wrinkled flame could be replaced by a smooth one which travels at the turbulent flame speed, S(sub T), defined by S(sub T)/S(sub L) = A(sub L) /A(sub P) where A(sub L) is the total flame surface area and AP is the area projected onto the mean direction of propagation. This relation can be expected to be valid when the flame structure is modified only slightly by the turbulence. More recent approaches have attempted to relate the turbulent flame speed to turbulence intensity, u(sub '), which presumably, characterizes the wrinkling of the flame.

Bourlioux, A.↗
NASA NTRS1995

A new non-eddy viscosity subgrid-scale model and its application to channel flow

To date, most large-eddy simulations (LES) have been carried out with eddy viscosity subgrid scale (SGS) models, with only a few exceptions that used the mixed model. Even though the assumptions behind Smagorinsky's model are rather stringent, it has been applied successfully to a variety of turbulent flows. This success is attributed to the ability of eddy viscosity models to drain energy from large scales, thus simulating the dissipative nature of turbulence. Most SGS models are absolutely dissipative, i.e. they remove energy from the large scales at every instant. However, SGS stresses may transfer energy back to the large scales intermittently; this reverse transfer or backscatter is especially important in geophysical flows and in transition. In a fully developed channel flow, there is reverse flow of energy from small to large scales near the walls, but eddy viscosity models are unable to account for this important feature. The dynamic localization eddy viscosity model of Ghosal et al. (1995) allows backscatter by co-evolving an auxiliary equation for the SGS energy; however, the computational cost is considerably larger than for conventional SGS models (Cabot 1994). In this report, a new non-eddy viscosity model based on local approximation of total quantities in terms of filtered ones is introduced; the scale similarity model of Bardina (1983) is a special case of this model. This procedure does not require the assumption of homogeneity, permits backscatter of energy from small to large scales, and is readily implemented in finite difference codes. The results of applying the proposed model to second order finite volume simulation of plane channel flow at high Reynolds numbers (Re(sub b) = 38,000) is described in this report. Greater emphasis is placed on the high Reynolds number flow since it provides a more rigorous test of the SGS model and its potential application. The results are compared to ones produced by the conventional and dynamic Smagorinsky models and the spectral LES of Piomelli (1993).

Shah, K. B.↗
NASA NTRS1994

Simulation of a turbulent flame in a channel

The interaction between turbulent premixed flames and channel walls is studied. Combustion is represented by a simple irreversible reaction with a large activation temperature. Feedback to the flowfield is suppressed by invoking a constant density assumption. The effect of wall distance on local and global flame structure is investigated. Quenching distances and maximum wall heat fluxes computed in laminar cases are compared to DNS results. It is found that quenching distances decrease and maximum heat fluxes increase relative to laminar flame values. It is shown that these effects are due to large coherent structures which push flame elements towards to wall. The effect of wall strain is studied in flame-wall interaction in a stagnation line flow; this is used to explain the DNS results. It is also shown that 'remarkable' flame events are produced by interaction with a horseshoe vortex: burnt gases are pushed towards the wall at high speed and induce quenching and high wall heat fluxes while fresh gases are expelled from the wall region and form finger-like structures. Effects of the wall on flame surface density are investigated, and a simple model for flame-wall interaction is proposed; its predictions compare well with the DNS results.

Bruneaux, G.↗
NASA NTRS1993

Linear stability of the confined compressible reacting mixing layer

This paper investigates the linear stability of confined mixing layers with special emphasis on the effects of heat release and compressibility. The results show that reflection of supersonic disturbances by the walls makes the confined supersonic mixing layer more unstable than the unconfined free shear layer. Decreasing the distance between the walls makes the flow more unstable. However, subsonic disturbances are relatively unaffected by the walls. Heat release and Mach number hardly change the growth rates of supersonic disturbances. The most unstable supersonic disturbances are two-dimensional in rectangular channel flows, but three-dimensional in partially confined flows. Finally, the reactants are not strongly mixed by supersonic instabilities, which mainly disturb one side of the layer.

Shin, D. S.↗
NASA NTRS1993

Mixing in a stratified shear flow: Energetics and sampling

Direct numerical simulations of the time evolution of homogeneous stably stratified shear flows have been performed for Richardson numbers from 0 to 1 and for Prandtl numbers between 0.1 and 2. The results indicate that mixing efficiency R(sub f) varies with turbulent Froude number in a manner consistent with laboratory experiments performed with Prandtl numbers of 0.7 and 700. However, unlike the laboratory results, for a particular Froude number, the simulations do not show a clear dependence on the magnitude of R(sub f) on Pr. The observed maximum value of R(sub f) is 0.25. When averaged over vertical length scales of an order of magnitude greater than either the overturning or Ozmidov scales of the flow, the simulations indicate that the dissipation rate epsilon is only weakly lognormally distributed with an intermittency of about 0.01 whereas estimated values in the ocean are 3 to 7.

Ivey, G. N.↗
NASA NTRS1992

Direct simulation of the stably stratified turbulent Ekman layer

The Navier-Stokes equations and the Boussinesq approximation were used to compute a 3D time-dependent turbulent flow in the stably stratified Ekman layer over a smooth surface. The simulation data are found to be in very good agreement with atmospheric measurements when nondimensionalized according to Nieuwstadt's local scaling scheme. Results suggest that, when Reynolds number effects are taken into account, the 'constant Froud number' stable layer model (Brost and Wyngaard, 1978) and the 'shearing length' stable layer model (Hunt, 1985) for the dissipitation rate of turbulent kinetic energy are both valid. It is concluded that there is good agreement between the direct numerical simulation results and large-eddy simulation results obtained by Mason and Derbyshire (1990).

Coleman, G. N.↗
NASA NTRS1991

Stability of low-speed reacting flows

Small perturbation, linear, inviscid stability theory is developed and applied to coflowing chemically reacting jets. A necessary condition for the existence of axisymmetric and azimuthal instabilities is derived. The eigenfunction structure provides considerable insight into the underlying physical mechanism. In nonpremixed cases, high temperature tends to occur in regions where the shear layer tries to roll up; this reduces the available momentum, suppressing the growth of instabilties. It is found that chemical reaction during the growth of instabilities is unimportant, suggesting that this effect could be ignored in a viscous theory. Both nonpremixed and premixed flames in axisymmetric shear layers are generally more stable to small perturbations than in the corresponding cold case. However, while for nonpremixed flames the frequency of the most amplified wave decreases with increasing heat release, the opposite occurs for premixed cases.

Mahalingham, S.↗
NASA NTRS1991

Methods for direct simulation of transition in hypersonic boundary layers

An implicit numerical algorithm for the time accurate solution of the compressible Navier-Stokes equations is described. Results for steady flow past a finite flat plate are presented, together with preliminary results for the temporal simulation of second mode instability in a flat plate boundary layer at Mach 4.5.

Vandervegt, J. J. W.↗
NASA NTRS1991

Effect of heat release on the stability of compressible reacting mixing layer

Reacting free shear layers are of fundamental importance in many industrial systems including gas turbine combustors and rockets. Efficient propulsion systems are essential for air breathing supersonic ramjets in the high Mach number range. A limiting factor in these engines is the time for fuel and oxidizer to mix in the combustion chamber; for fast mixing, the flow must be vigorously turbulent which requires the laminar flow to be unstable. Understanding the stability characteristics of compressible reacting free shear layers is, therefore, very important and may allow one to control the flow. Low speed shear layers are highly unstable but, as chemical reaction and compressibility effects tend to stabilize them, it is important to investigate the stability of high speed reacting mixing layers. The latter consists of two fluid streams containing fuel and oxidizer respectively, and the conclusions are expected to apply, with quantitative modifications, to other shear flows, e.g., jets. Since low speed reacting cases have been studied earlier, we concentrate on the effects of Mach number and heat release. We are primarily interested in solving the stability problem over a large range of Mach number and heat release. In order to understand the effect of the heat release on the stability of this flow, one must first study the characteristics of the non-reacting flow. Inviscid theory is a reliable guide for understanding stability of compressible shear flows at moderate and large Reynolds numbers and is the basis for this work.

Shin, Dong-Shin↗
NASA NTRS1991

Stability of compressible reacting mixing layer

Linear instability of compressible reacting mixing layers is analyzed with emphasis on the effects of heat release and compressibility. Laminar solutions of the compressible boundary-layer equations are used as the base flows. The parameters of this study are the adiabatic flame temperature, the Mach number of the upper stream, frequency, wavenumber, and the direction of propagation of the disturbance wave. Stability characteristics of the flow are presented. Three groups of unstable modes are found when the Mach number and/or heat release are large. Finally, it is shown that the unstable modes are two-dimensional for large heat release even in highly compressible flow.

Shin, D. S.↗
NASA NTRS1990

Direct numerical simulations of stably-stratified sheared turbulence: Implications for oceanic mixing

Direct numerical simulations of the time evolution of homogeneous stably stratified turbulent shear flows have been performed for several Richardson numbers Ri and Reynolds numbers R(sub lambda) in earlier works. The results show excellent agreement with length scale models developed from laboratory experiments to characterize oceanic turbulence. When the Richardson number Ri is less than the stationary value Ri(sub s), the turbulence intensity grows at all scales, and the growth rate appears to be a function of Ri. The size of the vertical density inversions also increases. On the other hand, when Ri is greater than or equal to Ri(sub s) the largest turbulent eddies become vertically constrained by buoyancy when the Ellison (turbulence) scale L(sub E) and the Ozmidov (buoyancy) scale L(sub O) are equal. At this point, the mixing efficiency is maximal and corresponds to a flux Richardson number R(sub f) = 0.20. The vertical mass flux becomes counter-gradient when epsilon = 19(nu)N(exp 2) and vertical density overturns are suppressed in less than half a Brunt-Vaisala period. The results of the simulations were also recast in terms of the Hydrodynamic Phase Diagram introduced in fossil turbulence models. The so-called point of fossilization occurs when epsilon = 4DCN(exp 2); Gibson proposed 13DCN(exp 2). This value is in agreement with indirect laboratory observations and field observations. Finally, the validity of the steady-state models to estimate vertical eddy diffusivities in the oceanic thermocline is discussed.

Itsweire, E. C.↗