0 and with the values m u(x, t) lying in an open subset, G, of R , the state space. The state space G arises because physical quantities such as the density or total energy should always be positive, thus the values of u are often con strained to an open set G. 2) The flux functions appearing in these balance laws are idealized through prescribed nonlinear functions, F.(u), mapping G into J j = 1, ... , N while source terms are defined by S(u, x, t) with S a given smooth function of these arguments with values in Rm. In parti- lar, the detailed microscopic effects of diffusion and dissipation are ignored. 3) A generalized version of the principle of virtual work is applied (see Antman [1]). The formal result of applying the three steps (1)-(3) is that the m physical quantities u define a weak solution of an m x m system of conservation laws, o I + N(Wt'u + r W ·F.(u) + W·S(u, x, t))dxdt (1.1) R xR j=l Xj J for all W E C~(RN x R+), W(x, t) E Rm." />