磁流体方腔槽道流整体线性稳定性数值方法研究

计算物理 ›› 2016, Vol. 33 ›› Issue (4): 379-390.

• 研究论文 • 下一篇

磁流体方腔槽道流整体线性稳定性数值方法研究

胡军¹, 刘婵², 张年梅², 倪明玖²

1. 北京应用物理与计算数学研究所, 北京 100094;
2. 中国科学院大学物理学院, 北京 100049

收稿日期:2015-06-03 修回日期:2015-10-08 出版日期:2016-07-25 发布日期:2016-07-25
作者简介:胡军(1978-),男,博士,副研究员,从事流体力学稳定性和数值模拟不确定度量化研究,E-mail:hu_jun@iapcm.ac.cn
基金资助:
国家自然科学基金重点项目（50936066，11125212）及国家自然科学基金面上项目（11172049，11472060）资助

Numerical Methods for Linear Global Stability of Magnetohydrodynamic Duct Flows

HU Jun¹, LIU Chan², ZHANG Nianmei², NI Mingjiu²

1. Institute of Applied Physics and Computational Mathematics, Beijing 100094, China;
2. University of Chinese Academy of Sciences, Beijing 100049, China

Received:2015-06-03 Revised:2015-10-08 Online:2016-07-25 Published:2016-07-25

摘要/Abstract

摘要： 将Chebyshev谱配置法和基于非均匀网格的高阶FD-q差分格式运用于磁流体方腔槽道流整体线性稳定性研究，比较两类数值方法的优缺点.Chebyshev谱配置法收敛快且精度高，但需要构造非常庞大的满矩阵，存储量和计算开销巨大；高阶FD-q差分格式采用了基于Kosloff-Tal-Ezer变换的Chebyshev谱配置点作为离散网格，在保持较高网格收敛精度的同时，差分格式可以采用稀疏矩阵进行存储，显著降低了存储量和计算开销.相比传统的谱配置法，基于非均匀网格的高阶FD-q差分格式计算效率得到显著的提升，将高阶FD-q差分格式运用于非正则模线性最优瞬态增长的计算，计算效果良好.

关键词: 磁流体方腔槽道流, 整体线性稳定性, 线性最优瞬态增长, Hunt槽道流

Abstract: Spectral Chebyshev collocation method and high-order FD-q finite difference method are used for global instability analysis of magetohydrodynamic(MHD) duct flows and compared for their merits and drawbacks. Spectral Chebyshev collocation method has faster convergence rate and high-order accuracy, while it needs to construct full general eigenvalue matrix which would consume large memory storage and a great deal of computational cost. High-order FD-q finite difference method adopts modified Chebyshev collocation points as discretization mesh grids based on Kosloff-Tal-Ezer transformation. FD-q method can maintain high convergence rate of mesh grids, and resulted general eigenvalue matrix is very sparse and can be stored with sparse matrix, which greatly reduces computational resource. In contrast to traditional spectral collocation method, non-uniform mesh based FD-q method obtains remarkable progress on computational efficiency, which is further demonstrated by computation of linear optimal transient growth for MHD duct flows.

Key words: magnetohydrodynamic duct flows, linear global stability, linear optimal transient growth, Hunt flows

中图分类号:

O361

胡军, 刘婵, 张年梅, 倪明玖. 磁流体方腔槽道流整体线性稳定性数值方法研究[J]. 计算物理, 2016, 33(4): 379-390.

HU Jun, LIU Chan, ZHANG Nianmei, NI Mingjiu. Numerical Methods for Linear Global Stability of Magnetohydrodynamic Duct Flows[J]. CHINESE JOURNAL OF COMPUTATIONAL PHYSICS, 2016, 33(4): 379-390.

参考文献

[1] DAVIDSON P A. Magnetohydrodynamics in materials processing[J]. Annu Rev Fluid Mech, 1999, 31:273-300.
[2] THOMAS B G, ZHANG L. Mathematical modeling of fluid flow in continuous casting:a review[J]. ISIJ Intl, 2001,41:1181-1193.
[3] BÜHLER L. Liquid metal magnetohydrodynamics for fusion blankets[M]//MOLOKOV S, MOREAU R, MOFFATT H K, ed. Magnetohydrodynamics:historical evolution and trends, Springer, 2007:171-194.
[4] 倪明玖.磁约束核聚变关键能量转换部件的磁流体力学研究[J].力学与实践,2010, 32(5):1-9.
[5] KAKUTANI T. The hydrodynamic stability of the modified plane Couette flow in the presence of a transverse magnetic field[J]. J Phys Soc Japan, 1964, 19:1041-1057.
[6] LEHNERT B. On the behaviour of an electrically conductive liquid in a magnetic field[J]. Ark Fys,1952, 5:69-90.
[7] HUNT J C R. Magnetohydrodynamic flow in rectangular ducts[J]. J Fluid Mech, 1965, 21:577-590.
[8] MOFFATT H K.On the suppression of turbulence by a uniform magnetic field[J]. J Fluid Mech, 1967, 28:571-592.
[9] DAVIDSON P A. Magnetic damping of jets and vortices[J]. J Fluid Mech, 1995, 299:153-186.
[10] SHATROV V, GERBETH G. Marginal turbulent magnetohydrodynamic flow in a square duct[J]. Physics of Fluids, 2010, 22:084101.
[11] KINET M, KNAEPEN B, MOLOKOV S. Instabilities and transition in magnetohydrodynamic flows in ducts with electrically conducting walls[J]. Physical Review Letters, 2009, 103(15):154501.
[12] PRIEDE J, ALEKSANDROVA S, MOLOKOV S. Linear stability of
Hunt's flow[J]. Journal of Fluid Mechanics, 2010, 649:115.
[13] PRIEDE J, ALEKSANDROVA, MOLOKOV S. Linear stability of magnetohydrodynamic flow in a perfectly conducting rectangular duct[J]. Journal of Fluid Mechanics, 2012, 708:111-127.
[14] KRASNOV D, ZIKANOV O, ROSSI M, et al. Optimal linear growth in magnetohydrodynamic duct flow[J]. Journal of Fluid Mechanics, 2010, 653:273-299.
[15] PAREDES P, HERMANNS M, CLAINCHE S L, THEOFILIS V. Order 10⁴ speedup in global linear instability analysis using matrix formation[J]. Comput Methods Appl Mech Engrg, 2013, 253:287-304.
[16] SCHMID P J, HENNINGSON D S. Stability and transition in shear flows[M]. New York:Springer-Verlag,2001,142.
[17] SCHMID P J. Nonmodal Stability Theory[J]. Ann Rev Fluid Mech, 2007, 39:129-162.
[18] TUCKERMAN L S, BARKLEY D. Bifurcation analysis for time-steppers[M]//DOEDEL E,TUCKERMAN L S, ed. Numerical methods for bifurcation problems and large-scale dynamical systems, New York:Springer, 2000.
[19] THEOFILIS V. Advances in global linear instability analysis of nonparallel and three-dimensional flows[J]. Prog Aero Sci, 2003, 39:249-315.
[20] HERMANNS M, HERNÀNDEZ J A. Stable high-order finite-difference methods based on non-uniform grid point distributions[J]. Int J Numer Meth Fluids, 2008, 56:233-255.
[21] KOSLOFF D, TAL-EZER H. A modified Chebyshev pseudospectral method with an O(N-1) time step restriction[J]. J Comput Phys, 1993, 104(2):457-469.
[22] THEOFILIS V. Global linear instability[J]. Annu Rev Fluid Mech, 2011, 43:319-352.
[23] DAVIS T A. Algorithm 832:UMFPACK, an unsymmetric-pattern multifrontal method[J]. ACM Transactions on Mathematical Software, 2004,30(2):196-199.
[24] AMESTORY P, DUFF I, KOSTER J. A fully asynchronous multifrontal solver using distributed dynamic scheduling[J]. SIAM J Matrix Anal, 2001, A 1:15-41.
[25] LEHOUCQ R B, SORENSEN D C, YANG C. ARPACK users guide:Solution of large-scale eigenvalue problems with implicitly restarted Arnoldi methods[M]. Philadelphia:SIAM,1998.
[26] BIAU D, SOUEIDA H, BPTTARO A. Transition to turbulence in duct flow[J]. J Fluid Mech, 2008, 596:133-142.
[27] BERRUT J P, TREFETHEN LN. Barycentric Lagrange interpolation[J]. SIAM Review, 2004, 46:501-517.
[28] HIGHAM N. The numerical stability of barycentric Lagrange interpolation[J]. IMA Journal of Numerical Analysis, 2004, 24:547-556.

磁流体方腔槽道流整体线性稳定性数值方法研究

Numerical Methods for Linear Global Stability of Magnetohydrodynamic Duct Flows

RichHTML

PDF

可视化

摘要/Abstract

引用本文

使用本文

参考文献

相关文章 10

编辑推荐

Metrics

作者中心

审稿中心

期刊浏览

期刊介绍

[1]	严春晖, 肖波, 王刚华, 陆禹, 李平. 一种基于磁通和电磁能流的磁扩散问题求解格式[J]. 计算物理, 0, (): 1-.
[2]	毛洁, 相凯, 王彦利, 王浩. 导电突扩管中磁流体管流的数值模拟[J]. 计算物理, 2018, 35(5): 597-605.
[3]	杨云. 磁流体数值模拟中一种消去磁场散度的方法[J]. 计算物理, 2018, 35(4): 437-442.
[4]	李源, 罗喜胜. 理想磁流体中激波与矩形密度界面相互作用的数值研究[J]. 计算物理, 2014, 31(6): 659-667.
[5]	李杰, 彭勇波. 磁流变液动态屈服行为的随机多尺度建模[J]. 计算物理, 2012, 29(1): 95-100.
[6]	吴昊. 一个磁流体力学中心型数值通量的耗散控制器[J]. 计算物理, 2010, 27(3): 361-368.
[7]	于翠影, 董志伟, 王玉芝, 王贵荣, 杨显俊, 赵强. 提高级联爆磁压缩发生器功率放大能力的新途径[J]. 计算物理, 2007, 24(6): 683-687.
[8]	朱自强, Islam Z, 朱一锟, 薛晓春. 流场/电磁场的数值优化计算[J]. 计算物理, 1998, 15(6): 648-653.
[9]	何启兵, 汤宇, 牟宗泽. 低杂波电流驱动模拟的数值方法[J]. 计算物理, 1996, 13(2): 159-166.
[10]	胡双辉, 黄林, 邱孝明. 湍流非线性问题的一种数值处理[J]. 计算物理, 1994, 11(2): 154-160.