We establish estimates on the error made by the Ritz method for quadratic energies on the space H 1 (Ω) in the approximation of the solution of variational problems with different boundary conditions. Special attention is paid to the case of Dirichlet boundary values which are treated with the boundary penalty method. We consider arbitrary and in general non linear classes V ⊆ H 1 (Ω) of ansatz functions and estimate the error in dependence of the optimisation accuracy, the approximation capabilities of the ansatz class and -in the case of Dirichlet boundary values -the penalisation strength λ. For non-essential boundary conditions the error of the Ritz method decays with the same rate as the approximation rate of the ansatz classes. For the boundary penalty method we obtain that given an approximation rate of r in H 1 (Ω) and an approximation rate of s in L 2 (∂Ω) of the ansatz classes, the optimal decay rate of the estimated error is min(s/2, r) ∈ [r/2, r] and achieved by choosing λn ∼ n s . We discuss how this rate can be improved, the relation to existing estimates for finite element functions as well as the implications for ansatz classes which are given through ReLU networks. Finally, we use the notion of Γ-convergence to show that the Ritz method converges for a wide class of energies including nonlinear stationary PDEs like the p-Laplace.
We provide convergence guarantees for the Deep Ritz Method for abstract variational energies. Our results cover non-linear variational problems such as the p-Laplace equation or the Modica-Mortola energy with essential or natural boundary conditions. Under additional assumptions, we show that the convergence is uniform across bounded families of right-hand sides.
We analyse the difference in convergence mode using exact versus penalised boundary values for the residual minimisation of PDEs with neural network type ansatz functions, as is commonly done in the context of physics informed neural networks. It is known that using an L 2 boundary penalty leads to a loss of regularity of 3/2 meaning that approximation in H 2 yields a priori estimates in H 1/2 . These notes demonstrate how this loss of regularity can be circumvented if the functions in the ansatz class satisfy the boundary values exactly. Furthermore, it is shown that in this case, the loss function provides a consistent a posteriori error estimator in H 2 norm made by the residual minimisation method. We provide analogue results for linear time dependent problems and discuss the implications of measuring the residual in Sobolev norms.
We compare different training strategies for the Deep Ritz Method for elliptic equations with Dirichlet boundary conditions and highlight the problems arising from the boundary values. We distinguish between an exact resolution of the boundary values by introducing a distance function and the approximation through a Robin Boundary Value problem. However, distance functions are difficult to obtain for complex domains. Therefore, it is more feasible to solve a Robin Boundary Value problem which approximates the solution to the Dirichlet Boundary Value problem, yet the naïve approach to this problem becomes unstable for large penalizations. A novel method to compensate this problem is proposed using a small penalization strength to pre-train the model before the main training on the target penalization strength is conducted. We present numerical and theoretical evidence that the proposed method is beneficial.
We present a three dimensional, time dependent model for bone regeneration in the presence of porous scaffolds to bridge critical size bone defects. Our approach uses homogenized quantities, thus drastically reducing computational cost compared to models resolving the microstructural scale of the scaffold. Using abstract functional relationships instead of concrete effective material properties, our model can incorporate the homogenized material tensors for a large class of scaffold microstructure designs. We prove an existence and uniqueness theorem for solutions based on a fixed point argument. We include the cases of mixed boundary conditions and multiple, interacting signalling molecules, both being important for application. Furthermore we present numerical simulations showing good agreement with experimental findings.
scite is a Brooklyn-based organization that helps researchers better discover and understand research articles through Smart Citations–citations that display the context of the citation and describe whether the article provides supporting or contrasting evidence. scite is used by students and researchers from around the world and is funded in part by the National Science Foundation and the National Institute on Drug Abuse of the National Institutes of Health.
hi@scite.ai
10624 S. Eastern Ave., Ste. A-614
Henderson, NV 89052, USA
Copyright © 2024 scite LLC. All rights reserved.
Made with 💙 for researchers
Part of the Research Solutions Family.