Numerical Analysis

Course "Numerical Analysis -- Special Topics" (2017-2018) -- Matlab Code

Here you can download two Matlab-based packages:

Code for linear (P1-P1) elements to solve Biot's problem.

Code for Taylor Hood (P2-P1) elements to solve Biot's problem.

After having unpacked the zip-files, one can start the session in Matlab by running

SpecialTopicsMesh ---- to generate the mesh
SpecialTopicsComp ---- to run the computational part (FEM)
SpecialTopicsPost ---- to make various output plots.

This is done for both the linear elements code and the Taylor-Hood code.

Regarding stabilisation, this should be built into the element matrices and vectors of the code for the linear elements. The last line in the file InputParameters.m, where beta is the stabilisation parameter, should be uncommented.

Last modified: February 1, 2018.