【UDF贴壁气泡的动量源项,添加就出错.另：x,y两个方向源项,用解释形式分别添加,只显示后添加的源项.垂直壁上附着一气泡,考虑tensionsurfaceforce和阿基米德力后,建立动量源项./#include"udf.h"#incl】

UDF贴壁气泡的动量源项,添加就出错.另：x,y两个方向源项,用解释形式分别添加,只显示后添加的源项.

垂直壁上附着一气泡,考虑tensionsurfaceforce和阿基米德力后,建立动量源项.

/

#include"udf.h"

#include"sg_mphase.h"/*包括体积分数宏C_VOF(C,T)*/

DEFINE_SOURCE(ymom_source,c,pri_th,dS,index)

{

realvof_g;/*定义单元格气相体积分数变量*/

Thread*mix_th;/*定义计算区域的指针*/

realc_yar,yar;

face_tf;

realNV_VEC(A);

realv_m,v0,v0_m,bv_m;/*bv边界区气相体积*/

realsource;

mix_th=THREAD_SUB_THREAD(pri_th,1);/*指向混合区主相即气相的指针*/

begin_c_loop(c,mix_th)

{

vof_g=C_VOF(c,mix_th);

if(vof_g>0.&&vof_g