CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Running, Solving & CFD

Divergence while using cyclic boundary conditions in twoPhaseEulerFoam

Register Blogs Members List Search Today's Posts Mark Forums Read

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   April 25, 2019, 13:01
Default Divergence while using cyclic boundary conditions in twoPhaseEulerFoam
  #1
Member
 
Join Date: Apr 2016
Posts: 30
Rep Power: 7
shanvach is on a distinguished road
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
shanvach is offline   Reply With Quote

Reply

Tags
cyclic bc, diverged, peqn.h, twophaseeulerfoam, ueqn.h

Thread Tools Search this Thread
Search this Thread:

Advanced Search
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
My radial inflow turbine Abo Anas CFX 27 May 11, 2018 01:44
Wrong multiphase flow at rotating interface Sanyo CFX 14 February 7, 2017 17:19
Boundary Conditions on Divergence in Marker and Cell (MAC) FluidFox Main CFD Forum 9 August 2, 2016 12:59
Wrong flow in ratating domain problem Sanyo CFX 17 August 15, 2015 06:20
Error finding variable "THERMX" sunilpatil CFX 8 April 26, 2013 07:00


All times are GMT -4. The time now is 21:43.