An optimization-based domain decomposition method for parabolic equations (Q2493734)

An optimization-based nonoverlapping domain decomposition method is presented for the 2D heat conduction equation of the form: \[ u_{t} - \triangle u = f, \quad \text{in}\quad \Omega \times (0,T), \] subject to the respective boundary and initial conditions \( u = 0\) on \(\Gamma \times (0,T),\) and \(u| _{t=0} = u_0\) in \(\Omega\). This method was proposed and analyzed earlier for the Poisson equation by \textit{M. D. Gunzburger, J. S. Peterson} and \textit{H. Kwon} [Comput. Math. Appl. 37, No. 10, 77--93 (1999; Zbl 0941.65123)]. The crux of the method is a constrained minimization problem for which the objective functional measures the jump in the dependent variables across the common boundaries between subdomains. An optimality system for the optimal solution is derived and used to define a gradient method. Convergence results are obtained for the gradient method. The present method is implemented to a test problem, and the results for rate of convergence are presented.