Mathematical programming techniques to solve biharmonic problems by a recursive projection algorithm (Q2641093)

The biharmonic equation \(\Delta^ 2u=0\) is considered on a star-shaped domain \(\Omega\) in the plane, with \(u=f\) and \(\partial_ nu=g\) prescribed on the boundary \(\partial \Omega\). The numerical method consists of global approximation by biharmonic polynomials. A finite number of points \(\theta_ j\) is selected on \(\partial \Omega\), and a recursive projection algorithm is used to minimize the sum \(\sum_{j}\{| u(\theta_ j)-f(\theta_ j)|^ 2+| \partial_ nu(\theta_ j)-g(\theta_ j)|^ 2+| \partial_ tu(\theta_ j)-\partial_ tf(\theta_ j)|^ 2\}.\) Here, \(\partial_ n\) denotes the normal derivative on \(\partial \Omega\), and \(\partial_ t\) the tangential derivative. Computations are presented with varying numbers of basis elements and evaluation points \(\theta_ j\).