Name
Title | ||
|---|---|---|
A spectral Galerkin method for the coupled Orr-Sommerfeld and induction equations for free-surface MHD |
Abstract | ||
|---|---|---|
We develop and test spectral Galerkin schemes to solve the coupled Orr-Sommerfeld and induction equations for parallel, incompressible MHD in free-surface and fixed-boundary geometries. The schemes' discrete bases consist of Legendre internal shape functions, supplemented with nodal shape functions for the weak imposition of the stress and insulating boundary conditions. The orthogonality properties of the basis polynomials solve the matrix-coefficient growth problem, and eigenvalue-eigenfunction pairs can be computed stably at spectral orders at least as large as p=3000 with p-independent roundoff error. Accuracy is limited instead by roundoff sensitivity due to non-normality of the stability operators at large hydrodynamic and/or magnetic Reynolds numbers (Re,Rm@?4x10^4). In problems with Hartmann velocity and magnetic-field profiles we employ suitable Gauss quadrature rules to evaluate the associated exponentially weighted sesquilinear forms without error. An alternative approach, which involves approximating the forms by means of Legendre-Gauss-Lobatto quadrature at the 2p-1 precision level, is found to yield equal eigenvalues within roundoff error. As a consistency check, we compare modal growth rates to energy growth rates in nonlinear simulations and record relative discrepancy smaller than 10^-^5 for the least stable mode in free-surface flow at Re=3x10^4. Moreover, we confirm that the computed normal modes satisfy an energy conservation law for free-surface MHD with error smaller than 10^-^6. The critical Reynolds number in free-surface MHD is found to be sensitive to the magnetic Prandtl number Pm, even at the Pm=O(10^-^5) regime of liquid metals. |
Year | DOI | Venue |
|---|---|---|
2009 | 10.1016/j.jcp.2008.10.016 | J. Comput. Physics |
Keywords | Field | DocType |
free-surface mhd,roundoff sensitivity,induction equation,incompressible mhd,eigenvalue problems,free-surface flow,legendre internal shape function,spectral galerkin method,65l60,eigenvalue problems spectral galerkin method hydrodynamic stability orr-sommerfeld equations free-surface mhd,p-independent roundoff error,76e05,roundoff error,76e17,orr–sommerfeld equations,modal growth rate,76e25,65l15,matrix-coefficient growth problem,hydrodynamic stability,energy growth rate,liquid metal,free surface,galerkin method,boundary condition,reynolds number,magnetic field,quadrature rule,normal modes,satisfiability,energy conservation,sesquilinear form,eigenvalues | Reynolds number,Hydrodynamic stability,Round-off error,Mathematical analysis,Spectral method,Gaussian quadrature,Magnetic Reynolds number,Mathematics,Legendre function,Induction equation | Journal |
Volume | Issue | ISSN |
228 | 4 | Journal of Computational Physics |
Citations | PageRank | References |
0 | 0.34 | 3 |
Authors | ||
3 | ||
Name | Order | Citations | PageRank |
|---|---|---|---|
| Dimitrios Giannakis | 1 | 5 | 2.60 |
| Paul F. Fischer | 2 | 181 | 23.03 |
| Robert Rosner | 3 | 25 | 2.76 |