ÇóÖú£¬ÎÒ×Ô¼º±àµÄudf³ÌÐò£¬Ê¹ÓÃÁËdefine adjustºÍudsÔ´ÏÔÚfluentÖбàÒëbuildºÍ¼ÓÔØload¶¼Ã»ÓÐÎÊÌ⣬µ«ÊǼÆËãÒ»¿ªÊ¼¾ÍÏÔʾuds-0·¢É¢ÁË£¬²»ÖªµÀÊÇʲôÔÒò£¬²¿·Ö´úÂëÈçÏ£º
#include "udf.h"
#define PI 3.14159 /* ¶¨Òå¦Ð*/
#define K 1.38e-23 /*¶¨Òå²£¶û×ÈÂü³£Êý */
#define M 2.99e-26 /*¶¨ÒåË®·Ö×ÓÖÊÁ¿*/
#define R 461.5 /*¶¨ÒåË®ÕôÆøÆøÌå³£Êý*/
#define f_v 0.4 /*Ïà¶Ôʪ¶È*/
real S1, G, J, pv1, ps1, T, P, sur_tension, density_c, density_v;/*²»±¥ºÍ¶È£¬¼ª²¼Ë¹ìÊ£¬³ÉºËÂÊ£¬Ë®ÕôÆø·ÖѹÁ¦£¬Ë®ÕôÆø±¥ºÍѹÁ¦£¬Î¶ȣ¬Ñ¹Á¦*/
/*¼ÆËã³ÉºËÂÊ£¬ÒºµÎÉú³¤Ëٶȣ¬ÁÙ½ç°ë¾¶*/
DEFINE_ADJUST(chenghelv, d)
{
Thread *t;
cell_t c;
thread_loop_c (t,d)
{
begin_c_loop_all(c,t)
{
T = C_T(c,t);
P = C_P(c,t);
ps1 = -2.88934 + 0.37966 * T - 0.00925 * pow(T, 2) + 1.54e-4 * pow(T, 3);/*±¥ºÍѹÁ¦ÄâºÏ*/
density_v = -0.00134 + 0.00136 * T - 3.27e-5 * pow(T, 2) + 7.76e-7 * pow(T, 3);/*ÕôÆøÃܶÈÄâºÏ*/
density_c = 1000.07904 + 0.01231 * T - 0.00586 * pow(T,2) + 1.6e-5 * pow(T,3);/*Ë®ÃܶÈÄâºÏ*/
sur_tension = 0.07566 - 1.41e-4 * T - 2.68e-7 * pow(T,2);/*±íÃæÕÅÁ¦ÄâºÏ*/
S1 = P / ps1;
pv1 = P * f_v;
G = 16 / 3 * PI * pow((M / (density_c * log(S1) * K * T)),2) * pow(sur_tension, 3); /*¼ÆË㼪²¼Ë¹×ÔÓÉìÊ*/
C_UDMI(c,t,0) = sqrt(2 / PI * sur_tension) * pow(M ,-3/2) * pow(density_v, 2) / density_c * exp( -G / K * T);/*³ÉºËÂÊ*/
end_c_loop_all(c,t)
}
}
}
/*µÚÒ»¸öÊäËÍ·½³ÌQ0Ô´Ïî*/
DEFINE_SOURCE(q0_source,c,t,dS,eqn)
{
return C_UDMI(c,t,0);
}
/*µÚ¶þ¸öÊäËÍ·½³ÌQ1Ô´Ïî*/
DEFINE_SOURCE(q1_source,c,t,dS,eqn)
{
return C_UDMI(c,t,0) * pow(C_UDMI(c,t,2),2) + C_R(c,t) * C_UDSI(c,t,0) * C_UDMI(c,t,1);
}
/*µÚÈý¸öÊäËÍ·½³ÌQ2Ô´Ïî*/
DEFINE_SOURCE(q2_source,c,t,dS,eqn)
{
return C_UDMI(c,t,0) * pow(C_UDMI(c,t,2),2) + 2 * C_R(c,t) * C_UDSI(c,t,1) * C_UDMI(c,t,1);
}
/*µÚËĸöÊäËÍ·½³ÌgµÄÔ´Ïî*/
DEFINE_SOURCE(g_source,c,t,dS,eqn)
{
return 4 * PI * density_c * (1/3 * C_UDMI(c,t,0) * pow(C_UDMI(c,t,2),3) + C_R(c,t) * C_UDSI(c,t,2) * C_UDMI(c,t,1));
}
/*ÄÜÁ¿·½³ÌÔ´Ïî*/
DEFINE_SOURCE(energy_source,c,t,dS,eqn)
{
real dt, temp, L;
dt = RP_Get_Real("physical-time-step"
temp = C_T(c,t);
L = 2502 - 2.42169 * temp;
return L * C_UDSI(c,t,3) / dt;
} |