24小时热门版块排行榜

返回列表

ersu

铁虫 (初入文坛)

应助: 0 (幼儿园)
金币: 4
帖子: 3
在线: 6.3小时
虫号: 5317928
注册: 2016-12-05
性别: MM
专业: 水力学与水信息学

[求助] 利用fluent模拟水库水温分层

我做的是水库的温度分层模拟，水面表层单元需计入水气热交换，现在编写表面传热及透过水面的净短波辐射（源项）的udf，我知道源项是在边界条件的fluid中导入，但是表面传热及水库水面的边界条件布置如何设置，以下是我编写的udf，希望各位大神帮助帮助，水面边界条件该怎么设置。另外我的udf是否有问题，导入进行计算总是不计算几步就停了。谢谢各位啦
（1）表面传热
#include "udf.h"
#define S 182.3 /* 晴天到达地面的太阳总辐射量 */
#define GAMMA 0.03 /* 水面反射率 */
#define BETA 0.65 /* 水体表面吸收率 */
#define EW 0.97 /* 水体辐射系数 */
#define ETA 0.5 /* 太阳辐射在水体中的衰减系数 */
#define SIGMA 5.67e-08 /* Stefan-bolzaman常数 */
#define TA 25.8 /* 水面以上2m处的大气温度 */
#define CR 0.6 /* 云层覆盖率 */
#define W 1.6 /* 水面上10m处的风速 */
#define TD 18.0 /* 露点温度 */
DEFINE_PROFILE(surface_heat_transfer,thread,index)
{
   real ac,a,es,ea,fw,x[ND_ND];
   real an,sn,br,e,c,ts;
   face_t f;
   ac=1-0.216*exp(-0.000074*TA*TA);
   a=ac*(1+0.017*CR*CR);
   ts= F_T(f,thread)-273.15; /* 水表面温度 */
   fw=9.2+0.46*W*W;/* 风函数 ,反映自由对流和强迫对流对蒸发的影响*/
   es=6.11*pow(10,7.5*ts/(ts+273.3));/*相应于水面温度TS的紧靠水面的空气饱和蒸发压力*/
   ea=6.11*pow(10,7.5*TD/(TD+273.3));/*水面上空气的蒸发压力*/
sn=S*(1-GAMMA)*BETA;/* 太阳短波辐射 */
an=(1-0.03)*SIGMA*a*pow(273.3+TA,4);/*大气长波辐射*/
br=SIGMA*0.97*pow(273.3+ts,4);/*水体长波的返回辐射*/
e=fw*(es-ea);/*水面蒸发热损失*/
c=0.47*fw*(ts-TA);/*热传导通量*/
begin_f_loop(f,thread)
  {
F_CENTROID(x,f,thread);

F_PROFILE(f, thread, index) = sn+an-br-e-c; /* 将通过水面进入水体的热通量存贮在F_PROFILE中 */

  }
  end_f_loop(f,thread)
}
（2）透过水面的净短波辐射
#include "udf.h"
DEFINE_SOURCE(cell_x_source,cell,thread,dS,eqn)
{
  real x[ND_ND];
  real z;
  real sn,source;
  C_CENTROID(x,cell,thread);
  z=x[2];
  sn=S*(1-GAMMA)*BETA;/* 太阳短波辐射 */
  source=(1-BETA)*sn*exp(-ETA*z);
  return source;
}

回复此楼