Simone Scacchi

Paola F. Antonietti, Lourenco Beirao da Veiga, Simone Scacchi et al.

In this paper we develop an evolution of the $C^1$ virtual elements of minimal degree for the approximation of the Cahn-Hilliard equation. The proposed method has the advantage of being conforming in $H^2$ and making use of a very simple set of degrees of freedom, namely 3 degrees of freedom per vertex of the mesh. Moreover, although the present method is new also on triangles, it can make use of general polygonal meshes. As a theoretical and practical support, we prove the convergence of the semi-discrete scheme and investigate the performance of the fully discrete scheme through a set of numerical tests.

Paola F. Antonietti, Matteo Bruggi, Simone Scacchi et al.

It is well known that the solution of topology optimization problems may be affected both by the geometric properties of the computational mesh, which can steer the minimization process towards local (and non-physical) minima, and by the accuracy of the method employed to discretize the underlying differential problem, which may not be able to correctly capture the physics of the problem. In light of the above remarks, in this paper we consider polygonal meshes and employ the virtual element method (VEM) to solve two classes of paradigmatic topology optimization problems, one governed by nearly-incompressible and compressible linear elasticity and the other by Stokes equations. Several numerical results show the virtues of our polygonal VEM based approach with respect to more standard methods.

Edoardo Centofanti, Giovanni Ziarelli, Simone Scacchi et al.

Solving inverse problems in cardiac electrophysiology consists in the recovery of physiological parameters from surface electrocardiogram (ECG) measurements, a task which is often computationally unfeasible due to the severe ill-posedness and the prohibitive computational complexity of PDE-constrained optimization. In this work, we introduce a data-driven framework leveraging Latent Dynamics Networks (LDNets) to construct efficient surrogate models of the forward operator. By mapping low-dimensional parameters, representing ectopic activation sites or ischemic region descriptors, to the ECG signals via latent dynamics governed by neural ordinary differential equations, our approach circumvents the computational burden of evaluating high-fidelity cardiac models during iterative parameter estimation. The surrogate is trained offline on high-fidelity data, enabling rapid and robust inversion. We validate the proposed framework through rigorous numerical experiments with synthetic data across both 2d and 3d geometries. Results show that the LDNet-based surrogate achieves precise reconstruction of cardiac parameters while drastically reducing computational overhead, thereby enabling near real-time clinical applications.