[Virtual-iCFD]

Hi everyone,

Below is my setup for controlDict to run interDyMFoam.

My concern is now (after 10 iterations) the deltaT is in the order of e-05 and continue to drops. I have tried to tighten pcorr to 1e-12 and with Gauss upwind for convection and also try Gauss linear limited 0.7 for laplacian. Also, decrease underelaxation for k and epsilon to 0.5. Still does not see the light. In this case, I also do the checkMesh and all look okay as indicated below. My PISO setting for Courant no. is also provided below. Considered this stage, my concern is time is not marching any way due to extremely small deltaT.

Any suggestion will be greatly appreciated.

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

application interDyMFoam;

startFrom latestTime;

startTime 0;

stopAt endTime;

endTime 0.1;

deltaT 0.005;

writeControl adjustableRunTime;

writeInterval 0.01;

purgeWrite 0;

writeFormat ascii;

writePrecision 12;

writeCompression uncompressed;

timeFormat general;

timePrecision 10;

runTimeModifiable yes;

adjustTimeStep yes;

maxCo 0.5;
maxAlphaCo 0.5;
maxDeltaT 0.005;

libs
(
"libincompressibleRASModels.so"
"libfvMotionSolvers.so"
"libforces.so"
);

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

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

Create polyMesh for time = 0

Time = 0

Mesh stats
points: 599829
faces: 1728840
internal faces: 1677599
cells: 565109
boundary patches: 6
point zones: 0
face zones: 0
cell zones: 0

Overall number of cells of each type:
hexahedra: 556015
prisms: 629
wedges: 0
pyramids: 0
tet wedges: 20
tetrahedra: 0
polyhedra: 8445

Checking topology...
Boundary definition OK.
Point usage OK.
Upper triangular ordering OK.
Face vertices OK.
Number of regions: 1 (OK).

Checking patch topology for multiply connected surfaces ...
Patch Faces Points Surface topology
sidewall 16384 16770 ok (non-closed singly connected)
inlet 4096 4225 ok (non-closed singly connected)
outlet 4096 4225 ok (non-closed singly connected)
floor 8192 8385 ok (non-closed singly connected)
ceiling 8192 8385 ok (non-closed singly connected)
floatingObject 10281 11433 ok (closed singly connected)

Checking geometry...
Overall domain bounding box (-2 -1.5 -2) (4 2 2)
Mesh (non-empty, non-wedge) directions (1 1 1)
Mesh (non-empty) directions (1 1 1)
Boundary openness (2.96151165684e-16 4.64197183787e-17 9.36902725343e-15) OK.
Max cell openness = 2.53757529159e-16 OK.
Max aspect ratio = 4.49379234381 OK.
Minumum face area = 2.09352045e-05. Maximum face area = 0.00341828125. Face area magnitudes OK.
Min volume = 4.3804331189e-07. Max volume = 0.0001602417. Total volume = 83.9563743795. Cell volumes OK.
Mesh non-orthogonality Max: 51.5433563696 average: 3.4337367223
Non-orthogonality check OK.
Face pyramids OK.
Max skewness = 3.44827063315 OK.

Mesh OK.

End


#------------

PISO
{
momentumPredictor no;
nCorrectors 3;
nNonOrthogonalCorrectors 2;
nAlphaCorr 1;
nAlphaSubCycles 3;
cAlpha 1;
correctPhi yes;
}