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

MRF Issues when using SimpleReactingParcelFoam

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

Like Tree1Likes
  • 1 Post By sebastien_F1

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   August 10, 2017, 16:16
Unhappy MRF Issues when using SimpleReactingParcelFoam
  #1
New Member
 
Christian
Join Date: Oct 2016
Posts: 2
Rep Power: 0
NotSoFatRabbit is on a distinguished road
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 is offline   Reply With Quote

Old   August 10, 2017, 21:24
Default
  #2
New Member
 
Christian
Join Date: Oct 2016
Posts: 2
Rep Power: 0
NotSoFatRabbit is on a distinguished road
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.
Attached Images
File Type: png ColdKiln.png (26.9 KB, 22 views)
NotSoFatRabbit is offline   Reply With Quote

Old   June 19, 2019, 05:54
Default
  #3
New Member
 
sebastien vilfayeau
Join Date: Feb 2012
Posts: 14
Rep Power: 10
sebastien_F1 is on a distinguished road
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.
sebastien_F1 is offline   Reply With Quote

Reply

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
Fan-assisted Natural Ventilation Simulation - Issues with results for an MRF case Tellur OpenFOAM Running, Solving & CFD 0 July 15, 2016 02:04
MRF solving issues Carno OpenFOAM Running, Solving & CFD 2 May 30, 2016 05:54
Possibly serious MRF implementation issue Ali Blues OpenFOAM Bugs 1 December 16, 2015 06:04
MRF setup andreas0209@hotmail.com OpenFOAM Pre-Processing 1 August 6, 2015 09:36
MRF and topoSet problem- Rotating volume doesn't rotate andreas0209@hotmail.com OpenFOAM 0 March 5, 2015 11:29


All times are GMT -4. The time now is 06:11.