| 查看: 2046 | 回复: 3 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
[求助]
利用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; } |
回复此楼» 猜你喜欢
基金申请
已经有44人回复
-
CSC与新西兰维多利亚大学PhD奖学金项目
已经有0人回复
-
物理学I论文润色/翻译怎么收费?
已经有117人回复
-
新西兰Robinson研究所 招聘CSC公派访问人员
已经有0人回复
-
帮我的英语口语老师找学生
已经有0人回复
-
什么时候开奖？
已经有13人回复
-
散金币祈福
已经有94人回复
-
青基已中
已经有1人回复
-
散金币祈福
已经有107人回复
-
求助，如何提取ELK的rt-TDDFT在某一时刻的自旋密度分布和电子密度分布
已经有1人回复
3楼2018-04-25 09:52:16
4楼2018-04-27 18:28:02
