CFD Online Discussion Forums

CFD Online Discussion Forums (https://www.cfd-online.com/Forums/)
-   OpenFOAM Running, Solving & CFD (https://www.cfd-online.com/Forums/openfoam-solving/)
-   -   Error in reactingFoam (https://www.cfd-online.com/Forums/openfoam-solving/147016-error-reactingfoam.html)

okstatecheme January 11, 2015 22:24
Error in reactingFoam
 
1 Attachment(s)
I'm trying to run a species tracking simulation using reactingFoam with noComubustion reaction model. My thermophysicalProperties dictionary appears to be correct but I keep getting this error:

Code:
/*---------------------------------------------------------------------------*\
| =========                |                                                |
| \\      /  F ield        | OpenFOAM: The Open Source CFD Toolbox          |
|  \\    /  O peration    | Version:  2.3.1                                |
|  \\  /    A nd          | Web:      www.OpenFOAM.org                      |
|    \\/    M anipulation  |                                                |
\*---------------------------------------------------------------------------*/
Build  : 2.3.1-bcfaaa7b8660
Exec  : reactingFoam
Date  : Jan 11 2015
Time  : 21:17:45
Host  : "brett-ubuntu"
PID    : 6271
Case  : /home/brett/Desktop/stack2
nProcs : 1
sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster
allowSystemOperations : Allowing user-supplied system call operations

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time

Create mesh for time = 0


Reading g
Creating reaction model

Selecting combustion model noCombustion<psiThermoCombustion>


--> FOAM FATAL IO ERROR:
"ill defined primitiveEntry starting at keyword 'mu' on line 55 and ending at line 55"

file: /home/brett/Desktop/stack2/constant/thermophysicalProperties at line 55.

    From function primitiveEntry::readEntry(const dictionary&, Istream&)
    in file lnInclude/IOerror.C at line 132.

FOAM exiting

brett@brett-ubuntu:~/Desktop/stack2$
Any ideas on what is going on? Attached is the thermophysicalProperties dict.

Code:
/*--------------------------------*- C++ -*----------------------------------*\
| =========                |                                                |
| \\      /  F ield        | OpenFOAM: The Open Source CFD Toolbox          |
|  \\    /  O peration    | Version:  2.3.1                                |
|  \\  /    A nd          | Web:      www.OpenFOAM.org                      |
|    \\/    M anipulation  |                                                |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version    2.0;
    format      ascii;
    class      dictionary;
    location    "constant";
    object      thermophysicalProperties;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

thermoType
{
    type            hePsiThermo;
    mixture        multiComponentMixture;
    transport      const;
    thermo          hConst;
    energy          sensibleInternalEnergy;
    equationOfState perfectGas;
    specie          specie;
}

species
(
    N2
    O2
    CO2
    H2S
    CH4
    H2O
);

inertSpecie        ;

O2
{
    specie
    {
        nMoles          1;
        molWeight      31.9988;
    }
    thermodynamics
    {
        Cp            918;
    Hf        0;
    }
    transport
    {
        mu              2.02-05;
        Pr              0.7;
    }
}

H2O
{
    specie
    {
        nMoles          1;
        molWeight      18.0153;
    }
    thermodynamics
    {
        Cp        3985;
    Hf        0;
    }
    transport
    {
        mu              1.8-05;
        Pr              0.7;
    }
}

CH4
{
    specie
    {
        nMoles          1;
        molWeight      16.0428;
    }
    thermodynamics
    {
        Cp        1000;
    Hf        0;
    }
    transport
    {
        mu              1.8-05;
        Pr              0.7;
    }
}

CO2
{
    specie
    {
        nMoles          1;
        molWeight      44.01;
    }
    thermodynamics
    {
        Cp            3640;
    Hf        0;
    }
    transport
    {
        mu              1.48-05;
        Pr              0.7;
    }
}

N2
{
    specie
    {
        nMoles          1;
        molWeight      28.0134;
    }
    thermodynamics
    {
        Cp        1040;
    Hf        0;
    }
    transport
    {
        mu              1.75-05;
        Pr              0.7;
    }
}

H2S
{
    specie
    {
        nMoles          1;
        molWeight      34.0758;
    }
    thermodynamics
    {
        Cp        1176;
    Hf        0;
    }
    transport
    {
        mu              1.13-05;
        Pr              0.7;
    }
}


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

Linse January 12, 2015 06:47
Entries for mu
 
I guess the problem will be solved by simply adding an "e" to the entries for mu, e.g. "2.02e-05" instead of "2.02-05".

Probably you are even lucky you did not have mu-values above 5e-05, otherwise you would have received results of a completely non-physical nature. ;-)

Cheers,
Bernhard

shenzhou1987 January 13, 2015 08:22
Hi, Do you solve your problem? I'm simulating free jet flow of a nozzle which contains two species, gas and air. There is no reaction between these two species. So can you tell me how to turn off combustion model and reaction model using reactingFoam? Hope for your reply. Thank you.

Linse January 13, 2015 08:54
Dear Shenzhou,

In future, for keeping the threads more comprehensible/logical AND for receiving results quicker, please
1.) Run a search via the search function of the forum. Many questions already are asked at some other place, and sometimes even I am astonished about the detail of some questions already asked elsewhere.
2.) If your question is not yet solved in a thread, please open a new thread with a descriptive headline, such as "Switching off reactions in reactingFoam". Then you do not have to hijack another thread which might or might not be solved yet. Readers looking for a solution for the problem specified in the thread title will be grateful! ;-)

Concerning your question: The first post in http://www.cfd-online.com/Forums/ope...urned-off.html (found by searching for "reactingFoam") already in the entry post provides the question to your answer.

Cheers,
Bernhard


PS: Please do not feel offended by 1.) and 2.), yours just was the 5th post hijacking a thread I read today...

okstatecheme January 13, 2015 21:37
Thank you Bernhard! I can't believe I missed that.

Linse January 14, 2015 08:26
Well, in recent years I have grown pedantic about spelling and such things... ;-)
Cheers,
Bernhard


All times are GMT -4. The time now is 01:39.