IMPROVED ENTROPY-ULTRA-BEE SCHEME FOR THE EULER SYSTEM OF GAS DYNAMICS  

作　　者：Rongsan Chen Dekang Mao 

机构地区：[1]School of Mathematics and Physics, China University of Geosciences, Wuhan, China [2]Department of Mathematics, Shanghai University, Shanghai, China

出　　处：《Journal of Computational Mathematics》2017年第2期121-151,共31页计算数学（英文）

基　　金：The research is supported by the National Natural Science Foundation of China No.10971132, No.11176015, No.11201436, No.11201435, the Leading Academic Discipline Project of Shanghai Municipal Education Commission, No.J50101, the National Basic Research Program of China（973 program） No.2011CB71064, No.2011CBT1065, and the Fundamental Research Funds for the Central Universities, China University of Geosciences（Wuhan）.

摘　　要：The Entropy-Ultra-Bee scheme was developed for the linear advection equation and extended to the Euler system of gas dynamics in [13]. It was expected that the technology be applied only to the second characteristic field of the system and the computation in the other two nonlinear fields be implemented by the Godunov scheme. However, the numerical experiments in [13] showed that the scheme, though having improved the wave resolution in the second field, produced numerical oscillations in the other two nonlinear fields. Sophisticated entropy increaser was designed to suppress the spurious oscillations by increasing the entropy when there are waves in the two nonlinear fields presented. However, the scheme is then not efficient neither robust with problem-related parameters. The purpose of this paper is to fix this problem. To this end, we first study a 3 × 3 linear system and apply the technology precisely to its second characteristic field while maintaining the computation in the other two fields be implemented by the Godunov scheme. We then follow the discussion for the linear system to apply the Entropy-Ultra-Bee technology to the second characteristic field of the Euler system in a linearlized field-by- field fashion to develop a modified Entropy-Ultra-Bee scheme for the system. Meanwhile a remark is given to explain the problem of the previous Entropy-Ultra-Bee scheme in [13]. A reference solution is constructed for computing the numerical entropy, which maintains the feature of the density and flats the velocity and pressure to constants. The numerical entropy is then computed as the entropy cell-average of the reference solution. Several limitations are adopted in the construction of the reference solution to further stabilize the scheme. Designed in such a way, the modified Entropy-Ultra-Bee scheme has a unified form with no problem-related parameters. Numerical experiments show that all the spurious oscillations in smooth regions are gone and the results are better than that of the previous Entropy-UltrThe Entropy-Ultra-Bee scheme was developed for the linear advection equation and extended to the Euler system of gas dynamics in [13]. It was expected that the technology be applied only to the second characteristic field of the system and the computation in the other two nonlinear fields be implemented by the Godunov scheme. However, the numerical experiments in [13] showed that the scheme, though having improved the wave resolution in the second field, produced numerical oscillations in the other two nonlinear fields. Sophisticated entropy increaser was designed to suppress the spurious oscillations by increasing the entropy when there are waves in the two nonlinear fields presented. However, the scheme is then not efficient neither robust with problem-related parameters. The purpose of this paper is to fix this problem. To this end, we first study a 3 × 3 linear system and apply the technology precisely to its second characteristic field while maintaining the computation in the other two fields be implemented by the Godunov scheme. We then follow the discussion for the linear system to apply the Entropy-Ultra-Bee technology to the second characteristic field of the Euler system in a linearlized field-by- field fashion to develop a modified Entropy-Ultra-Bee scheme for the system. Meanwhile a remark is given to explain the problem of the previous Entropy-Ultra-Bee scheme in [13]. A reference solution is constructed for computing the numerical entropy, which maintains the feature of the density and flats the velocity and pressure to constants. The numerical entropy is then computed as the entropy cell-average of the reference solution. Several limitations are adopted in the construction of the reference solution to further stabilize the scheme. Designed in such a way, the modified Entropy-Ultra-Bee scheme has a unified form with no problem-related parameters. Numerical experiments show that all the spurious oscillations in smooth regions are gone and the results are better than that of the previous Entropy-Ultr

关 键 词：Entropy-Ultra-Bee scheme Step-reconstruction Characteristic field Referencesolution. 

分 类 号：O[理学]