Next: Simulation
Up: Doubled Grid Density Across
Previous: Doubled Grid Density Across
If we are interested in using a grid of doubled density over a particular region of the problem domain (in order to surround a particular feature or an irregular part of the boundary), then we are faced with corners, and must develop special scattering junctions for them. An example of an irregular partitioning of the problem domain into two regions,
and
, is shown in Figure 4.36(a).
Figure 4.36: (a) A particular grid arrangement between a rectilinear mesh
and one of doubled grid density
, (b) a scattering junction at a corner point (labelled
), and (c) a scattering junction at another type of corner point (labelled
).
![\begin{figure}\begin{center}
\begin{picture}(540,330)
% graphpaper(0,0)(540,330...
...\put(425,217){\tiny {p}}
\end{picture} \end{center} \vspace{0.2in}
\end{figure}](img2005.png) |
Boundary points (labelled
) were treated previously, and it was found that waveguides connected to points
which lie along the boundary must have their admittances set to one-half what they would be in the interior of region
(i.e., to
, where
is the inductance at the center of the particular boundary waveguide).
There are two types of corners which can arise in an irregular domain decomposition of this kind: those which are concave with respect to region
(the grid point at such a corner is labelled
in Figure 4.36(a)), and those which are concave with respect to region
(labelled
). There are obviously four possible orientations for each type of corner, though, by symmetry, we need only treat the type shown. Because the admittances of the boundary waveguides (represented by thick lines connecting boundary points) are now prescribed, then for a corner junction the linking admittances are fixed; only the self-loop admittance may be varied. Scattering junctions corresponding to corners of type
and
are shown in Figure 4.36(b) and (c).
Suppose we have a corner point of type
located at coordinates
. Then, from the results previously given in this section, the admittances of the five waveguides connecting this corner to its neighbors will be
and the difference equation relating the corner junction voltage
to those of its neighbors is
where we have yet not specified
or
.
A Taylor expansion about
gives, in terms of the continuous variable
and neglecting higher-order terms,
We can thus conclude that the differencing occurring at the corner junction is not consistent with the lossless source-free parallel-plate system, regardless of our choice of the self-loop admittance. The coefficients of the spatial first derivative terms are incorrect, and there is an extra mixed-derivative term; neither vanishes in the limit as
(although if
is constant in a neighborhood surrounding the corner, then the first derivative terms vanish).
The corner junction is, however, still lossless, as are all junctions in this waveguide network, and hence a simulation of the parallel-plate system using such a mesh will be stable, regardless of inconsistencies at the corners. It is possible to argue, loosely speaking, that if the number of corners in the interface does not grow as the grid spacing is decreased (an example of this would be an enclosed rectangular doubled density grid, for which the number of corners will be four, independently of the grid spacing), then the error at the corners will become negligible. Although we will make no attempt to prove this, the simulation which we will present later in this section concurs readily with this assertion; indeed, in all tests we have run, any anomalous scattering at the corners is certainly far less important than the first-order scattering error (i.e., numerical reflection) along the interface itself.
In the interest, however, of making any scattering error at the corner junctions as small as possible, we should set, at a corner with coordinates
(in order that the wave speed at the corner is correct, in a gross sense)
where
and
are the inductance and capacitance at a corner point
. The stability requirement at such a corner is then
which is, like the stability condition at boundary points
, only marginally worse than the CFL bound, and mitigated by the somewhat worse bound to be found in region
(due to the decreased inter-junction spacing). See the discussion earlier in this section.
A similar argument follows for points
, and we also ideally have
constant in the neighborhood of such points. The setting of the self-loop admittance
should be
The resulting stability bound is
Next: Simulation
Up: Doubled Grid Density Across
Previous: Doubled Grid Density Across
Stefan Bilbao
2002-01-22