A data-driven shock capturing approach for discontinuous Galekin methods

We propose a data-driven artificial viscosity model for shock capturing in discontinuous Galerkin methods. The proposed model trains a multi-layer feedforward network to map from the element-wise solution to a smoothness indicator, based on which the artificial viscosity is computed. The data set for the training of the network is obtained using canonical functions. The compactness of the data set, which is critical to the success of training the network, is ensured by normalization and the adjustment of the range of the smoothness indicator. The network is able to recover the expected smoothness much more reliably than the original averaged modal decay model. Several smooth and non-smooth test cases are considered to investigate the performance of this approach. Convergence tests show that the proposed model recovers the accuracy of the corresponding linear schemes for smooth regions. For non-smooth flows, the model is observed to suppress spurious oscillations well.

[1]  J. Templeton,et al.  Reynolds averaged turbulence modelling using deep neural networks with embedded invariance , 2016, Journal of Fluid Mechanics.

[2]  Jun Zhu,et al.  Hermite WENO Schemes and Their Application as Limiters for Runge-Kutta Discontinuous Galerkin Method, III: Unstructured Meshes , 2009, J. Sci. Comput..

[3]  Huayong Liu,et al.  Geometric conservation law and applications to high-order finite difference schemes with stationary grids , 2011, J. Comput. Phys..

[4]  Michael S. Gashler,et al.  Training Deep Fourier Neural Networks to Fit Time-Series Data , 2014, ICIC.

[5]  Sanjiva K. Lele,et al.  An artificial nonlinear diffusivity method for supersonic reacting flows with shocks , 2005, J. Comput. Phys..

[6]  Brendan D. Tracey,et al.  A Machine Learning Strategy to Assist Turbulence Model Development , 2015 .

[7]  Lilia Krivodonova,et al.  Limiters for high-order discontinuous Galerkin methods , 2007, J. Comput. Phys..

[8]  Rainald Löhner,et al.  A Hermite WENO-based limiter for discontinuous Galerkin method on unstructured grids , 2007, J. Comput. Phys..

[9]  J. Hesthaven Numerical Methods for Conservation Laws: From Analysis to Algorithms , 2017 .

[10]  Jean-Luc Guermond,et al.  Entropy viscosity method for nonlinear conservation laws , 2011, J. Comput. Phys..

[11]  Qiqi Wang,et al.  Quantification of structural uncertainties in the k − ω turbulence model , 2010 .

[12]  Ngoc Cuong Nguyen,et al.  Dilation‐based shock capturing for high‐order methods , 2016 .

[13]  Chi-Wang Shu,et al.  A Comparison of Troubled-Cell Indicators for Runge-Kutta Discontinuous Galerkin Methods Using Weighted Essentially Nonoscillatory Limiters , 2005, SIAM J. Sci. Comput..

[14]  Pierre Sagaut,et al.  A problem-independent limiter for high-order Runge—Kutta discontinuous Galerkin methods , 2001 .

[15]  Andrew W. Cook,et al.  Short Note: Hyperviscosity for shock-turbulence interactions , 2005 .

[16]  Gianluca Iaccarino,et al.  Modeling of structural uncertainties in Reynolds-averaged Navier-Stokes closures , 2013 .

[17]  J. Templeton Evaluation of machine learning algorithms for prediction of regions of high Reynolds averaged Navier Stokes uncertainty , 2015 .

[18]  P. Woodward,et al.  The numerical simulation of two-dimensional fluid flow with strong shocks , 1984 .

[19]  Chunlei Liang,et al.  Computation of flows with shocks using the Spectral Difference method with artificial viscosity, I: Basic formulation and application , 2014 .

[20]  J. Hesthaven,et al.  Nodal Discontinuous Galerkin Methods: Algorithms, Analysis, and Applications , 2007 .

[21]  Jianxian Qiu,et al.  Hermite WENO schemes and their application as limiters for Runge-Kutta discontinuous Galerkin method II: Two dimensional case , 2005 .

[22]  Jean-Luc Guermond,et al.  Implementation of the entropy viscosity method with the discontinuous Galerkin method , 2013 .

[23]  Heng Xiao,et al.  Quantifying and reducing model-form uncertainties in Reynolds-averaged Navier-Stokes simulations: A data-driven, physics-informed Bayesian approach , 2015, J. Comput. Phys..

[24]  Sai Hung Cheung,et al.  Bayesian uncertainty analysis with applications to turbulence modeling , 2011, Reliab. Eng. Syst. Saf..

[25]  Eitan Tadmor,et al.  Solution of two‐dimensional Riemann problems for gas dynamics without Riemann problem solvers , 2002 .

[26]  H. T. Huynh,et al.  A Flux Reconstruction Approach to High-Order Schemes Including Discontinuous Galerkin Methods , 2007 .

[27]  Michael Dumbser,et al.  Runge-Kutta Discontinuous Galerkin Method Using WENO Limiters , 2005, SIAM J. Sci. Comput..

[28]  Parviz Moin,et al.  Suitability of artificial bulk viscosity for large-eddy simulation of turbulent flows with shocks , 2009, J. Comput. Phys..

[29]  Soshi Kawai,et al.  Localized artificial diffusivity scheme for discontinuity capturing on curvilinear meshes , 2008, J. Comput. Phys..

[30]  Zhi Jian Wang,et al.  A unifying lifting collocation penalty formulation including the discontinuous Galerkin, spectral volume/difference methods for conservation laws on mixed grids , 2009, J. Comput. Phys..

[31]  Chi-Wang Shu,et al.  The Runge-Kutta Discontinuous Galerkin Method for Conservation Laws V , 1998 .

[32]  Wai-Sun Don,et al.  Explicit discontinuous spectral element method with entropy generation based artificial viscosity for shocked viscous flows , 2017, J. Comput. Phys..

[33]  G. Lewicki,et al.  Approximation by Superpositions of a Sigmoidal Function , 2003 .

[34]  Chi-Wang Shu,et al.  A simple weighted essentially nonoscillatory limiter for Runge-Kutta discontinuous Galerkin methods , 2013, J. Comput. Phys..

[35]  Karthik Duraisamy,et al.  A paradigm for data-driven predictive modeling using field inversion and machine learning , 2016, J. Comput. Phys..

[36]  J. Peraire,et al.  Sub-Cell Shock Capturing for Discontinuous Galerkin Methods , 2006 .

[37]  Chi-Wang Shu,et al.  Runge-Kutta Discontinuous Galerkin Method Using WENO Limiters , 2005, SIAM J. Sci. Comput..

[38]  Jun Zhu,et al.  Runge-Kutta discontinuous Galerkin method using a new type of WENO limiters on unstructured meshes , 2013, J. Comput. Phys..

[39]  O. Friedrich,et al.  Weighted Essentially Non-Oscillatory Schemes for the Interpolation of Mean Values on Unstructured Grids , 1998 .

[40]  W. Cabot,et al.  A high-wavenumber viscosity for high-resolution numerical methods , 2004 .

[41]  Norbert Kroll ADIGMA - A European Initiative on the Development of Adaptive Higher-Order Variational Methods for Aerospace Applications : Results of a collaborative research project funded by the European Union, 2006-2009 , 2010 .

[42]  P. A. Durbin,et al.  Transition modeling using data driven approaches , 2014 .

[43]  J. S. Hesthaven,et al.  Viscous Shock Capturing in a Time-Explicit Discontinuous Galerkin Method , 2011, 1102.3190.