Concluding Remarks
Part I chapters give the necessary knowledge associated with cell behavior, especially when the cell reacts to various stimuli, particularly mechanical stresses. The magnitude and direction of these mechanical stresses applied by the flowing blood on the wetted surface of the endothelium as well as within the vessel wall varies during the cardiac cycle. The heart generates unsteady flows in a network of blood vessels characterized by complicated architecture and variable geometry both in space and time. Vessel geome- try varies over short distances. The vascular network of curved blood ves- sels is composed of successive geometrical singularities, mainly branchings. Three-dimensional blood flow is thereby developing, being conveyed in en- trance length (App. C.7). Moreover, blood vessels are deformable. Changes in transmural pressure (the pressure difference between the pressure at the wet- ted surface of the lumen applied by the moving blood on the vessel wall and the pressure at the external wall side, which depends on the activity on the neighbor organs) can also affect the shape of the vessel cross-section, especially when it becomes negative. More generally, the change in cross-section shape is due to taper, possible prints of adjacent organs with more or less progressive constriction and enlargment, and adaptation to branching (transition zone). These changes can induce three-dimensional blood motions displayed by vir- tual transverse currents, even if the vessel is considered straight (Part II). Local changes in the direction of stress components are caused by flow sep- aration and flow reversal during the cardiac cycle. Flow separation is set by an adverse pressure gradient when inertia forces and blood vorticity are high enough, especially in branching vascular segments. Due to the time-dependent feature of blood flow, the flow separation region spreads over a variable length during the cardiac cycle and can move. The location and variable size of the flow separation region depends on the flow distribution between the branches, which can vary during the cardiac cycle. Flow reversal occurs during the dias- tole of the left ventricle in elastic arteries, such as the aorta, and most of the muscular arteries, such as brachial and femoral arteries (but not in the carotid arteries). Flow reversal can be observed either in a region near the wall, more 490 Concluding Remarks or less wide with respect to the position of the local center of vessel curvature, or in the entire lumen. This book also provides a survey of flow modeling and simulations in the cardiovascular system. The basic knowledge of fluid mechanics necessary to understand blood flow behavior is given in part II. However, the Navier-Stokes equations, the set of partial differential equations (PDE) developed in the framework of continuum mechanics that governs the fluid flow is introduced in the following. Like any mathematical model, Navier-Stokes equations link the physical quantities that are involved in the studied process (both independent and dependent variables, time, space coordinates, temperature, pressure, ve- locity, energy, cross-section area, etc.) using equations which describe the sys- tem behavior. The mathematical model has a value only when realizable mea- surements validate the model, i.e., the model generates verifiable predictions. The Navier-Stokes equations correspond to the mathematical formulation of basic physical principles, especially the conservation of mass and momen- tum. Any fluid motion is indeed characterized by three principles: (1) the matter is neither lost nor created; (2) the variation rate of momentum of any fluid part is equal to the resulting forces (sum of all applied forces) on this fluid part; and (3) the energy is neither lost nor created.
Conservation Equations
Mass and momentum conservation equations are obtained from the analysis of the evolution of involved quantities in an infinitesimal control volume, the so-called fluid particle. The following differential operators are used: (1) the ∇ gradient operator =(∂/∂x1,∂/∂x2,∂/∂x 3), (2) the divergence operator ∇· ∇2 3 2 2 ∇ , and (3) the Laplace operator = i=1 ∂ /∂xi . The gradient p of scalar p is a vector of component (∇p)i = ∂ip. The gradient ∇v of vector v is a second order tensor with component (∇v)ij = ∂ivj. The first basic postulate states that matter is neither destroyed nor cre- ated, meaning that any matter finite volume neither disappers nor becomes infinite. The second basic postulate states that matter is impenetrable; any matter element does not share the same location with another element (con- tinuity axiom). ∂ The elements in presence are: (1) the mass change rate ρdV,and ∂t Ve (2) the mass flux across the surface Γe of the infinitesimal control volume 1 Ve, of unit normal n outward oriented − ρv · n dA. Using the divergence Γe theorem (Gauss formula),2 the equation of mass conservation, or the so-called equation of continuity, is obtained: 1 · = 3 The scalar product of two vectors provides a scalar: v n i=1 vini. 2 ∇·(ρv) dV = (ρv) · n dA. V Γ Concluding Remarks 491
∂tρ + ∇·(ρv)=0. (10.2)
The equation of momentum conservation of a moving fluid of mass density ρ, of dynamic viscosity µ, of kinematic viscosity ν = µ/ρ, which is conveyed with a velocity v(x,t) (x: Eulerian position, t: time) is obtained from the balance of surface and remote forces applied to the fluid particle. The elements in presence are: (1) the body forces f dV , (2) the pressure Ve forces − pn dA, (3) the viscous forces Cn dA, and (4) the inertia forces Γe Γe Dv Dv 3 4 ρ dV ,where = ∂tv +(v ·∇)v. Using the Gauss formulas, the Ve Dt Dt momentum conservation is obtained:
ρ(∂tv +(v ·∇)v)=f −∇p + ∇·C. (10.3)
The Navier–Stokes equations governing an unsteady flow of an incompress- ible (ρ constant) fluid in the absence of body forces are then given by: