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 *domai
来源:学生作业帮助网 编辑:作业帮 时间:2024/07/06 20:44:38
![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 *domai](/uploads/image/z/4868095-31-5.jpg?t=Fluent+%E7%94%A8UDF%E6%B1%82%E5%87%BA%E5%8F%A3%E5%B9%B3%E5%9D%87%E6%B8%A9%E5%BA%A6%E7%9A%84%E9%97%AE%E9%A2%98.%E4%B8%80%E4%B8%AA%E5%9C%86%E7%AE%A1%E4%B8%AD%E4%B8%80%E7%AB%AF%E6%B5%81%E5%85%A5%E4%B8%80%E7%AB%AF%E6%B5%81%E5%87%BA%2C%E5%A3%81%E9%9D%A2%E5%8A%A0%E7%83%AD%2C%E8%A6%81%E6%8A%8A%E5%87%BA%E5%8F%A3%E7%9A%84%E5%B9%B3%E5%9D%87%E6%B8%A9%E5%BA%A6%E8%B5%8B%E7%BB%99%E5%85%A5%E5%8F%A3udf%E5%A6%82%E4%B8%8B%EF%BC%9Areal+NV_VEC%28A%29%3Breal+avg_temp%3Breal+sum_T_A%3D0.0%3Breal+sum_A%3D0.0%3BThread+%2Athread_out%3Bface_t+f%3BDomain+%2Adomai)
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 *domai
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);
这句话搞不清楚是什么意?
面积乘以温度吗?
了表寸心.
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 *domai
计算平均值的方法一般有几何平均、算术平均以及加权平均,而在以上代码中,计算的平均值算是加权平均,其基本思想是:
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)