Skip to main content

ORIGINAL RESEARCH article

Front. Mater., 26 July 2023
Sec. Colloidal Materials and Interfaces
This article is part of the Research Topic Thermodynamics of Newtonian and non-Newtonian Nanofluids with Recent Advancements View all 13 articles

Flow and heat transfer in a meandering channel

Syed Ibrahim,
Syed Ibrahim1,2*Dil Nawaz Khan MarwatDil Nawaz Khan Marwat1Naeem UllahNaeem Ullah1Kottakkaran Sooppy Nisar,Kottakkaran Sooppy Nisar3,4 Kamran Kamran1
  • 1Department of Mathematics, Islamia College Peshawar, University Campus, Peshawar, Khyber Pakhtunkhwa, Pakistan
  • 2Riphah International University, Chakdara, Khyber Pakhtunkhwa, Pakistan
  • 3Department of Mathematics, College of Science and Humanities in Alkharj, Prince Sattam Bin Abdulaziz University, Alkharj, Saudi Arabia
  • 4School of Technology, Woxsen University, Hyderabad, India

Fluid flows occur due to internal or external forces such as wind, gravity, pressure gradients, side-wall motion, MHD, and free convection. This study examines how meanders impact heat transfer by studying the behavior of viscous fluid flow with streamwise vortices in a sinusoidal wavy meandering channel of non-uniform radius. The study simplifies the motion and energy equations governing the fluid flow using novel transformations and a regular perturbation method. By plotting graphs for different parameter values, such as Pr, Re, and Ec, it reveals that decreasing the wavelength leads to flow separation near the channel surface. However, the stream moves forward with a sudden meander disturbance, causing the flow to become rectilinear and independent of vertex-generating centrifugal forces. The study identifies a stream function using standard and established relations. The fluid flow patterns and temperature distribution behavior are shown in various plots, highlighting the significant impact of meanders on fluid flow.

1 Introduction

There are infinitely many possible types of meanders, and here we consider only a special type of meandering channel, i.e., the sinusoidal meandering channel. These channels are influenced by various agents that cause fluid motion, with the most common being the pressure gradient. In some cases, the pressure gradient may be mechanical in nature. The flows within these channels are predominantly induced by the pressure gradient, and they have numerous practical applications (Webb and Bergles, 1981; Bergles and Webb, 1985; Jensen et al., 1997; Ligrani et al., 2003). The partition of escalation structures can be performed on those that are responsible for necessary changes in the thermal and physical properties of fluids under consideration and the ones marked by perfection in fluid maxing (Fiebig, 1995a; Fiebig, 1995b; Jacobi and Shah, 1995). Mixing with improved mechanics is the point of attention in engineering devices functioning at small Reynolds numbers (Re) with laminar patterns (Fiebig, 1998). Commonly, it is understood that improved mixing can be done via compelling the fluid from the laminar state into a turbulent one or by the formation vortex generators. Vortex generators are functional but have significant pressure in their drag drawback (Patera and Mikic, 1986; Fiebig and Chen, 1999).

For nearly a century, centrifugal instability has been observed in shear layers. Previous studies focused on simple geometries and canonical flows, making it easy to calculate the curvature or meander of flow patterns. One example is the flow of fluid motion between rotating cylinders, where (Floryan, 1991) investigated the stability states of the distribution of rotation in an inviscid mechanism. In (Mohammadi and Floryan, 2013), addressed a highly viscous problem was addressed, and the critical conditions of the secondary flow onset were determined. In (Xu et al., 2016), comparable instability in curved channels was analyzed, and (Rayleigh, 1917) presented a case of boundary layer flows on concave sheets, taking into account centrifugal instability. In (Taylor, 1923), it was found that the special type of instability is prominent for fluid flow over concave and convex sheets, subject to the condition that the velocity components are non-monotonic in nature. In (Ghalambaz et al., 2016; Hayat et al., 2017; Srinivasacharya and Sibanda, 2020; Mahmud and Uddin, 2021), numerical simulation and analysis of different types of nanofluids flowing past various types of surfaces in the context of fluid flow and heat transfer in different channels was performed (Khan and Ahmed, 2015; Alsaedi et al., 2016; Nadeem et al., 2016; Izadi and Pourmehran, 2017; Khan et al., 2017; Sheikholeslami et al., 2017). These studies explore the impact of various factors such as thermal radiation, magnetic fields, bioconvection, entropy optimization, and chemical reactions on flow and heat transfer characteristics. As such, these articles are highly relevant for those interested in studying the behavior of nanofluids. All of these results apply to constant meanders and establish a clear relationship between streamline meanders and the wall meander. Consequently, the demonstration of this special type of stability revealed that it is represented by a single parameter, which is one possible way to describe this meandering with a single parameter.

The main objective of this study is to investigate the flow of a viscous fluid in a meandering channel with wavy walls driven by a pressure gradient. The equations of motion, including the conservation of the mass and momentum equations and the energy equation, are utilized to analyze the fluid flow and temperature distribution within the channel. The boundary conditions are imposed on the channel walls and at the center of the channel to maintain symmetry. Dimensionless variables are used to transform the governing equations into a dimensionless form for easy tracking of units. Additionally, the study explores special types of instabilities and their significance. In the first part, an assumption of fixed pressure gradient is imposed, assuming direct channel flow with parallel plates, and sinusoidal channel flow is driven by an identical pressure gradient. Flow rate variation is used to supplement the flow losses associated with the meandering channel. The generalized model is solved using perturbation method, and the stream function is calculated and examined using other techniques. The temperature distribution behavior is illustrated through various plots depicting different dimensionless parameters, such as Pr, Re, and Ec.

2 Geometry of the problem

Consider the flow of a viscous fluid in a meandering channel, which consists of two wavy walls separated by a fixed gap. The classical model for laminar flows in channels, tubes, and ducts is provided by Poiseuille. The flow in such a channel is primarily induced by a pressure gradient. There are various types of flows that pass through ducts and channels, and they have significant practical importance. The classification of scaling structures can be conducted based on those that have variations in the effective thermophysical properties of fluids, and the pursuit of optimal fluid mixing can lead to an infinite number of possible meander types. In this study, we will focus on the simplest sinusoidal channel, whose geometry is illustrated in Figure 1, where the superscript H denotes the upper and G denotes the lower walls of the channel. The channel extends along the x-axis from negative infinity to positive infinity, and the flow is driven by a pressure gradient.

FIGURE 1
www.frontiersin.org

FIGURE 1. Geometry of meandering channel and flow field under consideration.

3 Governing equations and their non-dimensional form

The general form of continuity equation for incompressible viscous flow is:

divV=0(1)

The velocity vector, denoted as V, has three components for three-dimensional flow. Similarly, the general form of the Navier–Stokes equations applies to steady and viscous flows of constant viscosity.

ρVV=P+μ2V,(2)

The variables ρ, μ, and P represent density, viscosity, and pressure of the fluid, respectively. To define the dimensionless variables, we use asterisks and reference length L and velocity U, as follows:

*=L,V*=VU,andP*=PρU2,(3)

After substituting the non-dimensional variables, which are defined in Eq. 3, into Eqs 1, 2, we obtain the dimensionless continuity and Navier–Stokes equations:

divV*=0,(4)
ρU2LV**V*=ρU2L*P*+μUL2*2V*,(5)

The above equation is multiplied by LρU2, and when simplified we get:

V**V*=*P*+1Re*2V*.(6)

where Re=ρULμ represents the Reynold’s number.

4 Modeling of the problem

Let us consider that for flow in a straight channel and for two-dimensional flow, the velocity is represented by the velocity V0, pressure by P0, the stream function of the flow is denoted by ψ0, and Q0 is the flow rate representative. Then, the continuity equation in two dimensional is transformed to:

u0x+v0y=0,(7)

The x- component of the Navier–Stoke equations is:

ρu0t+u0u0x+v0u0y=P0x+μ2u0x2+2u0y2,(8)

The y- component of the Navier–Stoke equations is:

ρv0t+u0v0x+v0v0y=P0y+μ2v0x2+2v0y2,(9)

The boundary conditions for the problem are established by utilizing the no-slip condition, along with considering the geometry of the problem. Thus, the boundary conditions can be written as:

u0=0aty=±1.(10)

Equations 79 represent the continuity and the x&y- momentum equations, respectively, for a one-dimensional straight channel. These equations can be greatly simplified by making certain assumptions, namely, that u00andv0=0. With this assumption, the solution for the governing equations for straight channel flow can be obtained by combining the axial velocity u0, normal velocity v0, and pressure P0, which can be expressed as follows:

V0=u0η,0=1η2,0P0=2ξRe+C1,(11)

Note that Eqs 111 for the fully developed flow in a straight channel has been reported in F.M. White [31] and Schlichting [32].

The volume flow rate for the channel flow of fixed width is:

Q0=11u0ηdη=43.(12)

The stream function ψ0 is simply evaluated as:

ψ0η=u0ηψ0=η33+η+C2.(13)

In a meandering channel, where the fluid flow moves in the positive x-axis direction, the Reynolds number is defined based on the maximum x-velocity and channel half-height.

The velocity field and other related field quantities for the fluid motion in the meandering channel can be expressed as follows:

Vb=V0+V1=u0η,0+u1ξ,η,v1ξ,η,Pbξ,η=P0ξ+P1ξ,η,ψb=ψ0η+ψ1ξ,η,Qb=Q0+Q1,(14)

In the context of the meandering channel geometry, the velocity V1, pressure P1, stream function ψ1, and flow rate Q1 are dimensionless quantities that are adjusted accordingly. Additionally, the velocity vector Vb=ub,vb is decomposed. These components are further decomposed as follows:

ub=u0η,u1ξ,ηandvb=0,v1ξ,η],(15)

This is achieved by substituting the new variables for the velocity field in the dimensionless governing equations, i.e., Eqs 4, 6, and eliminating pressure term.

The continuity equation Eq. 4 in the new variable is:

ubξ+vbη=0,(16)

and now using the decomposed forms of ub and vb from Eq. 15, we have:

u1ξ+v1η=0,(17)

Note that u0 and v0 are defined in Eq. 15, and some derivatives vanish:

u0ξ=v0ξ=0,

The dimensionless x momentum Eq. 6 for new variables becomes:

η(u0u1ξ)+η(u1u1ξ)+η(v1u0η)+η(v1u1η=2P1ξη+Re13u1ξ2η+3u1η3,(18)

The dimensionless y momentum Eq. 6 for new variables becomes:

ξ(u0v1ξ)+ξ(u1v1ξ)+ξ(v1v1η=2P1ξη+Re13v1ξ3+3v1ξη2,(19)

The pressure term is eliminated by expending the terms

η(u0u1ξ),η(v1u1η),ξ(u0v1ξ).

in Eqs 18, 19, and then by subtracting Eq. 19 from Eq. 18, we get:

u0ηu1ξ+u02u1ξη+η(u1u1ξ+v1u1η)+v1ηu0η+v12u0η2u02v1ξ2ξ(u1v1ξ+v1v1η)}=Re13u1ξ2η+3u1η33v1ξ33v1ξη2,(20)

We define a stream function (ψb) such that:

ub=ψbηandvb=ψbξu0+u1=ηψ0+ψ1andv1=ψ1η.(21)

Substituting in Eq. 20, we get:

u0ξ2ψ12u0η2ψ1ξRe14ψ1=ξ(u1v1ξ+v1v1η)η(u1u1ξ+v1u1η),(22)
where2ψ1=2ψ1η2+2ψ1ξ2&4ψ1=4ψ1η4+24ψ1ξ2η2+4ψ1ξ4.

From the above continuity Eq. 22, we have:

u1ξ+v1η=0u1u1ξ+v1η=v1u1ξ+v1η=0.(23)

By substituting in Eq. 22, we have:

u0ξ2ψ12u0η2ψ1ξRe14ψ1=ξξu1v1+ηv1v1ηξu1u1+ηu1v1.(24)

The no slip boundary conditions:

u0yUξ+u1ξ,yUξ=0,v1ξ,yUξ=0,u0yLξ+u1ξ,yLξ=0,v1ξ,yLξ=0.(25)

Reconsidering the total stream function:

ψ0yUξ+ψ1ξ,yUξ=34.ψ0yLξ+ψ1ξ,yLξ=0.(26)

Floryan (1997) used the Fourier expansion for the simplification of the above equation:

ψ1ξ,η=n=n=ϕnηeinαξ,(27)

where ϕn=ϕn* represents the conditions, and the star identifies the complex conjugate.

By substituting Eq. 27 into Eq. 24, we get:

D4ϕ0+2αReImk=1k=kD2ϕk*Dϕk=0,(28)

where D=η,Im represents the imaginary part. Integrating Eq. 28 w.r.t. ‘η’ we get:

D3ϕ0+2αReImk=1k=kDϕk*Dϕk=ARe,whereA=p1ξmean(29)

The equations of the model are valid for a constant pressure gradient, where both the straight and sinusoidal channels are driven by the same gradient. The flow rate variation is taken into account and measured due to the extra flow losses that occur in the meandering channel. Additionally, it is assumed that in Eq. 26, A is equal to zero, and the correction for the flow rate, Q1, can be easily calculated as per Floryan’s work in 1997.

5 Solution of the modified equations

Let us assume that the wavelength is a large quantity, denoted as α0. We can then obtain a solution for the problem. To regularize the flow domain, transformation θ=η2Scosσ is applied, where σ=αξ represents a slow scale. With these new variables, the governing equations in dimensionless form, denoted as Eqs 4, 6, are transformed into the following set of field equations given as:

αubσ+2Sαsinσubθ+vbθ=0,(30)
1+4S2α2sin2σ2ubθ2α22ubσ24Sα2sinσ2ubθσ2Sα2cosσubθ+2SαResinσububθ+Revbubθ+αReububσ+αRepbσ+2SαResinσpbθ=0,(31)
1+4S2α2sin2σ2vbθ2α22vbσ24Sα2sinσ2vbθσ2Sα2cosσvbθ+2SαResinσubvbθ+Revbvbθ+αReubvbσ+Repbθ=0,(32)

The no-slip boundary and free-stream conditions with the total flow rate given in Eq. 25 is reduced to the following simplest form:

ub=vb=0atθ=±1,ub=1atθ=0,Qb=11ubσ,θdη=43.(33)

In Eqs 3032 are three unknown quantities, i.e., ub,vb, and Pb that we need to determine by using the perturbation method. These unknown quantities are expanded in terms of α by using the following series:

ub=u^0+αu^1+α2u^2+α3u^3+Oα4,vb=v^0+αv^1+α2v^2+α3v^3+Oα4,Pb=α1P^1+P^0+αP^1+α2P^2+Oα3.(34)

Substituting the values of ub,vb, and Pb from Eq. 34 into Eqs 3032, we have:

ασu^0+αu^1+α2u^2+α3u^3+2Sαsinσθu^0+αu^1+α2u^2+α3u^3+θv^0+αv^1+α2v^2+α3v^3=0,(35)
1+4S2α2sin2σ2θ2u^0+αu^1+α2u^2+α3u^3α22σ2u^0+αu^1+α2u^2+α3u^34Sα2sinσ2θσu^0+αu^1+α2u^2+α3u^32Sα2cosσθu^0+αu^1+α2u^2+α3u^3+2SαResinσu^0+αu^1+α2u^2+α3u^3θu^0+αu^1+α2u^2+α3u^3+Rev^0+αv^1+α2v^2+α3v^3θu^0+αu^1+α2u^2+α3u^3+αReu^0+αu^1+α2u^2+α3u^3σu^0+αu^1+α2u^2+α3u^3+αReσα1P^1+P^0+αP^1+α2P^2+2SαResinσθα1P^1+P^0+αP^1+α2P^2=0,(36)
1+4S2α2sin2σ2θ2αv^1+α2v^2+α3v^3α22σ2αv^1+α2v^2+α3v^34Sα2sinσ2θσαv^1+α2v^2+α3v^32Sα2cosσθαv^1+α2v^2+α3v^3+2SαResinσu^0+αu^1+α2u^2+α3u^3θαv^1+α2v^2+α3v^3+Rev^0+αv^1+α2v^2+α3v^3θαv^1+α2v^2+α3v^3+αReu^0+αu^1+α2u^2+α3u^3σαv^1+α2v^2+α3v^3+Reθα1P^1+P^0+αP^1+α2P^2=0.(37)

6 Analysis of heat transfer in a meandering channel

The analysis of heat transfer in a meandering channel involves studying the heat transfer mechanisms that occur due to meanders in the channel. Heat transfer occurs in three modes: conduction, convection, and radiation. In the present study, a case of heat transfer through conduction and convection is taken into account. Conduction occurs when heat flows from regions of high temperature to regions of low temperature through a solid material. In a meandering channel, heat is conducted through the walls of the channel. The rate of heat transfer through conduction is proportional to the temperature gradient, the thermal conductivity of the material, and the cross-sectional area of the channel. Convection occurs when heat is transferred by the movement of fluids, either liquids or gases. In a meandering channel, heat is transferred through convection due to the flow of fluid through the channel. The rate of heat transfer through convection is proportional to the temperature difference between the fluid and the channel wall, the velocity of the fluid, and the heat transfer coefficient of the fluid. The heat transfer is analyzed by using the energy equation, which relates the rate of heat transfer to the temperature distribution in the channel. The energy equation is taken into account to study the conduction and convection mechanisms of heat transfer and the thermal properties of the meandering channel in the flow of fluids.

7 Formulation of the problem

The problem is modeled using the Navier–Stokes equations, along with the continuity and energy equations. To obtain a non-linear, simple partial differential equation, a defined set of transformations is applied. The regular perturbation technique is utilized to expand the heat transfer, Tb, in terms of " α " in the meandering channel.

The energy equation is:

ρcpDTDt=k2T+Φ,(38)

where T is the temperature distribution, ρ is density, cp is specific heat, k is thermal conductivity, and Φ is the dissipation term. The boundaries condition of the problem is at the upper wall, T=T1, and at the lower wall ,T=T0.

Consider the following dimensionless variables in order to transform the energy equation into dimensionless form:

T*=TT0T1T0,Φ*=L2ΦU2μ,(39)

The dimensionless numbers that appear in Eq. 38 are the Reynolds and Prandtl numbers, which are defined as:

Re=ρULμ,Pr=μcpk,Re.Pr=ρULcpk,1RePr=kρULcp,(40)
u*T*x*+v*T*y*=1RePr2T*x*2+2T*y*2+EcReΦ*,(41)

Equation 41 is the dimensionless form of the energy equation. On the other hand, the boundaries conditions for the temperature distribution in the dimensionless form are obtained, and the dimensional boundary conditions are given by taking:

T*x,y=TT0T1T0,(42)

The temperature at the upper wall is T=T1, so using definition from Eq. 42:

y=1,T=T1,T*x,y=TT0T1T0=T1T0T1T0,T*x,y=1,(43)

The temperature at the lower wall is T=T0, so using definition from Eq. 42:

y=1,T=T0,T*x,y=TT0T1T0=T0T0T1T0,T*x,y=0,(44)

The temperature of a flow in a straight channel is denoted by T^0, and so we have:

u0T^0x+v0T^0y=1RePr2T^0x2+2T^0y2+EcRe2u0x2+2v0y2+v0x+u0y2,(45)

The boundaries condition for the temperature at the upper wall, T^0=1, and at the lower wall, T^0=0.

Putting Eq. 11 into Eq. 45, after simplification we have:

2T^0y2=4PrEcy2,(46)

Integrating Eq. 46 twice w.r.t y, we get:

T^0y=PrEcy43+C1y+C2,(47)

where C1 and C2 are the constant of integration, which can be determined by using the boundary conditions imposed at the upper wall of the channel, T^01=1, and at lower wall of the channel, T^01=0.

T^0y=PrEc31y4+121+y,(48)

As the VbX=V0X+V1X and TbX=T0X+T1X, where VbX=ubi^+vbj^ and TbX are the velocity and temperature of the total flow, respectively, and V1X=u1i^+v1j^ and T1X are the velocity and temperature induced by the channel geometry, respectively.

It is assumed that α is a small quantity, and for this choice of α we expand the unknown quantities of T^b by using the regular perturbations to find T^b. The dimensionless form of the energy equation is transformed and expanded using the following transformation and expansion:

η=y2scasξ,ξ=αx,(49)
T^b=T^0+αT^1+α2T^2+Oα3,(50)

We get:

u02SαsinξT^0η+αT^1η+α2T^2η+αT^0ξ+αT^1ξ+α2T^2ξ+α2u22SαsinξT^0η+αT^1η+α2T^2η+αT^0ξ+αT^1ξ+α2T^2ξ+αv1T^0η+αT^1η+α2T^2η=1RePr[4S2α2sin2ξ2T^0η2+α2T^1η2+α22T^2η2+α22T^0ξ2+α2T^1ξ21+α22T^2ξ2+4α2Ssinξ2T^0ηξ+α2T^1ηξ+α22T^2ηξ+2Sα2cosξT^0η+αT^1η+α2T^2η+2T^0η2+α2T^1η2+α22T^2η2]+EcRe[u0η2+α2(8S2sin2ξu0η2+2u0ηu2η+4Ssinξu0ηv1η+2v1η2+8Ssinξu0ηu0ξ+2u0ξ+2u0ηv1ξ)+Oα3],(51)

The boundaries condition on the upper wall is y=h=1, and the boundaries condition on the lower wall is y=h=1 for the small α.

T^01=1,T^11=0,T^21=0,(52)
T^01=0,T^11=0,T^21=0,(53)

8 Results and discussion

With the solution represented as a power expansion in terms of σ and by equating the coefficients, we have:

u^0θ=1θ2,u^1σ,θ=0,u^2σ,θ=43Sθcosσ1θ2,
u^3σ,θ=2315ReSθsinσθ21θ4+21Sθθ24cosσ+θ234,
v^0σ,θ=0,v^1=2Ssinσ1θ2,v^2σ,θ=0,
v^3σ,θ=13Sθ21θ21+8Sθcosσsinσ,
P^1σ=2σRe+C0,P^0=C1,
P^1σ,θ=4Re1S2SσθsinσSsin2σ,
P^2σ,θ=2Scosσθ2θ33+θ55+1115S2cos2σ,
T^0η=PrEc31η4+121+η,
T^1η=11260sinξPrRe1+S1+η21055+η2+8EcPrη1111η2+10η4,

And

T^2=11260[1052pr2Re21+Scosξη4+1+1176EcPrSη5+η)cosξ+4EcPr3Re21+Sη7ηcosξ+109EcPr3Re21+Sη9+ηcosξ+88EcPr3Re21+S2η21sin2ξ+420Pr2Re21+S2η3ηsin2ξ+43EcPr11Pr2Re21+S2+3780S2η4+1sin2ξ+56EcPr3Re21+S2η6+1sin2ξ+20Pr2Re21+S2η7ηsin2ξ+50EcPr3Re21+S2η81sin2ξ+1129EcPr3Re21+S2η10+1sin2ξ+84PrRe1+Sη5+ηsinξScosξ+2PrRe1+Ssinξ+72Pr2Re1+Sη61cosξRe32EcSsinξ+43Prη3+ηcosξEc11Pr2Re21+S2100S+  630Re1+SSsinξ+  12η2+1cosξ1055Pr2Re21+S+24S352EcPr2Re1+SSsinξ+  20EcPr2Re1+SSη8+1sin2ξ].

Stream function ψσ,θ

The stream function ψσ,θ is used to define the velocity components such that:

ub=ψσ,θθ,vb=2αSsinσψσ,θθ+ψσ,θσ,(54)

By integrating ub with ‘θ’ and after finding the constant of integration, we get:

ψσ,θ=11260(420Sα21+θ22cosσ+reSα31+θ2267+2θ2+θ4sinσ+41+θ21052+θ+  reS2α31020θ+30θ212θ36θ4+3θ5sin2σ).(55)

The streamlines in the area adjacent to the lower wall are shown in Figures 2A, B, where S=5 amplitude), (Reynolds number) =1 , α=10, while Figures 2C is plotted for the whole channel, which determines that a diminution of the meandering wavelength results in the establishment of separation zones in the troughs. Figure 3 illustrates the variations of the pressure losses as a function of S and α. From Figure 2, the area adjacent to the lower wall is complex and exhibits a variety of behaviors. The fluid experiences both pressure and shear forces. The pressure forces result from the changes in the height of the channel, while the shear forces arise due to the velocity gradient near the channel walls. As a result of these forces, the streamlines in the area adjacent to the lower wall exhibit the behaviors of separation and reattachment. As the fluid flows over the wave crest, the pressure forces cause the flow to separate from the lower wall, resulting in a region of recirculating flow or eddies. As the fluid flows over the wave trough, the pressure forces cause the flow to reattach to the lower wall, resulting in a region of reversed flow. In addition to the primary flow along the channel axis, the pressure and shear forces are also induced the secondary flows perpendicular to the channel axis. These secondary flows create complex patterns of vortices and recirculation zones, as shown in Figure 2, which affect the overall flow behavior. The complex flow behavior in the area adjacent to the wall also leads to flow instability, such as the formation of Kelvin–Helmholtz instability waves, which can cause the flow to break down into smaller-scale vortices. Overall, the streamlines in the area adjacent to the lower wall of a meandering channel in a viscous fluid flow exhibit a variety of complex behaviors, which have an important implication for fluid transport and mixing in microfluidic and other applications.

FIGURE 2
www.frontiersin.org

FIGURE 2. The streamlines in the area adjacent to the lower wall are shown in (A,B) where S = 5 (amplitude), Re = 1, and α = 10, while (C) is plotted for the whole channel.

FIGURE 3
www.frontiersin.org

FIGURE 3. Variations of the pressure gradient correction Re[(dP_1)/dx] = -0.002.

9 Conclusion

The periodicity of the meandering channel affects the flow pattern and turbulence intensity of the fluid. A higher periodicity leads to a more regular flow pattern, while a lower periodicity leads to a more chaotic flow pattern. The presence of periodicity also leads to the development of secondary flows, such as Dean vortices, which affect the mixing and heat transfer characteristics of the flow. The Prandtl number relates the momentum diffusivity to the thermal diffusivity of a fluid. In a meandering channel, a higher Prandtl number results in a thicker thermal boundary layer, which affects the heat transfer characteristics of the flow. A lower Prandtl number, on the other hand, results in a thinner thermal boundary layer and a more efficient heat transfer. As shown in the Figures 48, the Eckert number relates the kinetic energy of a fluid to its thermal energy. In a meandering channel, a higher Eckert number results in a more energetic flow, which can lead to an increase in turbulence intensity and mixing. A lower Eckert number results in a less energetic flow, which leads to a more laminar flow pattern and reduced mixing. The Reynolds number represents the inertial forces to the viscous forces in a fluid. In a meandering channel, a higher Reynolds number results in a more turbulent flow pattern, with increased mixing and heat transfer. A lower Reynolds number results in a more laminar flow pattern, with reduced mixing and heat transfer. The effects of periodicity and the Prandtl number, Eckert number, and Reynolds number on fluid flow in a meandering channel are complex and interrelated.

FIGURE 4
www.frontiersin.org

FIGURE 4. Temperature distribution is uniform and smooth from 0 to 1 in the direction of η, whereas its behavior in ξ direction is periodic in nature. For fixed values of S = 0.0125, Re = 10, α = 0.1, Pr = 1, and Ec = 1.

FIGURE 5
www.frontiersin.org

FIGURE 5. Temperature distribution Tη,ξ is plotted against η for fixed values of ξ=0, S=0.0125, Re=1, Pr=2, and Ec=0.0001 and different values of α. Note that the Prandtl and Eckert numbers are the choice for water. The temperature profile is increasing with α, whereas each profile changes uniformly. Moreover, the profile is either linear or close to linear for small values of α.

FIGURE 6
www.frontiersin.org

FIGURE 6. Temperature distribution Tη,ξ is plotted against η for fixed values of ξ=0, S=0.0125, Re=1, α=0.25, Br=0.01, and different values of Pr. The range of Pr in water is 2Pr7, and Br=PrEc, the different values for the Prandtl numbers is Pr, the Eckert number becomes Ec=0.0085,0.0056,0.00425,0.0034. Note that the Prandtl and Eckert numbers are the choice for water. The temperature profile is increasing in Pr, whereas each profile changes uniformly. Moreover, the profile is either linear or close to linear for small values of Pr.

FIGURE 7
www.frontiersin.org

FIGURE 7. Temperature distribution Tη,ξ is plotted against η for fixed values of ξ=0, α=0.5, Re=1, Pr=2, and Ec=0.0001 and different values of S. Note that the Prandtl and Eckert numbers are the choice for water. The temperature profile is increased in S, whereas each profile changes uniformly. Moreover, the profile is either linear or close to linear for small values of S.

FIGURE 8
www.frontiersin.org

FIGURE 8. Temperature distribution Tη,ξ is plotted against η for fixed values of ξ=0, α=0.5, S=0.125, Pr=2, and Ec=0.000, and different values of Re. Note that the Prandtl and Eckert numbers are the choice for water. The temperature profile is increased in Re, whereas each profile changes uniformly. Moreover, the profile is either linear or close to linear for small values of Re.

In future, the flow and heat transfer in wavy meandering channels could include investigating the thermal profiles, skin friction, and Nusselt number calculations for various flow parameters. This could provide insights into the heat transfer characteristics of the flow and help optimize the design of such channels for specific applications. Additionally, examining the effects of different geometries and materials on the flow and heat transfer could yield valuable results.

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

DK proposed the mathematical model. SI performed all numerical computations and created the corresponding graphs. DK and SI provided the discussion of the graphs and their physical interpretation. NU conducted the literature review and compared the present simulations with classical data. KN and K carried out the final review and made amendments to the manuscript.

Funding

This study was supported by 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

Alsaedi, A., Rashad, A. M., and Hayat, T. (2016). Bioconvective Darcy–Forchheimer flow of the Ree–Eyring nanofluid through a stretching sheet with velocity and thermal slips. J. Mol. Liq. 220, 803–811.

Google Scholar

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

Google Scholar

Fiebig, Martin, and Chen, Yuwen (1999). Heat transfer enhancement by wing-type longitudinal vortex generators and their application to finned oval tube heat exchanger elements. Heat Transf. Enhanc. Heat Exch. 42, 79–105.

CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

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

CrossRef Full Text | Google Scholar

Ghalambaz, M., Rostami, B., and Rashidi, M. M. (2016). Thermal analysis of a radiative nanofluid over a stretching/shrinking cylinder with viscous dissipation. J. Mol. Liq. 224, 87–93.

Google Scholar

Hayat, T., Ahmad, B., Alsaedi, A., and Rasool, N. (2017). Theoretical analysis of the thermal characteristics of Ree–Eyring nanofluid flowing past a stretching sheet due to bioconvection. J. Mol. Liq. 242, 758–764.

Google Scholar

Izadi, M., and Pourmehran, O. (2017). Formation of hydrogen bonding during the Homann flow over a cylindrical disk of variable visco-elastic nano-materials in presence of activation energy. Int. J. Hydrogen Energy 42 (28), 18003–18016.

Google Scholar

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 (3), 295–309. doi:10.1016/0894-1777(95)00066-u

CrossRef Full Text | Google Scholar

Jensen, Michael K., Bergles, Arthur E., and Shome., Biswadip (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

CrossRef Full Text | Google Scholar

Khan, M., Saleem, S., and Rehman, K. U. (2017). Energy and mass transport through hybrid nanofluid flow passing over an extended cylinder with the magnetic dipole using a computational approach. J. Magnetism Magnetic Mater. 427, 107–114.

Google Scholar

Khan, W. A., and Ahmed, S. (2015). Significance of activation energy and entropy optimization in radiative stagnation point flow of nanofluid with cross-diffusion and viscous dissipation. Results Phys. 5, 301–311.

Google Scholar

Ligrani, Phil M., Oliveira, Mauro M., and Tim Blaskovich, (2003). Comparison of heat transfer augmentation techniques. AIAA J. 41 (3), 337–362. doi:10.2514/2.1964

CrossRef Full Text | Google Scholar

Mahmud, S., and Uddin, M. J. (2021). Numerical simulations for optimised flow of second-grade nanofluid due to rotating disk with nonlinear thermal radiation: Chebyshev spectral collocation method analysis J. Therm. Analysis Calorim. 145 (2), 1087–1100.

Google Scholar

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

CrossRef Full Text | Google Scholar

Nadeem, S., Akbar, N. S., and Lee, C. (2016) Numerical analysis of a time-dependent aligned MHD boundary layer flow of a hybrid nanofluid over a porous radiated stretching/shrinking surface J. Mol. Liq. 220, 945–954.

Google Scholar

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

CrossRef Full Text | Google Scholar

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

Google Scholar

Schlichting, H., and Kestin, J. (1961). Boundary layer theory. New York, NY, USA: McGraw-Hill.

Google Scholar

Sheikholeslami, M., Rokni, H. B., and Ganji, D. D. (2017). Numerical investigation of Darcy–Forchheimer hybrid nanofluid flow with energy transfer over a spinning fluctuating disk under the influence of chemical reaction and heat source. J. Therm. Analysis Calorim. 128 (2), 913–924.

Google Scholar

Srinivasacharya, D., and Sibanda, P. (2020). Investigation of hydromagnetic bioconvection flow of Oldroyd-B nanofluid past a porous stretching surface. J. Therm. Analysis Calorim. 139 (1), 237–247.

Google Scholar

Taylor, Geoffrey Ingram. "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. (1923).

Google Scholar

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, IA, USA: Iowa State University College of Engineering.

Google Scholar

White, F. M., and Majdalani, J. (2006). Viscous fluid flow. New York, NY, USA: McGraw-Hill, 433–434.

Google Scholar

Xu, Minghai, 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

CrossRef Full Text | Google Scholar

Keywords: stream wise vortices, meandering channel, flow separation, instability, nonuniform radius

Citation: Ibrahim S, Khan Marwat DN, Ullah N, Nisar KS and Kamran (2023) Flow and heat transfer in a meandering channel. Front. Mater. 10:1183175. doi: 10.3389/fmats.2023.1183175

Received: 09 March 2023; Accepted: 29 May 2023;
Published: 26 July 2023.

Edited by:

Noor Saeed Khan, University of Education Lahore, Pakistan

Reviewed by:

M. Riaz Khan, Quaid-i-Azam University, Pakistan
Nilankush Acharya, Jadavpur University, India

Copyright © 2023 Ibrahim, Khan Marwat, Ullah, Nisar and Kamran. 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

Disclaimer: 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.