复杂流场的多块并行数值模拟.pdfVIP

复杂流场的多块并行数值模拟 杨树池乔志德 西北工业大学翼型研究中心 摘要：本文着重介绍了三维多块并行数值计算程序EU3DMBP的多块网格计算结构，及在sol的 种复杂外形进行计算的结果，并与实验结果做了对比． 1 引言 计算机的飞速发展及计算格式的不断改进，使得CFD(计算流体力学)越来越广泛的应 用于飞行器性能分析及飞行器设计领域中。在实际应用中，CFD也不断克服新困难，不断 发展完善。但是，系统的来说，一套CFD软件能否有效的解决实际问题．应用于工程设 计，除了前后置处理之外，流场计算软件自身也要满足一定的要求： 速流动，能准确捕捉激波。 2)高效率(Effective)对定常问题，能在较短的时间内收敛到稳态解。 角、侧滑角，能适应各种复杂的外形。 计算格式的不断改进．使得计算结果的准确程度和计算程序的适应性等方面有了较大 程度的提高。但计算机硬件发展远远不能满足CFD工作人员的需求，因而近些年发展起 来的并行化网络数值计算为CFD发展提供了新的动力，利用现有的计算机条件，在网络 环境进行并行数值计算，可大大提高计算效率，缩短计算时间。 本文就是利用现有的计算条件进行并行计算的一个尝试，在SGI工作站上，分别采 P 用两个，三个和四个C U进行了计算．并就其并行效率问题进行了讨论。 2 控制方程及多块网格结构 本文所介绍的三维多块流场计算程序EU3DMBP的发展本身就是面向工程应用而开发 内外流动的无粘流场，计算格式空间采用中心有限体积格式进行离散．时间采用Runge． Kutta简化多步格式积分。 245 为提高计算格式的精度及正确捕捉激波，格式的人工粘性除包含二四阶混合导数的 Scalar Model之外，还包含矩阵粘性MatrixScheme及满足LED Dissipation Dissipation (LocalExtremum and Diminishing)条件的CUSP(ConvectiveUpstreamSplitPressue)格式 及各种改进措旌。同时为提高EU3DMBP计算效率，采用了当地时间步长、焓阻尼、变系 “柚月啊jb 数的隐式残值光顺等多种加速收敛措施。 正如前面所提到的，为使该软件能分折计算各种复杂流外形场，EU3DMBP除了采用 有限体积离散外，还采用了多块网格结构，即将整个流场划分成一系列任意大小的区域， 在保证各区域交界面上网格连续的情况下，各区域可独立生成任意拓扑结构的网格，以准 确模拟外形、真实反映流动特征，同时避免了单块网格的拓扑结构问题。多块网格计算过 程中主要的问题是相邻块之间的信息传递，由于各块网格可采用任意的拓扑结构．其边界 各异t包含物面边界、远场边界、对称边界、割缝边界、奇性边界、进气边界、排气边界、 出口边界等，其中相邻面的割缝连接方式非常复杂．有八种之多．见附图1。EU3DMBP 格与相邻块的单元相重合，适时更新各圈网格的物理量，可以独立求解每一块流场。 3 并行计算 EU3DMBP的另外一个显著特点在于其采用了并行处理，利用SGI的POWER FORTRAN语言．对原有的计算程序加以修改，可以在POWERCHALLEDGE工作站上oq。q瑚馕Ⅺ堋油门镧H罾!涠用墒j，…—j0■Ⅶ羽坷—Ⅻ“一曲i，㈡0：!女o Data)方法，由 用多CPU并行计算，EU3DMBP采用的是SPMD(SingleProgramMultiple 于POWER 部进行。 MIPS POWER 的个数，将每一个DO循环分成相应的几部分，由不同的CPU分别执行不同的部分，以下 面的求和为例： SUM；0．0 DO I=【．N SUM=