Peristaltic biofluids flow through vertical porous human vessels using third-grade non-Newtonian fluids model

In this paper, the heat and flow characteristic of third-grade non-Newtonian biofluids flow through a vertical porous human vessel due to peristaltic wall motion are studied. The third-grade model can describe shear thinning (or shear thickening) and normal stress differences, which is acceptable for biofluids modeling. In order to solve the governing equations, the assumption of long-wavelength approximation is utilized. This hypothesis emphasizes that the wavelength of the peristaltic wall motion is large in comparison with the radius of the human vessel, which is widely acceptable in biological investigations. The analytical perturbation method is employed to solve the governing equations. Consequently, analytical expressions for the velocity profile, shear stress, temperature field, and biofluid flow rate are obtained. In addition, the effects of the governing parameters such as the third-grade non-Newtonian parameter, Grashof Number, Eckert number, and porosity, on the results are examined.     

Steady flow of couple stress fluid through tubes of slowly varying cross-sections--application to blood flows

Steady laminar flow of a non-Newtonian fluid based on couple stress fluid theory, through narrow tubes of varying cross-sections has been studied theoretically. Asymptotic solutions are obtained for the basic equations and the expressions for the velocity field and the wall shear stress are derived for a general cross-section. Computation and discussions are carried out for the geometries which occur in the context of physiological flows or in particular blood flows. The tapered tubes and constricted tubes are of special importance. It is observed that increase in certain parameters results in erratic flow behaviour proximal to the constricted areas which is further enhanced by the increase in the geometric parameters. This elucidates the implications of the flow in the development of vascular lesions.     

Numerical simulation of peristaltic flow of a biorheological fluid with shear-dependent viscosity in a curved channel

Peristaltic motion of a non-Newtonian Carreau fluid is analyzed in a curved channel under the long wavelength and low Reynolds number assumptions, as a simulation of digestive transport. The flow regime is shown to be governed by a dimensionless fourth-order, nonlinear, ordinary differential equation subject to no-slip wall boundary conditions. A well-tested finite difference method based on an iterative scheme is employed for the solution of the boundary value problem. The important phenomena of pumping and trapping associated with the peristaltic motion are investigated for various values of rheological parameters of Carreau fluid and curvature of the channel. An increase in Weissenberg number is found to generate a small eddy in the vicinity of the lower wall of the channel, which is enhanced with further increase in Weissenberg number. For shear-thinning bio-fluids (power-law rheological index, n < 1) greater Weissenberg number displaces the maximum velocity toward the upper wall. For shear-thickening bio-fluids, the velocity amplitude is enhanced markedly with increasing Weissenberg number.     

Peristaltic flow of a couple stress fluid in an asymmetric channel

This paper presents an analysis of the peristaltic flow of a couple stress fluid in an asymmetric channel. The asymmetric nature of the flow is introduced through the peristaltic waves of different amplitudes and phases on the channel walls. Mathematical modelling corresponding to a two-dimensional flow has been carried out. The flow analysis is presented under long wavelength and low Reynolds number approximations. Closed form solutions for the axial velocity, stream function and the axial pressure gradient are given. Numerical computations have been carried out for the pressure rise per wavelength, friction forces and trapping. It is noted that there is a decrease in the pressure when the couple stress fluid parameter increases. The variation of the couple stress fluid parameter with the size of the trapped bolus is also similar to that of pressure. Furthermore, the friction force on the lower channel wall is greater than that on the upper channel wall.     

Numerical investigation of the non-Newtonian blood flow in a bifurcation model with a non-planar branch

The non-Newtonian fluid flow in a bifurcation model with a non-planar daughter branch is investigated by using finite element method to solve the three-dimensional Navier–Stokes equations coupled with a non-Newtonian constitutive model, in which the shear thinning behavior of the blood fluid is incorporated by the Carreau–Yasuda model. The objective of this study is to investigate the influence of the non-Newtonian property of fluid as well as of curvature and out-of-plane geometry in the non-planar daughter vessel on wall shear stress (WSS) and flow phenomena. In the non-planar daughter vessel, the flows are typified by the skewing of the velocity profile towards the outer wall, creating a relatively low WSS at the inner wall. In the downstream of the bifurcation, the velocity profiles are shifted towards the flow divider. The low WSS is found at the inner walls of the curvature and the lateral walls of the bifurcation. Secondary flow patterns that swirl fluid from the inner wall of curvature to the outer wall in the middle of the vessel are also well documented for the curved and bifurcating vessels. The numerical results for the non-Newtonian fluid and the Newtonian fluid with original Reynolds number and the corresponding rescaled Reynolds number are presented. Significant difference between the non-Newtonian flow and the Newtonian flow is revealed; however, reasonable agreement between the non-Newtonian flow and the rescaled Newtonian flow is found. Results of this study support the view that the non-planarity of blood vessels and the non-Newtonian properties of blood are an important factor in hemodynamics and may play a significant role in vascular biology and pathophysiology.     

Comparison of blood rheological models for physiological flow simulation

The present study investigates the flow effects that different blood constitutive equations induce when employed in numerical simulations in the framework of computational hemodynamics. In accord with experimental studies on the rheological behavior of blood, three blood constitutive equations namely the Casson, Power-Law and Quemada models were used for simulating the shear flow behavior of blood. The case studied is the flow in a channel with a moving part of the boundary and was selected because it reproduces the flow phenomena occurring in realistic arterial conditions. Flow simulation for every model is carried out assuming the same flow rate at the inlet of the channel and different Strouhal numbers reflecting different intensities of the boundary movement. Results show that the modeling of blood as non-Newtonian fluid has marked qualitative and quantitative effects on both the flow field and the wall shear stress whereas comparison of the different models shows good agreement between the flow effects by the Casson and Quemada models.     

Numerical 3D-simulation of pulsatile wall shear stress in an arterial T-bifurcation model

The structure of pulsatile blood flow and wall shear stress in a 90° T-bifurcation model is analysed numerically. The nonlinear Navier-Stokes equations for time-dependent incompressible Newtonian fluid flow are approximated using a newly developed pressure correction, finite element method. The wall shear stress is calculated from the finite element velocity field. The investigation shows viscous flow phenomena such as flow separation and stagnation and the distribution of high and low wall shear stress during the pulse cycle. Furthermore, the effect of a sharp corner the bifurcation edge on the wall shear stress is analysed. Detailed local flow investigation is required to examine fluid dynamic contribution to the development of arterial diseases such as atherosclerosis and thrombosis.     

A mathematical simulation of the ureter: effects of the model parameters on ureteral pressure/flow relations

Ureteral peristaltic mechanism facilitates urine transport from the kidney to the bladder. Numerical analysis of the peristaltic flow in the ureter aims to further our understanding of the reflux phenomenon and other ureteral abnormalities. Fluid-structure interaction (FSI) plays an important role in accuracy of this approach and the arbitrary Lagrangian-Eulerian (ALE) formulation is a strong method to analyze the coupled fluid-structure interaction between the compliant wall and the surrounding fluid. This formulation, however, was not used in previous studies of peristalsis in living organisms. In the present investigation, a numerical simulation is introduced and solved through ALE formulation to perform the ureteral flow and stress analysis. The incompressible Navier-Stokes equations are used as the governing equations for the fluid, and a linear elastic model is utilized for the compliant wall. The wall stimulation is modeled by nonlinear contact analysis using a rigid contact surface since an appropriate model for simulation of ureteral peristalsis needs to contain cell-to-cell wall stimulation. In contrast to previous studies, the wall displacements are not predetermined in the presented model of this finite-length compliant tube, neither the peristalsis needs to be periodic. Moreover, the temporal changes of ureteral wall intraluminal shear stress during peristalsis are included in our study. Iterative computing of two-way coupling is used to solve the governing equations. Two phases of nonperistaltic and peristaltic transport of urine in the ureter are discussed. Results are obtained following an analysis of the effects of the ureteral wall compliance, the pressure difference between the ureteral inlet and outlet, the maximum height of the contraction wave, the contraction wave velocity, and the number of contraction waves on the ureteral outlet flow. The results indicate that the proximal part of the ureter is prone to a higher shear stress during peristalsis compared with its middle and distal parts. It is also shown that the peristalsis is more efficient as the maximum height of the contraction wave increases. Finally, it is concluded that improper function of ureteropelvic junction results in the passage of part of urine back flow even in the case of slow start-up of the peristaltic contraction wave.     

Pulsatile non-Newtonian haemodynamics in a 3D bifurcating abdominal aortic aneurysm model

Numerical prediction of non-Newtonian blood flow in a 3D abdominal aortic aneurysm bifurcating model is carried out. The non-Newtonian Carreau model is used to characterise the shear thinning behaviour of the human blood. A physical inlet velocity waveform incorporating a radial velocity distribution reasonably representative of a practical case configuration is employed. Case studies subject to both equal and unequal outlet pressures at iliac bifurcations are presented to display convincingly the downstream pressure influences on the flow behaviour within the aneurysm. Simulations indicate that the non-Newtonian aspects of the blood cannot at all be neglected or given a cursory treatment. The wall shear stress (WSS) is found to change significantly at both the proximal and distal ends of the aneurysm. At the peak systole, the WSS is peak around the bifurcation point, whereas the WSS becomes zero in the bifurcation point. Differential downstream pressure fields display significant effects regarding the flow evolution in the iliac arteries, whereas little or no effects are observed directly on the flow details in the aneurysm.     

Modeling pulsatile flow in aortic aneurysms: effect of non-Newtonian properties of blood

Pulsatile flow in an axisymmetric rigid-walled model of an abdominal aorta aneurysm was analyzed numerically for various aneurysm dilations using physiologically realistic resting waveform at time-averaged Reynolds number of 300 and peak Reynolds number of 1607. Discretization of the governing equations was achieved using a finite element scheme based on the Galerkin method of weighted residuals. Comparisons with previously published work on the basis of special cases were performed and found to be in excellent agreement. Our findings indicate that the velocity fields are significantly affected by non-Newtonian properties in pathologically altered configurations. Non-Newtonian fluid shear stress is found to be greater than Newtonian fluid shear stress during peak systole. Further, the maximum shear stress is found to occur near the distal end of AAA during peak systole. The impact of non-Newtonian blood flow characteristics on pressure compared to Newtonian model is found insignificant under resting conditions. Viscous and inertial forces associated with blood flow are responsible for the changes in the wall that result in thrombus deposition and dilation while rupture of AAA is more likely determined by much larger mechanical stresses imposed by pulsatile pressure on the wall of AAA.     

Dynamic response of wall shear stress on the stenosed artery

The present study deals with an appropriate mathematical model of an artery in the presence of constriction in which the generated wall shear stress due to blood flow is analysed. The geometry of the stenosed arterial segment in the diseased state, causing malfunction of the cardiovascular system, is formed mathematically. The flowing blood contained in the stenosed artery is treated as non-Newtonian and the flow is considered to be two-dimensional. The motion of the arterial wall and its effect on local fluid mechanics is not ruled out from the present pursuit. The flow analysis applies the time-dependent, two-dimensional incompressible nonlinear Navier–Stokes equations for non-Newtonian fluids. The flow-field can be obtained primarily following the radial coordinate transformation, using the appropriate boundary conditions and finally adopting a suitable finite difference scheme numerically. The influences of flow unsteadiness, the arterial wall distensibility and the presence of stenosis on the flow-field and the wall shear stresses are quantified in order to indicate the susceptibility to atherosclerotic lesions and thereby to validate the applicability of the present theoretical model.     

Effects of the non-Newtonian viscosity of blood on flows in a diseased arterial vessel. Part 1: Steady flows

Effects of the non-Newtonian viscosity of blood on a flow in a coronary arterial casting of man were studied numerically using a finite element method. Various constitutive models were examined to model the non-Newtonian viscosity of blood and their model constants were summarized. A method to incorporate the non-Newtonian viscosity of blood was introduced so that the viscosity could be calculated locally. The pressure drop, wall shear stress and velocity profiles for the case of blood viscosity were compared for the case of Newtonian viscosity (0.0345 poise). The effect of the non-Newtonian viscosity of blood on the overall pressure drop across the arterial casting was found to be significant at a flow of the Reynolds number of 100 or less. Also in the region of flow separation or recirculation, the non-Newtonian viscosity of blood yields larger wall shear stress than the Newtonian case. The origin of the non-Newtonian viscosity of blood was discussed in relation to the viscoelasticity and yield stress of blood.     

The influence of the non-Newtonian properties of blood on the flow in large arteries: steady flow in a carotid bifurcation model.

Laser Doppler anemometry experiments and finite element simulations of steady flow in a three dimensional model of the carotid bifurcation were performed to investigate the influence of non-Newtonian properties of blood on the velocity distribution. The axial velocity distribution was measured for two fluids: a non-Newtonian blood analog fluid and a Newtonian reference fluid. Striking differences between the measured flow fields were found. The axial velocity field of the non-Newtonian fluid was flattened, had lower velocity gradients at the divider wall, and higher velocity gradients at the non-divider wall. The flow separation, as found with the Newtonian fluid, was absent. In the computations, the shear thinning behavior of the analog blood fluid was incorporated through the Carreau-Yasuda model. The viscoelastic properties of the fluid were not included. A comparison between the experimental and numerical results showed good agreement, both for the Newtonian and the non-Newtonian fluid. Since only shear thinning was included, this seems to be the dominant non-Newtonian property of the blood analog fluid under steady flow conditions.     

Pulsatile non-Newtonian flow characteristics in a three-dimensional human carotid bifurcation model.

Numerical analysis of flow phenomena and wall shear stresses in the human carotid artery bifurcation has been carried out using a three-dimensional geometrical model. The primary aim of this study is the detailed discussion of non-Newtonian flow velocity and wall shear stress during the pulse cycle. A comparison of non-Newtonian and Newtonian results is also presented. The applied non-Newtonian behavior of blood is based on measured dynamic viscosity. In the foreground of discussion are the flow characteristics in the carotid sinus. The investigation shows complex flow patterns especially in the carotid sinus where flow separation occurs at the outer wall throughout the systolic deceleration phase. The changing sign of the velocity near the outer sinus wall results in oscillating shear stress during the pulse cycle. At the outer wall of the sinus at maximum diameter level the shear stress ranges from -1.92 N/m2 to 1.22 N/m2 with a time-averaged value of 0.04 N/m2. At the inner wall of the sinus at maximum diameter level the shear stress range is from 1.16 N/m2 to 4.18 N/m2 with a mean of 1.97 N/m2. The comparison of non-Newtonian and Newtonian results indicates unchanged flow phenomena and rather minor differences in the basic flow characteristics.     

Wall shear stress distribution in the human carotid siphon during pulsatile flow

Wall shear stress distribution in the carotid siphon, which is a multiple curved segment of the internal carotid artery, is investigated numerically under physiological flow conditions. The computer simulation of flow through the model segment is based on the time-dependent, three-dimensional Navier-Stokes equations, solved numerically with a finite element method. The study shows the behavior of the wall shear stress-vector field and identifies the zones of high and low wall shear stress values during the cardiac cycle.     

Wall shear stress in backward-facing step flow of a red blood cell suspension.

An experimental investigation of the wall shear stress distribution downstream of a backward-facing step is carried out. The wall shear stress distribution was determined by measuring the deformation of a gel layer, attached to the wall downstream of the step. Speckle pattern interferometry was applied to measure the deformation of the gel layer. The measured deformation, combined with the properties of the gel layer, served as an input for a finite element solid mechanics computation to determine the stress distribution in the gel layer. The wall shear stress, required to generate the measured deformation of the gel layer, was determined from these computations. A Newtonian buffer solution and a non-Newtonian red blood cell suspension were used as measuring fluids. The deformation of the gel layer was determined for a Newtonian buffer solution to evaluate the method and to obtain the properties of the gel layer. Subsequently, the wall shear stress distribution for the non-Newtonian red blood cell suspension was determined for three different flow rates. The inelastic non-Newtonian Carreau-Yasuda model served as constitutive model for the red blood cell suspension. Using this model, the velocity and wall shear stress distribution were computed by means of a finite element fluid mechanics computation. From the comparison between the numerical and the experimental results, it can be concluded that wall shear stresses, induced by the red blood cell suspension, can be modeled accurately by employing a Carreau-Yasuda model.     

Mathematical modelling of ciliary propulsion of an electrically-conducting Johnson-Segalman physiological fluid in a channel with slip

Bionic systems frequently feature electromagnetic pumping and offer significant advantages over conventional designs via intelligent bio-inspired properties. Complex wall features observed in nature also provide efficient mechanisms which can be utilized in biomimetic designs. The characteristics of biological fluids are frequently non-Newtonian in nature. In many natural systems super-hydrophobic slip is witnessed. Motivated by these phenomena, in this paper, we discussed a mathematical model for the cilia-generated propulsion of an electrically-conducting viscoelastic physiological fluid in a ciliated channel under the action of magnetic field. The rheological behavior of the fluid is simulated with the Johnson-Segalman constitutive model which allows internal wall slip. The regular or coordinated movement of the ciliated edges (which line the internal walls of the channel) is represented by a metachronal wave motion in the horizontal direction which generates a two-dimensional velocity profile. This mechanism is imposed by a periodic boundary condition which generates propulsion in the channel flow. Under the classical lubrication approximation, the boundary value problem is non-dimensionalized and solved analytically with a perturbation technique. The influence of the geometric, rheological (slip and Weissenberg number) and magnetic parameters on velocity, pressure gradient and the pressure rise (evaluated via the stream function in symbolic software) are presented graphically and interpreted at length.     

A model for blood flow through a stenotic tube

This article deals with the introduction of the modified Casson's fluid model as the true representation for the blood for the steady laminar flow through a small diameter artery with axi-symmetric identical double stenoses in series. The governing equations are solved by using the finite element method. The results for the velocity profiles, the pressure and the wall shear stress distributions in addition to the location and length of the flow reversal zones have been brought out and discussed in reference to the severity of the disease. It has been observed that the non-Newtonian nature of the blood helps in reducing the magnitude of the peak wall shear stress at the throat and the length of the reversed flow regions in the post stenotic dilatation.     

Non-Newtonian Blood Flow in Left Coronary Arteries with Varying Stenosis: A Comparative Study

This paper presents Computational fluid dynamic (CFD) analysis of blood flow in three different 3-D models of left coronary artery (LCA). A comparative study of flow parameters (pressure distribution, velocity distribution and wall shear stress) in each of the models is done for a non-Newtonian (Carreau) as well as the Newtonian nature of blood viscosity over a complete cardiac cycle. The difference between these two types of behavior of blood is studied for both transient and steady states of flow. Additionally, flow parameters are compared for steady and transient boundary conditions considering blood as non-Newtonian fluid. The study shows that the highest wall shear stress (WSS), velocity and pressure are found in artery having stenosis in all the three branches of LCA. The use of Newtonian blood model is a good approximation for steady as well as transient blood flow boundary conditions if shear rate is above 100 s-1. However, the assumption of steady blood flow results in underestimating the values of flow parameters such as wall shear stress, pressure and velocity.