Inviscid Supersonic Flow Simulation
I am currently undertaking a Phd in hypersonics, and I am required to perform some supersonic CFD simulations in OpenFOAM. I am quite new to this, so I have a few (and I am sure I will have plenty more) questions.
Initially I am simulating a wedge in a supersonic flow. I have based my initial model on the wedge15Ma5 tutorial case in rhopSonicFoam.
My first question is when playing around with the wedge tutorial, if I increase the Mach number the solution rapidly diverges (courant number increases exponentially until the solver crashes). I altered the mach number in the 0/U file (changed it from 5 to 10, and also tried other values in between this range). Any ideas why this is happening?
My second question is what are the differences between rhoSonicFoam and rhopSonicFoam? I know that rhop is density pressure solver, but what does this mean? Am I using the right solvers in this case? What are their limitations etc?
If anyone could point me in the right direction as to where I could find this information, that would be great. Thanks,
Did you reduce the time step or enable the automatic time stepping to respect the CFL condition? I run it with Co = 1. However for better results stay under 0.5.
You might want to consider also rhoCentralFoam, which is based on a Riemann-free type of scheme developed by Turganov and Tadmor (the reference to the paper is in the code). I tried it on the same tutorial, and with the same changes I suggested above it runs OK.
Thanks Alberto, I had tried adjusting the time step, but it still diverged. Now that I have put adaptive time step on, it doesn't diverge. Thanks again for that.
I'm having a bit of trouble with a mesh I am generating based on this tutorial. I am essentially trying to model a cone in a supersonic flow. My blockMeshDict file is:
Any ideas as to what I am doing wrong:confused:
Thanks again very much for your help!
|All times are GMT -4. The time now is 10:23.|