Convex ENO high order multi-dimensional schemes without field by field decomposition or staggered grids

被引:141
作者
Liu, XD [1 ]
Osher, S
机构
[1] Univ Calif Santa Barbara, Dept Math, Santa Barbara, CA 93106 USA
[2] Univ Calif Los Angeles, Dept Math, Los Angeles, CA 90095 USA
基金
美国国家科学基金会;
关键词
D O I
10.1006/jcph.1998.5937
中图分类号
TP39 [计算机的应用];
学科分类号
081203 ; 0835 ;
摘要
Second order accurate (first order at extrema) cell averaged based approximations extending the Lax-Friedrichs central scheme, using component-wise rather than field-by-field limiting, have been found to give surprisingly good results for a wide class of problems involving shocks (see H. Nessyahu and E. Tadmor, J. Comput. Phys. 87, 408, 1990). The advantages of component-wise limiting compared to its counterpart, field-by-field limiting, are apparent: (1) No complete set of eigenvectors is needed and hence weakly hyperbolic systems can be solved. (2) Componentwise limiting is faster than field-by-field limiting. (3) The programming is much simpler, especially for complicated coupled systems of many equations. However, these methods are based on cell-averages in a staggered grid and are thus a bit complicated to extend to multiple dimensions. Moreover the staggering causes slight difficulties at the boundaries. In this work we modify and extend this component-wise central differencing based procedure in two directions: (1) Point values, rather than cell averages are used, thus removing the need for staggered grids, and also making the extension to multi-dimensions quite simple. We use TVD Runge-Kutta time discretizations to update the solution. (2) A new type of decision process, which follows the general ENO philosophy is introduced and used. This procedure enables us to extend our method to a third order component-wise central ENO scheme, which apparently works well and is quite simple to implement in multi-dimensions. Additionally, our numerical viscosity is governed by the local magnitude of the maximum eigenvalue of the Jacobian, thus reducing the smearing in the numerical results. We found a speed up of a factor of 2 in each space dimension, on a SGI O-2 workstation, over methods based on field-by-field decomposition limiting. The new decision process leads to new, "convex" ENO schemes which, we believe, are of interest in a more general setting. Our numerical results show the value of these new methods. (C) 1998 Academic Press.
引用
收藏
页码:304 / 330
页数:27
相关论文
共 22 条