# Investigation of fluid flow pattern in a 3D meandering tube

^{1}Department of Mathematics, Islamia College Peshawar, University Campus, Peshawar, Pakistan^{2}Department of Mathematics, College of Science and Humanities in Alkharj, Prince Sattam Bin Abdulaziz University, Alkharj, Saudi Arabia^{3}School of Technology, Woxsen University-Hyderabad, Hyderabad, India

Several types of meandering channels and their mathematical simulation have been proposed and discussed widely in the open literature. In the present study, the impact of a novel meandering tube geometry on streamwise vortices and pressure losses have been determined. Using a simplified Poiseuille flow simulation approach with a sinusoidal wavy meandering tube of non-uniform radius, the onset flow separation, vertex formation, and the impact of Reynolds number on field variables and stream function has been analyzed. Moreover, the linear stability theory has been implemented to trace the vertex formation. A decrease in wavelength leads to flow separation near the tube’s surface, but the flow becomes rectilinear with a sudden disturbance caused by the meander, becoming independent of vertex generating centrifugal forces. Novel insights are provided on the impact of meandering tube geometry on fluid flow and potential applications for enhancing flow conditions are suggested.

## 1 Introduction

The escalation process of heat and species transport in different flow phases with practical significance is an interesting and important subject in thermodynamics and fluid mechanics. Both the formation of fluid mixtures and changes in the functional thermophysical characteristics of the fluids fall within the category of escalation processes (Bergles and Webb, 1985; Jensen et al., 1997; Ligrani et al., 2003). Integrated heat exchangers that operate at low Reynolds numbers in case of laminar flows are treated specially for improved mixing (Webb and Bergles, 1981). It is often accepted that improved mixing may be achieved by driving a laminar-turbulent transformation under difficult conditions or by inserting vortex generators, which are efficient but have a substantial drag cost (Fiebig, 1995a; Fiebig, 1995b; Jacobi and Shah, 1995; Fiebig, 1998; Fiebig and Chen, 1999). In most cases, the hydrodynamic stabilities have been used for the transition of laminar flow to turbulent flow and, therefore, transverse grooves are used. However, a novel model has been demonstrated with the most accurate simulations, i.e., to shift laminar states without transiting to turbulent ones by using hydrodynamic instabilities. In this novel model, separating local flow is demonstrated well; the flow is driven by a slight oscillatory component and enters the resonant and stable separated shear layer (Patera and Mikic, 1986) Another type of bifurcation, i.e., centrifugal instability, has been taken into account and it produces simple vortices without suffering from the maximum drag cost associated with conventional vortex generators. Instability can improve vortex generators in limiting vortices decay, which decreases the minimum number of such generators even when the flow is only locally asymptotically stable. Both directly (by establishing a transverse transition) and indirectly (by creating a bypass transition), streamwise vortices can boost heat transfer. In most cases, a three-dimensional flow field has been developed at the end and transverse shear layers with inflection points and rapidly expanding secondary instability are formed (Floryan, 1991). Recent investigations have reported an alternate approach being taken into account of a fluid of higher thermal conductivity used with properly shaped channels while, at the same time, flow pressure being reduced immediately (Mohammadi and Floryan, 2013). This technique has provided an alternative solution to the most common development programs used on macro channels (Xu et al., 2016).

It has been known for nearly 100 years that the rotating shear layers are subjected to centrifugal instabilities. This instability has been taken into account for the situation of simple geometry and canonical flow, providing an example where it is possible to forecast the curvature of lines with ease. Using the flow state between rotating cylinders, Rayleigh introduced the inviscid technique in 1920 and determined the necessary stabilization condition through circulating distribution (Rayleigh, 1917). In 1923, Taylor included the full viscous problem and identified the crucial conditions that emerge as a result of the secondary flow (Taylor, 1923). Similar instability in curved channels was examined by Dean (Dean, 1928). The occurrence of centrifugal instability in the context of boundary layers on concave surfaces was proved by Görtler (Görtler, 1941). If the streamwise velocity distribution is not monotonic (Floryan, 1986), demonstrated that the instability is active in flows over concave as well as convex surfaces. A clear relationship between the streamline curvature and the wall curvature was provided by all of these investigations where the wall curvature was either constant or had been approximated as a constant. As a result, the critical stability condition might be expressed in terms of one parameter.

There are very few theoretical investigations that can pinpoint the starting conditions for flows in complex geometries where the wall curvature varies spatially. However, there are large number of numerical models and experimental studies that offer qualitative data (Gschwind et al., 1995). used a meandering amplitude of the same order of magnitude as the channel height, and (Nishimura et al., 1990; Tatsuo et al., 1990) used a large amplitude compared to the channel height to demonstrate the existence of streamwise vortices attributed to the centrifugal instability in sinusoidal channels. Rush et al. (Rush et al., 1999) has qualitatively identified the additional types of instability for corresponding geometries. Theoretical investigations which employed two-dimensional models were able to identify the effect of flow separation on heat transfer (Metwally and Manglik, 2004; Zhang et al., 2004) and identify the conditions that led to either single or double Hopf bifurcations or self-sustained oscillations (Guzmán et al., 2009). These models failed to identify the formation of the vortices. In their 2012 study, (Sui et al., 2012), took into account a three-dimensional rectangular channel with a very large meandering amplitude and applied numerical simulations to find a complex pattern of Dean’s vortices that changed over time and space. In the case of turbulent flow, consistent structures were found by (Pham et al., 2008). Other types of geometries, such as boundary layer flows over wavy surfaces (Saric and Ali, 1991) and Couette flows over wavy walls (Floryan, 2002), have also been found to demonstrate centrifugal instability. According to recent findings, streamwise vortices may be created by the placement of different triangular surface obstacles (Floryan and Asai, 2011) The transport of heat and fluid flow in various mediums and the studies considered the impact of magnetic fields, chemical reactions, porous media on fluid flow and thermal transport. The studies considered the impact of magnetic fields, chemical reactions, and porous media on fluid flow and thermal transport The results suggest that the addition of tri-hybrid nanoparticles and magnetic dipoles can enhance thermal transportation in Carreau Yasuda liquid, but may decrease the flow profile. This study also investigated the effects of different parameters on the peristaltic motion of hyperbolic tangent fluid in a curved compliant channel, which has potential applications in explaining blood transport dynamics. Numerical solutions and perturbation techniques were used to analyze and evaluate the results (Javed et al., 2021; Naseem et al., 2021; Wang et al., 2022). Different approaches related to the studies of fluid flow in microfluidic systems for biomedical engineering. This study analyzed the behavior of different types of fluids in different channel geometries by taking into account, convective conditions, thermal deposition effect, and chemical reactions. They investigated the impact of various parameters on flow quantities such as velocity, temperature, and concentration. This study also suggest the viability of electro-osmotic pumps for fluid flow in large osteoarticular implants (Hayat et al., 2015; Yasmin et al., 2020a; Yasmin et al., 2020b; Mehmood et al., 2020; Yasmin and Iqbal, 2021; Alyousef et al., 2023).

The main objective of the current analysis is to investigate the three-dimensional structure of the meandering geometries that cause the centrifugal force mechanism to produce streamwise vortices at the lowest possible cost as measured by pressure losses and without the interference of the traveling wave instability. This is the first comprehensive analysis of three-dimensional flow in a meandering tube that takes all potential instabilities into consideration. The new information would provide an accurate mathematical simulation for engineers and be a reasonable foundation upon which they may construct small heat exchangers that operate in the laminar flow domain. The study is divided into three main parts, specifically: i) mathematical modeling of flow in the meandering tube; ii) calculating the flow losses related to the meandering tube; and iii) identifying the geometric and flow characteristics that cause the centrifugal instability to predominate. The source in engineering procedures might be a mechanical pressure gradient. In this study, we look at the pressure gradient-driven (Poiseuille flow) flow of viscous fluid in a meandering tube with waves that is made up of fixed walls. A lot of interest has been shown in the flow of viscous fluid in a wavy meandering tube due to its applications in engineering and biological sciences, including regarding the development of muddy waves in river channels, the generation of wind waves on water and sandbanks in deserts, the movement of melting slides, rocket boosters, and the evaporation of film in burning chambers. Furthermore, physiologists and technicians have often attempted to create and explain blood and urine flow in terms of meandering channels (tubes).

## 2 Geometry of the problem

The flow of a viscous fluid in a meandering tube, whose geometry is shown in Figure 1, is the main focus of this paper.

**FIGURE 1**. Geometry of the meandering tube under consideration and flow regime with coordinate

The tube of variable radius is taken, whose radius is determined by

## 3 Governing equations

In order to express the current problem in proper coordinates, it is appropriate to choose a suitable coordinate system for the simulated problem. The well-established relationship between the Cartesian coordinates

where

The governing equations are non-dimensionalized by introducing dimensionless variables and the length is non-dimensionalized by

Note that the asterisk “*” represents the dimensional quantities and

**Continuity** equation,

The no slip boundary conditions at the wall and the symmetry conditions at the center of the meandering tube are given as:

where

By using the dimensionless variables as defined above, the continuity Eq. 1 is transformed as:

From above definitions of dimensionless variables, we have:

where

Where Eqs. 10–13 represent the dimensionless form of **continuity** and **components of Navier-Stokes** equations, respectively.

The dimensionless form of no slip boundary conditions at the wall and the symmetry conditions at the center of the meandering tube are obtained as:

where

## 4 Modal problem

The velocity vector *,* and

Further, the velocity vector

## 5 Solution of the problem

Next, our aim is to determine the solution of dimensionless Eqs. (10 - 16). As these equations have a parameter

Similarly, the dimensionless Eq. 16 has transformed (by introducing

In the case of small amplitude waviness, i.e.,

where

The values of

## 6 Results and discussion

The solution of the zeroth-order system, which is obtained by putting value from Eq. 23 into Eqs 18–21 and equating like powers of

Note that the solution in Eq. 24 for the fully developed flow in a straight duct has been reported in F.M. White (White and Majdalani, 2006) and Schlichting (Schlichting and Kestin, 1961).

Similarly, the first-order system is obtained by equating like powers of

where the functions

Moreover, the coefficients

For

Furthermore, the coefficients of the second series in the expression of

The recursive formula for

For

The recursive formula for

For

And the recursive formula for

For

Figure 2, Figure 3, Figure 4, Figure 5 describe the impact of meanders on the formation and behavior of vortices in the flow, producing complex patterns of stream contours. Vortices form near the wall of the tube and interact with the fluid flow. In Figure 6, fluid flow patterns are affected in several ways with variation in the number of helixes. As the number of helixes increases, the amplitude of the meanders increases, leading to more complex flow patterns with multiple recirculation zones near the wall of the tube. This also causes an increase in the pressure drop and the overall mixing in the flow.

**FIGURE 2**. The stream lines are drawn by using *ß* = 1, z = 1, and k = 10 in the domain 0 ≤ θ ≤ 2, −1 ≤ η ≤ 0 at the lower portion of the tube.

**FIGURE 3**. The stream lines are drawn by using *ß* = 1, z = 1, and k = 10 in the domain 0 ≤ θ ≤ 2, 0 ≤ η ≤ 1 at the upper portion of the tube.

**FIGURE 4**. The stream lines are drawn by using *ß* = 1, z = 1, and k = 0.1 in the domain 0 ≤ θ ≤, 0 ≤ η ≤ 1 at the upper portion of the tube.

**FIGURE 5**. The stream lines are drawn by using *ß* = 1, z = 1, and k = 10 in the domain 0 ≤ θ ≤ 2, −1 ≤ η ≤ 1 of the whole tube.

**FIGURE 6**. The stream lines are drawn by using **(A)** *ß* = 2, **(B)** *ß* = 5, and **(C)** *ß* = 0.1 in the domain 0 ≤ *z* ≤ 2, −1 ≤ η ≤ 1 of the whole tube.

## 7 Conclusion

In a meandering tube, streamwise vortices can improve the transportation of heat and species mass in a transverse direction. To create these vortices, a specific meandering wavelength is required, which can be measured using linear instability theory. Short wavelengths cause flow separation, while long wavelengths result in a rectilinear stream resembling a flat plate. Shear-driven instability also contributes to vortex formation. The effect of a steady laminar flow in a meandering tube with small amplitude wavy walls was studied. The governing equations were constructed using the continuity equation and Navier-Stokes equations with no slip boundary conditions. Suppositions were made to simplify the complex non-linear problem, including linear instabilities, uniform thermal characteristics, and laminar flow conditions. The equations were further simplified using available dimensionless variables and new transformations. The perturbation and power series approaches were used to solve the equations with the help of Mathematica, and the velocity profiles and stream contours were graphed using standard codes and definitions.

The flow in a meandering tube becomes unstable due to centrifugal impact, which results in the formation of streamwise vortices. This research has focused on identifying the lowest meandering amplitude that can create these vortices with minimal pressure loss. Two primary forms of vortex instability have been studied, with the most effective parameter being better for creating vortex instability than centrifugal instability. The meandering geometry that is most successful in producing vortices does not encourage traveling wave instability, which can delay the onset of laminar-turbulent transition. By using the reduced geometry model approach, the findings can be applied to various types of meandering tubes.

Further investigations can be conducted through experiments, simulations, and new applications. Experiments can provide valuable insights into the complex nature of the flow using advanced techniques like PIV, LDV, and HWA to measure velocities, turbulence, and vortices. Numerical simulations using CFD techniques can also be performed to investigate flow behavior, providing detailed information on velocity, pressure, and vortices, and help optimize tube design.

## Data availability statement

The original contributions presented in the study are included in the article/supplementary material, further inquiries can be directed to the corresponding author.

## Author contributions

The mathematical model has been proposed by DK, all the numerical computations and their graphs have been carried out by SI. The discussion of graphs and their physical interpretation has been given by DK and SI. The literature review and comparison of the present simulations with the classical data has been established by NU. The final review and amendments in the manuscript has been carried out by KN.

## Funding

This study is supported via funding from Prince Sattam bin Abdulaziz University project number (PSAU/2023/R/1444).

## Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

## Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

## References

Alyousef, H. A., Yasmin, H., Shah, R., Shah, N. A., El-Sherif, L. S., and El-Tantawy, S. A. (2023). Mathematical modeling and analysis of the steady electro-osmotic flow of two immiscible fluids: A biomedical application. *Coatings* 13 (1), 115. doi:10.3390/coatings13010115

Bergles, A. E., and Webb, R. L. (1985). *A guide to the literature on convective heat transfer augmentation*. Ames, Iowa: Iowa State University College of Engineering.

Dean, W. R. (1928). Fluid motion in a curved channel. *Proc. R. Soc. Lond. Ser. A, Contain. Pap. a Math. Phys. Character* 121, 402–420.

Fiebig, M., and Chen, Y. (1999). “Heat transfer enhancement by wing-type longitudinal vortex generators and their application to finned oval tube heat exchanger elements,” in *Heat transfer enhancement of heat exchangers* (Dordrecht: Springer), 79–105.

Fiebig, M. (1995). Embedded vortices in internal flow: Heat transfer and pressure loss enhancement. *Int. J. Heat Fluid Flow* 16, 376–388. doi:10.1016/0142-727x(95)00043-p

Fiebig, M. (1995). Vortex generators for compact heat exchangers. *J. Enhanc. Heat Transf.* 2, 1–20. doi:10.1615/jenhheattransf.v24.i1-6.10

Fiebig, M. (1998). Vortices, generators and heat transfer. *Chem. Eng. Res. Des.* 76, 108–123. doi:10.1205/026387698524686

Floryan, J. M., and Asai, M. (2011). On the transition between distributed and isolated surface roughness and its effect on the stability of channel flow. *Phys. Fluids* 23, 104101. doi:10.1063/1.3644694

Floryan, J. M. (2002). Centrifugal instability of Couette flow over a wavy wall. *Phys. fluids* 14, 312–322. doi:10.1063/1.1416185

Floryan, J. M. (1986). Görtler instability of boundary layers over concave and convex walls. *Phys. fluids* 29, 2380–2387. doi:10.1063/1.865531

Floryan, J. M. (1991). On the Görtler instability of boundary layers. *Prog. Aerosp. Sci.* 28, 235–271. doi:10.1016/0376-0421(91)90006-p

Görtler, H. (1941). Instabilität laminarer Grenzschichten an konkaven Wänden gegenüber gewissen dreidimensionalen Störungen. *ZAMM-Journal Appl. Math. Mechanics/Zeitschrift für Angewandte Math. und Mech.* 21, 250–252. doi:10.1002/zamm.19410210408

Gschwind, P., Regele, A., and Kottke, V. (1995). Sinusoidal wavy channels with Taylor-Goertler vortices. *Exp. Therm. fluid Sci.* 11, 270–275. doi:10.1016/0894-1777(95)00056-r

Guzmán, A. M., Cárdenas, M. J., Urzúa, F. A., and Araya, P. E. (2009). Heat transfer enhancement by flow bifurcations in asymmetric wavy wall channels. *Int. J. Heat Mass Transf.* 52, 3778–3789. doi:10.1016/j.ijheatmasstransfer.2009.02.026

Hayat, T., Iqbal, M., Yasmin, H., Alsaadi, F. E., and Gao, H. (2015). Simultaneous effects of Hall and convective conditions on peristaltic flow of couple-stress fluid in an inclined asymmetric channel. *Pramana* 85, 125–148. doi:10.1007/s12043-014-0888-1

Jacobi, A. M., and Shah, R. K. (1995). Heat transfer surface enhancement through the use of longitudinal vortices: A review of recent progress. *Exp. Therm. Fluid Sci.* 11, 295–309. doi:10.1016/0894-1777(95)00066-u

Javed, M., Imran, N., Arooj, A., and Sohail, M. (2021). Meta-analysis on homogeneous-heterogeneous reaction effects in a sinusoidal wavy curved channel. *Chem. Phys. Lett.* 763, 138200. doi:10.1016/j.cplett.2020.138200

Jensen, M. K., Bergles, A. E., and Shome., B. (1997). The literature on enhancement of convective heat and mass transfer. *J. Enhanc. Heat Transf.* 4, 1–6. doi:10.1615/jenhheattransf.v4.i1.10

Ligrani, P. M., Oliveira, M. M., and Blaskovich, T. (2003). Comparison of heat transfer augmentation techniques. *AIAA J.* 41, 337–362. doi:10.2514/2.1964

Mehmood, O. U., Qureshi, A. A., Yasmin, H., and Uddin, S. (2020). Thermo-mechanical analysis of non Newtonian peristaltic mechanism: Modified heat flux model. *Phys. A Stat. Mech. its Appl.* 550, 124014. doi:10.1016/j.physa.2019.124014

Metwally, H. M., and Manglik, R. M. (2004). Enhanced heat transfer due to curvature-induced lateral vortices in laminar flows in sinusoidal corrugated-plate channels. *Int. J. Heat Mass Transf.* 47, 2283–2292. doi:10.1016/j.ijheatmasstransfer.2003.11.019

Mohammadi, A., and Floryan, J. M. (2013). Groove optimization for drag reduction. *Phys. Fluids* 25, 113601. doi:10.1063/1.4826983

Naseem, T., Nazir, U., and Sohail, M. (2021). Contribution of Dufour and Soret effects on hydromagnetized material comprising temperature-dependent thermal conductivity. *Heat. Transf.* 50 (7), 7157–7175. doi:10.1002/htj.22222

Nishimura, T., Yano, K., Yoshino, T., and Kawamura, Y. (1990). Occurrence and structure of Taylor–Goertler vortices induced in two-dimensional wavy channels for steady flow. *J. Chem. Eng. Jpn.* 23, 697–703. doi:10.1252/jcej.23.697

Patera, A. T., and Mikic, B. B. (1986). Exploiting hydrodynamic instabilities. Resonant heat transfer enhancement. *Int. J. heat mass Transf.* 29, 1127–1138. doi:10.1016/0017-9310(86)90144-4

Pham, M. V., Plourde, F., and Doan, S. K. (2008). Turbulent heat and mass transfer in sinusoidal wavy channels. *Int. J. Heat Fluid Flow* 29 (5), 1240–1257. doi:10.1016/j.ijheatfluidflow.2008.04.002

Rayleigh, L. (1917). On the dynamics of revolving fluids. *Proc. R. Soc. Lond. Ser. A, Contain. Pap. a Math. Phys. Character* 93, 148–154.

Rush, T. A., Newell, T. A., and Jacobi, A. M. (1999). An experimental study of flow and heat transfer in sinusoidal wavy passages. *Int. J. heat mass Transf.* 42, 1541–1553. doi:10.1016/s0017-9310(98)00264-6

Saric, W. S., and Ali, B. (1991). “Görtler vortices with periodic curvature,” in *Boundary layer stability and transition to turbulence* (Berlin, Germany: Springer), 37–42.

Sui, Y., Teo, C. J., and Lee, P. S. (2012). Direct numerical simulation of fluid flow and heat transfer in periodic wavy channels with rectangular cross-sections. *Int. J. Heat Mass Transf.* 55, 73–88. doi:10.1016/j.ijheatmasstransfer.2011.08.041

Tatsuo, N., Shinichiro, M., Shingho, A., and Yuji, K. (1990). Flow observations and mass transfer characteristics in symmetrical wavy-walled channels at moderate Reynolds numbers for steady flow. *Int. J. Heat Mass Transf.* 33, 835–845. doi:10.1016/0017-9310(90)90067-5

Taylor, G. I. (1923). VIII. Stability of a viscous liquid contained between two rotating cylinders. *Philosophical Trans. R. Soc. Lond. Ser. A, Contain. Pap. a Math. or Phys. Character* 223, 605–615.

Wang, F., Sohail, M., Nazir, U., El-Zahar, E. R., Park, C., and Jabbar, N. (2022). An implication of magnetic dipole in Carreau Yasuda liquid influenced by engine oil using ternary hybrid nanomaterial. *Nanotechnol. Rev.* 11 (1), 1620–1632. doi:10.1515/ntrev-2022-0100

Webb, R. L., and Bergles, A. E. (1981). *Performance evaluation criteria for selection of heat transfer surface geometries used in low Reynolds number heat exchangers*. Ames, Iowa: Iowa State University College of Engineering.

Xu, M., Lu, H., Gong, L., Chai, J. C., and Duan, X. (2016). Parametric numerical study of the flow and heat transfer in microchannel with dimples. *Int. Commun. Heat Mass Transf.* 76, 348–357. doi:10.1016/j.icheatmasstransfer.2016.06.002

Yasmin, H., and Iqbal, N. (2021). Convective mass/heat analysis of an electroosmotic peristaltic flow of ionic liquid in a symmetric porous microchannel with soret and dufour. *Math. Problems Eng.* 2021, 1–14. doi:10.1155/2021/2638647

Yasmin, H., Iqbal, N., and Hussain, A. (2020). Convective heat/mass transfer analysis on johnson-segalman fluid in a symmetric curved channel with peristalsis: Engineering applications. *Symmetry* 12 (9), 1475. doi:10.3390/sym12091475

Yasmin, H., Iqbal, N., and Tanveer, A. (2020). Engineering applications of peristaltic fluid flow with hall current, thermal deposition and convective conditions. *Mathematics* 8 (10), 1710. doi:10.3390/math8101710

Zhang, J., Kundu, J., and Manglik, M. (2004). Effect of fin waviness and spacing on the lateral vortex structure and laminar heat transfer in wavy-plate-fin cores. *Int. J. Heat Mass Transf.* 47, 1719–1730. doi:10.1016/j.ijheatmasstransfer.2003.10.006

## Nomenclature

Keywords: stream wise vortices, meandering tube, centrifugal forces, flow separation, instability

Citation: Ibrahim S, Khan Marwat DN, Ullah N and Nisar KS (2023) Investigation of fluid flow pattern in a 3D meandering tube. *Front. Mater.* 10:1187986. doi: 10.3389/fmats.2023.1187986

Received: 16 March 2023; Accepted: 18 April 2023;

Published: 09 June 2023.

Edited by:

Noor Saeed Khan, University of Education Lahore, PakistanReviewed by:

Humaira Yasmin, King Faisal University, Saudi ArabiaMuhammad Sohail, Khwaja Fareed University of Engineering and Information Technology (KFUEIT), Pakistan

Copyright © 2023 Ibrahim, Khan Marwat, Ullah and Nisar. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Syed Ibrahim, syed.ibrahim@riphah.edu.pk