[hnemati]

Dears
I am doing Direct Numerical simulation at high Pr numbers.
I have question regarding the advection of scalar. Can I use flux limiter; like Koren limiter; for scalar?
I have to say that If I don't use the limiter, the scalar gives strange number which are out of bounds?
Thanks

[FMDenaro]

Originally Posted by hnemati

Dears
I am doing Direct Numerical simulation at high Pr numbers.
I have question regarding the advection of scalar. Can I use flux limiter; like Koren limiter; for scalar?
I have to say that If I don't use the limiter, the scalar gives strange number which are out of bounds?
Thanks

In DNS you should have such a fine grid to resolve all physical wavenumbers of your fields. Therefore, you should be able to solve your problem without any limiter. What is more, using a limiter it should not produce relevant effect if the relevant small scales are well resolved.

However, if you are solving a pure advection of scalar, that is Df/Dt=0, using a sharp discontinuities then no grid exists that ensures the absence of oscillation (apart the case of first order scheme) and a limiter is necessary for numerical reasons.