Related papers: The Leray-G{\aa}rding method for finite difference…
In many applications of practical interest, solutions of partial differential equation models arise as critical points of an underlying (energy) functional. If such solutions are saddle points, rather than being maxima or minima, then the…
We introduce the Local Increasing Regularity Method (LIRM) which allows us to get from \emph{local} a priori estimates, on solutions $u$ of a linear equation $\displaystyle Du=\omega ,$ \emph{global} ones. As an application we shall prove…
We present an augmented Lagrangian trust-region method to efficiently solve constrained optimization problems governed by large-scale nonlinear systems with application to partial differential equation-constrained optimization. At each…
A linear implicit finite difference method is proposed for the approximation of the solution to a periodic, initial value problem for a Schrodinger-Hirota equation. Optimal, second order convergence in the discrete $H^1-$norm is proved,…
A scheme is presented for accurately propagating the gravitational field constraints in finite difference implementations of numerical relativity. The method is based on similar techniques used in astrophysical magnetohydrodynamics and…
This paper is concerned with the initial-boundary value problem for a nonlinear hyperbolic system of conservation laws. We study the boundary layers that may arise in approximations of entropy discontinuous solutions. We consider both the…
In this paper, we introduce a Lagrange multiplier approach to construct linearly implicit energy-preserving schemes of arbitrary order for general Hamiltonian PDEs. Unlike the widely used auxiliary variable methods, this novel approach does…
We consider the leapfrog algorithm by Noakes for computing geodesics on Riemannian manifolds. The main idea behind this algorithm is to subdivide the original endpoint geodesic problem into several local problems, for which the endpoint…
For finite element approximations of transport phenomena, it is often necessary to apply a form of limiting to ensure that the discrete solution remains well-behaved and satisfies physical constraints. However, these limiting procedures are…
Fourier multiplier analysis is developed for nonlocal peridynamic-type Laplace operators, which are defined for scalar fields in $\mathbb{R}^n$. The Fourier multipliers are given through an integral representation. We show that the integral…
Topology optimization problems generally support multiple local minima, and real-world applications are typically three-dimensional. In previous work [I. P. A. Papadopoulos, P. E. Farrell, and T. M. Surowiec, Computing multiple solutions of…
A novel trust region method for solving linearly constrained nonlinear programs is presented. The proposed technique is amenable to a distributed implementation, as its salient ingredient is an alternating projected gradient sweep in place…
In this note we design a cut finite element method for a low order divergence free element applied to a boundary value problem subject to Stokes' equations. For the imposition of Dirichlet boundary conditions we consider either Nitsche's…
In the second part of this series, we use the Lagrange multiplier approach proposed in the first part \cite{CheS21} to construct efficient and accurate bound and/or mass preserving schemes for a class of semi-linear and quasi-linear…
In this paper we apply an augmented Lagrange method to a class of semilinear elliptic optimal control problems with pointwise state constraints. We show strong convergence of subsequences of the primal variables to a local solution of the…
We are concerned with the convergence of spectral method for the numerical solution of the initial-boundary value problem associated to the Korteweg-de Vries-Kawahara equation (in short Kawahara equation), which is a transport equation…
The existing paradox between theory and computational experiment for weak solutions of systems of conservation laws in higher space dimensions is arguably resolved. Apparently successful computations are identified with underlying…
The stability theory for hyperbolic initial boundary value problems relies most of the time on the Laplace transform with respect to the time variable. For technical reasons, this usually restricts the validity of stability estimates to the…
Finding the optimal embedding of networks into low-dimensional hyperbolic spaces is a challenge that received considerable interest in recent years, with several different approaches proposed in the literature. In general, these methods…
This paper proposes a partially inexact alternating direction method of multipliers for computing approximate solution of a linearly constrained convex optimization problem. This method allows its first subproblem to be solved inexactly…