Erik Burman, Johnny Guzman, Manuel A. Sanchez et al.
We prove an optimal error estimate for the flux variable for a stabilized unfitted Nitsche finite element method applied to an elliptic interface problem with discontinuous constant coefficients. Our result shows explicitly that this error estimate is totally independent of the diffusion coefficients
Johnny Guzman, Manuel A. Sanchez, Marcus Sarkis
We define a new finite element method for a steady state elliptic problem with discontinuous diffusion coefficients where the meshes are not aligned with the interface. We prove optimal error estimates in the $L^2$ norm and $H^1$ weighted semi-norm independent of the contrast between the coefficients. Numerical experiments validating our theoretical findings are provided.
Johnny Guzman, Alexandre Madureira, Marcus Sarkis et al.
We derive error estimates for the piecewise linear finite element approximation of the Laplace--Beltrami operator on a bounded, orientable, $C^3$, surface without boundary on general shape regular meshes. As an application, we consider a problem where the domain is split into two regions: one which has relatively high curvature and one that has low curvature. Using a graded mesh we prove error estimates that do not depend on the curvature on the high curvature region. Numerical experiments are provided.
Sijing Liu, Marcus Sarkis, Yi Zhang et al.
We propose a two-scale neural network method for optimal control problems governed by convection-dominated convection-diffusion-reaction equations. Building on two-scale architectures developed for singularly perturbed forward problems, we augment the spatial input with suitably rescaled features that become increasingly important as the diffusion coefficient becomes small. The approach employs separate neural networks for the state and adjoint state variables of the optimality system, reflecting the fact that these quantities develop sharp layers in different parts of the domain due to opposite convection fields. By choosing different center points for the two networks, the architecture naturally aligns with the layer location of each variable. We present two formulations of the method, one based on the first-order optimality conditions and another using penalization of the PDE constraint, and combine them with a successive training strategy that gradually decreases the diffusion coefficient toward its target value. Numerical experiments on benchmark problems illustrate the effectiveness and behavior of the proposed approach.
Kyle Dunn, Roger Lui, Marcus Sarkis
In this paper we introduce a finite element method for the Stokes equations with a massless immersed membrane. This membrane applies normal and tangential forces affecting the velocity and pressure of the fluid. Additionally, the points representing this membrane move with the local fluid velocity. We design and implement a high-accuracy cut finite element method (CutFEM) which enables the use of a structured mesh that is not aligned with the immersed membrane and then we formulate a time discretization that yields an unconditionally energy stable scheme. We prove that the stability is not restricted by the parameter choices that constrained previous finite element immersed boundary methods and illustrate the theoretical results with numerical simulations.
Eduardo Abreu, Ciro Diaz, Juan Galvis et al.
A new high-order conservative finite element method for Darcy flow is presented. The key ingredient in the formulation is a volumetric, residual-based, based on Lagrange multipliers in order to impose conservation of mass that does not involve any mesh dependent parameters. We obtain a method with high-order convergence properties with locally conservative fluxes. Furthermore, our approach can be straightforwardly extended to three dimensions. It is also applicable to highly heterogeneous problems where high-order approximation is preferred.
Johnny Guzman, Manuel A. Sanchez, Marcus Sarkis
We present higher-order piecewise continuous finite element methods for solving a class of interface problems in two dimensions. The method is based on correction terms added to the right-hand side in the standard variational formulation of the problem. We prove optimal error estimates of the methods on general quasi-uniform and shape regular meshes in maximum norms. In addition, we apply the method to a Stokes interface problem, adding correction terms for the velocity and the pressure, obtaining optimal convergence results.