{"title":"球上不可压缩磁流体力学的时空Lie-Poisson离散化","authors":"Klas Modin, Michael Roop","doi":"10.1093/imanum/draf024","DOIUrl":null,"url":null,"abstract":"We give a structure preserving spatio-temporal discretization for incompressible magnetohydrodynamics (MHD) on the sphere. Discretization in space is based on the theory of geometric quantization, which yields a spatially discretized analogue of the MHD equations as a finite-dimensional Lie–Poisson system on the dual of the magnetic extension Lie algebra $\\mathfrak{f}=\\mathfrak{su}(N)\\ltimes \\mathfrak{su}(N)^{*}$. We also give accompanying structure preserving time discretizations for Lie–Poisson systems on the dual of semidirect product Lie algebras of the form $\\mathfrak{f}=\\mathfrak{g}\\ltimes \\mathfrak{g^{*}}$, where $\\mathfrak{g}$ is a $J$-quadratic Lie algebra. The time integration method is free of computationally costly matrix exponentials. We prove that the full method preserves a modified Lie–Poisson structure and corresponding Casimir functions, and that the modified structure and Casimirs converge to the continuous ones. The method is demonstrated for two models of magnetic fluids: incompressible MHD and Hazeltine’s model.","PeriodicalId":56295,"journal":{"name":"IMA Journal of Numerical Analysis","volume":"26 1","pages":""},"PeriodicalIF":2.3000,"publicationDate":"2025-04-26","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":"0","resultStr":"{\"title\":\"Spatio-temporal Lie–Poisson discretization for incompressible magnetohydrodynamics on the sphere\",\"authors\":\"Klas Modin, Michael Roop\",\"doi\":\"10.1093/imanum/draf024\",\"DOIUrl\":null,\"url\":null,\"abstract\":\"We give a structure preserving spatio-temporal discretization for incompressible magnetohydrodynamics (MHD) on the sphere. Discretization in space is based on the theory of geometric quantization, which yields a spatially discretized analogue of the MHD equations as a finite-dimensional Lie–Poisson system on the dual of the magnetic extension Lie algebra $\\\\mathfrak{f}=\\\\mathfrak{su}(N)\\\\ltimes \\\\mathfrak{su}(N)^{*}$. We also give accompanying structure preserving time discretizations for Lie–Poisson systems on the dual of semidirect product Lie algebras of the form $\\\\mathfrak{f}=\\\\mathfrak{g}\\\\ltimes \\\\mathfrak{g^{*}}$, where $\\\\mathfrak{g}$ is a $J$-quadratic Lie algebra. The time integration method is free of computationally costly matrix exponentials. We prove that the full method preserves a modified Lie–Poisson structure and corresponding Casimir functions, and that the modified structure and Casimirs converge to the continuous ones. The method is demonstrated for two models of magnetic fluids: incompressible MHD and Hazeltine’s model.\",\"PeriodicalId\":56295,\"journal\":{\"name\":\"IMA Journal of Numerical Analysis\",\"volume\":\"26 1\",\"pages\":\"\"},\"PeriodicalIF\":2.3000,\"publicationDate\":\"2025-04-26\",\"publicationTypes\":\"Journal Article\",\"fieldsOfStudy\":null,\"isOpenAccess\":false,\"openAccessPdf\":\"\",\"citationCount\":\"0\",\"resultStr\":null,\"platform\":\"Semanticscholar\",\"paperid\":null,\"PeriodicalName\":\"IMA Journal of Numerical Analysis\",\"FirstCategoryId\":\"100\",\"ListUrlMain\":\"https://doi.org/10.1093/imanum/draf024\",\"RegionNum\":2,\"RegionCategory\":\"数学\",\"ArticlePicture\":[],\"TitleCN\":null,\"AbstractTextCN\":null,\"PMCID\":null,\"EPubDate\":\"\",\"PubModel\":\"\",\"JCR\":\"Q1\",\"JCRName\":\"MATHEMATICS, APPLIED\",\"Score\":null,\"Total\":0}","platform":"Semanticscholar","paperid":null,"PeriodicalName":"IMA Journal of Numerical Analysis","FirstCategoryId":"100","ListUrlMain":"https://doi.org/10.1093/imanum/draf024","RegionNum":2,"RegionCategory":"数学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":null,"EPubDate":"","PubModel":"","JCR":"Q1","JCRName":"MATHEMATICS, APPLIED","Score":null,"Total":0}
Spatio-temporal Lie–Poisson discretization for incompressible magnetohydrodynamics on the sphere
We give a structure preserving spatio-temporal discretization for incompressible magnetohydrodynamics (MHD) on the sphere. Discretization in space is based on the theory of geometric quantization, which yields a spatially discretized analogue of the MHD equations as a finite-dimensional Lie–Poisson system on the dual of the magnetic extension Lie algebra $\mathfrak{f}=\mathfrak{su}(N)\ltimes \mathfrak{su}(N)^{*}$. We also give accompanying structure preserving time discretizations for Lie–Poisson systems on the dual of semidirect product Lie algebras of the form $\mathfrak{f}=\mathfrak{g}\ltimes \mathfrak{g^{*}}$, where $\mathfrak{g}$ is a $J$-quadratic Lie algebra. The time integration method is free of computationally costly matrix exponentials. We prove that the full method preserves a modified Lie–Poisson structure and corresponding Casimir functions, and that the modified structure and Casimirs converge to the continuous ones. The method is demonstrated for two models of magnetic fluids: incompressible MHD and Hazeltine’s model.
期刊介绍:
The IMA Journal of Numerical Analysis (IMAJNA) publishes original contributions to all fields of numerical analysis; articles will be accepted which treat the theory, development or use of practical algorithms and interactions between these aspects. Occasional survey articles are also published.