Search

In this paper, we demonstrate for the first time how the Integrated Finite Element Neural Network (I-FENN) framework, previously proposed by the authors [1,2] , can efficiently simulate the entire loading history of non-local gradient damage propagation. To achieve this goal, we first adopt a Temporal Convolutional Network (TCN) as the neural network of choice to capture the history-dependent evolution of the non-local strain in a coarsely meshed domain. The quality of the network predictions governs the computational performance of I-FENN, and therefore we perform an extended investigation aimed at enhancing them. We explore a data-driven vs. physics-informed TCN setup to arrive at an optimum network training, evaluating the network based on a coherent set of relevant performance metrics. We address the crucial issue of training a physics-informed network with input data that span vastly different length scales by proposing a systematic way of input normalization and output un-normalization. We then integrate the trained TCN within the nonlinear iterative FEM solver and apply I-FENN to simulate the damage propagation analysis. I-FENN is always applied in mesh idealizations different from the one used for the TCN training, showcasing the framework's ability to be used at progressively refined mesh resolutions. We illustrate several cases that I-FENN completes the simulation using either a modified or a full Newton–Raphson scheme, and we showcase its computational savings compared to both the classical monolithic and staggered FEM solvers. We underline that we satisfy very strict convergence criteria for every increment across the entire simulation, providing clear evidence of the robustness and accuracy of I-FENN. All the code and data used in this work will be made publicly available upon publication of the article. • Non-local damage propagation is modeled for the first time with a TCN-based I-FENN. • I-FENN is faster than monolithic and staggered FEM and savings scale with model size. • I-FENN's robustness is shown by satisfying very strict convergence criteria. • Input/output normalization methods are applied to enhance the network training. • I-FENN is always applied in different training and testing mesh idealizations. [ABSTRACT FROM AUTHOR]

The numerical simulation of cardiac electrophysiology is a highly challenging problem in scientific computing. The Bidomain system is the most complete mathematical model of cardiac bioelectrical activity. It consists of an elliptic and a parabolic partial differential equation (PDE), of reaction–diffusion type, describing the spread of electrical excitation in the cardiac tissue. The two PDEs are coupled with a stiff system of ordinary differential equations (ODEs), representing ionic currents through the cardiac membrane. Developing efficient and scalable preconditioners for the linear systems arising from the discretization of such computationally challenging model is crucial in order to reduce the computational costs required by the numerical simulations of cardiac electrophysiology. In this work, focusing on the Bidomain system as a model problem, we have benchmarked two popular implementations of the Algebraic Multigrid (AMG) preconditioner embedded in the PETSc library and we have studied the performance on the calibration of specific parameters. We have conducted our analysis on modern HPC architectures, performing scalability tests on multi-core and multi-GPUs settings. The results have shown that, for our problem, although scalability is verified on CPUs, GPUs are the optimal choice, since they yield the best performance in terms of solution time. [ABSTRACT FROM AUTHOR]

Interface problems have long been a major focus of scientific computing, leading to the development of various numerical methods. Traditional mesh-based methods often employ time-consuming body-fitted meshes with standard discretization schemes or unfitted meshes with tailored schemes to achieve controllable accuracy and convergence rate. Along another line, mesh-free methods bypass mesh generation but lack robustness in terms of convergence and accuracy due to the low regularity of solutions. In this study, we propose a novel method for solving interface problems within the framework of the random feature method (RFM). This approach utilizes random feature functions in conjunction with a partition of unity as approximation functions, and solves a linear least-squares system to obtain the approximate solution. In the context of interface problems, two innovative and crucial components are incorporated into the RFM. Firstly, we utilize two sets of random feature functions on each side of the interface, allowing for the inclusion of low regularity or even discontinuous behaviors in the solution. Secondly, the construction of the loss function is based on the assessment of the partial differential equation, initial/boundary conditions, and the interface condition on collocation points. This approach ensures that these conditions are equally satisfied. Consequently, the challenges arising from geometric complexity primarily manifest in the generation of collocation points, a task amenable to standard methods. Importantly, the proposed method retains its meshfree characteristics and robustness when addressing problems featuring intricate geometries. We validate our method through a series of linear interface problems with increasingly complex geometries, including two-dimensional elliptic and three-dimensional Stokes interface problems, a three-dimensional elasticity interface problem, a moving interface problem with topological change, a dynamic interface problem with large deformation, and a linear fluid–solid interaction problem with complex geometry. Our findings show that despite the solution often being only continuous or even discontinuous, our method not only eliminates the need for mesh generation but also maintains high accuracy, akin to the spectral collocation method for smooth solutions. Remarkably, for the same accuracy requirement, our method requires two to three orders of magnitude fewer degrees of freedom than traditional methods, demonstrating its significant potential for solving interface problems with complex geometries and predetermined intricate evolution. [ABSTRACT FROM AUTHOR]