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

dpm-DEFINE_DPM_LAW

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

Reply
 
LinkBack Thread Tools Display Modes
Old   September 9, 2011, 21:52
Default dpm-DEFINE_DPM_LAW
  #1
New Member
 
Join Date: Mar 2011
Posts: 13
Rep Power: 6
tomas is on a distinguished road
hi,
I would like to simulate particle heating in the plasma torch,The total
heat flux on the particle surface in the plasma Q is described using the heat transfer attributable to the plasma flow and the radiation loss as
Q=\pi*(dp)^2*h*(T-TP)-\epsilon\sigma((TP)^4-(Ta)^4)
h is the heat transfer coefficient T is plasma temperature ,TP is particle temperature,epsilonis the emissivity,sigma is Stefan–Boltzmann con-
stant,Ta is the ambient temperature.dp is particle diameter
Q can also be written for the particle temperature change and the phase changes such as melting
Q=\pi/6*(\rho*cp*(dp)^3*dTP/dt(tp<Tm,Tm<TP<Tb)................................ .........(1)
cp is specific heat,Tm is melting point ,Tb is boiling point

Q=\pi/6*(\rho*Hm*(dp)^3*dx/dt(tp=Tm)......................................... ....................(2)
Hm is the latent heat,rho is particle mass
Once the particle is completely molten (x = 1), its temperature is allowed to follow again Eq 1
My question is: I'm using UDFs to define the the heat of fusion during the transition,if I use DPM in fluent ,how can i use DEFINE_DPM_LAW Macro
to discribe this Eq ,
Thx
tomas is offline   Reply With Quote

Old   September 12, 2011, 07:28
Smile
  #2
New Member
 
Join Date: Mar 2011
Posts: 13
Rep Power: 6
tomas is on a distinguished road
UDF that uses DEFINE_DPM_SWITCH to
switch between DPM laws using a criterion. The UDF switches to DPM_LAW_USER_1 which refers to condenshumidlaw since only one user law has
been defined.if i have two ,Any suggestions of how to tackle this issue?
tomas

Last edited by tomas; September 14, 2011 at 09:04.
tomas is offline   Reply With Quote

Old   September 13, 2011, 06:46
Default
  #3
New Member
 
Join Date: Mar 2011
Posts: 13
Rep Power: 6
tomas is on a distinguished road
nobody help me!!!!!!!!!
tomas is offline   Reply With Quote

Old   October 5, 2011, 11:31
Default
  #4
New Member
 
Join Date: Mar 2011
Posts: 13
Rep Power: 6
tomas is on a distinguished road
#include "udf.h"
#include "dpm.h"
#include "mem.h"
#define Ta 300
#define Em 0.8
#define Tm 1983
#define sgb 0.0000000567
#define Tb 2500
#define Cpp 800
#define Hm 369000
#define Hv 12480000
DEFINE_DPM_LAW(HeatLaw,p,ci)
{
real area, rel_vel, Re, Pr, HTC, delta_temp,delta_x,delta_p;
float x = 0;
cphase_state_t *c = &(p->cphase);
area = 4.0 * M_PI * pow(P_DIAM(p),2.0);
Pr = c->sHeat * c->mu / c->tCond;/* Prandtl number */
HTC =c->tCond * (2.0 + 0.6*pow(p->Re,0.5)*pow(Pr,1./3.))/ P_DIAM(p);
if (P_T(p) < Tm)
{
delta_temp = P_DT(p) * (area * HTC * (c->temp-P_T(p))-area*Em*sgb*(pow(P_T(p),4.0)-pow(Ta,4.0))) /(P_MASS(p)*Cpp);
P_T(p) = P_T(p) + delta_temp;
}
else if (P_T(p) = Tm && x>O&& 1>x )
{
delta_x=P_DT(p) *6* (area * HTC * (c->temp-P_T(p))-area*Em*sgb*(pow(P_T(p),4.0)-pow(Ta,4.0)))/(P_RHO(p)*Hm);
x=x + delta_x;
}

else if((P_T(p) < Tb )&&(P_T(p) >Tm ))
{
delta_temp = P_DT(p) * (area * HTC * (c->temp-P_T(p))-area*Em*sgb*(pow(P_T(p),4.0)-pow(Ta,4.0))) /(P_MASS(p)*Cpp);
P_T(p) = P_T(p) + delta_temp;
x=1;
}
else
{
P_T(p)=Tb;
delta_p=-P_DT(p) *2* (area * HTC * (c->temp-P_T(p))-area*Em*sgb*(pow(P_T(p),4.0)-pow(Ta,4.0)))/(P_RHO(p)*Hv);

P_DIAM(p)= P_DIAM(p)+delta_p;


}

}
tomas is offline   Reply With Quote

Old   October 5, 2011, 11:39
Angry
  #5
New Member
 
Join Date: Mar 2011
Posts: 13
Rep Power: 6
tomas is on a distinguished road
why it showed no change in melting process (Tp=tm )
tomas is offline   Reply With Quote

Old   November 21, 2011, 09:11
Default
  #6
New Member
 
Join Date: Mar 2011
Posts: 13
Rep Power: 6
tomas is on a distinguished road
.................
tomas is offline   Reply With Quote

Old   November 21, 2011, 09:47
Default
  #7
Senior Member
 
ghost82's Avatar
 
Daniele
Join Date: Oct 2010
Location: Italy
Posts: 915
Rep Power: 15
ghost82 will become famous soon enough
Sorry I can't help you but can you check this?

else if (P_T(p) = Tm && x>O&& 1>x )
ghost82 is offline   Reply With Quote

Old   November 25, 2011, 20:49
Default
  #8
New Member
 
Join Date: Mar 2011
Posts: 13
Rep Power: 6
tomas is on a distinguished road
hank you very much ghost82. I do apprecaite your concern. All the best,

Tomas
tomas is offline   Reply With Quote

Old   August 29, 2013, 05:20
Default DPM in combination with plasma
  #9
New Member
 
Ian Maes
Join Date: Oct 2011
Posts: 10
Rep Power: 5
ian.maes is on a distinguished road
Hej Tomas!

This thread is already more than 2 years old and unfortunately I cannot help you with your problem (if that would still be necessary...), but I am trying to model plasma gasification using the DPM model.

I read that you also combine a plasma flow with the DPM model, so I was hoping you could help me with something.

I have written several UDFs to calculate the high-temperature thermodynamic and transport properties of the plasma flow, but now I find that this results in a conflict with using the DPM model...
As far as I know, it's not possible to use a UDF for density, but only the ideal-gas-law when using combusting paricles.

Did you manage to solve this problem in some way?

Thank you in advance!

Ian
ian.maes 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
DPM and UDS mighelone FLUENT 0 June 3, 2011 07:27
DPM modeling Mohsin FLUENT 0 March 31, 2010 22:11
DPM with UDF - Step-by-Step Procedure???? Prashanth FLUENT 3 April 3, 2009 17:45
DPM - do the particles affect the liquid? Nikhil Dani FLUENT 0 January 1, 2009 12:58
slow after mesh for enhanced wall treat DPM joshkemp FLUENT 0 May 1, 2007 17:20


All times are GMT -4. The time now is 04:16.