[akamal1099]

Hello my friends,
I am experiencing some issues with KOmegaSSTLM using simpleFoam.
I am trying to simulate flow over an airfoil with slotted flap.

I used the boundary conditions which is found in the T3A tutorial for KOmegaSSTLM but the solution diverges.

I tried some troubleshooting and I will write down below summary of it:

1- Everything is fine when there is no layers around the airfoil and flap.
2- When I added layers the solution starts diverging for unknown reason.
3- Same thing happens with KOmegaSST

My layer was made to target y+ of 1 and my mesh is fine per meshCheck.
I am suspecting and error with boundary conditions due to me using wallfunctions for nut, omega and k.


Thanks my friends
Following is the omega file as an example.
I edited it to be zeroGradient because I have y+ equal to 1 and I don't need wallfunctions but still don't work. (Same for k and nut)

also my mesh is attached

/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  v2012                                 |
|   \\  /    A nd           | Website:  www.openfoam.com                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       volScalarField;
    location    "0";
    object      omega;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

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

internalField   uniform 1.7520318313820673;

boundaryField
{
    inlet
    {
        type            fixedValue;
        value           $internalField;
    }

    outlet
    {
        type            zeroGradient;
    }

	airfoilZeroDegSurf
    {
        type            zeroGradient;
    }
	
	flapFiftyDegSurf
    {
        type            zeroGradient;
    }

    upperEmpty
    {
        type            empty;
    }
	
	lowerEmpty
    {
        type            empty;
    }
	
	front
    {
        type            symmetry;
    }
	back
    {
        type            symmetry;
    }
}

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

[akamal1099]

Guys any idea?

[farzadmech]

Hello
did you know how to solve it? my simulation is also diverges.


Thanks,
Farzad

Originally Posted by akamal1099

Guys any idea?