船舶与海洋工程数值水池研究.pdfVIP

船舶与海洋工程数值水池研究 万德成 (上海交通大学 船舶海洋与建筑工程学院 海洋工程国家重点实验室 200240) 摘要 本文介绍一个新的船舶与海洋工程数值水池计算模型。“船舶与海洋工程数值水池”软件集成系统以 源代码工具 OpenFOAM 和 PETSc 为基础进行编程。前处理系统采用移动网格技术，实现复杂三维网格生 成，并与 Overset 网格技术结合，实现自适应网格和动网格的生成。后处理系统采用开发源代码POV-Ray 系统提供的图形库和辅助工具，进行后处理编程，实现数据后处理和流场动态三维演示。 关键词 数值水池，多重网格虚拟边界法，多尺度，OpenFOAM，PETSc ，POV-Ray 引 言 船舶与海洋工程数值水池是指利用计算机程序数值仿真实现物理船舶与海洋工程水池的功能， 包括造波、造流、波－波、波－流、波－流－结构物相互作用的实验。船舶与海洋工程数值水池的 优点十分明显：它节省大量场地、人力和费用，重复性好，无外界干扰，信息量丰富等，日益受到 船舶与海洋工程届的重视。特别是近年，随着计算机发展和计算方法的提高，以及持续增加的海洋 资源开发的需求，船舶与海洋工程数值水池研究已成为国际船舶与海洋工程学术界的研究热点，如 欧盟于 2005 年提出并正在实施的 VIRTUE （The Virtual Tank Utility in Europe ）计划。随着船舶CFD 软件对船舶性能计算和船型优化的重要作用，目前商业 CFD 软件已经应用到船流场计算，如 Fluent ， CFX，FEFLO ，CFDSHIP，SHIPFLOW 等都有很好应用例子，特别是 SHIPFLOW 和 CFDSHIP 是直 接针对船流场计算而设计的 CFD 软件，SHIPFLOW 还充分利用势流和粘性流的特性，把船流场分 为势流、边界层和粘性三个区域分别进行处理，然后进行三个区域的耦合计算；而 CFDSHIP 则直接 从粘性流出发，利用结构块网格处理复杂船边界问题和 Level-set 方法处理自由面问题。但这些软件 都有各自的局限性，无论是网格处理、求解器设计、后处理，还是计算结果可靠性分析上由于源代 码不公开而存在修正的困难，更为重要的是我们只能使用别人的软件，始终没有具有自己知识产权 的船舶性能计算软件，这与我们造船大国和强国的发展方向不符合，也不符合自主创新的科技发展 要求。 本文针对船舶与海洋工程复杂水动力学问题的复杂几何边界、运动界面、运动物体、多尺度、 时域非定常的特点，介绍一个新的船舶与海洋工程数值水池计算模型：利用多重网格虚拟边界法 (MFBM) 、嵌入边界/ 界面方法(Immersed boundary/interface method) 、ALE(Arbitrary Lagrangian- Eulerian)方法和移动网格技术、自适应时间步技术、分块嵌套技术(Overset)和并行技术，形成多尺度 网格嵌套流体数值解器(multi-scale mesh-nested CFD solver) 。对高雷诺数流动问题，采用多尺度大涡 模拟方法(multi-scale large eddy simulation)处理，即把速度场分解为多个不同尺度的运动，不同的尺 度对应一个或几个不同的物理机制，不同的物理机制采用不同的亚格子模型，不同的亚格子模型构 成的大涡模拟方程采用不同的数值方法和计算网格，这样多尺度大涡模拟方法就可以考虑多个不同 尺度上多个不同的物理量过程的耦合。为了在船舶与海洋工程数值水池中，反映实际海洋的风－浪 －流的真实环境，利用海洋中几个实际测点记录的风、浪、流数据，作为船舶与海洋工程数值水池 软件系统的初始边界条件，通过非线性非规则波相求解模拟方法(phase-resolved simulation method) 和多层优化迭代(multiple-level iterative optimization)，在空间域和时间域上重构真实海洋风－浪－流 的演化过程。建立这一海洋环境重构技术，不仅可以反映真实的海洋环境条件，而且通过利用物理 水池实验得到的部分测试数据作为数值水池的初始计算条件，进行数值水池软件计算，就可以大大 节约计算时间，也可以很好延续物理水池没有完成的实验。“船舶与海洋工程数值水池”软件集成系 统将