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.     

A Mathematical Study on Three Layered Oscillatory Blood Flow Through Stenosed Arteries

A mathematical model is constructed to examine the characteristics of three layered blood flow through the oscillatory cylindrical tube (stenosed arteries).The proposed model basically consists three layers of blood (viscous fluids with different viscosities) named as core layer (red blood cells),intermediate layer (platelets/white blood cells) and peripheral layer (plasma).The analysis was restricted to propagation of small-amplitude harmonic waves,generated due to blood flow whose wave length is larger compared to the radius of the arterial segment.The impacts of viscosity of fluid in peripheral layer and intermediate layer on the interfaces,average flow rate,mechanical efficiency,trapping and reflux are discussed with the help of numerical and computational results.This model is the generalized form of the preceding models.On the basis of present discussion,it is found that the size of intermediate and peripheral layers reduces in expanded region and enhances in contracted region with the increasing viscosity of fluid in peripheral layer,whereas,opposite effect is observed for viscosity of fluid in intermediate layer.Final conclusion is that the average flow rate and mechanical efficiency increase with the increasing viscosity of fluid in both layers,however,the effects of the viscosity of fluid in both layers on trapping and reflux are opposite to each other.     

Heat Transfer Analysis for Stationary Boundary Layer Slip Flow of a Power-Law Fluid in a Darcy Porous Medium with Plate Suction/Injection

In this paper, we investigate the slip effects on the boundary layer flow and heat transfer characteristics of a power-law fluid past a porous flat plate embedded in the Darcy type porous medium. The nonlinear coupled system of partial differential equations governing the flow and heat transfer of a power-law fluid is transformed into a system of nonlinear coupled ordinary differential equations by applying a suitable similarity transformation. The resulting system of ordinary differential equations is solved numerically using Matlab bvp4c solver. Numerical results are presented in the form of graphs and the effects of the power-law index, velocity and thermal slip parameters, permeability parameter, suction/injection parameter on the velocity and temperature profiles are examined.     

Double-diffusion convective biomimetic flow of nanofluid in a complex divergent porous wavy medium under magnetic effects

We explore the physical influence of magnetic field on double-diffusive convection in complex biomimetic (peristaltic) propulsion of nanofluid through a two-dimensional divergent channel. Additionally, porosity effects along with rheological properties of the fluid are also retained in the analysis. The mathematical model is developed by equations of continuity, momentum, energy, and mass concentration. First, scaling analysis is introduced to simplify the rheological equations in the wave frame of reference and then get the final form of equations after applying the low Reynolds number and lubrication approach. The obtained equations are solved analytically by using integration method. Physical interpretation of velocity, pressure gradient, pumping phenomena, trapping phenomena, heat, and mass transfer mechanisms are discussed in detail under magnetic and porous environment. The magnitude of velocity profile is reduced by increasing Grashof parameter. The bolus circulations disappeared from trapping phenomena for larger strength of magnetic and porosity medium. The magnitude of temperature profile and mass concentration are increasing by enhancing the Brownian motion parameter. This study can be productive in manufacturing non-uniform and divergent shapes of micro-lab-chip devices for thermal engineering, industrial, and medical technologies.

     

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.     

Peristaltic transport of blood: Casson model--II

The problem of peristaltic transport of blood in a uniform and non-uniform tube has been investigated, under zero Reynolds number and long wavelength approximation. Blood is represented by a two-layered fluid model consisting of a central layer of suspension of all erythrocytes, etc., assumed to be a Casson fluid, and a peripheral layer of plasma as a Newtonian fluid. A comparison of results with those without peripheral layer shows that the magnitude of the pressure rise, under a given set of conditions is smaller in the case of model with peripheral layer. It is found that, for a given flow rate, the pressure rise decreases as the viscosity of the peripheral layer decreases, and for a given non zero pressure drop, the flow rate increases as the viscosity of the peripheral layer decreases. However, the flow is independent of the presence of the peripheral layer, for zero pressure rise. Further, the pressure rise in the case of non-uniform geometry is found much smaller than the corresponding value in the uniform geometry.     

Mass transport in a microchannel bioreactor with a porous wall

A two-dimensional flow model has been developed to simulate mass transport in a microchannel bioreactor with a porous wall. A two-domain approach, based on the finite volume method, was implemented. For the fluid part, the governing equation used was the Navier-Stokes equation; for the porous medium region, the generalized Darcy-Brinkman-Forchheimer extended model was used. For the porous-fluid interface, a stress jump condition was enforced with a continuity of normal stress, and the mass interfacial conditions were continuities of mass and mass flux. Two parameters were defined to characterize the mass transports in the fluid and porous regions. The porous Damkohler number is the ratio of consumption to diffusion of the substrates in the porous medium. The fluid Damkohler number is the ratio of the substrate consumption in the porous medium to the substrate convection in the fluid region. The concentration results were found to be well correlated by the use of a reaction-convection distance parameter, which incorporated the effects of axial distance, substrate consumption, and convection. The reactor efficiency reduced with reaction-convection distance parameter because of reduced reaction (or flux), and smaller local effectiveness factor due to the lower concentration in Michaelis-Menten type reactions. The reactor was more effective, and hence, more efficient with the smaller porous Damkohler number. The generalized results could find applications for the design of bioreactors with a porous wall.     

In vivo mimicking model for solid tumor towards hydromechanics of tissue deformation and creation of necrosis

The present work addresses transvascular and interstitial fluid transport inside a solid tumor surrounded by normal tissue (close to an in vivo mimicking setup). In general, biological tissues behave like a soft porous material and show mechanical behavior towards the fluid motion through the interstitial space. In general, forces like viscous drag that are associated with the fluid flow may compress the tissue material. On the macroscopic level, we try to model the motion of fluids and macromolecules through the interstitial space of solid tumor and the normal tissue layer. The transvascular fluid transport is assumed to be governed by modified Starling’s law. The poroelastohydrodynamics (interstitial hydrodynamics and the deformation of tissue material) inside the tumor and normal tissue regions is modeled using linearized biphasic mixture theory. Correspondingly, the velocity distribution of fluid is coupled to the displacement field of the solid phase (mainly cellular phase and extracellular matrix) in both the normal and tumor tissue regions. The corresponding velocity field is used within the transport reaction equation for fluids and macromolecules through interstitial space to get the overall solute (e.g., nutrients, drug, and other macromolecules) distribution. This study justifies that the presence of the normal tissue layer plays a significant role in delaying/assisting necrosis inside the tumor tissue. It is observed that the exchange process of fluids and macromolecules across the interface of the tumor and normal tissue affects the effectiveness factor corresponding to the tumor tissue.     

Peristaltic Transport of Carreau-Yasuda Fluid in a Curved Channel with Slip Effects

The wide occurrence of peristaltic pumping should not be surprising at all since it results physiologically from neuro-muscular properties of any tubular smooth muscle. Of special concern here is to predict the rheological effects on the peristaltic motion in a curved channel. Attention is focused to develop and simulate a nonlinear mathematical model for Carreau-Yasuda fluid. The progressive wave front of peristaltic flow is taken sinusoidal (expansion/contraction type). The governing problem is challenge since it has nonlinear differential equation and nonlinear boundary conditions even in the long wavelength and low Reynolds number regime. Numerical solutions for various flow quantities of interest are presented. Comparison for different flow situations is also made. Results of physical quantities are interpreted with particular emphasis to rheological characteristics.     

Analysis of a self-propelling sheet with heat transfer through non-isothermal fluid in an inclined human cervical canal

The present theoretical analysis deals with biomechanics of the self-propulsion of a swimming sheet with heat transfer through non-isothermal fluid filling an inclined human cervical canal. Partial differential equations arising from the mathematical modeling of the proposed model are solved analytically. Flow variables like pressure gradient, propulsive velocity, fluid velocity, time mean flow rate, fluid temperature, and heat-transfer coefficients are analyzed for the pertinent parameters. Striking features of the pumping characteristics are explored. Propulsive velocity of the swimming sheet becomes faster for lower Froude number, higher Reynolds number, and for a vertical channel. Temperature and peak value of the heat-transfer coefficients below the swimming sheet showed an increase by the increment of Brinkmann number, inclination, pressure difference over wavelength, and Reynolds number whereas these quantities decrease with increasing Froude number. Aforesaid parameters have shown opposite effects on the peak value of the heat-transfer coefficients below and above the swimming sheet. Relevance of the current results to the spermatozoa transport with heat transfer through non-isothermal cervical mucus filling an inclined human cervical canal is also explored.     

Mathematica numerical simulation of peristaltic biophysical transport of a fractional viscoelastic fluid through an inclined cylindrical tube

This paper studies the peristaltic transport of a viscoelastic fluid (with the fractional second-grade model) through an inclined cylindrical tube. The wall of the tube is modelled as a sinusoidal wave. The flow analysis is presented under the assumptions of long wave length and low Reynolds number. Caputo's definition of fractional derivative is used to formulate the fractional differentiation. Analytical solutions are developed for the normalized momentum equations. Expressions are also derived for the pressure, frictional force, and the relationship between the flow rate and pressure gradient. Mathematica numerical computations are then performed. The results are plotted and analysed for different values of fractional parameter, material constant, inclination angle, Reynolds number, Froude number and peristaltic wave amplitude. It is found that fractional parameter and Froude number resist the flow pattern while material constant, Reynolds number, inclination of angle and amplitude aid the peristaltic flow. Furthermore, frictional force and pressure demonstrate the opposite behaviour under the influence of the relevant parameters emerging in the equations of motion. The study has applications in uretral biophysics, and also potential use in peristaltic pumping of petroleum viscoelastic bio-surfactants in chemical engineering and astronautical applications involving conveyance of non-Newtonian fluids (e.g. lubricants) against gravity and in conduits with deformable walls.     

Peristaltic transport of multilayered power-law fluids with distinct viscosities: a mathematical model for intestinal flows

This paper is concerned with the theoretical study of two-dimensional peristaltic flow of power-law fluids in three layers with different viscosities. The analysis is carried out under low Reynolds number and long wavelength approximations. The shapes of the interfaces are described by a system of non-linear algebraic equations which are solved numerically as streamlines. It is found that the non-uniformity in the intermediate and peripheral layers diminishes when the viscosity of the intermediate layer is increased and that of the outermost layer is kept unaltered for both the pseudo-plastic and dilatant fluids. Similar are the observations when the viscosity of the outermost layer is increased and that of the intermediate layer is kept fixed. The flow rate increases with the viscosities of the peripheral and the intermediate layers but the viscosity of the outermost layer is more effective. However, the knowledge of the effect of the viscosity of the intermediate layer facilitates us to achieve the required flow rate without disturbing the outermost layer. An increase in the flow behaviour index too favours larger flow rates. The trapping limits increase with the viscosity of the intermediate layer but decrease with the viscosity of the outermost layer and the flow behaviour index. Thus, a medicinal intervention that creates a more viscous intermediate layer and reduces pseudo plasticity may reduce constipation.     

Peristaltic motion of Phan–Thien–Tanner fluid in the presence of slip condition

This investigation considers the peristaltic flow of a Phan–Thien–Tanner fluid in the presence of slip condition and induced magnetic field. By use of the long wavelength and low Reynolds number approximations, closed form series solutions for stream function, pressure gradient, magnetic force function, axial induced magnetic field, and current density were obtained. The pressure gradient and frictional forces per wavelength were computed by numerical integration. The velocity slip condition in terms of shear stress is taken into account. Graphical results show the comparison between no-slip and viscous fluid cases. Pumping and trapping phenomena are discussed.     

Multiscale Modelling of Fluid and Drug Transport in Vascular Tumours

A model for fluid and drug transport through the leaky neovasculature and porous interstitium of a solid tumour is developed. The transport problems are posed on a micro-scale characterized by the inter-capillary distance, and the method of multiple scales is used to derive the continuum equations describing fluid and drug transport on the length scale of the tumour (under the assumption of a spatially periodic microstructure). The fluid equations comprise a double porous medium, with coupled Darcy flow through the interstitium and vasculature, whereas the drug equations comprise advection–reaction equations; in each case the dependence of the transport coefficients on the vascular geometry is determined by solving micro-scale cell problems.     

Interaction of peristaltic flow with pulsatile flow in a circular cylindrical tube

The effect of pulsatile flow on peristaltic transport in a circular cylindrical tube is analysed. The flow of a Newtonian viscous incompressible fluid in a flexible circular cylindrical tube on which an axisymmetric travelling sinusoidal wave is imposed, is considered. The initial flow in the tube is induced by an arbitrary periodic pressure gradient. A perturbation solution with amplitude ratio (wave amplitude/tube radius) as a parameter is obtained when the frequency of the travelling wave and that of the imposed pressure gradient are equal. The interaction effects of periodic wall induced flow and periodic pressure imposed flow are visualized through the presence of substantially different components of steady and higher harmonic oscillating flow in the first order flow solution. Numerical results show a strong variation of steady state velocity profiles with boundary wave number and Reynolds number and a strong phase shift behaviour of the flow in the radial direction.     

Density gradient system: II. A 50 channel programmable undulating diphragm peristaltic pump

A special peristaltic pump is described that has functioned as part of a system for density gradient formation and fractionation. Twenty-five pumping tubes are actuated on both the top and bottom of the pump. A Mylar diaphragm interposed between the rollers and the pumping tubes filters out the horizontal, stretching component of the forces imparted to the tubes. This greatly prolongs tube life, increases the allowable pressures that can be achieved with such a pump, and thus permits accurate delivery of viscous solutions.An explanation of the cause of the pulsations produced by peristaltic pumps is presented and the virtual elimination of these pulsations is demonstrated.Both velocity and direction of flow of the pump are controlled. By means of two independent bidirectional digital counters, preset volumes of fluid can be delivered and the total volume of liquid determined.Studies demonstrating the relative independence of fluid volume delivered at a preset count versus flow velocity and composition are presented.Other possibilities for use of the pump in automating density gradient analysis are discussed. Possibilities for employment of the pump for autoanalysis and in artificial organs are indicated.     

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.     

Soret and Dufour Effects on MHD Peristaltic Flow of Jeffrey Fluid in a Rotating System with Porous Medium

The objective of present paper is to examine the peristaltic flow of magnetohydrodynamic (MHD) Jeffrey fluid saturating porous space in a channel through rotating frame. Unlike the previous attempts, the flow formulation is based upon modified Darcy''s law porous medium effect in Jeffrey fluid situation. In addition the impacts due to Soret and Dufour effects in the radiative peristaltic flow are accounted. Rosseland’s approximation has been utilized for the thermal radiative heat flux. Lubrication approach is implemented for the simplification. Resulting problems are solved for the stream function, temperature and concentration. Graphical results are prepared and analyzed for different parameters of interest entering into the problems.     

The Onset of Double Diffusive Convection in a Viscoelastic Fluid-Saturated Porous Layer with Non-Equilibrium Model

The onset of double diffusive convection in a viscoelastic fluid-saturated porous layer is studied when the fluid and solid phase are not in local thermal equilibrium. The modified Darcy model is used for the momentum equation and a two-field model is used for energy equation each representing the fluid and solid phases separately. The effect of thermal non-equilibrium on the onset of double diffusive convection is discussed. The critical Rayleigh number and the corresponding wave number for the exchange of stability and over-stability are obtained, and the onset criterion for stationary and oscillatory convection is derived analytically and discussed numerically.     

Hydromagnetic Flow and Heat Transfer over a Porous Oscillating Stretching Surface in a Viscoelastic Fluid with Porous Medium

An analysis is carried out to study the heat transfer in unsteady two-dimensional boundary layer flow of a magnetohydrodynamics (MHD) second grade fluid over a porous oscillating stretching surface embedded in porous medium. The flow is induced due to infinite elastic sheet which is stretched periodically. With the help of dimensionless variables, the governing flow equations are reduced to a system of non-linear partial differential equations. This system has been solved numerically using the finite difference scheme, in which a coordinate transformation is used to transform the semi-infinite physical space to a bounded computational domain. The influence of the involved parameters on the flow, the temperature distribution, the skin-friction coefficient and the local Nusselt number is shown and discussed in detail. The study reveals that an oscillatory sheet embedded in a fluid-saturated porous medium generates oscillatory motion in the fluid. The amplitude and phase of oscillations depends on the rheology of the fluid as well as on the other parameters coming through imposed boundary conditions, inclusion of body force term and permeability of the porous medium. It is found that amplitude of flow velocity increases with increasing viscoelastic and mass suction/injection parameters. However, it decreases with increasing the strength of the applied magnetic field. Moreover, the temperature of fluid is a decreasing function of viscoelastic parameter, mass suction/injection parameter and Prandtl number.