Neumann–Neumann methods

In mathematics, Neumann–Neumann methods are domain decomposition preconditioners named so because they solve a Neumann problem on each subdomain on both sides of the interface between the subdomains. Just like all domain decomposition methods, so that the number of iterations does not grow with the number of subdomains, Neumann–Neumann methods require the solution of a coarse problem to provide global communication. The balancing domain decomposition is a Neumann–Neumann method with a special kind of coarse problem.

More specifically, consider a domain Template:Math, on which we wish to solve the Poisson equation

$$-\mathrm{\Delta }u=f,u{|}_{\partial \mathrm{\Omega }}=0$$

for some function f. Split the domain into two non-overlapping subdomains Template:Math and Template:Math with common boundary Template:Math and let Template:Math and Template:Math be the values of Template:Math in each subdomain. At the interface between the two subdomains, the two solutions must satisfy the matching conditions

$$u_1=u_2,{\partial }_{n_1}{u}_1={\partial }_{n_2}{u}_2$$

where $n_i$ is the unit normal vector to Γ in each subdomain.

An iterative method with iterations Template:Math for the approximation of each Template:Math (Template:Math) that satisfies the matching conditions is to first solve the Dirichlet problems

$$\begin{array}{cc}-& \mathrm{\Delta }{u}_i^{(k)}=f_i\text{in}{\mathrm{\Omega }}_i,\\ & {u_i^{(k)}|}_{\partial \mathrm{\Omega }}=0,{u_i^{(k)}|}_{\mathrm{\Gamma }}={\lambda }^{(k)}\end{array}$$

for some function Template:Math on Template:Math, where Template:Math is any inexpensive initial guess. We then solve the two Neumann problems

$$\begin{array}{cc}-& \mathrm{\Delta }{\psi }_i^{(k)}=0\text{in}{\mathrm{\Omega }}_i,\\ & {{\psi }_i^{(k)}|}_{\partial \mathrm{\Omega }}=0,{{\partial }_{n_i}{\psi }_i^{(k)}|}_{\mathrm{\Gamma }}=\omega ({\partial }_{n_1}{u}_1^{(k)}+{\partial }_{n_2}{u}_2^{(k)}).\end{array}$$

We then obtain the next iterate by setting

$${\lambda }^{(k+1)}={\lambda }^{(k)}-\omega ({\theta }_1{\psi }_1^{(k)}+{\theta }_2{\psi }_2^{(k)})\text{on}\mathrm{\Gamma }$$

for some parameters Template:Math, Template:Math and Template:Math.

This procedure can be viewed as a Richardson iteration for the iterative solution of the equations arising from the Schur complement method.

This continuous iteration can be discretized by the finite element method and then solved—in parallel—on a computer. The extension to more subdomains is straightforward, but using this method as stated as a preconditioner for the Schur complement system is not scalable with the number of subdomains; hence the need for a global coarse solve.

See also

• Neumann–Dirichlet method