First, I found a very nice program, called sumo, from a swedish firm http://www.larosterna.com/.
It is a simple CAD environment focused in designing airframes, and also provides options to generate surface and volume meshes for CFD simulation.
I decided to choose one of the default provided models, a small delta interceptor, and ended up with a surface mesh like this:
The volume mesh tool uses tetgen mesher http://tetgen.berlios.de/, and the mesh can be saved as .ele,.face and .node files, which can then be imported in OF via tetgenToFoam utility.
The result is a nice mesh of the airframe in a spherical boundary:
The OF solver of choice for the purpose, compressible, steady state and subsonic is rhoSimpleFoam.
The parameters for the simulation were uniform distributions with velocity 200m/s in the x direction, pressure 1e5 Pa, and temperature 300° K.
I have to say it is not so easy to obtain a convergent simulation in OF with tetrehedral unstructured meshes, it has been a difficult but instructive experience to get it.
My convergence criterion was 1e-4, and the process made it in 102 iterations, below the gnuplot out of foamLog:
Finally, a paraFoam screenshot of pressure distribution:
The lowest bound for pressure seems too low, it refers to a few cells spotting value greatly differing from their neighbours.
My next steps will be investigating these anomalies, and maybe see if rhoSimpleFoam can be used for supersonic flows.
