A reconstruction method for gappy and noisy arterial flow data

Efficient Cardiovascular Parameters Estimation for Fluid-Structure Simulations Using Gappy Proper Orthogonal Decomposition

As full-scale detailed hemodynamic simulations of the entire vasculature are not feasible, numerical analysis should be focused on specific regions of the cardiovascular system, which requires the identification of lumped parameters to represent the patient behavior outside the simulated computational domain. We present a novel technique for estimating cardiovascular model parameters using gappy Proper Orthogonal Decomposition (g-POD). A POD basis is constructed with FSI simulations for different values of the lumped model parameters, and a linear operator is applied to retain information that can be compared to the available patient measurements. Then, the POD coefficients of the reconstructed solution are computed either by projecting patient measurements or by solving a minimization problem with constraints. The POD reconstruction is then used to estimate the model parameters. In the first test case, the parameter values of a 3-element Windkessel model are approximated using artificial patient measurements, obtaining a relative error of less than 4.2%. In the second case, 4 sets of 3-element Windkessel are approximated in a patient's aorta geometry, resulting in an error of less than 8% for the flow and less than 5% for the pressure. The method shows accurate results even with noisy patient data. It automatically calculates the delay between measurements and simulations and has flexibility in the types of patient measurements that can handle (at specific points, spatial or time averaged). The method is easy to implement and can be used in simulations performed in general-purpose FSI software.

Proper Orthogonal Decomposition Analysis of the Flow Downstream of a Dysfunctional Bileaflet Mechanical Aortic Valve

PURPOSE

Aortic valve replacement remains the only viable solution for symptomatic patients with severe aortic valve stenosis. Despite their improved design and long history of successful operation, bileaflet mechanical heart valves are still associated with post-operative complications leading to valve dysfunction. Thus, the flow dynamics can be highly disturbed downstream of the dysfunctional valve.

METHODS

In this in vitro study, the flow dynamics downstream of healthy and dysfunctional bileaflet mechanical heart valves have been investigated using particle image velocimetry measurements. Proper orthogonal decomposition of the velocity field has been performed in order to explore the coherent flow features in the ascending aorta in the presence of a dysfunctional bileaflet mechanical heart valve.

RESULTS

The ability of proper orthogonal decomposition derived metrics to differentiate between heathy and dysfunctional cases is reported. Moreover, reduced-order modeling using proper orthogonal decomposition is thoroughly investigated not only for the velocity field but also for higher order flow characteristics such as time average wall shear stress, oscillatory shear index and viscous energy dissipation.

CONCLUSION

Considering these results, proper orthogonal decomposition can provide a rapid binary classifier to evaluate if the bileaflet mechanical valve deviates from its normal operating conditions. Moreover, the study shows that the size of the reduced-order model depends on which flow parameter is required to be reconstructed.

Dynamic Denoising and Gappy Data Reconstruction Based on Dynamic Mode Decomposition and Discrete Cosine Transform

Dynamic Mode Decomposition (DMD) is a data-driven method to analyze the dynamics, first applied to fluid dynamics. It extracts modes and their corresponding eigenvalues, where the modes are spatial fields that identify coherent structures in the flow and the eigenvalues describe the temporal growth/decay rates and oscillation frequencies for each mode. The recently introduced compressed sensing DMD (csDMD) reduces computation times and also has the ability to deal with sub-sampled datasets. In this paper, we present a similar technique based on discrete cosine transform to reconstruct the fully-sampled dataset (as opposed to DMD modes as in csDMD) from sub-sampled noisy and gappy data using l 1 minimization. The proposed method was benchmarked against csDMD in terms of denoising and gap-filling using three datasets. The first was the 2-D time-resolved plot of a double gyre oscillator which has about nine oscillatory modes. The second dataset was derived from a Duffing oscillator. This dataset has several modes associated with complex eigenvalues which makes them oscillatory. The third dataset was taken from the 2-D simulation of a wake behind a cylinder at Re = 100 and was used for investigating the effect of changing various parameters on reconstruction error. The Duffing and 2-D wake datasets were tested in presence of noise and rectangular gaps. While the performance for the double-gyre dataset is comparable to csDMD, the proposed method performs substantially better (lower reconstruction error) for the dataset derived from the Duffing equation and also, the 2-D wake dataset according to the defined reconstruction error metrics.

Analyzing transient turbulence in a stenosed carotid artery by proper orthogonal decomposition

High-resolution three-dimensional simulations (involving 100 million degrees of freedom) were employed to study transient turbulent flow in a carotid arterial bifurcation with a stenosed internal carotid artery (ICA). The geometrical model was reconstructed from MRI images, and in vivo velocity measurements were incorporated in the simulations to provide inlet and outlet boundary conditions. Due to the high degree of the ICA occlusion and the variable flow rate, a transitional and intermittent flow between laminar and turbulent states was established. Time- and space-window proper orthogonal decomposition (POD) was applied to quantify the different flow regimes in the occluded artery. A simplified version of the POD analysis that utilizes 2D slices only--more appropriate in the clinical setting--was also investigated.