Find answers, ask questions, and connect with our <br>community around the world.

Home Forums OpenFOAM Forum training case-pitzDaily_with_fine mesh

• # training case-pitzDaily_with_fine mesh

Posted by Deleted User on April 30, 2023 at 2:42 pm

Selamlar Hocam,

When I refine the mesh for pitzDaily case it gives the error below. When I look the log file I saw Ux and Uy residuals suddenly diminished. In CFDonline website I found a comment on this and it point out B.C. may be wrong. What can I do about this? I also checked mesh and it did not give any error.( I rearranged it as you said in lecture with doubling numbers.)

Best Regards.

******************************************************************************************************************

Time = 11

smoothSolver: Solving for Ux, Initial residual = 0.00264449, Final residual = 0.00012588, No Iterations 1
smoothSolver: Solving for Uy, Initial residual = 0.00776596, Final residual = 0.00020724, No Iterations 2
GAMG: Solving for p, Initial residual = 0.0525709, Final residual = 7.82874e-13, No Iterations 1
time step continuity errors : sum local = 5.53286e+18, global = -7.01846e+11, cumulative = -7.01399e+11
smoothSolver: Solving for epsilon, Initial residual = 2.68046e-13, Final residual = 2.68046e-13, No Iterations 0
smoothSolver: Solving for k, Initial residual = 3.57079e-10, Final residual = 3.57079e-10, No Iterations 0
ExecutionTime = 5.25 s ClockTime = 5 s

Time = 12

smoothSolver: Solving for Ux, Initial residual = 5.09839e-15, Final residual = 5.09839e-15, No Iterations 0
smoothSolver: Solving for Uy, Initial residual = 4.4968e-15, Final residual = 4.4968e-15, No Iterations 0
GAMG: Solving for p, Initial residual = 0.0911287, Final residual = 0.00630416, No Iterations 3
time step continuity errors : sum local = 7.02927e+32, global = 4.10551e+24, cumulative = 4.10551e+24
smoothSolver: Solving for epsilon, Initial residual = 0.989926, Final residual = 0.0399397, No Iterations 3
bounding epsilon, min: -1.43232e+62 max: 3.32293e+62 average: 8.5483e+57
smoothSolver: Solving for k, Initial residual = 1.47055e-07, Final residual = 1.47055e-07, No Iterations 0
ExecutionTime = 5.64 s ClockTime = 5 s

Time = 13

smoothSolver: Solving for Ux, Initial residual = 3.72233e-10, Final residual = 3.72233e-10, No Iterations 0
smoothSolver: Solving for Uy, Initial residual = 2.59676e-09, Final residual = 2.59676e-09, No Iterations 0
GAMG: Solving for p, Initial residual = 1, Final residual = 4.59526e-07, No Iterations 1
time step continuity errors : sum local = 2.66826e+52, global = -8.12133e+49, cumulative = -8.12133e+49
smoothSolver: Solving for epsilon, Initial residual = 1, Final residual = 0.0490915, No Iterations 1
bounding epsilon, min: -1.31802e+108 max: 7.84859e+108 average: 1.68838e+104
smoothSolver: Solving for k, Initial residual = 0.999015, Final residual = 6.79696e-12, No Iterations 1
bounding k, min: -1.75137e+59 max: 1.20812e+63 average: 5.01688e+58
ExecutionTime = 6.06 s ClockTime = 6 s

Time = 14

smoothSolver: Solving for Ux, Initial residual = 1.39541e-15, Final residual = 1.39541e-15, No Iterations 0
smoothSolver: Solving for Uy, Initial residual = 3.96373e-15, Final residual = 3.96373e-15, No Iterations 0
#0 Foam::error::printStack(Foam::Ostream&) in ~/OpenFOAM/OpenFOAM-v2212/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so
#1 Foam::sigFpe::sigHandler(int) in ~/OpenFOAM/OpenFOAM-v2212/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so
#2 ? in /lib/x86_64-linux-gnu/libc.so.6
#3 Foam::GaussSeidelSmoother::smooth(Foam::word const&, Foam::Field<double>&, Foam::lduMatrix const&, Foam::Field<double> const&, Foam::FieldField<Foam::Field, double> const&, Foam::UPtrList<Foam::lduInterfaceField const> const&, unsigned char, int) in ~/OpenFOAM/OpenFOAM-v2212/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so
#4 Foam::GaussSeidelSmoother::scalarSmooth(Foam::Field<double>&, Foam::Field<double> const&, unsigned char, int) const in ~/OpenFOAM/OpenFOAM-v2212/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so
#5 Foam::GAMGSolver::Vcycle(Foam::PtrList<Foam::lduMatrix::smoother> const&, Foam::Field<double>&, Foam::Field<double> const&, Foam::Field<double>&, Foam::Field<double>&, Foam::Field<double>&, Foam::Field<double>&, Foam::Field<double>&, Foam::PtrList<Foam::Field<double> >&, Foam::PtrList<Foam::Field<double> >&, unsigned char) const in ~/OpenFOAM/OpenFOAM-v2212/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so
#6 Foam::GAMGSolver::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const in ~/OpenFOAM/OpenFOAM-v2212/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so
#7 Foam::fvMatrix<double>::solveSegregated(Foam::dictionary const&) in ~/OpenFOAM/OpenFOAM-v2212/platforms/linux64GccDPInt32Opt/lib/libfiniteVolume.so
#8 Foam::fvMatrix<double>::solveSegregatedOrCoupled(Foam::dictionary const&) in ~/OpenFOAM/OpenFOAM-v2212/platforms/linux64GccDPInt32Opt/lib/libfiniteVolume.so
#9 Foam::fvMesh::solve(Foam::fvMatrix<double>&, Foam::dictionary const&) const in ~/OpenFOAM/OpenFOAM-v2212/platforms/linux64GccDPInt32Opt/lib/libfiniteVolume.so
#10 ? in ~/OpenFOAM/OpenFOAM-v2212/platforms/linux64GccDPInt32Opt/bin/simpleFoam
#11 ? in /lib/x86_64-linux-gnu/libc.so.6
#12 __libc_start_main in /lib/x86_64-linux-gnu/libc.so.6
#13 ? in ~/OpenFOAM/OpenFOAM-v2212/platforms/linux64GccDPInt32Opt/bin/simpleFoam

******************************************************************************************************************

2 Members · 14 Replies
• 14 Replies
• ### BarÄ±ÅŸ Bicer

Moderator
April 30, 2023 at 9:11 pm

Selam,

did you get any error with course mesh setup? if not you only fined the mesh and run same case setup? can you give more details what did you do step by step and got this error. Yeah definitely your BCs are wrong thats why sim. diverges.

can you also post her checkMesh result?

• ### Deleted User

Deleted User
May 1, 2023 at 9:23 am

Yes I run the case coarse mesh without any error, then I rearranged blockMeshDict by doubling mesh division numbers that you showed us in lecture. Then againg run blockMesh command and it gave the above error. I did not change anything except blockMeshDict file.

1)result of checkMesh is Just below
2)my rearranging in blockMeshDict is second

I try to get result of finer mesh analysis in same running folder, is it wrong way?

Thanks.
Best regards.

********************************************************************************************************************
/*—————————————————————————*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 2212 |
| \\ / A nd | Website: http://www.openfoam.com |
| \\/ M anipulation | |
\*—————————————————————————*/
Build : _66908158ae-20221220 OPENFOAM=2212 version=v2212
Arch : “LSB;label=32;scalar=64”
Exec : checkMesh
Date : May 01 2023
Time : 12:19:12
PID : 17266
I/O : uncollated
nProcs : 1
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

Time = 0

Mesh stats
points: 148383
faces: 343420
internal faces: 243380
cells: 97800
faces per cell: 6
boundary patches: 5
point zones: 0
face zones: 0
cell zones: 0

Overall number of cells of each type:
hexahedra: 97800
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 Surface topology
inlet 120 183 ok (non-closed singly connected)
outlet 228 345 ok (non-closed singly connected)
upperWall 892 1341 ok (non-closed singly connected)
lowerWall 1000 1503 ok (non-closed singly connected)
frontAndBack 97800 98922 ok (non-closed singly connected)

Checking faceZone topology for multiply connected surfaces…
No faceZones found.

No cellZones found.

No pointZones found.

Checking geometry…
Overall domain bounding box (-0.0206 -0.0254 -0.0005) (0.29 0.0254 0.0005)
Mesh has 2 geometric (non-empty/wedge) directions (1 1 0)
Mesh has 2 solution (non-empty) directions (1 1 0)
All edges aligned with or perpendicular to non-empty directions.
Boundary openness (-1.24969e-18 -4.04459e-18 -5.1961e-17) OK.
Max cell openness = 3.05871e-16 OK.
Max aspect ratio = 8.20162 OK.
Minimum face area = 4.22924e-08. Maximum face area = 1.28801e-06. Face area magnitudes OK.
Min volume = 2.11462e-11. Max volume = 4.76368e-10. Total volume = 1.4516e-05. Cell volumes OK.
Mesh non-orthogonality Max: 5.96544 average: 1.44923
Non-orthogonality check OK.
Face pyramids OK.
Max skewness = 0.261236 OK.
Coupled point location match (average 0) OK.

Mesh OK.

End

******************************************************************************************************************
blocks
(
hex (0 3 4 1 11 14 15 12)
//(18 30 1)
(36 60 2)

hex (2 5 6 3 13 16 17 14)
//(180 27 1)
(360 54 2)
edgeGrading (4 4 4 4 \$negY 1 1 \$negY 1 1 1 1)

hex (3 6 7 4 14 17 18 15)
//(180 30 1)
(360 60 2)
edgeGrading (4 4 4 4 \$posY \$posYR \$posYR \$posY 1 1 1 1)

hex (5 8 9 6 16 19 20 17)
//(25 27 1)
(50 54 2)

hex (6 9 10 7 17 20 21 18)
//(25 30 1)
(50 60 2)
);

• ### BarÄ±ÅŸ Bicer

Moderator
May 2, 2023 at 8:47 pm

Strange. I didnt see any problem. Can you run case manual and paste here all content of log file please.

• ### Deleted User

Deleted User
May 3, 2023 at 5:08 pm

Sure. It is below.

/*—————————————————————————*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 2212 |
| \\ / A nd | Website: http://www.openfoam.com |
| \\/ M anipulation | |
\*—————————————————————————*/
Build : _66908158ae-20221220 OPENFOAM=2212 version=v2212
Arch : “LSB;label=32;scalar=64”
Exec : simpleFoam
Date : May 03 2023
Time : 20:02:04
PID : 25123
I/O : uncollated
nProcs : 1
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

SIMPLE: convergence criteria
field p tolerance 0.01
field U tolerance 0.001
field “(k|epsilon|omega|f|v2)” tolerance 0.001

Selecting incompressible transport model Newtonian
Selecting turbulence model type RAS
Selecting RAS turbulence model kEpsilon
RAS
{
RASModel kEpsilon;
turbulence on;
printCoeffs on;
Cmu 0.09;
C1 1.44;
C2 1.92;
C3 0;
sigmak 1;
sigmaEps 1.3;
}

No MRF models present

No finite volume options present

Starting time loop

streamLine streamlines:
Employing velocity field U
automatic track length specified through number of sub cycles : 5

Time = 1

smoothSolver: Solving for Ux, Initial residual = 1, Final residual = 0.0416487, No Iterations 2
smoothSolver: Solving for Uy, Initial residual = 1, Final residual = 0.0837667, No Iterations 2
GAMG: Solving for p, Initial residual = 1, Final residual = 0.0944102, No Iterations 32
time step continuity errors : sum local = 1.65198, global = -0.51073, cumulative = -0.51073
smoothSolver: Solving for epsilon, Initial residual = 0.999801, Final residual = 0.0942148, No Iterations 2
bounding epsilon, min: -1.82748e+06 max: 1.46833e+08 average: 296787
smoothSolver: Solving for k, Initial residual = 1, Final residual = 0.000788883, No Iterations 1
bounding k, min: -8.62234 max: 606.951 average: 32.6869
ExecutionTime = 1.26 s ClockTime = 1 s

Time = 2

smoothSolver: Solving for Ux, Initial residual = 0.842241, Final residual = 0.0681018, No Iterations 2
smoothSolver: Solving for Uy, Initial residual = 0.983095, Final residual = 0.067607, No Iterations 2
GAMG: Solving for p, Initial residual = 0.648314, Final residual = 0.0362816, No Iterations 1
time step continuity errors : sum local = 168.958, global = -1.66091, cumulative = -2.17164
smoothSolver: Solving for epsilon, Initial residual = 0.623736, Final residual = 0.00261434, No Iterations 1
bounding epsilon, min: -5.2233e+06 max: 8.73317e+08 average: 4.13364e+06
smoothSolver: Solving for k, Initial residual = 0.629454, Final residual = 0.0484153, No Iterations 3
bounding k, min: -45.8184 max: 27057.3 average: 108.758
ExecutionTime = 1.68 s ClockTime = 2 s

Time = 3

smoothSolver: Solving for Ux, Initial residual = 0.617302, Final residual = 0.0292187, No Iterations 2
smoothSolver: Solving for Uy, Initial residual = 0.824055, Final residual = 0.0319441, No Iterations 2
GAMG: Solving for p, Initial residual = 0.0927089, Final residual = 0.00226491, No Iterations 2
time step continuity errors : sum local = 830.457, global = 1.21319, cumulative = -0.958447
smoothSolver: Solving for epsilon, Initial residual = 0.825569, Final residual = 0.00276052, No Iterations 1
bounding epsilon, min: -1.80314e+08 max: 2.15906e+09 average: 2.15019e+06
smoothSolver: Solving for k, Initial residual = 0.612212, Final residual = 0.0602642, No Iterations 2
bounding k, min: -1759.75 max: 20104.2 average: 128.098
ExecutionTime = 2.07 s ClockTime = 2 s

Time = 4

smoothSolver: Solving for Ux, Initial residual = 0.622281, Final residual = 0.0202116, No Iterations 2
smoothSolver: Solving for Uy, Initial residual = 0.824997, Final residual = 0.0795482, No Iterations 1
GAMG: Solving for p, Initial residual = 0.0790279, Final residual = 0.0048971, No Iterations 1
time step continuity errors : sum local = 1128.96, global = 2.51482, cumulative = 1.55637
smoothSolver: Solving for epsilon, Initial residual = 0.364604, Final residual = 0.0173408, No Iterations 1
bounding epsilon, min: -6.81332e+07 max: 1.09642e+10 average: 3.51909e+06
smoothSolver: Solving for k, Initial residual = 0.647961, Final residual = 0.0460089, No Iterations 1
bounding k, min: -176.667 max: 15233.4 average: 152.512
ExecutionTime = 2.49 s ClockTime = 2 s

Time = 5

smoothSolver: Solving for Ux, Initial residual = 0.591028, Final residual = 0.0271342, No Iterations 2
smoothSolver: Solving for Uy, Initial residual = 0.712959, Final residual = 0.0315978, No Iterations 2
GAMG: Solving for p, Initial residual = 0.073891, Final residual = 0.00724603, No Iterations 1
time step continuity errors : sum local = 1503.5, global = 2.04138, cumulative = 3.59775
smoothSolver: Solving for epsilon, Initial residual = 0.172608, Final residual = 0.0136951, No Iterations 2
bounding epsilon, min: -6.53154e+07 max: 7.88421e+09 average: 3.60464e+06
smoothSolver: Solving for k, Initial residual = 0.468569, Final residual = 0.0308663, No Iterations 2
bounding k, min: -314.157 max: 11628.9 average: 176.361
ExecutionTime = 2.97 s ClockTime = 3 s

Time = 6

smoothSolver: Solving for Ux, Initial residual = 0.516519, Final residual = 0.032357, No Iterations 2
smoothSolver: Solving for Uy, Initial residual = 0.613141, Final residual = 0.0392549, No Iterations 3
GAMG: Solving for p, Initial residual = 0.0810695, Final residual = 0.00308887, No Iterations 2
time step continuity errors : sum local = 630.03, global = 3.01534, cumulative = 6.6131
smoothSolver: Solving for epsilon, Initial residual = 0.5873, Final residual = 0.0131082, No Iterations 2
bounding epsilon, min: -4.73139e+09 max: 1.57686e+10 average: 6.12968e+06
smoothSolver: Solving for k, Initial residual = 0.315194, Final residual = 0.0281303, No Iterations 2
bounding k, min: -1745.15 max: 27997.6 average: 208.695
ExecutionTime = 3.49 s ClockTime = 3 s

Time = 7

smoothSolver: Solving for Ux, Initial residual = 0.495717, Final residual = 0.0256797, No Iterations 1
smoothSolver: Solving for Uy, Initial residual = 0.564044, Final residual = 0.0281303, No Iterations 2
GAMG: Solving for p, Initial residual = 0.0793183, Final residual = 0.00274383, No Iterations 2
time step continuity errors : sum local = 384.698, global = 1.42477, cumulative = 8.03787
smoothSolver: Solving for epsilon, Initial residual = 0.0734618, Final residual = 0.00650051, No Iterations 1
bounding epsilon, min: -1.23568e+09 max: 1.24391e+10 average: 5.90641e+06
smoothSolver: Solving for k, Initial residual = 0.997605, Final residual = 0.0764715, No Iterations 9
bounding k, min: -3908.64 max: 35809.3 average: 283.371
ExecutionTime = 3.96 s ClockTime = 4 s

Time = 8

smoothSolver: Solving for Ux, Initial residual = 1, Final residual = 0.035915, No Iterations 2
smoothSolver: Solving for Uy, Initial residual = 1, Final residual = 0.0328787, No Iterations 2
GAMG: Solving for p, Initial residual = 0.0650991, Final residual = 0.00189336, No Iterations 2
time step continuity errors : sum local = 3.89264e+11, global = 343.936, cumulative = 351.974
smoothSolver: Solving for epsilon, Initial residual = 1, Final residual = 0.0699936, No Iterations 2
bounding epsilon, min: -3.66921e+18 max: 5.69226e+19 average: 1.63666e+16
smoothSolver: Solving for k, Initial residual = 1, Final residual = 0.0133596, No Iterations 2
bounding k, min: -1.99581e+17 max: 4.44301e+18 average: 3.45236e+14
ExecutionTime = 4.46 s ClockTime = 4 s

Time = 9

smoothSolver: Solving for Ux, Initial residual = 0.00691155, Final residual = 0.000673623, No Iterations 2
smoothSolver: Solving for Uy, Initial residual = 0.0241889, Final residual = 0.00164334, No Iterations 2
GAMG: Solving for p, Initial residual = 0.446579, Final residual = 0.0406775, No Iterations 6
time step continuity errors : sum local = 5.67849e+11, global = 3.00112e+08, cumulative = 3.00113e+08
smoothSolver: Solving for epsilon, Initial residual = 5.19515e-06, Final residual = 5.19515e-06, No Iterations 0
smoothSolver: Solving for k, Initial residual = 0.52088, Final residual = 0.0515585, No Iterations 8
bounding k, min: -9.97804e+22 max: 3.91975e+24 average: 5.68961e+21
ExecutionTime = 4.96 s ClockTime = 5 s

Time = 10

smoothSolver: Solving for Ux, Initial residual = 0.0117616, Final residual = 0.000952165, No Iterations 1
smoothSolver: Solving for Uy, Initial residual = 0.0196859, Final residual = 0.00184292, No Iterations 1
GAMG: Solving for p, Initial residual = 1.48649e-08, Final residual = 1.48649e-08, No Iterations 0
time step continuity errors : sum local = 4.10086e+14, global = 1.46481e+08, cumulative = 4.46594e+08
smoothSolver: Solving for epsilon, Initial residual = 0.0024778, Final residual = 5.39404e-05, No Iterations 2
bounding epsilon, min: -2.67538e+24 max: 1.37203e+40 average: 2.08385e+36
smoothSolver: Solving for k, Initial residual = 0.0331729, Final residual = 0.00258841, No Iterations 9
bounding k, min: -2.82918e+19 max: 1.0179e+29 average: 5.07749e+25
ExecutionTime = 5.45 s ClockTime = 5 s

Time = 11

smoothSolver: Solving for Ux, Initial residual = 0.00264449, Final residual = 0.00012588, No Iterations 1
smoothSolver: Solving for Uy, Initial residual = 0.00776596, Final residual = 0.00020724, No Iterations 2
GAMG: Solving for p, Initial residual = 0.0525709, Final residual = 7.82874e-13, No Iterations 1
time step continuity errors : sum local = 5.53286e+18, global = -7.01846e+11, cumulative = -7.01399e+11
smoothSolver: Solving for epsilon, Initial residual = 2.68046e-13, Final residual = 2.68046e-13, No Iterations 0
smoothSolver: Solving for k, Initial residual = 3.57079e-10, Final residual = 3.57079e-10, No Iterations 0
ExecutionTime = 5.87 s ClockTime = 6 s

Time = 12

smoothSolver: Solving for Ux, Initial residual = 5.09839e-15, Final residual = 5.09839e-15, No Iterations 0
smoothSolver: Solving for Uy, Initial residual = 4.4968e-15, Final residual = 4.4968e-15, No Iterations 0
GAMG: Solving for p, Initial residual = 0.0911287, Final residual = 0.00630416, No Iterations 3
time step continuity errors : sum local = 7.02927e+32, global = 4.10551e+24, cumulative = 4.10551e+24
smoothSolver: Solving for epsilon, Initial residual = 0.989926, Final residual = 0.0399397, No Iterations 3
bounding epsilon, min: -1.43232e+62 max: 3.32293e+62 average: 8.5483e+57
smoothSolver: Solving for k, Initial residual = 1.47055e-07, Final residual = 1.47055e-07, No Iterations 0
ExecutionTime = 6.31 s ClockTime = 6 s

Time = 13

smoothSolver: Solving for Ux, Initial residual = 3.72233e-10, Final residual = 3.72233e-10, No Iterations 0
smoothSolver: Solving for Uy, Initial residual = 2.59676e-09, Final residual = 2.59676e-09, No Iterations 0
GAMG: Solving for p, Initial residual = 1, Final residual = 4.59526e-07, No Iterations 1
time step continuity errors : sum local = 2.66826e+52, global = -8.12133e+49, cumulative = -8.12133e+49
smoothSolver: Solving for epsilon, Initial residual = 1, Final residual = 0.0490915, No Iterations 1
bounding epsilon, min: -1.31802e+108 max: 7.84859e+108 average: 1.68838e+104
smoothSolver: Solving for k, Initial residual = 0.999015, Final residual = 6.79696e-12, No Iterations 1
bounding k, min: -1.75137e+59 max: 1.20812e+63 average: 5.01688e+58
ExecutionTime = 6.86 s ClockTime = 7 s

Time = 14

smoothSolver: Solving for Ux, Initial residual = 1.39541e-15, Final residual = 1.39541e-15, No Iterations 0
smoothSolver: Solving for Uy, Initial residual = 3.96373e-15, Final residual = 3.96373e-15, No Iterations 0

• ### BarÄ±ÅŸ Bicer

Moderator
May 3, 2023 at 9:25 pm

ok your initial condition looks problematic for turbulence paramaters. Can you run case laminar without turbulence. if it works then you should check the initial conditions for k and epsilon

• ### Deleted User

Deleted User
May 4, 2023 at 12:57 pm

Selamlar Hocam,

I ran the case as laminar and it finished succesfully, but value of residuals is not less as coarse mesh.(When coarse mesh-turbulence residuals is about .e-5, laminar case has value about .e-3).
When I checked turbulence value as you said and it was same and there is no any change (already I havent any changes already)

end of laminar log file is below.

******************************************************************************************************************

/*—————————————————————————*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 2212 |
| \\ / A nd | Website: http://www.openfoam.com |
| \\/ M anipulation | |
\*—————————————————————————*/
Build : _66908158ae-20221220 OPENFOAM=2212 version=v2212
Arch : “LSB;label=32;scalar=64”
Exec : simpleFoam
Date : May 04 2023
Time : 13:24:01
PID : 29179
I/O : uncollated
nProcs : 1
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

SIMPLE: convergence criteria
field p tolerance 0.01
field U tolerance 0.001
field “(k|epsilon|omega|f|v2)” tolerance 0.001

Selecting incompressible transport model Newtonian
Selecting turbulence model type laminar
Selecting laminar stress model Stokes
No MRF models present

No finite volume options present

Starting time loop

streamLine streamlines:
Employing velocity field U
automatic track length specified through number of sub cycles : 5

Time = 1

smoothSolver: Solving for Ux, Initial residual = 1, Final residual = 0.00074411, No Iterations 1
smoothSolver: Solving for Uy, Initial residual = 0.290984, Final residual = 8.63701e-05, No Iterations 1
GAMG: Solving for p, Initial residual = 1, Final residual = 0.0973048, No Iterations 32
time step continuity errors : sum local = 1.70263, global = -0.525404, cumulative = -0.525404

.
.
.
.
.
.
.
ExecutionTime = 1.08 s ClockTime = 1 s
Time = 1997

smoothSolver: Solving for Ux, Initial residual = 0.0353998, Final residual = 0.00220838, No Iterations 4
smoothSolver: Solving for Uy, Initial residual = 0.0782125, Final residual = 0.00582116, No Iterations 4
GAMG: Solving for p, Initial residual = 0.0950265, Final residual = 0.00876123, No Iterations 2
time step continuity errors : sum local = 7.40133, global = -1.2019, cumulative = -876.347
ExecutionTime = 469.88 s ClockTime = 472 s

Time = 1998

smoothSolver: Solving for Ux, Initial residual = 0.0351596, Final residual = 0.00218638, No Iterations 4
smoothSolver: Solving for Uy, Initial residual = 0.0771471, Final residual = 0.00573258, No Iterations 4
GAMG: Solving for p, Initial residual = 0.0942146, Final residual = 0.00861078, No Iterations 2
time step continuity errors : sum local = 7.40341, global = -1.25315, cumulative = -877.6
ExecutionTime = 470.05 s ClockTime = 472 s

Time = 1999

smoothSolver: Solving for Ux, Initial residual = 0.0346741, Final residual = 0.00215438, No Iterations 4
smoothSolver: Solving for Uy, Initial residual = 0.0763012, Final residual = 0.00562262, No Iterations 4
GAMG: Solving for p, Initial residual = 0.0926336, Final residual = 0.00841476, No Iterations 2
time step continuity errors : sum local = 7.42033, global = -1.28666, cumulative = -878.887
ExecutionTime = 470.23 s ClockTime = 472 s

Time = 2000

smoothSolver: Solving for Ux, Initial residual = 0.0344585, Final residual = 0.00213323, No Iterations 4
smoothSolver: Solving for Uy, Initial residual = 0.0750271, Final residual = 0.00547603, No Iterations 4
GAMG: Solving for p, Initial residual = 0.0894754, Final residual = 0.00816894, No Iterations 2
time step continuity errors : sum local = 7.39083, global = -1.29565, cumulative = -880.182
ExecutionTime = 470.63 s ClockTime = 473 s

streamLine streamlines write:
seeded 10 particles
Tracks:10
Total samples:19642
End

• ### Deleted User

Deleted User
May 4, 2023 at 1:02 pm

Turbulence files ss is below.

• ### Deleted User

Deleted User
May 4, 2023 at 1:03 pm

Turbulence files ss is below.

• ### Deleted User

Deleted User
May 4, 2023 at 1:03 pm

• ### BarÄ±ÅŸ Bicer

Moderator
May 4, 2023 at 4:22 pm

selamlar,

I copied the case and run it as coarse mesh. I didnt get any error with turbulence model. Also I checked your blockMesh where you divided Z directhion which is empty(since the case is 2D) into 2. It should be 1. not 2. Recareate your mesh again and run the case please.

• ### Deleted User

Deleted User
May 4, 2023 at 6:21 pm

Ok, Hocam.

Thanks.

• ### BarÄ±ÅŸ Bicer

Moderator
May 4, 2023 at 8:06 pm

welcome

• ### Deleted User

Deleted User
May 4, 2023 at 8:45 pm

You are right Hocam, as you said I change third dimension as 1 in blockMeshDict because of the case 2D and it worked and converged.

Thanks again.

Best regards.

• ### BarÄ±ÅŸ Bicer

Moderator
May 5, 2023 at 10:29 am

good to know that.