Hemodynamics of the Hepatic Venous Three-Vessel Confluences Using Particle Image Velocimetry

A novel potential mechanism for the development of portal vein thrombosis in cirrhosis based on portal hemodynamics

BACKGROUND

Marked changes in hemodynamics have been suggested to be a potential contributing factor to portal vein thrombosis (PVT) development. This study investigated the effect of portal hemodynamics based on the anatomical structure of the portal venous system on PVT development.

METHODS

The morphological features of portal venous system in patients with PVT and those without PVT subgroups were compared. In addition, idealized PV models were established to numerically evaluate the effect of the variation in the angulation of superior mesenteric vein (SMV) and splenic vein (SV) on the hemodynamics of portal venous system.

RESULTS

The angle α (angulation of SMV and SV) in patients with PVT was lower than that in patients without PVT (p < 0.0001), which was the only independent risk factor (odds ratio (OR), 0.90 (95% CI 0.84-0.95); p < 0.0001) for the presence of PVT. With the change in angle α, the flow pattern of blood flow changed greatly, especially the helical flow. When α = 80°, helical flow only appeared at the local PV near the intersection of SMV and SV. When α = 120°, most regions were occupied by the helical flow. In addition, the h₂ gradually increased with increasing α, when α = 80°, h₂ = 12.6 m/s²; when α = 120°, h₂ = 29.3 m/s².

CONCLUSIONS

The angulation of SV and SMV was closely associated with PVT development. Helical flow changed following the varying angulation of SV and SMV. Therefore, angulation of SV and SMV may help to identify high-risk cohorts for future PVT development earlier.

Application of Magnetic Resonance Imaging in Liver Biomechanics: A Systematic Review

MRI-based biomechanical studies can provide a deep understanding of the mechanisms governing liver function, its mechanical performance but also liver diseases. In addition, comprehensive modeling of the liver can help improve liver disease treatment. Furthermore, such studies demonstrate the beginning of an engineering-level approach to how the liver disease affects material properties and liver function. Aimed at researchers in the field of MRI-based liver simulation, research articles pertinent to MRI-based liver modeling were identified, reviewed, and summarized systematically. Various MRI applications for liver biomechanics are highlighted, and the limitations of different viscoelastic models used in magnetic resonance elastography are addressed. The clinical application of the simulations and the diseases studied are also discussed. Based on the developed questionnaire, the papers' quality was assessed, and of the 46 reviewed papers, 32 papers were determined to be of high-quality. Due to the lack of the suitable material models for different liver diseases studied by magnetic resonance elastography, researchers may consider the effect of liver diseases on constitutive models. In the future, research groups may incorporate various aspects of machine learning (ML) into constitutive models and MRI data extraction to further refine the study methodology. Moreover, researchers should strive for further reproducibility and rigorous model validation and verification.

Three-tissue microphysiological system for studying inflammatory responses in gut-liver Axis

The interaction between the gut and the liver, often known as the gut-liver axis, play crucial roles in modulating the body's responses to the xenobiotics as well as progression of diseases. Dysfunction of the axis can cause metabolic disorders as well as obesity, diabetes, and fatty liver disease. During the progression of such diseases, inflammatory responses involving the immune system also play an important part. In this study, we developed a three-tissue microphysiological system (MPS) that can accommodate three different cell types in separated compartments connected via fluidic channels in a microfluidic device. Using computational fluid dynamics, geometry of fluidic channels and flow conditions were optimized for seeding and culturing different cell types in the three-tissue MPS. Caco-2 (gut), RAW264.7 (immune), and HepG2 (liver) cells were seeded and cultured in the chip. Stimulation of the gut cells in the MPS with lipopolysaccharide (LPS) resulted in induction of inflammatory response and production of nitric oxide (NO) in all connected chambers. The anti-inflammatory effect of luteolin was demonstrated. Our study demonstrates that the three-tissue MPS can recapitulate the inflammatory responses involving the gut, liver and immune cells.

MRI-based modeling of spleno-mesenteric confluence flow

Characterization of hepatic blood flow magnitude and distribution can lead to a better understanding of the pathophysiology of liver disease. However, the underlying patterns and dynamics of hepatic flow, such as the helical flow structure that often develops following the spleno-mesenteric confluence (SMC) of the hepatic portal vein, have not yet been comprehensively studied. In this study, we used magnetic resonance image (MRI)-based computational models to study the effects of the helical flow structure and SMC geometry on portal blood flow distribution. Additionally, we examined these flow dynamics with four-dimensional (4D) flow MRI in a group of 12 cirrhotic patients and healthy subjects. A validation model was also created to compare computational data to particle image velocimetry (PIV) data. We found significant correlations between flow structure development, vessel geometry, and blood flow distribution in both virtually modified models and in healthy and cirrhotic subjects. However, the direction of these correlations varied among vessel configuration types. Nonetheless, validation model results displayed good qualitative agreement with computational model data.

Computational Pre-surgical Planning of Arterial Patch Reconstruction: Parametric Limits and In Vitro Validation

Surgical treatment of congenital heart disease (CHD) involves complex vascular reconstructions utilizing artificial and native surgical materials. A successful surgical reconstruction achieves an optimal hemodynamic profile through the graft in spite of the complex post-operative vessel growth pattern and the altered pressure loading. This paper proposes a new in silico patient-specific pre-surgical planning framework for patch reconstruction and investigates its computational feasibility. The proposed protocol is applied to the patch repair of main pulmonary artery (MPA) stenosis in the Tetralogy of Fallot CHD template. The effects of stenosis grade, the three-dimensional (3D) shape of the surgical incision and material properties of the artificial patch are investigated. The release of residual stresses due to the surgical incision and the extra opening of the incision gap for patch implantation are simulated through a quasi-static finite-element vascular model with shell elements. Implantation of different unloaded patch shapes is simulated. The patched PA configuration is pressurized to the physiological post-operative blood pressure levels of 25 and 45 mmHg and the consequent post-operative stress distributions and patched artery shapes are computed. Stress–strain data obtained in-house, through the biaxial tensile tests for the mechanical properties of common surgical patch materials, Dacron, Polytetrafluoroethylene, human pericardium and porcine xenopericardium, are employed to represent the mechanical behavior of the patch material. Finite-element model is experimentally validated through the actual patch surgery reconstructions performed on the 3D printed anatomical stenosis replicas. The post-operative recovery of the initially narrowed lumen area and post-op tortuosity are quantified for all modeled cases. A computational fluid dynamics solver is used to evaluate post-operative pressure drop through the patch-reconstructed outflow tract. According to our findings, the shorter incisions made at the throat result in relatively low local peak stress values compared to other patch design alternatives. Longer cut and double patch cases are the most effective in repairing the initial stenosis level. After the patch insertion, the pressure drop in the artery due to blood flow decreases from 9.8 to 1.35 mmHg in the conventional surgical configuration. These results are in line with the clinical experience where a pressure gradient at or above 50 mmHg through the MPA can be an indication to intervene. The main strength of the proposed pre-surgical planning framework is its capability to predict the intra-operative and post-operative 3D vascular shape changes due to intramural pressure, cut length and configuration, for both artificial and native patch materials.

Endoleak Assessment Using Computational Fluid Dynamics and Image Processing Methods in Stented Abdominal Aortic Aneurysm Models

Endovascular aortic aneurysm repair (EVAR) is a predominant surgical procedure to reduce the risk of aneurysm rupture in abdominal aortic aneurysm (AAA) patients. Endoleak formation, which eventually requires additional surgical reoperation, is a major EVAR complication. Understanding the etiology and evolution of endoleak from the hemodynamic perspective is crucial to advancing the current posttreatments for AAA patients who underwent EVAR. Therefore, a comprehensive flow assessment was performed to investigate the relationship between endoleak and its surrounding pathological flow fields through computational fluid dynamics and image processing. Six patient-specific models were reconstructed, and the associated hemodynamics in these models was quantified three-dimensionally to calculate wall stress. To provide a high degree of clinical relevance, the mechanical stress distribution calculated from the models was compared with the endoleak positions identified from the computed tomography images of patients through a series of imaging processing methods. An endoleak possibly forms in a location with high local wall stress. An improved stent graft (SG) structure is conceived accordingly by increasing the mechanical strength of the SG at peak wall stress locations. The presented analytical paradigm, as well as numerical analysis using patient-specific models, may be extended to other common human cardiovascular surgeries.

NUMERICAL SIMULATION OF HEMODYNAMICS IN PORTAL VEIN WITH THROMBOSIS BY COMPUTATIONAL FLUID DYNAMICS

Portal vein thrombosis (PVT) is an important complication that is associated with cirrhotic portal hypertension. The etiology is as yet unclear but could be closely related to the hemodynamics of the portal vein system. This paper investigated the hemodynamics in the portal vein model, both with and without thrombosis, as well as the effect of obstructions on the hemodynamics of the portal vein system using the computational fluid dynamics (CFD) method. PVT can probably develop in the inlets of the portal vein as well as the left/right branches of the portal vein because the distribution of wall shear stress satisfies the conditions for PVT formation based upon the simulation of the hemodynamics in the normal portal vein model. According to the above results, geometric models for a portal vein with a thrombus were constructed and the influence of different degrees (26%, 39%, 53% and 64%) of obstructions was studied. In the model with the maximum obstruction (64% blocked), the maximum velocity of portal vein (PV) increased up to twice than in the model without thrombosis, and the maximum wall shear stress of PV in the model with thrombosis (64% blocked) increased up to 9.4 Pa, whereas it was only 1.9 Pa in the model without thrombosis (nearly one fifth of the maximum wall shear stress). Excessive wall shear stress may cause mechanical damage to the blood vessels and induce physiological changes.

Effects of intraluminal thrombus on patient-specific abdominal aortic aneurysm hemodynamics via stereoscopic particle image velocity and computational fluid dynamics modeling

The pathology of the human abdominal aortic aneurysm (AAA) and its relationship to the later complication of intraluminal thrombus (ILT) formation remains unclear. The hemodynamics in the diseased abdominal aorta are hypothesized to be a key contributor to the formation and growth of ILT. The objective of this investigation is to establish a reliable 3D flow visualization method with corresponding validation tests with high confidence in order to provide insight into the basic hemodynamic features for a better understanding of hemodynamics in AAA pathology and seek potential treatment for AAA diseases. A stereoscopic particle image velocity (PIV) experiment was conducted using transparent patient-specific experimental AAA models (with and without ILT) at three axial planes. Results show that before ILT formation, a 3D vortex was generated in the AAA phantom. This geometry-related vortex was not observed after the formation of ILT, indicating its possible role in the subsequent appearance of ILT in this patient. It may indicate that a longer residence time of recirculated blood flow in the aortic lumen due to this vortex caused sufficient shear-induced platelet activation to develop ILT and maintain uniform flow conditions. Additionally, two computational fluid dynamics (CFD) modeling codes (Fluent and an in-house cardiovascular CFD code) were compared with the two-dimensional, three-component velocity stereoscopic PIV data. Results showed that correlation coefficients of the out-of-plane velocity data between PIV and both CFD methods are greater than 0.85, demonstrating good quantitative agreement. The stereoscopic PIV study can be utilized as test case templates for ongoing efforts in cardiovascular CFD solver development. Likewise, it is envisaged that the patient-specific data may provide a benchmark for further studying hemodynamics of actual AAA, ILT, and their convolution effects under physiological conditions for clinical applications.

Computational Hemodynamic Modeling of the Cardiovascular System

Computational methods and modeling are widely used in many fields to study the dynamic behaviour of different phenomena. Currently, the use of these models is an accepted practice in the biomedical field. One of the most significant efforts in this direction is applied to the simulation and prediction of pathophysiological conditions that can affect different systems of the human body. In this work, the design and development of a computational model of the human cardiovascular system is proposed. The structure of the model has been built from a physiological base, considering some of the mechanisms associated to the cardiovascular system. Thus, the aim of the model is the prediction, heartbeat by heartbeat, of some hemodynamic variables from the cardiovascular system, in different pathophysiological cardiac situations. A modular approach to development of the model has been considered in order to include new knowledge that could force the model's hemodynamic. The model has been validated comparing the results obtained with hemodynamic values published by other authors. The results show the usefulness and applicability of the model developed. Thus, different simulations of some cardiac pathologies and physical exercise situations are presented, together with the dynamic behaviors of the different variables considered in the model.