|May 30, 2011, 10:26||
Compressible calculation with Code_Saturne
Join Date: Apr 2011
Posts: 1Rep Power: 0
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.
|airfoil, code, compressible, saturne|
|Thread||Thread Starter||Forum||Replies||Last Post|
|compressible flow calculation error using rhoSimpleFoam solver||student4326||OpenFOAM||6||February 10, 2012 10:36|
|MRF and Heat transfer calculation||Susan YU||FLUENT||0||June 2, 2010 08:46|
|Compressible -> incompressible.||Jinwon||Main CFD Forum||6||November 23, 2007 22:07|
|Warning 097-||AB||CD-adapco||6||November 15, 2004 05:41|
|Compressible vs. Incompressible formulations||Fernando Velasco Hurtado||Main CFD Forum||3||January 7, 2000 17:51|