PimpleFoam freezes in first timestep when using turbulentDFSEMInlet - Jack, September 14, 2021

[JustinBonus]

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

• Jack, September 14, 2021

[JustinBonus]

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.

• Jack, September 14, 2021

[JustinBonus]

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?

• ajaybh, September 15, 2021

[JustinBonus]

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.

• Jack, September 30, 2021

[JustinBonus]

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

• ajaybh, October 06, 2021