Machine Learning Changes the Rules for Flux Limiters

Learning to integrate non-linear equations from highly resolved direct numerical simulations (DNSs) has seen recent interest due to the potential to reduce the computational load for numerical fluid simulations. Here, we focus on a specific problem in the integration of fluids: the determination of a flux-limiter for shock capturing methods. Focusing on flux limiters learned from data has the advantage of providing specific plug-and-play components for existing numerical methods. Since their introduction, a large array of flux limiters has been designed. Using the example of the coarse-grained Burgers’ equation, we show that fluxlimiters may be rank-ordered in terms of how well they integrate various coarse-grainings. We then develop theory to find an optimal flux-limiter and present a set of flux-limiters that outperform others tested for integrating Burgers’ equation on lattices with 2×, 3×, 4×, and 8× coarse-grainings. Our machine learned flux limiters have distinctive features that may provide new rules-of-thumb for the development of improved limiters. Additionally, we optimize over hyper-parameters, including coarse-graining, number of discretized bins, and diffusion parameter to find flux limiters that are best on average and more broadly useful for general applications. Our study provides a basis for understanding how flux limiters may be improved upon using machine learning methods.

[1]  B. V. Leer,et al.  Towards the ultimate conservative difference scheme V. A second-order sequel to Godunov's method , 1979 .

[2]  S. W. Schoombie,et al.  Exact analysis of nonlinear instability in a discrete Burgers' equation , 1991 .

[3]  P. Gaskell,et al.  Curvature‐compensated convective transport: SMART, A new boundedness‐ preserving transport algorithm , 1988 .

[4]  S. Osher,et al.  Weighted essentially non-oscillatory schemes , 1994 .

[5]  B. V. Leer,et al.  Towards the ultimate conservative difference scheme III. Upstream-centered finite-difference schemes for ideal compressible flow , 1977 .

[6]  Daniel Livescu,et al.  Application of a second-moment closure model to mixing processes involving multicomponent miscible fluids , 2011 .

[7]  Chris Benmore,et al.  Machine learning coarse grained models for water , 2019, Nature Communications.

[8]  J. Burgers A mathematical model illustrating the theory of turbulence , 1948 .

[9]  Barry Koren,et al.  A robust upwind discretization method for advection, diffusion and source terms , 1993 .

[10]  P. Woodward,et al.  The Piecewise Parabolic Method (PPM) for Gas Dynamical Simulations , 1984 .

[11]  P. Sweby High Resolution Schemes Using Flux Limiters for Hyperbolic Conservation Laws , 1984 .

[12]  S. Osher,et al.  Efficient implementation of essentially non-oscillatory shock-capturing schemes,II , 1989 .

[13]  S. Osher,et al.  Uniformly High-Order Accurate Nonoscillatory Schemes. I , 1987 .

[14]  A. Mohan,et al.  Compressed Convolutional LSTM: An Efficient Deep Learning framework to Model High Fidelity 3D Turbulence , 2019, 1903.00033.

[15]  G. D. van Albada,et al.  A comparative study of computational methods in cosmic gas dynamics , 1982 .

[16]  G. Pinton,et al.  Piecewise parabolic method for simulating one-dimensional shear shock wave propagation in tissue-mimicking phantoms , 2017 .

[17]  Michael M. McKerns,et al.  Building a Framework for Predictive Science , 2012, SciPy.

[18]  A. Harten High Resolution Schemes for Hyperbolic Conservation Laws , 2017 .

[19]  J. Cole On a quasi-linear parabolic equation occurring in aerodynamics , 1951 .

[20]  Liang Cheng,et al.  A review on TVD schemes and a refined flux-limiter for steady-state calculations , 2015, J. Comput. Phys..

[21]  Chang Ho Kim,et al.  MAPPED WENO SCHEMES BASED ON A NEW SMOOTHNESS INDICATOR FOR HAMILTON-JACOBI EQUATIONS , 2012 .

[22]  B. V. Leer,et al.  Towards the ultimate conservative difference scheme. II. Monotonicity and conservation combined in a second-order scheme , 1974 .

[23]  Ken Perlin,et al.  Accelerating Eulerian Fluid Simulation With Convolutional Networks , 2016, ICML.

[24]  Rainer Storn,et al.  Differential Evolution – A Simple and Efficient Heuristic for global Optimization over Continuous Spaces , 1997, J. Glob. Optim..

[25]  Fue-Sang Lien,et al.  Upstream monotonic interpolation for scalar transport with application to complex turbulent flows , 1994 .

[26]  Daniel Livescu,et al.  Leveraging Bayesian analysis to improve accuracy of approximate models , 2019, J. Comput. Phys..

[27]  P. Roe CHARACTERISTIC-BASED SCHEMES FOR THE EULER EQUATIONS , 1986 .