Atherosclerotic plaque rupture is the main cause of many cardiovascular diseases, and biomechanical factors play an important role in the process of plaque rupture. In the study of plaque biomechanics, there are relatively few studies based on fatigue fracture failure theory, and most of them mainly focus on the whole fatigue propagation process from crack initiation to plaque rupture, while there are few studies on the influence of crack on plaque rupture at a certain time in the process of fatigue propagation. In this paper, a two-dimensional plaque model with crack was established. Based on the theory of fracture mechanics and combined with the finite element numerical simulation method, the stress intensity factor (SIF) and related influencing factors at the crack tip in the plaque were studied. The SIF was used to measure the influence of crack on plaque rupture. The results show that the existence of crack can lead to local stress concentration, which increases the risk of plaque rupture. The SIF at the crack tip in the plaque was positively correlated with blood pressure, but negatively correlated with fibrous cap thickness and lipid pool stiffness. The effect of the thickness and angle of lipid pool on the SIF at the crack tip in the plaque was less than 4%, which could be ignored. This study provides a theoretical basis for the risk assessment of plaque rupture with cracks.
An implantable axial blood pump was designed according to the circulation assist requirement of severe heart failure patients of China. The design point was chosen at 3 L/min flow rate with 100 mm Hg pressure rise when the blood pump can provide flow rates of 2-7 L/min. The blood pump with good hemolytic and anti-thrombogenic property at widely operating range was designed by developing a structure that including the spindly rotor impeller structure and the diffuser with splitter blades and cantilevered main blades. Numerical simulation and particle image velocimetry (PIV) experiment were conducted to analyze the hydraulic, flow fields and hemolytic performance of the blood pump. The results showed that the blood pump could provide flow rates of 2-7 L/min with pressure rise of 60.0-151.3 mm Hg when the blood pump rotating from 7 000 to 11 000 r/min. After adding the splitter blades, the separation flow at the suction surface of the diffuser has been reduced efficiently. The cantilever structure changed the blade gap from shroud to hub that reduced the tangential velocity from 6.2 m/s to 4.3-1.1 m/s in blade gap. Moreover, the maximum scalar shear stress of the blood pump was 897.3 Pa, and the averaged scalar shear stress was 37.7 Pa. The hemolysis index of the blood pump was 0.168% calculated with Heuser’s hemolysis model. The PIV and simulated results showed the overall agreement of flow field distribution in diffuser region. The blood damage caused by higher shear stress would be reduced by adopting the spindle rotor impeller and diffuser with splitter blades and cantilevered main blades. The blood could flow smoothly through the axial blood pump with satisfactory hydraulics performance and without separation flow.
Due to their diverse types, complex causes, high incidence, and difficult treatment, lung diseases have become major killers threatening human life and health, and some lung diseases have a significant impact on alveolar morphology and histology. Numerical simulation of alveolar mechanical response, alveolar flow field information, multiphase flow, and material transport based on computational fluid dynamics is of great significance for lung disease diagnosis, clinical treatment, and in vitro experiments. Starting from the simplification and pathological differences of geometric and mechanical models, this paper analyzes and summarizes the conditions and application scenarios of the airflow dynamics calculation method in pulmonary alveoli, to provide a reference for further simulation and application of the alveolar region.
The current finite element analysis of vascular stent expansion does not take into account the effect of the stent release pose on the expansion results. In this study, stent and vessel model were established by Pro/E. Five kinds of finite element assembly models were constructed by ABAQUS, including 0 degree without eccentricity model, 3 degree without eccentricity model, 5 degree without eccentricity model, 0 degree axial eccentricity model and 0 degree radial eccentricity model. These models were divided into two groups of experiments for numerical simulation with respect to angle and eccentricity. The mechanical parameters such as foreshortening rate, radial recoil rate and dog boning rate were calculated. The influence of angle and eccentricity on the numerical simulation was obtained by comparative analysis. Calculation results showed that the residual stenosis rates were 38.3%, 38.4%, 38.4%, 35.7% and 38.2% respectively for the 5 models. The results indicate that the pose has less effect on the numerical simulation results so that it can be neglected when the accuracy of the result is not highly required, and the basic model as 0 degree without eccentricity model is feasible for numerical simulation.
The implantable miniaturized axial blood pump works at a high rotational speed, which increases the risk of blood damage. In this article, we aimed to reduce the possibility of hemolysis and thrombosis by designing a two-stage axial blood pump. Under the operation conditions of flow rate 5 L/min and outlet pressure of 100 mm Hg, we carried out the numerical simulation on the two-stage and single-stage blood pumps to compare the hemolysis and platelet activation state. The results turned out that the hemolysis index of two-stage axial blood pump was better while the platelet activation state was worse than those of single stage design. On the index of hemolysis level and platelet activation state, the design of the two-stage pump with the low and high-head impeller combination was better than the two-stage pump with the equal heads, or the high and low-head impeller combination. In terms of reducing the risk of blood damage for implantable miniaturized axial blood pump, the research result can provide some theoretical basis and new design ideas.
In order to investigate the application of lattice Boltzmann method (LBM) in the numerical simulation of computed tomography angiography-derived fractional flow reserve (FFRCT), an idealized narrowed tube model and two coronary stenosis arterymodels are studied. Based on the open source code library (Palabos), the relative algorithm program in the development environment (Codeblocks) was improved. Through comparing and analyzing the results of FFRCT which is simulated by LBM and finite element analysis software ANSYS, and the feasibility of the numerical simulation of FFRCT by LBM was verified . The results show that the relative error between the results of LBM and finite element analysis software ANSYS is about 1%, which vertifies the feasibility of simulating the coronary FFRCT by LBM. The simulation of this study provides technical support for developing future FFRCT application software, and lays the foundation for the calculation of clinical FFRCT.
The inhalation and deposition of particles in human pulmonary acinus region can cause lung diseases. Numerical simulation of the deposition of inhaled particles in the pulmonary acinus region has offered an effective gateway to the prevention and clinical treatment of these diseases. Based on some important affecting factors such as pulmonary acinar models, model motion, breathing patterns, particulate characteristics, lung diseases and ages, the present research results of numerical simulation in human pulmonary acinus region were summarized and analyzed, and the future development directions were put forward in this paper, providing new insights into the further research and application of the numerical simulation in the pulmonary acinus region.
Atherosclerosis is a complex and multi-factorial pathophysiological process. Researches over the past decades have shown that the development of atherosclerotic vulnerable plaque is closely related to its components, morphology, and stress status. Biomechanical models have been developed by combining with medical imaging, biological experiments, and mechanical analysis, to study and analyze the biomechanical factors related to plaque vulnerability. Numerical simulation could quantify the dynamic changes of the microenvironment within the plaque, providing a method to represent the distribution of cellular and acellular components within the plaque microenvironment and to explore the interaction of lipid deposition, inflammation, angiogenesis, and other processes. Studying the pathological mechanism of plaque development would improve our understanding of cardiovascular disease and assist non-invasive inspection and early diagnosis of vulnerable plaques. The biomechanical models and numerical methods may serve as a theoretical support for designing and optimizing treatment strategies for vulnerable atherosclerosis.
External support stent is a potential means for restricting the deformation and reducing wall stress of the vein graft, thereby improving the long-term patency of the graft in coronary artery bypass surgery. However, there still lacks a theoretical reference for choosing the size of stent based on the diameter of graft. Taking the VEST (venous external support) stent currently used in the clinical practice as the object of study, we constructed three models of VEST stents with different diameters and coupled them respectively to a model of the great saphenous vein graft, and numerically simulated the expansion-contraction process of the vein graft under the constraint of the stents to quantitatively evaluate the influence of stent size on the radial deformation and wall stress of the vein graft. The results showed that while the stent with a small diameter had a high restrictive effect in comparison with larger stents, it led to more severe concentration of wall stress and sharper changes in radial deformation along the axis of the graft, which may have adverse influence on the graft. In order to solve the aforementioned problems, we ameliorated the design of the stent by means of changing the cross-sectional shape of the thick and thin alloy wires from circle into rectangle and square, respectively, while keeping the cross-sectional areas of alloy wires and stent topology unchanged. Further numerical simulations demonstrated that the ameliorated stent evidently reduced the degrees of wall stress concentration and abrupt changes in radial deformation, which may help improve the biomechanical environment of the graft while maintaining the restrictive role of the stent.
A solid-liquid two-phase finite element model of articular cartilage and a microscopic finite element model of chondrocytes were established using the finite element software COMSOL in this study. The purpose of the study is to investigate the mechanics environment and the liquid flow field of the host cartilage chondrocytes in each layer by multi-scale method, under physiological load, with the different elastic modulus of artificial cartilage to repair cartilage defect. The simulation results showed that the uniform elastic modulus of artificial cartilage had different influences on the microenvironment of different layer chondrocytes. With the increase of the elastic modulus of artificial cartilage, the stress of the shallow surface layer and the intermediate layer chondrocytes increased and the stress of deep layer chondrocytes decreased. The flow field direction of the middle layer and the bottom layer of cartilage can also be changed by artificial cartilage implantation, as well as the ways of nourishment supply of the middle layer and underlying chondrocytes change. A barrier to underlying chondrocytes nutrition supply may be caused by this, thus resulting in the uncertainty of the repair results. With cross-scale finite element model simulation analysis of chondrocytes, we can quantitatively evaluate the mechanical environment of chondrocytes in each layer of the host cartilage. It is helpful to assess the clinical effect of cartilage defect reparation more accurately.