Error analysis for discretizations of parabolic problems using continuous finite elements in time and mixed finite elements in space.

Variational time discretization schemes are getting of increasing importance for the accurate numerical approximation of transient phenomena. The applicability and value of mixed finite element methods in space for simulating transport processes have been demonstrated in a wide class of works. We consider a family of continuous Galerkin-Petrov time discretization schemes that is combined with a mixed finite element approximation of the spatial variables. The existence and uniqueness of the semidiscrete approximation and of the fully discrete solution are established. For this, the Banach-Necas-Babuska theorem is applied in a non-standard way. Error estimates with explicit rates of convergence are proved for the scalar and vector-valued variable. An optimal order estimate in space and time is proved by duality techniques for the scalar variable. The convergence rates are analyzed and illustrated by numerical experiments, also on stochastically perturbed meshes.

Modeling of drug release from collagen matrices.

Drug release from collagen matrices is in most cases governed by diffusion from swollen matrices but also enzymatic matrix degradation or hydrophobic drug/collagen interactions may contribute. To reduce water uptake and to prolong the release, insoluble collagen matrices have been chemically or dehydrothermally crosslinked. Assuming Fickian diffusion a one-dimensional model was developed and tested that allows description of water penetration, swelling and drug release and that may be expanded considering a subsequent erosion process or interactions. Swelling is described by a volume balance. For dry collagen matrices crosslinked by thermal treatment the existence of a moving front separating the polymer from a gel phase was considered, and a convective term induced by the volume expansion was incorporated. The resulting moving boundary problem was solved using a method based on biquadratic finite elements in both space and time that is stable, shows high accuracy, and is suitable for solving problems with a singularity at the initial time point. The model was verified for insoluble collagen matrices at different crosslinking degrees for both chemical and thermal treatment. For constant diffusion coefficients a close form of the solution was derived yielding equivalent results to the numerical approach.