CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM

Problems rhoCentralFoam

Register Blogs Members List Search Today's Posts Mark Forums Read

Reply
 
LinkBack Thread Tools Display Modes
Old   October 30, 2013, 05:35
Default Problems rhoCentralFoam
  #1
Senior Member
 
Join Date: Dec 2010
Posts: 135
Rep Power: 6
eRzBeNgEl is on a distinguished road
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: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:
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 is offline   Reply With Quote

Old   October 30, 2013, 12:13
Default
  #2
Senior Member
 
Join Date: Dec 2010
Posts: 135
Rep Power: 6
eRzBeNgEl is on a distinguished road
any idea or help?
eRzBeNgEl is offline   Reply With Quote

Old   October 31, 2013, 16:35
Default
  #3
Member
 
Join Date: Jun 2012
Posts: 48
Rep Power: 5
maHein is on a distinguished road
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
maHein is offline   Reply With Quote

Reply

Thread Tools
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are On
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
Problems with rotational cyclic boundaries TReviol OpenFOAM Meshing & Mesh Conversion 8 July 11, 2014 03:45
[ICEM] Flow channel meshing problems StefanG ANSYS Meshing & Geometry 19 May 15, 2012 06:44
[ICEM] Problems with coedge curves and surfaces tommymoose ANSYS Meshing & Geometry 0 August 5, 2011 16:02
Two-phase air water flow problems by activating Wall Lubrication Force challenger85 CFX 5 November 5, 2009 06:44
Help required to solve Hydraulic related problems aero CFX 0 October 30, 2006 12:00


All times are GMT -4. The time now is 08:18.