[vinz]

Dear forumers,

I am currently trying to simulate the injection of different gases into air.
So far, I had no problem simulating the injection of gas1 at inlet1. No problem either when I try to inject gas2 at inlet2 of my problem.

I can also inject 50% of each gas at both inlets. So inlet1 is injecting 50% of gas1 and 50% of gas2. The same at inlet 2.

It is actually working fine for any proportion of gas1/gas2 EXCEPT 100% / 0% for one of the two gases. So I can't inject 100% gas1 at inlet1 and at the same time 100% of gas2 at inlet2.
For all these cases, I don't inject any air at the inlets (so fixedValue BC at 0.0 for specie air).

When I say that it does not work, it is that I get an fpe when the code is building the thermoPhysicalModel.
Even more strange, the behaviour is changing depending on the order of declaration of species in the thermophysicalDict file!

Here are some informations on the settings:
- thermophysicalProperties:

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

species
(
    air
    gas1
    gas2
);
 
inertSpecie air;

air
{
    specie
    {
        nMoles      1;
        molWeight   28.97;
    }
    thermodynamics
    {
        Cp          1005;
        Hf          0;
    }
    transport
    {
        mu          1.78e-05;
        Pr          0.7;
    }
}

gas1
{
    specie
    {
        nMoles      1;
        molWeight   29.56;
    }
    thermodynamics
    {
        Cp          2445;
        Hf          0;
    }
    transport
    {
        mu          9e-05;
        Pr          0.7;
    }
}

gas2
{
    specie
    {
        nMoles      1;
        molWeight   16.55;
    }
    thermodynamics
    {
        Cp          2894.53;
        Hf          0;
    }
    transport
    {
        mu          9e-05;
        Pr          0.7;
    }
}

- gas1

dimensions      [0 0 0 0 0 0 0];

internalField   uniform 0.0;

boundaryField
{
    "inlet1"
    {
        type            fixedValue;
        value           uniform 1.0;
    }
    "inlet2"
    {
        type            fixedValue;
        value           uniform 0.0;
    }
....
}

- gas2

dimensions      [0 0 0 0 0 0 0];

internalField   uniform 0.0;

boundaryField
{
    "inlet1"
    {
        type            fixedValue;
        value           uniform 0.0;
    }
    "inlet2"
    {
        type            fixedValue;
        value           uniform 1.0;
    }
....
}

- air

dimensions      [0 0 0 0 0 0 0];

internalField   uniform 1.0;

boundaryField
{
    "inlet1"
    {
        type            fixedValue;
        value           uniform 0.0;
    }

    "inlet2"
    {
        type            fixedValue;
        value           uniform 0.0;
    }
....
}

- Crash

Reading thermophysical properties

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

#0  Foam::error::printStack(Foam::Ostream&) in "/CFDBB/vincent/OpenFOAM//OpenFOAM-2.3.0/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#1  Foam::sigFpe::sigHandler(int) in "/CFDBB/vincent/OpenFOAM//OpenFOAM-2.3.0/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#2   in "/lib/x86_64-linux-gnu/libc.so.6"
#3  Foam::multiComponentMixture<Foam::constTransport<Foam::species::thermo<Foam::hConstThermo<Foam::perfectGas<Foam::specie> >, Foam::sensibleInternalEnergy> > >::patchFaceMixture(int, int) const in "/CFDBB/vincent/OpenFOAM//OpenFOAM-2.3.0/platforms/linux64GccDPOpt/lib/libreactionThermophysicalModels.so"
#4  Foam::heThermo<Foam::psiReactionThermo, Foam::SpecieMixture<Foam::multiComponentMixture<Foam::constTransport<Foam::species::thermo<Foam::hConstThermo<Foam::perfectGas<Foam::specie> >, Foam::sensibleInternalEnergy> > > > >::he(Foam::Field<double> const&, Foam::Field<double> const&, int) const in "/CFDBB/vincent/OpenFOAM//OpenFOAM-2.3.0/platforms/linux64GccDPOpt/lib/libreactionThermophysicalModels.so"
#5  Foam::heThermo<Foam::psiReactionThermo, Foam::SpecieMixture<Foam::multiComponentMixture<Foam::constTransport<Foam::species::thermo<Foam::hConstThermo<Foam::perfectGas<Foam::specie> >, Foam::sensibleInternalEnergy> > > > >::init() in "/CFDBB/vincent/OpenFOAM//OpenFOAM-2.3.0/platforms/linux64GccDPOpt/lib/libreactionThermophysicalModels.so"
#6  Foam::heThermo<Foam::psiReactionThermo, Foam::SpecieMixture<Foam::multiComponentMixture<Foam::constTransport<Foam::species::thermo<Foam::hConstThermo<Foam::perfectGas<Foam::specie> >, Foam::sensibleInternalEnergy> > > > >::heThermo(Foam::fvMesh const&, Foam::word const&) in "/CFDBB/vincent/OpenFOAM//OpenFOAM-2.3.0/platforms/linux64GccDPOpt/lib/libreactionThermophysicalModels.so"
#7  Foam::psiReactionThermo::addfvMeshConstructorToTable<Foam::hePsiThermo<Foam::psiReactionThermo, Foam::SpecieMixture<Foam::multiComponentMixture<Foam::constTransport<Foam::species::thermo<Foam::hConstThermo<Foam::perfectGas<Foam::specie> >, Foam::sensibleInternalEnergy> > > > > >::New(Foam::fvMesh const&, Foam::word const&) in "/CFDBB/vincent/OpenFOAM//OpenFOAM-2.3.0/platforms/linux64GccDPOpt/lib/libreactionThermophysicalModels.so"
#8  Foam::autoPtr<Foam::psiReactionThermo> Foam::basicThermo::New<Foam::psiReactionThermo>(Foam::fvMesh const&, Foam::word const&) in "/CFDBB/vincent/OpenFOAM//OpenFOAM-2.3.0/platforms/linux64GccDPOpt/lib/libreactionThermophysicalModels.so"
#9  Foam::psiReactionThermo::New(Foam::fvMesh const&, Foam::word const&) in "/CFDBB/vincent/OpenFOAM//OpenFOAM-2.3.0/platforms/linux64GccDPOpt/lib/libreactionThermophysicalModels.so"

I guess that the problem is really coming from some 0 value at one of the boundary conditions but I don't see any reason why.
Does someone has an idea on what's going on?

[shenzhou1987]

Hi, I have meet the similar problem. Can U email me your case, so that I can learn how to deal with boundary condition and initial conditon. In addition, my case includes 4 cellZones, and I want to patch different values for different cellZone as initial condition. Can U tell me how to do this?
Thank U very much!