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

rhoCentralFoam - Problem high difference pressure

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

Reply
 
LinkBack Thread Tools Display Modes
Old   June 11, 2012, 09:29
Default rhoCentralFoam - Problem high difference pressure
  #1
New Member
 
Join Date: May 2012
Posts: 1
Rep Power: 0
Nicolas_OF is on a distinguished road
Hi,

I'm new to OpenFoam (2.1.0). I try to simulate underexpanded axisymmetric supersonic flow of Argon. I used tutorial model : LadenburgJet60Psi.My simulation seems to work when difference between inlet and chamber pressure is small (factor of 10). But when this ratio becomes of the order of 1000, I get the famous error "floating exception point". I attempt different parameters, boundaries and I don't success to find where is my error.

Can Openfoam solve this kind of problem ?


Thank you and sorry for my English.




Here my parameters :

BlockMesh :

convertToMeters 0.001;

vertices
(
( 0 0 0)
(18 0 0)
( 0 1 -0.008726655121)
(18 1 -0.008726655121)
( 0 9 -0.017453310242)
(18 9 -0.017453310242)
( 0 1 0.008726655121)
(18 1 0.008726655121)
( 0 9 0.017453310242)
(18 9 0.017453310242)
);

blocks
(
hex (0 1 3 2 0 1 7 6) (60 10 1) simpleGrading (1 1 1)
hex (2 3 5 4 6 7 9 8) (60 10 1) simpleGrading (1 1 1)
);

edges
(
);

boundary
(
inlet
{
type patch;
faces
(
(0 2 6 0)
);
}

outlet
{
type patch;
faces
(
(1 3 7 1)
(3 5 9 7)
);
}

freestreamInlet
{
type patch;
faces
(
(2 4 8 6)
);
}

freestream
{
type patch;
faces
(
(4 8 9 5)
);
}

wedge1
{
type wedge;
faces
(
(0 2 3 1)
(2 4 5 3)
);
}

wedge2
{
type wedge;
faces
(
(0 1 7 6)
(6 7 9 8)
);
}
);

mergePatchPairs
(
);





ThermoPhysicalProperties :

thermoType ePsiThermo<pureMixture<sutherlandTransport<specieT hermo<hConstThermo<perfectGas>>>>>;

mixture
{
specie
{
nMoles 1;
molWeight 39.948;
}
thermodynamics
{
Cp 520;
Hf 30;
}
transport
{
As 2.054e-06;
Ts 162.4;
}
}


p :


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

internalField uniform 101;

boundaryField
{
inlet
{
type fixedValue;
value uniform 271724;
}

outlet
{
type waveTransmissive;
field p;
phi phi;
rho rho;
psi psi;
fieldInf 101;
gamma 1.67;
lInf 0.025;
value uniform 101;
}

freestream
{
type totalPressure;
value uniform 101;
p0 uniform 101;
U U;
phi phi;
rho none;
psi psi;
gamma 1.67;
}

freestreamInlet
{
type zeroGradient;
}

wedge1 {type wedge;}
wedge2 {type wedge;}
}


T :

dimensions [0 0 0 1 0 0 0];

internalField uniform 98.0;

boundaryField
{
inlet
{
type fixedValue;
value uniform 299.1;
}

outlet
{
type zeroGradient;
}

freestream
{
type totalTemperature;
value uniform 97;
T0 uniform 97;
U U;
phi phi;
rho none;
psi psi;
gamma 1.67;
}

freestreamInlet
{
type fixedValue;
value uniform 97.0;
}

wedge1 {type wedge;}
wedge2 {type wedge;}
}


U :

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

internalField uniform (0 0 0);

boundaryField
{
inlet
{
type fixedValue;
value uniform (515.6 0 0);
}

outlet
{
type inletOutlet;
inletValue uniform (0 0 0);
value uniform (0 0 0);
}

freestream
{
type zeroGradient;
}

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

wedge1 {type wedge;}
wedge2 {type wedge;}
}
Nicolas_OF 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
Extremely high Temperature needed to simulate pressurization with rhoCentralFoam shangzung OpenFOAM Running, Solving & CFD 0 August 6, 2010 06:35
Tutorial problem: "rhopSonicFoam - shockTube" with high pressure shangzung OpenFOAM 0 November 1, 2009 14:07
high pressure simulation mahi FLUENT 0 November 7, 2008 02:30
Does star cd takes reference pressure? monica CD-adapco 1 April 19, 2007 11:26
Multicomponent fluid Andrea CFX 2 October 11, 2004 05:12


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