The Gas Dynamics Equations

where is density, is volume velocity, is absolute pressure, and is total energy, internal plus kinetic. The three equations are not complete without a constitutive relation among the four dependent variables. The WDF people, in their treatment of

where is a constant which follows directly from thermodynamics (it is equal to the ratio of specific heats [203]).

Before proceeding any further, we mention the *scaling* of the dependent variables [16,49]; this is done, as in the linear problems discussed in Chapters 3 and 5, in order to optimize the stability condition on the resulting network. The variables are scaled as

The parameters and have dimensions of velocity and pressure, respectively, and nondimensionalize the system. will again become the space-step/time-step ratio in the numerical simulation routine, and plays a role similar to that of in the (1+1)D transmission line problem as discussed in §3.7 and follows directly from physical considerations.

Using the energy density definition (B.10), system (B.9) can be written in *non-conservative* [203] form (after some tedious algebraic manipulations) as

The problem here is that this system is not, in its present form, suitable for a circuit representation involving reciprocal elements^{}. Although it is not explicitly stated anywhere in the literature, the solution of Fettweis et al. has been to scale system (B.9) by left multiplication by the matrix diag
(though as previously mentioned, they work with the two-variable hydrodynamics system, or its analogues in higher dimensions). The scaled system takes the form:

While this scaling does not change smooth solutions to system (B.9), problems may occur if shocks are anticipated [181]. Such so-called

Finally, by employing the conservation of mass equation to simplify the other equations, the scaled system (B.12) can be written in *skew-selfadjoint form* [181] as:

where the symmetric matrices and are defined by

and is the state, , and . , in addition, will be positive definite if the density and pressure are positive everywhere. The significance of this skew-selfadjoint form can be seen by taking the inner product of (B.13) with , in which case we get

For the Cauchy problem (i.e., the problem is defined over the entire axis, so boundary conditions are effectively ignored), we may integrate over the domain to get

and thus is the global conserved quantity. It can be seen, by comparison between (B.14) and (B.15) with (3.3) and (3.5) (in the lossless case), that this skew-selfadjointness property is the natural extension of symmetric hyperbolicity to the nonlinear case. It is interesting that the generalized definitions of the inductor and capacitor, as per (B.1) and (B.3), are completely commensurate; we will see this in the next section. The theory of skew-selfadjoint forms has recently seen quite a bit of activity, in particular with regard to so-called

- MDKC and MDWDF for the Gas Dynamics Equations
- An Alternate MDKC and Scattering Network
- Entropy Variables