[yzf1215]

I'm running simpleFoam in OpenFOAM using k-omega-SST model on calculating flow around airfoil. But the solution blow up after several hundreds of iterations. Would any help me with that?

Here are the flow conditions:
U=139.4m/s
p=72319.7 Pa
Re=10e6
airfoil cord=1 m
computational domain length 40 m

Timestep=10e-6; should it be smaller or larger?
y+=1;



Here are the Boundary Conditions:

P boundary condition:



internalField uniform 72319.7;


AIRFOIL1
{
type zeroGradient;
}
INLET
{
type fixedValue;
value uniform 72319.7;
}
OUTLET
{
type zeroGradient;
}
UPPERWALL
{
type zeroGradient;
}
LOWERWALL
{
type zeroGradient;
}


U boundary condition:


internalField uniform (139.4 0 0);
AIRFOIL1
{
type fixedValue;
value uniform (0 0 0);
}
INLET
{
type fixedValue;
value (139.4 0 0);
}
OUTLET
{
type zeroGradient;
}
UPPERWALL
{
type zeroGradient;
}
LOWERWALL
{
type zeroGradient;
}

k boundary condition
internalField uniform 2.8982;

{
AIRFOIL1
{
type kqRWallFunction;
value uniform 0;
}
frontPlate
{
type symmetryPlane;
}
top
{
type symmetryPlane;
}
INLET
{
type fixedValue;
value uniform 2.8982;
}
OUTLET
{
type zeroGradient;
}
UPPERWALL
{
type kqRWallFunction;
value uniform 0;
}
LOWERWALL
{
type kqRWallFunction;
value uniform 0;
}

omega boundary condition
internalField uniform 697;

boundaryField
{
AIRFOIL1
{
type omegaWallFunction;
value uniform 0;

}
INLET
{
type fixedValue;
value uniform 697;
}
OUTLET
{
type zeroGradient;
}
UPPERWALL
{
type zeroGradient;
}
LOWERWALL
{
type zeroGradient;
}