140 research outputs found
IGS: an IsoGeometric approach for Smoothing on surfaces
We propose an Isogeometric approach for smoothing on surfaces, namely
estimating a function starting from noisy and discrete measurements. More
precisely, we aim at estimating functions lying on a surface represented by
NURBS, which are geometrical representations commonly used in industrial
applications. The estimation is based on the minimization of a penalized
least-square functional. The latter is equivalent to solve a 4th-order Partial
Differential Equation (PDE). In this context, we use Isogeometric Analysis
(IGA) for the numerical approximation of such surface PDE, leading to an
IsoGeometric Smoothing (IGS) method for fitting data spatially distributed on a
surface. Indeed, IGA facilitates encapsulating the exact geometrical
representation of the surface in the analysis and also allows the use of at
least globally continuous NURBS basis functions for which the 4th-order
PDE can be solved using the standard Galerkin method. We show the performance
of the proposed IGS method by means of numerical simulations and we apply it to
the estimation of the pressure coefficient, and associated aerodynamic force on
a winglet of the SOAR space shuttle
Effect of fibre orientation and bulk modulus on the electromechanical modelling of human ventricles
This work concerns the mathematical and numerical modeling of the heart. The aim is to enhance the understanding of the cardiac function in both physiological and pathological conditions. Along this road, a challenge arises from the multi-scale and multi-physics nature of the mathematical problem at hand. In this paper, we propose an electromechanical model that, in bi-ventricle geometries, combines the monodomain equation, the Bueno-Orovio minimal ionic model, and the Holzapfel-Ogden strain energy function for the passive myocardial tissue modelling together with the active strain approach combined with a model for the transmurally heterogeneous thickening of the myocardium. Since the distribution of the electric signal is dependent on the fibres orientation of the ventricles, we use a Laplace-Dirichlet Rule-Based algorithm to determine the myocardial fibres and sheets configuration in the whole bi-ventricle. In this paper, we study the influence of different fibre directions and incompressibility constraint and penalization on the compressibility of the material (bulk modulus) on the pressure-volume relation simulating a full heart beat. The coupled electromechanical problem is addressed by means of a fully segregated scheme. The numerical discretization is based on the Finite Element Method for the spatial discretization and on Backward Differentiation Formulas for the time discretization. The arising non-linear algebraic system coming from application of the implicit scheme is solved through the Newton method. Numerical simulations are carried out in a patient-specific biventricle geometry to highlight the most relevant results of both electrophysiology and mechanics and to compare them with physiological data and measurements. We show how various fibre configurations and bulk modulus modify relevant clinical quantities such as stroke volume, ejection fraction and ventricle contractility
Deep learning-based reduced order models in cardiac electrophysiology
Predicting the electrical behavior of the heart, from the cellular scale to
the tissue level, relies on the formulation and numerical approximation of
coupled nonlinear dynamical systems. These systems describe the cardiac action
potential, that is the polarization/depolarization cycle occurring at every
heart beat that models the time evolution of the electrical potential across
the cell membrane, as well as a set of ionic variables. Multiple solutions of
these systems, corresponding to different model inputs, are required to
evaluate outputs of clinical interest, such as activation maps and action
potential duration. More importantly, these models feature coherent structures
that propagate over time, such as wavefronts. These systems can hardly be
reduced to lower dimensional problems by conventional reduced order models
(ROMs) such as, e.g., the reduced basis (RB) method. This is primarily due to
the low regularity of the solution manifold (with respect to the problem
parameters) as well as to the nonlinear nature of the input-output maps that we
intend to reconstruct numerically. To overcome this difficulty, in this paper
we propose a new, nonlinear approach which exploits deep learning (DL)
algorithms to obtain accurate and efficient ROMs, whose dimensionality matches
the number of system parameters. Our DL approach combines deep feedforward
neural networks (NNs) and convolutional autoencoders (AEs). We show that the
proposed DL-ROM framework can efficiently provide solutions to parametrized
electrophysiology problems, thus enabling multi-scenario analysis in
pathological cases. We investigate three challenging test cases in cardiac
electrophysiology and prove that DL-ROM outperforms classical projection-based
ROMs.Comment: 28 page
La morfologia derivazionale tra analisi tradizionali e prospettive teoriche moderne
Dopo una breve panoramica dei percorsi degli studi di morfologia derivazionale all’interno della tradizione di linguistica indoeuropea, si prende in esame l’interazione tra questi studi e i modelli teorici e interpretativi sorti nell’ambito della linguistica teorica e tipologica. Si mostra come tale interazione sia proficua da un lato perché consente una migliore interpretazione dei dati offerti dalle lingue storiche, dall’altro perché accresce la consapevolezza delle caratteristiche peculiari del processo di ricostruzione della protolingua indoeuropea.After a brief overview of the studies in derivative morphology within the tradition of Indo-European linguistics, we examine the interaction between these studies and the theoretical and interpretative models that have emerged in the field of theoretical and typological linguistics. It is shown that this interaction is fruitful on the one hand because it allows a better interpretation of the data offered by historical languages, and on the other because it increases awareness of the peculiar characteristics of the process of reconstruction of the Indo-European protolanguage
Premessa
I contributi raccolti in questa sezione monografica I processi di derivazione tra linguistica indoeuropea e linguistica generale nascono da comunicazioni presentate durante una giornata di studio dall’omonimo titolo tenutasi il 17 giugno 2019 presso l’Università degli Studi di Milano e inquadrata all’interno di un progetto di ricerca diretto da chi scrive e finanziato dal Dipartimento di Studi letterari, filologici e linguistici
Mathematical analysis and numerical approximation of a general linearized poro-hyperelastic model
Abstract We describe the behavior of a deformable porous material by means of a poro-hyperelastic model that has been previously proposed in Chapelle and Moireau (2014) under general assumptions for mass and momentum balance and isothermal conditions for a two-component mixture of fluid and solid phases. In particular, we address here a linearized version of the model, based on the assumption of small displacements. We consider the mathematical analysis and the numerical approximation of the problem. More precisely, we carry out firstly the well-posedness analysis of the model. Then, we propose a numerical discretization scheme based on finite differences in time and finite elements for the spatial approximation; stability and numerical error estimates are proved. Particular attention is dedicated to the study of the saddle-point structure of the problem, that turns out to be interesting because velocities of the fluid phase and of the solid phase are combined into a single quasi-incompressibility constraint. Our analysis provides guidelines to select the componentwise polynomial degree of approximation of fluid velocity, solid displacement and pressure, to obtain a stable and robust discretization based on Taylor–Hood type finite element spaces. Interestingly, we show how this choice depends on the porosity of the mixture, i.e. the volume fraction of the fluid phase
Modeling cardiac muscle fibers in ventricular and atrial electrophysiology simulations
Since myocardial fibers drive the electric signal propagation throughout the
myocardium, accurately modeling their arrangement is essential for simulating
heart electrophysiology (EP). Rule-Based-Methods (RBMs) represent a commonly
used strategy to include cardiac fibers in computational models. A particular
class of such methods is known as Laplace-Dirichlet-Rule-Based-Methods (LDRBMs)
since they rely on the solution of Laplace problems. In this work we provide a
unified framework, based on LDRBMs, for generating full heart muscle fibers.
First, we review existing ventricular LDRBMs providing a communal mathematical
description and introducing also some modeling improvements with respect to the
existing literature. We then carry out a systematic comparison of LDRBMs based
on meaningful biomarkers produced by numerical EP simulations. Next we propose,
for the first time, a LDRBM to be used for generating atrial fibers. The new
method, tested both on idealized and realistic atrial models, can be applied to
any arbitrary geometries. Finally, we present numerical results obtained in a
realistic whole heart where fibers are included for all the four chambers using
the discussed LDRBMs
- …