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

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.9595e06, No Iterations 18
time step continuity errors : sum local = 3.56923e10, global = 4.85772e13, cumulative = 4.85772e13
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.68516e07, Final residual = 3.92466e10, No Iterations 1
Phase1 volume fraction = 0.812212 Min(alpha.water) = 0 Max(alpha.water) = 1
MULES: Correcting alpha.water
MULES: Correcting alpha.water
Phase1 volume fraction = 0.812212 Min(alpha.water) = 5.23392e10 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.8722e06, cumulative = 2.8722e06
GAMG: Solving for p_rgh, Initial residual = 0.00354171, Final residual = 8.81127e08, No Iterations 11
time step continuity errors : sum local = 9.5026e08, global = 3.78462e09, cumulative = 2.87599e06
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.54195e06, Final residual = 5.52944e09, No Iterations 2
Phase1 volume fraction = 0.812212 Min(alpha.water) = 9.59104e14 Max(alpha.water) = 1.00001
Applying the previous iteration compression flux
MULES: Correcting alpha.water
MULES: Correcting alpha.water
MULES: Correcting alpha.water
Phase1 volume fraction = 0.812212 Min(alpha.water) = 2.22188e06 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.16125e05, cumulative = 9.44885e05
GAMG: Solving for p_rgh, Initial residual = 0.116765, Final residual = 5.77483e08, No Iterations 17
time step continuity errors : sum local = 5.34217e05, global = 3.00548e05, cumulative = 6.44337e05
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.16352e09, 1
Smoothed flow time scale min/max = 3.16352e09, 5.52247e06
Damped flow time scale min/max = 3.16352e09, 5.52247e06
.
.
.
.
Time = 17PIMPLE: iteration 1
smoothSolver: Solving for alpha.water, Initial residual = 9.17273e08, Final residual = 5.06426e10, No Iterations 1
Phase1 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
Phase1 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.65821e08, 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.18341e104, 4.37143e87
Smoothed flow time scale min/max = 2.18341e104, 2.57979e101
Damped flow time scale min/max = 2.18341e104, 2.57979e101
Time = 18
PIMPLE: iteration 1
smoothSolver: Solving for alpha.water, Initial residual = 1.64103e07, Final residual = 1.10961e09, No Iterations 1
Phase1 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
Phase1 volume fraction = 0.812214 Min(alpha.water) = 0.000310351 Max(alpha.water) = 1.19992
———————————————————————————————————————————