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

Home Forums OpenFOAM Forum golf ball simulation

• # golf ball simulation

Posted by on May 20, 2022 at 10:44 pm

I did same simulation but my simulation is not converged .

————————-error —————————————-
Time = 1

smoothSolver: Â Solving for Ux, Initial residual = 1, Final residual = 0.00375137, No Iterations 1
smoothSolver: Â Solving for Uy, Initial residual = 0.999995, Final residual = 0.00293765, No Iterations 1
smoothSolver: Â Solving for Uz, Initial residual = 0.999993, Final residual = 0.00143062, No Iterations 1
GAMG: Â Solving for p, Initial residual = 1, Final residual = 0.032254, No Iterations 2
time step continuity errors : sum local = 0.161273, global = -0.137299, cumulative = -0.137299
smoothSolver: Â Solving for epsilon, Initial residual = 0.0569919, Final residual = 0.00552814, No Iterations 1
smoothSolver: Â Solving for k, Initial residual = 1, Final residual = 0.0993227, No Iterations 1
ExecutionTime = 7.79 s Â ClockTime = 8 s

Time = 2

smoothSolver: Â Solving for Ux, Initial residual = 0.273588, Final residual = 0.0259058, No Iterations 1
smoothSolver: Â Solving for Uy, Initial residual = 0.270249, Final residual = 0.0149572, No Iterations 1
smoothSolver: Â Solving for Uz, Initial residual = 0.269659, Final residual = 0.0150405, No Iterations 1
GAMG: Â Solving for p, Initial residual = 0.540361, Final residual = 0.0169022, No Iterations 2
time step continuity errors : sum local = 0.109318, global = -0.091879, cumulative = -0.229178
smoothSolver: Â Solving for epsilon, Initial residual = 0.000826276, Final residual = 4.55167e-05, No Iterations 1
bounding epsilon, min: -1.72876 max: 107343 average: 1576.14
smoothSolver: Â Solving for k, Initial residual = 0.539268, Final residual = 0.037327, No Iterations 1
ExecutionTime = 13.41 s Â ClockTime = 14 s

Time = 3

smoothSolver: Â Solving for Ux, Initial residual = 1, Final residual = 0.013802, No Iterations 1
smoothSolver: Â Solving for Uy, Initial residual = 1, Final residual = 0.013802, No Iterations 1
smoothSolver: Â Solving for Uz, Initial residual = 1, Final residual = 0.013802, No Iterations 1
GAMG: Â Solving for p, Initial residual = 1, Final residual = 0.0263338, No Iterations 1
time step continuity errors : sum local = 8.07886e+09, global = 2.74335e+09, cumulative = 2.74335e+09
smoothSolver: Â Solving for epsilon, Initial residual = 1, Final residual = 0.0187136, No Iterations 1
bounding epsilon, min: -4.53418e+11 max: 1.32632e+14 average: 4.03721e+09
smoothSolver: Â Solving for k, Initial residual = 1, Final residual = 0.0141131, No Iterations 1
bounding k, min: -6.27842e+06 max: 1.65333e+13 average: 1.788e+08
ExecutionTime = 18.88 s Â ClockTime = 20 s

Time = 4

smoothSolver: Â Solving for Ux, Initial residual = 1, Final residual = 0.0242878, No Iterations 1
smoothSolver: Â Solving for Uy, Initial residual = 1, Final residual = 0.0246007, No Iterations 1
smoothSolver: Â Solving for Uz, Initial residual = 1, Final residual = 0.0224357, No Iterations 1
GAMG: Â Solving for p, Initial residual = 1, Final residual = 0.0930323, No Iterations 1
time step continuity errors : sum local = 1.66715e+26, global = 5.26331e+25, cumulative = 5.26331e+25
smoothSolver: Â Solving for epsilon, Initial residual = 1, Final residual = 0.0958824, No Iterations 1
bounding epsilon, min: -3.2573e+36 max: 1.14152e+37 average: 1.49605e+32
smoothSolver: Â Solving for k, Initial residual = 1, Final residual = 0.0283157, No Iterations 1
bounding k, min: -1.7726e+35 max: 1.49759e+36 average: 2.29883e+31
ExecutionTime = 24.17 s Â ClockTime = 25 s

Time = 5

smoothSolver: Â Solving for Ux, Initial residual = 0.999998, Final residual = 0.0465079, No Iterations 3
smoothSolver: Â Solving for Uy, Initial residual = 0.999997, Final residual = 0.0465078, No Iterations 3
smoothSolver: Â Solving for Uz, Initial residual = 1, Final residual = 0.046508, No Iterations 3
GAMG: Â Solving for p, Initial residual = 1, Final residual = 0.09729, No Iterations 1
time step continuity errors : sum local = 1.11137e+64, global = 3.89077e+63, cumulative = 3.89077e+63
smoothSolver: Â Solving for epsilon, Initial residual = 1, Final residual = 0.0192271, No Iterations 2
bounding epsilon, min: -2.15672e+101 max: 7.06398e+102 average: 7.54211e+97
smoothSolver: Â Solving for k, Initial residual = 1, Final residual = 0.0288333, No Iterations 1
bounding k, min: -3.28193e+98 max: 2.34855e+101 average: 1.13643e+96
ExecutionTime = 30.16 s Â ClockTime = 32 s

Time = 6

smoothSolver: Â Solving for Ux, Initial residual = 0.999998, Final residual = 0.0783099, No Iterations 3
smoothSolver: Â Solving for Uy, Initial residual = 0.999994, Final residual = 0.0786014, No Iterations 3
smoothSolver: Â Solving for Uz, Initial residual = 0.999964, Final residual = 0.0550584, No Iterations 3
[0] [1] #0 Â Foam::error::printStack(Foam::Ostream&)#0 Â Foam::error::printStack(Foam::Ostream&) at ??:?
[1] #1 Â Foam::sigFpe::sigHandler(int) at ??:?
[0] #1 Â Foam::sigFpe::sigHandler(int) at ??:?
[1] #2 Â ? at ??:?
[0] #2 Â ? in /lib/x86_64-linux-gnu/libpthread.so.0
[1] #3 Â Foam::scalarProduct<double, double>::type Foam::sumProd<double>(Foam::UList<double> const&, Foam::UList<double> const&) in /lib/x86_64-linux-gnu/libpthread.so.0
[0] #3 Â Foam::scalarProduct<double, double>::type Foam::sumProd<double>(Foam::UList<double> const&, Foam::UList<double> const&) at ??:?
[1] #4 Â Foam::PCG::scalarSolve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const at ??:?
[0] #4 Â Foam::PCG::scalarSolve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const at ??:?
[1] #5 Â Foam::GAMGSolver::solveCoarsestLevel(Foam::Field<double>&, Foam::Field<double> const&) const at ??:?
[0] #5 Â Foam::GAMGSolver::solveCoarsestLevel(Foam::Field<double>&, Foam::Field<double> const&) const at ??:?
[1] #6 Â 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 at ??:?
[0] #6 Â 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 at ??:?
[1] #7 Â Foam::GAMGSolver::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const at ??:?
[0] #7 Â Foam::GAMGSolver::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const at ??:?
[0] #8 Â Foam::fvMatrix<double>::solveSegregated(Foam::dictionary const&) at ??:?
[1] #8 Â Foam::fvMatrix<double>::solveSegregated(Foam::dictionary const&) at ??:?
[1] #9 Â Foam::fvMatrix<double>::solveSegregatedOrCoupled(Foam::dictionary const&) at ??:?
[0] #9 Â Foam::fvMatrix<double>::solveSegregatedOrCoupled(Foam::dictionary const&) at ??:?
[1] #10 Â Foam::fvMesh::solve(Foam::fvMatrix<double>&, Foam::dictionary const&) const at ??:?
[0] #10 Â Foam::fvMesh::solve(Foam::fvMatrix<double>&, Foam::dictionary const&) const at ??:?
[0] #11 Â ? in /usr/lib/openfoam/openfoam2112/platforms/linux64GccDPInt32Opt/bin/simpleFoam
[0] #12 Â __libc_start_main at ??:?
[1] #11 Â ? in /lib/x86_64-linux-gnu/libc.so.6
[0] #13 Â  in /usr/lib/openfoam/openfoam2112/platforms/linux64GccDPInt32Opt/bin/simpleFoam
[1] #12 Â __libc_start_main? in /lib/x86_64-linux-gnu/libc.so.6
[1] #13 Â  in /usr/lib/openfoam/openfoam2112/platforms/linux64GccDPInt32Opt/bin/simpleFoam
[shivam-VirtualBox:03759] *** Process received signal ***
[shivam-VirtualBox:03759] Signal: Floating point exception (8)
[shivam-VirtualBox:03759] Signal code: Â (-6)
[shivam-VirtualBox:03759] [ 3] /usr/lib/openfoam/openfoam2112/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so(_ZN4Foam7sumProdIdEENS_13scalarProductIT_S2_E4typeERKNS_5UListIS2_EES8_+0x35)[0x7fba66b736a5]
[shivam-VirtualBox:03759] [ 4] ?/usr/lib/openfoam/openfoam2112/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so(_ZNK4Foam3PCG11scalarSolveERNS_5FieldIdEERKS2_h+0x857)[0x7fba66966cb7]
[shivam-VirtualBox:03759] [ 5] /usr/lib/openfoam/openfoam2112/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so(_ZNK4Foam10GAMGSolver18solveCoarsestLevelERNS_5FieldIdEERKS2_+0xbf)[0x7fba66990c1f]
[shivam-VirtualBox:03759] [ 6] /usr/lib/openfoam/openfoam2112/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so(_ZNK4Foam10GAMGSolver6VcycleERKNS_7PtrListINS_9lduMatrix8smootherEEERNS_5FieldIdEERKS8_S9_S9_S9_S9_S9_RNS1_IS8_EESD_h+0x620)[0x7fba66993580]
[shivam-VirtualBox:03759] [ 7] /usr/lib/openfoam/openfoam2112/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so(_ZNK4Foam10GAMGSolver5solveERNS_5FieldIdEERKS2_h+0x625)[0x7fba66995c15]
[shivam-VirtualBox:03759] [ 8] /usr/lib/openfoam/openfoam2112/platforms/linux64GccDPInt32Opt/lib/libfiniteVolume.so(_ZN4Foam8fvMatrixIdE15solveSegregatedERKNS_10dictionaryE+0x5c6)[0x7fba6aa73fd6]
[shivam-VirtualBox:03759] [ 9] /usr/lib/openfoam/openfoam2112/platforms/linux64GccDPInt32Opt/lib/libfiniteVolume.so(_ZN4Foam8fvMatrixIdE24solveSegregatedOrCoupledERKNS_10dictionaryE+0x3f8)[0x7fba6a2b4f78]
[shivam-VirtualBox:03759] [10] /usr/lib/openfoam/openfoam2112/platforms/linux64GccDPInt32Opt/lib/libfiniteVolume.so(_ZNK4Foam6fvMesh5solveERNS_8fvMatrixIdEERKNS_10dictionaryE+0x23)[0x7fba6a256ee3]
[shivam-VirtualBox:03759] [11] simpleFoam(+0x27a05)[0x5648bc8e3a05]
[shivam-VirtualBox:03759] [12] /lib/x86_64-linux-gnu/libc.so.6(__libc_start_main+0xe7)[0x7fba65169c87]
[shivam-VirtualBox:03759] [13] simpleFoam(+0x29c9a)[0x5648bc8e5c9a]
[shivam-VirtualBox:03759] *** End of error message ***
Â in /usr/lib/openfoam/openfoam2112/platforms/linux64GccDPInt32Opt/bin/simpleFoam
[shivam-VirtualBox:03760] *** Process received signal ***
[shivam-VirtualBox:03760] Signal: Floating point exception (8)
[shivam-VirtualBox:03760] Signal code: Â (-6)
[shivam-VirtualBox:03760] [ 3] /usr/lib/openfoam/openfoam2112/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so(_ZN4Foam7sumProdIdEENS_13scalarProductIT_S2_E4typeERKNS_5UListIS2_EES8_+0x35)[0x7f9ef84d36a5]
[shivam-VirtualBox:03760] [ 4] /usr/lib/openfoam/openfoam2112/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so(_ZNK4Foam3PCG11scalarSolveERNS_5FieldIdEERKS2_h+0x857)[0x7f9ef82c6cb7]
[shivam-VirtualBox:03760] [ 5] /usr/lib/openfoam/openfoam2112/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so(_ZNK4Foam10GAMGSolver18solveCoarsestLevelERNS_5FieldIdEERKS2_+0xbf)[0x7f9ef82f0c1f]
[shivam-VirtualBox:03760] [ 6] /usr/lib/openfoam/openfoam2112/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so(_ZNK4Foam10GAMGSolver6VcycleERKNS_7PtrListINS_9lduMatrix8smootherEEERNS_5FieldIdEERKS8_S9_S9_S9_S9_S9_RNS1_IS8_EESD_h+0x620)[0x7f9ef82f3580]
[shivam-VirtualBox:03760] [ 7] /usr/lib/openfoam/openfoam2112/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so(_ZNK4Foam10GAMGSolver5solveERNS_5FieldIdEERKS2_h+0x625)[0x7f9ef82f5c15]
[shivam-VirtualBox:03760] [ 8] /usr/lib/openfoam/openfoam2112/platforms/linux64GccDPInt32Opt/lib/libfiniteVolume.so(_ZN4Foam8fvMatrixIdE15solveSegregatedERKNS_10dictionaryE+0x5c6)[0x7f9efc3d3fd6]
[shivam-VirtualBox:03760] [ 9] /usr/lib/openfoam/openfoam2112/platforms/linux64GccDPInt32Opt/lib/libfiniteVolume.so(_ZN4Foam8fvMatrixIdE24solveSegregatedOrCoupledERKNS_10dictionaryE+0x3f8)[0x7f9efbc14f78]
[shivam-VirtualBox:03760] [10] /usr/lib/openfoam/openfoam2112/platforms/linux64GccDPInt32Opt/lib/libfiniteVolume.so(_ZNK4Foam6fvMesh5solveERNS_8fvMatrixIdEERKNS_10dictionaryE+0x23)[0x7f9efbbb6ee3]
[shivam-VirtualBox:03760] [11] simpleFoam(+0x27a05)[0x565378f8ca05]
[shivam-VirtualBox:03760] [12] /lib/x86_64-linux-gnu/libc.so.6(__libc_start_main+0xe7)[0x7f9ef6ac9c87]
[shivam-VirtualBox:03760] [13] simpleFoam(+0x29c9a)[0x565378f8ec9a]
[shivam-VirtualBox:03760] *** End of error message ***
————————————————————————–
mpirun noticed that process rank 1 with PID 0 on node shivam-VirtualBox exited on signal 8 (Floating point exception).
———————————error end—————————————————

if you can see after time 4 my k value is not stable .could you please tell me what will be the possible reason .

replied 2 years ago 3 Members · 17 Replies
• 17 Replies
• ### BarÄ±ÅŸ Bicer

Moderator
May 22, 2022 at 2:57 pm

Hi Shivam,

Your simulation diverges due to bounding k which leads to floating point error.

I think that your initial conditions for turbulence parameters (k and epsilon) are wrong posed and most probably assigned as 0.

Could you share these k and epsilon files located in /0 folder.

Best.

BarÄ±ÅŸ

• ### shivam.shivamchoubey

Member
May 23, 2022 at 8:58 am

/*——————————–*- C++ -*———————————-*\
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration | Website: https://openfoam.org
\\ / A nd | Version: 7
\\/ M anipulation |
\*—————————————————————————*/
FoamFile
{
version 2.0;
format ascii;
class volScalarField;
location “0”;
object k;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions [0 2 -2 0 0 0 0];

internalField uniform 1.5;

boundaryField
{
inlet
{
type fixedValue;
value \$internalField;
}

farfield
{

type inletOutlet;
inletValue \$internalField;
value \$internalField;
}

golfBall
{
type kqRWallFunction;
value \$internalField;
}

}

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

• ### shivam.shivamchoubey

Member
May 23, 2022 at 8:58 am

/*——————————–*- C++ -*———————————-*\
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration | Website: https://openfoam.org
\\ / A nd | Version: 7
\\/ M anipulation |
\*—————————————————————————*/
FoamFile
{
version 2.0;
format ascii;
class volScalarField;
location “0”;
object epsilon;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions [0 2 -3 0 0 0 0];

internalField uniform 4.31;

boundaryField
{
inlet
{
type fixedValue;
value \$internalField;
}

farfield
{

type inletOutlet;
inletValue \$internalField;
value \$internalField;
}

golfBall
{
type epsilonWallFunction;
value \$internalField;
}

}

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

• ### BarÄ±ÅŸ Bicer

Moderator
May 23, 2022 at 9:11 am

Hi again,

Did you change anything from the original setup what we provided for this case?

Your mesh is similar with us? can you try checkMesh utility and paste the result here please.

Thank you.

BarÄ±ÅŸ

• ### shivam.shivamchoubey

Member
May 23, 2022 at 10:18 am

I change nothing just ran snappyHexMesh.

———————-this is checkmesh output—————————
/*—————————————————————————*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 2112 |
| \\ / A nd | Website: http://www.openfoam.com |
| \\/ M anipulation | |
\*—————————————————————————*/
Build : _80318542-20220301 OPENFOAM=2112 patch=220310 version=2112
Arch : “LSB;label=32;scalar=64”
Exec : checkMesh
Date : May 23 2022
Time : 12:13:48
Host : shivam-VirtualBox
PID : 2189
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: 633204
faces: 1763758
internal faces: 1729893
cells: 569565
faces per cell: 6.13389
boundary patches: 3
point zones: 0
face zones: 0
cell zones: 0

Overall number of cells of each type:
hexahedra: 505003
prisms: 7104
wedges: 1030
pyramids: 0
tet wedges: 1903
tetrahedra: 16
polyhedra: 54509
Breakdown of polyhedra by number of faces:
faces number of cells
4 7526
5 6691
6 9875
7 5214
8 1865
9 16502
12 3810
15 2892
18 134

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 49 64 ok (non-closed singly connected)
farfield 1505 1564 ok (non-closed singly connected)
golfBall 32311 40421 ok (closed singly connected)

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

No cellZones found.

Checking geometry…
Overall domain bounding box (-0.99 -0.350001 -0.350002) (3.01 0.350002 0.350002)
Mesh has 3 geometric (non-empty/wedge) directions (1 1 1)
Mesh has 3 solution (non-empty) directions (1 1 1)
Boundary openness (3.7569e-17 -6.59073e-17 1.18903e-16) OK.
Max cell openness = 7.56116e-16 OK.
Max aspect ratio = 20.2217 OK.
Minimum face area = 3.31116e-09. Maximum face area = 0.0100045. Face area magnitudes OK.
Min volume = 7.12952e-13. Max volume = 0.00100027. Total volume = 1.95996. Cell volumes OK.
Mesh non-orthogonality Max: 64.8836 average: 8.3943
Non-orthogonality check OK.
Face pyramids OK.
Max skewness = 0.988576 OK.
Coupled point location match (average 0) OK.

Mesh OK.

End

• ### BarÄ±ÅŸ Bicer

Moderator
May 23, 2022 at 1:40 pm

Hi shivam again,

did you run case with a single CPU or parallel? Can you run the case with a single CPU and write output into log file and pls share with me.

thank you

• ### shivam.shivamchoubey

Member
May 23, 2022 at 2:13 pm

Hello Mr Bicer,

Thank you so much for the fast response. I found more things in the 0/U field value are saved inside the U file itself so it’s hard to open and edit .

————————-Here is my log file (Case run with one CPU) till 17 sec further in next comment ————————————————
/*—————————————————————————*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 2112 |
| \\ / A nd | Website: http://www.openfoam.com |
| \\/ M anipulation | |
\*—————————————————————————*/
Build : _80318542-20220301 OPENFOAM=2112 patch=220310 version=2112
Arch : “LSB;label=32;scalar=64”
Exec : simpleFoam
Date : May 23 2022
Time : 16:05:40
Host : shivam-VirtualBox
PID : 2756
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.0001
field U tolerance 0.0001
field “(k|omega|epsilon)” tolerance 0.0001

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

Time = 1

smoothSolver: Solving for Ux, Initial residual = 1, Final residual = 0.00989705, No Iterations 1
smoothSolver: Solving for Uy, Initial residual = 1, Final residual = 0.0663119, No Iterations 3
smoothSolver: Solving for Uz, Initial residual = 1, Final residual = 0.0656421, No Iterations 3
GAMG: Solving for p, Initial residual = 1, Final residual = 0.0304375, No Iterations 2
time step continuity errors : sum local = 0.144626, global = -0.114889, cumulative = -0.114889
smoothSolver: Solving for epsilon, Initial residual = 0.0553614, Final residual = 0.00537022, No Iterations 1
smoothSolver: Solving for k, Initial residual = 1, Final residual = 0.0176687, No Iterations 2
ExecutionTime = 13.38 s ClockTime = 14 s

functionObjects::Q Q writing field: Q
Time = 2

smoothSolver: Solving for Ux, Initial residual = 0.271847, Final residual = 0.00394546, No Iterations 2
smoothSolver: Solving for Uy, Initial residual = 0.27479, Final residual = 0.00614664, No Iterations 2
smoothSolver: Solving for Uz, Initial residual = 0.274415, Final residual = 0.00606402, No Iterations 2
GAMG: Solving for p, Initial residual = 0.412026, Final residual = 0.0121121, No Iterations 2
time step continuity errors : sum local = 0.113383, global = -0.0788744, cumulative = -0.193764
smoothSolver: Solving for epsilon, Initial residual = 0.0606698, Final residual = 0.00108669, No Iterations 2
bounding epsilon, min: -422167 max: 7.23235e+06 average: 24837.4
smoothSolver: Solving for k, Initial residual = 0.565944, Final residual = 0.0362071, No Iterations 1
bounding k, min: -26.112 max: 287.322 average: 5.69652
ExecutionTime = 22.71 s ClockTime = 24 s

functionObjects::Q Q writing field: Q
Time = 3

smoothSolver: Solving for Ux, Initial residual = 0.0377192, Final residual = 0.000855506, No Iterations 2
smoothSolver: Solving for Uy, Initial residual = 0.258764, Final residual = 0.0067388, No Iterations 2
smoothSolver: Solving for Uz, Initial residual = 0.337134, Final residual = 0.0088295, No Iterations 2
GAMG: Solving for p, Initial residual = 0.66364, Final residual = 0.0242389, No Iterations 2
time step continuity errors : sum local = 0.0939994, global = -0.0399866, cumulative = -0.23375
smoothSolver: Solving for epsilon, Initial residual = 0.021247, Final residual = 0.0008237, No Iterations 1
bounding epsilon, min: -1.02885e+07 max: 1.64674e+08 average: 329205
smoothSolver: Solving for k, Initial residual = 0.0656044, Final residual = 0.00114888, No Iterations 2
bounding k, min: -421.725 max: 17303.7 average: 15.7877
ExecutionTime = 31.43 s ClockTime = 32 s

functionObjects::Q Q writing field: Q
Time = 4

smoothSolver: Solving for Ux, Initial residual = 0.0456486, Final residual = 0.00162854, No Iterations 2
smoothSolver: Solving for Uy, Initial residual = 0.241119, Final residual = 0.0092247, No Iterations 2
smoothSolver: Solving for Uz, Initial residual = 0.21975, Final residual = 0.0103297, No Iterations 2
GAMG: Solving for p, Initial residual = 0.221276, Final residual = 0.00874707, No Iterations 2
time step continuity errors : sum local = 0.127948, global = -0.0211867, cumulative = -0.254937
smoothSolver: Solving for epsilon, Initial residual = 0.126194, Final residual = 0.0107563, No Iterations 1
bounding epsilon, min: -3.91161e+08 max: 2.8077e+10 average: 1.34829e+06
smoothSolver: Solving for k, Initial residual = 0.0600958, Final residual = 0.00147496, No Iterations 2
bounding k, min: -1068.46 max: 36875.2 average: 34.8867
ExecutionTime = 40.66 s ClockTime = 42 s

functionObjects::Q Q writing field: Q
Time = 5

smoothSolver: Solving for Ux, Initial residual = 0.0827382, Final residual = 0.00247526, No Iterations 2
smoothSolver: Solving for Uy, Initial residual = 0.216919, Final residual = 0.00624594, No Iterations 2
smoothSolver: Solving for Uz, Initial residual = 0.316906, Final residual = 0.00784881, No Iterations 2
GAMG: Solving for p, Initial residual = 0.030459, Final residual = 0.00141559, No Iterations 2
time step continuity errors : sum local = 0.154168, global = -0.0402135, cumulative = -0.295151
smoothSolver: Solving for epsilon, Initial residual = 0.115569, Final residual = 0.00786473, No Iterations 1
bounding epsilon, min: -6.38214e+10 max: 5.68063e+11 average: 8.76855e+06
smoothSolver: Solving for k, Initial residual = 0.336445, Final residual = 0.015374, No Iterations 1
bounding k, min: -32204.9 max: 398023 average: 65.2734
ExecutionTime = 51.47 s ClockTime = 52 s

functionObjects::Q Q writing field: Q
Time = 6

smoothSolver: Solving for Ux, Initial residual = 0.284195, Final residual = 0.01848, No Iterations 1
smoothSolver: Solving for Uy, Initial residual = 0.569072, Final residual = 0.0451165, No Iterations 1
smoothSolver: Solving for Uz, Initial residual = 0.495431, Final residual = 0.014658, No Iterations 2
GAMG: Solving for p, Initial residual = 0.0235749, Final residual = 0.00119702, No Iterations 2
time step continuity errors : sum local = 0.117156, global = -0.0655649, cumulative = -0.360715
smoothSolver: Solving for epsilon, Initial residual = 0.550587, Final residual = 0.0134038, No Iterations 1
bounding epsilon, min: -3.42406e+10 max: 3.35776e+12 average: 4.25111e+07
smoothSolver: Solving for k, Initial residual = 0.924521, Final residual = 0.0650929, No Iterations 2
bounding k, min: -6.79802e+06 max: 1.12786e+08 average: 716.561
ExecutionTime = 61.53 s ClockTime = 63 s

functionObjects::Q Q writing field: Q
Time = 7

smoothSolver: Solving for Ux, Initial residual = 0.397139, Final residual = 0.0392837, No Iterations 1
smoothSolver: Solving for Uy, Initial residual = 0.436024, Final residual = 0.0276398, No Iterations 1
smoothSolver: Solving for Uz, Initial residual = 0.577909, Final residual = 0.0366144, No Iterations 1
GAMG: Solving for p, Initial residual = 0.000121011, Final residual = 6.7821e-06, No Iterations 2
time step continuity errors : sum local = 0.128981, global = -0.0726883, cumulative = -0.433404
smoothSolver: Solving for epsilon, Initial residual = 0.254553, Final residual = 0.000661106, No Iterations 1
bounding epsilon, min: -1.75136e+11 max: 7.02792e+15 average: 3.73749e+10
smoothSolver: Solving for k, Initial residual = 0.50943, Final residual = 0.0229487, No Iterations 2
bounding k, min: -167947 max: 8.2005e+07 average: 2199.69
ExecutionTime = 73.9 s ClockTime = 77 s

functionObjects::Q Q writing field: Q
Time = 8

smoothSolver: Solving for Ux, Initial residual = 0.571546, Final residual = 0.0250485, No Iterations 2
smoothSolver: Solving for Uy, Initial residual = 0.599197, Final residual = 0.0576888, No Iterations 1
smoothSolver: Solving for Uz, Initial residual = 0.665602, Final residual = 0.0502761, No Iterations 1
GAMG: Solving for p, Initial residual = 3.01744e-06, Final residual = 2.08424e-07, No Iterations 2
time step continuity errors : sum local = 0.290697, global = -0.0532586, cumulative = -0.486662
smoothSolver: Solving for epsilon, Initial residual = 0.559445, Final residual = 0.0297444, No Iterations 1
bounding epsilon, min: -1.21069e+16 max: 4.60713e+16 average: 3.12325e+11
smoothSolver: Solving for k, Initial residual = 0.896683, Final residual = 0.013848, No Iterations 2
bounding k, min: -1.08803e+08 max: 3.33322e+09 average: 43563.9
ExecutionTime = 84.48 s ClockTime = 89 s

functionObjects::Q Q writing field: Q
Time = 9

smoothSolver: Solving for Ux, Initial residual = 0.493766, Final residual = 0.0394829, No Iterations 1
smoothSolver: Solving for Uy, Initial residual = 0.714466, Final residual = 0.0426581, No Iterations 1
smoothSolver: Solving for Uz, Initial residual = 0.790086, Final residual = 0.0247358, No Iterations 1
GAMG: Solving for p, Initial residual = 4.10875e-05, Final residual = 2.52148e-06, No Iterations 2
time step continuity errors : sum local = 0.269278, global = -0.0104, cumulative = -0.497062
smoothSolver: Solving for epsilon, Initial residual = 0.426421, Final residual = 0.0160667, No Iterations 1
bounding epsilon, min: -1.38272e+16 max: 5.78226e+17 average: 4.0044e+12
smoothSolver: Solving for k, Initial residual = 0.624362, Final residual = 0.0357542, No Iterations 2
bounding k, min: -2.55023e+09 max: 2.08613e+10 average: 500249
ExecutionTime = 94.15 s ClockTime = 99 s

functionObjects::Q Q writing field: Q
Time = 10

smoothSolver: Solving for Ux, Initial residual = 0.793844, Final residual = 0.0673995, No Iterations 1
smoothSolver: Solving for Uy, Initial residual = 0.786494, Final residual = 0.0580539, No Iterations 1
smoothSolver: Solving for Uz, Initial residual = 0.697902, Final residual = 0.0528468, No Iterations 1
GAMG: Solving for p, Initial residual = 2.06892e-06, Final residual = 1.12094e-07, No Iterations 1
time step continuity errors : sum local = 13.5857, global = 0.0958229, cumulative = -0.40124
smoothSolver: Solving for epsilon, Initial residual = 0.261632, Final residual = 0.00437305, No Iterations 1
bounding epsilon, min: -8.80803e+17 max: 5.7818e+19 average: 4.43787e+14
smoothSolver: Solving for k, Initial residual = 0.806056, Final residual = 0.0555003, No Iterations 1
bounding k, min: -2.79711e+11 max: 1.57161e+12 average: 9.6369e+06
ExecutionTime = 103.22 s ClockTime = 109 s

functionObjects::Q Q writing field: Q
Time = 11

smoothSolver: Solving for Ux, Initial residual = 0.23964, Final residual = 0.00895226, No Iterations 2
smoothSolver: Solving for Uy, Initial residual = 0.801511, Final residual = 0.0704146, No Iterations 1
smoothSolver: Solving for Uz, Initial residual = 0.292904, Final residual = 0.0201066, No Iterations 1
GAMG: Solving for p, Initial residual = 3.34711e-06, Final residual = 1.44611e-07, No Iterations 1
time step continuity errors : sum local = 244.216, global = -0.0443874, cumulative = -0.445627
smoothSolver: Solving for epsilon, Initial residual = 0.751309, Final residual = 0.027489, No Iterations 1
bounding epsilon, min: -3.33019e+19 max: 4.603e+21 average: 3.21218e+16
smoothSolver: Solving for k, Initial residual = 0.960702, Final residual = 0.0864602, No Iterations 2
bounding k, min: -3.98028e+12 max: 3.86245e+13 average: 8.23519e+08
ExecutionTime = 112.34 s ClockTime = 119 s

functionObjects::Q Q writing field: Q
Time = 12

smoothSolver: Solving for Ux, Initial residual = 0.494802, Final residual = 0.0381896, No Iterations 1
smoothSolver: Solving for Uy, Initial residual = 0.732431, Final residual = 0.0534336, No Iterations 1
smoothSolver: Solving for Uz, Initial residual = 0.365625, Final residual = 0.0123083, No Iterations 2
GAMG: Solving for p, Initial residual = 5.48267e-09, Final residual = 7.65741e-10, No Iterations 1
time step continuity errors : sum local = 445.903, global = 0.260896, cumulative = -0.184731
smoothSolver: Solving for epsilon, Initial residual = 0.245403, Final residual = 0.000626323, No Iterations 1
bounding epsilon, min: -1.68801e+21 max: 3.83857e+23 average: 7.66296e+18
smoothSolver: Solving for k, Initial residual = 0.350399, Final residual = 0.0202855, No Iterations 3
bounding k, min: -1.00178e+13 max: 4.9389e+15 average: 6.75376e+10
ExecutionTime = 122.31 s ClockTime = 129 s

functionObjects::Q Q writing field: Q
Time = 13

smoothSolver: Solving for Ux, Initial residual = 0.636521, Final residual = 0.0272121, No Iterations 1
smoothSolver: Solving for Uy, Initial residual = 0.460268, Final residual = 0.0327091, No Iterations 1
smoothSolver: Solving for Uz, Initial residual = 0.772591, Final residual = 0.0375421, No Iterations 1
GAMG: Solving for p, Initial residual = 2.00031e-08, Final residual = 9.84217e-10, No Iterations 1
time step continuity errors : sum local = 347941, global = -903.153, cumulative = -903.338
smoothSolver: Solving for epsilon, Initial residual = 0.000547182, Final residual = 9.44748e-23, No Iterations 1
bounding epsilon, min: -1.20432e+27 max: 1.3842e+28 average: 2.6464e+22
smoothSolver: Solving for k, Initial residual = 0.0402618, Final residual = 0.00291539, No Iterations 1
bounding k, min: -1.46515e+17 max: 7.00623e+18 average: 4.58365e+13
ExecutionTime = 131.2 s ClockTime = 139 s

functionObjects::Q Q writing field: Q
Time = 14

smoothSolver: Solving for Ux, Initial residual = 0.531102, Final residual = 0.0128858, No Iterations 2
smoothSolver: Solving for Uy, Initial residual = 0.347821, Final residual = 0.00822206, No Iterations 2
smoothSolver: Solving for Uz, Initial residual = 0.489354, Final residual = 0.027593, No Iterations 1
GAMG: Solving for p, Initial residual = 4.32136e-08, Final residual = 1.169e-09, No Iterations 1
time step continuity errors : sum local = 1.27291e+28, global = -9.8493e+24, cumulative = -9.8493e+24
smoothSolver: Solving for epsilon, Initial residual = 0.988829, Final residual = 0.0154013, No Iterations 1
bounding epsilon, min: -1.42513e+47 max: 1.05931e+48 average: 1.90237e+43
smoothSolver: Solving for k, Initial residual = 0.989608, Final residual = 0.0736308, No Iterations 1
bounding k, min: -8.52767e+38 max: 3.95141e+40 average: 6.28537e+35
ExecutionTime = 140.21 s ClockTime = 149 s

functionObjects::Q Q writing field: Q
Time = 15

smoothSolver: Solving for Ux, Initial residual = 0.0229298, Final residual = 0.00197612, No Iterations 1
smoothSolver: Solving for Uy, Initial residual = 0.00362096, Final residual = 0.000256112, No Iterations 1
smoothSolver: Solving for Uz, Initial residual = 0.00425882, Final residual = 0.000364833, No Iterations 1
GAMG: Solving for p, Initial residual = 0.985123, Final residual = 0.0127833, No Iterations 1
time step continuity errors : sum local = 3.75679e+38, global = 1.82726e+26, cumulative = 1.72876e+26
smoothSolver: Solving for epsilon, Initial residual = 0.11101, Final residual = 3.47106e-13, No Iterations 1
bounding epsilon, min: -2.11468e+73 max: 3.71333e+74 average: 3.40876e+69
smoothSolver: Solving for k, Initial residual = 1, Final residual = 0.0654468, No Iterations 3
bounding k, min: -2.30273e+64 max: 1.63014e+67 average: 9.57377e+61
ExecutionTime = 149.06 s ClockTime = 158 s

functionObjects::Q Q writing field: Q
Time = 16

smoothSolver: Solving for Ux, Initial residual = 1.18881e-11, Final residual = 1.18881e-11, No Iterations 0
smoothSolver: Solving for Uy, Initial residual = 4.06287e-12, Final residual = 4.06287e-12, No Iterations 0
smoothSolver: Solving for Uz, Initial residual = 1.57865e-12, Final residual = 1.57865e-12, No Iterations 0
GAMG: Solving for p, Initial residual = 1, Final residual = 0.0110236, No Iterations 2
time step continuity errors : sum local = 6.5558e+69, global = 6.01583e+69, cumulative = 6.01583e+69
smoothSolver: Solving for epsilon, Initial residual = 0.0107266, Final residual = 0.000115828, No Iterations 2
bounding epsilon, min: -2.88185e+105 max: 3.00966e+107 average: 1.99037e+103
smoothSolver: Solving for k, Initial residual = 1, Final residual = 0.00995004, No Iterations 2
bounding k, min: -9.18343e+100 max: 7.70188e+102 average: 9.19828e+97
ExecutionTime = 158.06 s ClockTime = 168 s

functionObjects::Q Q writing field: Q

• ### shivam.shivamchoubey

Member
May 23, 2022 at 2:14 pm

————————–terminal output —————————–

Time = 17

smoothSolver: Solving for Ux, Initial residual = 0.150307, Final residual = 0.0130307, No Iterations 2
smoothSolver: Solving for Uy, Initial residual = 0.00981227, Final residual = 0.000849715, No Iterations 2
smoothSolver: Solving for Uz, Initial residual = 0.0118886, Final residual = 0.00103051, No Iterations 2
#0 Foam::error::printStack(Foam::Ostream&) at ??:?
#1 Foam::sigFpe::sigHandler(int) at ??:?
#3 Foam::scalarProduct<double, double>::type Foam::sumProd<double>(Foam::UList<double> const&, Foam::UList<double> const&) at ??:?
#4 Foam::PCG::scalarSolve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const at ??:?
#5 Foam::GAMGSolver::solveCoarsestLevel(Foam::Field<double>&, Foam::Field<double> const&) const at ??:?
#6 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 at ??:?
#7 Foam::GAMGSolver::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const at ??:?
#8 Foam::fvMatrix<double>::solveSegregated(Foam::dictionary const&) at ??:?
#9 Foam::fvMatrix<double>::solveSegregatedOrCoupled(Foam::dictionary const&) at ??:?
#10 Foam::fvMesh::solve(Foam::fvMatrix<double>&, Foam::dictionary const&) const at ??:?
#11 ? in /usr/lib/openfoam/openfoam2112/platforms/linux64GccDPInt32Opt/bin/simpleFoam
#12 __libc_start_main in /lib/x86_64-linux-gnu/libc.so.6
#13 ? in /usr/lib/openfoam/openfoam2112/platforms/linux64GccDPInt32Opt/bin/simpleFoam

• ### BarÄ±ÅŸ Bicer

Moderator
May 23, 2022 at 2:41 pm

Hi shivam,

thank you very much for detail information. I got the error point. when I prepared case file I copied /0 folder instead of /0.orig. Therefore, you see the developed U velocity file with which the simulation starts and diverge.

Sandeep will change the case folder today.

Please download the new case folder with 0.orig file and copy it before run : cp -r 0.orig 0

then you can run the simulation.

If any problem happens let me know please.

thank you.

BarÄ±ÅŸ

• ### shivam.shivamchoubey

Member
May 23, 2022 at 2:45 pm

you are correct but I have already changed my velocity file ..
this is my velocity file …… .

………………………………… 0/U

/*——————————–*- C++ -*———————————-*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: v2112 |
| \\ / A nd | Website: http://www.openfoam.com |
| \\/ M anipulation | |
\*—————————————————————————*/
FoamFile
{
version 2.0;
format ascii;
class volVectorField;
object U;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions [0 1 -1 0 0 0 0];

internalField uniform (1 0 0);

boundaryField
{
inlet
{
type fixedValue;
value uniform (20 0 0);
}

farfield
{
type inletOutlet;
inletValue uniform (1 0 0);
value \$internalField;

}

golfBall
{
type noSlip;
}

}

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

• ### BarÄ±ÅŸ Bicer

Moderator
May 23, 2022 at 2:49 pm

So do you get still error?

Why did you give 1 m/s for internalFields??? It should be 0 at the beginning

• ### shivam.shivamchoubey

Member
May 23, 2022 at 2:57 pm

Yes I am getting same error .I start with 0 but just for checking i did 1 m/s but still i have same error .i will send you folder of setup.

• ### BarÄ±ÅŸ Bicer

Moderator
May 28, 2022 at 11:07 am

please refer to my new golf_ball post for soluiton.

Best.

• ### shivam.shivamchoubey

Member
May 30, 2022 at 4:10 pm

Hello Mr Bicer ,

Everything is working fine .I would like to know what is the reason behind errors

thanks and regards
Shivam

• ### shivam.shivamchoubey

Member
June 2, 2022 at 7:42 am

Hello Mr, Bicer

Everything is working fine now. I just want to know what was the problem so if in future the same problem will come I know how to solve it.

thanks and regards
Shivam

Page 1 of 2