Matteo Caldana,Paola F. Antonietti,Luca Dede'
Abstract:Finite element-based high-order solvers of conservation laws offer large accuracy but face challenges near discontinuities due to the Gibbs phenomenon. Artificial viscosity is a popular and effective solution to this problem based on physical insight. In this work, we present a physics-informed machine learning algorithm to automate the discovery of artificial viscosity models in a non-supervised paradigm. The algorithm is inspired by reinforcement learning and trains a neural network acting cell-by-cell (the viscosity model) by minimizing a loss defined as the difference with respect to a reference solution thanks to automatic differentiation. This enables a dataset-free training procedure. We prove that the algorithm is effective by integrating it into a state-of-the-art Runge-Kutta discontinuous Galerkin solver. We showcase several numerical tests on scalar and vectorial problems, such as Burgers' and Euler's equations in one and two dimensions. Results demonstrate that the proposed approach trains a model that is able to outperform classical viscosity models. Moreover, we show that the learnt artificial viscosity model is able to generalize across different problems and parameters.
What problem does this paper attempt to address?
### Problems the paper attempts to solve
This paper aims to solve the Gibbs phenomenon problem encountered when the discontinuous Galerkin (DG) method is used to solve conservation law equations. Specifically, although high - order finite element methods have high precision when dealing with conservation law equations, numerical oscillations will occur near the discontinuous regions of the solution due to the Gibbs phenomenon, and these oscillations will destroy the physical meaning of the solution. To overcome this problem, artificial viscosity (AV) is a commonly used solution. It restores the smoothness of the solution by adding an appropriate amount of dissipation in the discontinuous regions, thereby maintaining high - order accuracy.
However, existing artificial viscosity models have some limitations:
1. **Parameter - dependence**: Most artificial viscosity models require manual parameter adjustment, which is not only time - consuming but also requires in - depth professional knowledge.
2. **Computational cost**: Evaluating some artificial viscosity models (such as the entropy viscosity model) can be very expensive.
3. **Generalization ability**: Existing models are usually optimized for specific problems and are difficult to generalize to other problems or parameter settings.
To solve these problems, the paper proposes a hybrid learning framework based on physics - informed machine learning (PIML) for automatically discovering new and efficient artificial viscosity models. This framework combines the advantages of supervised learning and unsupervised learning, and uses the idea of reinforcement learning to train the neural network by minimizing the error between the reference solution and the current solution. This method not only avoids the need to construct large - scale data sets but also shows good generalization ability under different problems and parameter settings.
### Main contributions
1. **Hybrid learning framework**: A hybrid learning framework between supervised learning and unsupervised learning is proposed, using physical information to enhance the learning process.
2. **Automatically discover artificial viscosity models**: New and efficient artificial viscosity models are automatically discovered through neural networks without manual parameter tuning.
3. **Generalization ability**: The new model can perform well under different problems and parameter settings and has good generalization ability.
4. **Computational efficiency**: By optimizing input features and network structures, computational efficiency is improved and training costs are reduced.
### Mathematical models and methods
1. **Conservation law equations**: Consider the conservation law equations defined on a bounded open domain \(\Omega\subset\mathbb{R}^d\), where \(d\) can be 1 or 2, and time \(t\in(0, T]\), \(T > 0\).
2. **Discontinuous Galerkin method**: The discontinuous Galerkin method is used for spatial discretization, and the semi - discrete and algebraic forms of the equations are defined.
3. **Artificial viscosity models**: Classic artificial viscosity models are reviewed, and new models based on neural networks are introduced.
4. **Training algorithms**: A training algorithm based on physics - informed machine learning is proposed, which minimizes the error between the reference solution and the current solution through automatic differentiation techniques.
### Experimental results
The paper presents multiple numerical tests, including one - dimensional and two - dimensional scalar and vector problems (such as the Burgers equation and the Euler equation). The experimental results show that the proposed neural network artificial viscosity model can effectively reduce numerical oscillations, improve the accuracy of the solution, and perform well under different problems and parameter settings.
### Conclusion
The paper proposes a hybrid learning framework based on physics - informed machine learning for automatically discovering efficient artificial viscosity models. This method not only avoids the complexity of manual parameter tuning but also improves computational efficiency and generalization ability, providing a new and effective way to solve the Gibbs phenomenon in conservation law equations.