案例分享｜复杂仿真应用定制——TankSimulator液体灌箱仿真APP

1月前浏览1362

矩阵版头.gif

（一）问题及背景描述

罐箱是由箱体框架和压力容器罐体组成的设备，具有广泛用途和高效空间利用率，运输成本低廉。罐箱已被国际标准化，可用于公路、铁路、水路等交通方式。全球市场对罐箱需求巨大，其中我国是最大生产国，因此罐箱仿真模拟存在巨大市场需求。

对靖江某压力容器生产制造企业的产品设计制造及运行场景进行分析，总结其需要通过仿真解决的问题是，罐箱在运输过程中受到的冲击作用导致内部液体晃动并对罐箱结构造成冲击，易造成罐体结构损伤需要在设计之初即对相关场景进行仿真模拟，并优化设计。这是典型的流固耦合问题，内部液体冲击（流体部分还是气液两相流问题）可能引起应力集中、大变形以及与罐箱缺陷等结合问题导致失稳破坏，对安全性产生影响。

企业目前存在的问题主要有以下几点：

1.使用商业仿真前后处理软件、仿真求解器软件，需要支付高昂的License费用，企业前期投入会非常高。

2.目前没有专门用于罐体冲击仿真的软件，只能使用通用仿真软件，这就要求工程师掌握较多理论知识和软件使用经验，而流固耦合仿真需要熟悉结构和流体仿真，进一步提高了企业人力成本。

3.仿真过程需要硬件资源的支撑，流固耦合计算量也比常规仿真计算量大很多，企业自主采购搭建仿真集群成本高昂，且需要专人运维。

针对以上问题，国家超级计算无锡中心先进制造部基于一系列开源代码进行了二次开发和自主研发，以开展定制化应用开发。在本期文章中，就我们的定制应用开发平台深度集成开源仿真程序CalculiX、OpenFoam以及preCICE后形成的流固耦合罐箱设计TankSimulator APP进行详细介绍。

（二）解决方案

什么是流固耦合问题

流固耦合现象在工程中有许多实际案例，如航空领域的气动弹性问题，桥梁领域的颤振问题，都是流固耦合存在的场景。以模拟仿真的角度看，流-固耦合（以下简称FSI）是分别描述流体动力学和结构力学的定律之间的多物理场耦合。这种现象的特点是变形结构或运动结构与周围或内部的流体流动之间的相互作用，这种相互作用既可以是稳定的，也可以是振荡的[1]。

流固耦合问题的求解面临以下几个主要难点：

1.物理场的相互作用：流体和固体之间发生相互作用的过程往往十分复杂。

2.交界面确定及网格生成：FSI问题的求解需要确定流体和固体之间的交界面。这通常涉及到交界面的几何形状、位置和边界条件等问题。

3.耦合算法：FSI求解器需要设计有效的耦合算法来解决流体和固体之间的相互作用。这包括将流体和固体的方程耦合起来，处理界面条件，以及确保数值求解的稳定性和收敛性。

4.流体和固体求解在自身领域中已经具有一定的挑战性，而在FSI问题中更加复杂。流体求解中常见的问题包括数值不稳定性、边界条件的处理和网格质量等。而固体求解中需要考虑非线性材料行为、接触与摩擦等因素。

5.计算资源和计算效率：FSI问题通常涉及大规模计算和复杂的物理过程，对计算资源的需求较高。解决大规模FSI问题可能需要使用高性能计算平台，并优化计算算法和策略，以提高计算效率和减少计算成本。

我们的解决方案

在介绍解决方案之前，首先介绍一下TankSimulator APP的仿真场景。仿真对象是一个圆柱形罐箱，内部带有防波板，用于缓解流体晃动程度，几何模型示例如下图所示：

模拟的是压力容器在载具刹车过程中和刹停后内部液体自由晃动两个阶段中，内部液体的晃动过程及其对容器壁面冲击导致的容器壁面发生变形，这是一个双向FSI模拟仿真。

求解器选择及调试

结构固体力学求解器采用CalculiX，它是一个对标著名非线性结构分析商软ABAQUS的免费、开源的3D结构非线性分析程序[2]。本算例采用动力学求解模块，封头、容器壁、方波板采用壳单元，不考虑几何非线性。

流体力学求解器采用OpenFOAM，它是应用最广的开源计算流体力学软件，具有较好的二次开发特性。本算例采用的是interFoam求解器，基于VOF方法进行两相流仿真。

针对仿真场景，调试优化仿真参数设置，并通过实际算例进行验证，最终将与产品无关的仿真参数固化在求解器配置文件中，并实现全流程仿真自动化脚本开发。

耦合软件选择及开发

通过充分调研了两款仿真求解器的代码结构和运行特性，最终选择preCICE作为耦合软件，并基于官方适配器对求解器和适配器进行了二次开发，提高两个求解器的耦合效率。（适配器是用于在流体和固体之间实现耦合关系的组件或接口。它允许流体和结构两个求解器进行信息交换和数据传递，以模拟流体和固体之间的相互作用，包括了流体和固体的边界条件、位移、速度、压力等参数互相传递的标准接口。）

流固耦合方案决定了各求解器的解是如何组合的。一般分为串行/并行显式、串行/并行隐式四种耦合方案。压力容器自由晃动过程，属于几何结构简单、非稳态、较大时间步长的仿真问题，显式方案可以提供较好的计算效率和数值稳定性。因此，我们选择串行/并行显式耦合方案。

我们也对耦合器关键环节的参数及模式进行了定义：

1.数据映射。在多域仿真时，各个域通常独立的离散化，这样交接面处网格会出现不一致现象。因此，耦合器需要在不同网格之间进行数据插值、外推或者重构，以便流体和固体模拟可以在不同网格间进行信息传递。我们采用最近邻点插值方法将流体网格面上的信息传递到固体模拟中；采用最近点投影方法将固体模拟中的信息传递到流体网格面上。

2.数据交换。流体仿真和结构仿真需要建立通信，为了实现流固耦合的高效求解，需要超算集群多个计算节点之间的高效数据交换和协同计算，同时提高并行计算的效率和吞吐量，为此我们专门开发了一种优化后的分布式数据通信策略。

定制化界面开发

TankSimulator APP界面部分基于我们开发的定制应用开发框架快速搭建。对于框架本身，在前几期文章中已经有过详细介绍，这里就不展开介绍了。

定制应用开发案例之结构形状优化APP

前后处理平台介绍及二次开发案例

TankSimulator APP包括从几何模型导入的自动化前处理、仿真参数设置以及自动化后处理及报告生成。其中，根据实际工程和业务需要，我们设计了仿真参数输入界面，极大地简化了相关工程师进行繁杂的参数配置工作，降低了仿真的门槛。除此之外，自动化报告生成也使得仿真报告流程化和规范化。

（三）TankSimulator液体灌箱仿真APP 操作演示

1、登陆网址SimForge高性能仿真云进行账号注册

2、登陆平台后,找到对应图形应用图标(TankSimulator)，双击即可打开使用

神工坊2.0页面

下面是APP的界面展示，具体操作手册请登陆SimForge平台试用并查看。

TankSimulator APP界面展示

前处理及边界设置

一、固体/流体/仿真设置

Tanksimulator已为用户大幅简化了仿真流程，用户直接输入如下的参数即可

计算监控

在控制区选择队列和CPU核数。接下来需先点击计算准备按钮开始流体网格的生成并开启流体计算，如果是第一次生成流体网格,需要较长时间。观察到如下图的提示，说明已经准备就绪。

当计算准备过程完成后就可以点击开始计算开始流固双向耦合计算,计算过程会显示在中间及右边的文本框内,供观察计算进度。

Paraview后处理

计算完成后在计算监控面板中点击生成后处理结果。之后进入Paraview后处理标签页，相关页面已经准备好，用户可进行流体和固体的后处理。

自动生成报告

计算完成后，在计算监控面板标签页点击生成后处理结果，随后点击生成.docx报告，Tanksimulator将在后台自动提取分析数据，并绘图，自动生成完整的报告.docx文档。如下图所示，生成的报告可以从云桌面的文件系统中下载。

（四）结语

除TankSimulator液体灌箱仿真APP外，依托「SimForge HSF高性能数值模拟引擎」等技术，「神工坊」为不同行业、拥有特定仿真需求的客户快速定制复杂仿真应用，帮助客户实现设计研发、性能测试等过程的优化，从而降低成本、缩短周期、提高产品质量。

“仿真驱动创新，计算引领未来”，「神工坊」将始终践行“创新驱动、追求卓越、价值导向、开放合作”的价值观，继续保持对技术研发与创新的高度热情，与更多企业合作，探索更广泛、更深入的高性能数值模拟应用场景。

参考资料

[1]COMSOL.多物理场仿真百科[EB/OL].https://cn.comsol.com/multiphysics

SimForge 性能测试｜Fluent旋转机械稳态分析

前言CFD是工业仿真领域重要的分支之一，也是高性能计算的主要应用场景之一。本期选取了CFD领域的典型场景，稳态仿真计算案例——基于MRF方法的旋转机械流场分析，我们选用的软件是CFD领域最常用的仿真软件Fluent。我们来看下基于“神工坊”高性能工业仿真平台”的CFD稳态计算，和其他仿真云平台效率对比的情况。模拟与网格我们采用某品牌空调室外机作为稳态分析的仿真模型，如下图所示，左侧与后侧的进口流域，以及前侧的出口流域都考虑到计算中，并对空调内部结构简化后进行网格划分，最终网格单元数868万，其中，风扇叶片的旋转速度是850rpm。求解设置根据该款旋转机械的相关参数，经过理论计算得到该旋转机械的最大速度为25.6m/s，折合马赫数为0.075，为不可压缩流动，故选择压力基求解器，湍流模型选用了适用于旋转机械的k-εRealizable模型。对于动区域计算模型，本次稳态计算选择了网格静止不动的MRF旋转坐标系法，计算迭代步数400步，相关设置如下。仿真结果迭代完成之后仿真云图如下所示：仿真平台对比我们进行Fluent旋转机械稳态分析时,“神工坊”高性能工业仿真平台与其他两家仿真云平台的硬件参数如下表所示：计算过程中三个平台的一些输出日志如下图所示：本次仿真并行规模分别选取了16核、32核、64核、128核（受限于另外两个平台无法进行跨节点并行，并行规模无法进一步扩大），我们在“神工坊”平台进行了256核等更大规模的并行计算，结果显示计算用时会进一步缩短。“神工坊”高性能工业仿真平台与其他几家仿真云平台的计算时间如下图所示，其中，由于仿真云平台2最高只能64核并行使用，故图表中无仿真云平台2并行规模为128核的结果。可以发现，“神工坊”高性能工业仿真平台在进行稳态仿真分析时，表现出了绝对的速度优势。从16核到128核，其仿真计算速度都明显优于其他仿真云平台，且在相同并行核数下其仿真计算时间仅为其他仿真云平台的1/2不到，尤其是在64核并行时，其仿真计算时间更是只有仿真云平台1的1/3左右。我们以每个仿真云平台16核的计算时间为基本单位，计算各个平台的并行效率，结果如下图所示。我们可以发现“神工坊”高性能工业仿真平台的并行效率也是优于其他仿真云平台，且在每个核数下都保持着较高的并行效率。结论综上所述，SimForge高性能工业仿真云平台在进行Fluent稳态仿真分析时，无论是计算时间还是并行效率，均优于其他仿真云平台。（本文作者：郯俊建）“神工坊”SimForge高性能工业仿真平台以超算HPC集群作为硬件支撑，实现了跨节点大规模并行计算，对ANSYSFluent、Abaqus等仿真软件进行了CPU平台的高性能适配与优化，可以高效处理大规模网格模型以及复杂流场，大大缩短了企业仿真用时，提升工业设计效率，同时根据用户需求进行兼容性适配，力保每一核都能发挥出它的最大价值，欢迎广大工程师同行注册试用——1.专业GPU并行渲染，操作丝滑如本地仿真SimForge高性能仿真云平台实现了多GPU的分时共享，同时支持大规模仿真数据的多GPU服务端并行渲染，让10亿+网格可视化无压力。SimForge平台Fluent性能测试结果（注：仿真云平台2最高只能64核并行使用）2."超算级资源池+工业级软件栈"的垂直整合架构SimForge拥有亚洲最大的价值2亿的商业仿真软件授权，搭建了“传统商软+开源/国产软件+自研定制软件”的“工业级软件栈”，配合超算资源支持，单体软件并行最高可达2048核。3.打破工程仿真信息孤岛，在线协同范式升级实时同步操作界面：区别于大部分仿真云平台只有子母账号管理功能，SimForge在账号管理功能之外，是可以实现实时同步操作的。跨设备与跨账号无缝协同：关注到不同用户对账号信息保密需求不同，SimForge支持——“相同账号不同设备，同时登录协同操作”，“不同账号同时进入作业协同操作”两种使用情景。4.实时计费，精准结算，风险可控SimForge采用实时计费，精准结算的模式，让用户在使用过程中能够清晰了解费用情况，风险可控。无论是短期的集中计算任务，还是长期的研发项目，都能根据实际使用情况灵活计费，避免了资源浪费和不必要的成本支出，以及日常软硬件维护的烦恼。