Physics > Computational Physics
[Submitted on 8 Jul 2007 (v1), last revised 7 Jan 2008 (this version, v2)]
Title:An iterative semi-implicit scheme with robust damping
View PDFAbstract: An efficient, iterative semi-implicit (SI) numerical method for the time integration of stiff wave systems is presented. Physics-based assumptions are used to derive a convergent iterative formulation of the SI scheme which enables the monitoring and control of the error introduced by the SI operator. This iteration essentially turns a semi-implicit method into a fully implicit method. Accuracy, rather than stability, determines the timestep. The scheme is second-order accurate and shown to be equivalent to a simple preconditioning method. We show how the diffusion operators can be handled so as to yield the property of robust damping, i.e., dissipating the solution at all values of the parameter $\mathcal D\dt$, where $\mathcal D$ is a diffusion operator and $\dt$ the timestep. The overall scheme remains second-order accurate even if the advection and diffusion operators do not commute. In the limit of no physical dissipation, and for a linear test wave problem, the method is shown to be symplectic. The method is tested on the problem of Kinetic Alfvén wave mediated magnetic reconnection. A Fourier (pseudo-spectral) representation is used. A 2-field gyrofluid model is used and an efficacious k-space SI operator for this problem is demonstrated. CPU speed-up factors over a CFL-limited explicit algorithm ranging from $\sim20$ to several hundreds are obtained, while accurately capturing the results of an explicit integration. Possible extension of these results to a real-space (grid) discretization is discussed.
Submission history
From: Gregory W. Hammett [view email][v1] Sun, 8 Jul 2007 03:56:15 UTC (743 KB)
[v2] Mon, 7 Jan 2008 03:23:50 UTC (756 KB)
Current browse context:
physics.comp-ph
Change to browse by:
References & Citations
Bibliographic and Citation Tools
Bibliographic Explorer (What is the Explorer?)
Litmaps (What is Litmaps?)
scite Smart Citations (What are Smart Citations?)
Code, Data and Media Associated with this Article
alphaXiv (What is alphaXiv?)
CatalyzeX Code Finder for Papers (What is CatalyzeX?)
DagsHub (What is DagsHub?)
Gotit.pub (What is GotitPub?)
Hugging Face (What is Huggingface?)
Papers with Code (What is Papers with Code?)
ScienceCast (What is ScienceCast?)
Demos
Recommenders and Search Tools
Influence Flower (What are Influence Flowers?)
Connected Papers (What is Connected Papers?)
CORE Recommender (What is CORE?)
arXivLabs: experimental projects with community collaborators
arXivLabs is a framework that allows collaborators to develop and share new arXiv features directly on our website.
Both individuals and organizations that work with arXivLabs have embraced and accepted our values of openness, community, excellence, and user data privacy. arXiv is committed to these values and only works with partners that adhere to them.
Have an idea for a project that will add value for arXiv's community? Learn more about arXivLabs.