CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Fluent UDF and Scheme Programming

UDF for pellets

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

Reply
 
LinkBack Thread Tools Display Modes
Old   September 27, 2011, 10:35
Default UDF for pellets
  #1
New Member
 
Zsolt Vaszi
Join Date: Jul 2010
Posts: 8
Rep Power: 7
ZsoltiV is on a distinguished road
Hi everyone,

I want to simulate with eulerian multiphase model one gasifier. The gasifier is just a cylinder in whitch I have 3 phases. The first is the air need to gasification, wood pellets and ash. In this case I want to know just I the air is flows around the pellets. But I need an UDF for changeing the pellets to ash in time. I wrote this UDF. If some can please correct it:

#include "udf.h"


#define OMEGA 7.26E-3
#define K_F0 0.0257
#define K_S0 1.0
#define d_b 0.00683 /* mean diameter of bulk material */
real rho_b = 2250.; /* particle density */
real cp_b = 835.; /* specific heat capacity */
real A_V_sphere = 6.0/d_b; /* area-to-volume ratio of a sphere*/
real A_b_V; /* area-to-volume ratio */


real
k_solid(real e_gas)
{
real k_sol;
real a = K_S0/K_F0;
real b = 1.25*pow((1.-e_gas)/e_gas, 10./9.);
real t1 = (a-1.)/pow((1.-b/a),2.0)*(b/a)*log(a/b);
real t2 = -((b-1.)/(1.-b/a))-0.5*(b+1.);
real gm = (2./(1.-b/a))*(t1+t2);
real k_bs = K_F0*(OMEGA*a + (1.-OMEGA)*gm);
k_sol = k_bs/sqrt(1.-e_gas);
return (k_sol);
}

DEFINE_PROPERTY(conduct_gas,cell,thread)
{
real k_bf, k_gas;
k_bf = K_F0*(1.-sqrt(1.-C_VOF(cell,thread)));
k_gas = k_bf/C_VOF(cell, thread);
return (k_gas);
}


DEFINE_PROPERTY(conduct_solid,cell,thread)
{
real vol; /* cell volume */
real vol_m1; /* cell volume time step minus 1 */
real fac; /* factor in time derivative */
real fac_m1; /* factor at time step minus 1 */
real eps_b = 1.0 - C_POR(c,t); /*solid volume fraction*/

/* d(fac*V*phi)/dt = fac*V/dt*phi - fac_m1*V/dt*phi_m1
= apu * phi - su */
vol = C_VOLUME(c,t);
vol_m1 = vol;
fac = rho_b*cp_b*eps_b;
fac_m1 = fac;
*apu = -fac*vol / delta_t;
*su = fac_m1*vol_m1*C_UDSI_M1(c,t,i)/delta_t;
}
ZsoltiV 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
UDF parallel error: chip-exec: function not found????? shankara.2 Fluent UDF and Scheme Programming 1 January 16, 2012 23:14
How to add a UDF to a compiled UDF library kim FLUENT 3 October 26, 2011 21:38
Dynamic Mesh UDF Qureshi FLUENT 1 December 2, 2009 01:27
UDF...UDF...UDF...UDF Luc SEMINEL FLUENT 0 November 25, 2002 05:03
UDF, UDF, UDF, UDF Luc SEMINEL Main CFD Forum 0 November 25, 2002 05:01


All times are GMT -4. The time now is 07:53.