
We present a Virtual Element Method (VEM) for the solution of Dirichlet
problems for the quasilinear equation $\text{div} (k(u)\text{grad} u)=f$ with
essential boundary conditions. Within the VEM the nonlinear coefficient is
evaluated with the piecewise polynomial projection of the virtual element
ansatz.
Well posedness of the discrete problem and optimal order a priori error
estimates in the $H^1$ and $L^2$ norms are proven. In addition, the convergence
of fixed point iterations for the solution of the resulting nonlinear system is
established. Numerical examples confirm the convergence analysis.

Understanding how patterns and travelling waves form in chemical and
biological reactiondiffusion models is an area which has been widely
researched, yet is still experiencing fast development. Surprisingly enough, we
still do not have a clear understanding about all possible types of dynamical
regimes in classical reactiondiffusion models such as LotkaVolterra
competition models with spatial dependence. In this work, we demonstrate some
new types of wave propagation and pattern formation in a classical three
species cyclic competition model with spatial diffusion, which have been so far
missed in the literature. These new patterns are characterised by a high
regularity in space, but are different from patterns previously known to exist
in reactiondiffusion models, and may have important applications in improving
our understanding of biological pattern formation and invasion theory. Finding
these new patterns is made technically possible by using an automatic adaptive
finite element method driven by a novel a posteriori error estimate which is
proven to provide a reliable bound for the error of the numerical method. We
demonstrate how this numerical framework allows us to easily explore the
dynamical patterns both in two and three spatial dimensions.

We present a new $hp$version spacetime discontinuous Galerkin (dG) finite
element method for the numerical approximation of parabolic evolution equations
on general spatial meshes consisting of polygonal/polyhedral (polytopic)
elements, giving rise to prismatic spacetime elements. A key feature of the
proposed method is the use of spacetime elemental polynomial bases of
\emph{total} degree, say $p$, defined in the physical coordinate system, as
opposed to standard dGtimestepping methods whereby spatial elemental bases
are tensorized with temporal basis functions. This approach leads to a fully
discrete $hp$dG scheme using less degrees of freedom for each time step,
compared to standard dG timestepping schemes employing tensorized spacetime,
with acceptable deterioration of the approximation properties. A second key
feature of the new spacetime dG method is the incorporation of very general
spatial meshes consisting of possibly polygonal/polyhedral elements with
\emph{arbitrary} number of faces. A priori error bounds are shown for the
proposed method in various norms. An extensive comparison among the new
spacetime dG method, the (standard) tensorized spacetime dG methods, the
classical dGtimestepping, and conforming finite element method in space, is
presented in a series of numerical experiments.

An posteriori error analysis for the virtual element method (VEM) applied to
general elliptic problems is presented. The resulting error estimator is of
residualtype and applies on very general polygonal/polyhedral meshes. The
estimator is fully computable as it relies only on quantities available from
the VEM solution, namely its degrees of freedom and elementwise polynomial
projection. Upper and lower bounds of the error estimator with respect to the
VEM approximation error are proven. The error estimator is used to drive
adaptive mesh refinement in a number of test problems. Mesh adaptation is
particularly simple to implement since elements with consecutive coplanar
edges/faces are allowed and, therefore, locally adapted meshes do not require
any local mesh postprocessing.

We present the nonconforming Virtual Element Method (VEM) for the numerical
approximation of velocity and pressure in the steady Stokes problem. The
pressure is approximated using discontinuous piecewise polynomials, while each
component of the velocity is approximated using the nonconforming virtual
element space. On each mesh element the local virtual space contains the space
of polynomials of up to a given degree, plus suitable nonpolynomial functions.
The virtual element functions are implicitly defined as the solution of local
Poisson problems with polynomial Neumann boundary conditions. As typical in VEM
approaches, the explicit evaluation of the nonpolynomial functions is not
required. This approach makes it possible to construct nonconforming (virtual)
spaces for any polynomial degree regardless of the parity, for twoand
threedimensional problems, and for meshes with very general polygonal and
polyhedral elements. We show that the nonconforming VEM is infsup stable and
establish optimal a priori error estimates for the velocity and pressure
approximations. Numerical examples confirm the convergence analysis and the
effectiveness of the method in providing highorder accurate approximations.

We present in a unified framework new conforming and nonconforming Virtual
Element Methods (VEM) for general second order elliptic problems in two and
three dimensions. The differential operator is split into its symmetric and
nonsymmetric parts and conditions for stability and accuracy on their discrete
counterparts are established. These conditions are shown to lead to optimal
$H^1$ and $L^2$error estimates, confirmed by numerical experiments on a set
of polygonal meshes. The accuracy of the numerical approximation provided by
the two methods is shown to be comparable.

This work is concerned with the development of a spacetime adaptive
numerical method, based on a rigorous a posteriori error bound, for a
semilinear convectiondiffusion problem which may exhibit blowup in finite
time. More specifically, a posteriori error bounds are derived in the
$L^{\infty}(L^2)+L^2(H^1)$type norm for a first order in time
implicitexplicit (IMEX) interior penalty discontinuous Galerkin (dG) in space
discretization of the problem, although the theory presented is directly
applicable to the case of conforming finite element approximations in space.
The choice of the discretization in time is made based on a careful analysis of
adaptive time stepping methods for ODEs that exhibit finite time blowup. The
new adaptive algorithm is shown to accurately estimate the blowup time of a
number of problems, including one which exhibits regional blowup.

A discontinuous Galerkin (dG) method for the numerical solution of
initial/boundary value multicompartment partial differential equation (PDE)
models, interconnected with interface conditions, is presented and analysed.
The study of interface problems is motivated by models of mass transfer of
solutes through semipermeable membranes. More specifically, a model problem
consisting of a system of semilinear parabolic advectiondiffusionreaction
partial differential equations in each compartment, equipped with respective
initial and boundary conditions, is considered. Nonlinear interface conditions
modelling selective permeability, congestion and partial reflection are applied
to the compartment interfaces. An interior penalty dG method is presented for
this problem and it is analysed in the spacediscrete setting. The a priori
analysis shows that the method yields optimal a priori bounds, provided the
exact solution is sufficiently smooth. Numerical experiments indicate agreement
with the theoretical bounds and highlight the stability of the numerical method
in the advectiondominated regime.

This work is concerned with the derivation of a robust a posteriori error
estimator for a discontinuous Galerkin method discretisation of linear
nonstationary convectiondiffusion initial/boundary value problems and with
the implementation of a corresponding adaptive algorithm. More specifically, we
derive a posteriori bounds for the error in the $L^2(H^1)$type norm for an
interior penalty discontinuous Galerkin (dG) discretisation in space and a
backward Euler discretisation in time. An important feature of the estimator is
robustness with respect to the P\'{e}clet number of the problem which is
verified in practice by a series of numerical experiments. Finally, an adaptive
algorithm is proposed utilising the error estimator. Optimal rate of
convergence of the adaptive algorithm is observed in a number of test problems.

We consider a finite element method which couples the continuous Galerkin
method away from internal and boundary layers with a discontinuous Galerkin
method in the vicinity of layers. We prove that this consistent method is
stable in the streamline diffusion norm if the convection field flows
noncharacteristically from the region of the continuous Galerkin to the region
of the discontinuous Galerkin method. The stability properties of the coupled
method are illustrated by numerical experiments.

These pages review two families of mimetic finite difference methods: the
mixedtype methods presented in [Brezzi, Lipnikov, and Simoncini, M3AS, 2005]
and the nodal methods of [Brezzi, Buffa, and Lipnikov, M2AN, 2009]. The purpose
of this exercise it to highlight the similitudes underlying the construction of
the two families. The comparison prompts the definition of a piecewise linear
postprocessing of the nodal mimetic finite difference solution, as it was done
for the mixedtype method in [Cangiani and Manzini, CMAME, 2008].

In this short note we review deterministic simulation of biochemical
pathways, i.e. networks of biochemical reactions obeying the law of mass
action. It is meant as a basis for the MATLAB code, written by the author,
which permits easy input and simulation of general biochemical networks. This
work was carried out for the European Project `CardioWorkBench'.

We prove in an abstract setting that standard (continuous) Galerkin finite
element approximations are the limit of interior penalty discontinuous Galerkin
approximations as the penalty parameter tends to infinity. We apply this result
to equations of nonnegative characteristic form and the nonlinear, time
dependent system of incompressible miscible displacement. Moreover, we
investigate varying the penalty parameter on only a subset of a triangulation
and the effects of local superpenalization on the stability of the method,
resulting in a partly continuous, partly discontinuous method in the limit. An
iterative automatic procedure is also proposed for the determination of the
continuous region of the domain without loss of stability of the method.

For the stationary advectiondiffusion problem the standard continuous
Galerkin method is unstable without some additional control on the mesh or
method. The interior penalty discontinuous Galerkin method is stable but at the
expense of an increased number of degrees of freedom. The hybrid method
proposed in [5] combines the computational complexity of the continuous method
with the stability of the discontinuous method without a significant increase
in degrees of freedom. We discuss the implementation of this method using the
finite element library deal.ii and present some numerical experiments.