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

Home Forums OpenFOAM Forum Tutorial 4: Multiphase Flow

  • Tutorial 4: Multiphase Flow

    Posted by Sreeharsh on July 18, 2022 at 8:46 am

    Hi Dr.Baris,
    Need help in resolving error while running Tutorial 4: Multiphase flow. After 18 iterations, it is showing floating point error (I have attached the log file). The checkMesh quality showed fine. I also tired running the case single and parallel. But still it was showing the error. Also, why the max volume fraction of alpha.water was showing above 1 and min. is negative?

    // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
    Create time
    Create mesh for time = 0
    PIMPLE: Operating solver in PISO mode
    Using LTS
    Reading field p_rgh

    Reading field U

    Reading/calculating face flux field phi

    Reading transportProperties

    Selecting incompressible transport model Newtonian

    Selecting incompressible transport model Newtonian

    Selecting turbulence model type RAS

    Selecting RAS turbulence model kOmegaSST

    Selecting patchDistMethod meshWave

    RAS

    {

    RASModel kOmegaSST;

    turbulence on;

    printCoeffs on;

    alphaK1 0.85;

    alphaK2 1;

    alphaOmega1 0.5;

    alphaOmega2 0.856;

    gamma1 0.555556;

    gamma2 0.44;

    beta1 0.075;

    beta2 0.0828;

    betaStar 0.09;

    a1 0.31;

    b1 1;

    c1 10;

    F3 false;

    decayControl false;

    kInf 0;

    omegaInf 0;

    }

    Reading g

    Reading hRef

    Calculating field g.h

    No MRF models present

    No finite volume options present

    GAMGPCG: Solving for pcorr, Initial residual = 1, Final residual = 6.9595e-06, No Iterations 18

    time step continuity errors : sum local = 3.56923e-10, global = 4.85772e-13, cumulative = 4.85772e-13

    Starting time loop

    forces forces:

    Not including porosity effects

    Flow time scale min/max = 0.00017245, 1

    Smoothed flow time scale min/max = 0.00017245, 0.269039

    Time = 1

    PIMPLE: iteration 1

    smoothSolver: Solving for alpha.water, Initial residual = 0.00046962, Final residual = 7.73149e-10, No Iterations 1

    Phase-1 volume fraction = 0.810874 Min(alpha.water) = 0 Max(alpha.water) = 1

    MULES: Correcting alpha.water

    MULES: Correcting alpha.water

    Phase-1 volume fraction = 0.810868 Min(alpha.water) = -5.23392e-10 Max(alpha.water) = 1

    GAMG: Solving for p_rgh, Initial residual = 1, Final residual = 0.00966686, No Iterations 1

    time step continuity errors : sum local = 0.000521702, global = 6.839e-07, cumulative = 6.83901e-07

    GAMG: Solving for p_rgh, Initial residual = 0.0130115, Final residual = 6.20341e-08, No Iterations 14

    time step continuity errors : sum local = 1.63163e-08, global = 4.53136e-09, cumulative = 6.88432e-07

    smoothSolver: Solving for omega, Initial residual = 0.00441697, Final residual = 0.000118074, No Iterations 1

    smoothSolver: Solving for k, Initial residual = 1, Final residual = 0.000542951, No Iterations 1

    ExecutionTime = 7.74 s ClockTime = 8 s
    .
    .
    .
    Time = 17

    PIMPLE: iteration 1

    smoothSolver: Solving for alpha.water, Initial residual = 3.00196e-07, Final residual = 3.18514e-09, No Iterations 1

    Phase-1 volume fraction = 0.80847 Min(alpha.water) = -7.51046e-05 Max(alpha.water) = 1.18753

    Applying the previous iteration compression flux

    MULES: Correcting alpha.water

    MULES: Correcting alpha.water

    MULES: Correcting alpha.water

    Phase-1 volume fraction = 0.80847 Min(alpha.water) = -7.51046e-05 Max(alpha.water) = 1.18753

    GAMG: Solving for p_rgh, Initial residual = 0.99924, Final residual = 0.0081606, No Iterations 5

    time step continuity errors : sum local = 1.84912e+92, global = 2.76032e+90, cumulative = 2.76021e+90

    GAMG: Solving for p_rgh, Initial residual = 9.7503e-05, Final residual = 9.91133e-08, No Iterations 9

    time step continuity errors : sum local = 1.22538e+94, global = 1.43976e+93, cumulative = 1.44252e+93

    smoothSolver: Solving for omega, Initial residual = 0.264678, Final residual = 0.0258379, No Iterations 4

    bounding omega, min: -598725 max: 3.60428e+06 average: 11407

    smoothSolver: Solving for k, Initial residual = 0.210568, Final residual = 0.0205878, No Iterations 52

    bounding k, min: -0.254872 max: 1.1628 average: 0.000596272

    ExecutionTime = 135.84 s ClockTime = 140 s

    forces forces write:

    Sum of forces

    Total : (-1.9935e+200 -5.95023e+200 -1.73853e+201)

    Pressure : (-1.9935e+200 -5.95023e+200 -1.73853e+201)

    Viscous : (-1.51573e+104 -3.99209e+104 -2.78858e+105)

    Sum of moments

    Total : (-2.46138e+199 -4.35685e+201 1.48104e+201)

    Pressure : (-2.46138e+199 -4.35685e+201 1.48104e+201)

    Viscous : (6.2602e+102 -7.00594e+105 1.00472e+105)

    Flow time scale min/max = 4.8786e-104, 3.68663e-87

    Smoothed flow time scale min/max = 4.8786e-104, 5.71673e-101

    Damped flow time scale min/max = 4.8786e-104, 5.71673e-101

    Time = 18

    PIMPLE: iteration 1

    smoothSolver: Solving for alpha.water, Initial residual = 7.72657e-07, Final residual = 2.7603e-09, No Iterations 2

    Phase-1 volume fraction = 0.80847 Min(alpha.water) = -7.49532e-05 Max(alpha.water) = 1.18465

    Applying the previous iteration compression flux

    MULES: Correcting alpha.water

    MULES: Correcting alpha.water

    MULES: Correcting alpha.water

    Phase-1 volume fraction = 0.80847 Min(alpha.water) = -7.49532e-05 Max(alpha.water) = 1.18465

    [2] #0 Foam::error::printStack(Foam::Ostream&) at ??:?

    [2] #1 Foam::sigFpe::sigHandler(int) at ??:?

    [2] #2 ? in /lib/x86_64-linux-gnu/libc.so.6

    [2] #3 Foam::GAMGSolver::scale(Foam::Field&, Foam::Field&, Foam::lduMatrix const&, Foam::FieldField const&, Foam::UPtrList const&, Foam::Field const&, unsigned char) const at ??:?

    [2] #4 Foam::GAMGSolver::Vcycle(Foam::PtrList const&, Foam::Field&, Foam::Field const&, Foam::Field&, Foam::Field&, Foam::Field&, Foam::Field&, Foam::Field&, Foam::PtrList<Foam::Field >&, Foam::PtrList<Foam::Field >&, unsigned char) const at ??:?

    [2] #5 Foam::GAMGSolver::solve(Foam::Field&, Foam::Field const&, unsigned char) const at ??:?

    [2] #6 Foam::fvMatrix::solveSegregated(Foam::dictionary const&) at ??:?

    [2] #7 Foam::fvMatrix::solveSegregatedOrCoupled(Foam::dictionary const&) at ??:?

    [2] #8 Foam::fvMesh::solve(Foam::fvMatrix&, Foam::dictionary const&) const at ??:?

    [2] #9 ? in ~/OpenFOAM/OpenFOAM-v2112/platforms/linux64GccDPInt32Opt/bin/interFoam

    [2] #10 ? in /lib/x86_64-linux-gnu/libc.so.6

    [2] #11 __libc_start_main in /lib/x86_64-linux-gnu/libc.so.6

    [2] #12 ? in ~/OpenFOAM/OpenFOAM-v2112/platforms/linux64GccDPInt32Opt/bin/interFoam

    [harsh-VirtualBox:04835] *** Process received signal ***

    [harsh-VirtualBox:04835] Signal: Floating point exception (8)

    [harsh-VirtualBox:04835] Signal code: (-6)

    [harsh-VirtualBox:04835] Failing at address: 0x3e8000012e3

    [harsh-VirtualBox:04835] [ 0] /lib/x86_64-linux-gnu/libc.so.6(+0x42520)[0x7f513a334520]

    [harsh-VirtualBox:04835] [ 1] /lib/x86_64-linux-gnu/libc.so.6(pthread_kill+0x12c)[0x7f513a388a7c]

    [harsh-VirtualBox:04835] [ 2] /lib/x86_64-linux-gnu/libc.so.6(raise+0x16)[0x7f513a334476]

    [harsh-VirtualBox:04835] [ 3] /lib/x86_64-linux-gnu/libc.so.6(+0x42520)[0x7f513a334520]

    [harsh-VirtualBox:04835] [ 4] /home/harsh/OpenFOAM/OpenFOAM-v2112/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so(_ZNK4Foam10GAMGSolver5scaleERNS_5FieldIdEES3_RKNS_9lduMatrixERKNS_10FieldFieldIS1_dEERKNS_8UPtrListIKNS_17lduInterfaceFieldEEERKS2_h+0xfe)[0x7f513af67c8e]

    [harsh-VirtualBox:04835] [ 5] /home/harsh/OpenFOAM/OpenFOAM-v2112/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so(_ZNK4Foam10GAMGSolver6VcycleERKNS_7PtrListINS_9lduMatrix8smootherEEERNS_5FieldIdEERKS8_S9_S9_S9_S9_S9_RNS1_IS8_EESD_h+0x87d)[0x7f513af696ad]

    [harsh-VirtualBox:04835] [ 6] /home/harsh/OpenFOAM/OpenFOAM-v2112/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so(_ZNK4Foam10GAMGSolver5solveERNS_5FieldIdEERKS2_h+0x5b8)[0x7f513af6b9d8]

    [harsh-VirtualBox:04835] [ 7] /home/harsh/OpenFOAM/OpenFOAM-v2112/platforms/linux64GccDPInt32Opt/lib/libfiniteVolume.so(_ZN4Foam8fvMatrixIdE15solveSegregatedERKNS_10dictionaryE+0x5f2)[0x7f513e732c52]

    [harsh-VirtualBox:04835] [ 8] /home/harsh/OpenFOAM/OpenFOAM-v2112/platforms/linux64GccDPInt32Opt/lib/libfiniteVolume.so(_ZN4Foam8fvMatrixIdE24solveSegregatedOrCoupledERKNS_10dictionaryE+0x51a)[0x7f513db7a7ca]

    [harsh-VirtualBox:04835] [ 9] /home/harsh/OpenFOAM/OpenFOAM-v2112/platforms/linux64GccDPInt32Opt/lib/libfiniteVolume.so(_ZNK4Foam6fvMesh5solveERNS_8fvMatrixIdEERKNS_10dictionaryE+0x28)[0x7f513db20738]

    [harsh-VirtualBox:04835] [10] interFoam(+0x5299a)[0x563cb2b7899a]

    [harsh-VirtualBox:04835] [11] /lib/x86_64-linux-gnu/libc.so.6(+0x29d90)[0x7f513a31bd90]

    [harsh-VirtualBox:04835] [12] /lib/x86_64-linux-gnu/libc.so.6(__libc_start_main+0x80)[0x7f513a31be40]

    [harsh-VirtualBox:04835] [13] interFoam(+0x59065)[0x563cb2b7f065]

    [harsh-VirtualBox:04835] *** End of error message ***

    ————————————————————————–

    Primary job terminated normally, but 1 process returned

    a non-zero exit code. Per user-direction, the job has been aborted.

    Barış Bicer replied 1 year, 6 months ago 2 Members · 7 Replies
  • 7 Replies
  • Barış Bicer

    Moderator
    July 19, 2022 at 3:39 pm

    Hi Nair,

    Sınce your case is diverged alpha.water shows above 1 so it is wrong.

    Where can it come from?

    1. please paste here checkMesh results
    2. run case serial as laminar and paste also result here (diverged or not?)
    3. pls paste here system/fvSolution file

    and then lets see…

    Best

    Baris

  • Sreeharsh

    Member
    July 22, 2022 at 3:54 pm

    Hi Dr.Baris,

    As discussed, I am pasting the checkMesh results here. Also I ran the case as laminar, but there also it diverged. I am also pasting the warning which I received for setFields.

    1. Time = 0

    Mesh stats
    points: 919139
    faces: 2608507
    internal faces: 2538618
    cells: 845539
    faces per cell: 6.08739
    boundary patches: 7
    point zones: 0
    face zones: 0
    cell zones: 0

    Overall number of cells of each type:
    hexahedra: 788132
    prisms: 1511
    wedges: 0
    pyramids: 0
    tet wedges: 10
    tetrahedra: 0
    polyhedra: 55886
    Breakdown of polyhedra by number of faces:
    faces number of cells
    4 993
    5 827
    6 1285
    7 39811
    8 502
    9 12453
    10 14
    11 1

    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
    atmosphere 798 860 ok (non-closed singly connected)
    inlet 3192 3380 ok (non-closed singly connected)
    outlet 3192 3380 ok (non-closed singly connected)
    bottom 798 860 ok (non-closed singly connected)
    side 7056 7267 ok (non-closed singly connected)
    midPlane 27415 28561 ok (non-closed singly connected)
    hull 27438 28830 ok (non-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 (-26 -19 -16) (16 0 4)
    Mesh has 3 geometric (non-empty/wedge) directions (1 1 1)
    Mesh has 3 solution (non-empty) directions (1 1 1)
    Boundary openness (2.58769e-16 1.09768e-15 -1.12427e-16) OK.
    Max cell openness = 3.34969e-16 OK.
    Max aspect ratio = 75.9249 OK.
    Minimum face area = 7.7933e-06. Maximum face area = 1.00824. Face area magnitudes OK.
    Min volume = 1.07144e-07. Max volume = 0.933527. Total volume = 15958.8. Cell volumes OK.
    Mesh non-orthogonality Max: 69.9544 average: 6.61753
    Non-orthogonality check OK.
    Face pyramids OK.
    ***Max skewness = 5.23784, 1 highly skew faces detected which may impair the quality of the results
    <<Writing 1 skew faces to set skewFaces
    Coupled point location match (average 0) OK.

    Failed 1 mesh checks.
    End
    ————————————————————————————————————————————
    2. Laminar diverged case:
    PIMPLE: iteration 1

    smoothSolver: Solving for alpha.water, Initial residual = 9.17273e-08, Final residual = 5.06426e-10, No Iterations 1

    Phase-1 volume fraction = 0.812214 Min(alpha.water) = -0.000310776 Max(alpha.water) = 1.20042

    Applying the previous iteration compression flux

    MULES: Correcting alpha.water

    MULES: Correcting alpha.water

    MULES: Correcting alpha.water

    Phase-1 volume fraction = 0.812214 Min(alpha.water) = -0.000310776 Max(alpha.water) = 1.20042

    GAMG: Solving for p_rgh, Initial residual = 0.99959, Final residual = 0.00654842, No Iterations 4

    time step continuity errors : sum local = 9.52665e+91, global = 8.06293e+89, cumulative = 8.06127e+89

    GAMG: Solving for p_rgh, Initial residual = 0.000243327, Final residual = 7.65821e-08, No Iterations 11

    time step continuity errors : sum local = 2.25492e+93, global = 3.46299e+90, cumulative = 4.26912e+90

    ExecutionTime = 110.54 s ClockTime = 114 s

    forces forces write:

    Sum of forces

    Total : (-1.52811e+200 -3.71016e+200 -8.06535e+200)

    Pressure : (-1.52811e+200 -3.71016e+200 -8.06535e+200)

    Viscous : (-5.41544e+103 -1.42594e+104 -1.08155e+105)

    Sum of moments

    Total : (6.82809e+198 -2.00828e+201 9.16524e+200)

    Pressure : (6.82809e+198 -2.00828e+201 9.16524e+200)

    Viscous : (1.08138e+103 -2.71195e+105 3.58044e+104)

    Flow time scale min/max = 2.18341e-104, 4.37143e-87

    Smoothed flow time scale min/max = 2.18341e-104, 2.57979e-101

    Damped flow time scale min/max = 2.18341e-104, 2.57979e-101

    Time = 18

    PIMPLE: iteration 1

    smoothSolver: Solving for alpha.water, Initial residual = 1.64103e-07, Final residual = 1.10961e-09, No Iterations 1

    Phase-1 volume fraction = 0.812214 Min(alpha.water) = -0.000310351 Max(alpha.water) = 1.19992

    Applying the previous iteration compression flux

    MULES: Correcting alpha.water

    MULES: Correcting alpha.water

    MULES: Correcting alpha.water

    Phase-1 volume fraction = 0.812214 Min(alpha.water) = -0.000310351 Max(alpha.water) = 1.19992

    [2] #0 Foam::error::printStack(Foam::Ostream&) at ??:?

    [2] #1 Foam::sigFpe::sigHandler(int) at ??:?

    [2] #2 ? in /lib/x86_64-linux-gnu/libc.so.6

    [2] #3 Foam::GAMGSolver::scale(Foam::Field<double>&, Foam::Field<double>&, Foam::lduMatrix const&, Foam::FieldField<Foam::Field, double> const&, Foam::UPtrList<Foam::lduInterfaceField const> const&, Foam::Field<double> const&, unsigned char) const at ??:?

    [2] #4 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 ??:?

    [2] #5 Foam::GAMGSolver::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const at ??:?

    [2] #6 Foam::fvMatrix<double>::solveSegregated(Foam::dictionary const&) at ??:?

    [2] #7 Foam::fvMatrix<double>::solveSegregatedOrCoupled(Foam::dictionary const&) at ??:?

    [2] #8 Foam::fvMesh::solve(Foam::fvMatrix<double>&, Foam::dictionary const&) const at ??:?

    [2] #9 ? in ~/OpenFOAM/OpenFOAM-v2112/platforms/linux64GccDPInt32Opt/bin/interFoam

    [2] #10 ? in /lib/x86_64-linux-gnu/libc.so.6

    [2] #11 __libc_start_main in /lib/x86_64-linux-gnu/libc.so.6

    [2] #12 ? in ~/OpenFOAM/OpenFOAM-v2112/platforms/linux64GccDPInt32Opt/bin/interFoam

    [harsh-VirtualBox:06541] *** Process received signal ***

    [harsh-VirtualBox:06541] Signal: Floating point exception (8)

    [harsh-VirtualBox:06541] Signal code: (-6)

    [harsh-VirtualBox:06541] Failing at address: 0x3e80000198d

    [harsh-VirtualBox:06541] [ 0] /lib/x86_64-linux-gnu/libc.so.6(+0x42520)[0x7fe81bed0520]

    [harsh-VirtualBox:06541] [ 1] /lib/x86_64-linux-gnu/libc.so.6(pthread_kill+0x12c)[0x7fe81bf24a7c]

    [harsh-VirtualBox:06541] [ 2] /lib/x86_64-linux-gnu/libc.so.6(raise+0x16)[0x7fe81bed0476]

    [harsh-VirtualBox:06541] [ 3] /lib/x86_64-linux-gnu/libc.so.6(+0x42520)[0x7fe81bed0520]

    [harsh-VirtualBox:06541] [ 4] /home/harsh/OpenFOAM/OpenFOAM-v2112/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so(_ZNK4Foam10GAMGSolver5scaleERNS_5FieldIdEES3_RKNS_9lduMatrixERKNS_10FieldFieldIS1_dEERKNS_8UPtrListIKNS_17lduInterfaceFieldEEERKS2_h+0xea)[0x7fe81cb03c7a]

    [harsh-VirtualBox:06541] [ 5] /home/harsh/OpenFOAM/OpenFOAM-v2112/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so(_ZNK4Foam10GAMGSolver6VcycleERKNS_7PtrListINS_9lduMatrix8smootherEEERNS_5FieldIdEERKS8_S9_S9_S9_S9_S9_RNS1_IS8_EESD_h+0x87d)[0x7fe81cb056ad]

    [harsh-VirtualBox:06541] [ 6] /home/harsh/OpenFOAM/OpenFOAM-v2112/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so(_ZNK4Foam10GAMGSolver5solveERNS_5FieldIdEERKS2_h+0x5b8)[0x7fe81cb079d8]

    [harsh-VirtualBox:06541] [ 7] /home/harsh/OpenFOAM/OpenFOAM-v2112/platforms/linux64GccDPInt32Opt/lib/libfiniteVolume.so(_ZN4Foam8fvMatrixIdE15solveSegregatedERKNS_10dictionaryE+0x5f2)[0x7fe8202cec52]

    [harsh-VirtualBox:06541] [ 8] /home/harsh/OpenFOAM/OpenFOAM-v2112/platforms/linux64GccDPInt32Opt/lib/libfiniteVolume.so(_ZN4Foam8fvMatrixIdE24solveSegregatedOrCoupledERKNS_10dictionaryE+0x51a)[0x7fe81f7167ca]

    [harsh-VirtualBox:06541] [ 9] /home/harsh/OpenFOAM/OpenFOAM-v2112/platforms/linux64GccDPInt32Opt/lib/libfiniteVolume.so(_ZNK4Foam6fvMesh5solveERNS_8fvMatrixIdEERKNS_10dictionaryE+0x28)[0x7fe81f6bc738]

    [harsh-VirtualBox:06541] [10] interFoam(+0x5299a)[0x55b9824df99a]

    [harsh-VirtualBox:06541] [11] /lib/x86_64-linux-gnu/libc.so.6(+0x29d90)[0x7fe81beb7d90]

    [harsh-VirtualBox:06541] [12] /lib/x86_64-linux-gnu/libc.so.6(__libc_start_main+0x80)[0x7fe81beb7e40]

    [harsh-VirtualBox:06541] [13] interFoam(+0x59065)[0x55b9824e6065]

    [harsh-VirtualBox:06541] *** End of error message ***

    ————————————————————————–
    Primary job terminated normally, but 1 process returned
    a non-zero exit code. Per user-direction, the job has been aborted.
    ——————————————————————————————————————————–
    3. fvSolution file:
    solvers

    {

    “alpha.water.*”

    {

    nAlphaCorr 2;

    nAlphaSubCycles 1;

    cAlpha 1;

    icAlpha 0;

    MULESCorr yes;

    nLimiterIter 10;

    alphaApplyPrevCorr yes;

    solver smoothSolver;

    smoother symGaussSeidel;

    tolerance 1e-8;

    relTol 0;

    minIter 1;

    }

    “pcorr.*”

    {

    solver PCG;

    preconditioner

    {

    preconditioner GAMG;

    smoother GaussSeidel;

    tolerance 1e-5;

    relTol 0;

    };

    tolerance 1e-5;

    relTol 0;

    };

    p_rgh

    {

    solver GAMG;

    smoother DIC;

    tolerance 1e-7;

    relTol 0.01;

    };

    p_rghFinal

    {

    $p_rgh;

    relTol 0;

    }

    “(U|k|omega).*”

    {

    solver smoothSolver;

    smoother symGaussSeidel;

    nSweeps 1;

    tolerance 1e-7;

    relTol 0.1;

    minIter 1;

    };

    }

    PIMPLE

    {

    momentumPredictor no;

    nOuterCorrectors 1;

    nCorrectors 2;

    nNonOrthogonalCorrectors 0;

    maxCo 10;

    maxAlphaCo 5;

    rDeltaTSmoothingCoeff 0.05;

    rDeltaTDampingCoeff 0.5;

    nAlphaSpreadIter 0;

    nAlphaSweepIter 0;

    maxDeltaT 1;

    }

    relaxationFactors

    {

    equations

    {

    “.*” 1;

    }

    }

    cache

    {

    grad(U);

    }

    ————————————————————————————————————————–
    4. setFields warning:

    Create mesh for time = 0

    Reading setFieldsDict

    Setting field default values
    Setting internal values of volScalarField alpha.water

    Setting field region values
    Adding cells with centre within boxes 1((-999 -999 -999) (999 999 0.244))
    Setting internal values of volScalarField alpha.water
    Adding faces with centre within boxes 1((-999 -999 -999) (999 999 0.244))
    Setting patchField values of volScalarField alpha.water
    –> FOAM Warning :
    From bool setFaceFieldType(const Foam::word&, const Foam::fvMesh&, const labelList&, Foam::Istream&) [with Type = double; Foam::labelList = Foam::List<int>]
    in file setFields.C at line 277
    Ignoring internal face 0. Suppressing further warnings.
    On patch inlet set 1976 values
    On patch outlet set 1976 values
    On patch bottom set 798 values
    On patch side set 4368 values
    On patch midPlane set 16369 values
    On patch hull set 11268 values
    End

  • Barış Bicer

    Moderator
    July 22, 2022 at 4:13 pm

    Hi,

    Can you paste setFieldDict dictionary which is in system folder. Also, what do you see in alpha.water file located in /0 folder.

    Please copy here !

    • Sreeharsh

      Member
      July 22, 2022 at 4:27 pm

      Hi,

      setFields dictionary:
      defaultFieldValues

      (

      volScalarFieldValue alpha.water 0

      );

      regions

      (

      // Set cell values

      // (does zerogradient on boundaries)

      boxToCell

      {

      box (-999 -999 -999) (999 999 0.244);

      fieldValues

      (

      volScalarFieldValue alpha.water 1

      );

      }

      // Set patch values (using ==)

      boxToFace

      {

      box (-999 -999 -999) (999 999 0.244);

      fieldValues

      (

      volScalarFieldValue alpha.water 1

      );

      }

      );
      ———————————————————————————————-
      alpha.water in 0 folder (But I am unable to open the alpha.water file after I setfields):
      dimensions [0 0 0 0 0 0 0];

      internalField uniform 0;

      boundaryField

      {

      #includeEtc “caseDicts/setConstraintTypes”

      inlet

      {

      type fixedValue;

      value $internalField;

      }

      outlet

      {

      type variableHeightFlowRate;

      lowerBound 0;

      upperBound 1;

      value $internalField;

      }

      atmosphere

      {

      type inletOutlet;

      inletValue $internalField;

      value $internalField;

      }

      hull

      {

      type zeroGradient;

      }

      }
      ———————————————————————————

  • Barış Bicer

    Moderator
    July 23, 2022 at 8:07 am

    hı again,

    Ok I got the point. Because you run and manipulated alpha.water file already. delete the alpha.water file which you cant open and copy again from 0.org folder the original one. And then run setFields again. But before doing this you can make your mesh a bit finer and can try to run afterwards.

    Dont forget one thing that since it is multiphase flow and whenever you run the case alpha.water field is always initiated and not overwrite. therefore there is original 0.orh folder and thats why you have to copy alpha.water file from original folder and run setFields again if you wanna try a new condition.

    Best.

    Barış

    • Sreeharsh

      Member
      August 2, 2022 at 3:07 pm

      Hi Dr. Baris,

      As I was occupied with some official work last week, I was not able to reply.

      As you mentioned I have copied the alpha.water file from 0.orig folder to 0 folder whenever trying a new condition. But everytime the same error was coming. Even I tried making the mesh fine (which removed the skewness error which was there earlier), but this time also it got diverged and got crashed (see attached log file):

      —————————————————————————————————————–
      // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

      Create time

      Create mesh for time = 0

      PIMPLE: Operating solver in PISO mode

      Using LTS

      Reading field p_rgh

      Reading field U

      Reading/calculating face flux field phi

      Reading transportProperties

      Selecting incompressible transport model Newtonian

      Selecting incompressible transport model Newtonian

      Selecting turbulence model type laminar

      Selecting laminar stress model Stokes

      Reading g

      Reading hRef

      Calculating field g.h

      No MRF models present

      No finite volume options present

      GAMGPCG: Solving for pcorr, Initial residual = 1, Final residual = 6.9595e-06, No Iterations 18

      time step continuity errors : sum local = 3.56923e-10, global = 4.85772e-13, cumulative = 4.85772e-13

      Starting time loop

      forces forces:

      Not including porosity effects

      Flow time scale min/max = 0.00017245, 1

      Smoothed flow time scale min/max = 0.00017245, 0.269039

      Time = 1

      PIMPLE: iteration 1

      smoothSolver: Solving for alpha.water, Initial residual = 6.68516e-07, Final residual = 3.92466e-10, No Iterations 1

      Phase-1 volume fraction = 0.812212 Min(alpha.water) = 0 Max(alpha.water) = 1

      MULES: Correcting alpha.water

      MULES: Correcting alpha.water

      Phase-1 volume fraction = 0.812212 Min(alpha.water) = -5.23392e-10 Max(alpha.water) = 1

      GAMG: Solving for p_rgh, Initial residual = 1, Final residual = 0.00177646, No Iterations 1

      time step continuity errors : sum local = 0.000100447, global = 2.8722e-06, cumulative = 2.8722e-06

      GAMG: Solving for p_rgh, Initial residual = 0.00354171, Final residual = 8.81127e-08, No Iterations 11

      time step continuity errors : sum local = 9.5026e-08, global = 3.78462e-09, cumulative = 2.87599e-06

      ExecutionTime = 6.38 s ClockTime = 7 s

      forces forces write:

      Sum of forces

      Total : (-13077.5 24766.9 -62.7061)

      Pressure : (-13052.7 24740.8 -71.7739)

      Viscous : (-24.8188 26.1402 9.0678)

      Sum of moments

      Total : (1681.96 -38078.1 115996)

      Pressure : (1621.35 -36576.2 115732)

      Viscous : (60.6141 -1501.84 264.565)

      Flow time scale min/max = 0.000153987, 1

      Smoothed flow time scale min/max = 0.000153987, 0.232208

      Time = 2

      PIMPLE: iteration 1

      smoothSolver: Solving for alpha.water, Initial residual = 1.54195e-06, Final residual = 5.52944e-09, No Iterations 2

      Phase-1 volume fraction = 0.812212 Min(alpha.water) = -9.59104e-14 Max(alpha.water) = 1.00001

      Applying the previous iteration compression flux

      MULES: Correcting alpha.water

      MULES: Correcting alpha.water

      MULES: Correcting alpha.water

      Phase-1 volume fraction = 0.812212 Min(alpha.water) = -2.22188e-06 Max(alpha.water) = 1.00001

      GAMG: Solving for p_rgh, Initial residual = 0.0952834, Final residual = 0.000422383, No Iterations 4

      time step continuity errors : sum local = 0.00355688, global = 9.16125e-05, cumulative = 9.44885e-05

      GAMG: Solving for p_rgh, Initial residual = 0.116765, Final residual = 5.77483e-08, No Iterations 17

      time step continuity errors : sum local = 5.34217e-05, global = -3.00548e-05, cumulative = 6.44337e-05

      ExecutionTime = 11.46 s ClockTime = 12 s

      forces forces write:

      Sum of forces

      Total : (-8.47027e+08 -2.00935e+09 -3.27415e+09)

      Pressure : (-8.29023e+08 -1.98268e+09 -2.77119e+09)

      Viscous : (-1.80047e+07 -2.66618e+07 -5.02963e+08)

      Sum of moments

      Total : (3.00102e+08 -7.53815e+09 4.25943e+09)

      Pressure : (2.6415e+08 -6.26669e+09 4.18937e+09)

      Viscous : (3.59519e+07 -1.27146e+09 7.00646e+07)

      Flow time scale min/max = 3.16352e-09, 1

      Smoothed flow time scale min/max = 3.16352e-09, 5.52247e-06

      Damped flow time scale min/max = 3.16352e-09, 5.52247e-06

      .
      .
      .
      .
      Time = 17

      PIMPLE: iteration 1

      smoothSolver: Solving for alpha.water, Initial residual = 9.17273e-08, Final residual = 5.06426e-10, No Iterations 1

      Phase-1 volume fraction = 0.812214 Min(alpha.water) = -0.000310776 Max(alpha.water) = 1.20042

      Applying the previous iteration compression flux

      MULES: Correcting alpha.water

      MULES: Correcting alpha.water

      MULES: Correcting alpha.water

      Phase-1 volume fraction = 0.812214 Min(alpha.water) = -0.000310776 Max(alpha.water) = 1.20042

      GAMG: Solving for p_rgh, Initial residual = 0.99959, Final residual = 0.00654842, No Iterations 4

      time step continuity errors : sum local = 9.52665e+91, global = 8.06293e+89, cumulative = 8.06127e+89

      GAMG: Solving for p_rgh, Initial residual = 0.000243327, Final residual = 7.65821e-08, No Iterations 11

      time step continuity errors : sum local = 2.25492e+93, global = 3.46299e+90, cumulative = 4.26912e+90

      ExecutionTime = 110.6 s ClockTime = 116 s

      forces forces write:

      Sum of forces

      Total : (-1.52811e+200 -3.71016e+200 -8.06535e+200)

      Pressure : (-1.52811e+200 -3.71016e+200 -8.06535e+200)

      Viscous : (-5.41544e+103 -1.42594e+104 -1.08155e+105)

      Sum of moments

      Total : (6.82809e+198 -2.00828e+201 9.16524e+200)

      Pressure : (6.82809e+198 -2.00828e+201 9.16524e+200)

      Viscous : (1.08138e+103 -2.71195e+105 3.58044e+104)

      Flow time scale min/max = 2.18341e-104, 4.37143e-87

      Smoothed flow time scale min/max = 2.18341e-104, 2.57979e-101

      Damped flow time scale min/max = 2.18341e-104, 2.57979e-101

      Time = 18

      PIMPLE: iteration 1

      smoothSolver: Solving for alpha.water, Initial residual = 1.64103e-07, Final residual = 1.10961e-09, No Iterations 1

      Phase-1 volume fraction = 0.812214 Min(alpha.water) = -0.000310351 Max(alpha.water) = 1.19992

      Applying the previous iteration compression flux

      MULES: Correcting alpha.water

      MULES: Correcting alpha.water

      MULES: Correcting alpha.water

      Phase-1 volume fraction = 0.812214 Min(alpha.water) = -0.000310351 Max(alpha.water) = 1.19992
      ———————————————————————————————————————————-

  • Barış Bicer

    Moderator
    September 26, 2022 at 6:19 pm

    Hi again,

    Could you please decrease your time step to 0.1 and try it again.

    Thank you.

    Best regards.

    Baris

Log in to reply.

error: Content is protected !!