Numerical Methods for Linear Global Stability of Magnetohydrodynamic Duct Flows

Abstract

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

CLC Number:

O361

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.

References

[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.