黏性各向异性磁流体Kelvin-Helmholtz不稳定性: 二维数值研究

刘迎; 陈志华; 郑纯

doi:10.7498/aps.68.20181747

摘要

利用corner transport upwind和constrained transport算法求解非理想磁流体动力学方程组, 对匀强平行磁场作用下, 黏性各向异性等离子体自由剪切层中的Kelvin-Helmholtz不稳定性进行了数值模拟. 从流动结构、涡结构演化、磁场分布、横向磁压力、抗弯磁张力等角度对各向同性和各向异性黏性算例结果进行了讨论, 分析了黏性各向异性对Kelvin-Helmholtz不稳定性的影响. 结果表明, 黏性各向异性比黏性各向同性更利于流动的稳定. 其稳定性作用是由于磁感线方向上剪切速率降低导致界面卷起程度和圈数的降低, 并使卷起结构中小涡产生增殖、合并, 破坏了涡的常规增长, 从而导致流动的稳定. 黏性各向异性对横向磁压力的影响比对抗弯磁张力更大.

关键词:

Abstract

Kelvin-Helmholtz instability in anisotropic viscous fluid with uniform density in the presence of magnetic field is simulated through solving the non-ideal magnetohydrodynamic equations. The magnetic field is uniform and parallel to the stream. The magnetohydrodynamic equations are solved by corner transport upwind algorithm and constrained transport algorithm. In this paper, the influence of viscous anisotropy on Kelvin-Helmholtz instability is studied. The viscous anisotropy is embodied in the direction of the magnetic field, that is, viscosity parallel to the direction of the magnetic field line is much larger than that in the other directions. The results in the isotropic and the anisotropic viscous cases are compared from the aspects of flow structure, vortex evolution, and magnetic field distribution. It shows that the viscous anisotropy is more advantageous to the stability in a magnetized shear layer than the viscous isotropy. The flow structure evolves similarly in large scales but vortices evolve differently in small scales. Due to the decrease of the shear rate in the direction of the magnetic field lines, the rolling-up degree of interface and the number of laps decrease; the multiplication and merging of small vortices in the rolled-up structure destroy the regular growth of the vortex, which contributes to the stability of the flow. The increase of the magnetic energy at the sheared interface slows down effectively by the viscous anisotropy, which weakens the growth of the transverse magnetic pressure and anti-bending magnetic tension. However, viscous anisotropy shows much greater influence on the transverse magnetic pressure than on the anti-bending magnetic tension. The total enstrophy decreases slowly in viscous isotropy and anisotropy case. It increases quickly in late time in the former case, but is heavily suppressed in the latter case.

Keywords:

作者及机构信息

南京理工大学, 瞬态物理国家重点实验室, 南京　210094

通信作者: 陈志华, chenzh@mail.njust.edu.cn

基金项目: 国家自然科学基金(批准号: 11502117, 11272156)资助的课题.

Authors and contacts

Key Laboratory of Transient Physics, Nanjing University of Science and Technology, Nanjing 210094, China

Corresponding author: Chen Zhi-Hua, chenzh@mail.njust.edu.cn

Funds: Project supported by the National Natural Science Foundation of China (Grant Nos. 11502117, 11272156).

文章全文 : translate this paragraph

参考文献

[1]	Michikoshi S, Inutsuka, S I 2006 Astrophys. J. 641 1131 Google Scholar
[2]	Latter H N, Papaloizou J 2018 Mon. Not. R. Astron. Soc. 474 3110 Google Scholar
[3]	Masson A, Nykyri K 2018 Space Sci. Rev. 214 71 Google Scholar
[4]	Faganello M, Califano F 2017 J. Plasma Phys. 83 535830601 Google Scholar
[5]	Ma X Y, Otto A, Delamere P A 2014 J. Geophys. Res.: Space 119 781 Google Scholar
[6]	Srinivasan B, Tang X Z 2013 Phys. Plasmas 20 056307 Google Scholar
[7]	Wang L F, Ye W H, He X T, et al 2017 Sci. China: Phys. Mech. Astron. 60 055201 Google Scholar
[8]	Zhou Y 2017 Phys. Rep. 720-722 1 Google Scholar
[9]	Rao P A, Fuelberg H E 2000 Mon. Weather Rev. 128 3437 Google Scholar
[10]	Soloviev A V, Lukas R, Donelan M A, Haus B K, Ginis I 2017 J. Geophys. Res.: Oceans 122 10174 Google Scholar
[11]	Chandrasekhar S 1961 Hydrodynamic and Hydromagnetic Stability (Oxford: Clarendon Press)
[12]	Ho C M, Huerre P 1984 Annu. Rev. Fluid Mech. 16 365 Google Scholar
[13]	Winant C D, Browand F K 1974 J. Fluid Mech. 63 237 Google Scholar
[14]	Brown G L, Roshko A 1986 J. Fluid Mech. 170 499 Google Scholar
[15]	王立锋, 叶文华, 李英骏 2008 57 3038 Google Scholar Wang L F, Ye W H, Li Y J 2008 Acta Phys. Sin. 57 3038 Google Scholar
[16]	王立锋, 叶文华, 范征锋, 孙彦乾, 郑炳松, 李英骏 2009 58 6381 Google Scholar Wang L F, Ye W H, Fan Z F, Sun Y Q, Zheng B S, Li Y J 2009 Acta Phys. Sin. 58 6381 Google Scholar
[17]	Talwar S P 1965 Phys. Fluids 8 1295 Google Scholar
[18]	Duhau S, Gratton J 1975 J. Plasma Phys. 13 451 Google Scholar
[19]	Srivastava K M, Vyas D N 1979 Astrophys. Space Sci. 62 353 Google Scholar
[20]	Choudhury S R, Patel V L 1985 Phys. Fluids 28 3292 Google Scholar
[21]	Choudhury S R 1986 Phys. Fluids 29 1509 Google Scholar
[22]	Ruderman M S, Verwichte E, Erdélyi R, Goossens M 1996 J. Plasma Phys. 56 285 Google Scholar
[23]	Brown K G, Choudhury S R 2002 Q. Appl. Math. 60 601 Google Scholar
[24]	Prajapati R P, Chhajlani R K 2010 Phys. Plasmas 17 112108 Google Scholar
[25]	Gardiner T A, Stone J M 2008 J. Comput. Phys. 227 4123 Google Scholar
[26]	Liu C Q, Wang Y Q, Yang Y, Duan Z W 2016 Sci. China: Phys. Mech. Astron. 59 684711 Google Scholar
[27]	Liu Y, Chen Z H, Zhang H H, Lin Z Y 2018 Phys. Fluids 30 044102 Google Scholar

施引文献

图 1 计算模型

Fig. 1. Sketch of the computational model.

下载: 全尺寸图片幻灯片

图 2 不同时刻自由剪切层示踪剂浓度场　(a)各向同性黏性(Re = 10⁵); (b) 各向异性黏性(Re₀ = 10⁵, Re_|| = 100)

Fig. 2. Tracer concentration field at different times: (a) Isotropic viscous case (Re = 10⁵); (b) anisotropic viscous case (Re₀ = 10⁵, Re_|| = 100).

下载: 全尺寸图片幻灯片

图 3 不同时刻, 各向同性黏性(Re = 10⁵)算例涡结构云图及流线　(a) t = 5.5; (b) t = 7.0; (c) t = 7.5; (d) t = 8.0; (e) t = 9.0; (f) t = 10

Fig. 3. Rortex field with streamlines at different times in isotropic viscous fluid (Re = 10⁵): (a) t = 5.5; (b) t = 7.0; (c) t = 7.5; (d) t = 8.0; (e) t = 9.0; (f) t = 10.

下载: 全尺寸图片幻灯片

图 4 不同时刻各向异性黏性(Re₀ = 10⁵, Re_|| = 100)算例涡结构云图及流线　(a) t = 5.5; (b) t = 7.0; (c) t = 7.5; (d) t = 8.0; (e) t = 9.0; (f) t = 10

Fig. 4. Rortex field with streamlines at different times in anisotropic viscous fluid (Re₀ =10⁵, Re_|| = 100): (a) t = 5.5; (b) t = 7.0; (c) t = 7.5; (d) t = 8.0; (e) t = 9.0; (f) t = 10.

下载: 全尺寸图片幻灯片

图 5 不同时刻剪切层中磁场云图和磁感线　(a)各向同性黏性(Re = 10⁵); (b) 各向异性黏性(Re₀ = 10⁵, Re_|| = 100)

Fig. 5. Magnetic field with field lines at different times: (a) Isotropic viscous case (Re = 10⁵); (b) anisotropic viscous case (Re₀ = 10⁵, Re_|| = 100)

下载: 全尺寸图片幻灯片

图 6 纵向总动能随时间的变化

Fig. 6. Evolution of the longitudinal total kinetic energy E_ky .

下载: 全尺寸图片幻灯片

图 7 平均磁能密度放大倍数随时间的变化

Fig. 7. Evolution of the amplification factor of average magnetic energy.

下载: 全尺寸图片幻灯片

图 8 平均横向磁压力和抗弯磁张力随时间的变化　(a) 横向磁压力; (b) 抗弯磁张力

Fig. 8. Evolution of the average transverse magnetic pressure and anti-bending magnetic tension: (a) Transverse magnetic pressure; (b) anti-bending magnetic tension.

下载: 全尺寸图片幻灯片

图 9 流场总涡度拟能随时间的变化

Fig. 9. Evolution of the total enstrophy.

下载: 全尺寸图片幻灯片

map

[1]	Michikoshi S, Inutsuka, S I 2006 Astrophys. J. 641 1131 Google Scholar
[2]	Latter H N, Papaloizou J 2018 Mon. Not. R. Astron. Soc. 474 3110 Google Scholar
[3]	Masson A, Nykyri K 2018 Space Sci. Rev. 214 71 Google Scholar
[4]	Faganello M, Califano F 2017 J. Plasma Phys. 83 535830601 Google Scholar
[5]	Ma X Y, Otto A, Delamere P A 2014 J. Geophys. Res.: Space 119 781 Google Scholar
[6]	Srinivasan B, Tang X Z 2013 Phys. Plasmas 20 056307 Google Scholar
[7]	Wang L F, Ye W H, He X T, et al 2017 Sci. China: Phys. Mech. Astron. 60 055201 Google Scholar
[8]	Zhou Y 2017 Phys. Rep. 720-722 1 Google Scholar
[9]	Rao P A, Fuelberg H E 2000 Mon. Weather Rev. 128 3437 Google Scholar
[10]	Soloviev A V, Lukas R, Donelan M A, Haus B K, Ginis I 2017 J. Geophys. Res.: Oceans 122 10174 Google Scholar
[11]	Chandrasekhar S 1961 Hydrodynamic and Hydromagnetic Stability (Oxford: Clarendon Press)
[12]	Ho C M, Huerre P 1984 Annu. Rev. Fluid Mech. 16 365 Google Scholar
[13]	Winant C D, Browand F K 1974 J. Fluid Mech. 63 237 Google Scholar
[14]	Brown G L, Roshko A 1986 J. Fluid Mech. 170 499 Google Scholar
[15]	王立锋, 叶文华, 李英骏 2008 57 3038 Google Scholar Wang L F, Ye W H, Li Y J 2008 Acta Phys. Sin. 57 3038 Google Scholar
[16]	王立锋, 叶文华, 范征锋, 孙彦乾, 郑炳松, 李英骏 2009 58 6381 Google Scholar Wang L F, Ye W H, Fan Z F, Sun Y Q, Zheng B S, Li Y J 2009 Acta Phys. Sin. 58 6381 Google Scholar
[17]	Talwar S P 1965 Phys. Fluids 8 1295 Google Scholar
[18]	Duhau S, Gratton J 1975 J. Plasma Phys. 13 451 Google Scholar
[19]	Srivastava K M, Vyas D N 1979 Astrophys. Space Sci. 62 353 Google Scholar
[20]	Choudhury S R, Patel V L 1985 Phys. Fluids 28 3292 Google Scholar
[21]	Choudhury S R 1986 Phys. Fluids 29 1509 Google Scholar
[22]	Ruderman M S, Verwichte E, Erdélyi R, Goossens M 1996 J. Plasma Phys. 56 285 Google Scholar
[23]	Brown K G, Choudhury S R 2002 Q. Appl. Math. 60 601 Google Scholar
[24]	Prajapati R P, Chhajlani R K 2010 Phys. Plasmas 17 112108 Google Scholar
[25]	Gardiner T A, Stone J M 2008 J. Comput. Phys. 227 4123 Google Scholar
[26]	Liu C Q, Wang Y Q, Yang Y, Duan Z W 2016 Sci. China: Phys. Mech. Astron. 59 684711 Google Scholar
[27]	Liu Y, Chen Z H, Zhang H H, Lin Z Y 2018 Phys. Fluids 30 044102 Google Scholar

[1]	孙思杰, 蒋晗. 各向异性界面动力学对深胞晶生长形态稳定性的影响. , 2024, 73(11): 118101. doi: 10.7498/aps.73.20240362
[2]	张升博, 张焕好, 张军, 毛勇建, 陈志华, 石启陈, 郑纯. 激波与轻质气柱作用过程的磁场抑制特性. , 2024, 73(8): 084701. doi: 10.7498/aps.73.20231916
[3]	沈勇, 董家齐, 何宏达, 潘卫, 郝广周. 托卡马克理想导体壁与磁流体不稳定性. , 2023, 72(3): 035203. doi: 10.7498/aps.72.20222043
[4]	罗仕超, 吴里银, 常雨. 高超声速湍流流动磁流体动力学控制机理. , 2022, 71(21): 214702. doi: 10.7498/aps.71.20220941
[5]	范海龙, 陈明文. 磁场对二元合金凝固过程中糊状层稳定性的影响. , 2021, 70(6): 066401. doi: 10.7498/aps.70.20201748
[6]	丁明松, 江涛, 刘庆宗, 董维中, 高铁锁, 傅杨奥骁. 基于电流积分计算磁矢量势修正的低磁雷诺数方法. , 2020, 69(13): 134702. doi: 10.7498/aps.69.20200091
[7]	丁明松, 傅杨奥骁, 高铁锁, 董维中, 江涛, 刘庆宗. 高超声速磁流体力学控制霍尔效应影响. , 2020, 69(21): 214703. doi: 10.7498/aps.69.20200630
[8]	丁明松, 江涛, 董维中, 高铁锁, 刘庆宗, 傅杨奥骁. 热化学模型对高超声速磁流体控制数值模拟影响分析. , 2019, 68(17): 174702. doi: 10.7498/aps.68.20190378
[9]	董国丹, 郭则庆, 秦建华, 张焕好, 姜孝海, 陈志华, 沙莎. 不同磁场构型下Richtmyer-Meshkov不稳定性的数值研究及动态模态分解. , 2019, 68(16): 165201. doi: 10.7498/aps.68.20190410
[10]	董国丹, 张焕好, 林震亚, 秦建华, 陈志华, 郭则庆, 沙莎. 磁控条件下激波冲击三角形气柱过程的数值研究. , 2018, 67(20): 204701. doi: 10.7498/aps.67.20181127
[11]	王鹏, 薛纭, 楼智美. 黏性流体中超细长弹性杆的动力学不稳定性. , 2017, 66(9): 094501. doi: 10.7498/aps.66.094501
[12]	薛创, 丁宁, 孙顺凯, 肖德龙, 张扬, 黄俊, 宁成, 束小建. 脉冲功率驱动器与Z箍缩负载耦合的全电路数值模拟. , 2014, 63(12): 125207. doi: 10.7498/aps.63.125207
[13]	郎莎莎, 耿兴国, 臧渡洋. 八重准周期排列的短沟槽结构减阻机理分析. , 2014, 63(8): 084704. doi: 10.7498/aps.63.084704
[14]	袁永腾, 郝轶聃, 侯立飞, 涂绍勇, 邓博, 胡昕, 易荣清, 曹柱荣, 江少恩, 刘慎业, 丁永坤, 缪文勇. 流体力学不稳定性增长测量方法研究. , 2012, 61(11): 115203. doi: 10.7498/aps.61.115203
[15]	杨涓, 石峰, 杨铁链, 孟志强. 电子回旋共振离子推力器放电室等离子体数值模拟. , 2010, 59(12): 8701-8706. doi: 10.7498/aps.59.8701
[16]	吴翊, 荣命哲, 杨飞, 王小华, 马强, 王伟宗. 引入6波段P-1辐射模型的三维空气电弧等离子体数值分析. , 2008, 57(9): 5761-5767. doi: 10.7498/aps.57.5761
[17]	池长昀, 张金仓, 李领伟, 刘芬, 黎文峰, 敬超, 曹世勋, Miryala Muralidhar, 姚忻. 熔融织构(Nd-Eu-Gd)-Ba-Cu-O超导体的各向异性磁通跳跃及其不稳定性研究. , 2005, 54(5): 2307-2312. doi: 10.7498/aps.54.2307
[18]	李文飞, 张丰收. 非对称核物质的化学不稳定性与力学不稳定性. , 2001, 50(10): 1888-1895. doi: 10.7498/aps.50.1888
[19]	匡光力, G.WAIDMANN. TEXTOR托卡马克等离子体的磁流体动力学振荡特性. , 1994, 43(9): 1466-1475. doi: 10.7498/aps.43.1466
[20]	贾惟义, 张鹏翔. 磁晶各向异性场引起的YIG单晶微波器件温度不稳定性的最佳补偿. , 1976, 25(3): 254-264. doi: 10.7498/aps.25.254

计量

文章访问数: 8567
PDF下载量: 102
被引次数: 0

姓名
邮箱
手机号码
标题
留言内容
验证码

搜索

留言板

黏性各向异性磁流体Kelvin-Helmholtz不稳定性: 二维数值研究