亚洲免费av电影一区二区三区,日韩爱爱视频,51精品视频一区二区三区,91视频爱爱,日韩欧美在线播放视频,中文字幕少妇AV,亚洲电影中文字幕,久久久久亚洲av成人网址,久久综合视频网站,国产在线不卡免费播放

        ?

        Finite Difference Analysis of Time-Dependent Viscous Nanofluid Flow Between Parallel Plates

        2019-11-07 02:58:24SalmanAhmadHayatAlsaediKhanandWaleedAhmedKhan
        Communications in Theoretical Physics 2019年11期

        Salman Ahmad,T.Hayat,, A.Alsaedi, Z.H.Khan, and M.Waleed Ahmed Khan

        1Department of Mathematics, Quaid-I-Azam University 45320, Islamabad 44000, Pakistan

        2Nonlinear Analysis and Applied Mathematics(NAAM)Research Group,Department of Mathematics,Faculty of Science,King Abdulaziz University, P.O.Box 80257, Jeddah 21589, Saudi Arabia

        3State Key Laboratory of Hydraulics and Mountain River Engineering, College of Water Resource and Hydropower,Sichuan University, Chengdu 610065, China

        4Key Laboratory of Advanced Reactor Engineering and Safety, Ministry of Education, Tsinghua University, Beijing 100084, China

        Abstract This article computes effect of uniform magnetic field on nanofluid flow filling porous medium inside parallel sheets.Darcy’s law is used to characterize porous medium.Flow is caused by stretching of the lower sheet.The governing systems are computed using Finite difference technique.This scheme is used due to its more general and powerful nature to solve nonlinear problems.Results are obtained and discussed graphically.Nusselt number, skin friction, streamlines,velocity and temperature are emphasized.

        Key words:finite difference scheme, porous medium, magnetic field, nanofluid, parallel stretchable sheets

        Nomenclature

        1 Introduction

        Heat transfer rate has played a key role in saving energy and experimental costs.Presently technological fields need devices with more precision, accurate functioning,higher performance and more lifespan.That is why researchers are working for the thermal properties of materials.Nanofluids are the mixture of nanosized metallic as well as nonmetallic particles in some base fluids used for increasing thermal conductivity of fluids.In 1995 Choi[1]found the enhanced thermal properties of nanofluids as compared to simple fluids which compelled other researchers to give attention to such nanosized particles with extraordinary characteristics.Moreover, Buongiorno[2]discussed seven slip phenomena of nanofluid (i.e.inertia,Brownian diffusion, fluid drainage, diffusiophoresis, gravity, thermophoresis and Magnus).However he concluded that Brownian motion and thermophoresis are the key slip phenomena.Different types of particles e.g.copper (Cu),gold (Au), silver (Ag), and nickel (Ni), metal oxides like copper oxide(CuO),silicon dioxide(SiO2),oxide(Fe2O3),iron(III),titanium dioxide(TiO2)etc.[3?8]are extensively used by different researchers.Flow of nanofluid due to a rotating disk is examined by Turkyilmazoglu.[9]For cooling process, spraying of nanofluid by a rotating inclined disk is numerically investigated by Sheikholeslamiet al.[10]

        Magnetohydrodynamic flows over a stretching sheet have diverse applications in industrial processes.Examples of such processes include polymer sheet industries,cooling processes in many industries, glass sheets, parmachaology, biotechnology, nuclear research, glass sheet industries, drawing of plastic sheets, polymers extrusion,cooling processes of melting materials etc.In all such cases the flow situation is of vital importance.Quality of end product mainly depends upon the rate of heat transfer and skin friction.Such processes also depend on the boundary layer phenomenon along stretching surfaces as well as the mass transfer rate.Rajogopalet al.[11]initially studied viscoelastic fluid flow phenomenon caused by a stretching surface.MHD flow by a vertical plate was examined by Riley.[12]Fairbanks and Wike[13]analyzed the effect of chemical reaction in a uniform fluid flow over a horizontal sheet.Anderssonet al.[14]investigated chemically reactive flow over a stretching surface.The impacts of magnetic dipole and thermal radiation on Willamsion fluid flow are discussed by Hayatet al.[15]Some more applications are highlighted by Magyari and Keller.[16?17]

        Formulation of flow problems in the form of partial or ordinary differential equations is very important for precise modeling and simulation in engineering processes.These classes of differential equations are tackled both analytically as well as numerically.In this article numerical approach has been utilized to tackle the partial differential equations governing the problem of MHD flow of nanofluid in porous medium between two stretchable sheets.Applications of finite difference scheme[18?20]have been studied for solving the system of partial differential equations.Graphical results are obtained through various involved parameters for velocity, concentration and temperature.Quantities of engineering interests i.e.skin friction, mass and heat transfer rates are also discussed through graph.Main points of the present article are summarized in the concluding section.

        2 Mathematical Modeling

        In this section we model the flow of time-dependent viscous nanofluid inside two parallel sheets.Nanofluid occupies porous medium and Darcy’s law is used to characterized porous space.Influence of magnetic field is perpendicular to x-axis.Injection is taken on porous lower plate with stretching ratea.Upper plate is at a distanceLfrom lower plate.Temperature of lower plate isTwand upper plateTh.Figure 1 shows the flow diagram.Under the considered assumptions the flow expressions are:

        with constraints

        Fig.1 (Color online) Flow diagram.

        In order to transform the system into dimensionless form, consider the dimensionless variables

        Using dimensionless variables the flow expressions are reduced to

        3 Quantities of Engineering Interest

        Sherwood number, Nusselt number and skin friction are defined as

        whereqmrepresents mass flux,qwheat flux andτwshear stress.These are

        In dimensionless forms the Sherwood number, Nusselt number and skin friction are given by

        4 Constructing Finite Difference Scheme

        In order to solve the partial differential system given in Eqs.(8)–(12)with constraints(13),we used Finite difference scheme.The system contains four unknown function(U,V,θ,?)and three independent variables(τ,ξ,η).On the basis of finite difference technique, FD toolkit for partial derivatives are given by

        Using Eq.(17), partial differential equations system (Eqs.(7)–(10)) are reduced to

        with

        Sherwood and Nusselt numbers and skin friction become

        5 Results and Discussion

        Here we discuss the outcomes of velocity components(U(τ,ξ,η),V(τ,ξ,η)), concentration (?(τ,ξ,η)), temperature (θ(τ,ξ,η)), surface drag force (Cfx), heat transfer rate (Nux) and Sherwood number (Shx) with variation in pertinent parameters.Pertinent parameters included Reynolds number(M),Prandtl number(Pr),porosity parameter(λ),Hartmann number(Ha),injection parameter(β), thermophoresis variable (Nt) and Brownian variable(Nb).

        Fig.2 (Color online) Behaviors of β and ξ on U.

        Impacts ofβ,λ,MandHaon horizontal and vertical components of velocity are plotted graphically via Figs.2–7.

        Fig.3 (Color online) Behaviors of M and τ on U.

        Figure 2 shows the impact of injection variable on horizontal component of velocityU(τ,ξ,η).It shows that injection process tends to velocity enhancement due to an increase in boundary layer thickness.Figure 3 is sketched to show the behavior of velocityU(τ,ξ,η) with respect to time as well as Reynolds number (M).Increasing values of Reynolds number produces more turbulence in the flow.This is the reason for more disturbance inside the fluid particles which increases the velocity in x-direction.This figure also shows that velocity increases with time.Impacts of porosity parameter(λ)and Hartmann number(Ha) on velocity are portrayed in Fig.4.Since porosity parameter is inversely proportional to permeability of medium, so more porous medium decreases the speed of fluid particles.Hence velocity decreases.Similarly Hartmann number also decays the velocity due to Lorentz force produced inside the system, which slows the fluid particles.Injection parameter increases the velocity componentV(τ,ξ,η)due to injection of more fluid particles inside the system as shown in Fig.5.Vertical component of velocityV(τ,ξ,η) with respect to material parameterMand time is given in Fig.6.Stretching is produced in the lower plate, so less disturbance is felt at the upper plate.That is why velocity has a decreasing behavior due to material parameter, while it increases with respect to time.Behaviors of vertical velocity componentV(τ,ξ,η) due to Hartmann numberHaand porosity parameter are shown in Fig.7.Since magnetic field is parallel to velocity componentV(τ,ξ,η) so it has no effect on vertical velocityV(τ,ξ,η).However it decreases with respect to porosity parameter because of the decrease in the permeability of medium.

        Fig.4 (Color online) Behaviors of λ and Ha on U.

        Fig.5 (Color online) Behaviors of β and ξ on V.

        Effects of thermophoresis variable (Nt), Prandtl number (Pr) and Brownian variable (Nb) on temperatureθ(τ,ξ,η)are described in Figs.8–9.Figure 8 is delineated to show the impact ofPrand time on temperature.Since the lower plate is at higher temperature so temperature of the system gradually decreases with time.However it increases withPrdue to decreases in thermal diffusivity of fluid.As a result the system is heated up which consequently increases the temperature.Figure 9 shows the influences ofNtandNbon temperature.Here boundary layer thickness increases for largerNtandNband consequently temperature enhances.

        Fig.6 (Color online) Behaviors of τ and M on V.

        Fig.7 (Color online) Effects of λ and Ha on V.

        Fig.8 (Color online) Effects of Pr and τ on θ.

        Behaviors of thermophoresis variables (Nt), Schmidt number (Sc) and Brownian parameter (Nb) on concentration (?(τ,ξ,η)) are disclosed in Figs.10–11.Effects ofScand time on?(τ,ξ,η) are displayed in Fig.10.AsScis the relation of viscous diffusion rate and mass diffusion rate so largerScenhances the concentration field.With passing time the concentration show decreasing behavior.Likewise, Fig.11 is displayed to present the effects ofNtandNbon?(τ,ξ,η).Here concentration field increases forNbbut it decays for thermophoresis parameterNt.

        Fig.9 (Color online) Behaviors of Nt and Nb on θ.

        Fig.10 (Color online) Behaviors of Sc and τ on ?.

        Fig.11 (Color online) Behaviors of Nt and Nb on ?.

        Skin friction, mass and heat transfer rates with respect to pertinent parameters are discussed graphically in Figs.12–18.Skin friction increases with Hartmann numberHaand Reynolds numberM.Both parameters result in more resistance to fluid motion.Increasing behavior is observed as shown in Fig.12.Similarly porosity parameter also increases the drag force (see Fig.13).Effects ofNtandNbon heat transfer rate is discussed in Figs.15 and 16.Larger values ofNbandNtdecrease the temperature difference.That is whyNuxdecreases with both parameters.Similarly mass transfer rate for increasing values ofNbandNtare given in Figs.17 and 18.These figures show that mass transfer rate decreases with largerNbwhile it increases forNt.

        Fig.12 (Color online) Behaviors of M and Ha on Cfx.

        Fig.13 (Color online) Effect of β on Cfx.

        Fig.14 (Color online) Effect of λ on Cfx.

        Fig.15 (Color online) Behaviors of Pr and Nt on Nux.

        Fig.16 (Color online) Effect of Nb on Nux.

        Streamline for the current flow situation are sketched through Figs.20–23.The number of streamlines in these figures show the velocity of fluid.The more the streamlines, greater is the velocity of fluid.Figures 20–22 show that velocity decreases with Reynolds number(M),porosity parameter (λ) and Hartmann number (Ha) while it increases with injection parameter (β) (see Fig.23).

        Fig.17 (Color online) Effects of Nt and Sc on Shx.

        Fig.18 (Color online) Effects of Nb on Shx.

        In some limiting case, we draw the effect of magnetic parameter on skin friction in Fig.19.And the result is compared with Fig.5 in Ref.[20].Here we found a very good agreement.Therefore the present mathematical model and numerical method is verified.

        Fig.19 (Color online) Impact of magnetic parameter on skin friction.

        Fig.20 (Color online) Streamlines for M=1, 1.5, 3.

        Fig.21 (Color online) Streamlines for Ha=1, 2, 3.

        Fig.22 (Color online) Streamlines for λ=1, 2, 3.

        Fig.23 (Color online) Streamlines for β=1, 2, 5.

        6 Conclusions

        In this paper, we studied flow of viscous nanofluid between stretchable parallel plates in porous medium.The governing PDE’s system is transformed to dimensionless form by dimensionless parameters.For solution we used Numerical technique i.e.Finite difference scheme.From the results we obtained that velocity increases through Reynolds number and injection parameters while it reduces with magnetic and porosity parameters.Temperature (θ(τ,ξ,η))increases for larger Brownian variable (Nb), Prandtl number (Pr), and thermophoresis variable (Nt).Concentration(?(τ,ξ,η)) increases for larger Schmidt number and Brownian motion while it decays with thermophoresis parameter.Skin friction increases with Hartmann number, injection parameter and porosity parameter.Sherwood number increases withNtand decreases forNbandSc.

        亚洲国产精品无码久久久| 亚洲中文字幕高清视频| av一区二区三区综合网站| 中文字幕日韩精品有码视频| 无码少妇一区二区浪潮av| 久久精品国产亚洲不av麻豆| 亚洲视频在线中文字幕乱码| av在线观看免费天堂| 久久精品成人无码观看不卡| 国产午夜精品一区二区三区视频| 国产一区二区牛影视| 日本黄色高清视频久久| 蜜桃视频在线看一区二区三区| 欧美不卡一区二区三区| 日本免费人成视频播放| 久久er这里都是精品23| 久久精品国产69国产精品亚洲| 偷拍综合在线视频二区| 国产精品免费看久久久8 | 国产精品区一区二区三在线播放| 色综合视频一区中文字幕| 自拍偷拍亚洲一区| 国产在线精彩自拍视频| 国产精品会所一区二区三区| 久久精品免视看国产成人| 亚洲国产成人精品91久久久| 国产我不卡在线观看免费| 美女内射毛片在线看免费人动物| 成熟丰满熟妇高潮xxxxx视频| 亚洲成人免费网址| 成人免费毛片在线播放| 精品高朝久久久久9999| 国产精品无码一区二区三区免费| 国产婷婷丁香五月麻豆| 久久亚洲一区二区三区四区五| 亚洲中文字幕在线综合| 疯狂撞击丝袜人妻| 伊人婷婷色香五月综合缴激情| 国产av一区二区内射| 少妇真实被内射视频三四区| 夜夜爽一区二区三区精品|