13.07.2015 Views

Analyse expérimentale et modélisation du transfert de matière et du ...

Analyse expérimentale et modélisation du transfert de matière et du ...

Analyse expérimentale et modélisation du transfert de matière et du ...

SHOW MORE
SHOW LESS
  • No tags were found...

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

Annexe Chapitre 210864u'v' v' 220-2-4-6u'v' (DNS Rogers <strong>et</strong> al.[1])-8 v' 2 (DNS Rogers <strong>et</strong> al.[1])-10u'v' (RS Laun<strong>de</strong>r <strong>et</strong> al.[10])v' 2 (RS Laun<strong>de</strong>r <strong>et</strong> al.[10])-120 5 10 15StFigure 1: Development of the Reynolds stress tensor components of Reynolds stress tensor in a sheared turbulence(DNS by Rogers <strong>et</strong> al. [1])Once the Reynolds stress transport mo<strong>de</strong>lassessed, the analysis is then focussed onthe second or<strong>de</strong>r mo<strong>de</strong>lling of th<strong>et</strong>ransport equation of the turbulentpassive scalar. For this purpose theclosure of the pressure-scalar gradientcorrelations and of the <strong>de</strong>struction rat<strong>et</strong>ensor are evaluated, and the numericalresults are compared with the (DNS) dataof Rogers <strong>et</strong> al. [1].For this purpose different closures of thepressure-scalar gradient term were tested.In particular the three closures proposedby Wikström <strong>et</strong> al. [6] (mo<strong>de</strong>ls a, b, andc) were tested and an additionalformulation of the mo<strong>de</strong>l is proposed.This mo<strong>de</strong>l writes:DDt− C⎛' '' '⎜∂Ui ' 'uiθ= − ujθ+ uiuj⎝ ∂xjεu θ + Ck' 'θ1 i θ 4' ' ∂Θuiuj∂xj∂Θ ⎞⎟∂xj ⎠(15)In equation (15), the constants of theoriginal Wikström <strong>et</strong> al. [6] mo<strong>de</strong>l arefixedtor + 1C1= 1.3 ; C 35θ θ 4= 0.rC C = C 0 .θ 2=θ 3 θ 5=andIn these simulations the transportequation of the <strong>de</strong>struction rate ε θof thescalar variance q is computed using thefollowing constants ( C = 2. 4 , C = 1. 89 ,γ 1 γ 2*=*C 0.1 and C = 0. 95 ). The results ofγ 1 γ 2these simulations are reported in figure(2) and compared to the (DNS) data ofRogers <strong>et</strong> al. [1]. Figure (2) shows a goodagreement b<strong>et</strong>ween the DNS data ofRogers <strong>et</strong> al. [1] and the numerical results'of the turbulent passive scalar fluxes u'θ'and v'θ obtained with the second or<strong>de</strong>rmo<strong>de</strong>l as it is adjusted above.169

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!