[论文解读] EOF Library: Open-Source Elmer and OpenFOAM Coupler for Simulation of MHD With Free Surface

Juris Venčels, Andris Jakovičs+2 more|arXiv (Cornell University)|Jul 13, 2017

Magnetic confinement fusion research参考文献 2被引用 3

一句话总结

本文提出 EOF Library，这是一个开源的双向耦合框架，用于在 Elmer（有限元法）与 OpenFOAM（计算流体动力学）之间模拟具有自由表面的磁流体动力学（MHD）。该框架实现了求解器之间高效的数据插值与交换，通过一个二维轴对称 MHD 问题（含导电流体与电磁感应）进行了验证，结果与参考数据高度一致。

ABSTRACT

An open-source library called EOF which provides efficient two-way coupling between the finite element multi-physics simulation software Elmer and computational fluid dynamics software OpenFOAM is presented. Library controls interpolation and communication of fields between two codes. The main focus of this code is industrial magnetohydrodynamics with free surface, but it can be extended to other 2D and 3D multi-physics problems. Simulation results for 2D axisymmetric problem of conductive fluid with free surface surrounded by alternate electromagnetic field created by inductor are compared with analogous simulation results made by another author.

研究动机与目标

• 解决工业磁流体动力学（MHD）中自由表面模拟的精确性与效率需求，这是一个复杂的多物理场问题。
• 通过实现 Elmer 与 OpenFOAM 之间的无缝数据交换，克服独立求解器的局限性。
• 开发一个可重用的开源库，促进有限元与有限体积求解器在多物理场问题中的耦合。
• 实现二维与三维 MHD 问题的模拟，包含自由表面动力学及电磁场作用。
• 为金属铸造与液态金属加工等工业应用提供一个经过验证的框架。

提出的方法

• 使用共享内存与基于文件的通信方式，在 Elmer 与 OpenFOAM 之间实现双向数据耦合机制。
• 开发插值例程，实现在不同求解器的非结构化网格之间传递场量（如速度、磁场、压力）。
• 采用分区网格处理技术，支持并行计算并保持可扩展性。
• 将库作为插件集成到 Elmer 与 OpenFOAM 中，实现模块化与可扩展的耦合。
• 在 Elmer 中采用交错有限元法，在 OpenFOAM 中采用有限体积法，确保界面处耦合的一致性。
• 通过同步时间步长与使用保持守恒性与精度的插值方案，确保时间一致性。

实验结果

研究问题

• RQ1 能否在 Elmer 与 OpenFOAM 之间开发出一种高效、开源的双向耦合库，用于多物理场 MHD 模拟？
• RQ2 EOF Library 在电磁感应作用下对具有自由表面的 MHD 问题，其结果与已知参考结果的吻合度如何？
• RQ3 该耦合框架在二维轴对称 MHD 配置下的性能与稳定性如何？
• RQ4 该框架在 MHD 之外的其他二维与三维多物理场问题中可扩展到何种程度？
• RQ5 场插值方法对耦合系统中守恒性与收敛性的影响如何？

主要发现

• EOF Library 有效实现了 Elmer 与 OpenFOAM 之间的双向数据交换，支持复杂的多物理场模拟。
• 针对具有导电流体与电磁感应的二维轴对称 MHD 问题，模拟结果与另一项研究的参考数据高度一致。
• 通过一致的插值与时间同步，耦合框架保持了数值稳定性与精度。
• 该库具有可扩展性，可适用于 MHD 以外的其他多物理场问题。
• 库的开源特性促进了计算物理领域的结果可复现性与社区驱动的开发。
• 该实现证明了在涉及自由表面与电磁场的工业规模 MHD 模拟中具有可行性。