| ²é¿´: 1067 | »Ø¸´: 3 | ||
jimoxingzheгæ (³õÈëÎÄ̳)
|
[ÇóÖú]
Çó×é´óÏÀ¸øÎÒ¿´¿´±ÉÈ˱àдµÄUDF£¬¸Ð¼¤²»¾¡¡¡Ð¡µÜûÓнð±Ò°¡ÎØÎØ
|
|
#include "udf.h" #include "sg_mphase.h" #define A 16.0114 #define B 3456.80 #define C -8.67 #define L 0.00000125.0 #define POP 101325.0 #define K 0.1 #define C10H22_NUM 0 #define LIQ_PHASE_ID 0 #define GAS_PHASE_ID 1 DEFINE_SOURCE(liq_zhengkuiwan,c,liq_thread,dS,eqn) { int i; cell_t c; Thread*gas_thread,*mixture_thread; Domain*mixture_domain,*liq_domain,*gas_domain; Material *gas_material,*species; real c10h22_mole_fract,P_c10h22,total_mole,Mw[MAX_SPE_EQNS]; mixture_domain=Get_Domain(1); liq_domain=DOMAIN_SUB_DOMAIN(mixture_domain,LIQ_PHASE_ID); gas_domain=DOMAIN_SUB_DOMAIN(mixture_domain,GAS_PHASE_ID); //¸Ð¾õûʲôÓ㬵ÚÒ»´Î±àдudf¾Í¶¼Ð´ÉÏÁË mixture_thread=THREAD_SUPER_THREAD(liq_thread); gas_thread=THREAD_SUB_THREAD(mixture_thread,GAS_PHASE_ID); gas_material=THREAD_MATERIAL(gas_thread); real m_dot_liq,c_s,c_g,m_kuosan; real x[ND_ND],y=x[1]; real ps,v; //psΪҺµÎ±íÃæ±¥ºÍѹÁ¦ vΪµ¥¸öÍø¸ñµÄÌå»ý thread_loop_c(mixture_thread,mixture_domain) //ÔÚÕû¸ö¶þά¿ØÖÆÇø²éѯ»ìºÏÎïµÄµ¥ÔªÏß { begin_c_loop(c,mixture_thread) //ÔÚµ¥ÔªÏßÉϲéѯµ¥Ôª {if(0.01<=C_VOF(c,liq_thread)&&C_VOF(c,liq_thread)<1.0) //²éÑ¯ÆøÒº½çÃæ´¦µÄÍø¸ñ£¬ÎÒÑ¡ÁË0.01µ½1Ö®¼ä£¬ÒòΪµ±ÒºµÎµÄÌå»ý·ÖÊýСÓÚ0.01ʱ»á²úÉú·ÐÌÚÕô·¢ { ps=133*exp(A-B/(C_T(c,liq_th)+C); //AntoineÕôÆøÑ¹·½³ÌÇóµÃÒºµÎ±íÃæ±¥ºÍѹÁ¦ps c_s=ps/(C_RGAS(c,gas_thread)*C_T(c,liq_thread)); //ÇóµÃÒºµÎ±íÃæÅ¨¶È C_RGAS(c,gas_thread)ÆøÌå³£ÊýÕâÑù¿ÉÒÔÂ𣿣¿Ò²¿ÉÒÔÖ±½ÓÊäÈë8.3134 total_mole=0; mixture_species_loop(gas_material,species,i) { Mw=MATERIAL_PROP(species,PROP_mwi); total_mole+=C_YI(c,gas_thread,i)/Mw; //ÇóµÃÆøÒº½çÃæÍø¸ñÄÚÆøÌåµÄ×ÜĦ¶ûÊý } c10h22_mole_fract=(C_YI(c,gas_thread,C10H22_NUM)/Mw[C10H22_NUM])/total_mole; //ÇóµÃÆøÒº½çÃæÄÚc10h22µÄĦ¶û·ÖÊý P_c10h22=POP*c10h22_mole_fract; //ÇóµÃÆøÒº½çÃæÄÚc10h22µÄ·Öѹ c_g=P_c10h22/(C_RGAS(c,gas_thread)*C_T(c,liq_thread)); //ÇóµÃÆøÒº½çÃæÍø¸ñÄÚÆøÌåÖÐc10h22µÄŨ¶È m_kuosan=c_s-c_g; //ÆøÒº½çÃæÍø¸ñÄÚÒºµÎ±íÃæc10h22Ũ¶ÈºÍÆøÌåÖÐc10h22Ũ¶ÈµÄ²îÖµ£¬´Ó¶ø²úÉúÀ©É¢ v=L*L*L; m_dot_liq=142*K*y*m_kuosan*C_DIFF_EFF(c,gas_thread,i)/v*C_VOF(c,liq_th)*pow(1-C_VOF(c,liq_thread),2.0)*(1-(1-C_VOF(c,liq_thread))*c10h22_mole_fract); dS[eqn]= //×îÖÕÔ´Ïî dS»¹Ã»ÓÐÇ󣬱ȽÏÂé·³ } } end_c_loop(c,mixture_thread) return m_dot_liq; } [ Last edited by jimoxingzhe on 2013-4-9 at 16:14 ] |
» ²ÂÄãϲ»¶
340Çóµ÷¼Á
ÒѾÓÐ5È˻ظ´
ÉúÎïѧѧ˶Çóµ÷¼Á
ÒѾÓÐ10È˻ظ´
ÉϺ£µçÁ¦´óѧ²ÄÁÏ·À»¤ÓëвÄÁÏÖØµãʵÑéÊÒÕÐÊÕµ÷¼ÁÑо¿Éú£¨²ÄÁÏ¡¢»¯Ñ§¡¢µç»¯Ñ§£¬»·¾³£©
ÒѾÓÐ4È˻ظ´
²ÄÁÏѧÇóµ÷¼Á
ÒѾÓÐ6È˻ظ´
303Çóµ÷¼Á
ÒѾÓÐ5È˻ظ´
Ò»Ö¾Ô¸ÎäÀí085500»úеרҵ×Ü·Ö300Çóµ÷¼Á
ÒѾÓÐ7È˻ظ´
¿¼Ñе÷¼Á
ÒѾÓÐ4È˻ظ´
281Çóµ÷¼Á
ÒѾÓÐ4È˻ظ´
0805 316Çóµ÷¼Á
ÒѾÓÐ6È˻ظ´
085601Çóµ÷¼Á×Ü·Ö293Ó¢Ò»Êý¶þ
ÒѾÓÐ3È˻ظ´
» ±¾Ö÷ÌâÏà¹Ø¼ÛÖµÌùÍÆ¼ö£¬¶ÔÄúͬÑùÓаïÖú:
jimoxingzhe
гæ (³õÈëÎÄ̳)
- Ó¦Öú: 0 (Ó×¶ùÔ°)
- ½ð±Ò: 17.8
- Ìû×Ó: 4
- ÔÚÏß: 6.7Сʱ
- ³æºÅ: 2109486
- ×¢²á: 2012-11-06
2Â¥2013-04-09 17:20:17
jimoxingzhe
гæ (³õÈëÎÄ̳)
- Ó¦Öú: 0 (Ó×¶ùÔ°)
- ½ð±Ò: 17.8
- Ìû×Ó: 4
- ÔÚÏß: 6.7Сʱ
- ³æºÅ: 2109486
- ×¢²á: 2012-11-06
3Â¥2013-04-09 17:20:57
ÄäÃû
Óû§×¢Ïú (ÖøÃûдÊÖ)
- Ó¦Öú: 84 (³õÖÐÉú)
- ½ð±Ò: 5015.3
- É¢½ð: 1030
- ºì»¨: 35
- Ìû×Ó: 2148
- ÔÚÏß: 235.4Сʱ
- ³æºÅ: 0
- ×¢²á: 2011-10-29
- ÐÔ±ð: GG
- רҵ: ÒÁ˹À¼½Ì
4Â¥2013-04-13 22:51:05













»Ø¸´´ËÂ¥
