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

Temperature Dependant Porosity in Combustion

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

Reply
 
LinkBack Thread Tools Display Modes
Old   January 28, 2013, 12:48
Default Temperature Dependant Porosity in Combustion
  #1
Member
 
NormalVector's Avatar
 
NormalVector
Join Date: Oct 2010
Posts: 72
Rep Power: 6
NormalVector is on a distinguished road
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:
  • Pressure Based, Transient
  • Realizeable k-e
  • Finite Rate/Eddy Dissipation combustion
  • PISO pressure-velocity coupling (I've tried all of the other P-V couplings)
  • Standard pressure discretization (I've tried all other pressure discretizations)
  • the other discretizations are 1st order
Any ideas? I was originally using a full tetrahedral mesh but after I made the porous region structured it seemed to help. Not sure what the problem is.

Last edited by NormalVector; January 28, 2013 at 13:57.
NormalVector 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
monitoring point of total temperature rogbrito FLUENT 0 June 21, 2009 17:31
High temperature methane+air Peter FLUENT 5 January 26, 2009 19: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 01:34
Xoodles temperature change without combustion hannes OpenFOAM Running, Solving & CFD 5 October 17, 2005 16:31


All times are GMT -4. The time now is 02:19.