Inverse modelling of an aneurysm's stiffness using surrogate-based optimization of a three-dimensional fluid-structure interaction simulation

Characterization of the mechanical properties of arterial tissues is highly relevant, both as a direct diagnostic tool as well indirectly to provide material parameters for patientspecific computer simulations. This is especially true in the context of (aortic) aneurysms. In this work, we apply an inverse modelling approach to a model accounting for an aneurysm and the distal part of the circulation. The stiffness of the aneurysm itself and of the rest of the arterial system are modelled using two independent stiffness parameters. For given values of these parameters, the position of the arterial wall as a function of time is calculated using a forward simulation which takes the fluid-structure interaction (FSI) between the blood flow and the arterial wall into account. Using this forward simulation, the correct values of the stiffness parameters are obtained by minimizing a cost function, which is defined as the difference between the forward simulation and a measurement (in this case synthetic data from a simulation). The minimization is performed by means of surrogate-based optimization using a Kriging model combined with the expected improvement infill criterion. Results for a three-dimensional fluid-structure interaction model of the aneurysm show that the stiffness parameters converge to the correct values.

[1]  N. Cressie.  The origins of kriging , 1990 .

[2]  C. D. Perttunen,et al.  Lipschitzian optimization without the Lipschitz constant , 1993 .

[3]  Thomas J. Santner,et al.  Design and analysis of computer experiments , 1998 .

[4]  K. Bathe,et al.  Performance of a new partitioned procedure versus a monolithic procedure in fluid-structure interaction , 2009 .

[5]  Pierre-Yves Lagrée,et al.  An inverse technique to deduce the elasticity of a large artery , 2000 .

[6]  Fabio Nobile,et al.  Added-mass effect in the design of partitioned algorithms for fluid-structure problems , 2005 .

[7]  Andy J. Keane,et al.  Combining Global and Local Surrogate Models to Accelerate Evolutionary Optimization , 2007, IEEE Transactions on Systems, Man, and Cybernetics, Part C (Applications and Reviews).

[8]  Wolfgang A. Wall,et al.  Coupling strategies for biomedical fluid–structure interaction problems , 2010 .

[9]  B. Minasny,et al.  The Matérn function as a general model for soil variograms , 2005 .

[10]  P ? ? ? ? ? ? ? % ? ? ? ? , 1991 .

[11]  Donald R. Jones,et al.  Efficient Global Optimization of Expensive Black-Box Functions , 1998, J. Glob. Optim..

[12]  A. Keane,et al.  Evolutionary Optimization of Computationally Expensive Problems via Surrogate Modeling , 2003 .

[13]  Said Elias Francis,et al.  Continuous estimation of cardiac output and arterial resistance from arterial blood pressure using a third-order Windkessel model , 2007 .

[14]  Matteo Astorino,et al.  Fluid-structure interaction and multi-body contact. Application to aortic valves , 2009 .

[15]  N. M. Alexandrov,et al.  A trust-region framework for managing the use of approximation models in optimization , 1997 .

[16]  Jan Vierendeels,et al.  Stability analysis of Gauss-Seidel iterations in a partitioned simulation of fluid-structure interaction , 2010 .

[17]  R. Haftka.  Combining global and local approximations , 1991 .

[18]  Jan Vierendeels,et al.  Stability of a coupling technique for partitioned solvers in FSI applications , 2008 .

[19]  Piet Demeester,et al.  A Surrogate Modeling and Adaptive Sampling Toolbox for Computer Based Design , 2010, J. Mach. Learn. Res..

[20]  Wolfgang A. Wall,et al.  A computational strategy for prestressing patient‐specific biomechanical problems under finite deformation , 2010 .

[21]  Dick den Hertog,et al.  Maximin Latin Hypercube Designs in Two Dimensions , 2007, Oper. Res..

[22]  Jean-Frédéric Gerbeau,et al.  Parameter identification for a one-dimensional blood flow model , 2005 .

[23]  Hendrik Rogier,et al.  Surrogate-based infill optimization applied to electromagnetic problems , 2010 .

[24]  D. Duprez,et al.  A Comparison Between Systolic and Diastolic Pulse Contour Analysis in the Evaluation of Arterial Stiffness , 2001, Hypertension.