Eigenmode initialisation of 2D (magneto)hydrodynamic simulations
Fuente:
arXiv
Salvato in:
| Autori principali: | , |
|---|---|
| Natura: | Preprint |
| Pubblicazione: |
2026
|
| Soggetti: | |
| Accesso online: | |
| Tags: |
Aggiungi Tag
Nessun Tag, puoi essere il primo ad aggiungerne!!
|
| _version_ | 1866912930198454272 |
|---|---|
| author | De Jonghe, Jordi Russell, Alexander J. B. |
| author_facet | De Jonghe, Jordi Russell, Alexander J. B. |
| contents | The early evolution of unstable hydrodynamic and magnetohydrodynamic equilibria is often governed by a few dominant linear eigenmodes. We investigate whether initialising a simulation with a superposition of linear eigenmodes that contains the most unstable mode saves computation time, and how the selection of the included modes affects the non-linear evolution. Using the non-linear (magneto)hydrodynamic simulation code MPI-AMRVAC, the evolutions of a flow-sheared fluid interface, a Harris current sheet, and a flow-sheared plasma interface were simulated for various initial perturbations. The novel initial perturbations were linear eigenmodes of the equilibrium, or superpositions thereof, and calculated with the Legolas code. We benchmarked to initialisation with velocity noise and, in the case of the Harris sheet, initialisation with an analytic magnetic field perturbation. By perturbing an unstable equilibrium with a superposition of linear eigenmodes that contains the most unstable mode, significantly less computation time is spent on the linear stage of the evolution compared to traditional perturbation methods. In the best case, the simulation time needed to reach the non-linear stage is reduced by an order of magnitude. The inclusion or omission of certain modes in the initial perturbation is observed to affect the non-linear evolution to various degrees. The perturbation of equilibria with a superposition of linear eigenmodes that contains the most unstable mode allows simulations to reach a late-evolution stage faster, thus saving computation time. Additional eigenmodes can be included alongside the fastest growing mode to obtain further benefits, for example, to accelerate symmetry breaking in the non-linear stage, or to delay their effect on the non-linear evolution. Coupling spectroscopic codes with (magneto)hydrodynamic codes therefore offers significant advantages. |
| format | Preprint |
| id |
arxiv_https___arxiv_org_abs_2602_23849 |
| institution | arXiv |
| publishDate | 2026 |
| record_format | arxiv |
| spellingShingle | Eigenmode initialisation of 2D (magneto)hydrodynamic simulations De Jonghe, Jordi Russell, Alexander J. B. Instrumentation and Methods for Astrophysics Solar and Stellar Astrophysics The early evolution of unstable hydrodynamic and magnetohydrodynamic equilibria is often governed by a few dominant linear eigenmodes. We investigate whether initialising a simulation with a superposition of linear eigenmodes that contains the most unstable mode saves computation time, and how the selection of the included modes affects the non-linear evolution. Using the non-linear (magneto)hydrodynamic simulation code MPI-AMRVAC, the evolutions of a flow-sheared fluid interface, a Harris current sheet, and a flow-sheared plasma interface were simulated for various initial perturbations. The novel initial perturbations were linear eigenmodes of the equilibrium, or superpositions thereof, and calculated with the Legolas code. We benchmarked to initialisation with velocity noise and, in the case of the Harris sheet, initialisation with an analytic magnetic field perturbation. By perturbing an unstable equilibrium with a superposition of linear eigenmodes that contains the most unstable mode, significantly less computation time is spent on the linear stage of the evolution compared to traditional perturbation methods. In the best case, the simulation time needed to reach the non-linear stage is reduced by an order of magnitude. The inclusion or omission of certain modes in the initial perturbation is observed to affect the non-linear evolution to various degrees. The perturbation of equilibria with a superposition of linear eigenmodes that contains the most unstable mode allows simulations to reach a late-evolution stage faster, thus saving computation time. Additional eigenmodes can be included alongside the fastest growing mode to obtain further benefits, for example, to accelerate symmetry breaking in the non-linear stage, or to delay their effect on the non-linear evolution. Coupling spectroscopic codes with (magneto)hydrodynamic codes therefore offers significant advantages. |
| title | Eigenmode initialisation of 2D (magneto)hydrodynamic simulations |
| topic | Instrumentation and Methods for Astrophysics Solar and Stellar Astrophysics |
| url | https://arxiv.org/abs/2602.23849 |