next up previous
Next: Numerical Phase Velocity and Up: Finite Difference Interpretation Previous: Finite Difference Interpretation


MDWD Networks as Multi-step Schemes

Recall that the discretization step discussed in §3.5.3 consisted of the application of a spectral transformation of the form

$\displaystyle s_{j}\rightarrow\frac{2}{T_{j}}\frac{1-z_{j}^{-1}}{1+z_{j}^{-1}}$    

where $ s_{j}$ is the frequency domain transform variable corresponding to any MD-causal coordinate $ t_{j}$, and $ z_{j}^{-1}$ is the frequency domain unit shift in the same direction.

For spatially inhomogeneous problems, this spectral mapping is equivalent to the application of the trapezoid rule in direction $ t_{j}$. We can thus write, using operator notation,

$\displaystyle \frac{\partial}{\partial t_{j}}\rightarrow \frac{2}{T_{j}}\left(1+\delta_{j}\right)^{-1}\left(1-\delta_{j}\right)$ (3.69)

where $ \delta_{j}$ is a shift operator defined by

$\displaystyle \delta_{j}p(t_{1},\hdots, t_{j},\hdots, t_{k}) = p(t_{1},\hdots, t_{j}-T_{j},\hdots, t_{k})$    

when applied to any continuous function $ p({\bf t})$. Consider again the lossless (1+1)D transmission line equations
$\displaystyle \begin{eqnarray}l\frac{\partial i}{\partial t}+\frac{\partial u}{...
...\ c\frac{\partial u}{\partial t}+\frac{\partial i}{\partial x}=0 \end{eqnarray}$ (3.70a)

which can be written as
$\displaystyle \begin{eqnarray}\left(v_{0}l-r_{0}\right)\frac{\partial i_{1}}{\p...
...qrt{2}}\frac{\partial}{\partial t_{2}}\left(i_{2}-i_{1}\right)=0 \end{eqnarray}$ (3.71a)

under the application of coordinate transformation (3.19) and using scaled variables $ i_{1} = i$ and $ i_{2} = u/r_{0}$, as well as the scaled time variable $ t' = v_{0}t$. Under the substitution of (3.73), for $ j=1,2$, and using the generalized trapezoid rule in time, defined by

$\displaystyle \frac{\partial}{\partial t'}\rightarrow \frac{2}{T'}\left(1+\delta_{t'}\right)^{-1}\left(1-\delta_{t'}\right)$    

where $ \delta_{t'}$ is a shift in the scaled time direction $ t'$ of duration $ T'$, we get


$\displaystyle \left(v_{0}l-r_{0}\right)\frac{2}{T'}\left(1+\delta_{t'}\right)^{-1}\left(1-\delta_{t'}\right)i_{1}$ $\displaystyle +$ $\displaystyle \frac{\!\!\sqrt{2}r_{0}}{T_{1}}\left(1+\delta_{1}\right)^{-1}\left(1-\delta_{1}\right)\left(i_{1}+i_{2}\right)\notag$  
  $\displaystyle +$ $\displaystyle \frac{\!\!\sqrt{2}r_{0}}{T_{2}}\left(1+\delta_{2}\right)^{-1}\left(1-\delta_{2}\right)\left(i_{1}-i_{2}\right)=0$  
$\displaystyle \left(v_{0}cr_{0}^{2}-r_{0}\right)\frac{2}{T'}\left(1+\delta_{t'}\right)^{-1}\left(1-\delta_{t'}\right)i_{2}$ $\displaystyle +$ $\displaystyle \frac{\!\!\sqrt{2}r_{0}}{T_{1}}\left(1+\delta_{1}\right)^{-1}\left(1-\delta_{1}\right)\left(i_{2}+i_{1}\right)\notag$  
  $\displaystyle +$ $\displaystyle \frac{\!\!\sqrt{2}r_{0}}{T_{2}}\left(1+\delta_{2}\right)^{-1}\left(1-\delta_{2}\right)\left(i_{2}-i_{1}\right)=0$  

to second order in $ \Delta$. This can be rewritten as

$\displaystyle \begin{eqnarray}\left(1+\delta_{1}\right)\left(1+\delta_{2}\right...
...(1+\delta_{t'}\right)\left(\delta_{2}-\delta_{1}\right)i_{1} = 0 \end{eqnarray}$ (3.72a)

where we have used $ T'= v_{0}T = \Delta$, $ T_{1}=T_{2}=\sqrt{2}\Delta$, the fact that $ \delta_{1}\delta_{2} = \delta^{2}_{t'}$, and also the definitions of the port resistances of the MDWD network of Figure 3.14(b), given in (3.64)% latex2html id marker 81317
\setcounter{footnote}{2}\fnsymbol{footnote}. Upon replacing the quantities $ i_{1}$ and $ i_{2}$ by their respective grid functions $ I_{1,i}(n)$ and $ I_{2,i}(n)$, which take on values for $ n$ and $ i$ integer, (3.77a) and (3.77b) define recursions on a regular grid, of spacing $ \Delta$. (3.77a) can be written as
$\displaystyle \alpha_{i}I_{1,i}(n)$ $\displaystyle +$ $\displaystyle \beta_{i+1}I_{1,i+1}(n-1) + \beta_{i-1}I_{1,i-1}(n-1)+\gamma_{i}I_{1,i}(n-1)\notag$ (3.73)
  $\displaystyle -$ $\displaystyle I_{2,i-1}(n-1)+I_{2,i+1}(n-1)\notag$ (3.74)
  $\displaystyle -$ $\displaystyle \beta_{i+1}I_{1,i+1}(n-2) - \beta_{i-1}I_{1,i-1}(n-2)-\gamma_{i}I_{1,i}(n-2)\notag$ (3.75)
  $\displaystyle -$ $\displaystyle I_{2,i-1}(n-2)+I_{2,i+1}(n-2)\notag$ (3.76)
  $\displaystyle -$ $\displaystyle \alpha_{i}I_{1,i}(n-3)=0$ (3.77)

with

$\displaystyle \alpha_{i}\triangleq 1+\frac{R_{1}(i\Delta)}{R_{0}}\hspace{0.5in}...
...Delta)}{R_{0}}\hspace{0.5in}\gamma_{i}\triangleq 1-\frac{R_{1}(i\Delta)}{R_{0}}$ (3.78)

The recursion corresponding to (3.77b) is very similar, under the interchange of $ I_{1}$ and $ I_{2}$. Note that if $ l$ and $ c$ are constants, and if the difference scheme is operating at the CFL bound (so that $ R_{1} = R_{2} = 0$, then (3.78) can be simplified to

$\displaystyle I_{1,i}(n)-I_{1,i}(n-2) + I_{2,i+1}(n-1)-I_{2,i-1}(n-1) = 0$ (3.79)

which is a simple centered difference approximation to (3.74a) and which we will see again in the waveguide mesh context in §4.3.2. Unlike the case of the mesh however, away from the passivity bound we have a multi-step scheme [176] which involves three steps of ``look-back'' in order to update a grid variable at a particular location. The introduction of wave variables, then, can be considered to be a means of expanding the state of the system so that using the new state, the recursion (now in the form of the MDWDF of Figure 3.14) requires access only to wave quantities at the time step immediately preceding the current one.

In order to generate a scheme which operates on alternating interleaved grids (called offset sampling in [61]), it is possible to use a doubled time step of $ T' = 2\Delta$ in order to implement the generalized trapezoid rule applied to the time derivatives in (3.75a) and (3.75b), i.e.,

$\displaystyle \frac{\partial}{\partial t'}\rightarrow \frac{2}{T'}\left(1-\delt...
...rac{1}{\Delta}\left(1-\delta_{t'}^{2}\right)^{-1}\left(1+\delta_{t'}^{2}\right)$    

in which case we get, as an approximation to (3.74a),
$\displaystyle \alpha_{i}I_{i,1}(n)$ $\displaystyle +$ $\displaystyle \beta_{i+1}I_{1,i+1}(n-1)+\beta_{i+1}I_{1,i-1}(n-1)+ I_{2,i+1}(n-1)-I_{2,i-1}(n-1)\notag$ (3.80)
  $\displaystyle -$ $\displaystyle \beta_{i+1}I_{1,i+1}(n-3)-\beta_{i+1}I_{1,i-1}(n-3)+ I_{2,i+1}(n-3)-I_{2,i-1}(n-3)\notag$ (3.81)
  $\displaystyle -$ $\displaystyle \alpha_{i}I_{1,i}(n-4) = 0$ (3.82)

where $ \alpha$ and $ \beta$ are defined as per (3.79), but where $ R_{1}$ is now equal to $ \frac{1}{\Delta}\left(v_{0}l-r_{0}\right)$. This form also reduces to simple centered differences when $ l$ and $ c$ are constant, and when we are operating at the CFL bound.

The computational stencils corresponding to the two different schemes are shown in Figure 3.20; the top black dot in either picture represents the location of the grid variable currently being updated (either $ I_{1}$ or $ I_{2}$), and the other dots cover the discrete region of influence of the difference scheme. Notice in particular that each scheme has a width of only three grid points, corresponding to nearest-neighbor-only updating. Also, because these are multi-step methods, one might expect that we will have to take special care when initializing the scheme; we discuss this issue in §3.10. For the offset scheme of Figure 3.20(b), the stencil can be shifted one step to the left or right without any overlapping; thus such a scheme can subdivided into two mutually exclusive subschemes (operating only for $ n+i$ always even or always odd), one of which may be dropped from the calculating scheme entirely. This behavior appears in many of the difference schemes which we will come across subsequently; we will pay particular attention to such schemes during a spectral analysis of finite difference schemes in Appendix A.

Figure 3.20: Computational stencils of the equivalent multi-step schemes of MDWDFs for the (1+1)D transmission line equations-- (a) scheme (3.78) and (b) ``offset'' scheme (3.81).
\begin{figure}\begin{center}
\begin{picture}(400,230)
% graphpaper(0,0)(400,230...
...(a)}
\put(315,-40){(b)}
\end{picture} \end{center} \vspace{0.3in}
\end{figure}

One of the interesting features of the MDKC representation of a set of PDEs is that the same circuit can give rise to an entire family of MDWD networks, or, in other words, of difference methods, all of which are consistent with the original set of PDEs. In the case of the MDKC for the transmission line equations derived previously, although we have defined the directions of the various inductors (along which we will be integrating), at the circuit stage we have not as yet specified any spectral mapping which will determine the type of differencing to be applied. Any passivity-preserving mapping which is correct in the low frequency limit will give rise to a passive, consistent MDWD algorithm. We will examine the important implications of a more exotic type of mapping in §4.10, but it is also interesting to note that we can apply the trapezoid rule using different step sizes for all the reactive elements. The constraint on our choices of these step sizes is that all shifting operations refer, ultimately, to another grid point (for computability).

Finally, we note that in general, the determination of stability for a multi-step scheme can be quite difficult; even in the constant coefficient case, it will in general be necessary to perform Von Neumann analysis [176] (see Appendix A for such an analysis applied to difference schemes for the wave equation in (2+1)D and (3+1)D), which can be quite formidable. Here, however, we are ensured stability through the passivity condition on the network.


next up previous
Next: Numerical Phase Velocity and Up: Finite Difference Interpretation Previous: Finite Difference Interpretation
Stefan Bilbao 2002-01-22