some error using rhoCentralFoam
I have used rhoCentralFoam to calculate a supersonic nozzle. When I use the upwind scheme, it's all right. But when I change the interpolation scheme to the second order scheme, such as quadraticUpwind, there is a zone next to the wall boundary in which the velocity is lager than the main steam.See the attachment, the two figures show the mach number distribute in the nozzle.
Could someone tell me that why does this happen?
The following is the fvSchem file:
default Gauss linear;
div(tauMC) Gauss linear;
div(phi,k) Gauss Minmod phi 1;
default Gauss linear corrected;
interpolate(rho) quadraticUpwindFit phi 1;
reconstruct(rho) quadraticUpwindFit 1;
reconstruct(U) quadraticUpwindFit 1;
reconstruct(T) quadraticUpwindFit 1;
reconstruct(K) quadraticUpwindFit 1;
reconstruct(rho) linearUpwind phi;
reconstruct(U) linearUpwindV phi;
reconstruct(T) linearUpwind phi;
reconstruct(K) linearUpwind phi;
Is you mesh fine enough near wall?
Does it occur in linearUpwind too?
Or Have you found out the reason?
The case file is too large to upload. I'm not sure whether the mesh is fine enough. It occur in any other scheme except upwind. But the un-physical solution is not in the boundary area, and it develops from the throat. I'll go on to find out the reason. Best regards.
It seems as it linearUpwind scheme works well.
how did you make the mesh?
then the problem is resolved?:)
Sorry, there comes the latest news, the problem haven't been solved. The finest mesh is the one associated with the characteristc line.
If I using an apreciate discrete schem to solve the subsonic part of the flow well, then it will be fine, Because the high value region develops from this part.
|All times are GMT -4. The time now is 11:37.|