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.
Computational fluid dynamics was used to investigate the effect of the pathogenesis of membranous obstruction of inferior vena cava of Budd-Chiari syndrome with various angles between right hepatic vein and inferior vena cava. Mimics software was used to reconstruct the models from magnetic resonance imaging (MRI) angiograms of inferior vena cava, right hepatic vein, middle hepatic vein and left hepatic vein, and 3DMAX was used to construct the models of 30°, 60°, 90° and 120° angles between right hepatic vein and inferior vena cava, which was based on the reconstructed models.The model was conducted with clinical parameters in terms of wall shear stress distribution, static pressure distribution and blood velocity. The results demonstrated that the differences between wall shear stress and static pressure had statistical significance with various angles between right hepatic vein and inferior vena cava by SPSS. The pathogenesis of membranous obstruction of inferior vena cava had a correlation with the angles between right hepatic vein and inferior vena cava.
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 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.
Inhalable particles deposition in the human respiratory system is the main cause of many respiratory and cardiovascular diseases. It plays an important role in related disease prevention and treatment through establishing computer or external entity models to study rules of particle deposition. The paper summarized and analyzed the present research results of various inhalable particle deposition models of upper respiratory tract and pulmonary area, and expounded the application in the areas of disease inducement analysis, drug inhale treatment etc. Based on the review, the paper puts forward the problems and application limitations of present research, especially pointing out future emphasis in development directions. It will have a value of reference guidance for further systematic and in-depth study on the inhalable particle deposition simulation, experiment and application.
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.
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.
This study analyzed the inherent relation between arterial blood mass flow and muscle atrophy of residual limb to provide some necessary information and theoretical support for the clinical rehabilitation of lower limb amputees. Three-dimensional arterial model reconstruction was performed on both intact side and residual limb of a unilateral transfemoral amputee who is the subject. Then hemodynamic calculation was carried out to comparatively analyze the mass flow state at each arterial outlet of both lower extremities. The muscle atrophy ratio of residual limb was calculated by measuring the cross-sectional area of bilateral muscles. Based on the blood supply relationship, the correlation between arterial blood flow reduction ratio and muscle atrophy ratio was discussed. The results showed that the mass flow of superficial femoral arteries and lateral circumflex femoral arteries severely reduced. Meanwhile rectus femoris, vastus lateralis and vastus medialis which were fed by these arteries showed great atrophy too. On the contrary, the mass flow of deep femoral arteries and medial femoral circumflex arteries slightly reduced. Meanwhile gracilis, adductor longus, long head of biceps which were fed by these arteries showed mild atrophy too. These results indicated that there might be a positive and promotion correlation between the muscle atrophy ratio and the blood mass flow reduction ratio of residual limb during rehabilitation.
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.
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.