Hello Everyone

May you kindly assist me,

I am using viscoelasticFluidFoam solver  to simulate a flow in a pipe. And have modified it by adding temperature.

Now every time I run my simulation it gives me the below error and have tried to reduce delta T, Max Co. number and even check velocity value but nothing seems to be working

Courant Number mean: 8.19815e-06 max: 0.0016354 velocity magnitude: 0.336024
deltaT = 0.000121669
Time = 0.000730166

DILUPBiCG: Solving for Ux, Initial residual = 0.00475588, Final residual = 9.01518e-08, No Iterations 1
DILUPBiCG: Solving for Uy, Initial residual = 0.00477491, Final residual = 5.60245e-09, No Iterations 1
DILUPBiCG: Solving for Uz, Initial residual = 0.00153646, Final residual = 3.99944e-09, No Iterations 1
DICPCG: Solving for p, Initial residual = 0.202343, Final residual = 9.65702e-08, No Iterations 322
time step continuity errors : sum local = 1.74596e-13, global = -1.10905e-14, cumulative = -1.43122e-14
Model mode 1
DILUPBiCG: Solving for taufirstxx, Initial residual = 0.00504237, Final residual = 2.64016e-15, No Iterations 1
DILUPBiCG: Solving for taufirstxy, Initial residual = 0.00719778, Final residual = 1.69566e-15, No Iterations 1
DILUPBiCG: Solving for taufirstxz, Initial residual = 0.00372445, Final residual = 1.81474e-17, No Iterations 1
DILUPBiCG: Solving for taufirstyy, Initial residual = 0.00513511, Final residual = 2.97488e-15, No Iterations 1
DILUPBiCG: Solving for taufirstyz, Initial residual = 0.00372444, Final residual = 2.90091e-18, No Iterations 1
DILUPBiCG: Solving for taufirstzz, Initial residual = 0.00360069, Final residual = 9.85626e-17, No Iterations 1
Model mode 2
DILUPBiCG: Solving for tausecondxx, Initial residual = 0.00901374, Final residual = 4.5788e-16, No Iterations 1
DILUPBiCG: Solving for tausecondxy, Initial residual = 0.011415, Final residual = 3.71279e-16, No Iterations 1
DILUPBiCG: Solving for tausecondxz, Initial residual = 0.00373582, Final residual = 1.74955e-17, No Iterations 1
DILUPBiCG: Solving for tausecondyy, Initial residual = 0.00924354, Final residual = 6.99518e-16, No Iterations 1
DILUPBiCG: Solving for tausecondyz, Initial residual = 0.00373581, Final residual = 1.75101e-17, No Iterations 1
DILUPBiCG: Solving for tausecondzz, Initial residual = 0.00294321, Final residual = 5.88967e-17, No Iterations 1
Model mode 3
DILUPBiCG: Solving for tauthirdxx, Initial residual = 0.00929762, Final residual = 1.30313e-16, No Iterations 1
DILUPBiCG: Solving for tauthirdxy, Initial residual = 0.0117061, Final residual = 4.89053e-17, No Iterations 1
DILUPBiCG: Solving for tauthirdxz, Initial residual = 0.00373792, Final residual = 1.78278e-17, No Iterations 1
DILUPBiCG: Solving for tauthirdyy, Initial residual = 0.00953793, Final residual = 1.59044e-16, No Iterations 1
DILUPBiCG: Solving for tauthirdyz, Initial residual = 0.00373792, Final residual = 1.9548e-17, No Iterations 1
DILUPBiCG: Solving for tauthirdzz, Initial residual = 0.0026507, Final residual = 5.17534e-18, No Iterations 1
Model mode 4
DILUPBiCG: Solving for taufourthxx, Initial residual = 0.00940334, Final residual = 3.38134e-17, No Iterations 1
DILUPBiCG: Solving for taufourthxy, Initial residual = 0.0118159, Final residual = 1.53916e-17, No Iterations 1
DILUPBiCG: Solving for taufourthxz, Initial residual = 0.00373861, Final residual = 7.48391e-18, No Iterations 1
DILUPBiCG: Solving for taufourthyy, Initial residual = 0.00964937, Final residual = 8.37758e-17, No Iterations 1
DILUPBiCG: Solving for taufourthyz, Initial residual = 0.00373861, Final residual = 1.28578e-17, No Iterations 1
DILUPBiCG: Solving for taufourthzz, Initial residual = 0.00246207, Final residual = 8.77819e-17, No Iterations 1
DILUPBiCG: Solving for T, Initial residual = 0.45522, Final residual = 7.40274e-09, No Iterations 5
DILUPBiCG: Solving for Tw, Initial residual = 0.0828233, Final residual = 1.66062e-08, No Iterations 4
Floating point exception (core dumped)

CHECKMESH :

Create time

Create polyMesh for time = 0

Time = 0

Mesh stats
all points: 222101
live points: 222101
all faces: 654000
live faces: 654000
internal faces: 642000
cells: 216000
boundary patches: 3
point zones: 0
face zones: 0
cell zones: 0

Overall number of cells of each type:
hexahedra: 216000
prisms: 0
wedges: 0
pyramids: 0
tet wedges: 0
tetrahedra: 0
polyhedra: 0

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 Area [m^2] Surface topology
inlet 3600 3641 0.00196148 ok (non-closed singly connected)
outlet 3600 3641 0.00196148 ok (non-closed singly connected)
pipewall 4800 4880 0.392598 ok (non-closed singly connected)

Checking geometry...
This is a 3-D mesh
Overall domain bounding box (-0.025 -0.025 0) (0.025 0.025 2.5)
Mesh (non-empty, non-wedge) directions (1 1 1)
Mesh (non-empty) directions (1 1 1)
Mesh (non-empty, non-wedge) dimensions 3
Boundary openness (2.02337e-16 2.16504e-16 4.25812e-17) Threshold = 1e-06 OK.
Max cell openness = 2.74096e-16 OK.
Max aspect ratio = 195.315 OK.
Minumum face area = 3.67922e-07. Maximum face area = 8.17913e-05. Face area magnitudes OK.
Min volume = 1.53301e-08. Max volume = 2.95413e-08. Total volume = 0.00490369. Cell volumes OK.
Mesh non-orthogonality Max: 37.1451 average: 5.71625 Threshold = 70
Non-orthogonality check OK.
Face pyramids OK.
Max skewness = 0.569773 OK.

Mesh OK.

CONTROL DICT:
startFrom startTime;

startTime 0.0;

stopAt endTime;

endTime 2;

deltaT 1e-16;

writeControl adjustableRunTime;

writeInterval 1;

purgeWrite 0;

writeFormat ascii;

writePrecision 6;

writeCompression uncompressed;

timeFormat general;

timePrecision 6;

graphFormat raw;

runTimeModifiable yes;

adjustTimeStep on;

maxCo 0.25;

maxDeltaT 0.001;


I am not sure on how else to resolve this issue . May you kindly assist me or provide some guidance on how to resolve this problem


Thanks in advance