| 查看: 1699 | 回复: 2 | |||
[求助]
fluent uds方程对流项问题 用于求解level set方程及重新初始化方程 已有1人参与
|
|
大家好！我最近在尝试基于fluent实现level set方法追踪自由表面。不利用fluent自带的coupled level set and vof，而是希望通过uds求解level set方程、重新初始化方程。最后结合level set函数值求解物理量的控制方程。目前已利用uds求解level set函数，在求解重新初始化方程时遇到问题。 首先，在求解level set方程时，uds方程源项为0，扩散系数为0，对流项选择mass flow rate. 求解的方程如图1：  图1 利用UDS求解的Level Set方程 重新初始化方程如图2所示，对照UDS方程形式，认为对流项为w =sign phi0 *grad phi / |grad phi|  图2 重新初始化方程 对流项UDF编写时，令w0=sign phi0 /|grad phi|; w = w0*grad phi. 如下： DEFINE_UDS_FLUX(flux_reinil,f,t,i) { real w=0.0,w0; cell_t c0; Thread *t0; real NV_VEC(psi_vec),NV_VEC(A); c0=F_C0(f,t); t0=F_C0_THREAD(f,t); F_AREA(A,f,t); w0 = C_UDSI(c0,t0,0)/sqrt(C_UDSI(c0,t0,0)*C_UDSI(c0,t0,0)+0.01*0.01)/sqrt(C_UDSI_G(c0,t0,1)[0]*C_UDSI_G(c0,t0,1)[0]+C_UDSI_G(c0,t0,1)[1]*C_UDSI_G(c0,t0,1)[1]+C_UDSI_G(c0,t0,1)[2]*C_UDSI_G(c0,t0,1)[2]); NV_D(psi_vec,=,C_UDSI_G(c0,t0,1)[0],C_UDSI_G(c0,t0,1)[1],C_UDSI_G(c0,t0,1)[2]); NV_S(psi_vec,*=,w0); w = NV_DOT(psi_vec,A); return w; } 可以编译，可以计算，但是结果不能实现重新初始化。 问题多多:1）利用UDS解重新初始化方程的可行性；2）如果可行，对流项形式是否是w；3）DEFINE_UDS_FLUX的编写。这些都存在一些问题，烦请大家予以指导！或者推荐一些可供学习的资料、文献也行！非常感谢！！！  |
回复此楼» 猜你喜欢
filecode=后面第一个是大写字母
已经有7人回复
-
时间戳又变了8-15
已经有22人回复
-
欢迎发来filecode的Mz6后的代码验证其规律
已经有60人回复
-
有时候，自然基金真的不能太认真 （我的申报经验）
已经有10人回复
-
filecode
已经有12人回复
-
小木虫上这么多卖论文的，真有人买论文么？感觉没必要啊
已经有13人回复
-
哪位老哥知道今年的国自然具体哪一天放榜？
已经有7人回复
-
关于Filecode分析方法
已经有13人回复
-
各位道友，我要去昆明玩几天，回来见。
已经有8人回复
-
咱们一起用铁证分析2026国家社科基金中标与否
已经有22人回复
2楼2017-06-27 15:05:33
3楼2021-02-21 15:50:03
