[mdhfiz]

This is my first time testing a solver by using a significantly different geometry than the one provided in tutorials. Previously, I have fixed this issue (different solver) by changing the boundary conditions (the patch of U) of the test case. However, as of this solver, I am stuck. I'm not sure if it's because I have 2 inlets (an air inlet and a water inlet) or because of my boundary conditions. This is the out after I run twoPhaseEulerFoam:

/*---------------------------------------------------------------------------*\
  =========                 |
  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   \\    /   O peration     | Website:  https://openfoam.org
    \\  /    A nd           | Version:  7
     \\/     M anipulation  |
\*---------------------------------------------------------------------------*/
Build  : 7-f1e51ce8d094
Exec   : twoPhaseEulerFoam
Date   : Nov 15 2019
Time   : 08:34:29
Host   : "debian"
PID    : 7848
I/O    : uncollated
Case   : /home/hafiz/OpenFOAM/hafiz-7/run/jetSpray
nProcs : 1
sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster (fileModificationSkew 10)
allowSystemOperations : Allowing user-supplied system call operations

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time

Create mesh for time = 0


PIMPLE: No convergence criteria found


PIMPLE: No corrector convergence criteria found
        Calculations will do 3 corrections


PIMPLE: Operating solver in transient mode with 3 outer correctors



Reading g

Reading hRef
Creating twoPhaseSystem

Selecting thermodynamics package 
{
    type            heRhoThermo;
    mixture         pureMixture;
    transport       const;
    thermo          hConst;
    equationOfState perfectGas;
    specie          specie;
    energy          sensibleInternalEnergy;
}

Calculating face flux field phi.air
Selecting diameterModel for phase air: isothermal
Selecting turbulence model type RAS
Selecting RAS turbulence model mixtureKEpsilon
RAS
{
    RASModel        mixtureKEpsilon;
    turbulence      on;
    printCoeffs     on;
    Cmu             0.09;
    C1              1.44;
    C2              1.92;
    C3              1.92;
    Cp              0.25;
    sigmak          1;
    sigmaEps        1.3;
}

Selecting thermodynamics package 
{
    type            heRhoThermo;
    mixture         pureMixture;
    transport       const;
    thermo          eConst;
    equationOfState perfectFluid;
    specie          specie;
    energy          sensibleInternalEnergy;
}

Calculating face flux field phi.water
Selecting diameterModel for phase water: constant
Selecting turbulence model type RAS
Selecting RAS turbulence model mixtureKEpsilon
RAS
{
    RASModel        mixtureKEpsilon;
    turbulence      on;
    printCoeffs     on;
    Cmu             0.09;
    C1              1.44;
    C2              1.92;
    C3              1.92;
    Cp              0.25;
    sigmak          1;
    sigmaEps        1.3;
}

Selecting default blending method: linear
Selecting aspectRatioModel for (air in water): constant
Selecting aspectRatioModel for (water in air): constant
Selecting dragModel for (air and water): segregated
Selecting swarmCorrection for (air and water): none
Selecting dragModel for (air in water): SchillerNaumann
Selecting swarmCorrection for (air in water): none
Selecting dragModel for (water in air): SchillerNaumann
Selecting swarmCorrection for (water in air): none
Selecting virtualMassModel for (air in water): constantCoefficient
Selecting virtualMassModel for (water in air): constantCoefficient
Selecting heatTransferModel for (air in water): RanzMarshall
Selecting heatTransferModel for (water in air): RanzMarshall
Calculating field g.h

Reading field p_rgh

Creating field dpdt

Creating field kinetic energy K

No MRF models present

No finite volume options present
Courant Number mean: 0.0002 max: 0.0526316
Max Ur Courant Number = 0.088355
Calculating field DDtU1 and DDtU2


Starting time loop

fieldAverage fieldAverage1:
    Starting averaging at time 0

Courant Number mean: 0.0002 max: 0.0526316
Max Ur Courant Number = 0.088355
Time = 1e-06

PIMPLE: Iteration 1
MULES: Solving for alpha.air
MULES: Solving for alpha.air
alpha.air volume fraction = 1  Min(alpha.air) = 1  Max(alpha.air) = 1
Constructing momentum equations
smoothSolver:  Solving for e.air, Initial residual = 1, Final residual = 1.71606e-09, No Iterations 1
smoothSolver:  Solving for e.water, Initial residual = 1, Final residual = 2.2156e-16, No Iterations 1
min T.air 300
min T.water 300
GAMG:  Solving for p_rgh, Initial residual = 0.996029, Final residual = 1.67269e-11, No Iterations 2
#0  Foam::error::printStack(Foam::Ostream&) at ??:?
#1  Foam::sigFpe::sigHandler(int) at ??:?
#2  ? in "/lib/x86_64-linux-gnu/libc.so.6"
#3  ? at tensorField.C:?
#4  ? in "/home/hafiz/OpenFOAM/OpenFOAM-7/platforms/linux64GccDPInt64Opt/bin/twoPhaseEulerFoam"
#5  ? in "/home/hafiz/OpenFOAM/OpenFOAM-7/platforms/linux64GccDPInt64Opt/bin/twoPhaseEulerFoam"
#6  ? in "/home/hafiz/OpenFOAM/OpenFOAM-7/platforms/linux64GccDPInt64Opt/bin/twoPhaseEulerFoam"
#7  __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#8  ? in "/home/hafiz/OpenFOAM/OpenFOAM-7/platforms/linux64GccDPInt64Opt/bin/twoPhaseEulerFoam"
Floating point exception

My fail debugging approach to this problem was first, I checked the mesh:

checkMesh

and this is output:

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time

Create polyMesh for time = 0

Time = 0

Mesh stats
    points:           25704
    internal points:  0
    faces:            50501
    internal faces:   24799
    cells:            12550
    faces per cell:   6
    boundary patches: 5
    point zones:      0
    face zones:       0
    cell zones:       0

Overall number of cells of each type:
    hexahedra:     12550
    prisms:        0
    wedges:        0
    pyramids:      0
    tet wedges:    0
    tetrahedra:    0
    polyhedra:     0

Checking topology...
    Boundary definition OK.
 ***Total number of faces on empty patches is not divisible by the number of cells in the mesh. Hence this mesh is not 1D or 2D.
    Cell to face addressing 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                  
    waterInlet          1        4        ok (non-closed singly connected)  
    airInlet            50       102      ok (non-closed singly connected)  
    outlet              50       102      ok (non-closed singly connected)  
    walls               25500    25704    ok (non-closed singly connected)  
    defaultFaces        101      204      ok (non-closed singly connected)  

Checking geometry...
    Overall domain bounding box (0 0 0) (0.25 0.05 0.0001)
    Mesh has 1 geometric (non-empty/wedge) directions (1 0 0)
    Mesh has 1 solution (non-empty) directions (1 0 0)
    All edges aligned with or perpendicular to non-empty directions.
    Boundary openness (7.92192e-20 1.45684e-18 5.48246e-15) OK.
    Max cell openness = 1.05879e-16 OK.
    Max aspect ratio = 1 OK.
    Minimum face area = 5e-08. Maximum face area = 1e-06.  Face area magnitudes OK.
    Min volume = 5e-11. Max volume = 1e-10.  Total volume = 1.25e-06.  Cell volumes OK.
    Mesh non-orthogonality Max: 0 average: 0
    Non-orthogonality check OK.
    Face pyramids OK.
    Max skewness = 2.74172e-13 OK.
    Coupled point location match (average 0) OK.

Mesh OK.

End

I'm sure there's nothing wrong here, so I looked at my boundary conditions, and I'll share my U of air and water here:

U.air

/*--------------------------------*- C++ -*----------------------------------*\
  =========                 |
  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   \\    /   O peration     | Website:  https://openfoam.org
    \\  /    A nd           | Version:  7
     \\/     M anipulation  |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      binary;
    class       volVectorField;
    object      U.air;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions      [0 1 -1 0 0 0 0];

internalField   uniform (0 0 0);

boundaryField
{
    waterInlet
    {
        type               fixedValue;
        value              uniform (0 0 0);
    }
    
    airInlet
    {
        type            fixedValue;
        value           uniform (100 0 0);
    }
    
    outlet
    {
        type               pressureInletOutletVelocity;
        phi                phi.air;
        value              $internalField;
    }
    walls
    {
        type               fixedValue;
        value              uniform (0 0 0);
    }
}

// ************************************************************************* //

U.water

/*--------------------------------*- C++ -*----------------------------------*\
  =========                 |
  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   \\    /   O peration     | Website:  https://openfoam.org
    \\  /    A nd           | Version:  7
     \\/     M anipulation  |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      binary;
    class       volVectorField;
    object      U.water;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions      [0 1 -1 0 0 0 0];

internalField   uniform (0 0 0);

boundaryField
{
    waterInlet
    {
        type               fixedValue;
        value              uniform (0 176.71 0);
    }
    
    airInlet
    {
        type            fixedValue;
        value           uniform (0 0 0);
    }
    outlet
    {
        type               pressureInletOutletVelocity;
        phi                phi.water;
        value              $internalField;;
    }
    walls
    {
        type               fixedValue;
        value              uniform (0 0 0);
    }
}

// ************************************************************************* //

Are there anything wrong with my setup? The meshing (attached) seems to be fine.

PS: I am trying to simulate a water spray which comes from a small nozzle at the top left of the geometry and airflow comes from the inlet to further atomize the already present droplets. I know that lagrangian approach would be better than euler based solver but I read that some people have used it.

PSS: I won't be replying for another 6-7 hours because it's already midnight where I am (just in case somebody is replies to this thread)