磁场作用下液-气自由界面动态行为的两种数值处理方法.pdfVIP

磁场作用下液-气自由界面动态行为的 两种数值处理方法1 a,b,2 b b b a 雷作胜 , Yves FAUTRELLE , Roland ERNST , Jecqueline ETAY , 任忠鸣 a. 上海大学 （省部共建）现代冶金与材料制备重点实验室，上海大学，上海200072 b. Group EPM du laboratoir SIMAP, CNRS UMR 5266, 38402 St-MARTIN d ’HERES, France 摘要：本文结合两个简单例子，介绍了在磁场作用下液－气两相自由界面动态行为数值模拟方面开展的初步工作。第一个方 面，基于水平集方法 （Level Set Method ）建立了一个二维数学模型，计算了初始形状为圆形的单个气泡，在浮力作用下从 金属液中某一深度处释放后，上浮过程中气泡的变形、破碎等动态行为。垂直于这一二维平面施加一个静磁场，发现磁场对 金属液的制动效应减小了气泡的变形程度，降低了气泡的上浮速度，起到了控制气泡上浮行为的作用。第二个方面，基于 ALE （Arbitrary Lagrangian-Eulerian ）方法，研究了单个金属液滴在低频交变电磁场作用下，金属液自由界面的动态振动及 其稳定性行为。其物理过程为，在X-Y 平面内一个初始形状为圆形的金属液滴，当在Z 方向施加一个大小处处相等的低频 交变磁场时，在变化的电磁感应Lorentz 力和表面张力的共同作用下，金属液滴的自由界面变形、振荡，表现出十分复杂的 动态行为。为了处理界面形状变化时的动态表面张力，推导了基于 N-S 方程弱解形式的的自由界面边界条件。数值模拟结 果表明，在电磁场频率和磁感应强度这一参数空间内，金属液滴存在参数不稳定行为。在特定磁场频率下，当磁感应强度超 过某一临界值，在金属液滴中会激发不同模式的类多边形振动行为。通过数值模拟的结果给出了从模式2 到模式7 的稳定性 图，分析表明这一稳定性结构具有Mathieu 参数不稳定的典型特征。 关键词：磁场；金属液；气泡；自由界面；数值模拟 引言 在冶金和金属材料加工工艺中，很多过程都涉及到金属液-气体两相自由界面的动态行为，比较典型的 例子有：气泡在金属液中上浮时气泡的形态演化，金属悬浮熔炼、焊接过程金属熔池自由界面形貌，容器 内金属液自由表面的晃动，悬浮或者附着金属液滴的动态行为等等。通过施加形式各异的电磁场，来影响 或者控制液－气两相自由界面的行为，是近年来逐渐发展起来的一些技术手段，在这些工艺优化过程中都 涉及到一些磁流体力学和两相流问题。其中的核心问题是如何处理两相界面的动态演化，以及和界面形态 密切相关的曲率、表面张力等。 本文从两个简单的例子出发，介绍了在这方面开展的初步工作。第一个例子是基于水平集方法，数值 模拟了一个初始形状为圆形的单个气泡，从金属液中某一深度处释放后，上浮过程中气泡的变形、破碎等 动态行为，考查了静磁场对这一动态过程的影响。第二个例子是基于ALE （Arbitrary Lagrangian-Eulerian ） 方法，研究了在x-y 平面内的一个金属液滴在垂直于此平面的低频交变电磁场作用下，金属液自由界面的 动态振动及其稳定性行为。 1 国家自然科学基金和上海市科委(09QA1402200)资助项目 2 E_mail: lei_zsh@staff.shu.edu.cn 1 1 单个气泡上浮动态行为及静磁场对其的影响 1.1 问题描述 考虑如图 1 所示的物理过程，在x-y 平面内的金 属熔池某一深度处，释放一定大小的圆形气泡，在浮 y B ③ 力作用下气泡上浮并带动金属液流动，在流体的剪切 z gas