2　数值解法

本文采用混合有限分析法[5]对变换后的模型进行离散，用压力修正算法（SIMPLE方法）进行求解。

混合有限分析方法是在空间一维局部单元线性化常微分方程求精确解的基础上，按局部线性化算子迭加原理得到求解高维问题的计算模式。它避免了一般有限分析法系数中级数项出现，具有明显的自动迎风性质，数值频散与耗散要比有限差分方法低，计算稳定性好。

对于方程（5），将时变项并入源项中考虑，并取，f=J2/Γ，则式（5）写成

设变换平面上的ɑ和Z方向的步长分别为hi和kj，则式（6）的混合有限分析五点格式为

其中，。

流场初始条件取水库全场静止，水温水质分布参考实测值给定，采用ADI（交替方向隐式迭代）方法（TDMA算法）交替对两个方向扫描，两次全场扫描完成一轮迭代。