Author(s):  
Michael Dumbser ◽  
Claus-Dieter Munz

On Source Terms and Boundary Conditions Using Arbitrary High Order Discontinuous Galerkin SchemesThis article is devoted to the discretization of source terms and boundary conditions using discontinuous Galerkin schemes with an arbitrary high order of accuracy in space and time for the solution of hyperbolic conservation laws on unstructured triangular meshes. The building block of the method is a particular numerical flux function at the element interfaces based on the solution of Generalized Riemann Problems (GRPs) with piecewise polynomial initial data. The solution of the generalized Riemann problem, originally introduced by Toro and Titarev in a finite volume context, provides simultaneously a numerical flux function as well as a time integration method. The resulting scheme is extremely local since it integrates the PDE from one time step to the successive one in a single step using only information from the direct side neighbors. Since source terms are directly incorporated into the numerical flux via the solution of the GRP, our very high order accurate method is also able to maintain very well smooth steady-state solutions of PDEs with source terms, similar to the so-called well-balanced schemes which are usually specially designed for this purpose. Boundary conditions are imposed solving inverse generalized Riemann problems. Furthermore, we show numerical evidence proving that by using very high order schemes together with high order polynomial representations of curved boundaries, high quality solutions can be obtained on very coarse meshes.


Author(s):  
U. S. Vevek ◽  
B. Zang ◽  
T. H. New

AbstractA hybrid numerical flux scheme is proposed by adapting the carbuncle-free modified Harten-Lax-van Leer contact (HLLCM) scheme to smoothly revert to the Harten-Lax-van Leer contact (HLLC) scheme in regions of shear. This hybrid scheme, referred to as the HLLCT scheme, employs a novel, velocity-based shear sensor. In contrast to the non-local pressure-based shock sensors often used in carbuncle cures, the proposed shear sensor can be computed in a localized manner meaning that the HLLCT scheme can be easily introduced into existing codes without having to implement additional data structures. Through numerical experiments, it is shown that the HLLCT scheme is able to resolve shear layers accurately without succumbing to the shock instability.


2021 ◽  
Vol 36 (6) ◽  
pp. 313-336
Author(s):  
Ronald H. W. Hoppe ◽  
Youri Iliash

Abstract We are concerned with an Interior Penalty Discontinuous Galerkin (IPDG) approximation of the p-Laplace equation and an equilibrated a posteriori error estimator. The IPDG method can be derived from a discretization of the associated minimization problem involving appropriately defined reconstruction operators. The equilibrated a posteriori error estimator provides an upper bound for the discretization error in the broken W 1,p norm and relies on the construction of an equilibrated flux in terms of a numerical flux function associated with the mixed formulation of the IPDG approximation. The relationship with a residual-type a posteriori error estimator is established as well. Numerical results illustrate the performance of both estimators.


Author(s):  
Chanyoung Park ◽  
Samaun Nili ◽  
Justin Mathew ◽  
Frederick Ouellet ◽  
Rahul Koneru ◽  
...  

Abstract Uncertainty quantification (UQ) is an important step in the verification and validation of scientific computing. Validation is often inconclusive when uncertainty is larger than an acceptable range for both simulation and experiment. Therefore, uncertainty reduction (UR) is important to achieve meaningful validation. A unique approach in this paper is to separate model error from uncertainty such that UR can reveal the model error. This paper aims to share lessons learned from UQ and UR of a horizontal shock tube simulation, whose goal is to validate the particle drag force model for the compressible multiphase flow. Firstly, simulation UQ revealed the inconsistency in simulation predictions due to the numerical flux scheme, which was clearly shown using the parametric design of experiments. By improving the numerical flux scheme, the uncertainty due to inconsistency was removed, while increasing the overall prediction error. Secondly, the mismatch between the geometry of the experiments and the simplified 1D simulation model was identified as a lack of knowledge. After modifying simulation conditions and experiments, it turned out that the error due to the mismatch was small, which was unexpected based on expert opinions. Lastly, the uncertainty in the initial volume fraction of particles was reduced based on rigorous UQ. All these UR measures worked together to reveal the hidden modeling error in the simulation predictions, which can lead to a model improvement in the future. We summarized the lessons learned from this exercise in terms of empty success, useful failure, and deceptive success.


2015 ◽  
Vol 2015 ◽  
pp. 1-11 ◽  
Author(s):  
Thida Pongsanguansin ◽  
Montri Maleewong ◽  
Khamron Mekchay

A well-balanced scheme with total variation diminishing Runge-Kutta discontinuous Galerkin (TVD-RK DG) method for solving shallow water equations is presented. Generally, the flux function at cell interface in the TVD-RK DG scheme is approximated by using the Harten-Lax-van Leer (HLL) method. Here, we apply the weighted average flux (WAF) which is higher order approximation instead of using the HLL in the TVD-RK DG method. The consistency property is shown. The modified well-balanced technique for flux gradient and source terms under the WAF approximations is developed. The accuracy of numerical solutions is demonstrated by simulating dam-break flows with the flat bottom. The steady solutions with shock can be captured correctly without spurious oscillations near the shock front. This presents the other flux approximations in the TVD-RK DG method for shallow water simulations.


Sign in / Sign up

Export Citation Format

Share Document