Abstract

Abstract

In this work the NURBS-Enhanced Finite Element Method (NEFEM) is combined with a Discontinuous Galerkin (DG) formulation for the numerical solution of the Euler equations of gas dynamics. With the NEFEM approach numerical fluxes along curved boundaries are computed much more accurately due to the exact geometric representation of the computational domain. The proper implementation of the wall boundary condition provides accurate results even with a linear interpolation of the solution. A detailed comparison of the NEFEM in front of isoparametric finite elements (FE) is presented, demonstrating the superiority of the NEFEM approach for both linear and higher order computations.

Abstract

In this work, the NURBS-enhanced finite element method (NEFEM) is combined with a discontinuous Galerkin (DG) formulation for the numerical solution of Euler equations of gas dynamics. With NEFEM, numerical fluxes along curved boundaries are computed much more accurately due to the exact geometric representation of the computational domain. The proper implementation of the wall boundary condition and the exact geometry provide accurate results even with a linear approximation of the solution. A detailed comparison of NEFEM in front of isoparametric finite elements is presented, demonstrating the superiority of NEFEM approach for both linear and higher-order computations

Abstract

This article presents a novel shock-capturing technique for the discontinuous Galerkin (DG) method. The technique is designed for compressible flow problems, which are usually characterized by the presence of strong shocks and discontinuities. The inherent structure of standard DG methods seems to suggest that they are especially adapted to capture shocks because of the numerical fluxes based on suitable approximate Riemann solvers, which, in practice, introduces some stabilization. However, the usual numerical fluxes are not sufficient to stabilize the solution in the presence of shocks for large high-order elements. Here, a new basis of shape functions is introduced. It has the ability to change locally between a continuous or discontinuous interpolation depending on the smoothness of the approximated function. In the presence of shocks, the new discontinuities inside an element introduce the required stabilization because of numerical fluxes. Large high-order elements can therefore be used and shocks captured within a single element, avoiding adaptive mesh refinement and preserving the locality and compactness of the DG scheme. Several numerical examples for transonic and supersonic flows are studied to demonstrate the applicability of the proposed approach.

Abstract

Abstract

Abstract

This study is about the construction of a numerical scheme of the predictor-corrector type in conservative form for solving general systems of conservation laws in multiple space dimensions on unstructured meshes. The work is a generalization of the one-dimensional finite volume characteristics (FVC) scheme and is related to the work of Fayssal Benkhaldoun and Mohammed Sead. The construction of the intermediate state is based on the method of characteristics, while the corrective stage recovers the conservation equations. The scheme is accurate to first order, monotonic and entropic; it avoids Riemann solvers at each interface; it also allows for improved accuracy order in time and space on unstructured three-dimensional meshes in the framework of the finite volume method. The scheme's performance is evaluated through a series of test benchmarks for the three-dimensional version of the Euler equations..

Abstract

This paper aims at the development of a new stabilization formulation based on the Finite Calculus (FIC) scheme for solving the Euler equations using the Galerkin finite element method (FEM) on unstructured triangular grids. The FIC method is based on expressing the balance of fluxes in a space-time domain of finite size. It is used to prevent the creation of instabilities typically present in numerical solutions due to the high convective terms and sharp gradients. Two stabilization terms, respectively called streamline term and transverse term, are added via the FIC formulation to the original conservative equations in the spacetime domain. An explicit fourth order Runge-Kutta scheme is implemented to advance the solution in time. The presented numerical test examples for inviscid flows prove the ability of the proposed stabilization technique for providing appropriate solutions especially near shock waves. Although the derived methodology delivers precise results with a nearly coarse mesh, a mesh refinement technique is coupled to the solution process for obtaining a suitable mesh particularly in the high gradient zones.

Abstract