A study of blow-ups in the Keller-Segel model of chemotaxis

We study the Keller-Segel model of chemotaxis and develop a composite particle-grid numerical method with adaptive time stepping which allows us to accurately resolve singular solutions. The numerical findings (in two dimensions) are then compared with analytical predictions regarding formation and interaction of singularities obtained via analysis of the stochastic differential equations associated with the Keller-Segel model

An adaptive octree finite element method for PDEs posed on surfaces

The paper develops a finite element method for partial differential equations posed on hypersurfaces in $\mathbb{R}^N$, $N=2,3$. The method uses traces of bulk finite element functions on a surface embedded in a volumetric domain. The bulk finite element space is defined on an octree grid which is locally refined or coarsened depending on error indicators and estimated values of the surface curvatures. The cartesian structure of the bulk mesh leads to easy and efficient adaptation process, while the trace finite element method makes fitting the mesh to the surface unnecessary. The number of degrees of freedom involved in computations is consistent with the two-dimension nature of surface PDEs. No parametrization of the surface is required; it can be given implicitly by a level set function. In practice, a variant of the marching cubes method is used to recover the surface with the second order accuracy. We prove the optimal order of accuracy for the trace finite element method in $H^1$ and $L^2$ surface norms for a problem with smooth solution and quasi-uniform mesh refinement. Experiments with less regular problems demonstrate optimal convergence with respect to the number of degrees of freedom, if grid adaptation is based on an appropriate error indicator. The paper shows results of numerical experiments for a variety of geometries and problems, including advection-diffusion equations on surfaces. Analysis and numerical results of the paper suggest that combination of cartesian adaptive meshes and the unfitted (trace) finite elements provide simple, efficient, and reliable tool for numerical treatment of PDEs posed on surfaces

Space-time adaptive finite elements for nonlocal parabolic variational inequalities

This article considers the error analysis of finite element discretizations and adaptive mesh refinement procedures for nonlocal dynamic contact and friction, both in the domain and on the boundary. For a large class of parabolic variational inequalities associated to the fractional Laplacian we obtain a priori and a posteriori error estimates and study the resulting space-time adaptive mesh-refinement procedures. Particular emphasis is placed on mixed formulations, which include the contact forces as a Lagrange multiplier. Corresponding results are presented for elliptic problems. Our numerical experiments for $2$-dimensional model problems confirm the theoretical results: They indicate the efficiency of the a posteriori error estimates and illustrate the convergence properties of space-time adaptive, as well as uniform and graded discretizations.Comment: 47 pages, 20 figure