Abstract:We consider the design and analysis of numerical methods for approximating positive solutions to nonlinear geometric elliptic partial differential equations containing critical exponents. This class of problems includes the Yamabe problem and the Einstein constraint equations, which simultaneously contain several challenging features: high spatial dimension n >= 3, varying (potentially non-smooth) coefficients, critical (even super-critical) nonlinearity, non-monotone nonlinearity (arising from a non-convex energy), and spatial domains that are typically Riemannian manifolds rather than simply open sets in Rn. These problems may exhibit multiple solutions, although only positive solutions typically have meaning. This creates additional complexities in both the theory and numerical treatment of such problems, as this feature introduces both non-uniqueness as well as the need to incorporate an inequality constraint into the formulation. In this work, we consider numerical methods based on Galerkin-type discretization, covering any standard bases construction (finite element, spectral, or wavelet), and the combination of a barrier method for nonconvex optimization and global inexact Newton-type methods for dealing with nonconvexity and the presence of inequality constraints. We first give an overview of barrier methods in non-convex optimization, and then develop and analyze both a primal barrier energy method for this class of problems. We then consider a sequence of numerical experiments using this type of barrier method, based on a particular Galerkin method, namely the piecewise linear finite element method, leverage the FETK modeling package. We illustrate the behavior of the primal barrier energy method for several examples, including the Yamabe problem and the Hamiltonian constraint.

The Energy Technique for the Six-Step BDF Method

The weighted and shifted seven-step BDF method for parabolic equations

The variable two-step BDF method for parabolic equations

A NEW MULTI-STEP BDF ENERGY STABLE TECHNIQUE FOR THE EXTENDED FISHER-KOLMOGOROV EQUATION

Superconvergence analysis of an energy stable scheme for nonlinear reaction-diffusion equation with BDF mixed FEM

Stability and error estimates for the variable step-size BDF2 method for linear and semilinear parabolic equations

Stability and convergence of BDF2-ADI schemes with variable step sizes for parabolic equation

On energy stable, maximum-principle preserving, second order BDF scheme with variable steps for the Allen-Cahn equation

Energy-Stable Global Radial Basis Function Methods on Summation-By-Parts Form

A Positivity-Preserving, Energy Stable BDF2 Scheme with Variable Steps for the Cahn–Hilliard Equation with Logarithmic Potential

An enhanced and highly efficient semi-implicit combined Lagrange multiplier approach with preserving original energy law for dissipative systems

A revisit of the energy quadratization method with a relaxation technique

Decoupling multistep schemes for elliptic-parabolic problems

Von Neumann Stability Analysis for Multi-level Multi-step Methods

Barrier methods for critical exponent problems in geometric analysis and mathematical physics

Numerical and theoretical discussions for solving nonlinear generalized Benjamin–Bona–Mahony–Burgers equation based on the Legendre spectral element method

Structure-preserving BDF2 FE method for the coupled Schrödinger-Boussinesq equations

A comparison of duality and energy aposteriori estimates for L?(0,T;L2(Ω)) in parabolic problems

A Roadmap for Discretely Energy-Stable Schemes for Dissipative Systems Based on a Generalized Auxiliary Variable with Guaranteed Positivity

Unconditionally optimal error estimate of a linearized variable-time-step BDF2 scheme for nonlinear parabolic equations

Unconditional superconvergence analysis of a modified nonconforming energy stable BDF2 FEM for Sobolev equations with Burgers' type nonlinearity