The Waveguide Mesh in Radial Coordinates

We will look at waveguide meshes in general curvilinear coordinates in §4.8, but radial coordinates are an important special case, especially for musical instrument physical modelling applications (considering how many instruments exhibit some form of radial symmetry).

In terms of radial coordinates , where

the parallel-plate system (4.58) becomes

where we define radial and angular current densities by

and the effective radial material parameters and sources by

is a scaling coefficient which we will set, in anticipation of discretization, equal to , the ratio of the grid spacings in the and directions. We will allow these spacings to be, in general, different.

It is evident that system (4.82) has a form similar to its counterpart in rectilinear coordinates, apart from the extra factor of in (4.82c). The chief difference is that we now have different effective inductances and in the two coordinate directions, but, as we shall see, this anisotropy is easily taken care of (indeed, we could have defined anisotropic inductances and in the rectilinear case without greatly complicating matters). We can see immediately that when centered differences are applied to system (4.82), we will be able to operate on an interleaved grid in the coordinates. A version of Yee's algorithm in arbitrary curvilinear coordinates first appeared in [91], and is also discussed in [209]. The interleaved grid, viewed in rectilinear coordinates, is shown in Figure 4.28, where, as before, the dependent variable to be calculated at a particular grid point is indicated next to the point. Grey and white coloring of points indicates operation at alternating time steps.

Centered differencing yields a scheme nearly identical to (4.59), with, again, the difference that the inductance has a directional character. We can thus proceed directly to the waveguide mesh, and, furthermore, can use the same indexing as in the rectilinear case; now, the grid indices will refer to points . Due to the interleaved nature of the resulting difference approximations, we will have series junctions at locations and for , and integer (with associated junction currents and ) and parallel junctions at locations where we will calculate junction voltages , for , and integer (we return to the central grid point at later in this section). The computational molecule of the mesh is shown in Figure 4.29.
Referring to Figure 4.30, which gives the immittance nomenclature in the waveguide network, and where in addition we have the junction immittances defined by

for and integer, we can perform an analysis similar to the rectilinear case in order to determine that we must have

(The junction admittance will be dealt with shortly.) The source waves should be chosen as

where we may of course use the dual type of wave in regions where the loss parameters become small, as discussed in §4.3.7. Just as in the rectilinear case, these conditions define a family of waveguide networks which solve the radial transmission line equations. We here provide the impedance settings for voltage- and current-centered meshes, as well as stability bounds.

- Type I: Voltage-Centered Mesh
- Type II: Current-Centered Mesh
- Central Gridpoint in a Radial Mesh
- Simulation: Circular Region with Varying Inductance