RTD and FLOWs at steady and pulsatile flow in systemic circulation

RTD and FLOWs at steady and pulsatile flow in systemic circulation

Motivation•Dispersion of matter in systemic circulation. •Stimulus response method for system identification.•Flow distribution and flowrate evaluated from transit time distributions.•What is the effect of flow pulsation upon flowrate distribution in a branched system

Task prepared within the project FRVS 90/2010

Residence times

INLET OUTLET

Fast particle (short residence time)

Residence time = time spend by a particle inside a continuous system between inlet and outlet.

RTD Residence time distribution E() = histogram of residence times of all particles passing through outlet.

Applications of E(): Prediction of chemical reactions yield, transport of farmaceuticals, evaluation of active and stagnant volumes, diagnostics,… There are generally different RTDs for plasma and RTDs of blood particles. Different RTD exists at steady and pulsatile flow.

Stimulus response experiment

INLET OUTLET

All particles passing through inlet at time interval (tI,tI+dt)

are marked “red”

Unit IMPULSE stimulus

E()=impulse response

Experiments: “Red particles” – tracers (dyes, radionuclides, salts,…) are quickly injected to inlet. Detectors monitor concentration of tracer at outlet (impulse response).

In our case the tracer is KCl solution and reflective particles injected from syringe, detectors are electric conductivity probes.

Response by Convolution

E

cin

cout

t [s]

0

)()()( dtEctc inout

)(~

)(~)(~ fEfcfc inout Fourier transform of convolution

Convolution of stimulus and RTD function

cin(t), cout(t) concentration of tracer at inlet and outlet of a continuous system

Convolution enables to calculate response to an arbitrary stimulus function. It is possible to calculate impulse response E(t) by deconvolution from measured stimulus and response.

Mean time of response = mean time of stimulus + mean residence time

Variance of response = variance of stimulus + variance of E

RTD of pipe at convective regimeLaminar flow and short residence times: diffusion can be neglected and residence times can be calculated from velocity field. The diffusion free regime is characterized by

where Q-flow rate, L-pipe length, Dm-diffusion coefficient.

Example: C~1000 at aorta, arteries, vena cava.

1LD

QC

m

0

0.1

0.2

0.3

0.4

0.5

0 5 10 15t [s]

C,E

cin

cout

E()s8

sout 9

(s)First appearance time =4s

3

2

2)(

E

2

Analytical solution exists also for response to a pulse of a finite width (violet line). Time of the first appearance is not affected, only mean response time is increased.

Residence time distribution E() is the response to infinitely short impulse (delta function)

Laboratory model of circulation

Simulated part - pressure or conductivity transducers alternatively

Cresto Platinum wires

NODES of simulation model

Conductivity probes. Q-flowrate, C-response (defined in form of table)2

34

56

7

8

9

1011

12

13

14

1516

17

1819

20

21

22

23

2425

26

27

28

29

Data file defines the coordinates x,y,z of all 29 nodes

X [m]

Y [m]

1 -0.4 -0.3 -0.2 -0.1 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8

ELEMENTS of simulation model

P1(1,2,d1)

25

Data files define connectivity of

PIPES (2-node elements)

DIVISION (3-node elements)

P2(3,5,d2)P3(4,8,d2)

P4(8,11,d3)

P5(9,14,d3)P6(7,17,d3)

P7(6,20,d3)

V1(2,3,4,23,d1,d2,d2)

V2(5,6,7,24,d2,d3,d3)V3(8,9,10,25,d2,d3,d3)

V4(11,12,13,26,d3,d4,d4)

V5(14,16,15,27,d3,d4,d4) V6(17,18,19,28,d3,d4,d4)

V7(20,22,21,29,d3,d4,d4)

Nodal indices

RESPONSEs generated by PIPES

0

12 )()()( dtEctcx1,y1,z1 ,Q1, c1(t)

x2,y2,z2 Q2 ,c2(t)

d

3

21

2

2

))4/(()(

QLd

E

Convolution realized by FFT, and by inverse transformation.

RESPONSEs generated by WYES

0

114 )()()( dtEctc

0

242 )()()( dtEctc

0

343 )()()( dtEctcx1,y1,z1, Q1,c1(t)

x2,y2,z2 Q2,c2(t)

d1

x3,y3,z3 Q3,c3(t)

x4,y4,z4

Given c1(t)

Intermediate response

0

12 )()()( dtEctc mixed

final response

It is alternatively possible to characterize the whole wye as a perfectly mixed vessel. In this case the both responses c3 and c2 will be identical

PROBLEMS:

1. Derive impulse response of perfectly mixed vessel Emixed(t)

2. How to modify the responses in the case that detectors record not the tracer concentration but the flowrate of tracer (e.g. radiotracer – counter records rate of decays in the whole outlet).

Simulation in MATLABMATLAB M-files available at http:

Simulation proceeds according to M-file RTD.m in steps1. Reading input files xyzq,seq,cp,cv (txt) defining nodes, flowrate and elements2. Definition of time scale, time steps, and stimulus function at entry node 1.3. Evaluation of responses sequentially in all elements (sequence is defined in

vector seq). Responses are evaluated by FFT in functions

function cout=convpipe(t,n,dt,cin,d,l,q)

tmean=pi*d^2*l/(4*q);for i=1:n if t(i)<tmean/2

e(i)=0; else

e(i)=tmean^2/(2*t(i)^3); endendfcin=fft(cin);fe=fft(e);fcout=fcin.*fe;cout=ifft(fcout)*dt;end

Simulation in MATLABMATLAB M-files available at http:

Input data files:xyzq.txt x y z q (nodal coordinates and flowrates)

seq.txt e1 e2 … (sequence of evaluated elements, +pipes, -wyes)

cp.txt i j d (pipe -indices of nodes and diameter)

cv.txt i j k l d1 d2 d3 (wye - indices of nodes and diameters)

TASK:

1. Modify input files according to parameters of your experiment

2. Compare calculated and measured concentration responses

3. Evaluate flowrates in branches from the shortest/mean transit times

Experiments: UVP velocitiesUltrasound Doppler effect for measurement velocity profiles

1. Piezotransducer is transmitter as well as receiver of US pressure waves operating at frequency 4 or 8 MHz.

2. Short pulse of US waves is transmitted (repetition frequency 244Hz and more) and crystal starts listening received frequency reflected from particles in fluid.

3. Time delay of sampling (flight time) is directly proportional to the distance between transducer and the reflecting particle moving with the same velocity as liquid.

4. Received frequency differs from the transmitted frequency by Doppler shift Δf, that is proportional to the component of particle velocity in the direction of transducer axis.

PROBLEMS:

1. What is spatial resolution of velocity, knowing speed of sound in water (1400 m/s) and sampling frequency 8 MHz ?

2. Calculate flowrate in a circular pipe from recorded velocity profile (given angle )

http://biomechanika.cz

Experiments: flowrates

http://biomechanika.cz

Flowrate can be identified also from the stimulus response experiments

1. Mean residence time

2. Shortest residence time (time of the first appearance)

3. Cross-correlation of stimulus and response functions

inout ttV

Q

,

LABORATORY REPORT1. Front page: Title, authors, date

2. Content, list of symbols

3. Introduction, aims of project, references

4. Description of experimental setup

5. Experiments with injection of tracer at a steady state regime. Record responses and pressures.

6. Evaluate flowrate from recorded pressure drop. Check Re and flow regime.

7. Evaluate flowrate from recorded responses and transit times.

8. Evaluate flowrate from velocity profiles recorded by UVP monitor.

9. Simulation (steady state). Modify input files and M-files according to measured flowrates and concentration profile at inlet. Evaluate theoretical responses at nodes with conductivity detectors.

10. Adjust flowrate pulsation and repeat tracer experiments.

11.Evaluate flowrate distribution at pulsatile flow and compare relative flowrates with steady state case.

12.Conclusion (identify interesting results and problems encountered)