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

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

error: Content is protected !!