Modeling and Optimization of Peristaltic Fluid Transport in Axisymmetric, Porous Structures 1 Utilities
and Sustainable Engineering, The University of Trinidad and Tobago, Trinidad
and Tobago
1. INTRODUCTION Peristaltic transport is the movement of fluid achieved by alternating compression and relaxation along a distensible tube, which propels the fluid. The principle of peristaltic transport of fluid is extensively utilized in industrial and commercial pumps. Peristaltic pumps use rotating rollers pressed against a distensible tube, creating a pressurized flow. The fluid is moved through the tube with each rotating motion. The major advantage of a peristaltic pump is that the only contact the fluid being pumped has with the peristaltic pump is in interior of the distensible tube. This allows isolation of the fluid from the environment of the pump and is essential in the case of sterilized or hazardous fluid being transported Gupta & Gupta (2007). Peristaltic transport of fluid is common within the human body. Example of this include: urine flow from the kidney to the bladder, swallowing through the esophagus, movement of chyme in the gastro-intestinal tract, intra-uterine fluid motion, flow of spermatozoa in the ductus efferentes of the male reproductive tract, movement of ovum in the female fallopian tube, transport of lymph in the lymphatic vessels and the vasomotion of small blood vessels such as arterioles, venules and capillaries Shit et al. (2010). Earthworms also use peristalsis as a means of locomotion. The current study will focus on peristaltic transport in the human body, in particular transport in the gastro-intestinal tract and capillaries. Epithelial cells line major cavities and organs of the body. These include the stomach and small intestine, kidney, pancreas, and esophagus. Epithelial cells also form the structure of the lung. The epithelial cells that line capillaries are referred to as endothelial cells. Epithelial cells serve many functions which include transportation and absorption of nutrients. Absorption of nutrients by epithelial cells is accomplished by a process called transcytosis. Transcytosis occurs as membrane-bound carriers selectively transport materials between one part of the cell and another to maintain unique environments on either side of the cell Barnes et al. (2008). The selective absorption of nutrients or materials by epithelial cells facilitates modeling epithelial cells as porous media. This motivates the formulation of the current model. Early studies in peristaltic pumping was carried out by Shapiro et al. (1969). They investigated peristaltic pumping under small Reynolds number and a large tube aspect ratio to allow uniform pressure distribution over the tube cross section. They presented results for varying geometric parameters and introduced, explained, and illustrated reflux flow. They also referred to an experiment that confirmed their findings. They applied their model for peristaltic pumping to the gastro-intestinal system. Barton & Raynor (1968) explored the case where the wall disturbance wavelength is much larger than the average tube radius, and when the disturbance wavelength may be as small as the average radius. They used numerical techniques to determine the relation between average flow rate and pressure differential across a wavelength. Wei et al. (2003) investigated pressure driven flow of a Newtonian fluid in a two dimensional wavy, poro-elastic, cartesian walled channel. Biphasic mixing theory is employed for the poro-elastic layer and the flow is solved by a lubrication approximation using the aspect ratio of the channel. They obtained the velocity fields using a perturbation approach and discussed the fluid flow and geometrical parameters that develop and influence trapping or recirculation eddies within the Newtonian fluid. Mishra & Rao (2008) studied peristaltic transport in an axisymmetric tube, with a Brinkman porous peripheral region and a Newtonian fluid core region, using the lubrication approach. They allowed for deformation of the porous-fluid interface by determining the interface as part of the solution using the conservation of mass in both the porous medium and fluid medium separately. They used a shear stress jump condition at the interface between the peripheral and core regions. Pumping characteristics, trapping, and reflux phenomena are discussed for various parameters governing the flow. In this present study, peristaltic transport of an incompressible Newtonian fluid in an axisymmetric, porous, corrugated, or wavy tube is investigated. The porous layer is modelled as a Darcy region with a transition Brinkman layer. The analysis is performed in a wave frame of reference. Lubrication theory as well as perturbation analysis is used to solve the velocity field and pressure. The boundary conditions used at the Brinkman-fluid interface include jump in tangential stress, continuity or normal stresses and continuity of radial and axial velocities. The boundary conditions at the Brinkman-Darcy interface include the Beavers and Joseph boundary condition and continuity of normal stresses. The phenomena of reflux or recirculation is observed and its dependency on the fluid flow parameters of hydraulic resistivity, tangential stress jump and the Darcy slip parameter is studied from streamline and velocity field plots. 2. MODEL FORMULATION Peristaltic pumping is modeled as the motion of an infinite wave train along an axis symmetric porous cylindrical tube. The tube core is filled, and the porous layer saturated with the same incompressible viscous Newtonian fluid. The equation of the infinite wave train at the outer walls of the porous media, in cylindrical polar coordinates is given by: where The equation representing the Core - Brinkman interface is: The equation representing the Brinkman - Darcy interface is: Figure 1
Following the analysis of Shapiro
et al. (1969), the fixed frame of
reference coordinates In the fixed frame of reference, the outer wall of the
Darcy region, moves with velocity 2.1. Governing Equations - the Core and Porous region 2.1.1. Dimensionless Analysis The analysis is performed on the wave frame of reference. Following a similar analysis to Wei et al. (2003) the dimensionless quantities are: For the pressure p: Where The aspect ratio And Reynolds number The interfaces For simplicity, the hats are dropped from the dimensionless variables in the analysis that follows. 2.1.2. Dimensionless Core region In this region, we use the continuity equation: and the Navier Stokes equation: 2.1.3. Dimensionless Brinkman Region In the Brinkman region the fluid in the porous layer and The equations are: 2.1.4. Dimensionless Darcy region In the Darcy region 2.1.5. Dimensionless Boundary Conditions The boundary conditions for the interfaces At 1) The
axial core velocity At Jump in tangential stresses: where 2) Continuity of normal stresses: 3) Continuity of radial and axial velocities: At 1) Beavers and Joseph boundary condition on the axial velocity. where 2) Continuity of normal stresses: At 1) There is no penetration of the radial velocity at the outer wall. 2) There is a non zero negative axial velocity Shapiro et al. (1969) The total characteristic volumetric flow rate along the
tube is given by, Writing Eq3.3.24 in dimensionless form dropping the superscript hats: 2.1.6. Asymptotic Expansion The geometry of the tube and the wave train is such that
the average wavy variation of the wave train is small along the axial
direction. The aspect ratio The dimensionless governing equations and the
corresponding boundary conditions do not contain expressions in 2.1.7. Leading Order System Core region: Brinkman region: Darcy Region: Boundary
Conditions: 2.1.8. Solution to the leading order system Core region Brinkman region Darcy region where the terms: 2.1.9. Obtaining the Constants The boundary conditions are now applied and expressions
for the constants where: 3. RESULTS AND DISCUSSIONS The analytical expressions obtained for the stream functions and velocities for the different regions is now presented in their corresponding domains. These domains are complex two-dimensional regions, which comprises sinusoidal boundaries. The regions are meshed, using a triangular mesh, and the co-ordinates of each node substituted into the analytical expression for the stream function to obtain the velocity field and stream function plots that follow. Figure 2
The analysis that follows is done in the wave frame of
reference. The ratio, The values of the other constants used are: It should be noted that from In the streamline plots that follow, the fluid core region streamlines are red, the thin transition Brinkman layer is green and the Darcy region has purple streamlines. The color of the arrows in the velocity field plots represents velocity magnitudes. Dark blue represents the lowest velocities and red the highest velocities. Figure 3
Figure 4
Figure 5
3.1. Discussion Trapping and
Circulation as The phenomena of trapping was discussed by Shapiro et al. (1969). They defined trapping as the region of fluid having a mean speed of advance equal to the wave speed. The fluid is hence trapped in that it moves at least in the mean, with the wave itself. In the trapping region there is a recirculation of fluid particles. Trapping is observed in the fluid core region of the
streamlines and velocity field plots for varying hydraulic resistivity It is important to study the magnitude of the velocity fields, given by the colors of the field arrows. Within the circulation regions, the fluid velocity is at its slowest, this is illustrated by the blue color of the velocity field. Fluid or in the case of the human body, material, or nutrients, in this region, circulate at a slow speed or have a greater resident time. This allows for more efficient mixing and absorption of nutrients into the porous region underneath the peaks. Wei et al. (2003) commented that the circulation regions they observed inside the peaks of the poroelastic region allows a greater resident time for fluid particles near the adjacent trough. This is not the case for the core fluid in this model, in fact from the velocity field plots, the troughs adjacent to the circulating peak regions have a higher velocity than the circulating fluid. Geometric effects
on trapping and circulation. From Figure 4, the amplitude Affects of porous
layer thickness, Porous layer thickness 4. CONCLUSIONS And RECOMMENDATIONS This study has explored the intricate dynamics of peristaltic transport of an incompressible Newtonian fluid within an axisymmetric, porous, corrugated tube, utilizing a comprehensive analytical framework that incorporates porous layer theory and perturbation analysis. Our investigation reveals the critical role of hydraulic resistivity, geometric parameters of the wavy porous structure, and the porous layer thickness in shaping the transport characteristics, particularly in terms of trapping and circulation phenomena. The study demonstrates that the hydraulic resistivity parameter, α, plays a pivotal role in the formation of circulation regions within the fluid core, with significant implications for the efficiency of mixing and nutrient absorption processes. This finding not only corroborates the observations made by previous studies but also extends our understanding by highlighting the influence of the Brinkman model in capturing the nuances of fluid behavior in the core region adjacent to the porous peaks. Furthermore, the geometric configuration of the wavy porous layer, especially its amplitude and thickness, is shown to have a profound impact on the trapping and circulation regions. An increase in the amplitude of the wavy porous region amplifies the trapping phenomenon, while the thickness of the porous layer dictates the extent and positioning of circulation regions, thereby affecting the overall transport efficiency within the tube. This work contributes significantly to the body of knowledge on peristaltic pumping, particularly in the context of biological and industrial applications where the isolation and precise control of fluid transport are paramount. The insights garnered from our model underscore the potential of peristaltic transport mechanisms to optimize fluid flow and enhance the delivery and mixing of substances, whether in medical devices, industrial processes, or within the human body.
CONFLICT OF INTERESTS None. ACKNOWLEDGMENTS None. REFERENCES Barnes, P. J., Drazen, J. M., Rennard, S. I., & Thomson, N. C. (2008). Asthma and COPD: Basic Mechanisms and Clinical Management. Academic Press. Barton, C., & Raynor, S. (1968). Peristaltic Flow in Tubes. The Bulletin of Mathematical Biophysics, 30(4), 663–680. https://doi.org/10.1007/BF02476682 Chhabra, R. P., & Richardson, J. F. (2008). Non-Newtonian Flow and Applied Rheology: Engineering Applications. Butterworth-Heinemann. Goharzadeh, A., Khalili, A., & Jo̸rgensen, B. B. (2005). Transition Layer Thickness at a Fluid-Porous Interface. Physics of Fluids, 17(5), 057102. https://doi.org/10.1063/1.1894796 Gupta, A. K., & Gupta. (2007). Industrial Automation and Robotics. Firewall Media. Mishra, M., & Rao, A. R. (2008). Peristaltic Flow of a Two-Layer System in a Poroflexible Tube. Journal of Porous Media, 11(1), 51–71. https://doi.org/10.1615/JPorMedia.v11.i1.40 Pozrikidis, C. (2005). Axisymmetric Motion of a File of Red Blood Cells Through Capillaries. Physics of Fluids, 17(3), 031503-031503–031514. https://doi.org/doi:10.1063/1.1830484 Shapiro, A. H., Jaffrin, M. Y., & Weinberg, S. L. (1969). Peristaltic Pumping with Long Wavelengths at Low Reynolds Number. Journal of Fluid Mechanics, 37(04), 799–825. https://doi.org/10.1017/S0022112069000899 Shit, G. C., Roy, M., & Ng, E. Y. K. (2010). Effect of Induced Magnetic Field on Peristaltic Flow of a Micropolar Fluid in an Asymmetric Channel, 1007.0923. https://doi.org/10.1002/cnm.1397 Wei, H. H., Waters, S. L., Liu, S. Q., & Grotberg, J. B. (2003). Flow in a Wavy-Walled Channel Lined with a Poroelastic Layer. Journal of Fluid Mechanics, 492, 23–45. https://doi.org/10.1017/S0022112003005378
© IJETMR 2014-2024. All Rights Reserved. |