next up previous
Next: Higher-order Accuracy Revisited Up: Incorporating the DWN into Previous: Alternative MDKC for the

Alternative MDKC for (2+1)D Parallel-plate System

The same idea extends simply to the (2+1)D parallel-plate equations (3.67). In this case, we again look at the lossless source-free system, scaled by a factor of $ \Delta$, and can develop an MDKC along the lines of Figure 4.47(b), with $ u$ now treated as a voltage instead of a current as in Figure 3.17. We remind the reader that we use the notation $ D_{k}$, $ k=1,\hdots, 5$ in circuit diagrams to indicate directional derivatives in the directions $ t_{k}$ defined by the coordinate transformation (3.22). Also, we have written $ i_{1}$ and $ i_{2}$ for $ i_{x}$ and $ i_{y}$, in keeping with the WDF literature [62]. The alternative MDKC is shown in Figure 4.49, and the element values are given by

$\displaystyle L_{1} = L_{2} = \Delta\left(v_{0}l-r_{0}\right)\hspace{0.3in} C_{...
... \Delta\left(v_{0}c-\frac{2}{r_{0}}\right)\hspace{0.3in} L_{0} = \Delta r_{0}/2$    

$ r_{0}>0$ is, as before, a free parameter which can be scaled to achieve an optimal space step/time step ratio.

Figure 4.49: Alternative MDKC for the (2+1)D parallel-plate system in rectangular coordinates.
% graphpaper(0,0)(350,290...
\put(25,-15){\small {$A_{2}'$}}
\end{picture} \end{center} \end{figure}

The two two-ports with terminals $ A_{1}$, $ A_{1}'$, $ B$, $ B'$ and $ A_{2}$, $ A_{2}'$, $ B$, $ B'$ are both linear and shift-invariant, and their right-hand pairs of terminals are both connected in parallel with the capacitor. The hybrid matrices for these two-ports are

$\displaystyle {\bf K}_{1}(s_{1},s_{3})$ $\displaystyle = \frac{\Delta}{2}\begin{bmatrix}r_{0}\left(s_{1}+s_{3}\right)&s_{1}-s_{3}\\ s_{1}-s_{3}&\frac{1}{r_{0}}\left(s_{1}+s_{3}\right)\\ \end{bmatrix}$   $\displaystyle \mbox{{\rm Two-port$\hspace{0.1in}$\ $A_{1}A_{1}'BB'$}}$    
$\displaystyle {\bf K}_{2}(s_{2},s_{4})$ $\displaystyle = \frac{\Delta}{2}\begin{bmatrix}r_{0}\left(s_{2}+s_{4}\right)&s_{2}-s_{4}\\ s_{2}-s_{4}&\frac{1}{r_{0}}\left(s_{2}+s_{4}\right)\\ \end{bmatrix}$   $\displaystyle \mbox{{\rm Two-port$\hspace{0.1in}$\ $A_{2}A_{2}'BB'$}}$    

and under the spectral mappings
\begin{subequations}\begin{align}s_{1} &\rightarrow \frac{1}{T_{1}}\frac{(1-z_{1...
...z_{4}^{-1})(1+z_{2}^{-1})}{1+z_{2}^{-1}z_{4}^{-1}} \end{align}\end{subequations}

with $ T_{1}=T_{2}=T_{3}=T_{4}=\Delta/2$, we again have a decomposition of the discrete hybrid two-ports into pairs of multidimensional unit elements (impedance $ r_{0}$) connected in series/parallel, i.e.
$\displaystyle {\bf K}_{1}(z_{1}^{-1},z_{3}^{-1})$ $\displaystyle =$ $\displaystyle {\bf K}_{ue}(z_{1}^{-1},z_{3}^{-1})+{\bf K}_{ue}(-z_{3}^{-1},-z_{1}^{-1})$  
$\displaystyle {\bf K}_{2}(z_{2}^{-1},z_{4}^{-1})$ $\displaystyle =$ $\displaystyle {\bf K}_{ue}(z_{2}^{-1},z_{4}^{-1})+{\bf K}_{ue}(-z_{4}^{-1},-z_{2}^{-1})$  

The MDWD network, corresponding to ``checkerboard''-type sampling, is as shown in Figure 4.50, and where $ {\bf T_{1}},\hdots, {\bf T_{4}}$ are the shifts in the coordinate directions $ t_{1},\hdots, t_{4}$, and the port resistances are given by

$\displaystyle R_{1} = R_{2} = 2\left(v_{0}l-r_{0}\right)\hspace{0.3in}R_{3} = \frac{1}{2\left(v_{0}c-2/r_{0}\right)}\hspace{0.3in}R_{0} = r_{0}$    

Figure 4.50: MD network equivalent to type III DWN for the lossless source-free (2+1)D parallel-plate system under an alternative spectral mapping.
% graphpaper(0,0)(500,180...
...ut(488,66){\footnotesize {${\bf T}$}}
\end{picture} \end{center}\par\end{figure}

For the one-port inductors and the capacitor, we have again applied the trapezoid rule, with a step-size of $ T_{5} = \Delta$ (implying a delay of $ T$, written as $ {\bf T}$ in Figure 4.50).

When spatial dependence is expanded out, the signal flow graph is identical to the interleaved type III DWN for the parallel-plate equations shown in Figure 4.21, without the loss/source ports. As in (1+1)D, loss and source elements can be reintroduced into the MDKC of Figure 4.49 without difficulty.

It should also be possible to derive DWNs from MDKCs for the same system on alternative grids, such as the hexagonal and triangular grids mentioned in §4.6.1, though we have not investigated this in any detail. In this case, one would presumably begin from the MDKC under the appropriate coordinate transformation. One of these coordinate transformations (which generates a hexagonal grid under uniform sampling in the new coordinates), was discussed briefly in §3.3.3. The full MDKC for the parallel-plate equations in these coordinates is given in [62].

next up previous
Next: Higher-order Accuracy Revisited Up: Incorporating the DWN into Previous: Alternative MDKC for the
Stefan Bilbao 2002-01-22