| 查看: 1717 | 回复: 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的编写。这些都存在一些问题，烦请大家予以指导！或者推荐一些可供学习的资料、文献也行！非常感谢！！！  |
回复此楼» 猜你喜欢
人气不行了
已经有11人回复
-
2026国自然函评费到账
已经有22人回复
-
导师吐槽：我怎么摊上了这么个极品研究生！
已经有3人回复
-
2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战
已经有7人回复
-
今天基金会出结果吗？20260819
已经有18人回复
-
我面上完蛋了
已经有14人回复
-
范进中举一文的中心思想
已经有8人回复
-
只有每年这种时候来逛逛小木虫
已经有27人回复
-
明天应该可查了！？
已经有5人回复
-
如果此刻你正在为国基感到焦虑，不妨来听听这首《基金之外》
已经有7人回复
2楼2017-06-27 15:05:33
3楼2021-02-21 15:50:03
20