Home  << 1 >> 
Aylwin, R., & JerezHanckes, C. (2023). FiniteElement Domain Approximation for Maxwell Variational Problems on Curved Domains. SIAM J. Numer. Anal., 61(3), 1139–1171.
Abstract: We consider the problem of domain approximation in finite element methods for Maxwell equations on curved domains, i.e., when affine or polynomial meshes fail to exactly cover the domain of interest. In such cases, one is forced to approximate the domain by a sequence of polyhedral domains arising from inexact meshes. We deduce conditions on the quality of these approximations that ensure rates of error convergence between discrete solutions  in the approximate domains  to the continuous one in the original domain.

Aylwin, R., JerezHanckes, C., Schwab, C., & Zech, J. (2020). Domain Uncertainty Quantification in Computational Electromagnetics. SIAMASA J. Uncertain. Quantif., 8(1), 301–341.
Abstract: We study the numerical approximation of timeharmonic, electromagnetic fields inside a lossy cavity of uncertain geometry. Key assumptions are a possibly highdimensional parametrization of the uncertain geometry along with a suitable transformation to a fixed, nominal domain. This uncertainty parametrization results in families of countably parametric, Maxwelllike cavity problems that are posed in a single domain, with inhomogeneous coefficients that possess finite, possibly low spatial regularity, but exhibit holomorphic parametric dependence in the differential operator. Our computational scheme is composed of a sparse grid interpolation in the highdimensional parameter domain and an Hcurl conforming edge element discretization of the parametric problem in the nominal domain. As a steppingstone in the analysis, we derive a novel Strangtype lemma for Maxwelllike problems in the nominal domain, which is of independent interest. Moreover, we accommodate arbitrary small Sobolev regularity of the electric field and also cover uncertain isotropic constitutive or material laws. The shape holomorphy and edgeelement consistency error analysis for the nominal problem are shown to imply convergence rates for multilevel Monte Carlo and for quasiMonte Carlo integration, as well as sparse grid approximations, in uncertainty quantification for computational electromagnetics. They also imply expression rate estimates for deep ReLU networks of shapetosolution maps in this setting. Finally, our computational experiments confirm the presented theoretical results.

Duran, M., Godoy, E., RomanCatafau, E., & Toledo, P. A. (2022). Openpit slope design using a DtNFEM: Parameter space exploration. Int. J. Rock Mech. Min. Sci., 149, 104950.
Abstract: Given the sustained mineraldeposits oregrade decrease, it becomes necessary to reach greater depths when extracting ore by openpit mining. Steeper slope angles are thus likely to be required, leading to geomechanical instabilities. In order to determine excavation stability, mathematical modelling and numerical simulation are often used to compute the rockmass stressstate, to which some stability criterion needs to be added. A problem with this approach is that the volume surrounding the excavation has no clear borders and in practice it might be regarded as an unbounded region. Then, it is necessary to use advanced methods capable of dealing efficiently with this difficulty. In this work, a DtNFEM procedure is applied to calculate displacements and stresses in openpit slopes under geostatic stress conditions. This procedure was previously devised by the authors to numerically treat this kind of problems where the surrounding domain is semiinfinite. Its efficiency makes possible to simulate, in a short amount of time, multiple openpit slope configurations. Therefore, the method potentiality for openpit slope design is investigated. A regular openpit slope geometry is assumed, parameterised by the overallslope and benchface angles. Multiple geometrically admissible slopes are explored and their stability is assessed by using the computed stressfield and the MohrCoulomb failure criterion. Regions of stability and instability are thus explored in the parametric space, opening the way for a new and flexible designing tool for openpit slopes and related problems.
Keywords: DirichlettoNeumann map; Finite elements; Openpit; Slope design
