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

Selecting incompressible transport model Newtonian

Selecting incompressible transport model Newtonian

Selecting turbulence model type laminar

Selecting laminar stress model Stokes

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

error: Content is protected !!