| 查看: 2975 | 回复: 10 | ||||
[求助]
VOF体积分数的梯度为什么为0呢？？
|
|
分别使用如下两个UDF计算VOF中体积分数的梯度，两只情况下梯度值都是0？ 请高手指教，UDF是否有什么问题？ cas文件及UDF在附件中。 UDF2: #include "udf.h" DEFINE_ON_DEMAND(on_demand_calc) { Domain *subdomain; /* declare domain pointer since it is not passed as an argument to the DEFINE macro */ Thread *cell_thread; cell_t cell; float a1,a2,x[2]; subdomain=Get_Domain(2); /* returns phase with ID=2 domain pointer£¬meanning primary phase*/ /* Loop over all cell threads in the domain */ thread_loop_c(cell_thread,subdomain) { /* Loop over all cells */ begin_c_loop(cell,cell_thread) { C_UDSI(cell,cell_thread,0)=C_VOF(cell,cell_thread); } end_c_loop(cell,cell_thread) begin_c_loop(cell,cell_thread) { C_CENTROID(x,cell,cell_thread); C_UDMI(cell,cell_thread,0)=NV_MAG(C_UDSI_G(cell,cell_thread,0)); a1=C_VOF(cell,cell_thread); a2=C_UDMI(cell,cell_thread,0); printf("location:%f,%f,vof:%f,gradient:%f\n",x[0],x[1],a1,a2); } end_c_loop(cell,cell_thread) } } UDF3： #include "udf.h" DEFINE_ON_DEMAND(grad) { Domain *mixture_domain, *sub_domain; Thread *cell_thread; cell_t cell; int phase_domain_index; float a1,a2,x[2]; mixture_domain=Get_Domain(1); sub_domain_loop(sub_domain,mixture_domain,phase_domain_index) { /* loop if primary phase */ /* if (DOMAIN_ID(sub_domain) == 2) */ /* loop over all cell threads in the primary phase domain */ thread_loop_c(cell_thread,sub_domain) { /* loop over all cells in primary phase cell threads */ begin_c_loop_all(cell,cell_thread) { C_CENTROID(x,cell,cell_thread); C_UDSI(cell,cell_thread,0)=C_VOF(cell,cell_thread); C_UDMI(cell,cell_thread,0)=NV_MAG(C_UDSI_G(cell,cell_thread,0)); a1=C_VOF(cell,cell_thread); a2=C_UDMI(cell,cell_thread,0); printf("location:%f,%f,vof:%f,gradient:%f\n",x[0],x[1],a1,a2); } end_c_loop_all(cell,cell_thread) } } }  水-空气相分布 |
回复此楼» 本帖附件资源列表
-
欢迎监督和反馈：小木虫仅提供交流平台，不对该内容负责。
本内容由用户自主发布，如果其内容涉及到知识产权问题，其责任在于用户本人，如对版权有异议，请联系邮箱：xiaomuchong@tal.com - 附件 1 : FLUENTWORK.rar
2012-06-06 21:35:30, 35.79 K
» 收录本帖的淘帖专辑推荐
 udf相关 | 传热传质 |
» 猜你喜欢
职称评审没过，求安慰
已经有56人回复
-
最近几年招的学生写论文不引自己组发的文章
已经有5人回复
-
26申博自荐
已经有3人回复
-
A期刊撤稿
已经有4人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- VOF射流问题，为什么显示不出来两相啊
已经有11人回复
- 一般做复合材料，增强颗粒占基体的体积分数是多少？
已经有4人回复
- 利用CSF模型将气液相界面的表面张力添加到动量方程源项的相关问题
已经有17人回复
- 用VOF模型如何模拟不用流速的气液两相流
已经有4人回复
- fluent的vof两相流问题
已经有7人回复
- UDF中如何编写体积分数的梯度
已经有9人回复
- VOF模型静止壁面上速度不为0
已经有10人回复
- 求助三维vof模拟，最大网格尺寸0.5mm，结果气泡发散没有明显边界，
已经有20人回复
- FLUENT加入VOF模型后的一些问题
已经有15人回复
- Fluent欧拉模型中相的体积分数的设定
已经有16人回复
- 采用VOF模型对磨削加工中的流场进行仿真，仿真结果不对，摸索中，望高手相助！！
已经有28人回复
- fluent 6.3 13.0版本中多相流的比较
已经有20人回复
- vof模型中获取第二项的体积分数
已经有6人回复
- VOF模型中，体积分数云图中液体和气体之间总存在一个边界层，如何去掉？
已经有14人回复
- multiphase+相变UDF-by（vof-hrbeu)
已经有165人回复
- 管内气液两相流的模拟
已经有12人回复
- VOF，MIXTURE和双流体，在做多相流时这些模型怎么选择呢
已经有3人回复
- 【求助】多相流mixture模型计算液体两相流动边界条件中体积分数如何设置？【已解决】
已经有9人回复
2楼2012-06-07 10:05:06
3楼2012-06-07 14:55:03
4楼2012-06-07 14:56:17
5楼2012-06-07 16:19:49
6楼2012-06-07 16:20:36
7楼2012-06-07 21:55:13
【答案】应助回帖
★ ★
xiegangmai: 金币+2, 谢谢参与！ 2012-06-07 23:40:12
xiegangmai: 金币+2, 谢谢参与！ 2012-06-07 23:40:12
|
#define BOUNDARY_ID_fluid 2 /*全局变量，对整个流体区域*/ /*************/ DEFINE_ON_DEMAND(on_demand_calc) { Domain *subdomain; Thread *t; cell_t c; float a1,a2,x[ND_ND],xx,yy; FILE *fpc; domain=Get_Domain(1); /*基本相是不是默认1*/ t=Lookup_Thread(domain,2) /*混合相指针*/ t1=THREAD_SUB_THREAD(t,0); /*一相*/ t2=THREAD_SUB_THREAD(t,1);/*另一项*/ /* Loop over all cell threads in the domain */ thread_loop_c(c,t) { C_CENTROID(x,c,t); xx=x[0]; yy=x[1]; C_UDSI(c,t,0)=C_VOF(c,t); C_UDMI(c,t,0)=NV_MAG(C_UDSI_G(c,t,0)); a1=C_UDSI(c,t,0); fpc=fopen("gradient.txt","a"  ;fprinth(fpc,"%f %f %f\n", xx,yy,a1); fclose(fpc); } end_c_loop(c,t) |
8楼2012-06-07 22:29:13
9楼2012-06-07 22:31:45
10楼2014-04-24 21:01:24