next up previous
Next: Progressive Grid Density Doubling Up: Doubled Grid Density Across Previous: Corners

Simulation

To demonstrate the behavior of the interface discussed above, we will simulate the parallel-plate system (4.58), assuming no losses or sources, over a square region with side-length 1, with short-circuited boundary conditions (i.e., $ u=0$). Over a central square region with side lengths of $ \frac{1}{2}$ (bounded by the white square in Figure 4.37), we use a mesh of doubled density with respect to the outer region, where a simple rectilinear mesh is in place. The grid spacing in the outer region is $ \Delta = 0.01$. We set the capacitance per unit length $ c=1$ everywhere, and the inductance per unit length $ l$ is equal to 1 in the outer region. In the inner region, it also takes the value 1, except in the interior of the black circle in Figure 4.37(radius 0.1, center at $ x=0.6$, $ y=0.5$), where it rises, in a 2D raised cosine distribution, to a maximal value of $ l=11$.

Figure 4.37: Simulation of a parallel-plate system, with an inductive feature surrounded by a region of doubled grid density.
\begin{figure}\begin{center}
\begin{picture}(390,440)
\par\put(0,-1){\epsfig{fi...
...38.5){2}{\circle{32.8}}}
\end{picture} \end{center} \vspace{0.2in}
\end{figure}

The initial voltage distribution takes the form of a 2D raised cosine over a circular region of radius 0.07, centered at coordinates $ x=0.4$ and $ y=0.5$. The initial voltage takes a maximum of 1, and is zero everywhere outside this circle.

The voltage distribution is shown at three successive time instants, with light and dark coloring indicating areas of positive and negative plate voltage, respectively; the plots have been normalized and interpolated for better plotting results. At time $ t=0.15$, the voltage distribution has spread and begun to scatter from the inductive region inside the black circle. Note that the local wave speed decreases in the interior of the circle. By $ t=0.3$, it has progressed through the interface, with little visible numerical reflection, and by $ t=0.45$, the outward moving voltage distribution has reflected with inversion from the boundary at $ x=0$.


next up previous
Next: Progressive Grid Density Doubling Up: Doubled Grid Density Across Previous: Corners
Stefan Bilbao 2002-01-22