Modeling of the thermodilution curve for coronary flow assessment. Carel Stasse TU/e BMTE07.22

Similar documents
USING LOWER CLASS WEIGHTS TO CORRECT AND CHECK THE NONLINEARITY OF BALANCES

Technical Vibration - text 2 - forced vibration, rotational vibration

graph of unit step function t

ME 141. Engineering Mechanics

Addition & Subtraction of Polynomials

Reinforcement learning

() t. () t r () t or v. ( t) () () ( ) = ( ) or ( ) () () () t or dv () () Section 10.4 Motion in Space: Velocity and Acceleration

AB for hydrogen in steel is What is the molar flux of the hydrogen through the steel? Δx Wall. s kmole

Ans: In the rectangular loop with the assigned direction for i2: di L dt , (1) where (2) a) At t = 0, i1(t) = I1U(t) is applied and (1) becomes

X-Ray Notes, Part III

Homework 5 for BST 631: Statistical Theory I Solutions, 09/21/2006

IX.1.1 The Laplace Transform Definition 700. IX.1.2 Properties 701. IX.1.3 Examples 702. IX.1.4 Solution of IVP for ODEs 704

LAPLACE TRANSFORMS. 1. Basic transforms

Calculus 241, section 12.2 Limits/Continuity & 12.3 Derivatives/Integrals notes by Tim Pilachowski r r r =, with a domain of real ( )

f(x) dx with An integral having either an infinite limit of integration or an unbounded integrand is called improper. Here are two examples dx x x 2

Laplace Examples, Inverse, Rational Form

Circuits 24/08/2010. Question. Question. Practice Questions QV CV. Review Formula s RC R R R V IR ... Charging P IV I R ... E Pt.

Understanding small, permanent magnet brushed DC motors.

IX.1.1 The Laplace Transform Definition 700. IX.1.2 Properties 701. IX.1.3 Examples 702. IX.1.4 Solution of IVP for ODEs 704

Fractional Order Thermoelastic Deflection in a Thin Circular Plate

when t = 2 s. Sketch the path for the first 2 seconds of motion and show the velocity and acceleration vectors for t = 2 s.(2/63)

Positive and negative solutions of a boundary value problem for a

4.8 Improper Integrals

Motion on a Curve and Curvature

Degree of Approximation of a Class of Function by (C, 1) (E, q) Means of Fourier Series

Motion. ( (3 dim) ( (1 dim) dt. Equations of Motion (Constant Acceleration) Newton s Law and Weight. Magnitude of the Frictional Force

e t dt e t dt = lim e t dt T (1 e T ) = 1

can be viewed as a generalized product, and one for which the product of f and g. That is, does

EECE 301 Signals & Systems Prof. Mark Fowler

7 Wave Equation in Higher Dimensions

Sample Final Exam (finals03) Covering Chapters 1-9 of Fundamentals of Signals & Systems

On The Estimation of Two Missing Values in Randomized Complete Block Designs

Physics 15 Second Hour Exam

0 for t < 0 1 for t > 0

Science Advertisement Intergovernmental Panel on Climate Change: The Physical Science Basis 2/3/2007 Physics 253

LECTURE 5. is defined by the position vectors r, 1. and. The displacement vector (from P 1 to P 2 ) is defined through r and 1.

Physics 240: Worksheet 16 Name

ÖRNEK 1: THE LINEAR IMPULSE-MOMENTUM RELATION Calculate the linear momentum of a particle of mass m=10 kg which has a. kg m s

Bipartite Matching. Matching. Bipartite Matching. Maxflow Formulation

D zone schemes

Code_Aster. Elements of Fourier for axisymmetric structures

Consider a Binary antipodal system which produces data of δ (t)

ANALYSIS OF KINEMATICS AND KINETOSTATICS OF FOUR-BAR LINKAGE MECHANISM BASED ON GIVEN PROGRAM

More on ODEs by Laplace Transforms October 30, 2017

The Purpose of this talk The generation of the high-frequency resonant FEL wave by means of it s low-frequency wave as a pomp wave

The Production of Well-Being: Conventional Goods, Relational Goods and Status Goods

5.1-The Initial-Value Problems For Ordinary Differential Equations

Reinforcement learning

Network Flows: Introduction & Maximum Flow

u(t) Figure 1. Open loop control system

Primal and Weakly Primal Sub Semi Modules

ENGR 1990 Engineering Mathematics The Integral of a Function as a Function

Physics 120 Spring 2007 Exam #1 April 20, Name

Laplace Transform. Inverse Laplace Transform. e st f(t)dt. (2)

Fig. 1S. The antenna construction: (a) main geometrical parameters, (b) the wire support pillar and (c) the console link between wire and coaxial

To become more mathematically correct, Circuit equations are Algebraic Differential equations. from KVL, KCL from the constitutive relationship

Let. x y. denote a bivariate time series with zero mean.

Contraction Mapping Principle Approach to Differential Equations

Section P.1 Notes Page 1 Section P.1 Precalculus and Trigonometry Review

1. Kinematics of Particles

Chapter 6. Laplace Transforms

EE Control Systems LECTURE 2

U>, and is negative. Electric Potential Energy

CS4445/9544 Analysis of Algorithms II Solution for Assignment 1

EECE 301 Signals & Systems Prof. Mark Fowler

Physics 2A HW #3 Solutions

dt = C exp (3 ln t 4 ). t 4 W = C exp ( ln(4 t) 3) = C(4 t) 3.

Algorithmic Discrete Mathematics 6. Exercise Sheet

Discussion Session 2 Constant Acceleration/Relative Motion Week 03

Physics 201, Lecture 5

Representing Knowledge. CS 188: Artificial Intelligence Fall Properties of BNs. Independence? Reachability (the Bayes Ball) Example

5.2 GRAPHICAL VELOCITY ANALYSIS Polygon Method


One of the common descriptions of curvilinear motion uses path variables, which are measurements made along the tangent t and normal n to the path of

Control Volume Derivation

ANSWERS TO ODD NUMBERED EXERCISES IN CHAPTER

Chapter 2. Kinematics in One Dimension. Kinematics deals with the concepts that are needed to describe motion.

Then the number of elements of S of weight n is exactly the number of compositions of n into k parts.

Computer Propagation Analysis Tools

t s (half of the total time in the air) d?

Lecture 26: Leapers and Creepers

Lecture 22 Electromagnetic Waves

An Open cycle and Closed cycle Gas Turbine Engines. Methods to improve the performance of simple gas turbine plants

ENGI 4430 Advanced Calculus for Engineering Faculty of Engineering and Applied Science Problem Set 9 Solutions [Theorems of Gauss and Stokes]

Math 2142 Exam 1 Review Problems. x 2 + f (0) 3! for the 3rd Taylor polynomial at x = 0. To calculate the various quantities:

Design and Simulation of a Control System For the Bucket Wheel Excavator

Today - Lecture 13. Today s lecture continue with rotations, torque, Note that chapters 11, 12, 13 all involve rotations

ME 391 Mechanical Engineering Analysis

v T Pressure Extra Molecular Stresses Constitutive equations for Stress v t Observation: the stress tensor is symmetric

Lectures # He-like systems. October 31 November 4,6

PHYS PRACTICE EXAM 2

A LOG IS AN EXPONENT.

FM Applications of Integration 1.Centroid of Area

SIMPLIFIED ASSESSING OF SUBSTATION- ORIGINATED OUTAGES IN ANALYSIS OF TRANSMISSION SYSTEM RELIABILITY

Fourier-Bessel Expansions with Arbitrary Radial Boundaries

8.5 Circles and Lengths of Segments

( ) ( ) ( ) () () Signals And Systems Exam#1. 1. Given x(t) and y(t) below: x(t) y(t) (A) Give the expression of x(t) in terms of step functions.

Derivation of the differential equation of motion

Lecture 17: Kinetics of Phase Growth in a Two-component System:

Physics 232 Exam II Mar. 28, 2005

Transcription:

Eindhoven Univeiy of echnology Depen of Bioedicl Engineeing Meil echnology Modeling of he heodiluion cuve fo coony flow een Cel Se U/e 9-6-7 BME7. pof.d.i. F.N. vn de Voe

Conen H Inoducion 3 H Mehod 5. Inoducion 5. he ixing che 5.3 he ue 8.3. Conevion euion 8.3. He lnce fo he wll 9.3.3 heodiluion euion.4 Soluion of he pole.4. Pole definiion.4. Scling.4.3 Soluion 3 H3 eul 5 H4 Concluion nd dicuion 7 efeence 8

H Inoducion he coony ciculion upplie he h wih lood. Evey ie he h e, ou 5 pecen of he oke volue goe o he coony eie. In eing iuion, when he cdic oupu i ou 4.5 l/in, hi yield n vege coony lood flow of 5 l/in []. When coony ey i enoed, he flow hough h ey will decee, depending on he eveiy of he enoi. Since inevenion y Pecune nluinl Coony Angioply PCA o yp ugey in ode o evculie he coony ciculion hould only e pefoed in evee ce, i i ipon o e le o ccuely deeine coony lood flow in ode o dignoe he eveiy of enoi. Mehod o indice he degee of deceed lood flow, e he een of he Coony Flow eeve CF nd he Fcionl Flow eeve FF []. he CF give he io eween flow duing xiu vodilion hypeeic flow nd flow in e eline flow. he FF i he io eween xiu flow in he enoed ey nd he xiu flow in helhy ey. Boh io e ued o deeine he heodynicl e of he ey nd he peiphel vcul ed. o deeine he flow in n ey, he ehod of heodiluion cn e ued. hi ehod ue olu injecion of cold line ino he loode nd eue he epeue chnge locion dil of he injecion. An exple of he kind of epeue pofile which i eing eued, i picued elow. hi cuve, whee he eued epeue i ploed gin he ie, i clled heodiluion cuve figue.. figue.: he heodiluion cuve d fo N.H.J. Pijl, ued wih peiion 3

If he line injecion nd he heodiluion cuve would e iple del-funcion, he flow velociy would e deived y dividing he veled dince of he line y he ie ken o vel h dince. Since he cuve i no del-funcion, en ni ie, n, i inoduced. n i he veged ie i ke fo he cold line o vel fo he poin of injecion o he poin of epeue eueen. n i clculed uing []: n. Hee epeen he epeue dop given. Fo n nd he locion of he epeue eueen he flow velociy cn e deeined. he gol of hi udy i o odel he heodiluion cuve. hi i done in ode o e le o e he lood flow oe pecie nd o coe o n uoic flow deecion lgoih. In chpe he ehod ued o odel he heodiluion cuve e explined. Afe n inoducion, he fi p of he odel i peened. Hee he ixing of he cold line wih he lood in fon of he ennce o he coony ey i epeened y ixing che. he euling ixue hen vel ino he coony ey. he ey i odeled ue in ecion.3 nd he effec on he epeue y diffuion, convecion nd he exchnge hough he wll e clculed uing euion fo conevion of nd enegy. he euling heodiluion euion i hen olved in ecion.4. Chpe 3 conin he eul nd he epo i concluded in chpe 4. 4

Mehod. Inoducion In ode o odel he heodiluion cuve, wo copen e defined [3]. he fi copen i he cul injecion of he cold line ino he lood. he line i injeced hough chee which end ju in fon of he oic vlve nd he ennce o he igh o lef coony eie. Hee he line i ued o e ixed wih he lood efoe hi ixue vel ino he econd copen of he odel, he coony ey. In he coony ey, he ixue of line nd lood vel fo he ennce of he ey, p he epeue euing locion nd fuhe ino he coony nche.. he ixing che he locion whee he cold line i injeced i efeed o he ixing che. Hee he line ixe wih he lood nd hi ixue ene he coony ey. he ixing che i defined depiced in figue.:, lood flow line flow flow o o flow o coony ey, V, lood epeue line epeue ixue epeue V volue of he ixing e, figue.: Definiion of he ixing che Coplince of he ixing che i negleced, o conevion euie:. Fo conn nd eul deniy nd conn nd eul pecific he conn peue fo oh lood nd line, he he euion i: 5

V d. enging nd dividing y V give: d V V.3 o oin, he genel oluion fo fi ode diffeenil euion i ued: if y p x y x, wih h p x dx hen [ ] h h y x e e dx c Fo he ixing che diffeenil euion hi yield, wih : V V e e V.4 If he line injecion i ingle ue wve h egin nd end 4 figue., hn he coeponding epeue pofile, ed on he given pee vlue, i depiced in figue.3. Hee, edy ce i ued. o o figue.: he inpu flow; he injecion i epeened y he olid line. 6

76 l/ 4 l/ 7 l/ l/ V 8 l 37 C C figue.3: he ixing che oupu epeue he oupu epeue will e he inpu epeue fo he econd copen of he odel, which i he coony ey. 7

.3 he ue he epeue of he ixue, while veling hough he coony ey, chnge ove ie due o he convecion, diffuion nd he exchnge he wll. hee fco e ued o clcule he chnge in epeue of he ixue ove cein dince nd ie in n ey. o iplify he pole, he ey i defined ue wih pecificion given in figue.4. u d wih: u velociy in -diecion ue diu d wll hickne figue.4: Definiion of he ue.3. Conevion euion o clcule he chnge in epeue, conevion euion fo enegy nd e needed. he enegy euion wihou inenl he poducion nd vicou diipion uing cylindicl coodine i defined : ρc u θ λ[ ] θ θ p u u.5 wih lood deniy ρ, he cpciy of lood c p nd he diffuion coefficien of lood λ. he conevion euion i defined : ρ ρ v.6 ρ Blood i ued o e incopeile, o nd ρ. hi led o v.7 8

Wien in cylindicl coodine hi give: uθ u v u.8 θ Becue fully developed flow i ued, he velociie h e no in he diecion of u he longiudinl xi of he ue e eo: u u θ, o fo.8:. Suiuing hi ino he enegy euion.5 give: u, fo.9 λ wih:, he he euliing coefficien fo lood. ρ c p.3. He lnce fo he wll When uing.9 o deeine he epeue pofile in he wll of he ue, he he diffuion in -diecion cn e negleced, ecue << in he wll. hi led o he following euion fo he he lnce in he wll: w, fo < d. wih w he he euliing coefficien of he wll If he peneion deph of he he in he wll i uch lle hen he wll hickne, hen he e i eo. Since w i conn hi led o: d d d d. Inegion of. wih d he oue wll diee, whee i he ody epeue, nd he inne wll diee, whee : ln ln d. ln ln d wih: epeue in dil diecion epeue on he ouide of he wll epeue on he inide of he wll d wllhickne king he deivive of. give he e: 9

d ln ln.3 which cn e ued in.9..3.3 heodiluion euion Ineging.9 ove dil lice now give: d d d u d π π π π.4 ln ln d u π π π π.5 divide y π : ln ln d u.6 Now define ε : ln ln ln d c d p ρ λ ε.7 nd ince u i conn nd coepond wih he conn flow defined in he ixing che, A u. Hee A i he co-ecionl e of he ue. he euion fo he epeue chnge, o heodiluion euion, ecoe: o A ε, fo >, >.8

.4 Soluion of he pole.4. Pole definiion In ecion.3.3 he heodiluion euion i defined : A ε, fo >, >.9 o olve hi euion n iniil condiion cn e pecified:,. nd he ixing che oupu epeue, in ecion., i he oundy condiion fo : V V V e e F,. I i convenien o cle.9, in ode o iplify he pole..4. Scling, nd e he vile h need o e de dienionle, ince ll ohe e conn. hey e de dienionle uing:. Suiue hee in.9 nd enge: A ε.3 divide y nd uiue ε : A ε ε.4

Now uiue : εa εa.5 εa And wih D c, he euion ecoe: D c.6 o iplify he ue of he euion, he index i dopped nd he dienionle heodiluion euion ecoe: D c.7 Since he e in.7 e cled nd O, he vlue of D c i evlued uing [3]: ε. / A 7.e-6.4e-7 / 4e-6 3 / D c h vlue of ppoxiely 4.8e-8. hi en i i poile o dicd he econd ode e in.7:.8 In ode o olve.8 he vile nd e de dienionfull gin uing:, wih ε nd, wih εa I i convenien o leve dienionle fo now, u he index i einoduced o e le o dicen eween he dienionle nd dienionfull. Euion.8 now ed fe enging:

A ε.9 o iplify hi, he he euliing coefficien fo he wll, ε, i ued eo. hi leve he line dvecion euion: A.3 Suiuing he dienionle epeue, ccoding o only eul in he lo of he ccen:, ino.3 A.3.4.3 Soluion he oluion of.3 cn e oined wih he ehod of eping vile. wo new vile eplce, ccoding o, X Y. Suiuion in.3 nd enging yield: Y Y AX X.3 Euion.3 only hold if oh ide of he euion e eul o he e conn. In hi ce he conn i defined k nd wo epe euion e oined: Y Y k nd AX X k.33 Soluion of.33 nd einoducion of, led o: ln Y k C nd ln X k C A.34 Y k C e nd X A k C e.35, A k X Y Ce.36 he conn C nd k cn e olved uing he iniil nd oundy condiion fo. nd.. 3

4 Coining. nd.36 led o: k A Ce k A e C.37 Coining. nd.37 led o: k Ce F A k e F ln A F k.38 Now he vlue fo C nd k cn e uiued in.36: ln, A F e.39 o iplified: A A F,.4 wih: V V V e e F.4

3 eul Uing.4, epeue cn e clculed fo ny poiive nd. hi en h fo given dince fo he ennce of he ue, epeue pofile cn e oined. An exple i hown in figue 3., whee 5 c nd he ixing che oupu fo figue.3 i ued inpu. 76 l/ 4 l/ 7 l/ l/ V 8 l 37 C C figue 3.: he heodiluion oupu epeue he pofile h i oined i lo excly he e h of he ixing che oupu. Howeve, he pofile doe flen lighly owd 37 o C fo inceing eueen dince fo he ennce of he ue. Figue 3. how he ppoxiion of he pien d fo figue. y he heodiluion odel. he ppoxiion w de y viion of he pee fo.4. Mo of he pee e known nd hey cn e ued o deeine he wo eining unknown pee. hee unknown pee h need o e djued o oin pope ppoxiion e he coony flow,, nd he ixing che conol volue, V. he pee ued o fi he pien d in figue 3. e: 8 l/ 4 l/ 5 l/ 37 C C 5

5 l/ V 6 l figue 3.: heodiluion cuve i fied o pien d 6

H4 Concluion nd dicuion A odel w de o poduce heodiluion cuve. In he ixing che copen of he odel, he ixing of he injeced cold line wih he lood i odeled. In he econd copen, he npo of he ixue hough ue i odeled, euling in.8. he ue of.4 oundy condiion fo olving.8 finlly euled in.4. hi euion decie he epeue of he conveced line nd lood ixue. o e le o eie coony flow, evel eueen pecific pee ued in.4 nd.4 need o e defined. Uing hee pee, cuve cn e ched o pien d y djuing he unknown coony flow,, nd he ixing che conol volue, V. Howeve, hi ching pocedue need o e fuhe evlued efoe i cn e pplied. Fo he oined heodiluion cuve, i ppe h pien d cn e ppoxied wih he cuen odel. Whehe he ppoxiion yield ule vlue fo he coony flow h no een inveiged in hi udy. Alo, he odel w iplified o line dvecion pole, which poly educe he ccucy of he eiion of he coony flow. Howeve, if he npo of he ixue cn e decied y line dvecion euion, hen he inpu epeue pofile he ennce of he ue i highly defining fo he finl ppoxiion of he heodiluion cuve. hi en h fuue wok involving he heodiluion cuve pehp hould concene on oe dvnced odeling of he line injecion cuve inpu fo he ixing che in ed of he ep funcion ued in hi epo. hi could eul in ee heodiluion cuve ppoxiion nd coneuenly oe ccue coony flow clculion. 7

efeence [] A.C. Guyon, J.E. Hll, exook of Medicl Phyiology, W.B. Sunde Copny, [] B. De Buyne, N.H.J. Pijl, L. Sih, M. Wievegg, G.. Heyndickx, Coony heodiluion o Ae Flow eeve, Ciculion ; 4:3-6 [3] J.A.M. Schuun, D nlyi of heodiluion fo coony flow, Meepo M.6, 8