Compressible calculation with Code_Saturne
I am a new user of Code_Saturne and I am currently working on a project regarding the computation of the flow around an airfoil (NACA23012 here). The mesh was created using Salome and I successfully launched some incompressible cases with Code_Saturne using different meshes and turbulence models (without the GUI interface).
The mesh width is only 1 cell. The boundary conditions are the following : "wall" on the airfoil, "inlet" on all the farfield (which is infow or outflow depending on the direction) and "symmetry" at z=0 and z=1.
I am now trying to activate the compressible module in Code_Saturne and I am currently having trouble to set up the environment. I modified the following user files for my purpose : uscfcl.f90, uscfxi.f90, usini1.f90, usppmo.f90 and usproj.f90.
In particular, I added the "coumax" keyword in usini1.f90 and set it to 0.4 as adviced in the user manual (with idtvar=1 for variable time step).
Most of the parameters remain the same between incompressible and compressible cases (initialisation ...) and I tried to follow the user manual for specific parameters.
I compute the lift coefficient in usproj.f90 at each time step so I can follow its convergence and here is the problem : The calculation doesn't converge toward a fixed value but oscillate. I tried to change as many options as I can but the problem remains. Obviously, pressure and velocity fields as well as pressure coefficients around the airfoil are wrong.
Could somebody give me a lead to investigate the problem ? Or maybe a tutorial or some fortran files suitable for compressible calculations ?
I hope I am specific enough. Do not hesitate to answer for any further information.
The fortran user files and the mesh I used :
Thanks in advance for any help.
|All times are GMT -4. The time now is 00:37.|