基于OpenFOAM的投影法磁流体求解器开发与验证

doi:10.16568/j.0254-6086.201801003

核聚变与等离子体物理 ›› 2018, Vol. 2 ›› Issue (1): 15-20.DOI: 10.16568/j.0254-6086.201801003

基于OpenFOAM的投影法磁流体求解器开发与验证

毛洁1，王浩1，刘克2，王盛1，Claude B Reed3

(1. 杭州电子科技大学机械工程学院，杭州 621000； 2. 大连理工大学能源和动力工程学院，大连 116024； 3. Nuclear Engineering Division, Argonne National Laboratory, Argonne, IL, USA, 60439)

收稿日期:2017-06-02 修回日期:2017-10-18 出版日期:2018-03-15 发布日期:2018-03-15
作者简介:毛洁(1974-)，女，河北巨鹿人，博士，教授，主要从事磁流体动力学、计算流体力学数值模拟方法的研究。
基金资助:
国家自然科学基金(11375049)；国家磁约束核聚变能发展研究专项(2014GB125003)

Development and validation of magnetohydrodynamic solver with projection in OpenFOAM environment

MAO Jie, WANG Hao, LIU Ke, WANG Cheng, Claude B Reed

(1. School of Mechanical Engineering, Hangzhou Dianzi University, Hangzhou 310018; 2. School of Energy & Power Engineering, Dalian University of Technology, Dalian 116024; 3. Nuclear Engineering Division, Argonne National Laboratory, Argonne, IL, USA, 60439)

Received:2017-06-02 Revised:2017-10-18 Online:2018-03-15 Published:2018-03-15

摘要/Abstract

摘要：

在开源计算流体力学C++工具包OpenFOAM环境下开发了低磁雷诺数条件下的磁流体求解器，并进行了验证。采用投影算法求解动量方程和压力泊松方程；采用非结构网格同位相容守恒算法求解电势泊松方程、感应电流和洛伦兹力；采用边界耦合方法求解流固耦合电势场。通过对均匀磁场下导电方管和导电圆管内的完全发展磁流体层流的数值模拟和解析解的对比，对求解器进行了验证。进一步对非均匀强磁场作用下导电方管和导电圆管内完全发展磁流体层流进行了数值模拟，并与ALEX实验结果进行了比较。数值解和实验结果吻合良好。所开发的求解器可用于复杂结构强磁场作用下磁流体的数值模拟研究。

关键词: 磁流体, OpenFOAM, 投影算法, 相容守恒, 边界耦合, 验证

Abstract:

Magnetohydrodynamic flow solver with low magnetic Reynolds number in open source CFD environment of OpenFOAM was developed and validated. Four step projection methods were used to solve the momentum equation and the mass conservative equation. The consistent and conservative schemes based on collocated unstructured grid were applied to solve the electric potential Poisson equation, the induced electric current and the Lorentz force. A multi-region technique coupled by the boundary was also developed to solve magnetohydrodynamic flow with electric conducting walls. The solver was verified and validated by simulating duct and pipe flow with conducting walls under uniform magnetic field. Furthermore, the multi-region MHD solver was also validated by simulating liquid metal flows subject to a non-uniform magnetic field in an electrically conducting rectangular duct and pipe. The pressure distribution was compared with the ALEX experimental data. The results show that the solver has a good accuracy and can be used in the numerical study of MHD flows in fusion environment.

Key words: Magnetohydrodynamic, OpenFOAM, Projection method, Consistent and conservative scheme, Boundary couple, Validation

中图分类号:

TL61+2

毛洁1，王浩1，刘克2，王盛1，Claude B Reed3. 基于OpenFOAM的投影法磁流体求解器开发与验证[J]. 核聚变与等离子体物理, 2018, 2(1): 15-20.

MAO Jie, WANG Hao, LIU Ke, WANG Cheng, Claude B Reed. Development and validation of magnetohydrodynamic solver with projection in OpenFOAM environment[J]. NUCLEAR FUSION AND PLASMA PHYSICS, 2018, 2(1): 15-20.

参考文献

[1] Moreau R. Magnetohydrodynamics[M]. Kluwer Academic Publishers, 1990.

[2] Ni M J, Munipalli R, Morley N B, et al. A current density conservative scheme for incompressible MHD flows at a low magnetic Reynolds number, part I: on a rectangular collocated grid system[J]. J. Comp. Phys., 2007, 227(1): 174-204.

[3] Ni M J, Munipalli R, Huang P, et al. A current density conservative scheme for incompressible MHD flows at a low magnetic Reynolds number, part II: on an arbitrary collocated mesh[J]. J. Comp. Phys., 2007, 227(1): 205-228.

[4] Ni M J, Li J F. A consistent and conservative scheme for incompressible MHD flows at a low magnetic Reynolds number, part III: on a staggered mesh[J]. J. Comp. Phys., 2012, 231(2): 281-298.

[5] Chang N K. Computational analysis of a liquid metal magnetohydrodynamic flow in a manifold under a uniform magnetic field[J]. Fusion Science & Technology, 2013, 64(4): 787-799.

[6] Smolentsev S, Morley N, Abdou M. Code development for analysis of MHD pressure drop reduction in a liquid metal blanket using insulation technique based on a fully developed flow model[J]. Fusion Engineering & Design, 2005, 73(1): 83-93.

[7] Samad S A. The flow of conducting fluids through circular pipes having finite conductivity and finite thickness under uniform transverse magnetic fields
[J]. International Journal of Engineering Science, 1981, 19(9): 1221-1232.

[8] Hunt J C R. Magnetohydrodynamic flow in rectangular ducts[J]. Journal of Fluid Mechanics, 1965, 23(4): 563-581.

[9] Reed C B, Picologlou B F, Hua T Q, et al. ALEX results: A comparison of measurements from a round and a rectangular duct with 3-D code predictions[C]. Symposium on Fusion Engineering. Presented at the 12th Symposium on Fusion Engineering, Monterey, Calif. 12 Oct. 1987, 1987.

[10] Picologlou B F, Reed C B. Experimental investigation of 3-D MHD flows at high hartmann number and interaction parameter[M]. Liquid Metal Magnetohydrodynamics. Springer Netherlands, 1988. 239-241.

[11] Smolentsev S, Badia S, Bhattacharyay R, et al. An approach to verification and validation of MHD codes for fusion applications[J]. Fusion Engineering & Design, 2014, 100: 65-72.

[12] Albetschico X, Radhakrishnan H, Kassinos S, et al. Numerical simulation of a liquid-metal flow in a poorly conducting pipe subjected to a strong fringing magnetic field[J]. Physics of Fluids, 2011, 23(4): 171-194.

[13] Feng J, Chen H, He Q, et al. Further validation of liquid metal MHD code for unstructured grid based on OpenFOAM[J]. Fusion Engineering & Design, 2015, 100: 260-264.

[1]	董铠玮, 王俊, 王晓宇, 武兴华, 郭礼凯. 横向倾斜磁场作用下铁磁低活化钢矩形管道内磁流体流动的分析 [J]. 核工业西南物理研究院, 2025, 45(1): 98-103.
[2]	高正坤, 杨文军, 龚学余, 钟翊君. 磁场位形对内交换模不稳定性影响的研究[J]. 核工业西南物理研究院, 2025, 45(1): 104-110.
[3]	杨秋磊, 陈伟, 张洁, 程时葵, 梁绍勇, 张轶泼, HL- 实验团队. HL-3 装置欧姆等离子体中高频磁流体不稳定性的实验观测[J]. 核工业西南物理研究院, 2024, 44(3): 328-333.
[4]	李漫, 王先驱, 苏祥, 许宇鸿, 程钧, 刘海峰, 黄捷, 张欣, 刘海, 唐昌建 . 局部反磁剪切位形下高能量粒子模的线性不稳定性分析[J]. 核工业西南物理研究院, 2023, 43(4): 429-435.
[5]	闫家玮, 毛洁, 熊保龙. 平行壁变形导电管磁流体层流数值模拟[J]. 核聚变与等离子体物理, 2022, 42(4): 448-454.
[6]	蒙建朋, 贾瑞宝, 蔡强, 邱丽媛, 颉延风 . HL-2M 装置烘烤系统管路应力的仿真分析与优化[J]. 核聚变与等离子体物理, 2021, 41(s1): 388-391.
[7]	谭耀, 马骏, 项农 . 基于微扰守恒模型的磁流体模拟并行计算 [J]. 核聚变与等离子体物理, 2021, 41(1): 8-15.
[8]	王磊, 张秀杰, 潘传杰, 许增裕, 赵耀, 王凡. 聚变堆液态包层准二维磁流体湍流模型研究[J]. 核聚变与等离子体物理, 2021, 41(1): 21-25.
[9]	王坤，姜林村，张董，张阅，王天浩. 真空中金属丝电爆炸数值模拟研究[J]. 核聚变与等离子体物理, 2020, 40(4): 315-320.
[10]	陈豪杰，刘海华*，李亮玉，赵淘 . 微束等离子体弧焊电弧物理特性多场耦合分析[J]. 核聚变与等离子体物理, 2020, 40(3): 283-288.
[11]	徐超，毛洁*，尹跃广，王浩. 带引流条导电矩形管磁流体湍流大涡数值模拟[J]. 核聚变与等离子体物理, 2019, 39(2): 112-118.
[12]	刘佰奇1，阳倦成2，齐天煜1，任东伟1，倪明玖1. 水平磁场中液态锂膜流动实验装置及初步结果[J]. 核聚变与等离子体物理, 2019, 39(1): 62-67.
[13]	李俊龙1, 2，朱平2，高鹏3，任安业3. 空气域中铁磁流体动态控制研究[J]. 核聚变与等离子体物理, 2018, 2(4): 402-407.
[14]	王俊，张国书，张龙，王晓宇. 方管锂铅磁流体氚增殖剂中的氚输运模拟[J]. 核聚变与等离子体物理, 2018, 2(3): 254-259.
[15]	曹小岗1，魏强林1，韩磊2，马小春3，张志艳2，韦建军2，芶富均3. 液态锂综合实验平台中液态锂流面与氩等离子体相互作用[J]. 核聚变与等离子体物理, 2018, 2(3): 260-268.

基于OpenFOAM的投影法磁流体求解器开发与验证

Development and validation of magnetohydrodynamic solver with projection in OpenFOAM environment

PDF

可视化

摘要/Abstract

引用本文

使用本文

参考文献

相关文章 15

编辑推荐

Metrics