960 likes | 1.39k Views
醫用流體力學. Arterial Fluid Dynamics 邵耀華 台灣大學應用力學研究所. Physiological Fluid Dynamics. Evolution of Arterial Pressure Away from the heart. Systemic Arteries. Conduct blood flow from Left ventricle (LV) to peripheral organs Aortic valve Aortic arch (180° turn) Geometry changes :Tapering
E N D
醫用流體力學 Arterial Fluid Dynamics 邵耀華 台灣大學應用力學研究所 Physiological Fluid Dynamics
Evolution of Arterial Pressure Away from the heart
Systemic Arteries • Conduct blood flow from Left ventricle (LV) to peripheral organs • Aortic valve Aortic arch (180° turn) • Geometry changes :Tapering • Geometry changes : Branching • Mechanical properties changes
Incremental Elastic Modulus of human arteries change with age
Fluid Mechanics of Elastic Conduct • Mass Conservation • Conservation of momentum • Conservation of energy
Background • Fundamental VariablesPressure、 Flow • Geometrical VariablesSize、 Thickness 、 Length、 Curvature • Mechanical Properties Stiffness 、Visco-Elasticity
Equations of Viscous Pipe Flow • Consider a conduct filled with incompressible fluid of density and pressure p, let u be the only non-zero velocity component
Poiseuille’s Law (1840) • Assume steady flow, u= u(r ) with no body forces, the equation of motion
Laminar Poiseuillean flow • Rate of flow through the tube • Mean velocity of flow • Shear stress at the wall
Laminar Poiseuillean flow • Skin friction • Shear stress in terms of skin friction
Implication of Poiseuille’s Law • Q is proportional to the fourth power of the radius. • Q is directly proportional to the pressure difference. • Q is inversely proportional to the length of the tube. • If the arteries becomes constricted, the bloodpressure requires to supply the blood flow adequately will risesubstantially,leading to the state of hypertension.
Optimum design of Blood Vessel Bifurcation (Poiseuille’s formula) For a given pressure drop, 1% change in vessel radius results in a 4% changes in flow Murray (1926) Rosen (1967) Work done Metabolism Energy loss
Minimum cost function for optimum vessel configuration With respect to radiusa the optimum radius The optimum vessel radius is proportional flow to the 1/3 power, and
Optimum vessel bifurcation that with minimum cost function Minimize P at the bifurcation point B An optimum location B would be for arbitrary movements of B.
Let B displaced along A-B direction first The optimum is obtained when
Again, let B displaced in the C-B direction The optimum is obtained when Similarly, displaced B along D-B direction, we find
Similarly, displaced B along D-B direction, we find The continuity equation gives We find which is often referred to as Murray’s Law 37.5°
Let ao denotes the radius of the aorta, and assume equal bifurcation in all generation If the capillary blood vessel has a radius of 5 um and the radius of the aorta is 1.5 cm. We find n=30. The total number of blood vessel is about 230109. Note: in fact arteries rarely bifurcation symmetrically (a1=a2). For human, only one symmetric bifurcation. For dog, there are none.
Pulsatile Blood Flow • Consider pulsatile flow in a circular vessel, p=p(x, t) and u = u(r, t) • For a sinusoidal flow
Pulsatile Blood Flow(2) • The general solution of the ODE in the form involves Bessel functions of complex arguments U(r=a)=0 (non-slip) U(r=0)=finite
Pulsatile Blood Flow(3) • Introducing Womersley number • As 0, the velocity profile becomes parabolic. • As , viscosity is negligible U(r)=-i P/.
Analysis of Blood Flow using Elastic Theory • From Poiseuille’s Law the flux is proportional to the pressure difference (p1-p2). However, the blood flow in veins are remarkably non-linear. • The flow in elastic conduct gradually attains a maximum value as the pressure difference increases and then on longer increases.
Arterial Flow in Elastic Tube • Axial velocity, v • Lumen area, S
Pressure-Diameter relationship • Let T denotes the tension of the blood vessel per unit thickness, wall thickness h, vessel radius a • Let ro be the radius of zero tension state, the Hooke’s Law gives elastic constant E as
Poiseuille’s flow in elastic tube • Consider steady flow in elastic tube of length L, assume the tube is long and the pressure is function of axial coordinate z, let P1 and P2 denote the inlet and outlet pressure and the external pressure surrounding the tube is P0 • Assume the flow through the tube obey Poiseuille’s law, the flow becomes
Transmission of Pulse wave (Velocity) in elastic tube • Consider inviscid and incompressible fluid flow in elastic tube of lumen area A, • By linearizing the equations
Transmission of Pulse wave (Velocity) in elastic tube • Combining the continuity and momentum equations, • The wave equations Pulse Wave Velocity (PWV)
Analysis of Aortic Diastolic and Systolic Pressure Waveforms • Constitutive relationship between aortic volume and pressure where K is the volume elasticity of the aorta, and V0 is the end-systolic volume. • If the aorta is very soft (K is very small), let I(t) and Q(t) denote the inflow and outflow rates, we have
During diastole, the aortic valve is closed and there is no flow into the aorta. Hence I(t) =0. where is a non-invasive measure aortic volume elasticity. Let Td be the duration of diastolic phase, the aortic pressure (Pd) at the end of this phase or just prior to ejection is given by • The volume elasticity that depicts the exponential drop of aortic pressure is given by
Reynolds #; Strouhal #; Womersley • Reynolds number • Strouhal number • Womersley number
Wave propagation in Blood Vessel • Pulse wave propagation in arteries • A(x, t) depends on transmural pressure, Here c is the wave propagation velocity.
For thin walled elastic tube: • Consider the elasticity of the tube, arterial diameter blood pressure For a thick walled elastic tube:
Resonant vibration of flow in a circular tube • When the tube length is equal to the half wave length • This is called the fundamental frequency of the natural vibration. • Hemodynamics : Effects of Frequency on the Pressure-flow relationship of Arterial tree
Mean Velocity Profile (Dog Aorta)
Effect of sinusoidal pressure wave speed of various frequencies on the instantaneous aortic pressure