[SCHOLZC]

I am trying to setup a case which models a triplet injector. I have defined each of the species (reacting and products) in separate thermophysicalProperties files, and input the appropriate corresponding parameters for each (I am not sure if this is the best way to do it since I am relatively new to OpenFOAM). I am trying to model a Methane Reaction (CH4(liquid) & O2(liquid) reactants, CO2 & H2O biproducts), but as soon as the compilation (when attempting to run) reaches the O2 properties (which uses a reacting Mixutre instead of a pureMixture for the biproducts), it goes to reference the "reactions" file. I have looked at the BubbleColumnEvaporatingReacting tutorial and files, which have provided some guidance, but it seems that my model (which is running reactingMultiphaseEulerFoam) is requiring a different setup, or something is ill defined. It is only recognizing O2 as a viable species to refer to, and when including other species in the list it gives:

--> FOAM FATAL ERROR:
CO2 not found in table. Valid entries: 1(O2)

From function HashTable ...

If I remove all other species, it says that the "sum of mass fractions is zero for species 1(O2)". I have defined all of the species in phaseProperties & thermo.compressibleGas. Are there files I am missing or not attending to here? How can I resolve this issue?

Thank you! :-)

[SCHOLZC]

I was able to resolve the species not being recognized. I am now have issue with the following error:

Sum of mass fractions is zero for species

Does anyone know where this is resolved?

Thanks!

[avila.vc]

Hi, I have never dealt with this specific solver. But have you tried to modify just the phases names and proprieties from tutorials/multiphase/reactingMultiphaseEulerFoam/laminar/mixerVessel2D ? This case seems to be pretty well settled for 4 phases.

[sim2599]

Hi Scholzc,


How did you resolve the issue of the species not being recognized?

[kk415]

Hi All,

I am also facing the same error of "Sum of mass fractions is zero for species 2(N2 H2O)". I am using sprayFoam solver in OpenFoam-v1912.

My thermophysicalProperties file looks like this.

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

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

species         (N2 H2O);

N2
{
    specie
    {
        molWeight   28;
		Y			0.5;
    }
    thermodynamics
    {
        Cp          1000;
        Hf          0;
    }
    transport
    {
        mu          1.5e-5;
        Pr          0.7;
    }
}

H2O
{
    specie
    {
        molWeight   28;
		Y			0.5;
    }
    thermodynamics
    {
        Cp          4187;
        Hf          0;
    }
    transport
    {
        mu          1e-6;
        Pr          0.7;
    }
}
inertSpecie     N2;

liquids
{
    H2O
    {
        defaultCoeffs   yes;
    }
}

solids
{}

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

Any help will be deeply appreciated. I need to solve this problem soon.

[kk415]

It is resolved by adding the species file in the 0 folder.