[losiola]

Hello guys ,


i want to make a 2D airfoil simulation using the Kkl omega turbulence model using the pimpleFoam solver i have been looking through the tutorials and i ve found the elipskklomega tutorial.


Now i have 2 questions :


1/-once i run the tutorial the solutions diverges at 0.5X seconds , So why the tutorial is diverging if it is meant to be a model for us as users and learners of openFoam?


2/-Once i went through the files in the 0/ directory i ve found that they are not using any wall functions for kt kl or Omega !! does this mean that they are making a full wall resolution for the wall (elliptic shape ) ??!!


and these are the contents of the files


kl

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

internalField   uniform 0;

boundaryField
{
    inlet
    {
        type            fixedValue;
        value           uniform 0;
    }
    outlet
    {
        type            inletOutlet;
        inletValue      uniform 0;
        value           uniform 0;
    }
    up
    {
        type            symmetry;
    }
    hole
    {
        type            fixedValue;
        value           uniform 0;
    }
    frontAndBack
    {
        type            empty;
    }
 }

kt

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

internalField   uniform 0;

boundaryField
{
    inlet
    {
        type            fixedValue;
        value           uniform 0;
    }
    outlet
    {
        type            inletOutlet;
        inletValue      uniform 0;
        value           uniform 0;
    }
    up
    {
        type            symmetry;
    }
    hole
    {
        type            fixedValue;
        value           uniform 0;
    }
    frontAndBack
    {
        type            empty;
    }
}

nut

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

internalField   uniform 0;

boundaryField
{
    inlet
    {
        type           calculated;
        value          uniform 0;
    }
    outlet
    {
        type           calculated;
        value          uniform 0;
    }
    up
    {
        type           symmetry;
    }
    hole
    {
        type           nutLowReWallFunction;
        value          uniform 0;
    }
    frontAndBack
    {
        type            empty;
    }
}

Omega

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

internalField   uniform 1e-5;

boundaryField
{
    inlet
    {
        type            fixedValue;
        value           $internalField;
    }
    outlet
    {
        type            inletOutlet;
        inletValue      $internalField;
        value           $internalField;
    }
    up
    {
        type            symmetry;
    }
    hole
    {
        type            zeroGradient;
    }
    frontAndBack
    {
        type            empty;
    }
}

p

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

internalField   uniform 0;

boundaryField
{
    inlet
    {
       type             zeroGradient;
    }
    outlet
    {
        type            fixedValue;
        value           uniform 0;
    }
    up
    {
        type            symmetry;
    }
    hole
    {
        type            zeroGradient;
    }
    frontAndBack
    {
        type            empty;
    }
}

U

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

internalField   uniform (3.1 0 0);

boundaryField
{
    outlet
    {
       type             inletOutlet;
       inletValue       uniform (0 0 0);
       inlet            uniform (0 0 0);
    }
    inlet
    {
        type            fixedValue;
        value           uniform (3.1 0 0);
    }
    up
    {
        type            symmetry;
    }
    hole
    {
        type            noSlip;
    }
    frontAndBack
    {
        type            empty;
    }
}

Thank You