Working with an energymodel

Deranda · May 25, 2016, 06:56

Hi,

I start working with an energymodel for my simulation right now, just to get a boundary temperature for the bai-gosman collision model.
However, if I use the regular molecular weight of air (28.96 g/mol) the calculated density is 0.9 kg/m³. but it should be 8.238 kg/m³ (at 8.5bar and 86°C)
If I calculate the density with the thermodynamic equation of state
p*v=m*R*T
-> density = p/(T*R_L/M)
The calculated density is 8.238 kg/m³. Now I changed the molecular wight of air until the density fits. The result is 246 g/mol. To me it looks like the universal gas constant is different to 8.31 J/(mol*K)

Would be nice if you share your experiences with me

best regards

marmot · May 25, 2016, 07:11

Perhaps before changing MW, when you say density is that the cell-center value, using cell-center temperature, but perhaps you are looking at wall temperature and calculating a density based on this. Just an idea to check this,

Deranda · May 25, 2016, 11:02

well I've chosen 359.15K (86°C) as the inlet, outlet, boundary and initial temperature

marmot · May 25, 2016, 12:25

OK, did you set your reference pressure to 8.5bar and initialize with 0 static pressure, or did you initialize the pressure to be at 8.5bar, leaving the reference pressure at default, this could also be a reason,

May 25, 2016, 06:56	Working with an energymodel	#1
Deranda Member Nils Join Date: Nov 2015 Posts: 59 Rep Power: 11	Hi, I start working with an energymodel for my simulation right now, just to get a boundary temperature for the bai-gosman collision model. However, if I use the regular molecular weight of air (28.96 g/mol) the calculated density is 0.9 kg/m³. but it should be 8.238 kg/m³ (at 8.5bar and 86°C) If I calculate the density with the thermodynamic equation of state pv=mRT -> density = p/(TR_L/M) The calculated density is 8.238 kg/m³. Now I changed the molecular wight of air until the density fits. The result is 246 g/mol. To me it looks like the universal gas constant is different to 8.31 J/(mol*K) Would be nice if you share your experiences with me best regards

Similar Threads
Thread	Thread Starter	Forum	Replies	Last Post
Is the subsetMotionSolver working in parallel?	Arnoldinho	OpenFOAM Running, Solving & CFD	7	April 27, 2016 05:38
OpenFoam 1.7.1 in Ubuntu 10.04 won't start working	Leech	OpenFOAM Installation	8	August 10, 2014 21:42
DPM parallel is not working but serial is working	johnwinter	FLUENT	1	March 27, 2012 03:01
Working directory	mgonzalo	FLUENT	1	November 11, 2011 11:05
viewing cfx post while working on cfx solver manager	HMR	CFX	5	March 9, 2011 23:33

May 25, 2016, 07:11		#2
marmot Senior Member kevin alun Join Date: Sep 2011 Location: Germany Posts: 106 Rep Power: 15	Perhaps before changing MW, when you say density is that the cell-center value, using cell-center temperature, but perhaps you are looking at wall temperature and calculating a density based on this. Just an idea to check this,

May 25, 2016, 11:02		#3
Deranda Member Nils Join Date: Nov 2015 Posts: 59 Rep Power: 11	well I've chosen 359.15K (86°C) as the inlet, outlet, boundary and initial temperature

May 25, 2016, 12:25		#4
marmot Senior Member kevin alun Join Date: Sep 2011 Location: Germany Posts: 106 Rep Power: 15	OK, did you set your reference pressure to 8.5bar and initialize with 0 static pressure, or did you initialize the pressure to be at 8.5bar, leaving the reference pressure at default, this could also be a reason,