动静圆柱坐标系下N—S方程的分区算法及子域干涉处理研究.PDF

第25卷第2期 空气动力学学报 Vd．25。No．2 2007年06月 ACTAAERoIn烈AMⅡCASDⅡCA Jun．，2007 文章编号：0258．1825(2007)02．0211-05 动静圆柱坐标系下N—S方程的分区算法 及子域干涉处理研究 姜 健，刘 波，王掩刚，靳 军 (西北工业大学翼型、叶栅空气动力学国防科技重点实验室，陕西西安710072) 摘要：耦合了相对和绝对圆柱坐标系下三维Navier-Stokes方程，统一了动子域和静子域的控制方程；在子域计算 中，采用高分辨率的三阶ENN格式及LU—SGS隐式解法；在子域干涉处理中，研究并对比了动静子域干涉面的重 合处理与重叠处理；最终发展了一种新的既高效又准确的分区算法。通过对某带进口导叶的三级跨声速轴流压气 机叶栅内三维紊流流场的数值模拟，验证了该分区算法的有效性，计算结果与实验数据吻合较好。 关键词：分区算法；子域干涉；Navier-Stokes方程；圆柱坐标系 中图分类号：V211．3 文献标识码：A 模型和动静子域干涉面处理模型的时间推进调用实 0 引 言 现了一种新的串行分区算法，构成了一种既准确又高 对于相对和绝对圆柱坐标系共存计算问题，例如 多级叶轮机械内流场的数值模拟，Navier-Stokes方程数值求解方法。最后，数值模拟了某带进口导叶的三 的求解还局限在如下的方法上[1I，即在已有的单子域 级跨声速轴流压气机叶栅内三维紊流流场，验证了该 计算方法的基础上采用统一的网格，并且对动子域和 分区算法的有效性，计算结果与实验数据吻合较好。 静子域采用不同的坐标系，对不同的子域采用不同的 1动静圆柱坐标系下的N．S方程 方程，然后统一求解。显然这类方法要求大内存的大 型计算机，其发展也受到了限制。从而对分区算法的 本文将待求解速度量统一为绝对速度分量，从而 需求就日益迫切[2，3|。分区算法具有以下的优点：在 统一了相对和绝对圆柱坐标系下的控制方程。通过 每一子域内，容易生成高质量的网格；对于不同的子 张量形式的严格推导，并将相对圆柱坐标系(9，r，石) 域采用形式上统一的控制方程，以方便计算程序的模 块化设计；对于串行分区算法，在任何时刻计算机内 下，其形式为： 只需存放一个子域的计算数据，这就放宽了对计算机 aQaEaFaG aEvaF’， aGy= 内存的要求。 本文从上述理论出发，耦合了相对和绝对圆柱坐 (1) 标系下的三维Navier-Stokes方程，使得动静子域中使 甲 0 用统一的控制方程，而通过旋转角速度的选取来确定 rpu 一』D口M+r秽 1 是动子域还是静子域；研究对比了动静子域干涉面的 ，)一上 例 ； 五=了lD2+P—r仰 、一J J 重合处理和重叠处理，并发展为干涉面处理模型；在 rpw 0 子域计算中，采用高精度高分辨率的三阶ENN格式 rpe 0 以保证对