Question

Fluent 用UDF求出口平均温度的问题。

一个圆管中一端流入一端流出，壁面加热，要把出口的平均温度赋给入口
udf如下：
real NV_VEC(A);
real avg_temp;
real sum_T_A=0.0;
real sum_A=0.0;
Thread *thread_out;
face_t f;
Domain *domain;
DEFINE_ADJUST(adjust,d)
{domain=Get_Domain(1);
thread_out=Lookup_Thread(domain,4);/4是出口面的id/
begin_f_loop(f,thread_out)
{ F_AREA(A,f,thread_out);
sum_A+=NV_MAG(A);
sum_T_A+=NV_MAG(A)*F_T(f, thread_out);
}
end_f_loop(f,thread_out)

avg_temp=sum_T_A/sum_A;

begin_f_loop(f,t)
{
F_PROFILE(f,t,i) = avg_temp;
}
end_f_loop(f,t)
}

其中的

sum_T_A+=NV_MAG(A)*F_T(f, thread_out);
这句话搞不清楚是什么意？
面积乘以温度吗？
求好心人解答，不胜感激。
奉上10分，了表寸心。

Answer 1

计算平均值的方法一般有几何平均、算术平均以及加权平均，而在以上代码中，计算的平均值算是加权平均，其基本思想是：
1）假设各界面面积标记为Ai，界面温度为Ti，其总面积元为A=sum(Ai)；
2）假设各面积元的加权量为xi=Ai/A;
3）按加权平均量计算的平均值为Tm=xi*Ti
那么，转化为代码而言，sum_T_A+=NV_MAG(A)*F_T(f, thread_out)表示：
sum_T_A=sum_T_A+NV_MAG(A)*F_T(f, thread_out)

追问

wind_cool2009：
您好，非常感谢给于的帮助。
另外，我想做个壁面气泡在不同流体速度下的形态的研究，用fluent怎么把气面做在壁面上呢，单气泡的模型。
我之前直接patch 一个和壁面相交的气泡，貌似fluent没认出来，只是按照自由气泡显示了半个气泡而已，
可以给点意见吗？

追答

看了你的提问，估计你主要是追踪在不同流体速度下，气液界面分布
fluent本身提供了很多界面追踪方法，比如VOF,Level-set方式，如果从这个角度来考虑的话，可以在patch里面设置成第二相，但是这样有个问题——组分之间的扩散问题的考虑，最终算的结果可能是气泡融合到流体里面去了而无法表现去形态特征；
我建议的考虑方法：
为了研究气泡形态，忽略扩散问题，假设气液分界面为弹性界面，当在不同流体速度下时，其受力不同因而其形态也不同。根据此思路，我们可以采用动网格技术，初始的时候，把气液界面设为壁面，当流体作用时，简化为界面受力移动，然后到达平衡。
具体问题具体看，希望对你有帮助，Good luck!