Adaptive Algorithm for Linear Parabolic Problems for 2D Polygonal Domain

Algorithm solves linear parabolic problem with space-time dependent coefficients adaptively in time

Description

Solves FEM analogue to parabolic equation:

d/dt(u)-div(D(x,t)u) +C*Nabla(u)+R*u= f(x,t) u=u(x,t) scalar fct

u(x,t) = u_d(t) on Dirichlet-edge
D(x,t)*Nabla(u)*n(x) = g(x,t) on Neumann-edge.

1)Data can be space and time dependent.
2)Spatial and temporal grid is refined and coarsened adaptively and indepenently of the user.
3)Data is plugged-in via an easy to adapt driver file.
4)Extensive use of vectorisation methods.
5)Input files: finite element grid(coordinates.dat,elements.dat,dirichlet.dat,neumann.dat); driver file
6)Output file: Solution.mat
7)See doc/SpaceTimeAdaptiveAlgorithmForLinearParabolicProblems.pdf for documentation
8) demo.m for demo; help solveP.m for direct help

• MATLAB release: MATLAB 7.4 (R2007a)