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 shivam.shivamchoubey 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] Failing at address: 0x3e800000eaf
    [shivam-VirtualBox:03759] [ 0] /lib/x86_64-linux-gnu/libpthread.so.0(+0x12980)[0x7fba65552980]
    [shivam-VirtualBox:03759] [ 1] /lib/x86_64-linux-gnu/libpthread.so.0(raise+0xc7)[0x7fba65552817]
    [shivam-VirtualBox:03759] [ 2] /lib/x86_64-linux-gnu/libpthread.so.0(+0x12980)[0x7fba65552980]
    [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] Failing at address: 0x3e800000eb0
    [shivam-VirtualBox:03760] [ 0] /lib/x86_64-linux-gnu/libpthread.so.0(+0x12980)[0x7f9ef6eb2980]
    [shivam-VirtualBox:03760] [ 1] /lib/x86_64-linux-gnu/libpthread.so.0(raise+0xc7)[0x7f9ef6eb2817]
    [shivam-VirtualBox:03760] [ 2] /lib/x86_64-linux-gnu/libpthread.so.0(+0x12980)[0x7f9ef6eb2980]
    [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 .

    Barış Bicer 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—————————
    shivam@shivam-VirtualBox:~/Downloads/Lecture_24_GolfBall_case$ checkMesh
    /*—————————————————————————*\
    | ========= | |
    | \\ / 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
    Case : /home/shivam/Downloads/Lecture_24_GolfBall_case
    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.
    Cell to face addressing 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.

    Checking basic cellZone addressing…
    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
    Case : /home/shivam/Downloads/Lecture_24_GolfBall_case
    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

    Reading field p

    Reading field U

    Reading/calculating face flux field phi

    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 ??:?
    #2 ? in /lib/x86_64-linux-gnu/libpthread.so.0
    #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
    shivam@shivam-VirtualBox:~/Downloads/Lecture_24_GolfBall_case$

  • 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

Log in to reply.

error: Content is protected !!