Guardado en:
| Autores principales: | , , |
|---|---|
| Formato: | Preprint |
| Publicado: |
2025
|
| Materias: | |
| Acceso en línea: | https://arxiv.org/abs/2508.20255 |
| Etiquetas: |
Agregar Etiqueta
Sin Etiquetas, Sea el primero en etiquetar este registro!
|
| _version_ | 1866912557572292608 |
|---|---|
| author | Kirby, Robert C. MacLachlan, Scott P. Brubeck, Pablo D. |
| author_facet | Kirby, Robert C. MacLachlan, Scott P. Brubeck, Pablo D. |
| contents | Irksome is a library based on the Unified Form Language (UFL) that automates the application of Runge-Kutta time-stepping methods for finite element spatial discretizations of partial differential equations (PDEs). This paper describes recent updates to Irksome that allow users to express semidiscrete forms of PDEs that contain second-order temporal derivatives, whence it generates stage-coupled variational problems to be solved at each time step for Runge-Kutta-Nyström methods. Firedrake then generates code for these variational problems and provides a rich interface to PETSc for solving them. Directly discretizing second-order time derivatives with Runge-Kutta-Nyström methods provides several advantages relative to discretizing a rewritten first-order system with a standard Runge-Kutta method. Besides working with an interface closer to the problem formulation in UFL, avoiding these auxiliary variables means that Runge-Kutta-Nyström methods lead to smaller algebraic systems and better run-time. Our numerical results indicate that, with effective preconditioning, fully implicit Runge-Kutta-Nyström methods can be made competitive with more traditional explicit methods for wave equations. They are also (essentially) required to discretize wave-type equations with higher-order spatial derivatives. We also provide numerical experiments for fully dynamic poroelasticity, a system of mixed temporal order, where our time-stepping and algebraic solvers perform effectively even as we approach the incompressible limit. |
| format | Preprint |
| id |
arxiv_https___arxiv_org_abs_2508_20255 |
| institution | arXiv |
| publishDate | 2025 |
| record_format | arxiv |
| spellingShingle | Automated Runge-Kutta-Nyström time stepping for finite element methods in Irksome Kirby, Robert C. MacLachlan, Scott P. Brubeck, Pablo D. Numerical Analysis Irksome is a library based on the Unified Form Language (UFL) that automates the application of Runge-Kutta time-stepping methods for finite element spatial discretizations of partial differential equations (PDEs). This paper describes recent updates to Irksome that allow users to express semidiscrete forms of PDEs that contain second-order temporal derivatives, whence it generates stage-coupled variational problems to be solved at each time step for Runge-Kutta-Nyström methods. Firedrake then generates code for these variational problems and provides a rich interface to PETSc for solving them. Directly discretizing second-order time derivatives with Runge-Kutta-Nyström methods provides several advantages relative to discretizing a rewritten first-order system with a standard Runge-Kutta method. Besides working with an interface closer to the problem formulation in UFL, avoiding these auxiliary variables means that Runge-Kutta-Nyström methods lead to smaller algebraic systems and better run-time. Our numerical results indicate that, with effective preconditioning, fully implicit Runge-Kutta-Nyström methods can be made competitive with more traditional explicit methods for wave equations. They are also (essentially) required to discretize wave-type equations with higher-order spatial derivatives. We also provide numerical experiments for fully dynamic poroelasticity, a system of mixed temporal order, where our time-stepping and algebraic solvers perform effectively even as we approach the incompressible limit. |
| title | Automated Runge-Kutta-Nyström time stepping for finite element methods in Irksome |
| topic | Numerical Analysis |
| url | https://arxiv.org/abs/2508.20255 |