Systems Models of the Circula4on BENG 230C Lecture 2
|
|
- Opal Jefferson
- 6 years ago
- Views:
Transcription
1 Systems Models of the Circula4on BENG 230C Lecture 2 Why modeling Enhance insight in physiology Hypothesis genera5on Clinical applica5ons diagnosis training pla7orms for surgeons predict outcomes of surgical interven5ons predict outcomes of therapies
2 Example: Strain in failing hearts Circulatory models
3 Circula5on model elements Compliant vessel qin [ml/s] V [ml] qout [ml/s] Conservation of mass: dv/dt = qin - qout 5 Circula5on model elements Compliant vessel V [ml] p [kpa] C [ml/kpa] Pressure in vessel: p = V/C 6
4 Circula5on model elements Resistance R [kpa*s/ml] pin [kpa] pout [kpa] q [ml/s] Flow over resistance: q = (pin - pout)/r 7 Circula5on model elements Bifurcation q2 q1 q3 Conservation of mass: q1 = q2 + q3 8
5 Circula5on model elements Compliant vessel + Resistance pin qin V R pout qout C pin = V/C qout = (pin-pout)/r dv/dt = qin - qout 9 Circula5on model elements Compliant vessel + Resistance pin qin V R pout qout C By substitution of pin and qout in dv/dt: qin = C*dpin/dt + (pin-pout)/r With: pout = 0 pin = plv qin = qao R = Rper C = Cart this leads to... 10
6 Two- element windkessel model (ODo Frank, 1899) Compliant vessel + Resistance plv qao V Rper qout Cart qao = Cart*dplv/dt + plv/rper 11 windkessel models windkessel = airkettle used in plunger pumps to ensure a steady flow
7 Two- element windkessel model Electric equivalent plv qao V Rper qout Cart I 1 qao I 3 qout U plv I 2 dv/dt Cart Rper 13 Two- element windkessel model Electric equivalent plv qao V Rper qout Cart The 2wk model produces unrealistic aortic pressure waves, especially at higher frequencies 14
8 Three- element windkessel model (Broemser, 1930) Compliant vessel + two Resistances plv qao Zao V part Rper qart Cart p art = V/C art q ao = (p lv -p art )/Z ao q art = p art /R per dv/dt = q ao - q art p lv : LV pressure p art : arterial pressure q ao : aortic flow Z ao : aortic impedance R per : peripheral resistance C art : arterial compliance V: volume of aorta 15 Three- element windkessel model (Broemser, 1930) Compliant vessel + two Resistances plv qao Zao V part Rper qart Cart Produces realistic aortic pressure and flow waves 16
9 Circula5on model elements The LeO Ventricle Time-varying Compliant Cavity Vlv [ml] plv [kpa] Pressure in left ventricle: plv = Vlv/Clv(t) Clv(t) [ml/kpa] Taking into account that V 0 when p=0: plv = (Vlv-Vlv,0)/Clv(t) Or, more commonly, with time-varying elastance E=1/C: plv = Elv(t)*(Vlv-Vlv,0) 17 Circula5on model elements The LeO Ventricle Time-varying Elastance Vlv [ml] plv [kpa] Clv(t) [ml/kpa] We only need to find a function that describes E(t): During diastole: E = Emin At peak systole: E = Emax E(t) = y(t)*(emax - Emin) + Emin When t<dur_systole: y(t) = -cos(2*pi*t/dur_systole)/ else: y(t) = 0 18
10 Circula5on model elements Simple Valve (open/closed) q pin R pout Flow over valve: q = (pin - pout)/r when pin>pout q = 0 when pin<pout 19 Circula5on model elements The LeO Ventricle Time-varying compliant cavity with valves pla qmitral Rmitral Vlv [ml] plv [kpa] Zaorta paorta qaorta Elv(t) [kpa/ml] plv = Elv(t)*(Vlv-Vlv,0) qmitral = (pla - plv)/rmitral qmitral = 0 when pla>plv when pla<plv qaorta = (plv - paorta)/zaorta when plv>paorta qaorta = 0 dvlv/dt = qmitral - qaorta when plv<paorta 20
11 A simple closed-loop circulatory model p lv q ao p art V 1,C 1 Z ao q mit q art R 1 V 2,C 2 R 2 pvein LV connected to 2 compliant vessels: a total of 3 compartments A simple closed-loop circulatory model p lv q ao p art V 1,C 1 Z ao q mit q art R 1 V 2,C 2 1 p art = V 1 /C 1 R 2 if p lv > p art q ao = (p lv (t) - p art )/Z ao else q ao = 0 q art = (p art -p vein )/R 1 dv 1 /dt = q ao - q art p vein 2 p vein = V 2 /C 2 if p vein > p lv q mit = (p vein -p lv (t))/r 2 else q mit = 0 dv 2 /dt = q art - q mit 3 dv LV /dt = q mit - q ao
12 A simple closed-loop circulatory model This leads to a set of 3 ordinary differen5al equa5ons (ODEs): dv 1 /dt = q ao - q art dv 2 /dt = q art - q mit dv LV /dt = q mit - q ao These can be solved numerically using an ODE- solver (e.g. in Matlab) 23 A more detailed closed- loop circulatory model, Circulation model by Lu et al (AJP 2001)
13 Using Matlab ode23 func5on to solve ODEs [t,statevars]=ode23( your_odes,5mespan,y0,op5ons,parameters); input 5mespan: [0 t_end] or [0:dt:t_end] y0: array with ini5al condi5ons op5ons = odeset( RelTol,1e- 4, AbsTol,1e- 6); parameters: array with model parameter values name of func5on your_odes with system of ODEs output: t: 5me array statevars: array with solu5on of state variables Call ode23 from a main program Using Matlab your_odes create M- file your_odes.m with in it: func5on [dy,variables]=your_odes(t,statevars,flag,parameters) your equa5ons go here Variable names same as previous slide, and ignore flag
14 Example of 3- element windkessel model Program wk3_main: % wk3_main dt = 0.01; % 5me step [sec] t_end = 0.6; % end 5me of simula5on [sec] q lv Zao = 7e- 3; % aor5c impedance [kpa*sec/ml] p art Cao = 32.5; % aor5c compliance [ml/kpa] Rper = 0.25; % peripheral resistance [kpa*sec/ml] p lv Z ao % ini5al condi5ons y0(1,1) = 300; % ini5al aor5c volume [ml] % put parameters in array parameters(1) = Zao; parameters(2) = Cao; parameters(3) = Rper; Cao V R per q art op5ons = odeset('reltol',1e- 4,'AbsTol',1e- 3); % Solve the system of ODEs [t,statevars]=ode23( wk3_ode,[0:dt:t_end],y0,op5ons,parameters); % AOer solu5on is obtained, get rest of (non- state) variables values: for i=1:length(t), [dummy,variables(i,:)]=wk3_ode(t(i),statevars(i,:),0,parameters); end Example of 3- element windkessel model Program wk3_ode: func5on [dy,variables]=wk3_ode(t, statevars, flag, parameters) % Retrieve state variables Vao = statevars(1); % Retrieve parameter values Rao=parameters(1); Cao=parameters(2); Rper=parameters(3); % Make up some left ventricular pressure pmax = 16; %[kpa] plv=pmax*(-cos(2*pi*t/0.6)/2+0.5); part = Vao/Cao; if plv > part, qlv = (plv-part)/rao; % ventricular outflow else qlv = 0.0; % one-way valve end qart = part/rper; % arterial flow dy(1,1) = qlv-qart; %dvao/dt % Store non-state variables in array variables, so we have access and can plot these too variables(1) = qlv; variables(2) = qart; variables(3) = plv; variables(4) = part;
15 Example of 3- element windkessel model: results figure(1); clf subplot(311) plot(t,variables(:,[1 2])); ylabel('flows [ml/sec]','fontsize',16); legend('aortic','arterial','fontsize',16); subplot(312) plot(t,variables(:,[3 4])); ylabel('pressures [kpa]','fontsize',16); legend('lv','aortic'); subplot(313) plot(t,statevars(:,1)) xlabel('time [sec]','fontsize',16); ylabel('aortic volume [ml]','fontsize',16); Homework assignment Design and solve this model of the circula5on in Matlab Rap Cap Cvp Rvs Rvp Era Rtric Rpa Rao Ela Rmit Cvs Erv Elv Cas Ras
16 Parameter values % ====================== solution control parameters =============== dt = 0.001; % time step [sec] tend = 3.000; % end time of simulation [sec] % ======================= Heart ===================================== bcl = 0.600; % Basic cycle length [sec] twitchperiod = 0.300; % Duration of systole [sec], so here systole is half the time of the cycle % Left atrium t_av = 0.120; Emaxla = 0.078; Eminla = 0.071; restvlad = 14; % atrioventricular activation delay % maximum elastance of the left atrium % minimum elastance of left atrium % unloaded volume of left atrium % Right atrium Emaxra = 0.03; Eminra = 0.027; restvrad = 14; % maximum elastance of the left atrium % minimum elastance of the left atrium % unloaded volume of left atrium % Left ventricle Emaxlv = 1.5; % kpa/ml; Eminlv = 1/11; % kpa/ml; restvlvd = 26.1; % ml; % Maximum elastance of left ventricle % Minimum elastance of left ventricle % unloaded volume of left ventricle % Right ventricle Emaxrv = 1/3.01; % kpa/ml; Eminrv = 1/32.3; % kpa/ml; restvrvd = 22.3; % ml; % Maximum elastance of right ventricle % Minimum elastance of right ventricle % unloaded volume of right ventricle Parameter values (Cont d) and ini5al condi5ons % ==================== Systemic circulation ========================= R1s = 0.247; % kpa*sec/ml; % Resistance of vessel 1 R2s = 0.051; % kpa*sec/ml; % Resistance of vessel 2 Rlv = 7e-3; % kpa*sec/ml; % Aortic impedance Rmit = 5e-4; % kpa*sec/ml; % Mitral valve resistance C1s = 32.5; % ml/kpa; % Compliance of vessel 1 C2s = 432; % ml/kpa; % Compliance of vessel 2 % ======================= Pulmonary circulation ====================== R1p = 0.005; % kpa*sec/ml; % Resistance of vessel 1 R2p = 0.005; % kpa*sec/ml; % Resistance of vessel 2 Rrv = 2e-3; % kpa*sec/ml; % Arterial impedance of pulmonary artery Rtric = 5e-4; % kpa*sec/ml; % Resistance of tricuspid valve C1p = 41.7; % ml/kpa; % Compliance of vessel 1 C2p = 50; % ml/kpa; % Compliance of vessel 2 % =========================== Initial conditions ====================== i=1; statevar_init(i,1)=45.4;i=i+1; % LV statevar_init(i,1)=352;i=i+1;! % systemic arterial volume statevar_init(i,1)=1.04e3;i=i+1; % systemic venous volume statevar_init(i,1)=35;i=i+1; % RA statevar_init(i,1)=34.7;i=i+1; % RV statevar_init(i,1)=84.9;i=i+1; % pulmonary arterial volume statevar_init(i,1)=93.0;i=i+1; % pulmonary venous volume statevar_init(i,1)=37.3;i=i+1; % LA
17 Parameter values (Cont d) and ini5al condi5ons % ============= Put parameters together in one array ============== i=1; parameters(i)=bcl;i=i+1; parameters(i)=twitchperiod;i=i+1; parameters(i)=t_av;i=i+1; parameters(i)=emaxla;i=i+1; parameters(i)=eminla;i=i+1; parameters(i)=restvlad;i=i+1; parameters(i)=emaxra;i=i+1; parameters(i)=eminra;i=i+1; parameters(i)=restvrad;i=i+1; parameters(i)=r1s;i=i+1; parameters(i)=r2s;i=i+1; parameters(i)=rlv;i=i+1; parameters(i)=rmit;i=i+1; parameters(i)=c1s;i=i+1; parameters(i)=c2s;i=i+1; parameters(i)=restvlvd;i=i+1; parameters(i)=emaxlv;i=i+1; parameters(i)=eminlv;i=i+1; parameters(i)=r1p;i=i+1; parameters(i)=r2p;i=i+1; parameters(i)=rrv;i=i+1; parameters(i)=rtric;i=i+1; parameters(i)=c1p;i=i+1; parameters(i)=c2p;i=i+1; parameters(i)=restvrvd;i=i+1; parameters(i)=emaxrv;i=i+1; parameters(i)=eminrv;i=i+1; options = odeset('reltol',1e-4,'abstol',1e-6); % Solve the system of Ordinary Differential Equations [t,statevars]=ode23('circulation_statevardot',[0:dt:tend],statevar_init,options,parameters); % After solution is obtained, get rest of (non-state) variables values: for i=1:length(t), [dummy,variables(i,:)]=circulation_statevardot(t(i),statevars(i,:),0,parameters); end Programming the ODEs in a different file Program circula-on_statevardot.m function [dy,variables]=circulation_statevardot(t, statevar, flag, parameters) % Retrieve state variables i=1; Vlv = statevar(i); i=i+1; V1s = statevar(i); i=i+1; V2s = statevar(i); i=i+1; Vra = statevar(i); i=i+1; Vrv = statevar(i); i=i+1; V1p = statevar(i); i=i+1; V2p = statevar(i); i=i+1; Vla = statevar(i); i=i+1; % Retrieve parameter values i=1; bcl=parameters(i);i=i+1; twitchperiod=parameters(i);i=i+1; t_av=parameters(i);i=i+1; Emaxla=parameters(i);i=i+1; Eminla=parameters(i);i=i+1; restvlad=parameters(i);i=i+1; Emaxra=parameters(i);i=i+1; Eminra=parameters(i);i=i+1; restvrad=parameters(i);i=i+1; R1s=parameters(i);i=i+1; R2s=parameters(i);i=i+1; Rlv=parameters(i);i=i+1; Rmit=parameters(i);i=i+1; C1s=parameters(i);i=i+1; C2s=parameters(i);i=i+1; restvlvd=parameters(i);i=i+1; Emaxlv=parameters(i);i=i+1; Eminlv=parameters(i);i=i+1; R1p=parameters(i);i=i+1; R2p=parameters(i);i=i+1; Rrv=parameters(i);i=i+1; Rtric=parameters(i);i=i+1; C1p=parameters(i);i=i+1; C2p=parameters(i);i=i+1; restvrvd=parameters(i);i=i+1; Emaxrv=parameters(i);i=i+1; Eminrv=parameters(i);i=i+1; 34
18 ODEs con5nued % ======================= Equations ========================= % ==================== Cardiac elastance ======================= % ========================= Atria ============================ time_a = mod(t+t_av,bcl); % Time for atria, reset to zero for each new heartbeat % Systolic atrial activation function if (time_a<twitchperiod), ya = -cos(2.0*pi*time_a/twitchperiod)/ ;% systole else ya = 0.0; % diastole end % Left atrium Ela = ya*(emaxla-eminla) + Eminla; Pla = Ela*(Vla-restVlad); % Right atrium... right atrial elastance equations go here... % ========================== Ventricles ========================... ventricular elastance plus rest of circulation equations go here More specifically (1) Solve the model from t=0 to t=3 sec in Matlab with ode23 Use subplot to plot the following as a func5on of 5me in figure 1: Plot LV and RV pressures in subplot(311) Plot LV and RV volumes in subplot(312) Plot LV and RV in- and out- flows in subplot(313) In figure 2, plot the pressure- volume diagram for the LV (for all 5 beats) What is the stroke volume and ejec5on frac5on for the last beat?
19 More specifically (2) Next, re- run the model, but at t=3 sec we ll induce an acute leo ventricular infarct. Simulate this by reducing LV maximum elastance by 30%. Stop the simula5on at t=10 sec. In figure 3, plot the pressure- volume diagram for the LV (for all beats from t=0 to t=10 sec). What is the stroke volume and ejec5on frac5on for the last beat? Discuss what you observe. Remember to label your plots appropriately Hand in on paper (due date Tuesday April 12th, 2011) Matlab scripts (40 points) 3 figures (10 points each, total 30) Stroke volumes and ejec5on frac5ons for the normal and infarcted leo ventricle (5 points each, total 20) Discussion for figure 3 (10 points) Also your Matlab scripts to Amran (aasadi@ucsd.edu)
Simulating ventricular elastance with a heart-arterial interaction model
Simulating ventricular elastance with a heart-arterial interaction model Anita Gerstenmayer 1, Bernhard Hametner 2, Stephanie Parragh 1,2, Thomas Weber 3, Siegfried Wassertheurer 2 1 Department for Analysis
More informationCARDIOVASCULAR SIMULATOR
Harvard-MIT Division of Health Sciences and Technology HST.542J: Quantitative Physiology: Organ Transport Systems Instructors: Roger Mark and Jose Venegas CARDIOVASCULAR SIMULATOR Revision: January 4,
More informationMathematical physiology. 1.1 Derive a suitably scaled form of the Michaelis-Menten model for the reaction E + P,, λ = k 2
Problem sheet 1. 1.1 Derive a suitably scaled form of the Michaelis-Menten model for the reaction S + E and show that it depends on the parameters k 1 k 1 C k 2 E + P K = k 1 + k 2 k 1 S 0 λ = k 2 k 1
More informationCardiovascular Fluid Mechanics - lecture notes 8W090 -
Cardiovascular Fluid Mechanics - lecture notes 8W9 - F.N. van de Vosse (23) Eindhoven University of Technology department of Biomedical Engineering Preface i Preface As cardiovascular disease is a major
More informationStructural Identifiability Analysis of a Cardiovascular System Model
Preprints of the 19th World Congress The International Federation of Automatic Control Structural Identifiability Analysis of a Cardio System Model Antoine Pironet Pierre C. Dauby J. Geoffrey Chase James
More informationS#ff ODEs and Systems of ODEs
S#ff ODEs and Systems of ODEs Popula#on Growth Modeling Let the number of individuals in a given area at time t be. At time the number is so that is the number of individuals that have arrived in the area
More informationBiomechanics. Soft Tissue Biomechanics
Biomechanics cross-bridges 3-D myocardium ventricles circulation Image Research Machines plc R* off k n k b Ca 2+ 0 R off Ca 2+ * k on R* on g f Ca 2+ R0 on Ca 2+ g Ca 2+ A* 1 A0 1 Ca 2+ Myofilament kinetic
More informationParameter Identification Methods in a Model of the Cardiovascular System
Parameter Identification Methods in a Model of the Cardiovascular System Antoine Pironet Thomas Desaive Pierre C. Dauby J. Geoffrey Chase Paul D. Docherty GIGA-Cardiovascular Sciences, University of Liège,
More informationMathematical Modelling of the Cardiovascular System Haemodynamics
Mathematical Modelling of the Cardiovascular System Haemodynamics Anabel Hernández-Ramírez 1 Andrés Fraguela-Collar 1 Rafael Lemuz-López 2 1 Benemérita Universidad Autónoma de Puebla Physical and Mathematical
More informationMathematical Model. M. Umar Qureshi, Mitchel J. Colebank, and Mette S. Olufsen
Mathematical Model M. Umar Qureshi, Mitchel J. Colebank, and Mette S. Olufsen Department of Mathematics, North Carolina State University, Raleigh, North Carolina 27695 Friday 7 th September, 2018 The 1D
More informationMECH : a Primer for Matlab s ode suite of functions
Objectives MECH 4-563: a Primer for Matlab s ode suite of functions. Review the fundamentals of initial value problems and why numerical integration methods are needed.. Introduce the ode suite of numerical
More informationBME 419/519 Hernandez 2002
Vascular Biology 2 - Hemodynamics A. Flow relationships : some basic definitions Q v = A v = velocity, Q = flow rate A = cross sectional area Ohm s Law for fluids: Flow is driven by a pressure gradient
More information2803/01 Transport January 2005 Mark Scheme
2803/01 Transport January 2005 ADVICE TO EXAMINERS ON THE ANNOTATION OF SCRIPTS 1. Please ensure that you use the final version of the. You are advised to destroy all draft versions. 2. Please mark all
More informationA RULE-BASED CONTROLLER BASED ON SUCTION DETECTION FOR ROTARY BLOOD PUMPS
A RULE-BASED CONTROLLER BASED ON SUCTION DETECTION FOR ROTARY BLOOD PUMPS by Antonio Luiz S. Ferreira Master in Electrical Engineering Universidade Federal do Maranhão - Brazil, 998 Submitted to the Graduate
More informationA Validation Study of Mathematical Models of the Human Cardiovascular System for Short-term Regulation. Marwan Ahmad Chabbani, Bsc, Msc (Eng)
A Validation Study of Mathematical Models of the Human Cardiovascular System for Short-term Regulation Marwan Ahmad Chabbani, Bsc, Msc (Eng) A thesis submitted for the Degree of Doctor of Philosophy of
More informationThe Time Varying Elastance Model used as a Boundary Condition in Arterial Network Simulations
The Time Varying Elastance Model used as a Boundary Condition in Arterial Network Simulations Knut Petter Maråk Mechanical Engineering Submission date: July 2013 Supervisor: Leif Rune Hellevik, KT Co-supervisor:
More informationREVIEW SIMILITUDE IN THE CARDIOVASCULAR SYSTEM OF MAMMALS
The Journal of Experimental Biology 204, 395 407 (2001) Printed in Great Britain The Company of Biologists Limited 2001 JEB2982 395 REVIEW SIMILITUDE IN THE CARDIOVASCULAR SYSTEM OF MAMMALS THOMAS H. DAWSON*
More informationReceding horizon controller for the baroreceptor loop in a model for the cardiovascular system
SpezialForschungsBereich F 32 Karl Franzens Universität Graz Technische Universität Graz Medizinische Universität Graz Receding horizon controller for the baroreceptor loop in a model for the cardiovascular
More informationDoppler Imaging: The Basics
Doppler Imaging: The Basics Brian D. Coley, MD, FACR Cincinna< Children s Hospital Medical Center Cincinna
More informationESTIMATION AND IDENTIFICATION OF PARAMETERS IN A LUMPED CEREBROVASCULAR MODEL
MATHEMATICAL BIOSCIENCES doi:10.3934/mbe.2009.6.93 AND ENGINEERING Volume 6, Number 1, January 2009 pp. 93 115 ESTIMATION AND IDENTIFICATION OF PARAMETERS IN A LUMPED CEREBROVASCULAR MODEL Scott R. Pope
More informationKamran Ghalili, M.D. Rosamund Irwin, M.D.
Giffen, Todd Michael (29 yrs [3/13/1985]) Sex.M Procedure Notes Procedures filed by Ghalili, Kamran, MD at 10/08/08 1228 Author: Ghalili, Kamran, MD Service: (none) Author Physician Type: Filed: 10/08/08
More informationIntroduction to Signals and Systems General Informations. Guillaume Drion Academic year
Introduction to Signals and Systems General Informations Guillaume Drion Academic year 2017-2018 SYST0002 - General informations Website: http://sites.google.com/site/gdrion25/teaching/syst0002 Contacts:
More informationBIOE 110: Biomedical Physiology for Engineers Spring 2013 Midterm I Solutions Key
BIOE 110: Biomedical Physiology for Engineers Spring 2013 Midterm I Solutions Key QUESTION 1 Consider a chamber (at 25 C) consisting of two 1L solutions of mannose (MW 180 g/mol) separated by a semipermeable
More informationCopyright 2015 by Christina Battista. All Rights Reserved
ABSTRACT BATTISTA, CHRISTINA. Parameter Estimation of Viscoelastic Wall Models in a One-Dimensional Circulatory Network. (Under the direction of Mette S. Olufsen and Mansoor A. Haider.) Flow and pressure
More informationM. DAROWSKI 1, W. KLONOWSKI 1,2, M. KOZARSKI 1, G. FERRARI 3, K. ZIELIŃSKI 1, R. STEPIEN 1,2
M. DAROWSKI 1, W. KLONOWSKI 1,2, M. KOZARSKI 1, G. FERRARI 3, K. ZIELIŃSKI 1, R. STEPIEN 1,2 1 Polish Academy of Sciences (IBIB PAN) Institute of Biocybernetics and Biomedical Engineering Warsaw, Poland,
More informationApplied Mathematical Models in Human Physiology
IN TRANQUILLO MORS IN FLUCTU VITA Applied Mathematical Models in Human Physiology September 12, 2003 J.T. Ottesen, M.S. Olufsen, and J.K. Larsen UNIVERSITAS ROSKILDENSIS BioMath-Group Department of Mathematics
More informationInternational Journal of Scientific & Engineering Research, Volume 4, Issue 7, July-2013 ISSN
648 Local sensitivity analysis of cardiovascular system parameters R. Gul and S. Bernhard, Fachbereich Mathematik, FU Berlin, Germany. Dept. of Electrical Engineering and Information Technology, Pforzheim
More informationA connecting tube, which enables the movement of the liquid between two positions
Chapter 4: Flow In vivo, transportation of materials is controlled by the flow of a fluid such as liquid or gas. When the control of this flow does not work, life activity stagnates. In this chapter, we
More informationMAT 275 Laboratory 4 MATLAB solvers for First-Order IVP
MAT 275 Laboratory 4 MATLAB solvers for First-Order IVP In this laboratory session we will learn how to. Use MATLAB solvers for solving scalar IVP 2. Use MATLAB solvers for solving higher order ODEs and
More informationHemodynamics II. Aslı AYKAÇ, PhD. NEU Faculty of Medicine Department of Biophysics
Hemodynamics II Aslı AYKAÇ, PhD. NEU Faculty of Medicine Department of Biophysics Laplace s Law Relates the pressure difference across a closed elastic membrane on liquid film to the tension in the membrane
More informationIEEE TRANSACTIONS ON CIRCUITS AND SYSTEMS I: REGULAR PAPERS, VOL. 53, NO. 11, NOVEMBER
IEEE TRANSACTIONS ON CIRCUITS AND SYSTEMS I: REGULAR PAPERS, VOL. 53, NO. 11, NOVEMBER 2006 2459 Cycle-Averaged Models of Cardiovascular Dynamics Tushar A. Parlikar, Student Member, IEEE, Thomas Heldt,
More informationA numerical study of the effect of aortic wall compliance and blood flow through the ascending aorta
A numerical study of the effect of aortic wall compliance and blood flow through the ascending aorta Helen Lorna Morrissey Thesis presented in fulfillment blah blah blah for the Degree of BSc in Mechanical
More information10 0 = 1 1 second. Chapter 1
Chapter 1 10 0 = 1 1 second The word second is derived from the Latin word secundus or gradus secundus, which means second step or next step. The Romans divided the daylight time into 12 hours. As a further
More informationToni Lassila, Cristiano Malossi, Matteo Astorino, and Simone Deparis August 3, 2011
Geometrical multiscale model of an idealized left ventricle with fluid-structure interaction effects coupled to a one-dimensional viscoelastic arterial network Toni Lassila, Cristiano Malossi, Matteo Astorino,
More informationBayesian Networks for Cardiovascular Monitoring
Proceedings of the 8th IEEE EMBS Annual International Conference New York City, USA, Aug 3-Sept 3, 6 WeC5.3 Bayesian Networks for Cardiovascular Monitoring Jennifer M. Roberts, ushar A. Parlikar, homas
More informationUNIVERSITÀ DEGLI STUDI DI TRENTO
UNIVERSITÀ DEGLI STUDI DI TRENTO DEPARTMENT OF MATHEMATICS MASTER DEGREE IN MATHEMATICS A computational study of the effect of the guide wire on the pressure drop in coronary stenosis Supervisors: Prof.
More informationProject TOUCAN. A Study of a Two-Can System. Prof. R.G. Longoria Update Fall ME 144L Prof. R.G. Longoria Dynamic Systems and Controls Laboratory
Project TOUCAN A Study of a Two-Can System Prof. R.G. Longoria Update Fall 2009 Laboratory Goals Gain familiarity with building models that reflect reality. Show how a model can be used to guide physical
More informationChapter 9b: Numerical Methods for Calculus and Differential Equations. Initial-Value Problems Euler Method Time-Step Independence MATLAB ODE Solvers
Chapter 9b: Numerical Methods for Calculus and Differential Equations Initial-Value Problems Euler Method Time-Step Independence MATLAB ODE Solvers Acceleration Initial-Value Problems Consider a skydiver
More informationA Delay Recruitment Model of the Cardiovascular Control System
A Delay Recruitment Model of the Cardiovascular Control System A.C. Fowler M.J. McGuinness Submitted to Journal of Mathematical Biology, June 2004. Revised December 2004. Copyright will be owned by Springer.
More informationAn introduction to wave intensity analysis
Med Biol Eng Comput (9) 47:175 188 DOI 1.17/s11517-9-439-y SPECIAL ISSUE - ORIGINAL ARTICLE An introduction to wave intensity analysis Kim H. Parker Received: 17 July 8 / Accepted: 9 January 9 / Published
More informationPressure Estimation in the Systemic Arteries Using a Transfer Function. Carl-Johan Thore
Pressure Estimation in the Systemic Arteries Using a Transfer Function Carl-Johan Thore Division of Mechanics Master Thesis Department of Management and Engineering LIU - IEI - TEK - A - - 07 / 0069 -
More informationModelling and control of the human cardiovascular system
Research Collection Master Thesis Modelling and control of the human cardiovascular system Author(s): Gisler, Stefan Publication Date: 2011 Permanent Link: https://doi.org/10.3929/ethz-a-007207574 Rights
More informationA Note Regarding the Mathematical Treatment of a Class of Steady-State Compartmental Models of the Circulation
Montana Tech Library Digital Commons @ Montana Tech Center for Advanced Mineral and Metallurgical Processing (CAMP) Faculty Scholarship 9-1-2016 A Note Regarding the Mathematical Treatment of a Class of
More informationGCE AS/A level 1072/01 BIOLOGY BY2
Surname Centre Number Candidate Number Other Names 2 GCE AS/A level 1072/01 BIOLOGY BY2 P.M. TUESDAY, 14 January 2014 1 hour 30 minutes For s use Question Maximum Mark 1. 4 2. 12 Mark Awarded 1072 010001
More informationNumerical methods for cardiovascular problems: computational electrocardiology and fluid dynamics in moving domains
POLITECNICO DI MILANO Dipartimento di Matematica F. Brioschi Ph. D. course in Mathematical Engineering XXII cycle Numerical methods for cardiovascular problems: computational electrocardiology and fluid
More information24, B = 59 24, A = 55
Math 128a - Homework 8 - Due May 2 1) Problem 8.4.4 (Page 555) Solution: As discussed in the text, the fourth-order Adams-Bashforth formula is a formula of the type x n+1 = x n + h[af n + Bf n 1 + Cf n
More informationChapter 17 Current and Resistance
Chapter 17 Current and Resistance Current Practical applications were based on static electricity. A steady source of electric current allowed scientists to learn how to control the flow of electric charges
More informationPolitecnico di Milano. Characterization of delayed after-depolarization in extended FitzHugh-Nagumo models. Chiara Lelli
Politecnico di Milano Dipartimento di Matematica F. Brioschi Dottorato in Modelli e Metodi Matematici per l Ingegneria - XXIV Ciclo Characterization of delayed after-depolarization in extended FitzHugh-Nagumo
More informationEdexcel (B) Biology A-level
Edexcel (B) Biology A-level Topic 4: Exchange and Transport Notes Single celled organisms do not have specialised transport systems. Substances can enter the cell by passive transport as diffusion distance
More informationMAT 275 Laboratory 4 MATLAB solvers for First-Order IVP
MATLAB sessions: Laboratory 4 MAT 275 Laboratory 4 MATLAB solvers for First-Order IVP In this laboratory session we will learn how to. Use MATLAB solvers for solving scalar IVP 2. Use MATLAB solvers for
More informationAnimal structure and function
Animal structure and function The nervous system Parts of the nervous system 43C, 44B, 45D Brain structure and function Eyes Retina Neurons: How neurons communicate: Resting potential: The resting
More informationA differentiable, periodic function for pulsatile cardiac output based on heart rate and stroke volume
See discussions, stats, and author profiles for this publication at: https://www.researchgate.net/publication/1895386 A differentiable, periodic function for pulsatile cardiac output based on heart rate
More informationToni Lassila, A. Cristiano I. Malossi, Matteo Astorino, and Simone Deparis
Proceedings of the ECCOMAS Thematic International Conference on Simulation and Modeling of Biological Flows (SIMBIO 211) September 21 23, 211, VUB, Brussels, Belgium Geometrical multiscale model of an
More informationMAT 275 Laboratory 4 MATLAB solvers for First-Order IVP
MAT 75 Laboratory 4 MATLAB solvers for First-Order IVP In this laboratory session we will learn how to. Use MATLAB solvers for solving scalar IVP. Use MATLAB solvers for solving higher order ODEs and systems
More informationExploration, Processing and Visualization of Physiological Signals from the ICU
Exploration, Processing and Visualization of Physiological Signals from the ICU By Carlos A. Renjifo Submitted to the Department of Electrical Engineering and Computer Science in Partial Fulfillment of
More informationa:!ii PACKARD Learning Qualitative Models from Physiological Signals ft3hewlett
ft3hewlett a:!ii PACKARD Learning Qualitative Models from Physiological Signals David T. Hau, Enrico W. Coiera Intelligent Networked Computing Laboratory HP Laboratories Bristol HPL-95-31 March, 1995 patient
More informationME352 Project 4. Internal Combustion Engine Analysis. Jeremy Bourque Lab Division 9 April 24 th, 2018
ME352 Project 4 Internal Combustion Engine Analysis Jeremy Bourque Lab Division 9 April 24 th, 2018 Summary of Report This project aims to model an internal combustion engine from a lawnmower as a slidercrank
More informationModeling radiocarbon in the Earth System. Radiocarbon summer school 2012
Modeling radiocarbon in the Earth System Radiocarbon summer school 2012 Outline Brief introduc9on to mathema9cal modeling Single pool models Mul9ple pool models Model implementa9on Parameter es9ma9on What
More informationRank-Deficient Nonlinear Least Squares Problems and Subset Selection
Rank-Deficient Nonlinear Least Squares Problems and Subset Selection Ilse Ipsen North Carolina State University Raleigh, NC, USA Joint work with: Tim Kelley and Scott Pope Motivating Application Overview
More informationInvestigation Of Multiscale Fluid Structure Interaction Modeling Of Flow In Arterial Systems
University of Central Florida Electronic Theses and Dissertations Masters Thesis (Open Access) Investigation Of Multiscale Fluid Structure Interaction Modeling Of Flow In Arterial Systems 2013 Sebastian
More informationChapter 17. Current and Resistance
Chapter 17 Current and Resistance Electric Current The current is the rate at which the charge flows through a surface Look at the charges flowing perpendicularly through a surface of area A I av The SI
More informationTotal arterial inertance as the fourth element of the windkessel model
Total arterial inertance as the fourth element of the windkessel model NIKOS STERGIOPULOS, 1 BEREND E. WESTERHOF, 2 AND NICO WESTERHOF 3 1 Biomedical Engineering Laboratory, Swiss Federal Institute of
More informationArterial Macrocirculatory Hemodynamics
Arterial Macrocirculatory Hemodynamics 莊漢聲助理教授 Prof. Han Sheng Chuang 9/20/2012 1 Arterial Macrocirculatory Hemodynamics Terminology: Hemodynamics, meaning literally "blood movement" is the study of blood
More information14. FLUID DYNAMICS OF THE HEART AND ITS VALVES
14. FLUID DYNAMICS OF THE HEART AND ITS VALVES Charles S. Peskin and David M. McQueen 14.1 INTRODUCTION Blood is a viscous incompressible fluid which is propelled through the arteries, capillaries, and
More informationUniversity of Washington Department of Chemistry Chemistry 453 Winter Quarter 2013
Lecture 1 3/13/13 University of Washington Department of Chemistry Chemistry 53 Winter Quarter 013 A. Definition of Viscosity Viscosity refers to the resistance of fluids to flow. Consider a flowing liquid
More informationA COMPARISON OF WAVELET AND SIMPLICITY-BASED HEART SOUND AND MURMUR SEGMENTATION METHODS. A Thesis. presented to
A COMPARISON OF WAVELET AND SIMPLICITY-BASED HEART SOUND AND MURMUR SEGMENTATION METHODS A Thesis presented to the Faculty of California Polytechnic State University, San Luis Obispo In Partial Fulfillment
More informationGain-Induced Oscillations in Blood Pressure
arxiv:math/90811v1 [math.ds] 16 Aug 199 Gain-Induced Oscillations in Blood Pressure Roselyn M. Abbiw-Jackson William F. Langford September 14, 199 Abstract Mayer waves are long-period (6 to 1 seconds)
More informationThese videos and handouts are supplemental documents of paper X. Li, Z. Huang. An Inverted Classroom Approach to Educate MATLAB in Chemical Process
These videos and handouts are supplemental documents of paper X. Li, Z. Huang. An Inverted Classroom Approach to Educate MATLAB in Chemical Process Control, Education for Chemical Engineers, 9, -, 7. The
More information*AB121* *28AB12101* Biology. Assessment Unit AS 2 assessing [AB121] THURSDAY 23 JUNE, MORNING. Organisms and Biodiversity *AB121* TIME
Centre Number ADVANCED SUBSIDIARY (AS) General Certificate of Education 2016 Candidate Number Biology Assessment Unit AS 2 assessing Organisms and Biodiversity *AB121* [AB121] *AB121* THURSDAY 23 JUNE,
More informationSynthesising robust and optimal parameters for cardiac pacemakers using symbolic and evolutionary computation techniques
Synthesising robust and optimal parameters for cardiac pacemakers using symbolic and evolutionary computation techniques Marta Kwiatkowska 1, Alexandru Mereacre 1, Nicola Paoletti 1, and Andrea Patanè
More informationAngela Lungu. School of Medicine and Biomedical Sciences Department of Cardiovascular Science Medical Physics Group. Submitted for the degree of PhD
MRI IMAGE BASED MEASUREMENT, MODELLING AND DIAGNOSTIC INTERPRETATION OF PRESSURE AND FLOW IN THE PULMONARY ARTERIES: APPLICATIONS IN PULMONARY HYPERTENSION Angela Lungu School of Medicine and Biomedical
More informationLocalizing Cardiac Structures in Fetal Heart Ultrasound Video
Localizing Cardiac Structures in Fetal Heart Ultrasound Video Christopher P. Bridge 1, Christos Ioannou 2, and J. Alison Noble 1 1 Institute of Biomedical Engineering, University of Oford, Oford, UK, 2
More informationModelling of the Baroreflex-Feedback Mechanism with Time-Delay
Modelling of the Baroreflex-Feedback Mechanism with Time-Delay Johnny T. Ottesen IMFUFA, Roskilde University Postbox 260, DK-4000 Roskilde E-mail: Johnny@mmf.ruc.dk 20th February 1997 Running title: Baroreflex-Feedback
More informationPhysics 584 Computational Methods
Physics 584 Computational Methods Introduction to Matlab and Numerical Solutions to Ordinary Differential Equations Ryan Ogliore April 18 th, 2016 Lecture Outline Introduction to Matlab Numerical Solutions
More informationExample 3.3 Moving-average filter
3.3 Difference Equa/ons for Discrete-Time Systems A Discrete-/me system can be modeled with difference equa3ons involving current, past, or future samples of input and output signals Example 3.3 Moving-average
More information31545 Medical Imaging systems
31545 Medical Imaging systems Lecture 5: Blood flow in the human body Jørgen Arendt Jensen Department of Electrical Engineering (DTU Elektro) Biomedical Engineering Group Technical University of Denmark
More informationNonlinear Drug Infusion System
31 8 Biotechnol. Prog. 1995, 11, 31 8-332 Issues in the Design of a Multirate Model-Based Controller for a Nonlinear Drug Infusion System Ravi Gopinath,t,* B. Wayne Bequette,",? R. J. Roy? and H. Kaufmanll
More informationNoninvasive Estimation of Pulmonary Artery Pressure Using Heart Sound Analysis
Brigham Young University BYU ScholarsArchive All Theses and Dissertations 009-1-07 Noninvasive Estimation of Pulmonary Artery Pressure Using Heart Sound Analysis Aaron W. Dennis Brigham Young University
More informationfailure patients. There are two aspects of this project: one is the statistical analysis of the variables taken from heart failure patients, the other
Mathematical Approaches to Predictive Health Monitoring for Heart Failure Patients Yi-Ju Chao 127 Vincent Hall, Minnesota Center of Industrial Mathematics Abstract This article introduces both statistical
More informationLab 3A: Modeling and Experimentation: Two-Can System
Lab 3A: Modeling and Experimentation: Two-Can System Prof. R.G. Longoria Department of Mechanical Engineering The University of Texas at Austin June 26, 2014 1 Introduction 2 One-Can 3 Analysis 4 Two-Can
More informationAPPM 2360 Project 3 Mathematical Investigation of Cardiac Dynamics
APPM 2360 Project 3 Mathematical Investigation of Cardiac Dynamics Due: Thursday, December 6, 2018 by 4:59 p.m. Submit as a PDF to Assignments on Canvas 1 Introduction Cardiac Arrhythmia, or irregular
More informationIntroduction to System Dynamics
SE1 Prof. Davide Manca Politecnico di Milano Dynamics and Control of Chemical Processes Solution to Lab #1 Introduction to System Dynamics Davide Manca Dynamics and Control of Chemical Processes Master
More informationUse of mathematical modeling to study pressure regimes in normal and Fontan blood flow circulations
arxiv:86.246v [math.ap] 6 Jun 28 Use of mathematical modeling to study pressure regimes in normal and Fontan blood flow circulations J.P. Keener, M. Chugunova 2, R.M. Taranets 3, M.G. Doyle 4 Department
More informationDetermination of pressure data in aortic valves
Determination of pressure data in aortic valves Helena Švihlová a, Jaroslav Hron a, Josef Málek a, K.R.Rajagopal b and Keshava Rajagopal c a, Mathematical Institute, Charles University, Czech Republic
More informationLesson 9: Predator-Prey and ode45
Lesson 9: Predator-Prey and ode45 9.1 Applied Problem. In this lesson we will allow for more than one population where they depend on each other. One population could be the predator such as a fox, and
More informationECE 450 Homework #3. 1. Given the joint density function f XY (x,y) = 0.5 1<x<2, 2<y< <x<4, 2<y<3 0 else
ECE 450 Homework #3 0. Consider the random variables X and Y, whose values are a function of the number showing when a single die is tossed, as show below: Exp. Outcome 1 3 4 5 6 X 3 3 4 4 Y 0 1 3 4 5
More informationBENG 186B Winter 2014 Quiz 3. March 5, NAME (Last, First): This quiz is closed book and closed notes. You may use a calculator for algebra.
BENG 186B Winter 2014 Quiz 3 March 5, 2014 NAME (Last, First): This quiz is closed book and closed notes. You may use a calculator for algebra. Circle your final answers in the space provided; show your
More informationContrast magnetic resonance imaging for cardiovascular measurements an vitro study
Eindhoven University of Technology MASTER Contrast magnetic resonance imaging for cardiovascular measurements an vitro study Amilburu, R. Award date: 2005 Link to publication Disclaimer This document contains
More informationEE 4BD4 Lecture 26. Plethysmography
EE 4BD4 Lecture 26 Plethysmography 1 Definition Measure volume change over time Flow = dvolume/dt 2 Sequence Occlude venous return (raise venous cuff pressure to 50 mm Hg Allow limb segment to fill from
More informationENBE 415 Example Problems Dr. Arthur T. Johnson. Example Calculate the expected times for men swimming 500 and 600 m in competition.
ENBE 415 Example Problems Dr. Arthur T. Johnson Example 1..1 Calculate the expected times for men swimming 500 and 600 m in competition. The Riegel equation (1..1) will be used: t = ax b From Table 1..1,
More informationTAC1 TAC4 TAC7 TAC14 TAC21
Table S1 Gene qrt-pcr primer sequence Amplification efficiency α-sma (FW) 5 - GCCAGTCGCTGTCAGGAACCC -3 (RV) 5 - AGCCGGCCTAGAGCCCA -3 Procollagen-I (FW) 5 - AAGACGGGAGGGCGAGTGCT -3 (RV) 5 - AACGGGTCCCCTTGGGCCTT
More informationStress-Strain Analysis of Abdominal Aortic Wall: A Case of 3D Geometry Simulation
Energy Research Journal 1 (2): 165-170, 2010 ISSN 1949-0151 2010 Science Publications Stress-Strain Analysis of Abdominal Aortic Wall: A Case of 3D Geometry Simulation P. Khamdaengyodtai, P. Sakulchangsatjatai
More informationNumerical Integration of Ordinary Differential Equations for Initial Value Problems
Numerical Integration of Ordinary Differential Equations for Initial Value Problems Gerald Recktenwald Portland State University Department of Mechanical Engineering gerry@me.pdx.edu These slides are a
More informationModelling and simulation of an active fibre for cardiac muscle
Modelling and simulation of an active fibre for cardiac muscle P. Krejci J. Sainte-Marie M. Sorine J.M. Urquiza CRM-38 January 6 Mathematical Institute, Academy of Sciences of the Czech Republic, Žitná
More informationTHE ELASTIC PROPERTIES OF ARTERIES IN RELATION TO THE PHYSIOLOGICAL FUNCTIONS OF THE ARTERIAL SYSTEM
G.-\ STRO ENTEROLOGY Copy right 1967 by The Williams & Wilkins Co. Vol. 52, No.2, Part 2 Printed in U.S.A. THE ELASTIC PROPERTIES OF ARTERIES IN RELATION TO THE PHYSIOLOGICAL FUNCTIONS OF THE ARTERIAL
More informationBasic mechanisms of arrhythmogenesis and antiarrhythmia
EHRA EDUCATIONAL REVIEW AND PREPARATORY COURSE ON INVASIVE CARDIAC ELECTROPHYSIOLOGY EUROPEAN HEART HOUSE, February 2011 Basic mechanisms of arrhythmogenesis and antiarrhythmia Antonio Zaza Università
More informationVector Fields and Solutions to Ordinary Differential Equations using Octave
Vector Fields and Solutions to Ordinary Differential Equations using Andreas Stahel 6th December 29 Contents Vector fields. Vector field for the logistic equation...............................2 Solutions
More informationIntroduction and Background
Chapter 1 1 Introduction and Background The heart is a robust pump capable of beating rhythmically for over 2 ½ billion times in a lifetime. At very early stages of development, the embryonic heart is
More informationB4 Organising animals and plants. Student Book answers. B4.1 The blood. Question Answer Marks Guidance
B4. The blood Any three from: 3 transport of blood cells, transport of dissolved gases, transport of food, transport of hormones, removal of waste products, defence against infection, preventing blood
More information