Change search
ReferencesLink to record
Permanent link

Direct link
Gryphon - a Module for Time Integration of Partial Differential Equations in FEniCS
Norwegian University of Science and Technology, Faculty of Information Technology, Mathematics and Electrical Engineering, Department of Mathematical Sciences.
2012 (English)MasteroppgaveStudent thesis
Abstract [en]

This thesis aims to implement time integrators in the FEniCS framework. More specifically, the thesis focuses on selecting suitable time integrators, implement these and verify that the implementation works by applying them to various relevant test problems. This work resulted in a module for FEniCS, named Gryphon. The thesis is divided into four parts. The first part builds a theoretical framework which will motivate why singly diagonally implicit Runge-Kutta methods with an explicit first stage (ESDIRKs) should be considered for solving stiff ordinary differential equations (ODEs). It will also be shown how an ESDIRK method can be utilized to solve time dependent partial differential equations (PDEs) by solving the semidiscretized system arising from first applying a finite element method. We will restrict our attention to PDEs which either give rise to a pure ODE system or a DAE (differential-algebraic equation) system of index 1. The second part discusses the implementation of Gryphon, focusing on why such a module is useful and how the source code is structured. The third part is devoted to numerical experiments on the ESDIRK solvers implemented in Gryphon. The experiments will establish convergence and give some run-time statistics for various ESDIRK schemes. We will also see that L-stability is a favorable trait when working with stiff equations, by comparing an ESDIRK method to the trapezoidal rule. It will also be verified that the step size selectors implemented in Gryphon behaves as expected. As test problems we consider the heat equation, the Fisher-Kolmogorov equation, the Gray-Scott equations, the Fitzhugh-Nagumo equations and the Cahn-Hilliard equations. The fourth part is a user manual for Gryphon. All the parameters which can be changed by the user are explained. The manual also includes example code for solving the heat equation, the Gray-Scott equations and the Cahn-Hilliard equation, to get the reader starting on solving their own problems.

Place, publisher, year, edition, pages
Institutt for matematiske fag , 2012. , 73 p.
Keyword [no]
ntnudaim:7469, MTFYMA fysikk og matematikk, Industriell matematikk
URN: urn:nbn:no:ntnu:diva-18689Local ID: ntnudaim:7469OAI: diva2:566188
Available from: 2012-11-08 Created: 2012-11-08

Open Access in DiVA

fulltext(4498 kB)855 downloads
File information
File name FULLTEXT01.pdfFile size 4498 kBChecksum SHA-512
Type fulltextMimetype application/pdf
cover(184 kB)32 downloads
File information
File name COVER01.pdfFile size 184 kBChecksum SHA-512
Type coverMimetype application/pdf
attachment(289 kB)36 downloads
File information
File name ATTACHMENT01.zipFile size 289 kBChecksum SHA-512
Type attachmentMimetype application/zip

By organisation
Department of Mathematical Sciences

Search outside of DiVA

GoogleGoogle Scholar
Total: 855 downloads
The number of downloads is the sum of all downloads of full texts. It may include eg previous versions that are now no longer available

Total: 89 hits
ReferencesLink to record
Permanent link

Direct link