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

Question

UDF 贴壁气泡的动量源项,添加就出错.另：x,y两个方向源项,用解释形式分别添加,只显示后添加的源项.
垂直壁上附着一气泡,考虑tension surface force 和阿基米德力后,建立动量源项.
/
#include"udf.h"
#include "sg_mphase.h" /*包括体积分数宏C_VOF(C,T)*/
DEFINE_SOURCE(ymom_source,c,pri_th,dS,index)
{
real vof_g; /*定义单元格气相体积分数变量*/
Thread *mix_th; /*定义计算区域的指针*/
real c_yar,yar;
face_t f;
real NV_VEC(A);
real v_m,v0,v0_m,bv_m; /*bv边界区气相体积*/
real source;
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

Answer 1

首先,在程序本身就存在错误,v_m表述不正确,v0_m是密度乘以体积,而bv是体积分数乘以体积
标识/*42*/中,v未定义!
然后一个,你没有理解我上次的建议,如果按照现在的表达思路的话,应该是气泡运动以及在运动过程中掺混!
而我上次的建议则考虑的是以追踪界面变化,气泡在运动过程中的变形或者运动
具体的可以加我QQ,我没做过这方面,也想扩展一下应用技巧,交流一下,2208774663