Adaptive computation of fourth-order problems

Carsten Carstensen: Monday, 10th September at 11:00 Abstract: The popular (piecewise) quadratic schemes for the fourth-order plate bending problems based on triangles are the nonconforming Morley finite element, two discontinuous Galerkin, the C0 interior penalty, and the WOPSIP schemes. The first part of the presentation discusses recent applications to the linear bi-Laplacian and to semi-linear fourth-order problems like the stream function vorticity formulation of incompressible 2D Navier-Stokes problem and the von Kármán plate bending problem. The role of a smoother is emphasised and reliable and efficient a posteriori error estimators give rise to adaptive mesh-refining strategies that recover optimal rates in numerical experiments. The last part addresses recent developments on adaptive multilevel Argyris finite element methods. The presentation is based on joint work with B. Gräßle, (Humboldt) and N. Nataraj (IITB, Mumbai) partly reflected in the references below. REFERENCES [1] C. Carstensen, B. Gräßle,, and N. Nataraj. Unifying a posteriori error analysis of five piecewise quadratic discretisations for the biharmonic equation, J. Numer. Math., volume 32, pp. 77–109, 2024, arXiv:2310.05648. [2] C. Carstensen, B. Gräßle,, and N. Nataraj. A posteriori error control for fourth-order semilinear problems with quadratic nonlinearity, SIAM J. Numer. Anal., volume 62, pp. 919–945, 2024. [3] C. Carstensen, Jun Hu. Hierarchical Argyris finite element method for adaptive and multigrid algorithms, Comput. Methods Appl. Math., volume 21, pp. 529–556, 2021. [4] C. Carstensen, N. Nataraj. A Priori and a Posteriori Error Analysis of the Crouzeix–Raviart and Morley FEM with Original and Modified Right-Hand Sides, Comput. Methods Appl. Math., volume 21, pp. 289–315, 2021. [5] C. Carstensen, N. Nataraj, G.C. Remesan, D. Shylaja. Lowest-order FEM for fourth-order semi-linear problems with trilinear nonlinearity, Numerische Mathematik 154, pp. 323–368, 2023. [6] C. Carstensen, N. Nataraj. Lowest-order equivalent nonstandard finite element methods for biharmonic plates, ESAIM: Mathematical Modelling and Numerical Analysis, 56(1),…

Continue reading

Computable reliable bounds for Poincaré–Friedrichs constants via Čech–de-Rham complexes

Martin Licht: Thursday, 25th April at 11:00 Abstract: We derive computable and reliable upper bounds for Poincaré–Friedrichs constants of classical Sobolev spaces and, more generally, Sobolev de-Rham complexes. The upper bounds are in terms of local Poincaré–Friedrichs constants over subdomains and the smallest singular value of a finite-dimensional operator that is easily assembled from the geometric setting. Thus we reduce the computational effort when computing the Poincaré–Friedrichs constant of finite de-Rham complexes, and we provide computable reliable bounds even for the original Sobolev de-Rham complex. The reduction to a finite-dimensional system uses diagram chasing within a Čech–de-Rham complex.

Continue reading

Homogeneous multigrid for hybrid discretizations: application to HHO methods

Andreas Rupp: Tuesday, 2nd April at 14:00 Abstract: We prove the uniform convergence of the geometric multigrid V-cycle for hybrid high-order (HHO) and other discontinuous skeletal methods. Our results generalize previously established results for HDG methods, and our multigrid method uses standard smoothers and local solvers that are bounded, convergent, and consistent. We use a weak version of elliptic regularity in our proofs. Numerical experiments confirm our theoretical results.

Continue reading

A localized orthogonal decomposition strategy for hybrid discontinuous Galerkin methods

Roland Maier: Thursday, 4th April at 11:00 Abstract: We formulate and analyze a multiscale method for an elliptic problem with an oscillatory coefficient based on a skeletal (hybrid) formulation. More precisely, we employ hybrid discontinuous Galerkin approaches and combine them with the localized orthogonal decomposition methodology to obtain a coarse-scale skeletal method that effectively includes fine-scale information. This work is a first step to reliably merge hybrid skeletal formulations and localized orthogonal decomposition and unite the advantages of both strategies. Numerical experiments are presented to illustrate the theoretical findings.

Continue reading

Modeling some biological phenomena via the porous media approach

Zoubida Mghazli: Thursday, 18rd Jan at 14:00 Abstract: Many biological systems can be modeled by the ”porous medium approach”, such as the diffusion of nutrients and other macromolecules through and in biological tissues. In this presentation, after a brief introduction to the ”porous medium approach”, we present some biological systems viewed through this approach. This will mainly concern the process of biodegradation of household waste, the Trichoderma fungi and the flow of water in the plant root.

Continue reading