4 results
Search Results
Now showing 1 - 4 of 4
- ItemReduced-order unscented Kalman filter in the frequency domain: Application to computational hemodynamics(Berlin : Weierstraß-Institut für Angewandte Analysis und Stochastik, 2018) Müller, Lucas O.; Caiazzo, Alfonso; Blanco, Pablo J.Objective: The aim of this work is to assess the potential of the reduced order unscented Kalman filter (ROUKF) in the context of computational hemodynamics, in order to estimate cardiovascular model parameters when employing real patient-specific data. Methods: The approach combines an efficient blood flow solver for one-dimensional networks (for the forward problem) with the parameter estimation problem cast in the frequency space. Namely, the ROUKF is used to correct model parameter after each cardiac cycle, depending on the discrepancies of model outputs with respect to available observations properly mapped into the frequency space. Results: First we validate the filter in frequency domain applying it in the context of a set of experimental measurements for an in vitro model. Second, we perform different numerical experiments aiming at parameter estimation using patient-specific data. Conclusion: Our results demonstrate that the filter in frequency domain allows a faster and more robust parameter estimation, when compared to its time domain counterpart. Moreover, the proposed approach allows to estimate parameters that are not directly related to the network but are crucial for targeting inter-individual parameter variability (e.g., parameters that characterize the cardiac output). Significance: The ROUKF in frequency domain provides a robust and flexible tool for estimating parameters related to cardiovascular mathematical models using in vivo data.
- ItemMultiscale coupling of one-dimensional vascular models and elastic tissues(Berlin : Weierstraß-Institut für Angewandte Analysis und Stochastik, 2021) Heltai, Luca; Caiazzo, Alfonso; Müller, Lucas O.We present a computational multiscale model for the efficient simulation of vascularized tissues, composed of an elastic three-dimensional matrix and a vascular network. The effect of blood vessel pressure on the elastic tissue is surrogated via hyper-singular forcing terms in the elasticity equations, which depend on the fluid pressure. In turn, the blood flow in vessels is treated as a one-dimensional network. The pressure and velocity of the blood in the vessels are simulated using a high-order finite volume scheme, while the elasticity equations for the tissue are solved using a finite element method. This work addresses the feasibility and the potential of the proposed coupled multiscale model. In particular, we assess whether the multiscale model is able to reproduce the tissue response at the effective scale (of the order of millimeters) while modeling the vasculature at the microscale. We validate the multiscale method against a full scale (three-dimensional) model, where the fluid/tissue interface is fully discretized and treated as a Neumann boundary for the elasticity equation. Next, we present simulation results obtained with the proposed approach in a realistic scenario, demonstrating that the method can robustly and efficiently handle the one-way coupling between complex fluid microstructures and the elastic matrix.
- ItemComputational haemodynamics in stenotic internal jugular veins(Berlin : Weierstraß-Institut für Angewandte Analysis und Stochastik, 2013) Caiazzo, Alfonso; Montecinos, Gino; Müller, Lucas O.; Toro, Eleuterio F.; Haacke, E. MarkStenosis in internal jugular veins (IJVs) are frequently associated to pathological venous circulation and insufficient cerebral blood drainage. In this work, we set up a computational framework to assess the relevance of IJV stenoses through numerical simulation, combining medical imaging, patient-specific data and a mathematical model for venous occlusions. Coupling a three-dimensional (3D) description of blood flow in IJVs with a reduced one-dimesional model (1D) for major intracranial veins, we are able to model different anatomical configurations, an aspect of importance to understand the impact of IJV stenosis in intracranial venous haemodynamics. We investigate several stenotic configurations in a physiologic patient-specific regime, quantifying the effect of the stenosis in terms of venous pressure increase and wall shear stress patterns. Simulation results are in qualitative agreement with reported pressure anomalies in pathological cases. Moreover, they demonstrate the potential of the proposed multiscale framework for individualbased studies and computer-aided diagnosis.
- ItemAssessment of reduced order Kalman filter for parameter identification in one-dimensional blood flow models using experimental data(Berlin : Weierstraß-Institut für Angewandte Analysis und Stochastik, 2016) Caiazzo, Alfonso; Caforio, Federica; Montecinos, Gino; Müller, Lucas O.; Blanco, Pablo J.; Toro, Eleutero F.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.