我在水库三维模拟中,定义水温仅随z坐标变化,udf如下:
#include"udf.h"
DEFINE_INIT(start_temp,domain)
{
cell_t c;
Thread *thread;
real xc[ND_ND];
thread_loop_c(thread,domain)
{
begin_c_loop_all(c,thread)
{C_CENTROID(xc,c,thread);
if(xc[2]>=0&&xc[2]<=91)
C_T(c,thread)=6.53256+11.46744*exp(-0.04*(91-xc[2]))+10*exp(-0.018*(91-xc[2]))*cos(0.5236*(-2.15+1.3*exp(-0.005*(91-xc[2]))))+273.15;
}
end_c_loop_all(c,thread)
}
}
但是导入udf并初始化后,发现static temp在y=0平面(出口处)的分布如下图1,x=2502平面(出水管的截面)的温度分布如图2,好像受到了壁面温度的影响。而且最大值也不太对,按编译的公式,当z=91时,温度应该是300.1758,但是云图中最大值只有299.93,很奇怪,请大神看一下udf是否有问题?
后来我发现total temp的分层却很明显,如图3。
想请问一下,static temp是因为壁面影响才分布呈这样的吗????我该怎么调整,从哪里入手,可以将static temp初始化后变成分层的??
![水库分层三维模拟中,水温分布好像受到了壁面的影响,请问怎么解决?]()
1.png
![水库分层三维模拟中,水温分布好像受到了壁面的影响,请问怎么解决?-1]()
2.png
![水库分层三维模拟中,水温分布好像受到了壁面的影响,请问怎么解决?-2]()
3.png |