An optimized Schwarz waveform relaxation method for the unsteady convection diffusion equation in two dimensions (Q5916120)

The author is interested in solving the scalar convection diffusion equation in \(\Omega= \mathbb{R}^2\) \[ \begin{cases}{\partial w\over\partial t}+\vec{\mathbf b}\cdot\nabla w- \nu\Delta w+ cw= f\quad &\text{in }\Omega\times \mathbb{R}^+,\\ w(\cdot,\cdot,0)= w_0\quad &\text{in }\Omega,\end{cases} \] where \(\vec{\mathbf b}= (a,b)\) is a constant velocity field, \(\nu> 0\) is the viscosity and \(c> 0\). The purpose of the paper is to study the algorithm \[ \begin{cases}{\mathcal U}u^{n+1}= f\quad &\text{in }\Omega^-\times \mathbb{R}^+,\\ u^{n+1}(\cdot,\cdot,0)= w_0\quad &\text{in }\Omega^-,\\ {\mathcal B}^- u^{n+1}={\mathcal B}^- v^n\quad &\text{in }\Gamma_L\times\mathbb{R}^+,\end{cases}\;\begin{cases}{\mathcal L}v^{n+1}= f\quad &\text{in }\Omega^+\times\mathbb{R}^+,\\ v^{n+1}(\cdot,\cdot,0)= w_0\quad &\text{in }\Omega^+,\\{\mathcal B}^+ v^{n+1}={\mathcal B}^+ v^n\quad &\text{in }\Gamma_0\times\mathbb{R}^+,\end{cases} \] where \({\mathcal L}\equiv\partial_t+ \vec{\mathbf b}\cdot\nabla+ \nu\Delta+ c\), \(\Omega^-= (-\infty, L)\times \mathbb{R}\,(L>0)\), \(\Omega^+= (0,\infty)\times \mathbb{R}\), \(\Gamma_L= \{y\in\mathbb{R},x= L\}\), \(\Gamma_0= \{y\in\mathbb{R}, x= 0\}\). Conditions are given under which the algorithm converges. Numerical examples illustrate the efficiency of the method.