• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    Hydrodynamics and Heat Transfer Analysis of Airflow in a Sinusoidally Curved Channel

    2022-08-23 02:17:28AbidMemonAsifMemonKaleemullahBhattiThaninSitthiwiratthamandNichaphatPatanarapeelert
    Computers Materials&Continua 2022年6期

    Abid.A.Memon,M.Asif Memon,Kaleemullah Bhatti,Thanin Sitthiwiratthamand Nichaphat Patanarapeelert

    1Department of Mathematics and Social Sciences,Sukkur IBA University,Sukkur,65200,Sindh,Pakistan

    2Mathematics Department,Faculty of Science and Technology,Suan Dusit University,Bangkok,10300,Thailand

    3Department of Mathematics,Faculty of Applied Science,King Mongkut’s University of Technology North Bangkok,Bangkok,10800,Thailand

    Abstract:For heat transfer enhancement in heat exchangers,different types of channels are often tested.The performance of heat exchangers can be made better by considering geometry composed of sinusoidally curved walls.This research studies the modeling and simulation of airflow through a 2π units long sinusoidally curved wavy channel.For the purpose, two-dimensional Navier Stokes equations along with heat equations are under consideration.To simulate the fluid flow problem,the finite element-based software COMSOL Multiphysics 5.4 is used.The parametric study for Reynolds number from Re = 100 to Re = 1000 and the period of vibration P from 0 to 5 are observed.The surface plots, streamline patterns, contours, and graphs are presented for the velocity field magnitude, temperature, and pressure against the Reynolds number as well as period of vibration.The results are compared with various literature.It is found that due to the creation of periodic contraction regions the velocity magnitude of the flow is continuously increasing with the increase of Reynolds number,on the contrary the pressure is decreasing from inlet to outlet of the channel.Also,a periodic variation in the pressure distribution along the vibrating boundaries has been found with an average increase of 500%for the high Reynolds number.A novel work was done by expressing the rotation rate per second in terms of local Reynolds number for the recirculating regions found due to the periodic oscillation of the boundaries.The average temperature near the outlet where a fixed temperature is imposed initially is decreasing with an increase in Reynolds number.The convection process is weakened due to an increase of periodic vibration of boundaries.

    Keywords: Simulation of fluid flow; heat transfer; cosine function; wavy channel

    Nomenclature

    1 Introduction

    The wavy channels in which the contraction and expansion regions are created periodically is very essential in most of the devices for enhancing the heat distribution,pressure reduction,and maximizing the flow rate.Because of having the expansion and the contraction regions,the fluid passage through these channels attains an automatically created periodicity.Having the periodic variation in the crosssectional area,the fluid gains a regime of fully developed that is quite different from those channels having the constant cross-sectional area.Fluid flow problems in the wavy channels are often studied to examine the separated flows in which separated vortices are formed in the expansion regions,patterns of the flow velocity,pressure patterns,and the temperature distribution due to the periodic variation.Due to the periodic variation in the parameters of the geometry,the pattern of the fluid movement is largely affected by the flow velocities,loss in pressure,and heat enhancements[1–4].The wavy channels with the corrugated walls having the periodic expansion and contraction are widely studied in the many of the application like blood oxygenators[5,6],in the microfiltration cell[7],the wavy channels are also studied as the vortex generators to increase the thermal performance in solar panels[8,9].There is a lot of work focused upon the fluid flow through sinusoidally curved wavy channels,details are given below:

    In[10],the authors studied the laminar flow and thermal characteristics in the cavity composed of sinusoidal functions.It was discussed that due to the increase in the flow velocities, a large number of vortices are produced inside the cavities which might increase hydrodynamic and thermal enhancements in the channels.On examining the wavy channel [11] reported that when compared with the rectangular channel,the wavy types channel has the largest values of the friction factor and because of the reason the wavy channels are largely used for the cooling process.Investigating the flow pattern through the wavy channels [12] composed of sinusoidal functions with the parametric study for Reynolds number, it was reported that in the range of Re>700, a turbulent flow can be assumed in the wavy channel due to the creation of unsteady vortices in the expansion regions.An experimental investigation was performed[13–15]to investigate the laminar as well as the turbulent fluid flows and heat transfer with the use of a wide range of Reynolds numbers.It was reported that due to the presence of recirculation vortices in the regions,mixing of the fluid takes place quickly near the boundaries that will support enhancing the heat transfer in a wavy channel.Due to the curvy walls of the channel,the thermal boundary layers became destabilize and promote turbulence in the channel which might favor increasing the heat transfer.On investigating numerically[16]fluid flow and heat transfer through the sinusoidal wavy channel it was presented that the local Nusselt number possesses a direct relationship with the Reynolds number change periodically with the flow of direction.A twodimensional curved channel with a fully developed flow was studied[17]numerically along with the heat transfer.The research focused on the geometrical influence of Reynolds number,Prandtl number on flow distribution,and heat transfer.It was showed that with the increase of wall-to-wall distance and amplitude of the vibration of the wave, the heat transfer coefficient, streamlines contours are significantly affected.It was also endorsed that using the wavy geometry of any shape the pressure drop will be higher than that of a rectangular channel for the same dynamically and thermal parameters.The wavy channels are under consideration in many microfluidics devices for low Reynolds number.In many applications of[18–20]microfluidics it was stated that it is the only diffusion that is a mode of transport that took place in the wavy channels for low Reynolds number.Hence, to enhance the mixing in micro-channels it is essential to put Re<1 as described in [21,22].In the presence of a magnetic field, a slip flow was studied [23] for the fluid flow and heat transfer problem where the governing partial differential equation first converted in to fourth order differential equation and then a numerical solution obtained by Adomain decomposition.The two models for the hybrid nano-fluid were compared for a rotating channel with porous medium under the surface catalytic reaction[24].The single and multi-wall nanotubes were taken with the ethylene glycol as the base fluid.Another study was done for the purpose to compare the consequences heating with Newtonian impacts and homogeneous as well as heterogeneous reaction for the flow of water-silver nano-fluids [25].To see the better impacts of the study the thermal radiation and magnetic field were also added as the feature.Later the colloidal suspension with water for aluminum oxide and gamma-aluminum oxide was studied[26]in an elastic cylinder.The research becomes the beneficial due to the application of the thermal transport in the side of industries,aerodynamics and medical sciences.

    Besides all these works there is much effort devoted by researchers in last few years in modeling nanofluids flows through different curvy geometries with application in energy storage,drug delivery,cancer treatment, tissue generation, readers are referred to the works [27–29].In this paper, the hydrodynamics of airflow and heat transfer through a sinusoidally curved wavy channel is studied in detail.Section 2 presents the problem description.In Section 3,we present different mesh sizes and correlation analysis.The variation in flow variables including temperature distribution is observed by changing Reynolds number from 100 to 1000 and period of vibration of walls from 0 to 5.The problem is modeled in COMSOL Multiphysics 5.4 and solved using finite element method to obtain flow variables including temperature profile.Section 4 precisely discusses the results and presents a comparison analysis with the recently published work.A detailed parametric study using graphs is also carried out.In Section 5 we give the concluding remarks.

    2 Methodology

    2.1 Geometrical Construction of the Wavy Channel

    The current problem of the fluid flow and heat transfer is tested by making the open channel who’s upper and lower boundaries are shaped up by using the cosine function see Fig.1.Such geometries are creating the auto oscillation or vibration in a domain that generates hypes or contraction region in the channel which results increase the speed of the particles and reducing the dynamic pressure.The length of the channel is fixed of whereas the boundaries of the channel are vibrated with the period P.Here the boundary below of the channel is erected with the help of function and the lower boundary of the channel is erected by the function.WhereA= 0.3 cm is known as the amplitude of the vibration of the channel.Pis the period of the vibration of the channel.For the contemporary simulation, the periods of the vibration are tested from 0 to 5.By puttingP= 0, corresponds to the straight rectangular channel which is normally tested in many of the fluid flow problems.The total height of whereHis the length of the contracted region.An average inlet velocity will be forced from the left entrance of the channel while the right end of the channel is perceived as an outlet of the channel,where the pressure is getting as zero.The walls of the upper and lower walls of the channel are controlled by the slip boundary conditions.The common air as the testing fluid will be allowed to enter at the initial temperature of and the last vertical block of the channel will be kept at the temperature ofTheat=323.15K.

    Figure 1:Geometrical construction and boundary conditions of the problem

    2.2 Problem Formulation

    The modeling and simulation for the present fluid flow and heat transfer problem are examined by considering the steady-state, two-dimensional, and incompressible Newtonian fluid.We use COMSOL Multiphysics 5.4 for the simulation which implements streamline-upwind Petrov-Galerkin’s finite element method[15]to solve the coupling system of Navier Stokes equations and heat equation.In Fig.3, we present the work flow of the problem.Let=(u,v)denotes the velocity field of the fluid in two dimensions,whereuandvare respectively the horizontal and the vertical components of velocity field.The governing field equations are:

    wherepdenotes pressure and Re denotes the Reynolds number.Eqs.(1) and (2) can be written in component form as:

    The heat transfer for the fluids in the vector form and in component form can be written as by Eqs.(5)and(6)

    Following boundary conditions are set:

    (a) Slip boundary condition at upper and lower walls

    (b) Fully developed flow conditions at inlet and outlet of the channel are respectively

    and

    wherep0andTheatare the pressure at temperature applied at the outlet of the channel.Other computational parameters are:

    where Relocalis local Reynolds number,Nuxis the Nusselt number andthe averaged Nusselt number.

    3 Validation and Comparison Between Different Mesh Sizes

    3.1 Mesh Independent Study

    It is the describable quality of the finite element procedures to get the results up to the benchmark by performing a mesh independent study for the fluid flow problem and heat transfer.To approve the current fluid flow problem for a sinusoidally curved wavy channel a mesh independent study is performed along the channel.For this reason,five different types of meshes are practiced to estimate the ratio of the centerline horizontal component of velocityuto the inlet velocityuinat downstream of the channel.In Figs.2a–2e,five different meshes with different qualities are observed.Each of them is used to measure the ratiou/uinat the lower boundary of the channel see Fig.3.From the graph,it is very much explicit that with the implementation of fine to extremely fine meshes the problem of the flow is mesh-independent i.e.,the results achieved through the extremely fine meshes are conferring more accuracy to match the benchmark problem when studied with such type of the channel.

    Figure 2:Triangular meshing(a)normal(b)fine(c)finer(d)Extra fine and(e)Extremely fine

    Figure 3:The ratio uuin at the lower boundary against the x-axis for a different number of elements

    Statistics for extremely fine mesh as well different other meshes are given by the following Tab.1.

    Table 1: Mesh statistics

    3.2 Association with Correlation

    After doing the mesh independent study,we are adjusting our intention to approve the results for the local Nusselt number downstream of the channel by exploring the correlation available through different pieces of literature.According to the literature,the local Nusselt number for a flat plate in the state of forced convection is provided by the equation correlating the local Reynolds number and Prandtl number [7].However, our channel of concern is only flat for the periodP= 0 and gives satisfactory results when contrasted with the Nusselt number given by literature and the present work see Figs.4a–4c.Furthermore,the results are also adequate when we compare the results forP=3 see Figs.4d–4f.

    Figure 4:Correlation-based confirmation of the model for the local Nusselt number at the downstream of the channel at P=0(a)-(c),and P=3(d)-(f)

    4 Results and Discussion

    The laminar and Newtonian fluid flow onward with heat transfer was tested with the use of the Reynolds number Re from 100 to 278 K for a wavy channel with arcs whose boundaries are created by using the cosine functions.The boundaries are vibrated with the fixed amplitude of vibrationA=0.3 cm and with a different number of periodsP=0 toP=5.Designing such types of geometries based on the contraction and expansion regions with equal space of areas.Such type of geometries will create the automatic periodic flows along the channel and will be in favor of increasing velocities and decreasing pressures in the domain.Besides,for the current convection heat transfer problem,we are working to test the heat distribution due to the fluctuation of the boundaries at the exit of the channel.For this purpose, the air as a testing fluid is allowed to enter from the inlet of the domain with the temperature of 278 K,and the outlet of the wavy channel is fixed with the temperature of 323.15 K.The whole simulation was preparing on the finite element based commercial software of COMSOL MultiphysicsP.The outcomes are revealed with the surface plots,streamlines,and the contours for the velocity field,pressure,and temperature.The graphs are relocated to show various impacts of the Reynolds number as well as the periods of the vibrationPon the velocity field magnitude,pressure,and temperature near the outlet of the channel as well as for the contraction regions in the domain.

    4.1 Surface Plots of the Velocity Field with Streamlines and Axial Velocity

    The streamlines pattern with the surface plots is dispatched in Fig.5.For the straight rectangular channel, the flow of particles in the domain remains consistent i.e., there was no increment or declination in fluids velocities.Moreover, there was no formation of vortices in the domain.Next,we check the pattern of the velocity field in the vibrating channel by putting a period ofP=1,which produces only one contracted region of the flow.On passing the fluid in the region,the total volume of the fluid is reduced in the region which might increase the flow rate of the particles.With the increase in periods P,the number of contraction regions is increased where the maximum velocity can be deducted.Fluid particles after crossing the contraction region sliding down in the curved cavities where the recirculation will be produced.Due to the recirculation, the speed of the fluid particles reduced and then optimized after reaching the contraction region.In simple words,we would say the maximum velocity magnitude of the flow can be seen in the middle of the channel in terms of periods.Although with enhancing the Reynolds number a chaotic flow can be seen at the middle of the channel for example see forP= 5 at Re = 500 and 1000.The area of the chaotic region is increased with the increase in Reynolds number.Near the contraction region,the velocity magnitudes are calculated for theP= 5 at the particular Reynolds number see Figs.6a–6c.ForP= 5 the contraction regions are created from 0.5 to 5.5 cm with intervals of 0.5 cm.It is evident from Figs.7a–7c,the first contraction region atx= 0.5 cm has the normal velocity of the particles.Later, when the flow is passing other regions of the contraction the flow velocity is considerably increasing see Figs.6b and 6c.It is obvious that for larger the Reynolds number Re the increment in maximum flow velocity magnitude is larger when observing these regions of the contraction.

    Figure 5:Streamlines of the velocity field with surface plots for Re = 100,Re = 500,and Re = 1000(a)P=0(b)P=3(c)P=4(d)P=5(e)P=6

    Figure 6:Speed of the fluid particles near the contraction region(a)Re = 100(b)Re = 500,and(c)Re=1000

    Figure 7:Maximum flow rate at the outlet of the channel against Re

    In Fig.7,the maximum velocity magnitude at the outlet of the channel is calculated against the Reynolds number at different values of the periodPof the channel.It is obvious from the graph that constructing a channel with vibration or oscillating boundaries will produce the maximum flow rate at the outlet of the channel.Further,with the increase in Reynolds number,the maximum velocity at the outlet possessed a linear relationship with the Reynolds number.Except for the rectangular channel,all periods of the vibration of the channel are producing unique maximum velocity at the outlet of the channel for the particular Reynolds number.

    Applying a linear Regression procedure, we would make the correlation that will connect the maximum velocity at the outlet of the channel in terms of Reynolds number.Following equation is valid for the period of vibrationP >0:

    Fig.8 shows the computation of the axial velocity at the outlet of the channel for Re = 100 and Re = 300 for the period of the vibration.The graph is also produced to see the pattern of the flow at the exit of the channel.The graph can be compared with the pattern of the flow indicated in the literature[23]where the problem of the fluid flow and heat transfer was conducted with the channel consists of two period of vibration.The current problem has been extended for theP=3-P=5.In the current article,we are going to disclose the pattern of the axial velocity component at the outlet of the channel with the increase of periods of the vibration in the channel as well as Reynolds number.Mostly and actually,it is known that the construction of the wavy channel is most important in the applications where the fluid acquires speed automatically due to zig-zag motion.The more vibration given to the upper and lower boundary of the channel fluid gets stronger in terms of speed.Figs.9a–9c is showing the pattern of the axial velocity at the outlet of the channel.It can be distinguished that due to an increase in the value ofPfor the same value of Reynolds number,the axial velocity at the outlet of the channel increases significantly.Moreover, it is further predicted that with increasing value of Reynolds number by keeping same value ofPthe velocity is increasing but for some values ofPit is showing different patterns,see Figs.9a and 9c.

    Figure 8:Axial velocity profile at the outlet for the period P = 2 of the vibration of the channel and for(a)Re=100 and(b)Re=300

    Figure 9:Axial velocity profile at the outlet for the indicated period of vibration of the channel for(a)Re=100(b)2.5 and(c)Re=1000

    4.2 Pressure

    Pressure contours are obtained in Fig.11 using two different values of the Reynolds number,Re=500 and Re=1000 together with different values ofP.The contour presentation intimates that creating such geometry can be used as pressure-reducing pumps.The geometry is created through the special cosine function creating the expansion and the region of contraction due to the vibration of the boundaries.In the contraction region,the minimum pressure can be observed see Figs.10d–10f.Although with the rise in Reynolds number the pressure at the contracting channel reduced sufficiently.In the expansion region where the small vortices are creating which sustain the pressure to largely minimize.

    Figure 10:Pressure presentation in the domain with the periods of the vibration P=1(a and d),P=3(b and e),and P=5(c and f)

    Figure 11:Estimated results of the pressure(Pa)near all contracted regions for P=5 at(a)Re=100(b)Re=500 and(c)Re=1000

    In Figs.11a–11c,we have bestowed the estimation of dynamic pressure near the contracted regions against the y-axis for a particular Reynolds number with the only observation of periodsP= 5.It is clear from the graph as we go forward in the domain for the observation of the pressure after each contracted region, we found that the pressure is reducing continuously.With the extension of more periodsPof the vibration the pressure further can be subdued.Although,the increase in inlet velocities at the inlet of the channel supports the more loss of the pressure at the contracted regions.

    Pressure drop at the outlet of the channel is the function of Reynolds number as well as that of the period of vibration.By increasing the values ofPor Re the minimum pressure at the outlet is decreasing in a non-linear way as shown in Fig.12.

    Mostly in the industrial area,the wavy/curvy channels are constructed due to the need of optimum pressure in the region of interest.When the elements of any fluid passing through the contracted region the volumes are reduced which results a significant increase in the pressure.In Figs.13a–13c to 15a–15c we describe the pressure distribution along the length of the channel.The pattern can be compared with the patterns reported in the references[24]and[25].It is obvious from Figs.13a–13c while the fluid came into contact at the contracted regions due to the reduction of volume the pressure increases at the region.The pressure in the channel shows a periodic variation in the channel, see Figs.15a–15c.

    Figure 12:Minimum pressure at the outlet against the Reynolds number

    Figure 13: Pressure distribution along the length of the curvy channel with the period P = 2 (a)Re=100,(b)Re=500 and(c)Re=1000

    Also,with the increase in Reynolds number the range of periodicity increases for example Fig.14a the range of the amplitude for pressure can be seen at 2.5.While observing Fig.14c we can see the range of amplitude for the pressure is 60.It is also concluded that due to the increase in the periodic vibration of the boundaries of the channel an insignificant increase in the pressure occurs for moderate Reynolds and a significant increase for high Reynolds number.In our problem, the optimum pressures occur for Re=1000 with vibration periods ofP=5.It has been already cleared that and also described in many articles that adding the period of vibration at the boundaries of the channel producing maximum and minimum pressure at a time due to the oscillatory pattern of the pressure.Although the present problem has been simulated with the zero pressure at the inlet of the curvy channel it could be seen that for Re = 1000 atP= 5,the pressure declines about 600%with an average declination of 500%.Using a curve-fitting procedure we can estimate some sinusoidal regression curves forP= 5 which are given below.The curves give a better estimate of where major data points are lying for the pressure related to x-coordinate.

    Figure 14: Pressure distribution along the length of the curvy channel with the period P = 3 (a)Re=100,(b)Re=500 and(c)Re=1000

    Figure 15: Pressure distribution along the length of the curvy channel with the period P = 5 (a)Re=100,(b)Re=500 and(c)Re=1000

    4.3 The Wall Shear Stress

    In the field of fluid dynamics,the wall shear stress can be regarded as the rate at which the velocity inclines while the fluid is moving along the wall.In common, shear stresses are occurring when the elements of fluid acquiring the motion relative to one another.For a simple pipe, zero velocity is assumed at the pipe wall whereas it is increasing while the fluid element is moving towards the center of the pipe.For a non-smooth pipe,the determination of the wall shear stress is essential.Here,we are going to present the graphs for wall shear stress against the length of the pipe see Figs.16a–16c to 18a–18c.The pattern of the graphs can be matched to the pattern of the graphs reported in[24]where the channel was considered for the two period of vibration and here the problem is being focused on more period of vibration of the channel.From the graphs,it is evident that with the increase in Reynolds number the wall shear stress at the inlet of the channel increases whereas the wall shear stress at the outlet of the channel is zero for all the cases due to the zero pressure condition.The pattern of behavior of the wall shear stress along the length of the channel is not very much clear for the moderate number of periodsP= 2 andP= 3.Even the Reynolds number is changing from lower to higher values to observe the pattern see Figs.16a–16c and 17a–17c.Whereas the wall shear stress at the contracted region was found maximum in these cases.With the increase in the period of vibration of the walls as shown in Figs.18a–18c,the wall shear stress is decreasing after each contraction region.For a low Reynolds number Re=100,the wall shear stress shows declination like a step function of the length.For a particular Reynolds number,the wall shear stress is increasing with the increase in periodsPof the vibration i.e.,compare Fig.16a with Fig.17a and 18a.

    Figure 16:Wall shear stress against the length of channel for P = 2(a)Re = 100,(b)Re = 500 and(c)Re=1000

    Figure 17:Wall shear stress against the length of channel for P = 3(a)Re = 100,(b)Re = 500 and(c)Re=1000

    Figure 18:Wall shear stress against the length of channel for P = 5 and(a)Re = 100,(b)Re = 500 and(c)Re=1000

    4.4 Measurement of Rotating Regions

    It is often can be seen in the fluid flow problems which are observed under the presence of the bluff bodies,cavities and might be they can be observed due to the periodic motion of the boundaries of the channel,the fluid often rotates along its axis.The rotational regions becomes clearer when the fluid is observed under the presence of bluff bodies or obstacles rather than without of them.When fluid with any material slides down upstream of the channel,the rotating regions can be observed along the sides of the downstream.The length and intensity of these rotating regions depend upon the velocity of the fluids.The fluid needs some velocity magnitude to form the vortex along with the downstream of the region.While observing the wavy channel in the current problem the rotating regions were observed in Figs.5d–5f.In Figs.19a–19f, we are presenting graphs where we are measuring the rotation rate downstream of these rotating regions against the local Reynolds number.

    Figure 19:Relationship of rotation rate(s-1)with local Reynolds number for Re=1000 and(a)P=0(b)P=1(c)P=2(d)P=3(e)P=4 and(f)P=5

    Although, there was no rotational region found for the periodP= 0 f the vibration even we can measure the rotation rate at the end of the wavy channel.From Fig.1a,we found the maximum rotation per second for a zero period is 30 rotations per second.Whereas the maximum rotation for periodP=1 andP=2 of the vibration,we found 50 and 220 rotation per second respectively.These are the cases in which the rotating regions can be observed in the presentation of streamlines.Now look at the cases forP= 3 toP= 5,in these cases,the rotating regions are obvious for Re = 1000 see Figs.5a–5f.On observing these rotation regions forP= 3,4,5 we found the maximum rotation rates per second are 500, 900, and 1500 respectively.One more thing is that the rotation rate of the fluid is very place to place for example we can see Figs.19d–19f the maximum rotation per second is optimum at expansion regions.

    4.5 Temperature

    The heat distribution in the channel is investigated by imposing a temperature at the outlet of the channel and the air is allowed to enter the domain with an initial temperature.Our concern is to lessen the temperature near to the outlet of the channel by searching such a Reynolds number and the period of the vibration of the boundaries.The temperature distributions with the surface plots in the channel are presented in Figs.20a–20l for the particular Reynolds number as well as the Periods of the vibration of the channel.On following Fig.10 for Re = 100 the average temperature near the outlet is increasing with the increase in the period of the vibration of the boundaries.Similar cases can be perceived for the other Reynolds numbers of 500 and 1000 in Fig.20a–20l.In the contrast, with fixing the periods of the vibration,see Figs.20a–20c,an increase in the Reynolds number will decrease in the average temperature near the outlet of the channel.A similar pattern of the distribution near the outlet of the channel can be recognized for example Figs.20d–20l.

    Figure 20: Temperature distribution in the channel from cold place to hot place in the periods (a-c)P=0,(d-f)P=2,(g-i)P=4 and(j-l)P=5

    Therefore,the graphical presentation(see Fig.21a)of the average temperature near the outlet of the channel against the Reynolds number is depicted for the particular periods.The graph is showing the same consequences as described above in Figs.20a–20l.It is also apparent from the graph that the declination of the average temperature persists consistently after Re = 300 i.e., there will be no declination of the temperature near the outlet of the channel.The average Nusselt number at the outlet of the channel is estimated for all Reynolds numbers at the particular period of the vibration, see Fig.21b.It appears from the graph that the convection process near the outlet is growing stronger with the increase in Reynolds number as well as the period of the vibration of the boundaries.

    Figure 21: (a) Average temperature against the Reynolds number near the outlet of the channel (b)Average Nusselt number against Reynolds number near the outlet of the channel

    5 Conclusion

    The present work encountered airflow and temperature distribution through a sinusoidally curved wavy channel of length 2π.Carrying out parametric study by varying Reynolds number in the range from 100 to 1000 and period of vibration from 0 to 5.COMSOL Multiphysic 5.4 was used to solve the nonlinear equations and to simulate the results.A mesh independent study of the finite element method was performed with the use of the different configurations of meshing by computing u/uinat the downstream of the channel.To approve the results,the local Nusselt number at the lower boundary of the channel was compared with the correlation available in the literature.The velocity magnitude,pressure,temperature,maximum velocity,minimum pressure,average temperature at the outlet,and the local Nusselt number at the outlet was estimated in terms of Reynolds number at particular period of vibration of the boundary.This study is important because of its applications in various engineering fields.One such problem is in the study of heat exchangers.A sinusoidally curved wavy channel which satisfies the need of rapid heat transfer is always needed.The results were presented graphically and following conclusions are drawn:

    1.Maximum velocity magnitude is found at the contraction regions.At the expansion regions,small recirculating vortices are created which resist the speed of the flow while the fluid is sliding down the region.The maximum velocity magnitude at the outlet is increasing linearly with the increase in Reynolds number but it is very low forP=0.

    2.The pressure in the channel is continuously decreasing after passing through contraction regions for all the Reynolds numbers and the period of vibration of the boundaries.It is suggested that using the wavy channel with an equal number of periods can be used as the pressure reducing pump in most engineering applications.

    3.The pressure distribution along the length of the channel is having cyclic variation along the channel.Due to cyclic fluctuation in the pressure,it is obvious to find the maximum pressure at the expansion region and minimum pressure at the contraction region.Due to considering the periodic variation,we can obtain a regression equation involving a sinusoidal function that can cover the major points along the length of the channel.In the article,we have recommended some equations for lower,moderate,and high Reynolds numbers.

    4.It was also found that for low and moderate Reynolds numbers,the shear stress at the lower wall of the channel shows a periodic variation from inlet to outlet of the channel.Whereas,for high Reynolds number, the wall shear stress at the bottom wall is decreasing after each contraction region.This means the velocity gradient after each contraction region declines than the previous.

    5.By measuring the rotation rate in terms of period of vibration of the channel it was seen that forP= 0 (straight channel) only some rotation rate was measured near the outlet of the channel.Whereas,for the other periods ofP= 1 toP= 3,maximum rotation rates were deducted at the expansion regions of the channel.Moreover,the rotation rate is increasing with the increase in Reynolds number as well as the period of vibration.

    6.A flow of heat from a cold body to a hot body was tested for the wavy channel.It was concluded that the average temperature near the hot body is decreasing for Re<300 and remains constant for 300<Re<1000.Moreover, an increase in the period of vibration of boundaries will increase the average temperature near the outlet of the channel.

    Acknowledgment:The authors are highly indebted to anonymous referees for their helpful suggestions and comments which led to improvements of our original manuscript.

    Funding Statement:This research was funded by King Mongkut’s University of Technology North Bangkok.Contract no.KMUTNB-63-KNOW-20.

    Conflicts of Interest:The authors declare that they have no conflicts of interest to report regarding the present study.

    久久久久久免费高清国产稀缺| 88av欧美| 男女做爰动态图高潮gif福利片| 欧美黑人精品巨大| 精品国产乱码久久久久久男人| 99riav亚洲国产免费| 久久久国产成人精品二区| 可以在线观看毛片的网站| 黑人巨大精品欧美一区二区mp4| 大型黄色视频在线免费观看| 午夜福利免费观看在线| 日本一区二区免费在线视频| 男男h啪啪无遮挡| 午夜两性在线视频| 深夜精品福利| 亚洲五月婷婷丁香| 亚洲欧美日韩高清在线视频| 不卡一级毛片| 88av欧美| 久9热在线精品视频| 美女扒开内裤让男人捅视频| 一级a爱视频在线免费观看| 亚洲美女黄片视频| 亚洲九九香蕉| 精品日产1卡2卡| 一本综合久久免费| 亚洲国产精品999在线| 视频区欧美日本亚洲| 欧美色视频一区免费| 日韩欧美国产在线观看| 亚洲专区国产一区二区| 亚洲,欧美精品.| 一个人免费在线观看的高清视频| 好男人在线观看高清免费视频 | 每晚都被弄得嗷嗷叫到高潮| 国产一区二区在线av高清观看| 嫁个100分男人电影在线观看| 色av中文字幕| 日本成人三级电影网站| 亚洲第一青青草原| 老鸭窝网址在线观看| 热99re8久久精品国产| 成人一区二区视频在线观看| 神马国产精品三级电影在线观看 | 91国产中文字幕| 一卡2卡三卡四卡精品乱码亚洲| 香蕉久久夜色| 亚洲国产毛片av蜜桃av| 51午夜福利影视在线观看| 三级毛片av免费| 真人做人爱边吃奶动态| 国产av在哪里看| 日韩视频一区二区在线观看| 成年免费大片在线观看| 曰老女人黄片| 草草在线视频免费看| 欧美zozozo另类| 精品欧美一区二区三区在线| 曰老女人黄片| 国产主播在线观看一区二区| 高清毛片免费观看视频网站| 一级毛片精品| 亚洲国产毛片av蜜桃av| 男女视频在线观看网站免费 | 国产亚洲精品一区二区www| 女人高潮潮喷娇喘18禁视频| 黄频高清免费视频| 美女扒开内裤让男人捅视频| 在线观看www视频免费| 国产av一区二区精品久久| av视频在线观看入口| 亚洲欧美精品综合久久99| 无限看片的www在线观看| 97人妻精品一区二区三区麻豆 | 中文字幕av电影在线播放| 亚洲精品中文字幕在线视频| 老汉色av国产亚洲站长工具| 可以在线观看的亚洲视频| 真人做人爱边吃奶动态| 18禁黄网站禁片免费观看直播| 日韩视频一区二区在线观看| 国产成人欧美在线观看| 一二三四在线观看免费中文在| 亚洲精品久久成人aⅴ小说| 满18在线观看网站| 美女扒开内裤让男人捅视频| 国产亚洲精品久久久久久毛片| 草草在线视频免费看| 母亲3免费完整高清在线观看| 女人被狂操c到高潮| 黑人操中国人逼视频| 神马国产精品三级电影在线观看 | 人成视频在线观看免费观看| 国产单亲对白刺激| 啦啦啦韩国在线观看视频| 亚洲精品一卡2卡三卡4卡5卡| 亚洲成人免费电影在线观看| 黄色视频不卡| 亚洲成国产人片在线观看| 中文字幕精品免费在线观看视频| 欧美在线黄色| 久久天躁狠狠躁夜夜2o2o| 伦理电影免费视频| 韩国精品一区二区三区| 欧美性猛交黑人性爽| 国产麻豆成人av免费视频| 一区二区三区国产精品乱码| 男女下面进入的视频免费午夜 | 9191精品国产免费久久| 国产亚洲精品第一综合不卡| 丝袜美腿诱惑在线| 亚洲欧美精品综合久久99| 777久久人妻少妇嫩草av网站| 男人舔女人下体高潮全视频| 听说在线观看完整版免费高清| 久久精品91无色码中文字幕| 99在线视频只有这里精品首页| tocl精华| xxxwww97欧美| 国产又黄又爽又无遮挡在线| av电影中文网址| 男人舔奶头视频| 少妇粗大呻吟视频| 精品国产美女av久久久久小说| 亚洲精品久久国产高清桃花| 欧美色欧美亚洲另类二区| 一级毛片女人18水好多| 欧美日韩黄片免| 欧美激情极品国产一区二区三区| 两性夫妻黄色片| 美女扒开内裤让男人捅视频| 国产亚洲精品av在线| 成人国产综合亚洲| 亚洲全国av大片| 午夜福利一区二区在线看| 久久精品影院6| 在线观看免费日韩欧美大片| 真人做人爱边吃奶动态| 久久精品aⅴ一区二区三区四区| a级毛片在线看网站| 天堂影院成人在线观看| 欧美又色又爽又黄视频| 欧美zozozo另类| 欧美大码av| 不卡一级毛片| 亚洲欧美一区二区三区黑人| 最新美女视频免费是黄的| 久久久久久免费高清国产稀缺| 欧美黑人欧美精品刺激| 久久久精品国产亚洲av高清涩受| 一级毛片精品| 国产高清激情床上av| 亚洲成av人片免费观看| 色综合婷婷激情| 国产精品久久久av美女十八| 91国产中文字幕| 丁香六月欧美| 日韩精品免费视频一区二区三区| 制服诱惑二区| 香蕉久久夜色| 国产精品一区二区三区四区久久 | 欧美激情久久久久久爽电影| 在线天堂中文资源库| 最近在线观看免费完整版| 国产亚洲精品av在线| 国产高清激情床上av| 后天国语完整版免费观看| 美国免费a级毛片| 一级a爱片免费观看的视频| 一个人观看的视频www高清免费观看 | 午夜视频精品福利| 欧美成人性av电影在线观看| 国产精品电影一区二区三区| 亚洲成人精品中文字幕电影| 欧美色视频一区免费| 国产精品亚洲美女久久久| 一进一出好大好爽视频| av片东京热男人的天堂| 欧美精品亚洲一区二区| 长腿黑丝高跟| 99热这里只有精品一区 | 在线看三级毛片| 欧美激情久久久久久爽电影| 亚洲男人天堂网一区| 一级作爱视频免费观看| 国产成+人综合+亚洲专区| 久久99热这里只有精品18| 亚洲人成电影免费在线| av免费在线观看网站| 啦啦啦 在线观看视频| 国产单亲对白刺激| 国内少妇人妻偷人精品xxx网站 | 夜夜看夜夜爽夜夜摸| 亚洲,欧美精品.| 亚洲七黄色美女视频| 成年女人毛片免费观看观看9| 午夜免费鲁丝| 日韩欧美 国产精品| 天天添夜夜摸| 成年免费大片在线观看| 国产免费男女视频| 精品久久久久久,| 成人免费观看视频高清| 在线观看免费午夜福利视频| 亚洲av中文字字幕乱码综合 | 欧美日韩乱码在线| 女人被狂操c到高潮| 亚洲自偷自拍图片 自拍| 成人国产一区最新在线观看| 国产精品九九99| 久99久视频精品免费| 老司机深夜福利视频在线观看| 国产不卡一卡二| 亚洲五月婷婷丁香| 亚洲国产精品合色在线| 国产精品自产拍在线观看55亚洲| 91国产中文字幕| 欧美中文综合在线视频| 久久精品国产综合久久久| 精品国产一区二区三区四区第35| 一本久久中文字幕| 精品久久久久久久人妻蜜臀av| 国产一区二区三区视频了| 久久精品国产清高在天天线| 亚洲国产精品成人综合色| 午夜免费鲁丝| 香蕉av资源在线| 最新在线观看一区二区三区| 成人av一区二区三区在线看| 一边摸一边抽搐一进一小说| 国产亚洲精品av在线| 欧美激情 高清一区二区三区| 国产极品粉嫩免费观看在线| 国产精品一区二区精品视频观看| 国产一区二区三区在线臀色熟女| 在线观看免费午夜福利视频| 久久精品国产综合久久久| 国内少妇人妻偷人精品xxx网站 | 天天添夜夜摸| 免费在线观看日本一区| 亚洲成人国产一区在线观看| 亚洲国产中文字幕在线视频| 色在线成人网| 在线观看免费午夜福利视频| 欧美日韩黄片免| 18禁美女被吸乳视频| 国内久久婷婷六月综合欲色啪| 757午夜福利合集在线观看| 国产高清激情床上av| 非洲黑人性xxxx精品又粗又长| 精品高清国产在线一区| 国产成人系列免费观看| 大型av网站在线播放| 国产亚洲欧美精品永久| 欧美精品亚洲一区二区| 亚洲av第一区精品v没综合| 中文字幕久久专区| 在线免费观看的www视频| 成人国语在线视频| 久久国产精品人妻蜜桃| 少妇粗大呻吟视频| 国产亚洲av嫩草精品影院| 免费在线观看视频国产中文字幕亚洲| 亚洲熟妇中文字幕五十中出| 国产aⅴ精品一区二区三区波| 午夜免费成人在线视频| 一级片免费观看大全| 啦啦啦韩国在线观看视频| 欧美黑人精品巨大| 色av中文字幕| 97碰自拍视频| 制服丝袜大香蕉在线| 国产av一区在线观看免费| 成人特级黄色片久久久久久久| 国产精品免费一区二区三区在线| 国产成人精品无人区| 国产精品久久电影中文字幕| 欧美精品亚洲一区二区| 日韩三级视频一区二区三区| 首页视频小说图片口味搜索| 777久久人妻少妇嫩草av网站| 亚洲,欧美精品.| 亚洲一卡2卡3卡4卡5卡精品中文| svipshipincom国产片| 动漫黄色视频在线观看| 国产三级黄色录像| 国产真实乱freesex| 丁香六月欧美| 波多野结衣高清无吗| 精品国产乱码久久久久久男人| 免费在线观看黄色视频的| 久久国产亚洲av麻豆专区| 日韩成人在线观看一区二区三区| a在线观看视频网站| 国产精品乱码一区二三区的特点| 午夜福利免费观看在线| 日韩大尺度精品在线看网址| 婷婷精品国产亚洲av在线| 琪琪午夜伦伦电影理论片6080| 日韩欧美国产一区二区入口| 国产精品免费一区二区三区在线| 久久午夜亚洲精品久久| 国产亚洲精品综合一区在线观看 | 中文字幕精品免费在线观看视频| 丁香欧美五月| av超薄肉色丝袜交足视频| 男女那种视频在线观看| 久久久久亚洲av毛片大全| 国产精品二区激情视频| 成人亚洲精品av一区二区| 欧美激情久久久久久爽电影| 久久精品影院6| av片东京热男人的天堂| 黑人操中国人逼视频| 亚洲国产精品999在线| 黄频高清免费视频| 精品无人区乱码1区二区| 国产午夜福利久久久久久| 中国美女看黄片| 亚洲午夜精品一区,二区,三区| 波多野结衣高清作品| 中文在线观看免费www的网站 | 免费高清在线观看日韩| 欧美最黄视频在线播放免费| 香蕉av资源在线| www.精华液| 国产av在哪里看| 欧美成狂野欧美在线观看| av免费在线观看网站| 大型av网站在线播放| 国产成+人综合+亚洲专区| 国产精品日韩av在线免费观看| 国产久久久一区二区三区| 在线视频色国产色| 亚洲欧美激情综合另类| 神马国产精品三级电影在线观看 | 看免费av毛片| 91麻豆精品激情在线观看国产| 99热6这里只有精品| 曰老女人黄片| www日本黄色视频网| 18禁美女被吸乳视频| 国产av又大| 亚洲精品av麻豆狂野| 动漫黄色视频在线观看| 精品久久久久久,| 91成年电影在线观看| 午夜福利视频1000在线观看| 国产伦一二天堂av在线观看| 黄色 视频免费看| 免费电影在线观看免费观看| 欧美人与性动交α欧美精品济南到| 亚洲第一电影网av| 亚洲国产高清在线一区二区三 | 亚洲成av人片免费观看| 亚洲一区二区三区不卡视频| 十八禁网站免费在线| 精品久久久久久成人av| 午夜福利在线在线| 两性夫妻黄色片| 国产亚洲精品av在线| 日韩 欧美 亚洲 中文字幕| 麻豆国产av国片精品| 国产亚洲精品一区二区www| 国产精品久久久久久亚洲av鲁大| 十八禁网站免费在线| 亚洲男人的天堂狠狠| 午夜影院日韩av| 中亚洲国语对白在线视频| 成人永久免费在线观看视频| 国产私拍福利视频在线观看| 1024手机看黄色片| 国产免费男女视频| 国产精品九九99| 黄片播放在线免费| 中文字幕人成人乱码亚洲影| 高清毛片免费观看视频网站| 天天一区二区日本电影三级| 国产精品精品国产色婷婷| 亚洲欧美日韩高清在线视频| 黄色女人牲交| e午夜精品久久久久久久| 国产亚洲精品久久久久久毛片| 1024香蕉在线观看| 国产97色在线日韩免费| 亚洲成人免费电影在线观看| 久久久久久久精品吃奶| 精品国产国语对白av| 中文字幕人妻丝袜一区二区| 精品日产1卡2卡| 久久欧美精品欧美久久欧美| 国产精品久久久久久人妻精品电影| 老司机福利观看| 老司机深夜福利视频在线观看| 一a级毛片在线观看| 日本撒尿小便嘘嘘汇集6| 变态另类成人亚洲欧美熟女| 一二三四社区在线视频社区8| 免费在线观看黄色视频的| 激情在线观看视频在线高清| a级毛片在线看网站| www日本在线高清视频| 少妇裸体淫交视频免费看高清 | 国产成人啪精品午夜网站| 午夜a级毛片| 国产aⅴ精品一区二区三区波| 国产高清激情床上av| 精品久久久久久成人av| 黄片小视频在线播放| 黄色丝袜av网址大全| 免费在线观看完整版高清| 法律面前人人平等表现在哪些方面| 精品久久久久久,| www国产在线视频色| www日本黄色视频网| 国产激情久久老熟女| 欧美在线黄色| 精品久久久久久成人av| 91麻豆精品激情在线观看国产| √禁漫天堂资源中文www| 成人三级黄色视频| 国产精品 欧美亚洲| 99re在线观看精品视频| 丰满人妻熟妇乱又伦精品不卡| 一a级毛片在线观看| 国产片内射在线| 草草在线视频免费看| 观看免费一级毛片| АⅤ资源中文在线天堂| 国产亚洲精品综合一区在线观看 | 免费在线观看黄色视频的| 久久久久久人人人人人| 黄色片一级片一级黄色片| 女人高潮潮喷娇喘18禁视频| 亚洲一区高清亚洲精品| 久热这里只有精品99| 成人18禁在线播放| 亚洲精品美女久久久久99蜜臀| 久久伊人香网站| 免费电影在线观看免费观看| 成人欧美大片| 俄罗斯特黄特色一大片| 亚洲色图 男人天堂 中文字幕| 欧美日韩精品网址| 国产精品久久久久久人妻精品电影| 欧美精品啪啪一区二区三区| 久久久久国内视频| 亚洲成人久久性| 人人妻人人澡人人看| 国产精品影院久久| 精品久久蜜臀av无| av欧美777| 观看免费一级毛片| 久久久国产精品麻豆| 亚洲av片天天在线观看| 亚洲七黄色美女视频| 亚洲五月色婷婷综合| 日韩大码丰满熟妇| 色哟哟哟哟哟哟| 亚洲欧美激情综合另类| 国产精品香港三级国产av潘金莲| 亚洲国产中文字幕在线视频| 深夜精品福利| 国产精品电影一区二区三区| 久久中文字幕人妻熟女| 久久国产精品影院| 精品欧美国产一区二区三| 97碰自拍视频| 欧美另类亚洲清纯唯美| 美女高潮喷水抽搐中文字幕| 热re99久久国产66热| 欧美激情高清一区二区三区| 婷婷精品国产亚洲av| 午夜精品在线福利| 亚洲真实伦在线观看| 精品久久久久久成人av| 久久精品影院6| 中文资源天堂在线| 中国美女看黄片| 12—13女人毛片做爰片一| 亚洲欧美激情综合另类| 在线观看一区二区三区| 一级毛片女人18水好多| 日韩欧美免费精品| 久久国产精品影院| 无人区码免费观看不卡| 久久狼人影院| 正在播放国产对白刺激| 国产一区二区激情短视频| 国产伦在线观看视频一区| 国产精品 国内视频| 夜夜看夜夜爽夜夜摸| 国产麻豆成人av免费视频| 国产一区在线观看成人免费| 国产av一区在线观看免费| 国产三级黄色录像| 女性被躁到高潮视频| 成人手机av| 老熟妇仑乱视频hdxx| 国产单亲对白刺激| 黄片小视频在线播放| 国产1区2区3区精品| or卡值多少钱| 欧美日韩精品网址| 欧美乱妇无乱码| 亚洲精品在线观看二区| 精品久久久久久久久久久久久 | 欧美另类亚洲清纯唯美| 亚洲av成人av| 在线十欧美十亚洲十日本专区| 99久久无色码亚洲精品果冻| 午夜福利高清视频| 91麻豆av在线| 亚洲久久久国产精品| 一本大道久久a久久精品| 两个人视频免费观看高清| 国产一级毛片七仙女欲春2 | 中文字幕人成人乱码亚洲影| 人妻丰满熟妇av一区二区三区| 91成年电影在线观看| 亚洲av五月六月丁香网| 可以在线观看的亚洲视频| 中文字幕人成人乱码亚洲影| 看免费av毛片| 久久中文字幕人妻熟女| 黑人巨大精品欧美一区二区mp4| 男人操女人黄网站| 波多野结衣高清无吗| 国产激情久久老熟女| 十八禁网站免费在线| 好看av亚洲va欧美ⅴa在| 高清毛片免费观看视频网站| 日韩免费av在线播放| 老司机福利观看| 久久久水蜜桃国产精品网| 中文字幕精品免费在线观看视频| 中国美女看黄片| 日韩精品青青久久久久久| 少妇的丰满在线观看| 在线永久观看黄色视频| 国产精品野战在线观看| 一个人观看的视频www高清免费观看 | 日本a在线网址| 色哟哟哟哟哟哟| 久久久久国产一级毛片高清牌| 国产一卡二卡三卡精品| 精品国产国语对白av| 最近最新中文字幕大全免费视频| 亚洲精品国产区一区二| 国产97色在线日韩免费| 在线永久观看黄色视频| 亚洲国产欧美网| 性色av乱码一区二区三区2| 久久精品91蜜桃| 国产亚洲欧美98| 欧美黄色片欧美黄色片| 精品国产国语对白av| 亚洲五月天丁香| 成人国产一区最新在线观看| 国产片内射在线| 美女高潮到喷水免费观看| 日韩精品中文字幕看吧| 女警被强在线播放| 日韩欧美一区二区三区在线观看| 久久国产亚洲av麻豆专区| 日本 欧美在线| 国产精品久久久久久亚洲av鲁大| 国产激情久久老熟女| 国产午夜精品久久久久久| 老司机在亚洲福利影院| 午夜激情福利司机影院| 久久精品国产清高在天天线| 熟妇人妻久久中文字幕3abv| 好男人在线观看高清免费视频 | 午夜成年电影在线免费观看| 成人亚洲精品一区在线观看| 精品福利观看| 免费在线观看成人毛片| 亚洲第一青青草原| 美国免费a级毛片| 亚洲国产欧美一区二区综合| 日韩欧美国产一区二区入口| 制服人妻中文乱码| 丝袜美腿诱惑在线| 午夜激情av网站| 99久久国产精品久久久| 亚洲精品久久成人aⅴ小说| 窝窝影院91人妻| 欧美黑人巨大hd| 婷婷精品国产亚洲av在线| 欧美久久黑人一区二区| 亚洲欧美一区二区三区黑人| 丝袜在线中文字幕| 精品欧美一区二区三区在线| 日韩三级视频一区二区三区| 国产精品1区2区在线观看.| 成年免费大片在线观看| 亚洲熟妇中文字幕五十中出| 欧美不卡视频在线免费观看 | 国产精品永久免费网站| 叶爱在线成人免费视频播放| 黑人巨大精品欧美一区二区mp4| 国产成人一区二区三区免费视频网站| 草草在线视频免费看| 亚洲一区二区三区不卡视频| 亚洲色图 男人天堂 中文字幕| 18禁黄网站禁片免费观看直播| 人人澡人人妻人| 日本a在线网址| 自线自在国产av| 久久伊人香网站| 午夜亚洲福利在线播放| 99久久久亚洲精品蜜臀av| 亚洲第一电影网av| 亚洲性夜色夜夜综合|