Design Innovation Papers

Implementation of an Automated Peripheral Resistance Device in a Mock Circulatory Loop With Characterization of Performance Values Using Simulink Simscape and Parameter Estimation

[+] Author and Article Information
Charles E. Taylor

Department of Biomedical Engineering,
Virginia Commonwealth University,
401 West Main Street RM1229,
Richmond, VA 23220
e-mail: taylorce@mymail.vcu.edu

Gerald E. Miller

Department of Biomedical Engineering,
Virginia Commonwealth University,
401 West Main Street RM1229,
Richmond, VA 23220
e-mail: gemiller@vcu.edu

Manuscript received January 6, 2012; final manuscript received August 7, 2012; published online October 11, 2012. Assoc. Editor: Paul A. Iaizzo.

J. Med. Devices 6(4), 045001 (Oct 11, 2012) (7 pages) doi:10.1115/1.4007458 History: Received January 06, 2012; Revised August 07, 2012

Accurate peripheral resistance simulation in a mock circulatory loop is critical to the evaluation of ventricular assist devices and heart valves. Implementation of an automated device that is capable of accurate resistance settings and precise reproduction of cardiovascular parameters allows for improved construction of experimental conditions within a mock circulatory loop. A mock circulatory loop resistor that employs a proportional valve design is proposed; a piston extending into the flow path to produce a resistance to flow. Real-time position feedback of the piston is used to determine orifice size, providing resolution in the change of resistance over time. Characterization of the physical system with The MathWorks SIMULINKSIMSCAPE™ block set allowed the determination of objective device parameters; the discharge coefficient and critical Reynolds number. The determination of these values was achieved utilizing the SIMULINK™ Parameter Estimation™ tool, experimental data, and a computational plant model of the experimental setup. With this information, an accurate computational model of the resistance device is presented for use in determining resistance settings in silico prior to implementation in the mock circulatory loop. Experimental in vitro trials verified the repeatability of the automated resistor performance by means of a staircase testing of piston position during several different continuous flow rates of a glycerin/water solution.

© 2012 by ASME
Your Session has timed out. Please sign back in to continue.


Widmaier, E. P.Raff, H., and Strang, K. T., eds., 2004, Vander, Sherman, and Luciano's Human Physiology: The Mechanisms of Body Function, 9th ed., McGraw-Hill, New York.
Fung, Y. C., 1993, Biomechanics: Mechanical Properties of Living Tissues, 2nd ed., Springer-Verlag, New York.
Fung, Y. C., 1997, Biomechanics: Circulation, 2nd ed., Springer-Verlag, New York.
Liu, Y., Allaire, P., Wu, Y., Wood, H., and Olsen, D., 2006, “Construction of an Artificial Heart Pump Performance Test System,” Cardiovasc. Eng., 6(4), pp. 151–158. [CrossRef] [PubMed]
Timms, D., Hayne, M., McNeil, K., and Galbraith, A., 2005, “A Complete Mock Circulation Loop for the Evaluation of Left, Right, and Biventricular Assist Devices,” Artif. Organs, 29(7), pp. 564–572. [CrossRef] [PubMed]
Legendre, D., Fonseca, J., Andrade, A., Biscegli, J. F., Manrique, R., Guerrino, D., Prakasan, A. K., Ortiz, J. P., and Lucchi, J. C., 2008, “Mock Circulatory System for the Evaluation of Left Ventricular Assist Devices, Endoluminal Prostheses, and Vascular Diseases,” Artif. Organs, 32(6), pp. 461–467. [CrossRef] [PubMed]
Zappe, R. W., 2004, Valve Selection Handbook, Elsevier, Boston, pp. 73–75.
Gelman, S., 1995, “The Pathophysiology of Aortic Cross-Clamping and Unclamping,” Anesthesiology, 82(4), pp. 1026–1060. [CrossRef] [PubMed]
Watton, J., 2009, Fundamentals of Fluid Power Control, Cambridge University Press, New York, pp. 74–79.
Liptak, B. G., 2005, Instrument Engineer's Handbook, CRC, Boca Raton, pp. 559–566.
Miller, S., 2009, “Modeling Physical Systems as Physical Networks With the Simscape Language,” 6th Vienna International Conference on Mathematical Modeling, Vienna, Austria, February 11–13.
The MathWorks, 2010, “ Simscape 3.0 User's Guide,” The MathWorks, Natick, MA.
Wu, Y., Allaire, P. E., Tao, G., and Olsen, D., 2005, “Modeling, Estimation and Control of Cardiovascular Systems With a Left Ventricular Assist Device,” IEEE Trans. Control Syst. Technol., 15(4), pp. 754–767. [CrossRef]
Johansson, P., and Andersson, B., 2008, “Comparison of Simulation Programs for Supercapacitor Modelling: Model Creation and Verification,” Master's thesis, Chalmers University of Technology, Gothenburg, Sweden.
Merritt, H. E., 1967, Hydraulic Control Systems, Wiley-Interscience, New York, p. 114.
Segur, J. B., and Oberstar, H. E., 1951, “Viscosity of Glycerol and Its Aqueous Solutions,” Ind. Eng. Chem., 43(9), pp. 2117–2120. [CrossRef]
Dow Chemical Corporation, 2011, “Viscosity of Aqueous Glycerine Solutions,” http://www.dow.com/glycerine/resources/table18.htm
Shankar, P. N., and Kumar, M., 1994, “Experimental Determination of the Kinematic Viscosity of Glycerol-Water Mixtures,” Math. Phys. Sci., 444(1922), pp. 573–581. [CrossRef]
Dow Chemical Corporation, 2011, “Density of Glycerine-Water Solutions,” http://www.dow.com/glycerine/resources/table4_3140.htm


Grahic Jump Location
Fig. 2

Cutaway view of the CAD model for the resistor assembly. Components of the device are as follows: (a) linear actuator lead screw, (b) stepper motor, (c) resistor motor frame, (d) rotation constraint bar, (e) piston, (f) O-rings, (g) threaded flow outlet, and (h) resistor block. This depiction is analogous to the annotated photograph of the assembly displayed as Fig. 3.

Grahic Jump Location
Fig. 1

A manually actuated clamp employed as a resistance mechanism by Legendre et al. [6] (left). An example of a Hass (New York, NY) proportional pinch valve as an automated solution to resistance actuation by Timms et al. [5] (right).

Grahic Jump Location
Fig. 3

Resistor as installed in the mock circulatory loop. The elements of the design are as follows: (1) vertical optical rail, (2) linear actuator lead screw, (3) stepper motor, (4) resistor motor frame, (5) rotation constraint bar, (6) piston, (7) flow inlet, (8) 1 of 4 mounting holes for optical rail attachment, (9) flow outlet, (10) potentiometer, and (11) connector with cable leading to driver board (11).

Grahic Jump Location
Fig. 5

Modified mock circulatory loop used in this analysis; the black arrow counter clockwise from pump indicates flow direction. The order of the pump, compliance chamber, resistor, and venous reservoir is consistent with a mock circulatory loop design. The purpose of pressure sensors downstream of the compliance chamber and resistor is to acquire differential pressure across the resistor. The flowmeter has been repositioned downstream of the resistor for better acuity in flow rate measurement for this segment of the loop. The upstream and downstream throttling valves provide fluid flow regulation and backpressure control, respectively [10].

Grahic Jump Location
Fig. 6

SIMULINK model of the experimental system utilizing simscape blocks to model the physical components of the mock circulatory loop. The simulated flow is from left to right, with components similar to those described in the experimental setup and in Fig. 2. The experimental flow sensor data is permuted as a flow source in the computational model; the model needs a volume source. The downstream pressure sensor experimental data is represented as a pressure source in the model, since there needs to be an exit condition for the simulated flow.

Grahic Jump Location
Fig. 7

The top graph depicts the tuning dataset with the result of the model using the tuned values to that of the experimental dataset for the upstream resistor. The flow rate graph (middle graph) indicates that an approximately constant flow was maintained through the resistor movement (bottom graph).

Grahic Jump Location
Fig. 4

LabVIEW interface used to control the resistor and collect data. The slider allows for a position to be selected and the radix below it enables a particular set point to be entered. The read back is the intermittently updated position read back from the position controller. The button at the bottom communicates the new set point to the controller, initiating the movement of the resistor. This panel can be inserted into other LabVIEW interfaces when constructing control panels for mock circulatory loops, where this device would be a subsystem.

Grahic Jump Location
Fig. 8

Performance of the tuned SIMULINK model using input data with different conditions to that of the tuning data. The model simulation at the lower flow rate showed signs of sensitivity to the fluctuations in flow, which was also seen in the tuning data. The similarity of the simulation to the experimental data provides verification that the parameter estimation was a success.




Some tools below are only available to our subscribers or users with an online account.

Related Content

Customize your page view by dragging and repositioning the boxes below.

Related Journal Articles
Related eBook Content
Topic Collections

Sorry! You do not have access to this content. For assistance or to subscribe, please contact us:

  • TELEPHONE: 1-800-843-2763 (Toll-free in the USA)
  • EMAIL: asmedigitalcollection@asme.org
Sign In