CFD Online Discussion Forums

CFD Online Discussion Forums (https://www.cfd-online.com/Forums/)
-   OpenFOAM (https://www.cfd-online.com/Forums/openfoam/)
-   -   NACA0012 compressible static simulation (https://www.cfd-online.com/Forums/openfoam/87231-naca0012-compressible-static-simulation.html)

paolopaolo88 April 14, 2011 10:21

NACA0012 compressible static simulation
 
1 Attachment(s)
Hi everybody,
I am a new user of openFoam, and I would like to use it to do some dynamic stall simulations for a naca 0012 for M=0.7-0.85.

I think that a good solver for this purpose is rhoPisoFoam with a kOmegaSST turbolence model with a later implementation of the rotation of the grid.
So I was trying to launch a first static simulation at M=0.5 with rhoPisoFoam and the kOmegaSST model at zero angle of attack to validate the solver with some experimental data.

Unfortunatly If I set case with the kOmegaSST model using these conditions
http://turbmodels.larc.nasa.gov/sst.html
I get this error, that lets me think that there is a problem with my thermoPhysicalProperties or the schemes I am using to evaluate them.
I post the case without the grid (that is not my property).

the length of a cell in the flow direction is approximately 0.0045 and in y direction near the wall is 0.0000055. As you can notice the grid is very fine.
Other info about the grid are:
Quote:

(
body
{
type wall;
physicalType wallAdiabatic;
nFaces 970;
startFace 385090;
}
defaultFaces
{
type empty;
nFaces 386060;
startFace 386060;
}
inlet-outlet
{
type patch;
nFaces 970;
startFace 772120;
}
)
The error I get is this one:

Quote:

--> FOAM FATAL ERROR:
Maximum number of iterations exceeded

From function specieThermo<thermo>::T(scalar f, scalar T0, scalar (specieThermo<thermo>::*F)(const scalar) const, scalar (specieThermo<thermo>::*dFdT)(const scalar) const) const
in file /home/opencfd/OpenFOAM/OpenFOAM-1.7.x/src/thermophysicalModels/specie/lnInclude/specieThermoI.H at line 67.

FOAM aborting

#0 Foam::error::printStack(Foam::Ostream&) in "/opt/openfoam171/lib/linux64GccDPOpt/libOpenFOAM.so"
#1 Foam::error::abort() in "/opt/openfoam171/lib/linux64GccDPOpt/libOpenFOAM.so"
#2 Foam::hPsiThermo<Foam::pureMixture<Foam::constTran sport<Foam::specieThermo<Foam::hConstThermo<Foam:: perfectGas> > > > >::calculate() in "/opt/openfoam171/lib/linux64GccDPOpt/libbasicThermophysicalModels.so"
#3 Foam::hPsiThermo<Foam::pureMixture<Foam::constTran sport<Foam::specieThermo<Foam::hConstThermo<Foam:: perfectGas> > > > >::correct() in "/opt/openfoam171/lib/linux64GccDPOpt/libbasicThermophysicalModels.so"
#4
in "/opt/openfoam171/applications/bin/linux64GccDPOpt/rhoPisoFoam"
#5 __libc_start_main in "/lib/libc.so.6"
#6
in "/opt/openfoam171/applications/bin/linux64GccDPOpt/rhoPisoFoam"
Aborted
I really hope someone can help me,
thank you very much for any suggestions.

Paolo

paolopaolo88 April 14, 2011 10:29

Another thing that can help is this one:

when I try to launch the case using a much less fine unstructured grid, adjusting the omega wall parameter and the deltaT to a much lower value, the case works (the results converge at least).

For the finer mesh case I am using a DT=1.5e-05 but the grid is so fine that I get a max Co Number of some hundreds.

Thank you again.
Paolo


All times are GMT -4. The time now is 07:39.