| ²é¿´: 885 | »Ø¸´: 5 | ||
äìÒ¶·ÉÑï½ð³æ (ÕýʽдÊÖ)
|
[ÇóÖú]
udf±àÒë´íÎó
|
|
ÕâÊÇÔ´´úÂ룺 /*************************************************************** UDF for customizing the default Syamlal drag law in Fluent ****************************************************************/ #include "udf.h" #define pi 4.*atan(1.) #define diam2 3.e-4 DEFINE_EXCHANGE_PROPERTY(custom_drag3,cell,mix_thread,s_col,f_col) { Thread *thread_g, *thread_s; real x_vel_g, x_vel_s, y_vel_g, y_vel_s, abs_v, slip_x, slip_y, rho_g, rho_s, mu_g, reyp, afac, bfac, void_g, vfac, fdrgs, taup, k_g_s; /* find the threads for the gas (primary) */ /* and solids (secondary phases) */ thread_g = THREAD_SUB_THREAD(mix_thread, s_col);/* gas phase */ thread_s = THREAD_SUB_THREAD(mix_thread, f_col);/* solid phase*/ /* find phase velocities and properties*/ x_vel_g = C_U(cell, thread_g); y_vel_g = C_V(cell, thread_g); x_vel_s = C_U(cell, thread_s); y_vel_s = C_V(cell, thread_s); rho_g = C_R(cell, thread_g); rho_s = C_R(cell, thread_s); mu_g = C_MU_L(cell, thread_g); /*compute slip*/ slip_x = x_vel_g - x_vel_s; slip_y = y_vel_g - y_vel_s; abs_v = sqrt(slip_x*slip_x + slip_y*slip_y); /*compute Reynold's number*/ reyp = rho_g*abs_v*diam2/mu_g; /* compute particle relaxation time */ taup = rho_s*diam2*diam2/18./mu_g; void_g = C_VOF(cell, thread_g);/* gas vol frac*/ /*compute drag and return drag coeff, k_g_s*/ afac = pow(void_g,4.14); if(void_g<=0.85) bfac = 0.281632*pow(void_g, 1.28); else bfac = pow(void_g, 9.076960); vfac = 0.5*(afac-0.06*reyp+sqrt(0.0036*reyp*reyp+0.12*reyp*(2.*bfac- afac)+afac*afac)); fdrgs = void_g*(pow((0.63*sqrt(reyp)/ vfac+4.8*sqrt(vfac)/vfac),2))/24.0; k_g_s = (1.-void_g)*rho_s*fdrgs/taup; return k_g_s; } µ¥»úËĺ˲¢ÐмÆËã±àÒëʱÌáʾ´íÎó£ºerror:C compiler cannot execute cpp-I¡°c:\fluent.inc\...", ÕâÖÖÇé¿ö¸ÄÔõôÐÞ¸ÄÄØ |
» ²ÂÄãϲ»¶
²ÄÁϵ÷¼Á
ÒѾÓÐ7È˻ظ´
»¯Ñ§¹¤³Ì085602 305·ÖÇóµ÷¼Á
ÒѾÓÐ10È˻ظ´
289Çóµ÷¼Á
ÒѾÓÐ15È˻ظ´
291 Çóµ÷¼Á
ÒѾÓÐ7È˻ظ´
274Çóµ÷¼Á
ÒѾÓÐ14È˻ظ´
±±¾©ÁÖÒµ´óѧ˶µ¼ÕÐÉú¹ã¸æ
ÒѾÓÐ3È˻ظ´
309Çóµ÷¼Á
ÒѾÓÐ5È˻ظ´
292Çóµ÷¼Á
ÒѾÓÐ8È˻ظ´
Çóµ÷¼Á
ÒѾÓÐ4È˻ظ´
Ò»Ö¾Ô¸ Î÷±±´óѧ ×Ü·Ö282 Ó¢ÓïÒ»62 Çóµ÷¼Á
ÒѾÓÐ3È˻ظ´

kylafree
ÖÁ×ðľ³æ (ÖªÃû×÷¼Ò)
- Ó¦Öú: 166 (¸ßÖÐÉú)
- ½ð±Ò: 17477
- É¢½ð: 883
- ºì»¨: 18
- Ìû×Ó: 6210
- ÔÚÏß: 1161.7Сʱ
- ³æºÅ: 1061079
- ×¢²á: 2010-07-20
- ÐÔ±ð: GG
- רҵ: Á÷ÌåÁ¦Ñ§
2Â¥2012-04-11 17:03:54
äìÒ¶·ÉÑï
½ð³æ (ÕýʽдÊÖ)
- Ó¦Öú: 2 (Ó×¶ùÔ°)
- ½ð±Ò: 566.5
- É¢½ð: 964
- ºì»¨: 5
- Ìû×Ó: 749
- ÔÚÏß: 849.6Сʱ
- ³æºÅ: 1152526
- ×¢²á: 2010-11-21
- ÐÔ±ð: MM
- רҵ: ÄÜÔ´»¯¹¤

3Â¥2012-04-11 18:26:26
kylafree
ÖÁ×ðľ³æ (ÖªÃû×÷¼Ò)
- Ó¦Öú: 166 (¸ßÖÐÉú)
- ½ð±Ò: 17477
- É¢½ð: 883
- ºì»¨: 18
- Ìû×Ó: 6210
- ÔÚÏß: 1161.7Сʱ
- ³æºÅ: 1061079
- ×¢²á: 2010-07-20
- ÐÔ±ð: GG
- רҵ: Á÷ÌåÁ¦Ñ§
4Â¥2012-04-11 19:15:17
äìÒ¶·ÉÑï
½ð³æ (ÕýʽдÊÖ)
- Ó¦Öú: 2 (Ó×¶ùÔ°)
- ½ð±Ò: 566.5
- É¢½ð: 964
- ºì»¨: 5
- Ìû×Ó: 749
- ÔÚÏß: 849.6Сʱ
- ³æºÅ: 1152526
- ×¢²á: 2010-11-21
- ÐÔ±ð: MM
- רҵ: ÄÜÔ´»¯¹¤

5Â¥2012-04-11 19:52:40
ÇïÉñºÃÀ÷º¦
гæ (СÓÐÃûÆø)
- Ó¦Öú: 0 (Ó×¶ùÔ°)
- ½ð±Ò: 323.5
- Ìû×Ó: 84
- ÔÚÏß: 62.4Сʱ
- ³æºÅ: 2413670
- ×¢²á: 2013-04-11
- ÐÔ±ð: GG
- רҵ: ´«ÈÈ´«ÖÊѧ
6Â¥2014-06-06 17:09:16













»Ø¸´´ËÂ¥
5