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.
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.
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.
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.
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.
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.
To solve the problem of stent malapposition of intravascular stents, explore the design method of intravascular body-fitted stent structure and to establish an objective apposition evaluation method, the support and apposition performance of body-fitted stent in the stenotic vessels with different degrees of calcified plaque were simulated and analyzed. The traditional tube-mesh-like stent model was constructed by using computational aided design tool SolidWorks, and based on this model, the body-fitted stent model was designed by means of projection algorithm. Abaqus was used to simulate the crimping-expansion-recoil process of the two stents in the stenotic vessel with incompletely calcified plaque and completely calcified plaque respectively. A comprehensive method for apposition evaluation was proposed considering three aspects such as separation distance, fraction of non-contact area and residual volume. Compared with the traditional stent, the separation distances of the body-fitted stent in the incompletely calcified plaque model and the completely calcified plaque model were decreased by 21.5% and 22.0% respectively, the fractions of non-contact areas were decreased by 11.3% and 11.1% respectively, and the residual volumes were decreased by 93.1% and 92.5% respectively. The body-fitted stent improved the apposition performance and was effective in both incompletely and completely calcified plaque models. The established apposition performance evaluation method of stent considered more geometric factors, and the results were more comprehensive and objective.
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.
A three-dimensional (3D) model of human anterior chamber is reconstructed to explore the effect of different corneal temperatures on the heat transfer in the chamber. Based on the optical coherence tomography imaging of the volunteers with normal anterior chamber, a 3D anterior chamber model was reconstructed by the method of UG parametric design. Numerical simulation of heat transfer and aqueous humor flow in the whole anterior chamber were analyzed by the finite volume methods at different corneal temperatures. The results showed that different corneal temperatures had obvious influence on the temperature distribution and the aqueous flow in the anterior chamber. The temperature distribution is linear and axial symmetrical around the pupillary axis. As the temperature difference increases, the symmetry becomes poorer. Aqueous floated along the warm side and sank along the cool side which forms a vortexing flow. Its velocity increased with the addition of temperature difference. Heat fluxes of cornea, lens andiris were mainly affected by the aqueous velocity. The higher the velocity, the bigger more absolute value of the above-mentioned heat fluxes became. It is practicable to perform the numerical simulation of anterior chamber by the optical coherence tomography imaging. The results are useful for studying the important effect of corneal temperature on the heat transfer and aqueous humor dynamics in the anterior chamber.
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.