Type I: Voltage-Centered Mesh

$$Y_{\rho^{+},i,j} = Y_{\rho^{-},i,j} = Y_{\theta^{+},i,j} = Y_{\theta^{-},i,j} = \frac{\Delta_{\rho}c_{u,i,j}}{2T}\hspace{0.3in}Y_{c,i,j} = 0$$ (4.97)

$$Z_{c,i+\frac{1}{2},j} = \frac{\Delta_{\rho}l_{\rho,i,j}}{T}-\frac{2T}{\Delta_{\rho}c_{u,i,j}} + \frac{\Delta_{\rho}l_{\rho,i+1,j}}{T}-\frac{2T}{\Delta_{\rho}c_{u,i+1,j}}$$ (4.98)

$$Z_{c,i,j+\frac{1}{2}} = \frac{\Delta_{\theta}l_{\theta,i,j}}{T}-\frac{2T}{\Delta_{\rho}c_... ...+ \frac{\Delta_{\theta}l_{\theta,i,j+1}}{T}-\frac{2T}{\Delta_{\rho}c_{u,i,j+1}}$$ (4.99)

The stability constraints (which follow from the requirement of positivity of $Z_{c}$ everywhere) are

$$\frac{\Delta_{\rho}}{T}\geq \max_{i,j}\sqrt{\frac{2}{l_{i,j}c_{i,... ...{T}\geq\max_{i,j}\left(\frac{1}{\rho_{i}}\sqrt{\frac{2}{l_{i,j}c_{i,j}}}\right)$$ (4.100)

There is thus a dependence on $\rho$ in the second condition (relating the angular spacing $\Delta_{\theta}$ to $T$), which we expect, since the spacing between the junctions at a given radius now varies linearly with the radius. Stability bounds are, for a radial mesh, necessarily more severe than in the rectilinear case, due to this variation in grid spacing.