MRF Issues when using SimpleReactingParcelFoam

NotSoFatRabbit · August 10, 2017, 17:16

Greetings one and all,

I have been trying to model coal combustion using SimpleReactingParcelFoam. Every time I run the solver I receive an error relating to MRF. To the best of my understanding there isn't a MRF model involved in this solver so I have no clue as to why an MRF related error seems to be causing issues. Below is the crashed solver:

Code:

/*---------------------------------------------------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  v1612+                                |
|   \\  /    A nd           | Web:      www.OpenFOAM.com                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
Build  : v1612+
Exec   : simpleReactingParcelFoam
Date   : Aug 10 2017
Time   : 15:55:52
Host   : "christian-XPS"
PID    : 8255
Case   : /home/christian/openFoamRUN/work/v1612+/burnerSimpleReactingParcelFoam1612
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.000000


SIMPLE: no convergence criteria found. Calculations will run for 15 steps.


Reading g
Creating combustion model

Selecting combustion model PaSR<rhoChemistryCombustion>
Selecting chemistry type 
{
    chemistrySolver ode;
    chemistryThermo rho;
}

Selecting thermodynamics package 
{
    type            heRhoThermo;
    mixture         reactingMixture;
    transport       sutherland;
    thermo          janaf;
    energy          sensibleEnthalpy;
    equationOfState perfectGas;
    specie          specie;
}

Selecting chemistryReader foamChemistryReader
    elements not defined in "/home/christian/openFoamRUN/work/v1612+/burnerSimpleReactingParcelFoam1612/constant/reactions"
chemistryModel: Number of species = 6 and reactions = 2
Selecting ODE solver seulex
    using integrated reaction rate
Creating component thermo properties:
    multi-component carrier - 6 species
    liquids - 1 components
    solids - 2 components

Reading field U

Reading/calculating face flux field phi

Creating turbulence model

Selecting turbulence model type RAS
Selecting RAS turbulence model kEpsilon
kEpsilonCoeffs
{
    Cmu             0.09;
    C1              1.44;
    C2              1.92;
    C3              -0.33;
    sigmak          1;
    sigmaEps        1.3;
}

Creating multi-variate interpolation scheme

No MRF models present

Selecting radiationModel P1
Selecting absorptionEmissionModel binaryAbsorptionEmission
Selecting absorptionEmissionModel constantAbsorptionEmission
Selecting absorptionEmissionModel cloudAbsorptionEmission
Selecting scatterModel cloudScatter
Selecting transmissivityModel none

Constructing reacting cloud
Constructing particle forces
    Selecting particle force sphereDrag
    Selecting particle force gravity
Constructing cloud functions
    Selecting cloud function patchPostProcessing1 of type patchPostProcessing
    Selecting cloud function particleTracks1 of type particleTracks
Constructing particle injection models
Creating injector: model1
Selecting injection model patchInjection
    Constructing 3-D injection
Selecting distribution model fixedValue
Selecting dispersion model stochasticDispersionRAS
Selecting patch interaction model standardWallInteraction
Selecting stochastic collision model none
Selecting surface film model none
Selecting U integration scheme Euler
Selecting heat transfer model RanzMarshall
Selecting T integration scheme analytical
Selecting composition model singleMixtureFraction
Selecting phase change model liquidEvaporation
Participating liquid species:
    H2O
Selecting devolatilisation model none
Selecting surface reaction model none
No finite volume options present


Starting time loop

Time = 0.000010



--> FOAM FATAL ERROR: 

    [U[0 1 -2 0 0 0 0] ] + [(rho*MRFZoneList:acceleration)[1 -2 -2 0 0 0 0] ]

    From function void Foam::checkMethod(const Foam::fvMatrix<Type>&, const Foam::DimensionedField<Type, Foam::volMesh>&, const char*) [with Type = Foam::Vector<double>]
    in file /opt/OpenFOAM/OpenFOAM-v1612+/src/finiteVolume/lnInclude/fvMatrix.C at line 1292.

FOAM aborting

#0  Foam::error::printStack(Foam::Ostream&) at ??:?
#1  Foam::error::abort() at ??:?
#2  void Foam::checkMethod<Foam::Vector<double> >(Foam::fvMatrix<Foam::Vector<double> > const&, Foam::DimensionedField<Foam::Vector<double>, Foam::volMesh> const&, char const*) at ??:?
#3  Foam::tmp<Foam::fvMatrix<Foam::Vector<double> > > Foam::operator+<Foam::Vector<double> >(Foam::tmp<Foam::fvMatrix<Foam::Vector<double> > > const&, Foam::tmp<Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh> > const&) at ??:?
#4  ? at ??:?
#5  __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#6  ? at ??:?
Aborted (core dumped)

If anyone has any idea as to why I am seeing this error or more insight as to what is causing this error please let me know

!! Thanks in advance

NotSoFatRabbit · August 10, 2017, 22:24

I have attached the rho, p, pReal and U files as well as an image of what the case should look like.

Code:

/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  2.3.0                                 |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       volScalarField;
    object      rho;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

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

internalField   uniform 1.23;

boundaryField
{
    CoolerInlet 
    {
        type           zeroGradient;
    }
    AxInlet
    {
        type           zeroGradient;
    }
    TxInlet
    {
        type           zeroGradient;
    }
    SxInlet
    {
        type           zeroGradient;
    }
    KilnFarFieldExit 
    {
        type           fixedValue;
        value          $internalField;
    }
    ".*Wall.*"
    {
        type           zeroGradient;
    }
}

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

Code:

/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  v1706                                 |
|   \\  /    A nd           | Web:      www.OpenFOAM.com                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       volVectorField;
    location    "0";
    object      U;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

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


internalField   nonuniform List<vector> 
71676
(
(7.868595544 -1.101043453 11.33870473)
(-4.300311962 -0.2047302783 36.8489932)
(-5.02353482 0.0603800697 36.48272214)
(-3.142003803 0.226653061 32.29823214)
(-2.466122667 0.1531160831 29.79640277)
(-2.461805154 0.06999445143 33.77437758)
(-1.868883172 0.04316328079 40.64620763)
(-1.107449094 0.02260406809 46.75560802)
(-0.5018793325 0.006878261595 50.98365586)
(-0.1118547005 0.001555653395 53.18263994)
(0.1424500393 0.006428377257 53.13565268)
(0.3174053273 0.0151999123 50.40918342)
(0.4287574214 0.02408269523 45.17495783)
(0.4744194907 0.02933808403 38.76933274)
(0.4642216332 0.02791861526 32.60334695)
(0.4077914053 0.02159288311 27.24203027)
(0.3809780182 0.0130547208 23.02424007)
(0.08085403086 -0.0002319090571 18.71387508)
(8.536119161 -1.248611293 9.492207445)
.....lots more reference points....
)
;
    }
}


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

Code:

/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  v1706                                 |
|   \\  /    A nd           | Web:      www.OpenFOAM.com                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       volScalarField;
    location    "0";
    object      p;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

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


internalField   nonuniform List<scalar> 
71676
(
100801.8473
100765.5019
101059.2122
101214.3135
101332.0604
101385.7618
101392.1996
101368.4076
101344.7738
101331.9319.........
)
;

boundaryField
{
    CoolerInlet
    {
        type            zeroGradient;
    }
    BurnerFaceWall
    {
        type            zeroGradient;
    }
    AxInlet
    {
        type            zeroGradient;
    }
    TxInlet
    {
        type            zeroGradient;
    }
    SxInlet
    {
        type            zeroGradient;
    }
    KilnWalls
    {
        type            zeroGradient;
    }
    KilnFarFieldExit
    {
        type            fixedValue;
        value           uniform 101325;
    }
}


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

Code:

/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  2.3.0                                 |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       volScalarField;
    object      p;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

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

internalField   uniform 101325;

boundaryField
{
    CoolerInlet 
    {
        type           zeroGradient;
    }
    AxInlet
    {
        type           zeroGradient;
    }
    TxInlet
    {
        type           zeroGradient;
    }
    SxInlet
    {
        type           zeroGradient;
    }
    KilnFarFieldExit 
    {
        type           fixedValue;
        value          $internalField;
    }
    ".*Wall.*"
    {
        type           zeroGradient;
    }
}

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

The attached image, shows a cross section of a coal burner with the inlets on the left side.

sebastien_F1 · June 19, 2019, 06:54

Hi,

It is an issue of units.

1/ You might have done a mistake in your units. Check your p_rgh units.

2/or have initialize your flow with potentialFoam. potentialFoam is an incompressible solver, and you cannot use directly this solution to initialize a compressible solution unless you delete the phi generated by potentialFoam and only keep the velocity field.

Best,
Sebastien

Similar Threads
Thread	Thread Starter	Forum	Replies	Last Post
MRF and topoSet problem- Rotating volume doesn't rotate	andreas0209@hotmail.com	OpenFOAM	1	April 4, 2021 14:35
Fan-assisted Natural Ventilation Simulation - Issues with results for an MRF case	Tellur	OpenFOAM Running, Solving & CFD	0	July 15, 2016 03:04
MRF solving issues	Carno	OpenFOAM Running, Solving & CFD	2	May 30, 2016 06:54
Possibly serious MRF implementation issue	Ali Blues	OpenFOAM Bugs	1	December 16, 2015 07:04
MRF setup	andreas0209@hotmail.com	OpenFOAM Pre-Processing	1	August 6, 2015 10:36

June 19, 2019, 06:54		#3
sebastien_F1 New Member sebastien vilfayeau Join Date: Feb 2012 Posts: 14 Rep Power: 14	Hi, It is an issue of units. 1/ You might have done a mistake in your units. Check your p_rgh units. 2/or have initialize your flow with potentialFoam. potentialFoam is an incompressible solver, and you cannot use directly this solution to initialize a compressible solution unless you delete the phi generated by potentialFoam and only keep the velocity field. Best, Sebastien uckmhnds likes this.