21.04.2014 Views

these simulation numerique et modelisation de l'ecoulement autour ...

these simulation numerique et modelisation de l'ecoulement autour ...

these simulation numerique et modelisation de l'ecoulement autour ...

SHOW MORE
SHOW LESS

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

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

The sensor used in the 2 nd or<strong>de</strong>r operator is also the maximum sensor ζ COL<br />

max .<br />

2 nd and 4 th or<strong>de</strong>r operators are then applied on each species.<br />

For energy we also build the modified coefficient based on ζ COL<br />

max<br />

A.4 The 4 mo<strong>de</strong>ls implemented in AVBP<br />

smu4 ′ = max(0, smu4 − ζ COL<br />

max smu2) (A.28)<br />

and we apply the 2 nd or<strong>de</strong>r AV with the same sensor.<br />

The subtilty comes when <strong>de</strong>aling with momentum. The 4 th or<strong>de</strong>r operator is not applied<br />

smu4 ′ = 0<br />

(A.29)<br />

and the 2 nd or<strong>de</strong>r operator uses the energy sensor ζ COL<br />

E<br />

instead of the maximum sensor ζ COL<br />

max .<br />

This mo<strong>de</strong>l is particularly <strong>de</strong>dicated to LES of reactive flows. The lack of 4 th or<strong>de</strong>r AV on momentum<br />

allows to keep many small scale structures, while damping the wiggles on energy and species.<br />

A.4.4<br />

“SLK” mo<strong>de</strong>l<br />

This is an improvement of the “Colin” mo<strong>de</strong>l. We have noticed that applying no 4 th or<strong>de</strong>r AV at all on<br />

momentum can y<strong>et</strong> lead to non-negligible wiggles in some cases (often <strong>de</strong>pending on the quality of the<br />

mesh) and this had to be avoi<strong>de</strong>d. This new “SLK” mo<strong>de</strong>l is very similar to the “Colin” mo<strong>de</strong>l, except<br />

that instead of s<strong>et</strong>ting the modified 4 th or<strong>de</strong>r coefficient to zero for momentum equation, it is s<strong>et</strong> to 10 %<br />

of the value used for the other equations. This gives :<br />

ζE SLK = ζΩ C j<br />

(ρE), ζY SLK = max<br />

k=1,neqs ζC Ω j<br />

(ρ k ) and ζmax<br />

SLK<br />

= max(ζE<br />

SLK , ζY SLK ) (A.30)<br />

smu4 ′ Y,E = max(0, smu4 − ζmax SLK smu2) and smu4 ′ U = 1 max(0, smu4 − ζSLK max smu2)<br />

10<br />

(A.31)<br />

This mo<strong>de</strong>l is very well suited for computations on poor quality meshes that exhibit velocity wiggles<br />

with the standard “Colin” mo<strong>de</strong>l.<br />

201

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

Saved successfully!

Ooh no, something went wrong!