Luca Gerardo-Giorda

Nicole Cusimano, Félix del Teso, Luca Gerardo-Giorda et al.

In this work, we propose novel discretizations of the spectral fractional Laplacian on bounded domains based on the integral formulation of the operator via the heat-semigroup formalism. Specifically, we combine suitable quadrature formulas of the integral with a finite element method for the approximation of the solution of the corresponding heat equation. We derive two families of discretizations with order of convergence depending on the regularity of the domain and the function on which the spectral fractional Laplacian is acting. Our method does not require the computation of the eigenpairs of the Laplacian on the considered domain, can be implemented on possibly irregular bounded domains, and can naturally handle different types of boundary constraints. Various numerical simulations are provided to illustrate performance of the proposed method and support our theoretical results.

Argyrios Petras, Massimiliano Leoni, Jose M. Guerra et al.

Radiofrequency catheter ablation (RFCA) is an effective treatment for cardiac arrhythmias. Although generally safe, it is not completely exempt from the risk of complications. The great flexibility of computational models can be a major asset in optimizing interventional strategies, if they can produce sufficiently precise estimations of the generated lesion for a given ablation protocol. This requires an accurate description of the catheter tip and the cardiac tissue. In particular, the deformation of the tissue under the catheter pressure during the ablation is an important aspect that is overlooked in the existing literature, that resorts to a sharp insertion of the catheter into an undeformed geometry. As the lesion size depends on the power dissipated in the tissue, and the latter depends on the percentage of the electrode surface in contact with the tissue itself, the sharp insertion geometry has the tendency to overestimate the lesion obtained, especially when a larger force is applied to the catheter. In this paper we introduce a full 3D computational model that takes into account the tissue elasticity, and is able to capture the tissue deformation and realistic power dissipation in the tissue. Numerical results in FEniCS-HPC are provided to validate the model against experimental data, and to compare the lesions obtained with the new model and with the classical ones featuring a sharp electrode insertion in the tissue.

Victorita Dolean, Martin Gander, Luca Gerardo-Giorda

Over the last two decades, classical Schwarz methods have been extended to systems of hyperbolic partial differential equations, and it was observed that the classical Schwarz method can be convergent even without overlap in certain cases. This is in strong contrast to the behavior of classical Schwarz methods applied to elliptic problems, for which overlap is essential for convergence. Over the last decade, optimized Schwarz methods have been developed for elliptic partial differential equations. These methods use more effective transmission conditions between subdomains, and are also convergent without overlap for elliptic problems. We show here why the classical Schwarz method applied to the hyperbolic problem converges without overlap for Maxwell's equations. The reason is that the method is equivalent to a simple optimized Schwarz method for an equivalent elliptic problem. Using this link, we show how to develop more efficient Schwarz methods than the classical ones for the Maxwell's equations. We illustrate our findings with numerical results.