A vector spherical harmonic spectral code for linearised magnetohydrodynamics

D. J. Ivers, C. G. Phillips


Linearised rotating magnetohydrodynamic stability code for the steady axisymmetric basic states of an electrically conducting fluid sphere is described. The code generates compact hybrid angular spectral forms of the magnetic induction, heat and Boussinesq Navier-Stokes equations, using toroidal and poloidal representations of the perturbation vector fields, and vector or scalar spherical harmonic expansions of all fields. The momentum equation may include inertial, Coriolis, buoyancy, viscous and magnetic Lorentz forces. Three subroutines evaluate the spectral interactions of products. There are only six radial functions, which are discretised using uniform second-order finite differences. The resulting large scale complex non-hermitian generalised eigen- and critical-value problems are solved using inverse and Newton-Raphson iteration methods, respectively. Test results are presented for several models.

Full Text:


DOI: https://doi.org/10.21914/anziamj.v44i0.689

Remember, for most actions you have to record/upload into this online system
and then inform the editor/author via clicking on an email icon or Completion button.
ANZIAM Journal, ISSN 1446-8735, copyright Australian Mathematical Society.