|
[Sponsors] |
January 28, 2013, 11:48 |
Temperature Dependant Porosity in Combustion
|
#1 |
Member
NormalVector
Join Date: Oct 2010
Posts: 71
Rep Power: 15 |
Hello,
I'm trying to model the temperature dependent porosity of a porous medium in a combustion application. I'm using a UDF to loop over the cells in the porous medium and assign a porosity value based on the cell temperature. I don't believe my UDF is wrong so I'm not posting in the UDF section of the forums but I'm including it here for completeness: __________________________________________________ __________ #include "udf.h" DEFINE_PROFILE(porosity_udf_test,t,nv) { real temperature; real a; cell_t c; begin_c_loop(c,t) { temperature = C_T(c,t); a = 0.5*((exp(2*temperature - 1000)-1)/(exp(2*temperature - 1000)+1)) + 0.5; C_PROFILE(c,t,nv) = a; } end_c_loop(c,t) } __________________________________________________ __________ If you notice, I'm trying to approximate a step in porosity (ie, below one temperature the porous medium is there and above it it is not) using a hyperbolic tangent function centered around 500. When I run my simulation, it's fine until the porous media heats up and then suddenly Pressure Correction diverges. Here are my relevant solver settings:
Last edited by NormalVector; January 28, 2013 at 12:57. |
|
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
monitoring point of total temperature | rogbrito | FLUENT | 0 | June 21, 2009 17:31 |
High temperature methane+air | Peter | FLUENT | 5 | January 26, 2009 18:04 |
Combustion temperature | yu_korogi | OpenFOAM Running, Solving & CFD | 1 | May 6, 2008 04:00 |
Shell Model for Low Temperature Combustion | Phil | Main CFD Forum | 0 | March 17, 2006 00:34 |
Xoodles temperature change without combustion | hannes | OpenFOAM Running, Solving & CFD | 5 | October 17, 2005 16:31 |