2023/02/21 by Siyu Cen, Cen, Siyu, Bangti Jin +5 · 2 citations
Environmental Science · Mathematics · Physics and Astronomy · #FOS: Mathematics #Groundwater flow and contamination studies #Model Reduction and Neural Networks #Numerical Analysis (math.NA) #Numerical methods in inverse problems
paper · pdf · doi:10.48550/arxiv.2302.10773
openalex publication_date 2023/02/21 · openalex created_date 2023/02/24 · openalex updated_date 2026/07/28
In this work we investigate the numerical identification of the diffusion coefficient in elliptic and parabolic problems using neural networks. The numerical scheme is based on the standard output least-squares formulation where the Galerkin finite element method (FEM) is employed to approximate the state and neural networks (NNs) act as a smoothness prior to approximate the unknown diffusion coefficient. A projection operation is applied to the NN approximation in order to preserve the physical box constraint on the unknown coefficient. The hybrid approach enjoys both rigorous mathematical foundation of the FEM and inductive bias / approximation properties of NNs. We derive a priori error estimates in the standard L2(Ω) norm for the numerical reconstruction, under a positivity condition which can be verified for a large class of problem data. The error bounds depend explicitly on the noise level, regularization parameter and discretization parameters (e.g., spatial mesh size, time step size, and depth, upper bound and number of nonzero parameters of NNs). We also provide extensive numerical experiments, indicating that the hybrid method is very robust for large noise when compared with the pure FEM approximation.