9512.net

# Wavelet Cross-Correlation Analysis of Turbulent Mixing from Large-Eddy-Simulations

arXiv:physics/0003029v1 [physics.flu-dyn] 14 Mar 2000

Wavelet Cross-Correlation Analysis of Turbulent Mixing from Large-Eddy-Simulations
S. Sello and J. Bellazzini
Enel Research A. Pisano 120, Pisa 56122, ITALY

Contact e-mail: sello@pte.enel.it

1

Introduction

The complex interactions existing between turbulence and mixing in a blu?body stabilised ?ame con?guration is investigated by means of a wavelet crosscorrelation analysis on Large Eddy Simulations. The combined approach allows to better point out typical features of unsteady turbulent ?ows with mixing through the characterisation of the processes involved both in time and scales. The wavelet cross-correlation analysis of the time signals of velocity and mixture fraction ?uctuations can be an an e?ective tool to study the processes involved in turbulent mixing ?ows which are of great interest in combustion problems.

2

Generalities on wavelet cross-correlation

The continuous wavelet transform of a function f (t) is de?ned as the convolution between f and a dilated function ψ called wavelet mother: 1 Wf (a, τ ) = √ a
+∞

f (t)ψ ? (
?∞

t?τ )dt, a

(1)

where a is the dilation parameter, which plays the same role as the frequency in Fourier analysis, and τ indicates the translation parameter corresponding to the position of the wavelet in the physical space. In the present study we use 2 the complex Morlet wavelet (ψ(t) = eiω0 t e?t /2 ) as wavelet mother. Let Wf (a, τ ) and Wg (a, τ ) be the continuous wavelet transforms of f (t) and g(t). We de?ne the wavelet cross-scalogram as
? Wf g (a, τ ) = Wf (a, τ )Wg (a, τ ),

(2)

2

where the symbol ? indicates the complex conjugate. When the wavelet mother is complex, the wavelet cross-scalogram Wf g (a, τ ) is also complex and can be written in terms of its real and imaginary parts: Wf g (a, τ ) = CoWf g (a, τ ) ? iQuadWf g (a, τ ).
+∞ +∞ +∞

(3)
2

It can be shown that the following equation holds if f (t), g(t) ∈ L (?) f (t)g(t)dt = 1/cψ
?∞ 0 ?∞

CoWf g (a, τ )dτ da,

(4)

where 1/cψ is a constant depending on the choice of the wavelet mother.

3

Cross wavelet coherence functions

The highly redundant information from a multiscale wavelet analysis of time series must be reduced by means of suitable selective procedures and quantities, in order to extract the main features correlated to an essentially intermittent dynamics. In this study, we analysed and compared the properties of two complementary wavelet local correlation coe?cents which are able to well evidence peculiar and anomalous local events associated to the vortex dynamics. More precisely, given two signals f (t) and g(t), we refer to the so-called Wavelet Local Correlation Coe?cent (Buresti et. al ), de?ned as: W LCC(a, τ ) = CoWf g (a, τ ) . | Wf (a, τ ) || Wg (a, τ ) | (5)

This quantity is essentially a measure of the phase coherence of the signals. Here we introduce the Cross Wavelet Coherence Function (CWCF) de?ned as: CW CF (a, τ ) = 2 | Wf g (a, τ ) |2 , | Wf (a, τ ) |4 + | Wg (a, τ ) |4 (6)

which is essentially a measure of the intensity coherence of the signals. Using the polar coordinates we can write the wavelet transforms of Wf (a, τ ), Wg (a, τ ) and Wf g (a, τ ) as: Wf (a, τ ) = ρf e?θf Wg (a, τ ) = ρg e?θg (7) (8)

Wf g (a, τ ) = ρf ρg e?(θg ?θf ) , and the Cross Wavelet Coherence Function can be written also as: CW CF (a, τ ) = 2ρ2 ρ2 f g . ρ4 + ρ4 g f

(9)

It is easy to observe the two basic properties of the function (6): CW CF (a, τ ) = 0 =? ρf = 0 0 ≤ CW CF ≤ 1 or ρg = 0 (10) (11)

? a, τ.

3

4

Numerical simulation

We considered a laboratory-scale axisymmetric ?ame of methane-air in a non con?ned blu?-body con?guration. More precisely, the burner consists of a 5.4 mm diameter methane jet located in the center of a 50 mm diameter cylinder. Air is supplied through a 100 mm outer diameter coaxial jet around the 50 mm diameter blu?-body. The Reynolds number of the central jet is 7000 (methane velocity =21 m/s) whereas the Reynolds number of the coaxial jet is 80000 (air velocity =25 m/s). This is a challenging test case for all the turbulence models, as well documented in the ERCOFTAC report (Chatou, 1994) . Moreover, due to the highly intermittent, unsteady dynamics involved and the high turbulence level, especially for the reactive case, the Large Eddy Simulation (LES) appears as the most adequate numerical approach (Sello et. al ).

5

Results and discussion

In this analysis we are mainly interested to relations existing between evolution of turbulence and mixing, for the reactive case. Previous DNS simulations on coaxial jets at di?erent Reynolds numbers, show the ability of the wavelet crosscorrelation analysis to better investigate the relations between mixing process and the dynamics of vorticity (Salvetti et. al ). Thus, the signals analysed here are velocity ?uctuations (for Reynolds stress contributions) and mixture fraction ?uctuations (for mixing evolution) from LES. As an example, Figure 1 shows the wavelet co-spectrum maps for a signi?cant time interval in the pseudo-stationary regime of motion. The main contributions to the Reynolds stress are evidenced by high intensity correlations (red) and anti-correlations (blue) regions, which evolve intermittently. The dominant frequencies involved are located around 130 Hz. For the mechanisms responsable of the evolution of mixing, we note that the same regions of high Reynolds stress correspond to high correlation, or cooperation, between velocity and mixture fraction ?uctuations, suggesting that, at the selected location, the same events of stretching and tilting of the vorticity layer, drive both Reynolds stress and mixing evolutions. Note that the large high value region located at low frequencies in the right map is statistically not signi?cant if we assume a proper red noise background spectrum. To better investigate the role of the high correlation regions, we performed a cross section in the wavelet map at the frequency 160 Hz. Figure 2 (left) shows the time behaviour of the coherence functions WLCC, eq.(5), and CWCF, eq.(6). Here the phase and intensity coherence of signals are almost equivalent, but we can clearly point out an important anomalous event occurred at around t=0.19 s, corresponding to a loss of both intensity and phase coherence, followed by a change of the correlation sign. The link between this event and the dynamics of vorticity is evidenced by Figure 2 (right), which displays the wavelet map of the related vorticity signal. The higher frequency signi?cant regions (≈ 730 Hz) result strongly intermittent, with a bifurcation to lower and higher values than

4

average, followed by a drop of activity, in phase with the anomalous event.

Figure 1: Cross-Wavelet co-spectrum maps for axial and radial velocity ?uctuations (left) and for axial velocity and mixture fraction ?uctuations (right) at a given spatial point near the edge of the central jet.

Figure 2: Coherence functions for axial velocity and mixture fraction ?uctuations (left) and wavelet map of vorticity time series (right). These few examples support the usefulness of the cross-wavelet analysis approach to better investigate turbulent mixing processes in real systems.

References
 G. Buresti and G. Lombardi. Application of continuous wavelet transforms to the analysis of experimental turbulent velocity signals. Proc. of the 1st Int. Symp. on Turb. Shear Flow Phen., S. Barbara USA, Sept. 1999.  EDF Direction des Etudes et Researches. 1st A.S.C.F. Workshop Final Results, Chatou, France, October 1994.  S. Sello and G. Mariotti. Large eddy simulation of a blu? body stabilised ?ame. Proc. of the 4th ETMM Int. Symp., Ajaccio France, May 1999.  M.V. Salvetti, G. Lombardi and F. Beux. Application of a wavelet crosscorrelation technique to the analysis of mixing. AIAA Jour., 37:1007–1009, 1999.

This figure "fig1a_sel.gif" is available in "gif" format from: http://arXiv.org/ps/physics/0003029v1

This figure "fig1b_sel.gif" is available in "gif" format from: http://arXiv.org/ps/physics/0003029v1

This figure "fig2a_sel.gif" is available in "gif" format from: http://arXiv.org/ps/physics/0003029v1

This figure "fig2b_sel.gif" is available in "gif" format from: http://arXiv.org/ps/physics/0003029v1