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

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

    Barış Bicer replied 11 months, 3 weeks ago 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.

    ********************************************************************************************************************
    kadir@kadir-VirtualBox:~/OpenFOAM/kadir-v2212/run/pitzDailysteady$ checkMesh
    /*—————————————————————————*\
    | ========= | |
    | \\ / 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
    Host : kadir-VirtualBox
    PID : 17266
    I/O : uncollated
    Case : /home/kadir/OpenFOAM/kadir-v2212/run/pitzDailysteady
    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.
    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 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.

    Checking basic cellZone addressing…
    No cellZones found.

    Checking basic pointZone addressing…
    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

    kadir@kadir-VirtualBox:~/OpenFOAM/kadir-v2212/run/pitzDailysteady$

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

    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)
    simpleGrading (2.5 1 1)

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

  • 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
    Host : kadir-VirtualBox
    PID : 25123
    I/O : uncollated
    Case : /home/kadir/OpenFOAM/kadir-v2212/run/pitzDailysteady
    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

    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

    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
    Host : kadir-VirtualBox
    PID : 29179
    I/O : uncollated
    Case : /home/kadir/OpenFOAM/kadir-v2212/run/pitzDailysteady
    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

    Reading field p

    Reading field U

    Reading/calculating face flux field phi

    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.

Log in to reply.

error: Content is protected !!