[shanvach]

Hi all,

I have been trying to use twoPhaseEulerFoam with cyclic boundary conditions.However I am getting nan everywhere as soon as it solves pEqn.H. Don't seem to know what's going wrong. I have added meanVelocityForce to drive the pipe flow.

fvOptions

momentumSource
{
type meanVelocityForce;
active yes;

meanVelocityForceCoeffs
{
selectionMode all;

fields (U.air);
Ubar (0.00125 0 0);
}
meanVelocityForceCoeffs
{
selectionMode all;

fields (U.water);
Ubar (0.00125 0 0);
}
}

I am getting the following error

/*---------------------------------------------------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  dev                                   |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
|*---------------------------------------------------------------------------*|
|* OpenFOAM for Windows 18.10 (v1)                                           *|
|* Built by CFD Support, www.cfdsupport.com (based on Symscape).             *|
\*---------------------------------------------------------------------------*/
Build  : dev-59be3e7
Exec   : D:\OpenFOAM_18.10\18.10\cygwin64\opt\OpenFOAM\OpenFOAM-dev\platforms\cygwin64mingw-w64DPInt32Opt\bin\twoPhaseEulerFoam.exe
Date   : Apr 25 2019
Time   : 17:51:14
Host   : "CKLAB20"
PID    : 11976
I/O    : uncollated
Case   : D:/OpenFOAM_18.10/18.10/cygwin64/OpenFOAM/OpenFOAM-dev/tutorials/multiphase/twoPhaseEulerFoam/laminar/cylinder_2D_10d_4k_2
nProcs : 1
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
        Calclations will do 3 corrections


Reading g

Reading hRef
test
Creating twoPhaseSystem

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

--> FOAM Warning :
    From function Foam::Field<Type>::Field(const Foam::word&, const Foam::dictionary&, Foam::label) [with Type = double; Foam::label = int]
    in file /opt/CFDSupportFOAM4.0/beta/OpenFOAM-dev/src/OpenFOAM/lnInclude/Field.C at line 324
    Reading "D:/OpenFOAM_18.10/18.10/cygwin64/OpenFOAM/OpenFOAM-dev/tutorials/multiphase/twoPhaseEulerFoam/laminar/cylinder_2D_10d_4k_2/0/T.air.boundaryField.WALL" from line 40 to line 41
    expected keyword 'uniform' or 'nonuniform', assuming deprecated Field format from Foam version 2.0.
Calculating face flux field phi.air
Selecting diameterModel for phase air: constant
Selecting turbulence model type laminar
Selecting laminar stress model Stokes
Selecting thermodynamics package
{
    type            heRhoThermo;
    mixture         pureMixture;
    transport       const;
    thermo          hConst;
    equationOfState rhoConst;
    specie          specie;
    energy          sensibleInternalEnergy;
}

Calculating face flux field phi.water
Selecting diameterModel for phase water: constant
Selecting turbulence model type laminar
Selecting laminar stress model Stokes
Selecting default blending method: none
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
test
Calculating field g.h

Reading field p_rgh

Creating field dpdt

Creating field kinetic energy K

No MRF models present

Creating finite volume options from "constant/fvOptions"

Selecting finite volume options model type meanVelocityForce
    Source: momentumSource
    - selecting all cells
    - selected 3486 cell(s) with volume 1.251e-011
    Initial pressure gradient = 0

Courant Number mean: 0.00996 max: 0.00996
Calculating field DDtU1 and DDtU2


Starting time loop

Courant Number mean: 0.00996 max: 0.00996
Max Ur Courant Number = 0
deltaT = 0.0001
Time = 0.0001

PIMPLE: Iteration 1
MULES: Solving for alpha.air
MULES: Solving for alpha.air
alpha.air volume fraction = 0.8  Min(alpha.air) = 0.8  Max(alpha.air) = 0.8
Constructing momentum equations
Selecting patchDistMethod meshWave
Selecting patchDistMethod meshWave
--> FOAM Warning :
    From function const volVectorField& Foam::wallDist::n() const
    in file fvMesh/wallDist/wallDist/wallDist.C at line 174
    n requested but 'nRequired' not specified in the wallDist dictionary
    Recalculating y and n fields.
sc=0
dyvi2=6.49167
fsid=(-0 -0 -0)
dis=2.5e-005
fl=(-0 -0 -0)
ncap=(0 1 0)
csid=77.4597
Pressure gradient source: uncorrected Ubar = 0.00125, pressure gradient = -3.22063e-008
smoothSolver:  Solving for e.air, Initial residual = 0.513552, Final residual = 0.513552, No Iterations 0
smoothSolver:  Solving for e.water, Initial residual = 0.514248, Final residual = 0.514248, No Iterations 0
min T.air 300
min T.water 300
nonePBiCGStab:  Solving for p_rgh, Initial residual = 0.999028, Final residual = nan, No Iterations 10000
Pressure gradient source: uncorrected Ubar = nan, pressure gradient = nan
PIMPLE: Iteration 2
MULES: Solving for alpha.air
MULES: Solving for alpha.air
alpha.air volume fraction = nan  Min(alpha.air) = nan  Max(alpha.air) = nan
Constructing momentum equations
Selecting patchDistMethod meshWave
Selecting patchDistMethod meshWave
--> FOAM Warning :
    From function const volVectorField& Foam::wallDist::n() const
    in file fvMesh/wallDist/wallDist/wallDist.C at line 174
    n requested but 'nRequired' not specified in the wallDist dictionary
    Recalculating y and n fields.
sc=nan
dyvi2=nan
fsid=(nan nan nan)
dis=2.5e-005
fl=(nan nan nan)
ncap=(0 1 0)
csid=77.4597
Pressure gradient source: uncorrected Ubar = nan, pressure gradient = nan
smoothSolver:  Solving for e.air, Initial residual = nan, Final residual = nan, No Iterations 1000
smoothSolver:  Solving for e.water, Initial residual = nan, Final residual = nan, No Iterations 1000
min T.air nan
min T.water nan
nonePBiCGStab:  Solving for p_rgh, Initial residual = nan, Final residual = nan, No Iterations 10000


--> FOAM FATAL IO ERROR:
wrong token type - expected Scalar, found on line 0 the word 'nan'

file: D:/OpenFOAM_18.10/18.10/cygwin64/OpenFOAM/OpenFOAM-dev/tutorials/multiphase/twoPhaseEulerFoam/laminar/cylinder_2D_10d_4k_2/system/data.solverPerformance.p_rgh at line 0.

    From function Foam::Istream& Foam::operator>>(Foam::Istream&, Foam::doubleScalar&)
    in file lnInclude/Scalar.C at line 93.

FOAM exiting