[vishalgarg474]

Create time

Create mesh for time = 0


PIMPLE: Operating solver in PISO mode

Reading field p_rgh

Calculating face flux field phi.waterdia1
Selecting diameterModel for phase waterdia1: constant
Calculating face flux field phi.waterdia2
Selecting diameterModel for phase waterdia2: constant
Calculating face flux field phi.air
Selecting diameterModel for phase air: constant
Selecting dragModel for phase waterdia2: blended
Selecting dragModel for phase waterdia2: SchillerNaumann
Selecting dragModel for phase air: SchillerNaumann
Selecting dragModel for phase waterdia1: blended
Selecting dragModel for phase waterdia1: SchillerNaumann
Selecting dragModel for phase waterdia2: SchillerNaumann
Selecting dragModel for phase waterdia1: blended
Selecting dragModel for phase waterdia1: SchillerNaumann
Selecting dragModel for phase air: SchillerNaumann
Selecting turbulence model type laminar

Reading g

Reading hRef
Calculating field g.h

No MRF models present

DICPCG:  Solving for pcorr, Initial residual = 1, Final residual = 4.39209e-06, No Iterations 27
time step continuity errors : sum local = 6.23677e-13, global = 6.57419e-14, cumulative = 6.57419e-14
Courant Number mean: 3.72941e-06 max: 0.000135479

Starting time loop

Courant Number mean: 3.72941e-06 max: 0.000135479
deltaT = 1.2e-08
Time = 1.2e-08

PIMPLE: iteration 1
MULES: Solving for alpha.waterdia1
waterdia1 volume fraction, min, max = 2.982e-08 0 0.3
MULES: Solving for alpha.waterdia2
waterdia2 volume fraction, min, max = 1.988e-08 0 0.2
MULES: Solving for alpha.air
air volume fraction, min, max = 1 0 1
Phase-sum volume fraction, min, max = 1 0 1
MULES: Solving for alpha.waterdia1
waterdia1 volume fraction, min, max = 5.964e-08 0 0.3
MULES: Solving for alpha.waterdia2
waterdia2 volume fraction, min, max = 3.976e-08 0 0.2
MULES: Solving for alpha.air
air volume fraction, min, max = 1 0 1
Phase-sum volume fraction, min, max = 1 0 1

Hey Guys, I am facing this error in multiphaseEulerFoam.

I am trying to analyze a water spray using this solver with air and two other phases for water with different droplet diameter. The source code of multiphaseEulerFoam is not completely executed and stops before zonePhaseVolumes.H (or Maybe it stops within the zonePhaseVolumes.H file).


As I understood, the fluid.solve() command is executed. Because the print statement in multiphasesystem.c solve and then solveAlphas is printed.



I am attaching my case folder as a zip file. Could you please help me understand the error or tell me why is the code stopping before the entire source code of multiphaseEulerFoam is executed?



Thank you in advance.

[vishalgarg474]

I found the problem.


The sum of volume fractions was not 1 at the boundary.

[Eslam Reda]

Boss