April 25, 2019, 13:01
|
Divergence while using cyclic boundary conditions in twoPhaseEulerFoam
|
#1
|
Member
Join Date: Apr 2016
Posts: 30
Rep Power: 10
|
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
Code:
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
Code:
/*---------------------------------------------------------------------------*\
| ========= | |
| \\ / 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
|
|
|