Longitudinal wall shear stress evaluation using centerline projection approach in the numerical simulations of the patient-based carotid artery

In this numerical study, areas of the carotid bifurcation and of a distal stenosis in the internal carotid artery are closely observed to evaluate the patient's current risks of ischemic stroke. An indicator for the vessel wall defects is the stress exerted by blood on the vessel tissue, typically expressed by the amplitude of the wall shear stress vector (WSS) and its oscillatory shear index. To detect negative shear stresses corresponding with reversal flow, we perform orientation-based shear evaluation. We investigate the longitudinal component of the wall shear vector, where tangential vectors aligned longitudinally with the vessel are necessary. However, resulting from imaging segmentation resolution of patients' computed tomography angiography scans and stenotic regions, the geometry model's mesh is non-smooth on its surface areas and the automatically generated tangential vector field is discontinuous and multi-directional, making an interpretation of our orientation-based risk indicators unreliable. We improve the evaluation of longitudinal shear stress by applying the projection of the vessel's centerline to the surface to construct smooth tangential field aligned longitudinally with the vessel. We validate our approach for the longitudinal WSS component and the corresponding oscillatory index by comparing them to results obtained using automatically generated tangents in both rigid and elastic vessel modeling and to amplitude-based indicators. We present the major benefit of our longitudinal WSS evaluation based on its directionality for the cardiovascular risk assessment, which is the detection of negative WSS indicating persistent reversal or transverse flow. This is impossible in the case of the amplitude-based WSS.

Integrated particle image velocimetry and fluid-structure interaction analysis for patient-specific abdominal aortic aneurysm studies

BACKGROUND

Understanding the hemodynamics of an abdominal aortic aneurysm (AAA) is crucial for risk assessment and treatment planning. This study introduces a low-cost, patient-specific in vitro AAA model to investigate hemodynamics using particle image velocimetry (PIV) and flow-simulating circuit, validated through fluid-structure interaction (FSI) simulations.

METHODS

In this study, 3D printing was employed to manufacture a flexible patient-specific AAA phantom using a lost-core casting technique. A pulsatile flow circuit was constructed using off-the-shelf components. A particle image velocimetry (PIV) setup was built using an affordable laser source and global shutter camera, and finally, the flow field inside the AAA was analyzed using open-source software. Fluid-structure interaction (FSI) simulations were performed to enhance our understanding of the flow field, and the results were validated by PIV analysis. Both steady-state and transient flow conditions were investigated.

RESULTS

Our experimental setup replicated physiological conditions, analyzing arterial wall deformations and flow characteristics within the aneurysm. Under constant flow, peak wall deformations and flow velocities showed deviations within - 12% to + 27% and - 7% to + 5%, respectively, compared to FSI simulations. Pulsatile flow conditions further demonstrated a strong correlation (Pearson coefficient 0.85) in flow velocities and vectors throughout the cardiac cycle. Transient phenomena, particularly the formation and progression of vortex structures during systole, were consistently depicted between experimental and numerical models.

CONCLUSIONS

By bridging high-fidelity experimental observations with comprehensive computational analyses, this study underscores the potential of integrated methodologies in enhancing our understanding of AAA pathophysiology. The convergence of realistic AAA phantoms, precise PIV measurements at affordable cost point, and validated FSI models heralds a new paradigm in vascular research, with significant implications for personalized medicine and bioengineering innovations.

Structural numerical analysis of a branched modular stent-graft for aneurysms encompassing all zones of the aortic arch

The development of stent-grafts for the total repair of aneurysms in the aortic arch is still a technical challenge due mainly to the anatomical complexity of this region. Research performed here structurally evaluated a modular branched stent-graft for aneurysms encompassing all zones of the aortic arch by means of numerical simulations using fluid-structure interaction. The geometric domain obtained by means of computed tomography was subjected to physiological boundary conditions. The blood was modelled as non-Newtonian by the Carreau model, and the arterial wall was modelled as anisotropic hyperelastic by the Holzapfel model. The material adopted for the stents was Nitinol, and expanded polytetrafluoroethylene (ePTFE) was used for the graft. A comparison of the structural behaviour of the aneurysmal aortic arch before and after stent-graft implantation was performed. The numerical flow model was experimentally verified in vitro on a representative test bench of blood flow in the aortic arch. The stent-graft was shown to minimally modify arterial wall dynamics and was not susceptible to migration and endoleak. Peak stresses and strains were found in the stents and graft, respectively, while the stresses in the aneurysm sac were significantly reduced, of the order of 97.5%, due to the isolation of the arterial wall by the stent-graft.

The impact of endotracheal intubation on oxygen delivery, trachea pressure and wall deformation

This paper concerns improving endotracheal tube (ETT) insertion through advanced computational science modelling. The study aims to better understand endotracheal intubation (ETI) and reduce medical errors in intensive and critical care units since ETT insertion is unique for each patient, depending on age, gender, size, physiology, and underlying health conditions. We have employed computational fluid dynamics and biomechanics modelling to investigate the effect of ETT for three ventilation modes on (a) local oxygen delivery to the lungs, (b) air pressure and wall shear stress at the tracheal walls, and (c) oscillatory elastic deformation of the tracheal tissues and muscle. For the first time, we reveal how the ventilation mode and ETT insertion in the trachea may induce major complications, especially in long periods of ETT. We show that rotating the ETT or displacing it by 2 mm only can induce a significant rise in the tracheal pressure up to 177 cmH2O. This study, for the first time, shows the vital role of computers in biology and medicine to provide enhanced decision-making-support to clinicians and medical doctors dealing with ETI.

Numerical Study of a Thrombus Migration Risk in Aneurysm After Coil Embolization in Patient Cases: FSI Modelling

PURPOSE

There are still many challenges for modelling a thrombus migration process in aneurysms. The main novelty of the present research lies in the modelling of aneurysm clot migration process in a realistic cerebral aneurysm, and the analysis of forces suffered by clots inside an aneurysm, through transient FSI simulations.

METHODS

The blood flow has been modelled using a Womersley velocity profile, and following the Carreau viscosity model. Hyperelastic Ogden model has been used for clot and isotropic linear elastic model for the artery walls. The FSI coupled model was implemented in ANSYS® software. The hemodynamic forces suffered by the clot have been quantified using eight different clot sizes and positions inside a real aneurysm.

RESULTS

The obtained results have shown that it is almost impossible for clots adjacent to aneurysm walls, to leave the aneurysm. Nevertheless, in clots positioned in the centre of the aneurysm, there is a real risk of clot migration. The risk of migration of a typical post-coiling intervention clot in an aneurysm, in contact with the wall and occupying a significant percentage of its volume is very low in the case studied, even in the presence of abnormally intense events, associated with sneezes or impacts.

CONCLUSIONS

The proposed methodology allows evaluating the clot migration risk, vital for evaluating the progress after endovascular interventions, it is a step forward in the personalized medicine, patient follow-up, and helping the medical team deciding the optimal treatment.

Numerical investigation of shock-induced bubble collapse dynamics and fluid-solid interactions during shock-wave lithotripsy

In this paper we investigate the bubble collapse dynamics under shock-induced loading near soft and rigid bio-materials, during shock wave lithotripsy. A novel numerical framework was developed, that employs a Diffuse Interface Method (DIM) accounting for the interaction across fluid-solid-gas interfaces. For the resolution of the extended variety of length scales, due to the dynamic and fine interfacial structures, an Adaptive Mesh Refinement (AMR) framework for unstructured grids was incorporated. This multi-material multi-scale approach aims to reduce the numerical diffusion and preserve sharp interfaces. The presented numerical framework is validated for cases of bubble dynamics, under high and low ambient pressure ratios, shock-induced collapses, and wave transmission problems across a fluid-solid interface, against theoretical and numerical results. Three different configurations of shock-induced collapse applications near a kidney stone and soft tissue have been simulated for different stand-off distances and bubble attachment configurations. The obtained results reveal the detailed collapse dynamics, jet formation, solid deformation, rebound, primary and secondary shock wave emissions, and secondary collapse that govern the near-solid collapse and penetration mechanisms. Significant correlations of the problem configuration to the overall collapse mechanisms were found, stemming from the contact angle/attachment of the bubble and from the properties of solid material. In general, bubbles with their center closer to the kidney stone surface produce more violent collapses. For the soft tissue, the bubble movement prior to the collapse is of great importance as new structures can emerge which can trap the liquid jet into induced crevices. Finally, the tissue penetration is examined for these cases and a novel tension-driven tissue injury mechanism is elucidated, emanating from the complex interaction of the bubble/tissue interaction during the secondary collapse phase of an entrapped bubble in an induced crevice with the liquid jet.

Onset and nature of flow-induced vibrations in cerebral aneurysms via fluid-structure interaction simulations

Clinical, experimental, and recent computational studies have demonstrated the presence of wall vibrations in cerebral aneurysms, thought to be induced by blood flow instability. These vibrations could induce irregular, high-rate deformation of the aneurysm wall, and potentially disrupt regular cell behavior and promote deleterious wall remodeling. In order to elucidate, for the first time, the onset and nature of such flow-induced vibrations, in this study we imposed a linearly increasing flow rate on high-fidelity fluid-structure interaction models of three anatomically realistic aneurysm geometries. Prominent narrow-band vibrations in the range of 100-500 Hz were found in two out of the three aneurysm geometries tested, while the case that did not exhibit flow instability did not vibrate. Aneurysm vibrations consisted mostly of fundamental modes of the entire aneurysm sac, with the vibrations exhibiting more frequency content at higher frequencies than the flow instabilities driving those vibrations. The largest vibrations occurred in the case which exhibited strongly banded fluid frequency content, and the vibration amplitude was highest when the strongest fluid frequency band was an integer multiple of one of the natural frequencies of the aneurysm sac. Lower levels of vibration occurred in the case which exhibited turbulent-like flow with no distinct frequency bands. The current study provides a plausible mechanistic explanation for the high-frequency sounds observed in cerebral aneurysms, and suggests that narrow-band (vortex-shedding type) flow might stimulate the wall more, or at least at lower flow rates, than broad-band, turbulent-like flow.

Peristaltic regimes in esophageal transport

A FLIP device gives cross-sectional area along the length of the esophagus and one pressure measurement, both as a function of time. Deducing mechanical properties of the esophagus including wall material properties, contraction strength, and wall relaxation from these data are a challenging inverse problem. Knowing mechanical properties can change how clinical decisions are made because of its potential for in-vivo mechanistic insights. To obtain such information, we conducted a parametric study to identify peristaltic regimes by using a 1D model of peristaltic flow through an elastic tube closed on both ends and also applied it to interpret clinical data. The results gave insightful information about the effect of tube stiffness, fluid/bolus density and contraction strength on the resulting esophagus shape through quantitive representations of the peristaltic regimes. Our analysis also revealed the mechanics of the opening of the contraction area as a function of bolus flow resistance. Lastly, we concluded that peristaltic driven flow displays three modes of peristaltic geometries, but all physiologically relevant flows fall into two peristaltic regimes characterized by a tight contraction.

Validation of the Reduced Unified Continuum Formulation Against In Vitro 4D-Flow MRI

We previously introduced and verified the reduced unified continuum formulation for vascular fluid-structure interaction (FSI) against Womersley's deformable wall theory. Our present work seeks to investigate its performance in a patient-specific aortic setting in which assumptions of idealized geometries and velocity profiles are invalid. Specifically, we leveraged 2D magnetic resonance imaging (MRI) and 4D-flow MRI to extract high-resolution anatomical and hemodynamic information from an in vitro flow circuit embedding a compliant 3D-printed aortic phantom. To accurately reflect experimental conditions, we numerically implemented viscoelastic external tissue support, vascular tissue prestressing, and skew boundary conditions enabling in-plane vascular motion at each inlet and outlet. Validation of our formulation is achieved through close quantitative agreement in pressures, lumen area changes, pulse wave velocity, and early systolic velocities, as well as qualitative agreement in late systolic flow structures. Our validated suite of FSI techniques offers a computationally efficient approach for numerical simulation of vascular hemodynamics. This study is among the first to validate a cardiovascular FSI formulation against an in vitro flow circuit involving a compliant vascular phantom of complex patient-specific anatomy.

The Impact of a Limited Field-of-View on Computed Hemodynamics in Abdominal Aortic Aneurysms: Evaluating the Feasibility of Completing Ultrasound Segmentations with Parametric Geometries

To improve abdominal aortic aneurysm (AAA) rupture risk assessment, a large, longitudinal study on AAA hemodynamics and biomechanics is necessary, using personalized fluid-structure interaction (FSI) modeling. 3-dimensional, time-resolved ultrasound (3D+t US) is the preferred image modality to obtain the patient-specific AAA geometry for such a study, since it is safe, affordable and provides temporal information. However, the 3D+t US field-of-view (FOV) is limited and therefore often fails to capture the inlet and aorto-iliac bifurcation geometry. In this study, a framework was developed to add parametric inlet and bifurcation geometries to the abdominal aortic aneurysm geometry by employing dataset statistics and parameters of the AAA geometry. The impact of replacing the patient-specific inlet and bifurcation geometries, acquired using computed tomography (CT) scans, by parametric geometries was evaluated by examining the differences in hemodynamics (systolic and time-averaged wall shear stress and oscillatory shear index) in the aneurysm region. The results show that the inlet geometry has a larger effect on the AAA hemodynamics (median differences of 7.5 to 18.8%) than the bifurcation geometry (median differences all below 1%). Therefore, it is not feasible to replace the patient-specific inlet geometry by a generic one. Future studies should investigate the possibilities of extending the proximal FOV of 3D+t US. However, this study did show the feasibility of adding a parametric bifurcation geometry to the aneurysm geometry. After extending the proximal FOV, the obtained framework can be used to extract AAA geometries from 3D+t US for FSI simulations, despite the absence of the bifurcation geometry.

Impact of cyclic bending on coronary hemodynamics

It remains unknown that the degree of bias in computational fluid dynamics results without considering coronary cyclic bending. This study aims to investigate the influence of different rates of coronary cyclic bending on coronary hemodynamics. To model coronary bending, a multi-ring-controlled fluid-structural interaction model was designed. A coronary artery was simulated with various cyclic bending rates (0.5, 0.75 and 1 s, corresponding to heart rates of 120, 80 and 60 bpm) and compared against a stable model. The simulated results show that the hemodynamic parameters of vortex Q-criterion, temporal wall shear stress (WSS), time-averaged WSS (TaWSS) and oscillatory shear index (OSI) were sensitive to the changes in cyclic rate. A higher heart rate resulted in higher magnitude and larger variance in the hemodynamic parameters. Whereas, the values and distributions of flow velocity and relative residence time (RRT) did not show significant differences between different bending periods. This study suggests that a stable coronary model is not sufficient to represent the hemodynamics in a bending coronary artery. Different heart rate conditions were found to have significant impact on the hemodynamic parameters. Thus, cyclic bending should be considered to mimic the realistic hemodynamics in future patient-specific coronary hemodynamics studies.

Flow-based method demonstrates improved accuracy for calculating wall shear stress in arterial flows from 4D flow MRI data

Four-dimensional flow magnetic resonance imaging (i.e., 4D flow MRI) has become a valuable tool for the in vivo assessment of blood flow within large vessels and cardiac chambers. As wall shear stress (WSS) has been correlated with the development and progression of cardiovascular disease, focus has been directed at developing techniques to quantify WSS directly from 4D flow MRI data. The goal of this study was to compare the accuracy of two such techniques - termed the velocity and flow-based methods - in the setting of simplified and complex flow scenarios. Synthetic MR data were created from exact solutions to the Navier-Stokes equations for the steady and pulsatile flow of an incompressible, Newtonian fluid through a rigid cylinder. In addition, synthetic MR data were created from the predicted velocity fields derived from a fluid-structure interaction (FSI) model of pulsatile flow through a thick-walled, multi-layered model of the carotid bifurcation. Compared to the analytical solutions for steady and pulsatile flow, the flow-based method demonstrated greater accuracy than the velocity-based method in calculating WSS across all changes in fluid velocity/flow rate, tube radius, and image signal-to-noise (p < 0.001). Furthermore, the velocity-based method was more sensitive to boundary segmentation than the flow-based method. When compared to results from the FSI model, the flow-based method demonstrated greater accuracy than the velocity-based method with average differences in time-averaged WSS of 0.31 ± 1.03 Pa and 0.45 ± 1.03 Pa, respectively (p <0.005). These results have implications on the utility, accuracy, and clinical translational of methods to determine WSS from 4D flow MRI.

Semi-Automated Construction of Patient-Specific Aortic Valves from Computed Tomography Images

This paper presents a semi-automatic method for the construction of volumetric models of the aortic valve using computed tomography angiography images. Although the aortic valve typically cannot be segmented directly from a computed tomography angiography image, the method described herein uses manually selected samples of an aortic segmentation derived from this image to inform the construction. These samples capture certain physiologic landmarks and are used to construct a volumetric valve model. As a demonstration of the capabilities of this method, valve models for 25 pediatric patients are created. A selected valve anatomy is used to perform fluid-structure interaction simulations using the immersed finite element/difference method with physiologic driving and loading conditions. Simulation results demonstrate this method creates a functional valve that opens and closes normally and generates pressure and flow waveforms that are similar to those observed clinically.

The pathogenesis of superior mesenteric artery dissection: An in-depth study based on fluid-structure interaction and histology analysis

OBJECTIVE

To explore the role of hemodynamic factors in the occurrence of superior mesenteric artery (SMA) dissection (SMAD) using a fluid-structure interaction (FSI) simulation method, and to identify histopathologic changes occurring in the wall of the SMA.

METHODS

A total of 122 consecutive patients diagnosed with SMAD and 122 controls were included in this study. Hemodynamic factors were calculated using a FSI simulation method. Additionally, SMA specimens obtained from 12 cadavers were stained for histological quantitative analysis.

RESULTS

The mean aortomesenteric angle (59.7° ± 21.4° vs 48.2° ± 16.8°; p < .001) and SMA maximum curvature (0.084 ± 0.078 mm^-1 vs 0.032 ± 0.023 mm^-1; p < .001) were higher in SMAD patients than the controls. Larger aortomesenteric angles and SMA curvatures were associated with higher and more concentrated wall shear stress at anterior wall of the SMA curve segment, co-located with the dissection origins. The mean thickness of media (325.18 ± 44.87 µm vs 556.92 ± 138.32 µm; p = .003) was thinner in the anterior wall of the SMA curve than in the posterior wall. The area fractions of elastin (17.96% ± 3.36% vs 27.06% ± 4.18%; p = .002) and collagen (45.43% ± 6.89% vs 55.57% ± 7.57%; p = .036) were lower in anterior wall of the SMA curve than in posterior wall.

CONCLUSION

Increased aortomesenteric angle and SMA curvature are risk factors for SMAD. Both of these factors can cause local hemodynamic abnormalities, which can lead to histopathologic changes in anterior wall of SMA.

Experimental and Mouse-Specific Computational Models of the Fbln4SMKO Mouse to Identify Potential Biomarkers for Ascending Thoracic Aortic Aneurysm

PURPOSE

To use computational methods to explore geometric, mechanical, and fluidic biomarkers that could correlate with mouse lifespan in the Fbln4SMKO mouse. Mouse lifespan was used as a surrogate for risk of a severe cardiovascular event in cases of ascending thoracic aortic aneurysm.

METHODS

Image-based, mouse-specific fluid-structure-interaction models were developed for Fbln4SMKO mice (n = 10) at ages two and six months. The results of the simulations were used to quantify potential biofluidic biomarkers, complementing the geometrical biomarkers obtained directly from the images.

RESULTS

Comparing the different geometrical and biofluidic biomarkers to the mouse lifespan, it was found that mean oscillatory shear index (OSImin) and minimum time-averaged wall shear stress (TAWSSmin) at six months showed the largest correlation with lifespan (r2 = 0.70, 0.56), with both correlations being positive (i.e., mice with high OSImean and high TAWSSmin tended to live longer). When change between two and six months was considered, the change in TAWSSmin showed a much stronger correlation than OSImean (r2 = 0.75 vs. 0.24), and the correlation was negative (i.e., mice with increasing TAWSSmin over this period tended to live less long).

CONCLUSION

The results highlight potential biomarkers of ATAA outcomes that can be obtained through noninvasive imaging and computational simulations, and they illustrate the potential synergy between small-animal and computational models.

Planktos: An Agent-Based Modeling Framework for Small Organism Movement and Dispersal in a Fluid Environment with Immersed Structures

Multiscale modeling of marine and aerial plankton has traditionally been difficult to address holistically due to the challenge of resolving individual locomotion dynamics while being carried with larger-scale flows. However, such problems are of paramount importance, e.g., dispersal of marine larval plankton is critical for the health of coral reefs, and aerial plankton (tiny arthropods) can be used as effective agricultural biocontrol agents. Here we introduce the open-source, agent-based modeling software Planktos targeted at 2D and 3D fluid environments in Python. Agents in this modeling framework are relatively tiny organisms in sufficiently low densities that their effect on the surrounding fluid motion can be considered negligible. This library can be used for scientific exploration and quantification of collective and emergent behavior, including interaction with immersed structures. In this paper, we detail the implementation and functionality of the library along with some illustrative examples. Functionality includes arbitrary agent behavior obeying either ordinary differential equations, stochastic differential equations, or coded movement algorithms, all under the influence of time-dependent fluid velocity fields generated by computational fluid dynamics, experiments, or analytical models in domains with static immersed mesh structures with sliding or sticky collisions. In addition, data visualization tools provide images or animations with kernel density estimation and velocity field analysis with respect to deterministic agent behavior via the finite-time Lyapunov exponent.

Cavitation bubble interaction with compliant structures on a microscale: A contribution to the understanding of bacterial cell lysis by cavitation treatment

Numerous studies have already shown that the process of cavitation can be successfully used for water treatment and eradication of bacteria. However, most of the relevant studies are being conducted on a macro scale, so the understanding of the processes at a fundamental level remains poor. In attempt to further elucidate the process of cavitation-assisted water treatment on a scale of a single bubble, the present paper numerically addresses interaction between a collapsing microbubble and a nearby compliant structure, that mechanically and structurally resembles a bacterial cell. A fluid-structure interaction methodology is employed, where compressible multiphase flow is considered and the bacterial cell wall is modeled as a multi-layered shell structure. Simulations are performed for two selected model structures, each resembling the main structural features of Gram-negative and Gram-positive bacterial cell envelopes. The contribution of two independent dimensionless geometric parameters is investigated, namely the bubble-cell distance δ and their size ratio ς. Three characteristic modes of bubble collapse dynamics and four modes of spatiotemporal occurrence of peak local stresses in the bacterial cell membrane are identified throughout the parameter space considered. The former range from the development of a weak and thin jet away from the cell to spherical bubble collapses. The results show that local stresses arising from bubble-induced loads can exceed poration thresholds of cell membranes and that bacterial cell damage could be explained solely by mechanical effects in absence of thermal and chemical ones. Based on this, the damage potential of a single microbubble for bacteria eradication is estimated, showing a higher resistance of the Gram-positive model organism to the nearby bubble collapse. Microstreaming is identified as the primary mechanical mechanism of bacterial cell damage, which in certain cases may be enhanced by the occurrence of shock waves during bubble collapse. The results are also discussed in the scope of bacteria eradication by cavitation treatment on a macro scale, where processes of hydrodynamic and ultrasonic cavitation are being employed.

Effects of membrane and flexural stiffnesses on aortic valve dynamics: identifying the mechanics of leaflet flutter in thinner biological tissues

Valvular pathologies that induce deterioration in the aortic valve are a common cause of heart disease among aging populations. Although there are numerous available technologies to treat valvular conditions and replicate normal aortic function by replacing the diseased valve with a bioprosthetic implant, many of these devices face challenges in terms of long-term durability. One such phenomenon that may exacerbate valve deterioration and induce undesirable hemodynamic effects in the aorta is leaflet flutter, which is characterized by oscillatory motion in the biological tissues. While this behavior has been observed for thinner bioprosthetic valves, the specific underlying mechanics that lead to leaflet flutter have not previously been identified. This work proposes a computational approach to isolate the fundamental mechanics that induce leaflet flutter in thinner biological tissues during the cardiac cycle. The simulations in this work identify reduced flexural stiffness as the primary factor that contributes to increased leaflet flutter in thinner biological tissues, while decreased membrane stiffness and mass of the thinner tissues do not directly induce flutter in these valves. The results of this study provide an improved understanding of the mechanical tissue properties that contribute to flutter and offer significant insights into possible developments in the design of bioprosthetic tissues to account for and reduce the incidence of flutter.

Pressure Difference Estimation in Non-stenotic Carotid Bifurcation Phantoms Using Vector Flow Imaging

Local pressure differences estimated using vector flow imaging (VFI) and direct catheterization in seven carotid bifurcation phantoms were compared with simulated pressure fields. VFI correlated strongly with simulated peak pressure differences (r = 0.99, p < 0.00001), with an average overestimation of 12.3 Pa (28.6%). The range between the lowest and highest pressure difference of VFI underestimated simulations by 4.6 Pa (8.06%; r = 0.99, p < 0.0001). The catheter method exhibited no correlation (r = -0.09, p = 0.85). Ten repeated measurements on one phantom revealed a small standard deviation (SD) for VFI (SD = 8.4%, mean estimated SD = 11.5%), but not for the catheter method (SD = 785.6%). An in vivo peak systolic pressure difference of 97.9 Pa (estimated SD = 30%) was measured using VFI in one healthy individual. This study indicates that VFI pressure difference estimation is feasible in phantoms and in vivo and realistic estimates of the SD can be attained from the data.

Targeted drug delivery of magnetic microbubble for abdominal aortic aneurysm: an in silico study

Targeted drug delivery (TDD) to abdominal aortic aneurysm (AAA) using a controlled and efficient approach has recently been a significant challenge. In this study, by using magnetic microbubbles (MMBs) under a magnetic field, we investigated the MMBs performance in TDD to AAA based on the amount of surface density of MMBs (SDMM) adhered to the AAA lumen. The results showed that among the types of MMBs studied in the presence of the magnetic field, micromarkers are the best type of microbubble with a -[Formula: see text] increase in SDMM adhered to the critical area of AAA. The results show that applying a magnetic field causes the amount of SDMM adhered to the whole area of AAA to increase -[Formula: see text] times compared to the condition in which the magnetic field is absent. This optimal and maximum value occurs for Definity MMBs with - 3.3 μm diameter. Applying a magnetic field also increases the adhesion surface density by - [Formula: see text], - [Formula: see text], and -[Formula: see text] times for the Micromarker, Optison, and Sonovue microbubbles, respectively, relative to the condition in which the magnetic field is absent. It was shown that using MBBs under magnetic field has the best performance in delivery to AAA for patients with negative inlet blood flow. Also, we have exposed that in an efficient TDD to AAA using MMBs, decreasing the density of MMBs increases drug delivery efficiency and performance. When density is - [Formula: see text], there is the highest difference (about - 75%) between the SDMM adhered to AAA in the presence of a magnetic field and in the absence of a magnetic field.

Investigation the effect of geometry and position of polymeric heart valves on hemodynamic with fluid-structure interaction numerical method

Computational modeling and numerical simulation of heart valve dynamics incorporating both fluid dynamics and valve structural replications has been challenging. In this study, we developed a double-coupled fluid-structure interaction (FSI) model using arbitrary lagrangian eulerian(ALE) and steered adaptive mesh(SAM). So we were looking to simulate transcatheter aortic valve (TAV) hemodynamic performance throughout entire cardiac cycles [1]. To reach this object, semi-real geometry of aorta and aortic polymeric valves has been created. At model inlet, left ventricular pressure and at the model outlet the elastic porous tube have been considered. Nonlinear finite element way and Sparse solver was utilized to couple fluid and solid equation. Consequently, extensive and comparative simulation were performed to investigate the impact of valve elasticity and valve positions on hemodynamics and solid parameters. Effective Orifice Area(EOA) also has been calculated [1]. The simulation results indicated that the lower of the elastic modulus cause to increase the EOA. Furthermore, the result of valve position showed, whenever the valve is closer to sinuses, a greater EOA and lower stresses impose on the leaflet are achievable.

Exploring the sensitivity in jellyfish locomotion under variations in scale, frequency, and duty cycle

Jellyfish have been called one of the most energy-efficient animals in the world due to the ease in which they move through their fluid environment, by product of their bell kinematics coupled with their morphological, muscular, material properties. We investigated jellyfish locomotion by conducting in silico comparative studies and explored swimming performance across different fluid scales (i.e., Reynolds Number), bell contraction frequencies, and contraction phase kinematics (duty cycle) for a jellyfish with a fineness ratio of 1 (ratio of bell height to bell diameter). To study these relationships, an open source implementation of the immersed boundary method was used (IB2d) to solve the fully coupled fluid-structure interaction problem of a flexible jellyfish bell in a viscous fluid. Thorough 2D parameter subspace explorations illustrated optimal parameter combinations in which give rise to enhanced swimming performance. All performance metrics indicated a higher sensitivity to bell actuation frequency than fluid scale or duty cycle, via Sobol sensitivity analysis, on a higher performance parameter subspace. Moreover, Pareto-like fronts were identified in the overall performance space involving the cost of transport and forward swimming speed. Patterns emerged within these performance spaces when highlighting different parameter regions, which complemented the global sensitivity results. Lastly, an open source computational model for jellyfish locomotion is offered to the science community that can be used as a starting place for future numerical experimentation.

Effect of plaque compositions on fractional flow reserve in a fluid-structure interaction analysis

Coronary artery disease involves the reduction of blood flow to the myocardium due to atherosclerotic plaques. The findings of myocardial ischemia may indicate severe coronary stenosis, but many studies have demonstrated a mismatch between lumen stenosis and fractional flow reserve (FFR). Recently, some clinical studies have found that the composition of atherosclerotic plaques may be a potential missing link between stenosis and ischemia. To investigate the relationship between myocardial ischemia and plaque composition, we have developed and adopted a new fluid-structure interaction (FSI) patient-specific coronary plaque model, based on computed tomography angiography data, to assess the impact on FFR as a biomechanical indicator of ischemia. A total of 180 analyses have been performed in 3D-FSI coronary artery disease models based on plaque compositions, plaque location, and stenosis degree. Hemodynamic analysis of simulation results and comparisons with other methods has been conducted to validate our models. Our results have successfully verified that the different compositions of plaques have resulted in differences in the calculated FFR. The mean FFR values with lipid plaques are [Formula: see text] as compared to the mean FFR values in lesions with fibrous plaques [Formula: see text] and calcified plaques [Formula: see text]. Besides, FFR differences between the three different plaque compositions have been shown to increase as the diameter stenosis increased. Plaque composition affects vascular stiffness and vascular dilation ability, and thereby affects the stenosis degree, resulting in abnormal FFR leading to myocardial ischemia. This interrelationship can help to diagnose the cause of high-risk coronary artery disease, leading to myocardial ischemia.

Fully coupled dynamic simulations of bioprosthetic aortic valves based on an embedded strategy for fluid-structure interaction with contact

AIMS

This work aims at presenting a fully coupled approach for the numerical solution of contact problems between multiple elastic structures immersed in a fluid flow. The key features of the computational model are (i) a fully coupled fluid-structure interaction with contact, (ii) the use of a fibre-reinforced material for the leaflets, (iii) a stent, and (iv) a compliant aortic root.

METHODS AND RESULTS

The computational model takes inspiration from the immersed boundary techniques and allows the numerical simulation of the blood-tissue interaction of bioprosthetic heart valves (BHVs) as well as the contact among the leaflets. First, we present pure mechanical simulations, where blood is neglected, to assess the performance of different material properties and valve designs. Secondly, fully coupled fluid-structure interaction simulations are employed to analyse the combination of haemodynamic and mechanical characteristics. The isotropic leaflet tissue experiences high-stress values compared to the fibre-reinforced material model. Moreover, elongated leaflets show a stress concentration close to the base of the stent. We observe a fully developed flow at the systolic stage of the heartbeat. On the other hand, flow recirculation appears along the aortic wall during diastole.

CONCLUSION

The presented FSI approach can be used for analysing the mechanical and haemodynamic performance of a BHV. Our study suggests that stresses concentrate in the regions where leaflets are attached to the stent and in the portion of the aortic root where the BHV is placed. The results from this study may inspire new BHV designs that can provide a better stress distribution.

Valveless pumping behavior of the simulated embryonic heart tube as a function of contractile patterns and myocardial stiffness

During development, the heart begins pumping as a valveless multilayered tube capable of driving blood flow throughout the embryonic vasculature. The mechanical properties and how they interface with pumping function are not well-defined at this stage. Here, we evaluate pumping patterns using a fluid-structure interaction computational model, combined with experimental data and an energetic analysis to investigate myocardial mechanical properties. Through this work, we propose that a myocardium modeled as a Neo-Hookean material with a material constant on the order of 10 kPa is necessary for the heart tube to function with an optimal pressure and cardiac output.

Fluid-structure interaction in a fully coupled three-dimensional mitral-atrium-pulmonary model

This paper aims to investigate detailed mechanical interactions between the pulmonary haemodynamics and left heart function in pathophysiological situations (e.g. atrial fibrillation and acute mitral regurgitation). This is achieved by developing a complex computational framework for a coupled pulmonary circulation, left atrium and mitral valve model. The left atrium and mitral valve are modelled with physiologically realistic three-dimensional geometries, fibre-reinforced hyperelastic materials and fluid–structure interaction, and the pulmonary vessels are modelled as one-dimensional network ended with structured trees, with specified vessel geometries and wall material properties. This new coupled model reveals some interesting results which could be of diagnostic values. For example, the wave propagation through the pulmonary vasculature can lead to different arrival times for the second systolic flow wave (S2 wave) among the pulmonary veins, forming vortex rings inside the left atrium. In the case of acute mitral regurgitation, the left atrium experiences an increased energy dissipation and pressure elevation. The pulmonary veins can experience increased wave intensities, reversal flow during systole and increased early-diastolic flow wave (D wave), which in turn causes an additional flow wave across the mitral valve (L wave), as well as a reversal flow at the left atrial appendage orifice. In the case of atrial fibrillation, we show that the loss of active contraction is associated with a slower flow inside the left atrial appendage and disappearances of the late-diastole atrial reversal wave (AR wave) and the first systolic wave (S1 wave) in pulmonary veins. The haemodynamic changes along the pulmonary vessel trees on different scales from microscopic vessels to the main pulmonary artery can all be captured in this model. The work promises a potential in quantifying disease progression and medical treatments of various pulmonary diseases such as the pulmonary hypertension due to a left heart dysfunction.

Impact of left atrial appendage location on risk of thrombus formation in patients with atrial fibrillation

Most strokes in patients with atrial fibrillation (AF) are thought to arise from thrombus formation in the left atrial appendage (LAA). Assessing the hemodynamics in LAA and left atrium (LA) may provide some insights in the evaluation of the risk of thrombus formation. This study aims to find out the impact of different LAA locations with respect of LA on the risk of thrombus formation within LAA in patients with AF. Three different LAA locations at LA were modeled and a fully coupled fluid-structure interaction analysis was performed. A discrete phase method was used for particle residence analysis to evaluate risk of the thrombus formation. The results showed that LAA positions on the LA affected the LAA flow velocity distribution, passive contraction ability, and particle residence. In particular, the left pulmonary veins (PVs) had a greater influence on the LAA hemodynamics than the right PVs. The LAA had the lowest contractibility when it was located between left superior and left inferior PVs, and in this case, a larger number of particles were resided, which indicated a higher risk of thrombus formation. The present work provides a quantitative way to evaluate the risk of thrombus formation within LAA in patients with AF.

In Vitro Clot Trapping Efficiency of the FDA Generic Inferior Vena Cava Filter in an Anatomical Model: An Experimental Fluid-Structure Interaction Benchmark

PURPOSE

Robust experimental data for performing validation of fluid-structure interaction (FSI) simulations of the transport of deformable solid bodies in internal flow are currently lacking. This in vitro experimental study characterizes the clot trapping efficiency of a new generic conical-type inferior vena cava (IVC) filter in a rigid anatomical model of the IVC with carefully characterized test conditions, fluid rheological properties, and clot mechanical properties.

METHODS

Various sizes of spherical and cylindrical clots made of synthetic materials (nylon and polyacrylamide gel) and bovine blood are serially injected into the anatomical IVC model under worst-case exercise flow conditions. Clot trapping efficiencies and their uncertainties are then quantified for each combination of clot shape, size, and material.

RESULTS

Experiments reveal the clot trapping efficiency increases with increasing clot diameter and length, with trapping efficiencies ranging from as low as approximately 42% for small 3.2 mm diameter spherical clots up to 100% for larger clot sizes. Because of the asymmetry of the anatomical IVC model, the data also reveal the iliac vein of clot origin influences the clot trapping efficiency, with the trapping efficiency for clots injected into the left iliac vein up to a factor of 7.5 times greater than that for clots injected into the right iliac (trapping efficiencies of approximately 10% versus 75%, respectively).

CONCLUSION

Overall, this data set provides a benchmark for validating simulations predicting IVC filter clot trapping efficiency and, more generally, low-Reynolds number FSI modeling.

Computational investigation of left ventricular hemodynamics following bioprosthetic aortic and mitral valve replacement

The left ventricle of the heart is a fundamental structure in the human cardiac system that pumps oxygenated blood into the systemic circulation. Several valvular conditions can cause the aortic and mitral valves associated with the left ventricle to become severely diseased and require replacement. However, the clinical outcomes of such operations, specifically the postoperative ventricular hemodynamics of replacing both valves, are not well understood. This work uses computational fluid-structure interaction (FSI) to develop an improved understanding of this effect by modeling a left ventricle with the aortic and mitral valves replaced with bioprostheses. We use a hybrid Arbitrary Lagrangian-Eulerian/immersogeometric framework to accommodate the analysis of cardiac hemodynamics and heart valve structural mechanics in a moving fluid domain. The motion of the endocardium is obtained from a cardiac biomechanics simulation and provided as an input to the proposed numerical framework. The results from the simulations in this work indicate that the replacement of the native mitral valve with a tri-radially symmetric bioprosthesis dramatically changes the ventricular hemodynamics. Most significantly, the vortical motion in the left ventricle is found to reverse direction after mitral valve replacement. This study demonstrates that the proposed computational FSI framework is capable of simulating complex multiphysics problems and can provide an in-depth understanding of the cardiac mechanics.

A Bayesian Framework to Estimate Fluid and Material Parameters in Micro-swimmer Models

To advance our understanding of the movement of elastic microstructures in a viscous fluid, techniques that utilize available data to estimate model parameters are necessary. Here, we describe a Bayesian uncertainty quantification framework that is highly parallelizable, making parameter estimation tractable for complex fluid-structure interaction models. Using noisy in silico data for swimmers, we demonstrate the methodology's robustness in estimating fluid and elastic swimmer parameters, along with their uncertainties. We identify correlations between model parameters and gain insight into emergent swimming trajectories of a single swimmer or a pair of swimmers. Our proposed framework can handle data with a spatiotemporal resolution representative of experiments, showing that this framework can be used to aid in the development of artificial micro-swimmers for biomedical applications, as well as gain a fundamental understanding of the range of parameters that allow for certain motility patterns.

Biomechanical effects of the novel series LVAD on the aortic valve

BACKGROUND AND OBJECTIVE

The series type of LVAD (i.e., BJUT-II VAD) is a novel left ventricular assist device, whose effects on the aortic valve remain unclear.

METHODS

The biomechanical effects of BJUT-II VAD on the aortic valve were investigated by using a fluid-structure interaction method. The geometric model of BJUT-II VAD was virtually implanted into the ascending aorta to generate the realistic flow pattern for the aortic valve (i.e., support). In addition, the biomechanical states of the aortic valve without BJUT-II VAD support was computed as control (i.e., control case).

RESULTS

Results demonstrated that the biomechanical effects of BJUT-II VAD were quite different from that resulting from traditional "bypass LVAD." Compared with those in the control case, BJUT-II VAD support could significantly reduce the stress load of the leaflet (maximum stress, 0.5 MPa in the control case vs. 0.12 MPa in the support case). Similarly, the rapid valve opening time (100 ms in the control case vs. 175 ms in the support case) and rapid valve closing time (50 ms in the control case vs. 150 ms in the support case) in the support case were obviously longer than those in the control case. Moreover, BJUT-II VAD support reduced retrograde blood flow during the diastolic phase and significantly changed the distribution of WSS of the leaflets.

CONCLUSIONS

In summary, while unloading the left ventricle, BJUT-II VAD could provide beneficial biomechanical states for the aortic leaflets, thereby reducing the risk of aortic valve disease.

A numerical simulation method of natural fragment formation and injury to human thorax

OBJECTIVE

Fragment injury is a type of blast injury that is becoming more and more common in military campaigns and terrorist attacks. Numerical simulation methods investigating the formation of natural fragments and injuries to biological targets are expected to be developed.

METHODS

A cylindrical warhead model was established and the formation process of natural fragments was simulated using the approach of tied nodes with failure through the explicit finite element (FE) software of LS-DYNA. The interaction between the detonation product and the warhead shell was simulated using the fluid-structure interaction algorithm. A method to simulate the injury of natural fragments to a biological target was presented by transforming Lagrange elements into smooth particle hydrodynamics (SPH) particles after the natural fragments were successfully formed. A computational model of the human thorax was established to simulate the injury induced by natural fragments by the node-to-surface contact algorithm with erosion.

RESULTS

The discontinuous velocities of the warhead shell at different locations resulted in the formation of natural fragments with different sizes. The velocities of natural fragments increased rapidly at the initial stage and slowly after the warhead shell fractured. The initial velocities of natural fragments at the central part of the warhead shell were the largest, whereas those at both ends of the warhead shell were the smallest. The natural fragments resulted in bullet holes that were of the same shape as that of the fragments but slightly larger in size than the fragments in the human thorax after they penetrated through. Stress waves propagated in the ribs and enhanced the injury to soft tissues; additionally, ballistic pressure waves ahead of the natural fragments were also an injury factor to the soft tissues.

CONCLUSION

The proposed method is effective in simulating the formation of natural fragments and their injury to biological targets. Moreover, this method will be beneficial for simulating the combined injuries of natural fragments and shock waves to biological targets.

Patient-Specific Bicuspid Aortic Valve Biomechanics: A Magnetic Resonance Imaging Integrated Fluid-Structure Interaction Approach

Congenital bicuspid aortic valve (BAV) consists of two fused cusps and represents a major risk factor for calcific valvular stenosis. Herein, a fully coupled fluid-structure interaction (FSI) BAV model was developed from patient-specific magnetic resonance imaging (MRI) and compared against in vivo 4-dimensional flow MRI (4D Flow). FSI simulation compared well with 4D Flow, confirming direction and magnitude of the flow jet impinging onto the aortic wall as well as location and extension of secondary flows and vortices developing at systole: the systolic flow jet originating from an elliptical 1.6 cm² orifice reached a peak velocity of 252.2 cm/s, 0.6% lower than 4D Flow, progressively impinging on the ascending aorta convexity. The FSI model predicted a peak flow rate of 22.4 L/min, 6.7% higher than 4D Flow, and provided BAV leaflets mechanical and flow-induced shear stresses, not directly attainable from MRI. At systole, the ventricular side of the non-fused leaflet revealed the highest wall shear stress (WSS) average magnitude, up to 14.6 Pa along the free margin, with WSS progressively decreasing towards the belly. During diastole, the aortic side of the fused leaflet exhibited the highest diastolic maximum principal stress, up to 322 kPa within the attachment region. Systematic comparison with ground-truth non-invasive MRI can improve the computational model ability to reproduce native BAV hemodynamics and biomechanical response more realistically, and shed light on their role in BAV patients' risk for developing complications; this approach may further contribute to the validation of advanced FSI simulations designed to assess BAV biomechanics.

A computational study of amoeboid motility in 3D: the role of extracellular matrix geometry, cell deformability, and cell-matrix adhesion

Amoeboid cells often migrate using pseudopods, which are membrane protrusions that grow, bifurcate, and retract dynamically, resulting in a net cell displacement. Many cells within the human body, such as immune cells, epithelial cells, and even metastatic cancer cells, can migrate using the amoeboid phenotype. Amoeboid motility is a complex and multiscale process, where cell deformation, biochemistry, and cytosolic and extracellular fluid motions are coupled. Furthermore, the extracellular matrix (ECM) provides a confined, complex, and heterogeneous environment for the cells to navigate through. Amoeboid cells can migrate without significantly remodeling the ECM using weak or no adhesion, instead utilizing their deformability and the microstructure of the ECM to gain enough traction. While a large volume of work exists on cell motility on 2D substrates, amoeboid motility is 3D in nature. Despite recent progress in modeling cellular motility in 3D, there is a lack of systematic evaluations of the role of ECM microstructure, cell deformability, and adhesion on 3D motility. To fill this knowledge gap, here we present a multiscale, multiphysics modeling study of amoeboid motility through 3D-idealized ECM. The model is a coupled fluid‒structure and coarse-grain biochemistry interaction model that accounts for large deformation of cells, pseudopod dynamics, cytoplasmic and extracellular fluid motion, stochastic dynamics of cell-ECM adhesion, and microstructural (pore-scale) geometric details of the ECM. The key finding of the study is that cell deformation and matrix porosity strongly influence amoeboid motility, while weak adhesion and microscale structural details of the ECM have secondary but subtle effects.

The effect of the entry and re-entry size in the aortic dissection: a two-way fluid-structure interaction simulation

Aortic dissection (AD) is one of the most catastrophic cardiovascular diseases. AD occurs when a layer inside the aorta is disrupted and gives rise to the formation of a true lumen and a false lumen. These lumens can be connected through tears in the intimal flap which are known as entries. Despite being known for about two centuries, the effects of many factors on the morbidity and mortality of this disease are still unknown. As the blood interaction with the aorta is crucial in the severity and the progression of the aortic dissection, a biomechanical approach is chosen to investigate the influence of different morphologies on the severity of this disease. Using the finite element method (FEM) and the fluid-structure interaction (FSI) approach, we have evaluated the blood flow characteristics along the diseased aorta, in conjunction with the deformation of the aortic wall. In this study, an idealized geometry of a dissected descending aorta (type B) with two entries has been studied. The values for the diameter of the entry tear were chosen to be 5 mm and 10 mm. Therefore, a total of four conditions were investigated. According to our results, the retrograde flow through the proximal tear is dependent on the size of the distal re-entry and vice versa. Our results revealed that when both entry and re-entry tears are 10 mm in diameter, the flow passes through the true and false lumens with smaller resistance, resulting in a smaller flutter of the intimal flap, and therefore more stable intimal flap. Major oscillation frequencies of 2.5 Hz and 7.4 Hz were observed for the oscillation of the intimal flap, and amplitudes of the waves with higher frequencies were negligible.

Biomechanical Assessment of Bicuspid Aortic Valve Phenotypes: A Fluid-Structure Interaction Modelling Approach

PURPOSE

Bicuspid aortic valve (BAV) is a congenital heart malformation with phenotypic heterogeneity. There is no prior computational study that assesses the haemodynamic and valve mechanics associated with BAV type 2 against a healthy tricuspid aortic valve (TAV) and other BAV categories.

METHODS

A proof-of-concept study incorporating three-dimensional fluid-structure interaction (FSI) models with idealised geometries (one TAV and six BAVs, namely type 0 with lateral and anterior-posterior orientations, type 1 with R-L, N-R and N-L leaflet fusion and type 2) has been developed. Transient physiological boundary conditions have been applied and simulations were run using an Arbitrary Lagrangian-Eulerian formulation.

RESULTS

Our results showed the presence of abnormal haemodynamics in the aorta and abnormal valve mechanics: type 0 BAVs yielded the best haemodynamical and mechanical outcomes, but cusp stress distribution varied with valve orifice orientation, which can be linked to different cusp calcification location onset; type 1 BAVs gave rise to similar haemodynamics and valve mechanics, regardless of raphe position, but this position altered the location of abnormal haemodynamic features; finally, type 2 BAV constricted the majority of blood flow, exhibiting the most damaging haemodynamic and mechanical repercussions when compared to other BAV phenotypes.

CONCLUSION

The findings of this proof-of-concept work suggest that there are specific differences across haemodynamics and valve mechanics associated with BAV phenotypes, which may be critical to subsequent processes associated with their pathophysiology processes.

Modelling secondary lymphatic valves with a flexible vessel wall: how geometry and material properties combine to provide function

A three-dimensional finite-element fluid/structure interaction model of an intravascular lymphatic valve was constructed, and its properties were investigated under both favourable and adverse pressure differences, simulating valve opening and valve closure, respectively. The shear modulus of the neo-Hookean material of both vascular wall and valve leaflet was varied, as was the degree of valve opening at rest. Also investigated was how the valve characteristics were affected by prior application of pressure inflating the whole valve. The characteristics were parameterised by the volume flow rate through the valve, the hydraulic resistance to flow, and the maximum sinus radius and inter-leaflet-tip gap on the plane of symmetry bisecting the leaflet, all as functions of the applied pressure difference. Maximum sinus radius on the leaflet-bisection plane increased with increasing pressure applied to either end of the valve segment, but also reflected the non-circular deformation of the sinus cross section caused by the leaflet, such that it passed through a minimum at small favourable pressure differences. When the wall was stiff, the inter-leaflet gap increased sigmoidally during valve opening; when it was as flexible as the leaflet, the gap increased more linearly. Less pressure difference was required both to open and to close the valve when either the wall or the leaflet material was more flexible. The degree of bias of the valve characteristics to the open position increased with the inter-leaflet gap in the resting position and with valve inflation pressure. The characteristics of the simulated valve were compared with those specified in an existing lumped-parameter model of one or more collecting lymphangions and used to estimate a revised value for the constant in that model which controls the rate of valve opening/closure with variation in applied pressure difference. The effects of the revised value on the lymph pumping efficacy predicted by the lumped-parameter model were evaluated.

A partition of unity approach to fluid mechanics and fluid-structure interaction

For problems involving large deformations of thin structures, simulating fluid-structure interaction (FSI) remains a computationally expensive endeavour which continues to drive interest in the development of novel approaches. Overlapping domain techniques have been introduced as a way to combine the fluid-solid mesh conformity, seen in moving-mesh methods, without the need for mesh smoothing or re-meshing, which is a core characteristic of fixed mesh approaches. In this work, we introduce a novel overlapping domain method based on a partition of unity approach. Unified function spaces are defined as a weighted sum of fields given on two overlapping meshes. The method is shown to achieve optimal convergence rates and to be stable for steady-state Stokes, Navier-Stokes, and ALE Navier-Stokes problems. Finally, we present results for FSI in the case of 2D flow past an elastic beam simulation. These initial results point to the potential applicability of the method to a wide range of FSI applications, enabling boundary layer refinement and large deformations without the need for re-meshing or user-defined stabilization.

Fluid-structure interaction simulations of patient-specific aortic dissection

Credible computational fluid dynamic (CFD) simulations of aortic dissection are challenging, because the defining parallel flow channels-the true and the false lumen-are separated from each other by a more or less mobile dissection membrane, which is made up of a delaminated portion of the elastic aortic wall. We present a comprehensive numerical framework for CFD simulations of aortic dissection, which captures the complex interplay between physiologic deformation, flow, pressures, and time-averaged wall shear stress (TAWSS) in a patient-specific model. Our numerical model includes (1) two-way fluid-structure interaction (FSI) to describe the dynamic deformation of the vessel wall and dissection flap; (2) prestress and (3) external tissue support of the structural domain to avoid unphysiologic dilation of the aortic wall and stretching of the dissection flap; (4) tethering of the aorta by intercostal and lumbar arteries to restrict translatory motion of the aorta; and a (5) independently defined elastic modulus for the dissection flap and the outer vessel wall to account for their different material properties. The patient-specific aortic geometry is derived from computed tomography angiography (CTA). Three-dimensional phase contrast magnetic resonance imaging (4D flow MRI) and the patient's blood pressure are used to inform physiologically realistic, patient-specific boundary conditions. Our simulations closely capture the cyclical deformation of the dissection membrane, with flow simulations in good agreement with 4D flow MRI. We demonstrate that decreasing flap stiffness from [Formula: see text] to [Formula: see text] kPa (a) increases the displacement of the dissection flap from 1.4 to 13.4 mm, (b) decreases the surface area of TAWSS by a factor of 2.3, (c) decreases the mean pressure difference between true lumen and false lumen by a factor of 0.63, and (d) decreases the true lumen flow rate by up to 20% in the abdominal aorta. We conclude that the mobility of the dissection flap substantially influences local hemodynamics and therefore needs to be accounted for in patient-specific simulations of aortic dissection. Further research to accurately measure flap stiffness and its local variations could help advance future CFD applications.

Immersogeometric fluid-structure interaction modeling and simulation of transcatheter aortic valve replacement

The transcatheter aortic valve replacement (TAVR) has emerged as a minimally invasive alternative to surgical treatments of valvular heart disease. TAVR offers many advantages, however, the safe anchoring of the transcatheter heart valve (THV) in the patients anatomy is key to a successful procedure. In this paper, we develop and apply a novel immersogeometric fluid-structure interaction (FSI) framework for the modeling and simulation of the TAVR procedure to study the anchoring ability of the THV. To account for physiological realism, methods are proposed to model and couple the main components of the system, including the arterial wall, blood flow, valve leaflets, skirt, and frame. The THV is first crimped and deployed into an idealized ascending aorta. During the FSI simulation, the radial outward force and friction force between the aortic wall and the THV frame are examined over the entire cardiac cycle. The ratio between these two forces is computed and compared with the experimentally estimated coefficient of friction to study the likelihood of valve migration.

Numerical analysis of aortic hemodynamics under the support of venoarterial extracorporeal membrane oxygenation and intra-aortic balloon pump

BACKGROUND AND OBJECTIVE

A gap still exists in the hemodynamic effect of intra-aortic balloon pump (IABP), venoarterial extracorporeal membrane oxygenation (VA-ECMO), and VA-ECMO plus IABP on the blood perfusion of the coronary artery, brain, and lower limb; the relation between heart flow and ECMO flow; and the wall stress of vessels.

METHODS

A finite-element model of the aorta, ECMO, and IABP was proposed to calculate the mechanical response via fluid-structure interaction. Heart failure (HF), IABP, ECMO, and ECMO plus IABP were utilized to study the effect of support models.

RESULTS

For the pressure curve, VA-ECMO weakened the dicrotic notch of pressure compared with HF and the pulsatile index (0.494 vs. 0.706 vs. 0.471 vs. 0.613). IABP, ECMO, and ECMO plus IABP increased the perfusion of the coronary, brain, and renal artery compared with HF. However, ECMO and ECMO plus IABP clearly reduced the blood flow of the left arteria femoralis compared to that of the right arteria femoralis (ECMO: 194.04 vs. 730.80 mL/min; ECMO plus IABP: 342.15 vs. 947.22 mL/min). In addition, the flow of ECMO accessed the renal artery more than the left ventricular flow. Greater ventricular flow perfused to the renal artery at a diastolic period for ECMO plus IABP, especially at the time points of 2.192 s and 2.304 s. Compared to the velocity distribution with ECMO, the flow of the right arteria femoralis was increased in the process of IABP-on. According to these four cases, the stress of the vascular wall was increased for ECMO support at the systolic period. The peak wall stress of ECMO is increased by 20% at 1.68 s.

CONCLUSIONS

ECMO plus IABP is more conducive to the blood supply than other cases from the result of numerical simulation. The location of blood intersection was generated in the region of the renal artery, which is estimated carefully.

Effects of Normal Variation in the Rotational Position of the Aortic Root on Hemodynamics and Tissue Biomechanics of the Thoracic Aorta

PURPOSE

Variation in the rotational position of the aortic root relative to the left ventricle is present in normal trileaflet aortic valves. Its impact on the resulting fluid mechanics of blood flow in the thoracic aorta and structural mechanics in the aortic wall are unknown. We aimed to determine the regional hemodynamic and biomechanical differences in different rotational positions of the normal aortic root (clockwise, central, and counterclockwise positions).

METHOD

Cardiac magnetic resonance imaging (CMR) data was acquired from a normal pediatric patient. These were used for reconstruction of the aortic valve and thoracic aorta 3D model. Fluid-structure interaction (FSI) simulations were employed to study the influence of the root rotation with a central position as compared to observed extreme variations. Patient-specific phase-encoding CMR data were used to assess the validity of computed blood flow. The 3D FSI model was coupled with Windkessel boundary conditions that were tuned for physiological pressures. A grid velocity function was adopted for the valve motion during the systolic period.

RESULTS

The largest wall shear stress level is detected in the clockwise positioned aortic root at the sinutubular junction. Two counter-rotating vortex cores are formed within the aortic root of both the central and extreme root configurations, however, in the clockwise root the vortex system becomes more symmetric. This also coincides with more entrainment of the valve jet and more turbulence production along the shear layer.

CONCLUSION

A clockwise rotational position of the aortic root imparts an increased wall shear stress at the sinutubular junction and proximal ascending aorta in comparison to other root rotation positions. This may pose increased risk for dilation of the sinutubular junction and ascending aorta in the patient with a clockwise positioned aortic root compared to other normal positional configurations.

A computational model of a network of initial lymphatics and pre-collectors with permeable interstitium

Initial lymphatic vessels are made up of overlapped endothelial cells that act as unidirectional valves enabling one-way drainage of tissue fluid into the lumen of the initial lymphatics when there is a favourable pressure gradient. Initial lymphatics subsequently drain this fluid into the collecting lymphatics. This paper describes a computational model for a network of passive rat mesenteric lymphatic vessels with sparse secondary valves. The network was simulated with the secondary valves both operational and non-operational. The effects on the cycle-mean outflow-rate from the network of both inflammation and the resistance of the surrounding interstitium were considered. The cycle-mean outflow-rate is sensitive to vessel stiffness. If the influence of primary-valve resistance is reduced relative to that of interstitial resistance and intravascular resistance, there is no absolute advantage of extrinsic pumping, since maximum outflow-rate occurs when vessels are rigid. However, there is relative advantage, in that the outflow-rate at intermediate stiffness is higher with the secondary valves functioning than when they are deactivated. If primary-valve resistance dominates, then extrinsic pumping of non-rigid vessels provides absolute advantage. The nonlinear relation between pressure drop and flow-rate of the endothelial primary valves, combined with downstream compliance and pulsatile external pressure, constitutes a separate mechanism of pumping. By enabling the consideration of interactions between multiple phenomena (primary valves, secondary valves, a real network geometry with multiple branches, deformable vessel walls, interstitial resistance and external pressures), the model offers a perspective for delineating physiological phenomena that have not yet been fully linked to the biomechanics of fluid flow through initial lymphatic networks.

Numerical study on the dynamics of primary cilium in pulsatile flows by the immersed boundary-lattice Boltzmann method

An explicit immersed boundary-lattice Boltzmann method is applied to numerically investigate the dynamics of primary cilium in pulsatile blood flows with two-way fluid-structure interaction considered. To well characterize the effect of cilium basal body on cilium dynamics, the cilium base is modeled as a nonlinear rotational spring attached to the cilium's basal end as proposed by Resnick (Biophys J 109:18-25, 2015. https://doi.org/10.1016/j.bpj.2015.05.031). After several careful validations, the fluid-cilium interaction system is investigated in detail at various pulsatile flow conditions that are characterized by peak Reynolds numbers ([Formula: see text]) and Womersley numbers ([Formula: see text]). The periodic flapping of primary cilium observed in our simulations is very similar to the in vivo ciliary oscillation captured by O'Connor et al. (Cilia 2:8, 2013. https://doi.org/10.1186/2046-2530-2-8). The cilium's dynamics is found to be closely related to the [Formula: see text] and [Formula: see text]. Increase the [Formula: see text] or decrease the [Formula: see text] bring to an increase in the cilium's flapping amplitude, tip angular speed, basal rotation, and maximum tensile stress. It is also demonstrated that by reducing the [Formula: see text] or enhancing the [Formula: see text] to a certain level, one can shift the flapping pattern of cilium from its original two-side one to a one-side one, making the stretch only happen on one particular side. During the flapping process, the location of the maximum tensile stress is not always found at the basal region; instead, it is able to propagate from time to time within a certain distance to the base. Due to the obstruction of the primary cilium, the distribution of wall shear stress no longer remains uniform as in the absence of cilia. It oscillates in space with the minimum magnitude which is always found near where the cilium is located. The presence of cilium also reduces the overall level of wall shear stress, especially at the region near the cilium's anchor point.

Clinical-pathological correlations of BAV and the attendant thoracic aortopathies. Part 1: Pluridisciplinary perspective on their hemodynamics and morphomechanics

Clinical BAV manifestations pertain to faulty aortic valve (AOV) function, the associated aortopathy, and other complications such as endocarditis, thrombosis and thromboembolism. BAV arises during valvulogenesis when 2 of the 3 leaflets/cusps of the AOV are fused together. Ensuing asymmetric BAV morphologies alter downstream ejection jet flow-trajectories. Based on BAV morphologies, ejection-flows exhibit different wall-impingement and scouring patterns in the proximal aorta, with excessive hydrodynamic wall-shear that correlates closely with mural vascular smooth muscle cell and extracellular matrix disruptions, revealing hemodynamic participation in the pathogenesis of BAV-associated aortopathies. Since the embryologic regions implicated in both BAV and aortopathies derive from neural crest cells and second heart field cells, there may exist a common multifactorial/polygenic embryological basis linking the abnormalities. The use of Electronic Health Records - encompassing integrated NGS variant panels and phenotypic data - in clinical studies could speed-up comprehensive understanding of multifactorial genetic-phenotypic and environmental factor interactions. This Survey represents the first in a 2-article pluridisciplinary work. Taken in toto, the series covers hemodynamic/morphomechanical and environmental (milieu intérieur) aspects in Part 1, and molecular, genetic and associated epigenetic aspects in Part 2. Together, Parts 1-2 should serve as a reference-milestone and driver for further pluridisciplinary research and its urgent translations in the clinical setting.

Simultaneous evaluation of plaque stability and ischemic potential of coronary lesions in a fluid-structure interaction analysis

The measurement of fractional flow reserve (FFR) and superficial wall stress (SWS) identifies inducible myocardial ischemia and plaque vulnerability, respectively. A simultaneous evaluation of both FFR and SWS is still lacking, while it may have a major impact on therapy. A new computational model of one-way fluid-structure interaction (FSI) was implemented and used to perform a total of 54 analyses in virtual coronary lesion models, based on plaque compositions, arterial remodeling patterns, and stenosis morphologies under physiological conditions. Due to a greater lumen dilation and more induced strain, FFR in the lipid-rich lesions (0.81 ± 0.15) was higher than that in fibrous lesions (0.79 ± 0.16, P = 0.001) and calcified lesions (0.79 ± 0.16, P = 0.001). Four types of lesions were further defined, based on the combination of cutoff values for FFR (0.80) and maximum relative SWS (30 kPa): The level of risk increased from (1) plaques with mild-to-moderate stenosis but negative remodeling for lipid-rich (Type A: non-ischemic, stable) to (2) lipid-rich plaques with mild-to-moderate stenosis and without-to-positive remodeling (Type B: non-ischemic, unstable) or plaques with severe stenosis but negative remodeling for lipid-rich (Type C: ischemic, stable) to (3) lipid-rich plaques with severe stenosis and without-to-positive remodeling (Type D: ischemic, unstable). The analysis of FSI to simultaneously evaluate inducible myocardial ischemia and plaque stability may be useful to identify coronary lesions at a high risk and to ultimately optimize treatment. Further research is warranted to assess whether a more aggressive treatment may improve the prognosis of patients with non-ischemic, intermediate, and unstable lesions.

A hemodynamic model with a seepage condition and fluid-structure interactions for blood flow in arteries with symmetric stenosis

To strengthen the detailed understanding of arterial stenosis, we construct a novel hemodynamic model. Frequently used symmetric stenosis is employed in this work. Being different from a traditional model, this numerical model adopts microcirculation resistance as an outlet boundary condition, which is called a seepage condition. Meanwhile, fluid-structure interactions are used in the numerical simulation considering the interrelationship of blood and arterial wall. Our results indicate that (i) the region upstream of stenosis experiences very high pressures during cardiac cycles, (ii) pressure drops much faster as the flow moves into the stenotic region, and (iii) high flow velocities and high shear stresses occur in the post-stenosis region. This work provides evidence that there is a strong effect of the function of microcirculation on stenosis. This contributes to evaluating potential stenotic behavior in arteries and is pivotal in guiding disease treatment.

Development of a Computational Method for Simulating Tricuspid Valve Dynamics

Computational modeling can be used to improve understanding of tricuspid valve (TV) biomechanics and supplement knowledge gained from benchtop and large animal experiments. The aim of this study was to develop a computational model of the TV using high resolution micro-computed tomography (μCT) imaging and fluid-structure interaction simulations. A three-dimensional TV model, incorporating detailed leaflet and chordal geometries, was reconstructed from μCT images of an excised porcine TV obtained under diastolic conditions. The leaflets were described using non-linear stress-strain relations and chordal properties were iteratively adjusted until valve closure was obtained. The leaflet coaptation zone obtained from simulation of valve closure was validated against μCT images of the TV captured at peak systole. The computational model was then used to simulate a regurgitant TV morphology and investigate changes in closure dynamics. Overall, the mean stresses in the leaflet belly region and the chordae tendinae of the regurgitant TV were 7% and 3% higher than the same regions of the normal TV. The maximum principal strain in the leaflet belly of the regurgitant TV was also 9% higher than the same regions of the normal TV. It is anticipated that this computational model can be used in future studies for further understanding of TV biomechanics and associated percutaneous repairs.

Investigating heartbeat-related in-plane motion and stress levels induced at the aortic root

Background

The axial motion of aortic root (AR) due to ventricular traction was previously suggested to contribute to ascending aorta (AA) dissection by increasing its longitudinal stress, but AR in-plane motion effects on stresses have never been studied. The objective is to investigate the contribution of AR in-plane motion to AA stress levels.

Methods

The AR in-plane motion was assessed on magnetic resonance imagining data from 25 healthy volunteers as the movement of the AA section centroid. The measured movement was prescribed to the proximal AA end of an aortic finite element model to investigate its influences on aortic stresses. The finite element model was developed from a patient-specific geometry using LS-DYNA solver and validated against the aortic distensibility. Fluid–structure interaction (FSI) approach was also used to simulate blood hydrodynamic effects on aortic dilation and stresses.

Results

The AR in-plane motion was 5.5 ± 1.7 mm with the components of 3.1 ± 1.5 mm along the direction of proximal descending aorta (PDA) to AA centroid and 3.0 ± 1.3 mm perpendicularly under the PDA reference system. The AR axial motion elevated the longitudinal stress of proximal AA by 40% while the corresponding increase due to in-plane motion was always below 5%. The stresses at proximal AA resulted approximately 7% less in FSI simulation with blood flow.

Conclusions

The AR in-plane motion was comparable with the magnitude of axial motion. Neither axial nor in-plane motion could directly lead to AA dissection. It is necessary to consider the heterogeneous pressures related to blood hydrodynamics when studying aortic wall stress levels.

Electronic supplementary material

The online version of this article (10.1186/s12938-019-0632-7) contains supplementary material, which is available to authorized users.

Fluid-structure interaction simulation of the brain-skull interface for acute subdural haematoma prediction

Traumatic brain injury is a leading cause of disability and mortality. Finite element-based head models are promising tools for enhanced head injury prediction, mitigation and prevention. The reliability of such models depends heavily on adequate representation of the brain-skull interaction. Nevertheless, the brain-skull interface has been largely simplified in previous three-dimensional head models without accounting for the fluid behaviour of the cerebrospinal fluid (CSF) and its mechanical interaction with the brain and skull. In this study, the brain-skull interface in a previously developed head model is modified as a fluid-structure interaction (FSI) approach, in which the CSF is treated on a moving mesh using an arbitrary Lagrangian-Eulerian multi-material formulation and the brain on a deformable mesh using a Lagrangian formulation. The modified model is validated against brain-skull relative displacement and intracranial pressure responses and subsequently imposed to an experimentally determined loading known to cause acute subdural haematoma (ASDH). Compared to the original model, the modified model achieves an improved validation performance in terms of brain-skull relative motion and is able to predict the occurrence of ASDH more accurately, indicating the superiority of the FSI approach for brain-skull interface modelling. The introduction of the FSI approach to represent the fluid behaviour of the CSF and its interaction with the brain and skull is crucial for more accurate head injury predictions.