Symmetric Hyperbolic Systems

In a distributed problem, the dependent variables are functions not only of time , but also of location within an -dimensional spatial domain
, with coordinates
. Such a problem is referred to as an D problem in the WDF literature [131]. Problems without spatial dependence will be called *lumped* problems. If the equations which define the problem include differential operators, we are faced with solving a set of partial differential equations (PDEs).

A particularly important family of PDE systems are the *symmetric hyperbolic* [74,82] systems of the form

Here, , the state, is a -element column vector defined over coordinates and . and , , are real symmetric matrices

Though it is possible to extend this definition to include cases where the matrices
may depend on , or even (in which case system (3.1) is nonlinear), this simpler form describes a wide variety of physical systems, from electromagnetics to string, membrane, beam, plate, shell, and elastic solid dynamics, to transmission line systems, to linear acoustics, etc. Symmetric hyperbolic systems are important because they form a subclass of *strongly hyperbolic* systems, for which the initial-value problem is well-posed [176]. Roughly speaking, to say that a system is well-posed is to say that the growth of its solution is bounded in a well-defined way; growth in an norm cannot be faster than exponential. This concept is elaborated in detail in [82,176]. We can examine this growth in the present case as follows.

First, assume that the problem is defined over an unbounded spatial domain , so that we can drop any consideration of boundary conditions, and also that the forcing function . We now take the inner product of (the transpose of ) with (3.1) to get

where we have replaced by its symmetric part . Due to the symmetry of and the , we can then write

Now, integrate (3.3) over , to get

where is the D differential volume element. The expression is easily seen to be the divergence of a vector field, and by the Divergence Theorem [174], the integral of this quantity can be replaced by a surface integral over the problem boundary--because we have assumed no boundary, this integral vanishes, and we are left with

The quantity

can be interpreted as the total energy of system (3.1) at time . Note that due to the positivity requirement on , it is a positive definite function of the state, . If is positive semi-definite, then we must have, from (3.5), that

which implies that

In other words, the energy of the system must decrease as time progresses.

In the MD circuit models that we will discuss, what we will be doing, in essence, is dividing this energy up among various reactive MD circuit elements. We will elaborate on this in the sections on the (1+1)D transmission line and (2+1)D parallel-plate system. The passivity condition is essentially equivalent to (3.7). Also, the symmetric nature of the systems will be reflected, in the circuit models, by the use of mainly *reciprocal* [12] circuit elements, though non-reciprocal elements (gyrators) will come into play if is not symmetric (it is not required to be, and note that system (3.1) is well-posed regardless of the form of [82]). We have not explored the application of passive circuit methods to systems which are more generally strongly hyperbolic, for which energy estimates such as (3.7) can also be derived [82]. This would appear to be a worthy direction of future research.