Hi Daniel,
Thank you so much for your reply.
Do you need to preserve a sharp interface or are you more interested
in a higher order scheme? The Van Leer term can be modified to not
have the limiting properties and just be a plain old QUICK scheme if
you prefer.
I just want to make sure the result is as accurate as possile, not due to
the method I choose. I used the periodic mesh, so "residual" or small
"turbulence left behind" will become a big error after several circles. It
may cause me have a wrong explanation for the influence of convection term.
I'm not sure what is a plain old QUICK scheme. Do you have more information
about it?
Another question is: can I describe the convection term as
convTerm=(convCoeff*var.getArithmeticFaceValue()).getDivergence()
convCoeff is a vector. What's the difference if I define the convection
term
by myself without using the (default) scheme.
You can certainly do that, but your scheme will be explicit and the
getDivergence operator doesn't know about boundary conditions. You can
of course include boundary conditions by modifying the term
appropriately.
I also had met a difficulty to define the boundary conditions for my scheme
(convTerm=(convCoeff*var.getArithmeticFaceValue()).getDivergence()).
Would please show me how to include the boundary conditions in my scheme, so
I may get a similar result as Van Leer or a higher order scheme?
Thanks again for your help.
Kundar
