Publications
Title : | Transparent boundary conditions for evolution equations in infinite periodic strips |
Year : | 2012 |
Type : | paper in peer-reviewed journal |
Authors : | J. Coatléven |
Abstract : | We consider the solution of a generic equation $\gamma\rho(\mathbf{x})\partial^p_tu(\mathbf{x},t)-\Delta u(\mathbf{x},t) +V(\mathbf{x})u(\mathbf{x},t) = f(\mathbf{x},t)$, $\mathbf{x} = (x,y)$, for $t>0$, $p=1,2$ in a domain $\Omega$ which is infinite in $x$ and bounded in $y$. We assume that $f(\cdot,t)$ is supported for all $t>0$ in $\Omega_0 = \{\mathbf{x} \in \Omega \; | \; -a_- < x < a_+\}$ and that $\rho(\mathbf{x})$ and $V(\mathbf{x})$ are x-periodic in $\Omega \setminus \Omega_0$. We consider the associated $\theta$-scheme in time, to obtain a semidiscretized problem. We then show how to obtain for each time step exact boundary conditions on the vertical segments, $\Gamma_0^- = \{\mathbf{x}\in \Omega\; | \; x=-a_-\}$ and $\Gamma_0^+ = \{\mathbf{x}\in \Omega \;| \; x=a_+\}$, that will enable us to find the solution on $\Omega_0 \cup \Gamma_0^+ \cup \Gamma_0^-$. Then the solution can be extended in $\Omega$ in a straightforward manner from the values on $\Gamma_0^-$ and $\Gamma_0^+$. The method is based on the solution of local problems on a single periodicity cell, solved during an initialization step. The exact boundary conditions as well as the extension operators can be obtained for each time step through elementary computations using the solution of these local cell problems. |
Themes : |
Artificial boundaries Transient waves |
Reference : | SIAM journal on Scientific Computing - vol. 34(3) (pp 1563-1583 ) |