Fluent 通量（ρAU）计算不守恒

kate-huzhizi

通量（ρAU）计算中遇到不守恒问题（A 是面积；U是面上的法向速度），困扰了一礼拜了，大伙帮我看看吧，不胜感激

我用fluent做了一个不可压缩流场稳态计算，离散格式用的是一阶迎风

收敛之后我想在流场里面取一个假想的体（长方体）来监测组成这个体的各面上的通量情况 （模型如图）；

步骤：
1.在fluent里面用surface-plane surface功能把组成这个长方体的各个面给定义出来
（一共三个面；其余三个面是建筑物表面和地板被定义成wall的没有质量交换因此不考虑）
2.然后用file-export-solution data 选择cell center 把各surface上的速度存出来
3. 求surface上个各个小网格面上的速度U
    小网格面上的速度由这个面相邻的两个cell中心的速度去计算：
对两边的cell中心的速度求算术平均（左右两边cell中心离面的距离相等）
4.对各面进行ρAU计算；
最后对各面的通量加和发现入流通量和出流通量不相等不知道什么问题？

大家给我支支招吧，谢谢谢谢

[ 本帖最后由 kate-huzhizi 于 2012-10-11 16:30 编辑 ]