174 research outputs found
Parameter identification for a one-dimensional blood flow model
International audienceThe purpose of this work is to use a variational method to identify some of the parameters of one-dimensional models for blood flow in arteries. These parameters can be fit to approach as much as possible some data coming from experimental measurements or from numerical simulations performed using more complex models. A nonlinear least squares approach to parameter estimation was taken, based on the optimization of a cost function. The resolution of such an optimization problem generally requires the efficient and accurate computation of the gradient of the cost function with respect to the parameters. This gradient is computed analytically when the one-dimensional hyperbolic model is discretized with a second order Taylor-Galerkin scheme. An adjoint approach was used. Some preliminary numerical tests are shown. In these simulations, we mainly focused on determining a parameter that is linked to the mechanical properties of the arterial walls, the compliance. The synthetic data we used to estimate the parameter were obtained from a numerical computation performed with a more accurate model: a three-dimensional fluid-structure interaction model. The first results seem to be promising. In particular, it is worth noticing that the estimated compliance which gives the best fit is quite different from the values that are commonly used in practice.Le but de ce travail est d'identifier certains des paramètres existant dans des modèles 1-d d'écoulement sanguin dans des artères. Ces paramètres peuvent permettre d'approcher autant que possible des configurations géométriques réalistes ou des données expérimentales. Une approche de l'estimation de paramètres par moindres carrés non-linéaires a été adoptée, basée sur l'optimisation d'une certaine fonction coût. La résolution d'un tel problème de minimisation requiert le calcul efficace et précis du gradient de la fonction coût par rapport aux paramètres. Le gradient est discrétisé analytiquement dans le cas d'une discrétisation du modèle hyperbolique 1-d par le schema de Taylor-Galerkin. Une approche par l'état adjoint a été employée. Des premiers résultats numériques sont fournis. Pour ces simulations, nous nous sommes concentrés sur la détermination d'un paramètre lié aux propriétés mécaniques de la paroi artérielle. Les données synthétiques utilisées pour l'estimation de ce paramètre ont été obtenues à partir d'un modèle beaucoup plus raffiné : un modèle 3-d d'interaction fluide-structure. Les résultats semblent intéressants car le paramètre estimé est assez différent de ce à quoi on s'attendrait a priori
Numerical method of characteristics for one-dimensional blood flow
Mathematical modeling at the level of the full cardiovascular system requires
the numerical approximation of solutions to a one-dimensional nonlinear
hyperbolic system describing flow in a single vessel. This model is often
simulated by computationally intensive methods like finite elements and
discontinuous Galerkin, while some recent applications require more efficient
approaches (e.g. for real-time clinical decision support, phenomena occurring
over multiple cardiac cycles, iterative solutions to optimization/inverse
problems, and uncertainty quantification). Further, the high speed of pressure
waves in blood vessels greatly restricts the time step needed for stability in
explicit schemes. We address both cost and stability by presenting an efficient
and unconditionally stable method for approximating solutions to diagonal
nonlinear hyperbolic systems. Theoretical analysis of the algorithm is given
along with a comparison of our method to a discontinuous Galerkin
implementation. Lastly, we demonstrate the utility of the proposed method by
implementing it on small and large arterial networks of vessels whose elastic
and geometrical parameters are physiologically relevant
Comparison of reduced models for blood flow using Runge-Kutta discontinuous Galerkin methods
One-dimensional blood flow models take the general form of nonlinear
hyperbolic systems but differ greatly in their formulation. One class of models
considers the physically conserved quantities of mass and momentum, while
another class describes mass and velocity. Further, the averaging process
employed in the model derivation requires the specification of the axial
velocity profile; this choice differentiates models within each class.
Discrepancies among differing models have yet to be investigated. In this
paper, we systematically compare several reduced models of blood flow for
physiologically relevant vessel parameters, network topology, and boundary
data. The models are discretized by a class of Runge-Kutta discontinuous
Galerkin methods
High-order fully well-balanced numerical methods for one-dimensional blood flow with discontinuous properties
In this paper, we are interested in the numerical study of the one-dimensional blood flow model with discontinuous mechanical and geometrical properties. We present the mathematical model together with its nondimensional form. We do an exhaustive investigation of all its stationary solutions and we propose high-order fully well-balanced numerical methods that are able to preserve all of them. They are based on the combination of the Generalized Hydrostatic Reconstruction and well-balanced reconstruction operators. These methods are able to deal with more than one discontinuous parameter. Several numerical tests are shown to prove its well-balanced and high-order properties, and its convergence to the exact solutions.The research of EPG and CP was partially supported by the Spanish Government (SG), the European Regional Development Fund (ERDF), the Regional Government of Andalusia (RGA), and the University of Málaga (UMA) through the projects of reference RTI2018-096064-B-C21 (SG-ERDF), UMA18-Federja-161 (RGA-ERDF-UMA), and P18-RT-3163 (RGA-ERDF). EPG was also financed by the European Union – NextGenerationEU. // Funding for open access charge: Universidad de Málaga / CBU
Assessment of reduced order Kalman filter for parameter identification in one-dimensional blood flow models using experimental data
This work presents a detailed investigation of
a parameter estimation approach based on the
reduced order unscented Kalman filter (ROUKF) in the context of one-dimensional blood flow models.
In particular, the main aims of this study are (i) to investigate the effect of
using real measurements vs. synthetic data (i.e., numerical results of the same in silico model, perturbed with white
noise) for the estimation and (ii) to identify
potential difficulties and limitations of the approach in clinically realistic applications in order
to assess the applicability of the filter to such setups.
For these purposes, our numerical study is based on the in vitro model of the arterial network
described by [Alastruey et al. 2011, J. Biomech. {\bf 44}], for
which experimental
flow and pressure measurements are available at few selected locations.
In order to mimic clinically relevant situations,
we focus on the estimation of terminal resistances and arterial wall parameters related to vessel mechanics
(Young's modulus and thickness) using few experimental observations (at most a single
pressure or flow measurement per
vessel).
In all cases, we first perform a theoretical identifiability analysis based on the generalized sensitivity
function, comparing then the results obtained with the ROUKF,
using either synthetic or experimental data, to results obtained using reference parameters and to available
measurements
Recommended from our members
Assessment of reduced order Kalman filter for parameter identification in one-dimensional blood flow models using experimental data
This work presents a detailed investigation of a parameter estimation
approach based on the reduced order unscented Kalman filter (ROUKF) in the
context of one-dimensional blood flow models. In particular, the main aims of
this study are (i) to investigate the effect of using real measurements vs.
synthetic data (i.e., numerical results of the same in silico model,
perturbed with white noise) for the estimation and (ii) to identify potential
difficulties and limitations of the approach in clinically realistic
applications in order to assess the applicability of the filter to such
setups. For these purposes, our numerical study is based on the in vitro
model of the arterial network described by [Alastruey et al. 2011, J.
Biomech. 44], for which experimental flow and pressure measurements are
available at few selected locations. In order to mimic clinically relevant
situations, we focus on the estimation of terminal resistances and arterial
wall parameters related to vessel mechanics (Youngs modulus and thickness)
using few experimental observations (at most a single pressure or flow
measurement per vessel). In all cases, we first perform a theoretical
identifiability analysis based on the generalized sensitivity function,
comparing then the results obtained with the ROUKF, using either synthetic or
experimental data, to results obtained using reference parameters and to
available measurements
A novel, FFT-based one-dimensional blood flow solution method for arterial network
In the present work, we propose an FFT-based method for solving blood flow equations in an arterial network with variable properties and geometrical changes. An essential advantage of this approach is in correctly accounting for the vessel skin friction through the use of Womersley solution. To incorporate nonlinear effects, a novel approximation method is proposed to enable calculation of nonlinear corrections. Unlike similar methods available in the literature, the set of algebraic equations required for every harmonic is constructed automatically. The result is a generalized, robust and fast method to accurately capture the increasing pulse wave velocity downstream as well as steepening of the pulse front. The proposed method is shown to be appropriate for incorporating correct convection and diffusion coefficients. We show that the proposed method is fast and accurate and it can be an effective tool for 1D modelling of blood flow in human arterial networks
- …