Lawrence Berkeley National Laboratory

We present a variational multiscale formulation for the numerical solution of one-dimensional systems of conservation laws. The key idea of the proposed formulation, originally presented by Hughes [Comput. Methods Appl. Mech. Engrg., 127 (1995) 387-401], is a multiple-scale decomposition into resolved grid scales and unresolved subgrid scales. Incorporating the e.ect of the subgrid scales onto the coarse scale problem results in a .nite element method with enhanced stability properties, capable of accurately representing the sharp features of the solution. In the formulation developed herein, the multiscale split is invoked prior to any linearization of the equations. Special attention is given to the choice of the matrix of stabilizing coe.cients and the discontinuity capturing difusion. The methodology is applied to the one-dimensional simulation of three-phase .ow in porous media, and the shallow water equations. These numerical simulations clearly show the potential and applicability of the formulation for solving highly nonlinear, nearly hyperbolic systems on very coarse grids. Application of the numerical formulation to multidimensional problems is presented in a forthcoming paper.