摘要
GDQ method is a kind of high order accurate numerical methods developed several years ago, which have been successfully used to simulate the solution of smooth engineering problems such as structure mechanics and incompressible fluid dynamics. In this paper, extending the traditional GDQ method, we develop a new kind of discontinuous GDQ methods to solve compressible flow problems of which solutions may be discontinuous. In order to capture the local features of fluid flows, firstly, the computational domain is divided into many small pieces of subdomains. Then, in each small subdomain, the GDQ method is implementedand some kinds of numerical flux limitation conditions will be required to keep the correct flow direction. At the boundary interface between subdomains, we also use some kind of flux conditions according to the flow direction. The numerical method obtained by the above steps has the advantages of high order accuracy and easy to treat boundary conditions. It can simulate perfectly nonlinear waves such as shock, rarefaction wave and contact discontinuity. Finally, the numerical experiments on one dimensional Burgers equation and Euler equations are given.The numerical results verify the validation of the method.
GDQ method is a kind of high order accurate numerical methods developed several years ago, which have been successfully used to simulate the solution of smooth engineering problems such as structure mechanics and incompressible fluid dynamics. In this paper, extending the traditional GDQ method, we develop a new kind of discontinuous GDQ methods to solve compressible flow problems of which solutions may be discontinuous. In order to capture the local features of fluid flows, firstly, the computational domain is divided into many small pieces of subdomains. Then, in each small subdomain, the GDQ method is implemented and some kinds of numerical flux limitation conditions will be required to keep the correct flow direction. At the boundary interface between subdomains, we also use some kind of flux conditions according to the flow direction. The numerical method obtained by the above steps has the advantages of high order accuracy and easy to treat boundary conditions. It can simulate perfectly nonlinear waves such as shock, rarefaction wave and contact discontinuity. Finally, the numerical experiments on one dimensional Burgers equation and Euler equations are given . The numerical results verify the validation of the method.
出处
《计算数学》
CSCD
北大核心
2004年第3期293-302,共10页
Mathematica Numerica Sinica
基金
航空科学基金项目(01A52003及02A52004)
"十五"国防预研项目资助.