基于PETSc的区域分解法数值模拟圆柱势流绕流问题.pdfVIP

第九届全国水动力学学术会议暨第二十二届全国水动力学研讨会文集 基于PETSc的区域分解 圆柱势流绕流问题 王吉飞 万德成 (．1：海交通人学船舶海洋‘j建筑T程学院海洋T程国家蕈点实验室200030) T00lkitfor scientific 行了并行数值计算，验证了网格的分辨率、计算精度、并行加速比，以及不同区域划分形式 对数值计算效果的影响，计算得到的流场速度分布和物体表面压力分布结果与解析解比较， 两者吻合很好。数值计算结果表明，基于PETsc的区域分解法和有限元法结合可以高效率实 现理想不可压流体圆柱势流绕流问题的并行数值计算，为今后开展复杂黏性流体流动问题的 并行计算奠定了基础． 关键词：势流；有限元法；钝体绕流；PETsc；区域分解；并行计算 l引言 圆柱绕流在实际j【程应用中非常普遍。在船舶与海洋-』：程中，洋流对柱形桥墩的冲击， 海洋平台立柱、立管绕流等均为圆柱绕流问题，因此研究圆柱绕流具有重要的实用意义。当 前计算能力不足是CFD发展所面临的一火挑战。为解决：f程应用方面的大规模设计计算问题， 并行化是大幅度提高计算效率的最有效手段。区域分解法是应用较广泛的一种粗粒度并行计 算方法。在流体力学有限元计算中，区域划分可在网格生成之前进行，也可在网格生成之后 进行，后者的优点是可以方便地对串行程序作并行改造而无需对网格生成形式进行修改，也 ToolkitforScientjfic 展的并行库，非常适合偏微分方程离散后人型稀疏代数方稗组的并行数值求解。 本研究基丁．有限元法和并行计算：L具箱PETsc，开发和编制区域分解法程序，数值求解 理想不可压流体中的圆柱绕流问题。首先对整个计算流域的网格进行区域划分，各子区域由 重替过渡区域连接，然后对各子区域网格用有限单元离散，单元形式为任意人面体单元，单 Inte而ce)库函数 (GMRES方法加lLU预估条件)进行求解，并使用MPl(Message P嬲sing 对各子区域间并行计算数据进行通信。本研究结构安排如下：首先描述圆柱势流绕流的区域 分解数值算法，然后给出数值模拟结果，并进行并行性能分析，最后给出结论。 ．270． 第九届全国水动力学学术会议暨第二十二届全国水动力学研讨会文集 2数值算法 2．1有限元空间离散 对于三维理想不可压流体的流动，可求解以速度势缈为未知数的Laplace方程边值问题： 。 p∥产7‘“(训㈡∈力 f鲁+窖+鲁：o I良2。勿2‘瑟2 { ^ (1) h项 矧。=g 其中Q为计算流体区域，S为本质边界，为自然边界。方程(1)的弱解积分表达式为： ㈤… 加。苏锄 却却 瑟瑟。 坞一7～ 蜷挈+等掣+警掣瑚=协啦 本研究使用任意六面体单元离散计算流域，采用八结点的La可绷ge线性插值函数进行单元 分析。进行单元分析时，为了使不同形状和尺寸的任意六面体单元均能使用规则的插值函数进 行单元分析，首先要进行坐标变换，将笛卡尔坐标(x，少，z)变换为无量纲的局部坐标(孝，7，f)， 值函数为： (f=1，…，8)