Numerical Methods for PDEs: Robustness and Efficiency

Report on Current Developments in the Research Area

General Direction of the Field

The recent developments in the research area are marked by a significant push towards enhancing the robustness, efficiency, and theoretical soundness of numerical methods for solving complex partial differential equations (PDEs). The field is witnessing a convergence of advanced computational techniques, rigorous error analysis, and innovative algorithmic frameworks, all aimed at addressing the intricacies of nonlinear and nonsmooth problems.

One of the key trends is the integration of adaptive mesh refinement techniques with advanced numerical methods, such as the virtual element method (VEM) and the finite element method (FEM). These methods are being tailored to handle quasilinear elliptic PDEs and other challenging problems, with a focus on deriving computable error estimators that drive adaptive refinement strategies. This approach not only improves the accuracy of numerical solutions but also optimizes computational resources by concentrating refinement efforts where the error is most significant.

Another notable direction is the development of globalized and inexact semismooth Newton methods for solving nonsmooth fixed-point equations and variational inequalities. These methods are being extended to Banach spaces and are shown to achieve q-superlinear convergence under certain contraction assumptions. The ability to handle inexact function evaluations and Newton steps, combined with globalization techniques, makes these methods highly versatile for a wide range of applications, including those involving semilinear PDEs with nondifferentiable operators.

The field is also seeing a strong emphasis on formalizing the mathematical foundations of numerical methods, particularly in the context of the finite element method. This trend reflects a growing interest in ensuring the highest level of confidence in numerical simulations by formalizing proofs in rigorous mathematical frameworks, such as Coq. This approach is crucial for establishing the soundness of numerical algorithms and for facilitating their adoption in critical applications.

Noteworthy Developments

  • Virtual Element Method (VEM) for Quasilinear Elliptic PDEs: The development of a computable error estimator for VEM, coupled with adaptive mesh refinement, represents a significant advancement in handling complex PDEs.

  • Globalized Inexact Semismooth Newton Method: The extension of this method to Banach spaces and its application to nonsmooth problems, including those involving semilinear PDEs, is a notable contribution to the field.

  • Formalization of Finite Element Method in Coq: The detailed formal proofs provided for the construction of Lagrange finite elements in Coq underscore the importance of rigorous mathematical foundations in numerical methods.

Sources

A posteriori error analysis of the virtual element method for second-order quasilinear elliptic PDEs

A Globalized Inexact Semismooth Newton Method for Nonsmooth Fixed-point Equations involving Variational Inequalities

On inf-sup stability and optimal convergence of the quasi-reversibility method for unique continuation subject to Poisson's equation

Finite element method. Detailed proofs to be formalized in Coq

Built with on top of