Computational Modeling and Analysis of Murmurs Generated by Modeled Aortic Stenoses

On the use of steady flow rates for approximating flow instabilities and vibrations in intracranial aneurysms

Recent computational and experimental studies of intracranial aneurysms have revealed potential mechanisms of aneurysm bruits and murmurs, driven by flow instabilities rather than by stable pulsatile flow. Some of these studies have been conducted under the assumption of constant flow rate (steady flow); however the validity of this assumption has not been evaluated for high-frequency flow instability, or vibrations from fluid-structure interaction (FSI) simulations. We evaluated the time-averaged wall shear stress, flow instability and vibration amplitude of steady flow simulations, performed at both cycle-averaged and peak-systolic flow rates, and compared these to recent pulsatile FSI simulations. Wall shear stress fields of pulsatile flow (time-averaged and peak values) were well-approximated by the respective steady-flow FSI simulations, and the spatial distribution and frequency content of flow instability and vibrations were reasonably approximated by the steady flow simulations at peak-systolic flow rates. However, the level of flow instability and vibration was generally over-predicted by the steady flow simulations at peak-systolic flow rates as flow remained unstable for longer than in the pulsatile simulation, while no flow instability was detected for steady flow at cycle-averaged flow rates. Additionally, the amplitude of flow instability and vibration fluctuated considerably in the steady flow simulations, while the pulsatile simulations exhibited consistent vibration amplitudes (less than 10 % variation at peak systole between cycles). Finally, steady flow simulations at peak-systolic conditions required 2-3x more compute time than the pulsatile simulations for the same time duration. Therefore, we recommend using pulsatile flow simulations when investigating vibrations and flow instabilities.

A 3D scaling law for supravalvular aortic stenosis suited for stethoscopic auscultations

In this study a frequency scaling law for 3D anatomically representative supravalvular aortic stenosis (SVAS) cases is proposed. The law is uncovered for stethoscopy's preferred auscultation range (70-120 Hz). LES simulations are performed on the CFD solver Fluent, leveraging Simulia's Living Heart Human Model (LHHM), modified to feature hourglass stenoses that range between 30 to 80 percent (mild to severe) in addition to the descending aorta. For physiological hemodynamic boundary conditions the Windkessel model is implemented via a UDF subroutine. The flow-generated acoustic signal is then extracted using the FW-H model and analyzed using FFT. A preferred receiver location that matches clinical practice is confirmed (right intercostal space) and a correlation between the degree of stenosis and a corresponding acoustic frequency is obtained. Five clinical auscultation signals are tested against the scaling law, with the findings interpreted in relation to the NHS classification of stenosis and to the assessments of experienced cardiologists. The scaling law is thus shown to succeed as a potential quantitative decision-support tool for clinicians, enabling them to reliably interpret stethoscopic auscultations for all degrees of stenosis, which is especially useful for moderate degrees of SVAS. Computational investigation of more complex stenotic cases would enhance the clinical relevance of this proposed scaling law, and will be explored in future research.

A CFD-FFT approach to hemoacoustics that enables degree of stenosis prediction from stethoscopic signals

In this paper, we identify a new (acoustic) frequency-stenosis relation whose frequencies lie within the recommended auscultation threshold of stethoscopy (< 120 Hz). We show that this relation can be used to extend the application of phonoangiography (quantifying the degree of stenosis from bruits) to widely accessible stethoscopes. The relation is successfully identified from an analysis restricted to the acoustic signature of the von Karman vortex street, which we automatically single out by means of a metric we propose that is based on an area-weighted average of the Q-criterion for the post-stenotic region. Specifically, we perform CFD simulations on internal flow geometries that represent stenotic blood vessels of different severities. We then extract their emitted acoustic signals using the Ffowcs Williams-Hawkings equation, which we subtract from a clean signal (stenosis free) at the same heart rate. Next, we transform this differential signal to the frequency domain and carefully classify its acoustic signatures per six (stenosis-)invariant flow phases of a cardiac cycle that are newly identified in this paper. We then automatically restrict our acoustic analysis to the sounds emitted by the von Karman vortex street (phase 4) by means of our Q-criterion-based metric. Our analysis of its acoustic signature reveals a strong linear relationship between the degree of stenosis and its dominant frequency, which differs considerably from the break frequency and the heart rate (known dominant frequencies in the literature). Applying our new relation to available stethoscopic data, we find that its predictions are consistent with clinical assessment. Our finding of this linear correlation is also unlike prevalent scaling laws in the literature, which feature a small exponent (i.e., low stenosis percentage sensitivity over much of the clinical range). They hence can only distinguish mild, moderate, and severe cases. Conversely, our linear law can identify variations in the degree of stenosis sensitively and accurately for the full clinical range, thus significantly improving the utility of the relevant scaling laws... Future research will investigate incorporating the vibroacoustic role of adjacent organs to expand the clinical applicability of our findings. Extending our approach to more complex 3D stenotic morphologies and including the vibroacoustic role of surrounding organs will be explored in future research to advance the clinical reach of our findings.

Computational Methods for Fluid-Structure Interaction Simulation of Heart Valves in Patient-Specific Left Heart Anatomies

Given the complexity of human left heart anatomy and valvular structures, the fluid–structure interaction (FSI) simulation of native and prosthetic valves poses a significant challenge for numerical methods. In this review, recent numerical advancements for both fluid and structural solvers for heart valves in patient-specific left hearts are systematically considered, emphasizing the numerical treatments of blood flow and valve surfaces, which are the most critical aspects for accurate simulations. Numerical methods for hemodynamics are considered under both the continuum and discrete (particle) approaches. The numerical treatments for the structural dynamics of aortic/mitral valves and FSI coupling methods between the solid Ωs and fluid domain Ωf are also reviewed. Future work toward more advanced patient-specific simulations is also discussed, including the fusion of high-fidelity simulation within vivo measurements and physics-based digital twining based on data analytics and machine learning techniques.

Computational Modeling of Aortic Stenosis With a Reduced Degree-of-Freedom Fluid-Structure Interaction Valve Model

In this study, a novel reduced degree-of-freedom (rDOF) aortic valve model is employed to investigate the fluid-structure interaction (FSI) and hemodynamics associated with aortic stenosis. The dynamics of the valve leaflets are determined by an ordinary differential equation with two parameters and this rDOF model is shown to reproduce key features of more complex valve models. The hemodynamics associated with aortic stenosis is studied for three cases: a healthy case and two stenosed cases. The focus of the study is to correlate the hemodynamic features with the source generation mechanism of systolic murmurs associated with aortic stenosis. In the healthy case, extremely weak flow fluctuations are observed. However, in the stenosed cases, simulations show significant turbulent fluctuations in the ascending aorta, which are responsible for the generation of strong wall pressure fluctuations after the aortic root mostly during the deceleration phase of the systole. The intensity of the murmur generation increases with the severity of the stenosis, and the source locations for the two diseased cases studied here lie around 1.0 inlet duct diameters (Do) downstream of the ascending aorta.