[Tscar]

Hello,
may I ask how can I fix following error? Unfortunately, I was not able to find right answer in already existing threads.

My case is a flow in vertical pipe. Lower inlet air has a higher temperature and by buyoncy the air goes up. I already modeled the case with buoyantBoussinesqPimpleFoam and now i would like to add heat transfer from air to the solid pipe (wich has uniform outer temperature).

I created multi region mesh in Salome and i used following commands to export the mesh to openfoam and to create solid and fluid regions

ideasUnvToFoam meshname.unv
splitMeshRegions -cellZones -overwrite

Create time

Create fluid mesh for region fluid for time = 0

Create solid mesh for region solid for time = 0

*** Reading fluid mesh thermophysical properties for region fluid

    Adding to thermoFluid

Selecting thermodynamics package 
{
    type            heRhoThermo;
    mixture         pureMixture;
    transport       const;
    thermo          hConst;
    equationOfState perfectGas;
    specie          specie;
    energy          sensibleEnthalpy;
}

#0  Foam::error::printStack(Foam::Ostream&) at ??:?
#1  Foam::sigFpe::sigHandler(int) at ??:?
#2  ? in "/lib/x86_64-linux-gnu/libc.so.6"
#3   Foam::heRhoThermo<Foam::rhoThermo,  Foam::pureMixture<Foam::constTransport<Foam::species::thermo<Foam::hConstThermo<Foam::perfectGas<Foam::specie>  >, Foam::sensibleEnthalpy> > > >::calculate() at ??:?
#4    Foam::rhoThermo::addfvMeshConstructorToTable<Foam::heRhoThermo<Foam::rhoThermo,   Foam::pureMixture<Foam::constTransport<Foam::species::thermo<Foam::hConstThermo<Foam::perfectGas<Foam::specie>  >, Foam::sensibleEnthalpy> > > > >::New(Foam::fvMesh  const&, Foam::word const&) at ??:?
#5   Foam::autoPtr<Foam::rhoThermo>  Foam::basicThermo::New<Foam::rhoThermo>(Foam::fvMesh const&,  Foam::word const&) at ??:?
#6  Foam::rhoThermo::New(Foam::fvMesh const&, Foam::word const&) at ??:?
#7  ? at ??:?
#8  __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#9  ? at ??:?

Problems appears to be in the thermophysicalProperties dictionary for fluid region.

Thank you.

[Tscar]

I forgot to add my termophysicalProperties file

FoamFile
{
    version     2.0;
    format      ascii;
    class       dictionary;
    location    "constant";
    object      thermophysicalProperties;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

thermoType
{
    type            heRhoThermo;    
    mixture         pureMixture;    
    transport       const;
    thermo          hConst;
    equationOfState perfectGas;
    specie          specie;
    energy          sensibleEnthalpy;
}

mixture
{
    specie
    {
        nMoles      1;
        molWeight   28.9;
    }
    thermodynamics
    {
        Cp          1004.5;
        Hf          274228.5;
    }
    transport
    {
        mu          1.71969e-05;    //13.3*1.293e-6
        Pr          0.715;
    }
}


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

[fabian_roesler]

Are you solving in parallel? If so you have to add values for temperature and pressure to all boundary fields, especially the processor fields. Otherwise you will have division by zero in the first iteration of the ideal gas law.

Cheers

Fabian

[manbo]

Hi!
I'm experiencing the same issue after almost seven years with version 8 of openfoam. Did you at last find what was wrong back then?