急！流场问题，请各位提宝贵建议！

tcsjk

本人接触CFD不久，我在做气流绕林带（气场内有林带）问题时，有一些问题一直困扰我，请斑竹和各位大虾提出宝贵意见和建议：
    建模时，我将流场设置为二维，林带设置为矩形固体元（有一定宽度和高度，并且有一定的疏透度，也就是说有一定的阻力和透气性，与单纯的绕平板流不同，），使用标准的k－e模型，现在问题来了：
     1、这应该属于气固两相问题吧，所以设置材料时，是不是既要设置气体材料，又要设置固体元材料才行？我看到有些外文文献（不知道自己理解对了，呵呵）意思是说当林带结构远远小于气流场的时候，可以不必考虑固体元太多，只要全局考虑气流场，把固体元的阻力影响以动量方程源项体现。问题一：是不是只要设置气体材料即可？
    2、问题二：其他边界条件还好设置，主要就是固体元的边界问题，是要设置壁面类型的wall，还是作为内部表面边界的多孔介质阶跃，抑或就作为内部界面（interior）不管？
    3、由于入口截面速度不是恒定值，K 动能和 E耗散初始值也不恒定，已经有如下的公式表示：

    我也知道要用UDF来定义，但有关UDF函数的资料不多，虽然C语言没有问题，看了一下总觉得不对点（不知道有没有跟我一样。呵呵），还有就是方程中的源项，S＝－1/2*Cd*A*U*U（Cd为阻力系数，A为单位面积密度，U为速度），我想也要用到UDF吧。
     问题三：哪位高手能不能就上面的等式写一个程序出来，或者提示一下也行，或者能告诉我哪里有UDF的手册（时间比较紧，中文的最好）。在此谢谢啦！

tcsjk

我写了一段代码，请大家看看对不对
#include "udf.h"
DEFINE_PROFILE(inlet_x_velocity, thread, index) /*定义初始速度u=u1/k*log(y/y0)*/  
{
  real x[ND_ND];   
  real y;
  face_t f;
  real u1=0.202
  real  k=0.42
  real  y0=0.00053
  begin_f_loop(f, thread)
    {
      F_CENTROID(x,f,thread);
      y = x[1];
    F_PROFILE(f, thread, index) = u1/k*log(y/y0);
      end_f_loop(f, thread)
}
DEFINE_PROFILE(K_velocity, thread, index) /*定义初始K值，它与初始速度有关k=1/2*u*1.5%*/  
{
  
   face_t f;
   Thread *t;
  begin_f_loop(f, thread)
   F_PROFILE(f, thread, index) = 1/2*F_U(f,t)*0.015
  end_f_loop(f, thread)
}