写出了多孔介质内不可压缩流体流动的基本方程组,利用流线迭代法对充满多孔介质的变截面圆管内的流体的流动进行数值计算,并通过算例对充满多孔介质管道中的流动情况进行分析.

采用三阶Ｔａｙｌｏｒ－Ｇａｌｅｒｋｉｎ（ＴＧ）有限元法离散Ｎ－Ｓ方程，并对离散格式中由高阶时间导数转化而来的附加项进行了缩并运算修正，使其起到了人工粘性的作用。由于该法隐含了流线迎风的耗散作用，可较好地应用于较高Ｒｅｙｎｏｌｄｓ数流动的计算。同时采用压力校正法求解不可压缩流场的各原始变量，并讨论了压力校正Ｐｏｉｓｓｏｎ方程的边界条件；最后给出了二维方腔拖带流的计算结果，证明了算法的正确性。