CFD Online Discussion Forums

CFD Online Discussion Forums (https://www.cfd-online.com/Forums/)
-   OpenFOAM Running, Solving & CFD (https://www.cfd-online.com/Forums/openfoam-solving/)
-   -   LES model on surfaceMountedCube tutorial (https://www.cfd-online.com/Forums/openfoam-solving/243526-les-model-surfacemountedcube-tutorial.html)

Teresa Sun June 23, 2022 08:31

LES model on surfaceMountedCube tutorial
 
The tutorial copied from the link below can not be excuted.
https://develop.openfoam.com/Develop...aceMountedCube

The problem is that the solver stops before time 0, neither continuing nor throwing error message

I tried to debug with following tests:
a. decrease the mesh resolution by twice
b. increase viscosity, decrease Reynolds number by ten times
c. add finite volume option
d. try different solver pimpleFoam and pisoFoam
e. try different LES model like Smagorinsky, spalartAllmaras and QR(developed by our group)

The result doesn't change, the solver still stops before time 0, neither continuing nor throwing error message

Teresa Sun June 23, 2022 10:07

The output in log.pimpleFoam
 
1 Attachment(s)
Code:

\*---------------------------------------------------------------------------*/
Build  : _b45f8f6f58-20200629 OPENFOAM=2006
Arch  : "LSB;label=32;scalar=64"
Exec  : pisoFoam -parallel
Date  : Jun 22 2022
Time  : 17:04:39
Host  : pg-node006
PID    : 9472
I/O    : uncollated
Case  : /scratch/p302694/sufaceMountedCube/spalart/fullCase
nProcs : 8
Hosts  :
(
    (pg-node006 1)
    (pg-node024 1)
    (pg-node029 1)
    (pg-node045 1)
    (pg-node050 1)
    (pg-node054 1)
    (pg-node058 1)
    (pg-node059 1)
reate time

Create mesh for time = 0


PISO: Operating solver in PISO mode

Reading field p


Pstream initialized with:
    floatTransfer      : 0
    nProcsSimpleSum    : 0
    commsType          : nonBlocking
    polling iterations : 0
trapFpe: Floating point exception trapping enabled (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster (fileModificationSkew 5, maxFileModificationPolls 20)
allowSystemOperations : Allowing user-supplied system call operations

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

Create time

Create mesh for time = 0


PISO: Operating solver in PISO mode

Reading field p


Reading field U

Turbulent DFSEM patch inlet: interpolating field R from "/scratch/p302694/sufaceMountedCube/spalart/fullCase/constant/boundaryData/inlet/0/R"
Turbulent DFSEM patch inlet: interpolating field L from "/scratch/p302694/sufaceMountedCube/spalart/fullCase/constant/boundaryData/inlet/0/L"
Turbulent DFSEM patch inlet: interpolating field U from "/scratch/p302694/sufaceMountedCube/spalart/fullCase/constant/boundaryData/inlet/0/U"
Reading/calculating face flux field phi

Selecting incompressible transport model Newtonian
Selecting turbulence model type LES
Selecting LES turbulence model Smagorinsky
Selecting LES delta type IDDESDelta
Selecting LES hmax type maxDeltaxyzCubeRoot
Selecting patchDistMethod meshWave
LES
{
    turbulence      on;
    LESModel        Smagorinsky;
    printCoeffs    on;
    delta          IDDESDelta;
    IDDESDeltaCoeffs
    {
        hmax            maxDeltaxyzCubeRoot;
        maxDeltaxyzCubeRootCoeffs
        {
        }
    }
    Ce              1.048;
    Ck              0.094;
}

No MRF models present

No finite volume options present

Starting time loop

fieldAverage fieldAverage1:
    Restarting averaging for fields:
        U: starting averaging at time 0
        p: starting averaging at time 0


chen112p June 25, 2022 08:10

works fine for me. maybe post your full case here. I'm suspecting issue comes from your parallel setting.


All times are GMT -4. The time now is 01:26.