SimCenter Forum

Research in Natural Hazards Engineering [Archived] => Wind Engineering (WE-UQ) => Topic started by: Jack on September 14, 2021, 10:37:00 PM

Title: PimpleFoam freezes in first timestep when using turbulentDFSEMInlet
Post by: Jack on September 14, 2021, 10:37:00 PM

I was using PimpleFoam with turbulentDFSEMInlet on a supercomputer and the simulation freezes at the first timestep. Has anyone experienced this before? Is there any way I can use to debug this?

'''
Pstream initialized with:
    floatTransfer      : 0
    nProcsSimpleSum    : 0
    commsType          : nonBlocking
    polling iterations : 0
sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster (fileModificationSkew 10)
allowSystemOperations : Allowing user-supplied system call operations

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time

Create mesh for time = 0


PIMPLE: No convergence criteria found


PIMPLE: Operating solver in transient mode with 1 outer corrector
PIMPLE: Operating solver in PISO mode


Reading field p

Reading field 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 cubeRootVol
LES
{
    LESModel        Smagorinsky;
    turbulence      on;
    printCoeffs     on;
    delta           cubeRootVol;
    dynamicKEqnCoeffs
    {
        filter          simple;
    }
    cubeRootVolCoeffs
    {
        deltaCoeff      1;
    }
    PrandtlCoeffs
    {
        delta           cubeRootVol;
        cubeRootVolCoeffs
        {
            deltaCoeff      1;
        }
        smoothCoeffs
        {
            delta           cubeRootVol;
            cubeRootVolCoeffs
            {
                deltaCoeff      1;
            }
            maxDeltaRatio   1.1;
        }
        Cdelta          0.158;
    }
    vanDriestCoeffs
    {
        delta           cubeRootVol;
        cubeRootVolCoeffs
        {
            deltaCoeff      1;
        }
        smoothCoeffs
        {
            delta           cubeRootVol;
            cubeRootVolCoeffs
            {
                deltaCoeff      1;
            }
            maxDeltaRatio   1.1;
        }
        Aplus           26;
        Cdelta          0.158;
    }
    smoothCoeffs
    {
        delta           cubeRootVol;
        cubeRootVolCoeffs
        {
            deltaCoeff      1;
        }
        maxDeltaRatio   1.1;
    }
    Ce              1.048;
    Ck              0.094;
}

No MRF models present

No finite volume options present
Courant Number mean: 0.007375 max: 4.72

Starting time loop

Courant Number mean: 0.007375 max: 4.72
deltaT = 0.01
Time = 0.01

Patch: inlet eddy box:
    volume    : 38.4
    maxSigmaX : 2
'''

Thanks,
Jack

Title: Re: PimpleFoam freezes in first timestep when using turbulentDFSEMInlet
Post by: Jack on September 14, 2021, 11:09:26 PM

Update:
The inlet turbulent generation process seems to be sensitive to the number of cells in the whole domain. I trimmed down the domain to the region near inlet plane and the program works. Now I am planning to store the inlet field and replace the turbulentDFSEMInlet with timeVaryingMappedFixedValue so I can run it on the whole domain.

Title: Re: PimpleFoam freezes in first timestep when using turbulentDFSEMInlet
Post by: ajaybh on September 15, 2021, 09:30:56 PM

Are you running it on single or multiple processor? One issue could be with the domain decomposition. Sometimes, there might not be any cells in one of the domains when partitioned. This could cause an issue for the simulation to stop after the first step. Could you please check this?


Could you also please share the OpenFOAM files?

Title: Re: PimpleFoam freezes in first timestep when using turbulentDFSEMInlet
Post by: Jack on September 30, 2021, 08:05:07 PM

Thanks Ajay. Sorry for the late reply.

I was running it on multiple processors. I retried it with the original full case and it seems to run without any issue now. This may be related to the decomposition you mentioned.

Title: Re: PimpleFoam freezes in first timestep when using turbulentDFSEMInlet
Post by: ajaybh on October 06, 2021, 08:07:16 PM

Sorry, I somehow missed this. Did you try with a new decomposition?