CFD Online Discussion Forums

CFD Online Discussion Forums (https://www.cfd-online.com/Forums/)
-   OpenFOAM (https://www.cfd-online.com/Forums/openfoam/)
-   -   Problems rhoCentralFoam (https://www.cfd-online.com/Forums/openfoam/125743-problems-rhocentralfoam.html)

eRzBeNgEl October 30, 2013 05:35

Problems rhoCentralFoam
 
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:
Quote:

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::printStack(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::psiThermo, Foam::pureMixture<Foam::constTransport<Foam::speci es::thermo<Foam::hConstThermo<Foam::perfectGas<Foa m::specie> >, Foam::sensibleInternalEnergy> > > >::init() in "/opt/openfoam222/platforms/linux64GccDPOpt/lib/libfluidThermophysicalModels.so"
#4 Foam::heThermo<Foam::psiThermo, Foam::pureMixture<Foam::constTransport<Foam::speci es::thermo<Foam::hConstThermo<Foam::perfectGas<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::perfectGas<Foam::s pecie> >, Foam::sensibleInternalEnergy> > > > >::New(Foam::fvMesh const&, Foam::word const&) in "/opt/openfoam222/platforms/linux64GccDPOpt/lib/libfluidThermophysicalModels.so"
#6 Foam::autoPtr<Foam::psiThermo> Foam::basicThermo::New<FosiThermo>(Foam::fvMesh const&, Foam::word const&) in "/opt/openfoam222/platforms/linux64GccDPOpt/lib/libfluidThermophysicalModels.so"
#7 Foam::psiThermo::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:
Quote:

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:
Quote:

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:
Quote:

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:
Quote:

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
Quote:

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 October 30, 2013 12:13

any idea or help?

maHein October 31, 2013 16:35

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


All times are GMT -4. The time now is 15:07.