[ibelunatic]

I am trying to compute lift and drag from 3 airfoils, set up in a triangle in 2D using simpleFoam (OpenFOAM 4.1) and the kkl-omega turbulence model. Upon running the simulation, time step continuity errors occur and my omega, kl and kt blow up. I do not know if I am using the most appropriate fvSolution solvers. I would appreciate if someone could have a look. I have also included my boundary conditions. My polymesh passed the checkMesh fyi, I can submit it if required.

fvSolution:

FoamFile
{
version 2.0;
format ascii;
class dictionary;
location "system";
object fvSolution;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

solvers
{
p
{
solver GAMG;
tolerance 1e-08;
relTol 1.e-3;
smoother GaussSeidel;
nPreSweeps 2;
nPostSweeps 2;
cacheAgglomeration true;
nCellsInCoarsestLevel 10;
agglomerator faceAreaPair;
mergeLevels 1;
}

"(U|kt|kl|omega)"
{
solver smoothSolver;
smoother GaussSeidel;
nSweeps 2;
tolerance 1e-10;
relTol 1.e-3;
}

}

SIMPLE
{
nNonOrthogonalCorrectors 2;
pRefCell 0;
pRefValue 0;
}

relaxationFactors
{
p 0.1;
U 0.1;
kl 0.1;
kt 0.1;
omega 0.1;
}


// ************************************************** *********************** //


Boundary conditions:

U:

FoamFile
{
version 2.0;
format ascii;
class volVectorField;
object U;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions [0 1 -1 0 0 0 0];

internalField uniform (10.0 0 0);

boundaryField
{
exit
{
type freestream;
freestreamValue uniform (10.0 0 0);
}
inlet
{
type freestream;
freestreamValue uniform (10.0 0 0);
}
top
{
type slip;
}
bottom
{
type slip;
}
airfoila
{
type fixedValue;
value uniform (0 0 0);
}
airfoilb
{
type fixedValue;
value uniform (0 0 0);
}
airfoilc
{
type fixedValue;
value uniform (0 0 0);
}
front
{
type empty;
}
back
{
type empty;
}


}

// ************************************************** *********************** //

p:

FoamFile
{
version 2.0;
format ascii;
class volScalarField;
object p;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions [0 2 -2 0 0 0 0];

internalField uniform 0;

boundaryField
{
inlet
{
type zeroGradient;
}
exit
{
type fixedValue;
value uniform 0;
}
top
{
type slip;
}
bottom
{
type slip;
}
airfoila
{
type zeroGradient;
}
airfoilb
{
type zeroGradient;
}
airfoilc
{
type zeroGradient;
}
front
{
type empty;
}
back
{
type empty;
}
}

// ************************************************** *********************** //

omega:

FoamFile
{
version 2.0;
format ascii;
class volScalarField;
location "0";
object omega;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions [0 0 -1 0 0 0 0];

internalField uniform 1.74964;

boundaryField
{
inlet
{
type freestream;
freestreamValue $internalField;
}
exit
{
type freestream;
freestreamValue $internalField;
}
top
{
type slip;
}
bottom
{
type slip;
}
airfoila
{
type fixedValue;
value uniform 1e-11;
}
airfoilb
{
type fixedValue;
value uniform 1e-11;
}
airfoilc
{
type fixedValue;
value uniform 1e-11;
}
front
{
type empty;
}
back
{
type empty;
}
}


// ************************************************** *********************** //

kl:

FoamFile
{
version 2.0;
format ascii;
class volScalarField;
location "0";
object kl;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions [0 2 -2 0 0 0 0];

internalField uniform 0;

boundaryField
{
inlet
{
type freestream;
freestreamValue $internalField;
}
exit
{
type freestream;
freestreamValue $internalField;
}
top
{
type slip;
}
bottom
{
type slip;
}
airfoila
{
type fixedValue;
value uniform 0;
}
airfoilb
{
type fixedValue;
value uniform 0;
}
airfoilc
{
type fixedValue;
value uniform 0;
}
back
{
type empty;
}
front
{
type empty;
}
}


// ************************************************** *********************** //

kt:

FoamFile
{
version 2.0;
format ascii;
class volScalarField;
location "0";
object kt;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions [0 2 -2 0 0 0 0];

internalField uniform 0.015;

boundaryField
{
inlet
{
type freestream;
freestreamValue $internalField;
}
exit
{
type freestream;
freestreamValue $internalField;
}
top
{
type slip;
}
bottom
{
type slip;
}
airfoila
{
type fixedValue;
value uniform 1e-11;
}
airfoilb
{
type fixedValue;
value uniform 1e-11;
}
airfoilc
{
type fixedValue;
value uniform 1e-11;
}
front
{
type empty;
}
back
{
type empty;
}
}


// ************************************************** *********************** //

nut:

FoamFile
{
version 2.0;
format ascii;
class volScalarField;
location "0";
object nut;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions [0 2 -1 0 0 0 0];

internalField uniform 0;

boundaryField
{
inlet
{
type calculated;
value uniform 0;
}
exit
{
type calculated;
value uniform 0;
}
top
{
type calculated;
value uniform 0;
}
bottom
{
type calculated;
value uniform 0;
}
airfoila
{
type nutLowReWallFunction;
value uniform 0;
}
airfoilb
{
type nutLowReWallFunction;
value uniform 0;
}
airfoilc
{
type nutLowReWallFunction;
value uniform 0;
}
front
{
type empty;
}
back
{
type empty;
}
}


// ************************************************** *********************** //