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.
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.

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] [ 0] /lib/x86_64-linux-gnu/libc.so.6(+0x42520)[0x7fe81bed0520]

[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] [ 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;

nAlphaSweepIter 0;

maxDeltaT 1;

}

relaxationFactors

{

equations

{

“.*” 1;

}

}

cache

{

}

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

Create mesh for time = 0

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 !!