[Zaffar]

Rob3rt,

Unfortunately, the problem has not been resolved yet.
And any help would be appreciated.

[Tian Jin]

Originally Posted by tidal_Tom

I do not know what is causing your errors. I am a relatively novice user so it is beyond me. below are my files I used for flm and fmm. I don't think fmm can be zero, as the calculation for turbulent viscosity divides by fmm, I set this to a value of 1. flm can be zero as it is the nominator.

Have you added flm and fmm in the fvSolution and fvSchemes? lines I added are below.

fvSolution:

flm
{
solver PBiCG;
preconditioner DILU;
tolerance 1e-05;
relTol 0;
}

fmm
{
solver PBiCG;
preconditioner DILU;
tolerance 1e-05;
relTol 0;
}

fvSchemes:

divSchemes
{
default none;
div(phi,U) Gauss linear;
div(phi,flm) Gauss limitedLinear 1;
div(phi,fmm) Gauss limitedLinear 1;
div(phi,B) Gauss limitedLinear 1;
div(B) Gauss linear;
div((nuEff*dev(T(grad(U))))) Gauss linear;
}

and I added:

dynLagrangianCoeffs
{
filter simple;
ce 1.048;
theta 1.5;
}

in the LESProperties dictionary.

flm:

/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 2.0.0 |
| \\ / A nd | Web: www.OpenFOAM.com |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class volScalarField;
object flm;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions [0 4 -4 0 0 0 0];

internalField uniform 0;

boundaryField
{
inlet
{
type fixedValue;
value uniform 0;
}

outlet
{
type zeroGradient;
}

outside
{
type fixedValue;
value uniform 0;
}

object
{
type fixedValue;
value uniform 0;
}

wedge1
{
type wedge;
}

wedge2
{
type wedge;
}
}

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

fmm:

/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 2.0.0 |
| \\ / A nd | Web: www.OpenFOAM.com |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class volScalarField;
object fmm;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions [0 4 -4 0 0 0 0];

internalField uniform 1;

boundaryField
{
inlet
{
type fixedValue;
value uniform 1;
}

outlet
{
type zeroGradient;
}

outside
{
type zeroGradient;
}

object
{
type zeroGradient;
}

wedge1
{
type wedge;
}

wedge2
{
type wedge;
}
}

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

Thank you very much. It can works. But when I run it in my case, I found that turbulent viscosity are 0 in the internalField. Do you know what's wrong with that?

[Tian Jin]

Originally Posted by anishtain4

Hi,

Since there is no clear example of dynamic Lagrangian LES model in the OpenFOAM tutorials, I attached a case I had to make work.
This is for OpenFOAM 3.0.1,
Let me know if you have any questions about it.

The codes to run it are as follow:

cp -r 0.org 0;
blockMesh;
perturbUChannel;
pimpleFoam;

One also needs

It's very kind of you. When I run your case, it gives me the error, do you know why it happened?
~/Desktop/test_solver/dynLag$ perturbUChannel
perturbUChannel: command not found