Other Waveguide Networks for Timoshenko's System

The difference between the two forms presented here is only in the coupling between the two lines, each type of which gives rise to a different type of offset sampling. For example, in the configuration shown in Figure 5.10, we calculate and at the same time instant, and at the same grid location, whereas in Figure 5.11, the updating of with respect to is staggered in both time and space.

There are many ways of setting the network immittances so as to approximate the Timoshenko system. Since both networks are simply just pairs of coupled (1+1)D transmission lines, the derivation of the constraints on the junction immittances is familiar. Because the two networks are so similar, we will state these constraints only for the network shown in Figure 5.11:

We use the bar notation again to indicate an approximation to at least second-order accuracy in the spatial step , and the tilde to distinguish between the immittances in the two transmission lines. Before we begin setting the immittances values, however, we must comment on the coupling waveguides which run between the parallel junctions which calculate and ; each such waveguide has a non-reciprocal character (as evidenced by the sign inversion in only one of the two signal paths). We have left the admittance of this coupling waveguide independent of position. The reason for this is seen from the difference equations in the junction variables. For example, in the configuration of Figure 5.11, we have, after following through familiar manipulations,

which will approximate (5.17b) only if . The same statement is true of the other configuration as well (where the admittance of each coupling waveguide is set to ).

It is convenient to set the immittances such that the self-loops at the parallel junctions (recall the transmission line) can be dropped entirely from the calculation; in fact, in this lossless case, we can drop the explicit calculation of and altogether if we wish, since when the self-loops at these junctions disappear, we are left with simple ``through'' junctions which do not scatter (recall Figure 4.14). Eventually, though, we may wish to have access to this variable because it is indeed the time derivative of , the primary dependent variable of interest. We thus set

and

It is quite interesting to examine the positivity conditions on the self-loop admittances. For example, will be positive if, for all integer , we have

which, using , is equivalent to

It is easy to see that, in contrast to the stability conditions we derived for the various transmission lines, for the ideal beam, and for Timoshenko's system in §5.2.1, for a passive realization there is now a

Assuming that this condition is met, we then must have

A positivity requirement on yields a similar set of conditions,

and |

In the limit as becomes small, both conditions on the spatial step size approach conventional CFL-type bounds (recall the wave speeds implied by Timoshenko's system, mentioned in the beginning of this section); both must be satisfied for the waveguide network to be *concretely passive*, though they may be stable in the Von Neumann sense (see Appendix A). It would appear that the special character of these passivity conditions on the scheme are the direct result of the approximation of the memoryless asymmetric line coupling by a reactive waveguide coupling.