Computational & Technology Resources
an online resource for computational,
engineering & technology publications 

CivilComp Proceedings
ISSN 17593433 CCP: 84
PROCEEDINGS OF THE FIFTH INTERNATIONAL CONFERENCE ON ENGINEERING COMPUTATIONAL TECHNOLOGY Edited by: B.H.V. Topping, G. Montero and R. Montenegro
Paper 4
An Adaptive Mesh Method for Hyperbolic Conservation Laws J. Felcman^{1} and P. Kubera^{2}
^{1}Faculty of Mathematics and Physics, Charles University, Prague, Czech Republic
J. Felcman, P. Kubera, "An Adaptive Mesh Method for Hyperbolic Conservation Laws", in B.H.V. Topping, G. Montero, R. Montenegro, (Editors), "Proceedings of the Fifth International Conference on Engineering Computational Technology", CivilComp Press, Stirlingshire, UK, Paper 4, 2006. doi:10.4203/ccp.84.4
Keywords: mesh adaptation, geometric conservation law, finite volume method, moving discontinuities, Euler equations.
Summary
The subjectmatter of this paper is the numerical simulation of the multidimensional inviscid
compressible transonic gas flow in channels, past airfoils and cascades of profiles.
The adaptive strategy is applied to the numerical solution of problem governed by hyperbolic
partial differential equations. An adaptive mesh is constructed in the framework of the cellcentred
finite volume scheme. An anisotropic mesh adaptation strategy is followed by a recovery of the
approximate solution on the new mesh. The main feature of the proposed method is to keep the
mass conservation of the numerical solution at each adaptation step. This allows the solution of the
nonstationary problems. The geometric conservation law is employed.
The numerical solution for the case of a nonstationary discontinuity wave is presented.
We consider the flow of an inviscid perfect gas in a bounded domain and time interval with . Here or 3 for two or threedimensional flow, respectively, and we suppose that is polygonal in two dimensions or polyhedral in three dimensions, respectively. Further we suppose, that the flow is adiabatic and we neglect the outer volume force. Our goal is to solve numerically the initialboundary value problem described by the Euler equations. As a discrete analogy of the continuous problem we use the ADER higherorder finite volume scheme based on a polynomial reconstruction of a piecewise constant finite volume solution. For the study of the numerical order of accuracy of the proposed finite volume scheme see [1,2] and references therein. In [3] a time marching finite volume method for nonstationary problems was proposed and tested in for onedimensional case. Here we present its further development and twodimensional results. The algorithm is generally formulated in three dimensions. It consists of three basic sections at each time step: the time evolution of the numerical solution, the mesh adaptation and the recomputing of the numerical solution from the mesh before the adaptation to the mesh after the adaptation. In one time step the finite volume scheme is evaluated twice. Firstly for the prediction how to adapt the mesh, further for the update of the numerical solution itself. In the prediction part, we forecast the evolution of the numerical solution and adapt the mesh. The anisotropic mesh adaptation (AMA) is applied. For its description see e.g. [4]. In [4] the necessary condition for the properties of the Nsimplicial mesh, on which the discretization error is below the prescribed tolerance, is formulated. It is shown, how to control this necessary condition by the interpolation error and the anisotropic mesh adaptation technique is applied. For two and threedimensional numerical examples see e.g. [5]. In the AMA technique, the equilateral mesh is constructed in the least squares sense. The length of an edge of an Nsimplicial mesh is measured in the numerical solution dependent Riemann norm. After the adaptive mesh is constructed it is necessary to recompute the solution on the old mesh to its recovery on the newly adapted mesh. According to [6] the geometric mass conservation law has to be satisfied in this computational step. The perturbation method from [6] is applied. The twodimensional numerical example employing the described higherorder method with proposed mesh adaptation is presented. The numerical aspects of the method are demonstrated with the example of the twodimensional channel flow described by the Euler equations equipped with the discontinuous initial condition. The mesh adaptation corresponding to the various time instants is depicted. The proposed nonstationary mesh adaptation can help to improve the numerical results of a cell centered finite volume method. Numerical tests in progress show that it prevents the numerical solution from spurious oscillations. The proposed adaptation is value for the sufficiently precise solution despite increased computational time. References
purchase the fulltext of this paper (price £20)
go to the previous paper 
