Transparent boundary conditions for evolution equations in infinite periodic strips

2012
Type de publication :
Article (revues avec comité de lecture)
Journal :
SIAM journal on Scientific Computing, vol. 34(3), pp. 1563-1583
HAL :
hal-00969313
Résumé :
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.
BibTeX :
@article{Coa-2012-1,
    author={Julien Coatléven },
    title={Transparent boundary conditions for evolution equations in 
           infinite periodic strips },
    doi={10.1137/110838030 },
    journal={SIAM journal on Scientific Computing },
    year={2012 },
    volume={34(3) },
    pages={1563--1583},
}