Abstract
We analytically and numerically solve the PCP equation, u(xx)u(yy) = 1, with homogeneous Dirichlet boundary conditions on the unit square. Chebyshev and Fourier spectral methods with low degree truncations yield moderate accuracy but the usual exponential rate of convergence of spectral methods is destroyed by the boundary singularities of the solution. In the sequel to this work, we will apply a variety of strategies including a change-of-coordinates and singular basis functions to recover spectral accuracy in spite of the boundary singularities. As preparation for this numerical study, we find explicit solutions to related problems to the two-dimensional PCP equation in a domain with a boundary that is an ellipse and the three-dimensional PCP equation in a cubic domain. We also analyze the boundary behavior of these solutions: all have complicated singularities with unbounded first derivatives. (C) 2021 Elsevier B.V. All rights reserved.
NSTL
Elsevier