ABSTRACT
Objective To explore the application of three parameter identification methods (impedance modulus curve method, impedance component method, and genetic algorithm) in solving parameter identification problem of the 11-element lumped parameter model in the circle of Willis. Methods Using the flow and pressure waveforms of the internal carotid arteries and vertebral arteries on both sides as inlet conditions, parameter values of the model under normal and bilateral vertebral artery stenosis conditions were calculated. The recognition algorithm was verified by using Simulink models, and finally the stability of the recognition algorithm was verified by adding a certain noise to the flow. Results Under normal circumstances, the proximal resistances obtained by the impedance modulus curve method were larger, and the resistances of the anterior communicating artery obtained by the impedance component method were larger. The genetic algorithm could obtain relatively reasonable model parameter values. In the case of vertebral artery stenosis on both sides, the impedance modulus curve method could obviously get the results of the increasement in proximal resistances of the posterior circulation, but the results obtained by the impedance component method and the genetic algorithm mainly lied in that the distal resistance had a larger increase. Conclusions There are still differences between the pressure data calculated by the parameters identified by the above three methods and the actual data, which are considered as modeling errors, source data errors and calculation errors. The impedance modulus curve method has a certain effect in distinguishing changes of the proximal and distal resistances, but there exist large errors in identification of some parameters. The impedance component method can identify the parameters, but this method is unstable with large calculation errors. Genetic algorithm can obtain a better approximate solution, but it has certain problems in distinguishing vertebral artery stenosis. The combination of impedance modulus curve method and genetic algorithm may play a better role in future application of this model for disease diagnosis.
ABSTRACT
We aimed to propose a novel computational approach to predict the electromechanical performance of pre- and post-mitral valve cerclage annuloplasty (MVCA). Furthermore, we tested a virtual estimation method to optimize the left ventricular basement tightening scheme using a pre-MVCA computer model. The present model combines the three-dimensional (3D) electromechanics of the ventricles with the vascular hemodynamics implemented in a lumped parameter model. 3D models of pre- and post-MVCA were reconstructed from the computed tomography (CT) images of two patients and simulated by solving the electromechanical-governing equations with the finite element method. Computed results indicate that reduction of the dilated heart chambers volume (reverse remodeling) appears to be dependent on ventricular stress distribution. Reduced ventricular stresses in the basement after MVCA treatment were observed in the patients who showed reverse remodeling of heart during follow up over 6 months. In the case who failed to show reverse remodeling after MVCA, more virtual tightening of the ventricular basement diameter than the actual model can induce stress unloading, aiding in heart recovery. The simulation result that virtual tightening of the ventricular basement resulted in a marked increase of myocardial stress unloading provides in silico evidence for a functional impact of MVCA treatment on cardiac mechanics and post-operative heart recovery. This technique contributes to establishing a pre-operative virtual rehearsal procedure before MVCA treatment by using patient-specific cardiac electromechanical modeling of pre-MVCA.
Subject(s)
Humans , Computer Simulation , Follow-Up Studies , Heart , Hemodynamics , Mechanics , Methods , Mitral ValveABSTRACT
Abstract Introduction A ventricular assist device (VAD) is an electromechanical pump used to treat heart failures. For designing the physiological control system for a VAD, one needs a mathematical model and its related parameters. This paper presents a characterization procedure for determining the model parameter values of the electrical, mechanical, and hydraulic subsystems of a pediatric Rotary Blood Pump (pRBP). Methods An in vitro test setup consisting of a pRBP prototype, a motor driver module, an acrylic reservoir, mechanical resistance and tubings, pressure and fluid flow sensors, and data acquisition, processing, and visualization system. The proposed procedure requires a set of experimental tests, and a parameter estimation algorithm for determining the model parameters values. Results The operating limits of the pRBP were identified from the steady-state data. The relationship between the pressure head, flow rate, and the rotational speed of the pRBP was found from the static tests. For the electrical and mechanical subsystems, the dc motor model has a viscous friction coefficient that varies nonlinearly with the flow. For the hydraulic subsystem, the pressure head is assumed to be a sum of terms related to the resistance, the inertance, the friction coefficient, and the pump speed. Conclusion The proposed methodology was successfully applied to the characterization of the pRBP. The combined use of static and dynamic tests provided a precise lumped parameter model for representing the pRBP dynamics. The agreement, regarding mean squared deviation, between experimental and simulated results demonstrates the correctness and feasibility of the characterization procedure.
ABSTRACT
Objective To analyze and compare hemodynamic features of two different options for modified B-T shunt (MBTS) surgery, namely end-to-side(ETS) and side-to-side (STS), so as to provide references for clinical treatment of single ventricle heart defect syndrome. MethodsThe real geometric model was reconstructed by medical images obtained from a patient with hypoplastic left heart syndrome (HLHS); MBTS surgery was simulated through virtual operations; a lumped parameter model (LPM) was constructed based on physiological data of the patient; the post-operational boundary conditions of computational fluid dynamics (CFD) models (namely STS model and ETS model) were predicted based on the LPM; numerical simulation was conducted on two CFD models by using finite volume method. Results Flow details and wall shear stress distributions were all obtained for two models. The mean oscillatory shear index (OSI) of ETS model and STS model in part of pulmonary arteries was 3.058×10-3 and 13.624×10-3, respectively, while the energy loss was 116.5 and 94.8 mW, respectively, and blood flow rate ratios of left pulmonary artery to right pulmonary artery (RRPA/LPA) were 0.8 and 1.72, respectively. Conclusions There were nearly no differences between two CFD models in energy loss, which led to a relatively small impact on the surgery. The STS model had a more balanced pulmonary artery blood perfusion and a smaller mean OSI in part of pulmonary arteries, therefore, the STS model was superior to the ETS model. This study provides an important theoretical support and reference for treating patients with HLHS.
ABSTRACT
Objective Modified B-T shunt (MBTS) and central shunt (CS) are two common surgical procedures for the treatment of tetralogy of fallot (TOF). The purpose is to analyze and compare the hemodynamic features of MBTS and CS. Methods 3D anatomy was reconstructed by medical images obtained from a patient with TOF, and two computational models were generated through virtual operations. A lumped parameter model was constructed to predict the post-operational boundary conditions. Computational fluid dynamics (CFD) was performed for the two models. Results A persistent pulmonary blood perfusion was observed in each model both during the systolic phase and diastolic phase, but the maximum velocities in the shunt were different for the two models. The pressure drop of the shunt in CS model was higher than that in MBTS model. The wall shear stress of the shunt in the MBTS model ranged unevenly from 0.025 to 340 Pa, while the wall shear stress in CS model ranged relatively evenly from 32.2 to 72.6 Pa. Conclusions Pulmonary artery blood was increased effectively for both options. The blood perfusion of right upper extremity was decreased in the MBTS model. More blood was directed into the pulmonary artery in CS model. Attention should be paid to the fact that the pressure gradient was large at the proximal anastomosis in both models in clinic. This study provides important theoretical references for surgeons to make choice from the surgery options in the treatment with TOF.
ABSTRACT
Objective Try to set up a nonlinear lumped parameter model of intra-aorta pump by studying the relationship between the pressure difference and the blood flow rate at the head of the pump, so as to predict the hemodynamic parameters of the pump. Methods Only the parameters of the pump, without hemodynamic parameters of the circulating system, were used to esfablish the model. It was composed of a speed-controlled current source representing the flow rate driven by impeller, an internal resistant representing the resistance of the radial clearance, an inductance denoting the inertance of the blood. Results The model could simulate the physiological status of the heart under all the situations from pulmonary congestion to ventricular collapse. The characteristic equation of the pump was derived with parameters determined by experimental data in vitro. Conclusions To verify the accuracy of the model, the prediction value calculated from the model was compared with the one recorded from experiment in vitro. The results showed that the error in between was less than 5%, which indicated that this model could predict the pressure difference of the pump accurately.