Running les for aerofoil

Project Support
1

I am trying to run les simulation for flow over Naca0012 using WALE subgrid model.I first created the 3-D domain and run the case for steady state case using simpleFoam and then i used the converged solution as initial solution for the les by renaming the folder (6000 to 1e-5) then i decomposed the domain. I am using pimpleFoam for les simulation but the problem is that the simulation is not converging and the pressure equation solver reach the maximum no fo iteration (1000)without reaching a converged solution.

Is there any examples or have anyone worked on les for aerofoil ?

This the RANS results which is set as initial value for les
the velocity fields

image
image950Ă—563 38.2 KB

pressure fields
image
image768Ă—464 43.3 KB

after running les simulation for few time step results seems to be not converged
the velocity field
image
image756Ă—456 23.1 KB

pressure fields
image
image758Ă—455 76 KB

The boundary conditions
u
/--------------------------------- C++ -----------------------------------
| ========= | |
| \ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \ / O peration | Version: 5.x |
| \ / A nd | Web: www.OpenFOAM.org |
| \/ M anipulation | |
*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class volVectorField;
location “0”;
object U;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

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

internalField uniform (71.3 0 0);

boundaryField
{
aerofoil
{
type fixedValue;
value uniform (0 0 0);
}
top
{
type symmetryPlane;
}
bottom
{
type symmetryPlane;
}
inlet
{
type fixedValue;
value uniform (71.3 0 0);
}
outlet
{
{
type freestream;
freestreamValue uniform (71.3 0 0);
value uniform (71.3 0 0);
}
}
front
{
type cyclic;
}
back
{
type cyclic;
}
}

// ************************************************************************* /
p
/--------------------------------- C++ -----------------------------------
| ========= | |
| \ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \ / O peration | Version: 5.x |
| \ / A nd | Web: www.OpenFOAM.org |
| \/ M anipulation | |
*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class volScalarField;
location “0”;
object p;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

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

internalField uniform 0;

boundaryField
{
aerofoil
{
type zeroGradient;
}
top
{
type symmetryPlane;
}
bottom
{
type symmetryPlane;
}
inlet
{
type zeroGradient;
}
outlet
{
type fixedValue;
value uniform 0;
}
front
{
type cyclic;
}
back
{
type cyclic;
}
}

// ************************************************************************* /
nut
/--------------------------------- C++ -----------------------------------
| ========= | |
| \ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \ / O peration | Version: 5.x |
| \ / A nd | Web: www.OpenFOAM.org |
| \/ M anipulation | |
*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class volScalarField;
location “0”;
object nut;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

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

internalField uniform 0.3;

boundaryField
{
aerofoil
{
type zeroGradient;
}
top
{
type symmetryPlane;
}
bottom
{
type symmetryPlane;
}
inlet
{
type calculated;
value uniform 0.3;
}
outlet
{
type calculated;
value uniform 0.3;
}
front
{
type cyclic;
}
back
{
type cyclic;
}
}

// ************************************************************************* //
turbulence properties
/--------------------------------- C++ -----------------------------------
| ========= | |
| \ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \ / O peration | Version: 2.2.0 |
| \ / A nd | Web: www.OpenFOAM.org |
| \/ M anipulation | |
*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class dictionary;
location “constant”;
object turbulenceProperties;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

simulationType LES;

LES
{
LESModel WALE;
turbulence on;
printCoeffs on;

delta vanDriest;

vanDriestCoeffs
{
delta cubeRootVol;
cubeRootVolCoeffs
{
deltaCoeff 1;
}

Aplus           26;
Cdelta          0.158;

}
}
// ************************************************************************* //
fvschemes
/--------------------------------- C++ -----------------------------------
| ========= | |
| \ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \ / O peration | Version: 2.2.0 |
| \ / A nd | Web: www.OpenFOAM.org |
| \/ M anipulation | |
*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class dictionary;
location “system”;
object fvSchemes;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

ddtSchemes
{
default backward;
}

gradSchemes
{
default Gauss linear;
grad(p) Gauss linear;
grad(U) Gauss linear;
}

divSchemes
{
default none;
div(phi,U) Gauss linear;
div(phi,k) Gauss limitedLinear 1;
div(phi,B) Gauss limitedLinear 1;
div(B) Gauss linear;
div(phi,nuTilda) Gauss limitedLinear 1;
div((nuEff*dev2(T(grad(U))))) Gauss linear;
}

laplacianSchemes
{
default none;
laplacian(nuEff,U) Gauss linear corrected;
laplacian((1|A(U)),p) Gauss linear corrected;
laplacian(DkEff,k) Gauss linear corrected;
laplacian(DBEff,B) Gauss linear corrected;
laplacian(DnuTildaEff,nuTilda) Gauss linear corrected;
}

interpolationSchemes
{
default linear;
interpolate(U) linear;
}

snGradSchemes
{
default corrected;
}

fluxRequired
{
default no;
p ;
}

// ************************************************************************* //
fvsolution
/--------------------------------- C++ -----------------------------------
| ========= | |
| \ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \ / O peration | Version: v1906 |
| \ / A nd | Web: www.OpenFOAM.com |
| \/ M anipulation | |
*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class dictionary;
location “system”;
object fvSolution;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

solvers
{
p
{
solver GAMG;
tolerance 0;
relTol 0.1;
smoother GaussSeidel;
}

pFinal
{
    $p;
    smoother        DICGaussSeidel;
    tolerance       1e-06;
    relTol          0;
}

"(U|k|nuTilda)"
{
    solver          smoothSolver;
    smoother        GaussSeidel;
    tolerance       1e-05;
    relTol          0.1;
    minIter         1;
}

"(U|k|nuTilda)Final"
{
    $U;
    tolerance       1e-05;
    relTol          0;
}

}

PIMPLE
{
nOuterCorrectors 10;
nCorrectors 2;
nNonOrthogonalCorrectors 1;
pRefCell 0;
pRefValue 0;

residualControl //(or residualControl, depends on the version of OF)
{
    p
   { 
      tolerance 1e-4; 
      relTol 0;
   } // I was said 1e-4 is enough
    U 
  { 
     tolerance 1e-4; 
     relTol       0;
  }
}

relaxationFactors
{
fields
{
p 0.7;
pFinal 1;
}
equations
{
“(U|k|nuTilda)” 0.7;
“(U|k|nuTilda)Final” 1;
}
}

}

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

2

Hi @zeinelserfy!

Could you share the project link with us? Tagging the @cfd_squad here who might guide you through the process and give you some useful tips.

Best,

Jousef

3

Project link please?

4

how to upload project link?
I can upload the files like that of openFoam ?

5

I quickly assumed that you already had a SimScale project that is not converging, if not, then I am not sure how much help you will get here unless one of our openfoam experts jumps in …

6

simscale is based on openFoam, right?

7

Yes but I for one, can’t help you unless you set that project up in SimScale… (I am not familiar with a lot of what you describe and I do not know if that setup is even possible in SimScale)

10

Hi @zeinelserfy!

We do not use external resources to help our users, would you be willing to share or re-post the thread and integrate it into this one? After that either edit your post with the link or remove it, up to you :slight_smile: From there we can then start investigating the issue(s).

Cheers and thanks!

Jousef

11

i have updated my post if anyone can help

13

Post closed! In case you want to run your simulation on SimScale, feel free to re-open a new post.

Cheers and thanks!

Jousef

1 Like