MPS-FEM数值分析带自由面的流固耦合问题

来源：计算力学学报杂志2017年第3期北京时间：

作者：张友林,邹璐,万德成

单位：上海交通大学 船舶海洋与建筑工程学院 海洋工程国家重点实验室 高新船舶与深海开发装备协同创新中心, 上海 20040,上海交通大学 船舶海洋与建筑工程学院 海洋工程国家重点实验室 高新船舶与深海开发装备协同创新中心, 上海 20040,上海交通大学 船舶海洋与建筑工程学院 海洋工程国家重点实验室 高新船舶与深海开发装备协同创新中心, 上海 20040

摘要：随着计算科学的发展，研究人员为探索流固耦合问题的物理机理而提出了众多的数值方法。其中，耦合的移动粒子半隐式方法MPS（Moving Particle Semi-Implicit method）和有限单元法FEM（Finite Element method）为流固耦合问题的数值仿真工作提供了新的途径。本文所有流场的数值模拟工作均采用课题组自主开发的无网格法求解器MLParticle-SJTU来完成。该求解器在原始的MPS法基础上，对核函数、压力梯度模型、压力泊松方程的求解和自由面判断方式等方面进行了改进。此外，在该求解器框架内，基于FEM法拓展了针对结构场进行求解的功能。首先，对MPS和FEM方法的理论模型及其耦合策略进行了介绍。然后，采用该自研MPS-FEM耦合求解器，数值模拟了溃坝流动对弹性结构的冲击及其相互作用的标准问题。通过将结构变形及自由面波型变化等结果与已发表结果进行对比，验证了该求解器在处理带自由面剧烈变化的粘性流体和柔性变形结构的耦合作用问题上的可行性。

关键词：无网格法,MPS方法,有限元法,MPS-FEM,流固耦合,自由面流动,MLParticle-SJTU求解器,

基金资助：国家自然科学基金（51379125，51490675，11432009，51579145）；长江学者奖励计划（2014099）；上海高校特聘教授（东方学者）岗位跟踪计划（2013022）；上海市优秀学术带头人计划（17XD1402300）；上海市船舶工程重点实验室基金（K2015-11）；工信部数值水池创新专项VIV/VIM项目（2016-23/09）资助项目