[eRzBeNgEl]

Hi guys,I have a strange problem with running rhoCentralFoam, this is what I did so far:
-Copied forwardstep tutorial - simulations runs the tutorial
-changed mesh and inlet conditions (see attachment) - simulation got an error!


Error Message:

Create time
Create mesh for time = 0
Reading thermophysical properties


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


#0 Foam::error:rintStack(Foam::Ostream&) in "/opt/openfoam222/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#1 Foam::sigFpe::sigHandler(int) in "/opt/openfoam222/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#2 in "/lib/x86_64-linux-gnu/libc.so.6"
#3 Foam::heThermo<Foam:siThermo, Foam:ureMixture<Foam::constTransport<Foam::speci es::thermo<Foam::hConstThermo<Foam:erfectGas<Foa m::specie> >, Foam::sensibleInternalEnergy> > > >::init() in "/opt/openfoam222/platforms/linux64GccDPOpt/lib/libfluidThermophysicalModels.so"
#4 Foam::heThermo<Foam:siThermo, Foam:ureMixture<Foam::constTransport<Foam::speci es::thermo<Foam::hConstThermo<Foam:erfectGas<Foa m::specie> >, Foam::sensibleInternalEnergy> > > >::heThermo(Foam::fvMesh const&, Foam::word const&) in "/opt/openfoam222/platforms/linux64GccDPOpt/lib/libfluidThermophysicalModels.so"
#5 FoamsiThermo::addfvMeshConstructorToTable<Foam::he PsiThermo<Foam:siThermo, FoaureMixture<Foam::constTransport<Foam::species:: thermo<Foam::hConstThermo<Foam:erfectGas<Foam::s pecie> >, Foam::sensibleInternalEnergy> > > > >::New(Foam::fvMesh const&, Foam::word const&) in "/opt/openfoam222/platforms/linux64GccDPOpt/lib/libfluidThermophysicalModels.so"
#6 Foam::autoPtr<Foam:siThermo> Foam::basicThermo::New<FosiThermo>(Foam::fvMesh const&, Foam::word const&) in "/opt/openfoam222/platforms/linux64GccDPOpt/lib/libfluidThermophysicalModels.so"
#7 Foam:siThermo::New(Foam::fvMesh const&, Foam::word const&) in "/opt/openfoam222/platforms/linux64GccDPOpt/lib/libfluidThermophysicalModels.so"
#8
in "/home/iagaxtma/OpenFOAM/iagaxtma-2.2.2/platforms/linux64GccDPOpt/bin/rhoCentralFoamRK4"
#9 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#10
in "/home/iagaxtma/OpenFOAM/iagaxtma-2.2.2/platforms/linux64GccDPOpt/bin/rhoCentralFoamRK4"
Floating point exception (core dumped)
OTE]

This is my U:

dimensions [0 1 -1 0 0 0 0];


internalField uniform (1 0 0);


boundaryField
{
in
{
type fixedValue;
value uniform (1 0 0);
}


out
{
type inletOutlet;
inletValue uniform (1 0 0);
value uniform (1 0 0);
}


sym1
{
type empty;
}


sym2
{
type empty;
}


cyl
{
type fixedValue;
value uniform (0 0 0);
}

T:

dimensions [0 0 0 1 0 0 0];
internalField uniform 1;

boundaryField
{
in
{
type fixedValue;
value uniform 1;
}
out
{
type inletOutlet;
inletValue uniform 1;
value uniform 1;
}
sym1
{
type empty;
}
sym2
{
type empty;
}
cyl
{
type zeroGradient;
}


}

and p:

dimensions [1 -1 -2 0 0 0 0];


internalField uniform 0;


boundaryField
{
in
{
type fixedValue;
value uniform 0;
}


out
{
type zeroGradient;
}


sym1
{
type empty;
}


sym2
{
type empty;
}


cyl
{
type zeroGradient;
}

and boundary file:

FoamFile
{
version 2.0;
format ascii;
class polyBoundaryMesh;
location "constant/polyMesh";
object boundary;
}


5
(
cyl
{
type wall;
nFaces 356;
startFace 248132;
}
in
{
type patch;
nFaces 178;
startFace 248488;
}
out
{
type patch;
nFaces 178;
startFace 248666;
}
sym1
{
type empty;
nFaces 124244;
startFace 248844;
}
sym2
{
type empty;
nFaces 124244;
startFace 373088;
}
)

thermofile

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


// Note: these are the properties for a "normalised" inviscid gas
// for which the speed of sound is 1 m/s at a temperature of 1K
// and gamma = 7/5
mixture
{
specie
{
nMoles 1;
molWeight 11640.3;
}
thermodynamics
{
Cp 2.5;
Hf 0;
}
transport
{
mu 0.1;
Pr 1;
}
}

Where is the problem?

[eRzBeNgEl]

any idea or help?

[maHein]

Hello eRzBeNgEl,

you set the pressure to zero. This leads to problems in the thermo class.

You should set it to one similar to the mentioned tutorial case.

Regards,

maHein