fluent中壁面怒谢尔特数的计算问题

songsh0915

最近在模拟管内流体层流流动传热的时候，需要查看壁面的局部怒谢尔特数，在默认的情况下，在查看wall flux下的Nu时，沿管长方向居然有严重的跳跃，与实际的单相流体的Nu严重不符。后来，查看相关资料，设了参考温度、管径参数后，局部怒谢尔特数沿管长方向的曲线倒是象那么回事，但我要看的是局部怒谢尔特数，在管长方向的每个位置的壁温和流体温度都是不一样的，参考温度该怎么设啊？
  简而言之，就是怎么看沿管长的壁面局部局部怒谢尔特数？
fluent中自身对Nu有定义：Nu=q.D/[(Tw-Tref).K],其中q为壁面热流密度，Tw为壁面温度；Tref为设置的参考温度，是个常数；K为流体导热系数。如果要求局部Nu的话，Tref应该是沿程处处不同的啊，也即应该是沿程流体的温度，此时Nux=q.D/[(Tw,x-Tf,x).K]，下标x表示沿管长的具体位置。所以我实际需要的局部Nu与fluent自身的Nu还是有一定差别，请问这个问题该如何解决？
Tf,x求起来实际比较麻烦的，理论上应该是每个截面流体的平均温度。
如果我要自定义的话，Tw,x和Tf,x改怎么定义呢？