Preserving the positivity of the deformation gradient determinant in intergrid interpolation by combining RBFs and SVD: Application to cardiac electromechanics

DOI10.1016/J.CMA.2023.116292MaRDI QIDQ6187630zbMATH OpenOpenAlexFDO

Authors Michele Bucelli, Francesco Regazzoni, Luca Dedé, Alfio Quarteroni

Publication date 15 January 2024

Published in Computer Methods in Applied Mechanics and Engineering (Search for Journal in Brave)

Copyright license Creative Commons Attribution 4.0 International

Full work available at URL https://arxiv.org/abs/2306.00457

zbMATH Keywords

singular value decompositionpositivity preservingradial basis function interpolationcardiac modelingmultiphysics modeling

Mathematics Subject Classification ID

Numerical analysis (65-XX) Mechanics of deformable solids (74-XX)

Abstract: The accurate robust and efficient transfer of the deformation gradient tensor between meshes of different resolution is crucial in cardiac electromechanics simulations. We present a novel method that combines rescaled localized Radial Basis Function (RBF) interpolation with Singular Value Decomposition (SVD) to preserve the positivity of the determinant of the deformation gradient tensor. The method involves decomposing the evaluations of the tensor at the quadrature nodes of the source mesh into rotation matrices and diagonal matrices of singular values; computing the RBF interpolation of the quaternion representation of rotation matrices and the singular value logarithms; reassembling the deformation gradient tensors at quadrature nodes of the destination mesh, to be used in the assembly of the electrophysiology model equations. The proposed method overcomes limitations of existing interpolation methods, including nested intergrid interpolation and RBF interpolation of the displacement field, that may lead to the loss of physical meaningfulness of the mathematical formulation and then to solver failures at the algebraic level, due to negative determinant values. The proposed method enables the transfer of solution variables between finite element spaces of different degrees and shapes and without stringent conformity requirements between different meshes, enhancing the flexibility and accuracy of electromechanical simulations. Numerical results confirm that the proposed method enables the transfer of the deformation gradient tensor, allowing to successfully run simulations in cases where existing methods fail. This work provides an efficient and robust method for the intergrid transfer of the deformation gradient tensor, enabling independent tailoring of mesh discretizations to the particular characteristics of the physical components concurring to the of the multiphysics model.

Recommendations

• An intergrid transfer operator using radial basis functions with application to cardiac electromechanics
• Robust radial basis function interpolation based on geodesic distance for the numerical coupling of multiphysics problems
• Efficient mesh deformation based on radial basis function interpolation by means of the inverse fast multipole method
• High-order conservative positivity-preserving DG-interpolation for deforming meshes and application to moving mesh DG simulation of radiative transfer
• Isogeometric mixed collocation of nearly-incompressible electromechanics in finite deformations for cardiac muscle simulations

Cites work

• 3D-0D closed-loop model for the simulation of cardiac biventricular electromechanics
• A cardiac electromechanical model coupled with a lumped-parameter model for closed-loop blood circulation
• A comprehensive and biophysically detailed computational model of the whole human heart electromechanics
• A computationally efficient physiologically comprehensive 3D-0D closed-loop model of the heart and circulation
• A fast cardiac electromechanics model coupling the eikonal and the nonlinear mechanics equations
• A matrix-free high-order solver for the numerical solution of cardiac electrophysiology
• A rescaled localized radial basis function interpolation on non-Cartesian and nonconforming grids
• Advanced Linear Algebra
• An intergrid transfer operator using radial basis functions with application to cardiac electromechanics
• Anatomically accurate high resolution modeling of human whole heart electromechanics: A strongly scalable algebraic multigrid solver method for nonlinear deformation
• Computational model of three-dimensional cardiac electromechanics
• Constitutive modelling of passive myocardium: a structurally based framework for material characterization
• Electro-fluid-mechanics of the heart
• Fast Radial Basis Function Interpolation via Preconditioned Krylov Iteration
• Fast fitting of radial basis functions: Methods based on preconditioned GMRES iteration
• INTERNODES: an accurate interpolation-based method for coupling the Galerkin solutions of PDEs on subdomains featuring non-conforming interfaces
• Integrated heart -- coupling multiscale and multiphysics models for the simulation of the cardiac function
• Isogeometric analysis of the electrophysiology in the human heart: numerical simulation of the bidomain equations on the atria
• Machine learning for fast and reliable solution of time-dependent differential equations
• Machine learning of multiscale active force generation models for the efficient simulation of cardiac electromechanics
• Mathematical cardiac electrophysiology
• Mathematical modelling of the human cardiovascular system. Data, numerical approximation, clinical applications
• Modeling cardiac muscle fibers in ventricular and atrial electrophysiology simulations
• Modeling the cardiac electromechanical function: A mathematical journey
• Multipatch isogeometric analysis for electrophysiology: simulation in a human heart
• On approximate cardinal preconditioning methods for solving PDEs with radial basis functions
• Parallel iterative methods for sparse linear systems
• Piecewise polynomial, positive definite and compactly supported radial functions of minimal degree
• Sensitivity analysis of a strongly-coupled human-based electromechanical cardiac model: effect of mechanical parameters on physiologically relevant biomarkers
• The INTERNODES method for applications in contact mechanics and dedicated preconditioning techniques
• The \textsc{deal.II} finite element library: design, features, and insights
• The \texttt{deal.II} library, version 9.4
• The significant effect of the choice of ionic current integration method in cardiac electro-physiological simulations
• The superconvergent patch recovery anda posteriori error estimates. Part 1: The recovery technique
• Thermodynamically consistent orthotropic activation model capturing ventricular systolic wall thickening in cardiac electromechanics
• p4est: scalable algorithms for parallel adaptive mesh refinement on forests of octrees

Cited in

(1)

• Robust radial basis function interpolation based on geodesic distance for the numerical coupling of multiphysics problems

This page was built for publication: Preserving the positivity of the deformation gradient determinant in intergrid interpolation by combining RBFs and SVD: Application to cardiac electromechanics

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q6187630)