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

Home Forums OpenFOAM Forum Tutorial 4: Multiphase Flow Reply To: Tutorial 4: Multiphase Flow

  • 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

error: Content is protected !!