GLIMM'S METHOD FOR GAS DYNAMICS
We investigate Glimm's method, a method for constructing approximate solutions to systems of hyperbolic conservation laws in one space variable by sampling explicit wave solutions. It is extended to several space variables by operator splitting. We consider two functional problems. 1) We propose a highly accurate form of the sampling procedure, in one space variable, based on the van der Corput sampling sequence. We test the improved sampling procedure numerically in the case of inviscid compressible flow in one space dimension and find that it gives high resolution results both in the smooth parts of the solution, as well as the discontinuities. 2) We investigate the operator splitting procedure by means of which the multidimensional method is constructed. An 0(1) error stemming from the use of this procedure near shocks oblique to the spatial grid is analyzed numerically in the case of the equations for inviscid compressible flow in two space dimensions. We present a hybrid method which eliminates this error, consisting of Glimm's method, used in continuous parts of the flow, and the nonlinear Godunov's method, used in regions where large pressure jumps are generated. The resulting method is seen to be a substantial improvement over either of the component methods for multidimensional calculations.
- Research Organization:
- Lawrence Berkeley National Lab. (LBNL), Berkeley, CA (United States)
- Sponsoring Organization:
- Computing Sciences Directorate; Physics Division
- DOE Contract Number:
- DE-AC02-05CH11231
- OSTI ID:
- 1077126
- Report Number(s):
- LBL-11532
- Country of Publication:
- United States
- Language:
- English
Similar Records
Analysis of the effect of operator splitting and of the sampling procedure on the accuracy of Glimm's method
Computer implementation of Glimm's method. [Subroutine GLIMM, for solution of nonlinear hyperbolic equations, in FORTRAN for CDC 7600 computer]