Numerical Flux Functions for RANS-k-omega Computations with a Line-Preconditioned p-Multigrid DG Solver

Author(s):  
Marcel Wallraff ◽  
Tobias Leicht ◽  
Markus Lange-Hegermann
Keyword(s):  
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.


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.


2018 ◽  
Vol 358 ◽  
pp. 1-20 ◽  
Author(s):  
Gianmarco Mengaldo ◽  
Daniele De Grazia ◽  
Rodrigo C. Moura ◽  
Spencer J. Sherwin

2019 ◽  
Vol 219 (3) ◽  
pp. 1900-1914 ◽  
Author(s):  
T Möller ◽  
W Friederich

SUMMARY An existing nodal discontinuous Galerkin (NDG) method for the simulation of seismic waves in heterogeneous media is extended to media containing fractures with various rheological behaviour. Fractures are treated as 2-D surfaces where Schoenberg’s linear slip or displacement discontinuity condition is applied as an additional boundary condition to the elastic wave equation which is in turn implemented as an additional numerical flux within the NDG formulation. Explicit expressions for the new numerical flux are derived by considering the Riemann problem for the elastic wave equation at fractures with varying rheologies. In all cases, we obtain further first order differential equations that fully describe the temporal evolution of the particle velocity jump at the fracture. Our flux formulation allows to separate the effect of a fracture from flux contributions due to simple welded interfaces enabling us to easily declare element faces as parts of a fracture. We make use of this fact by first generating the numerical mesh and then building up fractures by selecting appropriate element faces instead of adjusting the mesh to pre-defined fracture surfaces. The implementation of the new numerical fluxes into NDG is verified in 1-D by comparison to an analytical solution and in 2-D by comparing the results of a simulation valid in 1-D and 2-D. Further numerical examples address the effect of fracture systems on seismic wave propagation in 1-D and 2-D featuring effective anisotropy and coda generation. Finally, a study of the reflective and transmissive behaviour of fractures indicates that reflection and transmission coefficients are controlled by the ratio of signal frequency and relaxation frequency of the fracture.


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.


Sign in / Sign up

Export Citation Format

Share Document