A new set of consistent boundary conditions for Yee scheme approximations of
wave equations in one and two space dimensions are developed and analyzed.
We show how the classical staircase boundary conditions for hard reflections
or, in the electromagnetic case, conducting surfaces in many cases give O(1)
errors in 2D.
The proposed conditions keep the structure of the Yee scheme and are thus
well suited for high performance computing. The higher accuracy is achieved
by modifying the coefficients in the difference stencils near the boundary.
We study stability and convergence and we present numerical examples. Some
preliminary results for the threee dimensional case are included.
|