34.时域有限差分3.docVIP

IV. Mur’s Absorbing Boundary Condition ? Introduction In radiation and scattering problems, the domain in which the field will be computed is assumed unbounded. Since no computers can store an unlimited amount of data, a certain method has to be designed for limiting the domain in which the field is computed. This is, certain kind of lattice truncation has to be implemented. Proper truncation of the lattice requires that the size of the truncated lattice must be large enough to fully contain the antennas and scatters, and any outgoing waves disappear at the truncated lattice boundary without reflection during the continuous time stepping of the algorithm. It is desirable to impose a boundary condition on the outer surface of the truncated mash such that the unbounded surrounding is modeled as accurately as possible. Boundary conditions of this type are termed the absorbing boundary conditions (abbreviated to ABC). A number of methods for setting up ABC have already be put forward and satisfactory results are acquired. Among them, the method supposed by Mur can gain the approximation of higher order. If the boundary surface is perfectly conducting, a “hard” boundary condition instead of “soft” absorbing boundary condition on this surface is applied, where tangential components of and normal components of vanish. ? Mur’s Absorbing Boundary Conditions For the sake of simplicity, only cubic cells in free space for a FDTD lattice structure at hand are considered. The wave equation governing a time-varying field component is (22) (23) We will assume, without loss of generality, that the lattice truncation surface is the plane x=0, and the lattice is located in the region . Let , where the constants are called the inverse velocity components, and the solution to wave equation possesses the form of (24) where g denotes a function of time and space. [proof] Let then and similarly, , Hence On t