A Deep Double Ritz Method (D2RM) for solving Partial Differential Equations using Neural Networks

被引:14
作者
Uriarte, Carlos [1 ,2 ]
Pardo, David [1 ,2 ,3 ]
Muga, Ignacio [1 ,4 ]
Munoz-Matute, Judit [1 ,5 ]
机构
[1] Basque Ctr Appl Math BCAM, Alameda Mazarredo 14, E-48009 Bilbao, Spain
[2] Univ Basque Country UPV EHU, E-48940 Leioa, Spain
[3] Basque Fdn Sci Ikerbasque, Plaza Euskadi 5, E-48009 Bilbao, Spain
[4] Pontifical Catholic Univ Valparaiso PUCV, Ave Brasil 2950, Valparaiso, Chile
[5] Oden Inst Computat Engn & Sci OICES, 201 E 24th St, Austin, TX 78712 USA
关键词
Partial Differential Equations; Variational formulation; Residual minimization; Ritz Method; Optimal test functions; Neural networks; FRAMEWORK; ALGORITHM;
D O I
10.1016/j.cma.2023.115892
中图分类号
T [工业技术];
学科分类号
08 ;
摘要
Residual minimization is a widely used technique for solving Partial Differential Equations in variational form. It minimizes the dual norm of the residual, which naturally yields a saddle-point (min-max) problem over the so-called trial and test spaces. In the context of neural networks, we can address this min-max approach by employing one network to seek the trial minimum, while another network seeks the test maximizers. However, the resulting method is numerically unstable as we approach the trial solution. To overcome this, we reformulate the residual minimization as an equivalent minimization of a Ritz functional fed by optimal test functions computed from another Ritz functional minimization. We call the resulting scheme the Deep Double Ritz Method (D2RM), which combines two neural networks for approximating trial functions and optimal test functions along a nested double Ritz minimization strategy. Numerical results on different diffusion and convection problems support the robustness of our method, up to the approximation properties of the networks and the training capacity of the optimizers. (c) 2023 The Author(s). Published by Elsevier B.V. This is an open access article under the CC BY-NC-ND license (http://creativecommons.org/licenses/by-nc-nd/4.0/).
引用
收藏
页数:24
相关论文
共 55 条
[1]  
Abadi M., 2016, PREPRINT, DOI DOI 10.48550/ARXIV.1603.04467
[2]   Numerical solution of inverse problems by weak adversarial networks [J].
Bao, Gang ;
Ye, Xiaojing ;
Zang, Yaohua ;
Zhou, Haomin .
INVERSE PROBLEMS, 2020, 36 (11)
[3]  
Baydin AG, 2018, J MACH LEARN RES, V18
[4]  
Bochev PB, 2009, APPL MATH SCI, V166, P3, DOI 10.1007/b13382_1
[5]   A least-squares approach based on a discrete minus one inner product for first order systems [J].
Bramble, JH ;
Lazarov, RD ;
Pasciak, JE .
MATHEMATICS OF COMPUTATION, 1997, 66 (219) :935-955
[6]   Neural control of discrete weak formulations: Galerkin, least squares & minimal-residual methods with quasi-optimal weights [J].
Brevis, Ignacio ;
Muga, Ignacio ;
van der Zee, Kristoffer G. .
COMPUTER METHODS IN APPLIED MECHANICS AND ENGINEERING, 2022, 402
[7]   A machine-learning minimal-residual (ML-MRes) framework for goal-oriented finite element discretizations [J].
Brevis, Ignacio ;
Muga, Ignacio ;
van der Zee, Kristoffer G. .
COMPUTERS & MATHEMATICS WITH APPLICATIONS, 2021, 95 :186-199
[8]   ON THE STABILITY OF DPG FORMULATIONS OF TRANSPORT EQUATIONS [J].
Broersen, D. ;
Dahmen, W. ;
Stevenson, R. P. .
MATHEMATICS OF COMPUTATION, 2018, 87 (311) :1051-1082
[9]   (PARAMETRIZED) FIRST ORDER TRANSPORT EQUATIONS: REALIZATION OF OPTIMALLY STABLE PETROV-GALERKIN METHODS [J].
Brunken, Julia ;
Smetana, Kathrin ;
Urban, Karsten .
SIAM JOURNAL ON SCIENTIFIC COMPUTING, 2019, 41 (01) :A592-A621
[10]   Isogeometric Residual Minimization Method (iGRM) with direction splitting preconditioner for stationary advection-dominated diffusion problems [J].
Calo, V. M. ;
Los, M. ;
Deng, Q. ;
Muga, I. ;
Paszynski, M. .
COMPUTER METHODS IN APPLIED MECHANICS AND ENGINEERING, 2021, 373 (373)