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

    Effect of Magnetic Field on Forced Convection between Two Nanofluid Laminar Flows in a Channel

    2016-12-09 06:12:54AfrasiabRaisiandAhmadQanbary
    關(guān)鍵詞:節(jié)理繪制自動

    Afrasiab Raisi and Ahmad Qanbary

    ?

    Effect of Magnetic Field on Forced Convection between Two Nanofluid Laminar Flows in a Channel

    Afrasiab Raisi*and Ahmad Qanbary

    This paper provides a numerical study of forced convection between hot and cold nanofluid laminar flows that are separated by a thin membrane, in a horizontal channel. Outer surface of channels’ walls are thermally insulated and divide into two parts; namely NMP and MP. NMP is the channel’s wall from the entrance section to the middle section of channel that is not influenced by magnetic field. MP is the channel’s wall from the middle section to the exit section of channel which is influenced by a uniform-strength transverse magnetic field. The governing equations for both hot and cold flows are solved together using the SIMPLE algorithm. The effects of pertinent parameters, such as Reynolds number (10≤≤500), Hartman number (0≤≤60) and the solid volume fraction of copper nano-particles (0≤≤0.05), are studied. The results are reported in terms of streamlines, isotherms, velocity and temperature profiles and local and average Nusselt number. The results of the numerical simulation indicate that the increase in Reynolds number and the solid volume fraction lead to increase in Nusselt number. Meanwhile, the results also show that the rate of heat transfer between the flows increases as the Hartmann number increases, especially at higher values of the Reynolds number.

    magnetic field, forced convection, nanofluid, membrane, channel

    1 Introduction

    In order to save the energy consumption in various industries, including the automotive, electronics, aerospace, food processing and manufacturing industries, the heat exchangers that are used must be smaller and lighter and also with higher thermal performance. Due to restrictions on the size of the heat exchangers, the heat exchange surface cannot be increased. The use of fluids that have better heat transfer characteristics is a good way to increase the efficiency of such heat exchangers. Nano?uid is a suspension of solid nanoparticles in conventional base fluids such as water, oil or ethylene glycol, and has a higher thermal conductivity than the base fluid depending on the shape, size and type of the solid nanoparticles[1–5].

    In recent years, forced convection of nanofluids in different channels has been investigated by several researchers. Izadi, et al[6], numerically investigated laminar forced convection of alumina-water nanofluid in an annulus. Their findings showed that adding nanoparticles to the base fluid increased the convective heat transfer coefficient, but had a little effect on the velocity profile. Bianco, et al[7], numerically studied the developing laminar forced convection flow of a water–Al2O3nanofluid in a circular tube, submitted to a constant and uniform heat flux at the wall. They showed that the single and two-phase models led to the same results and the temperature dependent models presented higher values of heat transfer coefficient and Nusselt number. Santra, et al[8], examined laminar forced convective heat transfer of copper-water nanofluid in a horizontal parallel plate channel. For this purpose, the nanofluid was considered as Newtonian as well as non-Newtonian. They concluded that the thermal behavior of Newtonian and non- Newtonian nanofluid was more or less the same. Also, the sharp increase in the wall shear stress with increasing the solid volume fraction for non- Newtonian nanofluid, was one of their outcomes. Manca, et al[9], reported the results of an investigation on forced convection with Al2O3-water nanofluids in a two dimensional channel under uniform heat flux condition. In their report an increase in the heat transfer coefficient and required pumping power as solid volume fraction increases was confirmed.

    SAha, et al[10], presented a numerical investigation of turbulent forced convection flow of water based Al2O3and TiO2nanofluids flowing through a horizontal circular pipe under uniform heat flux boundary condition applied to the wall. Based on the results of this study, with increasing Reynolds number and solid volume fraction and decreasing size of nanoparticles, Nusselt number increased.

    In experimental research Peyghambarzadeh, et al[11]and Hussein, et al[12], experimentally examined the friction factor and forced convection heat transfer enhancement of a water based nanofluid with various nanoparticles in the automotive cooling system. In general, the results showed that the use of nanofluids increased the friction factor and forced convective heat transfer coefficient, so that the percentage increase in heat transfer coefficient was more than the percentage increase of the friction factor.

    In most of the existing equipment in various industries, a magnetic field is created due to the electric field. Thus the flow of electrically conducting fluid through a channel or a duct is exposed to a transverse magnetic field. The results of the research conducted have shown that the fluid flow subjected to a magnetic field, experience a Lorentz force. Lorentz force directly affects the flow field and consequently the heat transfer is also affected. Therefore, understanding of the flow and thermal behaviors of the channels that are cooled by electrically conducting fluids and are exposed to a magnetic field has long been of interest to researchers. Siegel[13]was among the first researchers who have studied the hydromagnetic flows. He studied convective heat transfer for the fully developed and laminar flow of an incompressible and electrically conducting fluid in a parallel plate channel submitted to a constant and uniform heat flux at the walls and exposed to a transverse magnetic field. The results showed that the flow field and the heat transfer in the channel were affected by the magnetic field. Alpher[14]continued and extended the earlier study done by Siegel[13]. In particular the earlier analysis of Siegel involving non-conducting plates were corrected and extended to plates of finite conductivity. Back[15]analytically studied the uniform laminar flow of an electrically conducting fluid in a constant cross-section, parallel-plate channel with isothermal, electrically nonconducting walls. In this study, the effect of the Joule heating parameter was also considered. The results showed that the thermal entrance lengths were proportional to Reynolds and Prandtl numbers, similar to the case of no applied magnetic field. Additionally, these lengths decreased with an increase in Joule heating parameter. Lahjomri, et al[16], analytically investigated the thermal developing laminar Hartmann flow through a parallel plate channel, including viscous dissipation, Joule heating and axial conduction with a step change in wall temperature. The results showed that in the absence of viscous dissipation and Joule heating, the Hartman effect led to a reduction of the bulk temperature. This study also showed that viscous dissipation and Joule heating played a significantly different role in the heat transfer, depending on the electrical condition of the walls. In another study, Lahjomri, et al[17], reiterated his earlier work under a constant heat flux boundary condition. This time, they showed that the heat transfer increased with the increasing Hartmann number until finally a saturated state was achieved.

    With the arrival of nanofluids to the technology arena, the interest to examine the hydro-magnetic flow of nanofluids has greatly increased. Aminossadati, et al[18], numerically examined the laminar forced convection of a water–Al2O3nanofluid flowing through a partially heated microchannel. They indicated that for any value of the Reynolds number, velocity and temperature along the centerline of the microchannel decreased as the Hartmann number increased. Their results also showed that the Nusselt number increased as the Reynolds and Hartman numbers increased. Mahian, et al[19], analytically investigated the effects of using TiO2/water nanofluid on the entropy generation in an isoflux vertical annulus under MHD flow effects. In this study the flow was assumed as laminar, steady and fully developed inside the annulus. The results showed that using TiO2/water nanofluid reduced the entropy generation in the annulus while an increase in the Hartmann number increased the entropy generation number. SyamSundar, et al[20], experimentally examined the forced convection heat transfer and friction factor in a tube with Fe3O4magnetic nanofluid. In this study, the heat transfer coefficient increased by 30.96% and friction factor by 10.01% at 0.6% volume concentration compared to flow of water at similar operating conditions. Ghofrani, et al[21], established an experimental investigation on laminar forced convection heat transfer of ferrofluids inside a circular copper tube under an alternating magnetic field. They investigated the effects of magnetic field, volume concentration and Reynolds number on the convective heat transfer coefficient. They found that increasing the frequency of the alternating magnetic field and the solid volume fraction, led to better heat transfer enhancement. Malvandi, et al[22], presented a theoretical study on the laminar flow and forced convective heat transfer of water/alumina nanofluid inside a parallel plate channel in the presence of a uniform magnetic field. They applied a modified two component model for water/alumina nanofluid and looked at the migration of nanoparticles. The results indicated that nanoparticles moved from the heated walls toward the core region of the channel and constructed a non-uniform nanoparticles distribution.

    It is worth mentioning that previous studies have solely investigated heat transfer in a single channel. However, in many engineering applications, such as heat exchangers, hot and cold fluids flow through a pair of adjacent channels, and forced convective heat transfer takes place between them. In fact, this system can be a simple model of a heat exchanger, which is exposed to a uniform magnetic field to increase the rate of heat transfer. The main aim of the present work is a numerical investigation of heat transfer performance of a Cu-water nanofluid in a pair of horizontal channels that is under the influence of a transverse magnetic field. It is assumed that the flows are separated by a thin plat (membrane) which its mass is negligible.

    2 Problem Description

    Fig. 1 shows the schematic diagram of a two-dimensional horizontal channel wherein hot and cold flows are separated by a membrane. Dimensionless height and length areandrespectively.

    Fig. 1. A schematic diagram of the physical model

    Both upper and lower walls of the channel are thermally insulated. The first part of the channel which is not influenced by magnetic field (Non-Magnetic part, i.e., NMP here after) has a dimensionless length ofand the second part influenced by magnetic field (magnetic part, i.e., MP from now on) has the same dimensionless length ofHot flow with velocityand temperatureand cold flow with velocityand temperatureenter lower and upper channels respectively. Both hot and cold nanofluid flows contain water and Cu nanoparticles and are assumed to be laminar, Newtonian and incompressible. Constant thermophysical properties are considered for the nanofluid, except for the density variation in the buoyancy forces that is determined by the Boussinesq approximation.

    The thermophysical properties of pure water (base fluid) and the Cu-nanoparticles, at a fixed temperature of 25℃, are presented in Table 1[23]. The Prandtl number of water is assumed to be

    Table 1. Thermophysical properties of water and Cu at 25℃

    3 Mathematical Formulation

    The non-dimensional governing equations of continuity, momentum, and energy for hot and cold flows can be written as follow (is for hot flow andis for cold flow:

    (2)

    (3)

    (5)

    The Reynolds, Hartman and Prandtl numbers are:

    The effective density, heat capacitance and the thermal diffusivity coefficient of the nanofluid are defined based on the properties of water and Cu[23]:

    (8)

    (9)

    The electrical conductivity of nanofluid was presented by Maxwell[24]as below:

    The effective dynamic viscosity of the nanofluid can be modeled by Brinkman[25]:

    (11)

    The thermal conductivity of the nanofluidkwas determined by Patel, et al[26]. For the two-component entity of spherical-particle suspension, the model gives:

    where

    (13)

    (15)

    In the above-mentioned equations,is the Boltzmann constantandis an experimental constant (=36 000) which was suggested by Patel, et al[26]. The molecular diameter of the solid nanoparticles is100 nm and the molecular diameter of the base fluid is2°A.

    The non-dimensional boundary conditions used to solve Eqs. (1)–(4) are as follows:

    Since the membrane is thin and its mass is negligible, the thermal boundary conditions on it are as follows:

    (17)

    4 Numerical Approach

    The non-dimensional governing equations for hot and cold flows (Eqs. (1)–(4)) with the boundary conditions given in Eqs. (16) and (17) are solved simultaneously by the control volume formulation using Patankar’s SIMPLE algorithm[27].

    The convection-diffusion terms are discretized by a power-law scheme, and the system is numerically modeled in FORTRAN. A regular rectangular domain with a uniform grid is used. The convergence criteria are to reduce the maximum mass residual of the grid control volume below 108.

    After solving both the hot and the cold governing equations forandother useful quantities such as the Nusselt number can be determined. The local Nusselt number () at the membrane can be expressed as:

    where 𝜆 is the convection heat transfer coefficient and is calculated as follows:

    (19)

    Thus the local Nusselt number by using the dimensionless parameters given in Eq. (5) is as follows:

    (21)

    Finally, the average Nusselt number () is determined by integratingalong the membrane separator:

    5 Grid Independency and Code Validation

    Table 2 indicates the solution grid independence for two different Hartmann numbers

    Table 2. Grid independency study at Re=100 f=0.03, D=1/2, VC1=1

    The average Nusselt number and the dimensionless bulk mean temperature of hot and cold flows at channels’ outlet, i.e., (1)outand (2)outare presented at different grid points when=0.03,1=1, and=1/2. The obtained results from Table 2 illustrates that the grid size of 81′400 satisfies grid independency.

    Fig. 2(a) compares the isotherms of this work with those of Santra, et al[8]to validate our code. Moreover, the variation of temperature along the centerline of the channel for=100 and=0.02 of Al2O3-Water Nanofluid is compared with that of Aminossadati, et al[18]in Fig. 2(b).

    Furthermore, in Fig. 2(c) the fully developed velocity profiles at three values of Hartman number contrasted with Back’s findings[15]presented by Eq. (23). Good agreements are seen in all Figs. 2(a), 2(b) and 2(c).

    Fig.2. Validation code against: (a) Isotherms inRef. [8], (b) Variation of temperature along the centerline of the channel inRef. [18], and (c) Analytical velocity profile modeled inRef. [15]

    6 Results

    The effects of the Reynolds number (10≤≤500), the Hartman number (0≤≤60) and the solid volume fraction (0≤≤0.05) on the thermal performance of the channel are studied and are presented in the following sections forand

    6.1 Effects of Reynolds and Hartman numbers

    In this section, the effects ofandon heat transfer, isotherms and streamlines are studied when=0.03. Streamlines and isotherms are shown in Fig. 3 at different values offor=0 and 40. In all values of, streamlines match together at NMP, whereas at=40 both the hot and the cold flows affected by magnetic field at MP where streamlines expand due to Lorentz force.

    Fig. 3. Streamlines (top) and isotherms (bottom) for Cu-water nanofluid (f=0.03) for different Re at Ha=40 and Ha=0

    Isotherms are not affected by the magnetic field at=10, and they are the same for=0 and 40. This is due to the small amounts of hot and cold flows’ velocities at=10, so both flows reach a balanced temperature very soon. Other findings show that heat transfer occurs throughout the channel length at=100 and 500. Thus the temperature patterns do penetrate into the exit section and at=40, isotherms are influenced by magnetic field at MP. The comparison between the isotherms of the nanofluid near the membrane shows that the heat transfer increases as theincreases.

    Fig. 4 presents the variations of dimensionless velocities along the centerline of hot or cold flows at=10 100 and 500. It is observed that both flows have the same velocity variations along the flow centerline in all values ofbecause of their similar geometric and hydrodynamic properties. In this figure, it can be seen that at=10 centerline velocity stabilizes much sooner in comparison with other values of. In addition, the length of hydrodynamic entrance increases as theincreases. However,maxdecreases due to the magnetic field effect by entering hot and cold flows to the MP. The gradient of this decrease depends on.

    Fig. 4. Variation of horizontal velocity along the centerlineof the channel for various Reynolds numbers atHa=40, f=0.03

    In Fig. 5,maxvariations are presented at four different values of Hartman number (=0, 20, 40 and 60) and=100. In the absence of magnetic field (=0), the maximum of horizontal flow velocity is constant at fully developed conditions, due to Lorentz force effect, however,maxdecreases at0. The highest decrease ofmaxachieves for the largest amount of(see Eq.(2)).

    Fig. 6 shows the velocity profiles at four different cross sections of channel; namely=7.5, 15, 22.5 and 30,for hot and cold nanofluid flows at=40 and=100. According to this figure, the velocity profile at=7.5 is parabolic and fully developed. The magnetic field starts at=15 where the shape of velocity profile begins to change. In two sections of=22.5 and 30, magnetic field yields a similar flat velocity profile for both sections.

    Reynolds and Hartman numbers’ effects onNuare shown in Figs. 7(a) and 7(b). Highercorresponds tohigherNu. Regardless of the values ofand, the maximum amount ofNuis obtained at the channel inlet due to the highest temperature gradient between hot and cold flows. After heat exchange between two nanofluid flows along the membrane length,Nustarts to decrease. For=10, two flows reach a thermal equilibrium immediately. Thus, in this case, magnetic field has no significant effect on isotherms’ variations andNu. For=100and 500, an increasing jump can be observed inNudiagram at the entrance of MP. This happens because in the existence of magnetic field, the velocity of two flows near the walls and membrane increases.

    節(jié)理等密圖的自動繪制………………………………………………………………………………………………李玉華(1.21)

    Fig. 5. Variation of horizontal velocity along the centerline of the channel for various Hartman numbers at=100,=0.03

    Fig. 6. Horizontal velocity profiles across various sections of channel at=0.03,=100,=40

    The effects ofon Nusselt number can be better understood by studying values presented in Table 3. As mentioned before, the increase inleads to the increase inNuat higher values of. After running the computer code for different values of the Reynolds and Hartmann numbers, the following correlation was found between the average Nusselt, Reynolds and Hartmann numbers when the solid volume fraction=0.03.

    where

    (25)

    Fig. 7. Variation of local Nusselt number along the channel length at=0.03, and at (a) various Reynolds numbers at=40, (b) various Hartman number at=100

    Table 3. Effects of Ha on the Num and its relative increase percentage (%Incr.) at various Reand f=0.03

    Note: %means relative increase percentage,

    6.2 Effect of solid volume fraction

    Figs. 8(a) and 8(b) show that the dimensionless bulk mean temperatures of the hot and the cold flows start to decrease and increase respectively as soon as they enter the channel. The rate of temperature change is higher at low values of, so that for=10, the dimensionless bulk mean temperatures of the two flows become equal before reaching the exit section for both the nanofluid and pure water while for=100 and 500 these temperatures are substantially different from each other. For all values of, the changes in dimensionless bulk mean temperature are more rapid for nanofluid flow than that of pure water. Figs. 8(a) and 8(b) also show that in all values of, the hot and the cold flows have the same hydrodynamic conditions, and variations of flows’ temperature along the channel are symmetrical.

    Fig. 8. Variation ofdimensionless balk mean temperature along the channel length for nanofluid (=0.03), and pure water at various Reynolds numbers and=40, for (a) hot flow and (b) cold flow

    Nuat different values of copper nanoparticles volume fraction (0≤≤0.05) for=0 and 40 is presented in Table 4 at different values of. In all cases, higher values ofNucan be achieved at higher amounts of. As a result, it can be said that the effects ofandon increasingNuare more significant for higher values of. For example, at=500 and=40, an increase of 3% in the volume fraction of copper nanoparticles results in approximately 18.84% increase in theNu, whereas this values is 1.0% for=10 and=40.

    After running the computer code for different values of the Reynolds numbers and solid volume fractions, the following correlation was found between the average Nusselt and Reynolds numbers and solid volume fractions when the Hatrmann number is=0 and 50≤≤500.

    where

    (27)

    Table 4. Effects of f on the Num and its relative increase percentage for different values of Re and Ha

    Note: %means relative increase percentage,

    7 Conclusions

    The flows were becoming fully developed later at higher values of the Reynolds number. Streamliners were expanded at MP due to applying the magnetic field. The Hartman number had no effect on heat transfer rate at=10 because two flows were reached an equal temperature immediately. This is while at higher values of Reynolds number, the average Nusselt number increased as the Hartman number increased. Moreover, a more noticeable increase in the Nusselt number with solid volume fraction was evident at higher values of Reynolds number.

    Nomenclature

    B0—Magnetic field strength, c—Experimental constant, Eq. (12), Cp—Specific heat, J ? kg–1? K–1, d—Distance between channel’s bottom wall and membrane, m, D—Dimensionless distance between channel’s bottom wall and membrane, D=dh, h—Channel height, m H—Dimensionless channel height, H=h/h=1, Ha—Hartman number, k—Thermal conductivity, W ? m–1? K–1, l—Length of channel, m, L—Dimensionless length of channel, L1—Dimensionless length of first region, without of magnetic field, L2—Dimensionless length of second region, with of magnetic field MP—Part of channel’s wall that is influenced by magnetic field, (magnetic part) NMP—Part of channel’s wall that is not influenced by magnetic field, (non-magnetic part), Nux—Local Nusselt number on the membrane, Num—Average Nusselt number along the membrane, p1—Hot flow pressure, Pa, p2—Cold flow pressure, Pa, —Modified pressure, P1—Dimensionless pressure for hot flow, P2—Dimensionless pressure for cold flow, Pe—Peclet number, Pe= Pr—Prandtl number, Pr= Re—Reynolds number, Re= T1, T2—Temperature of hot and cold flows respectively, K, Tc—Inlet temperature of cold flow, K, Th—Inlet temperature of hot flow, K, us—Brownian motion velocity, m/s, u, v—Velocity components in x, y directions, m/s, U1, V1—Dimensionless velocity of hot flow, U2, V2—Dimensionless velocity of cold flow, VC1—Velocity coefficient, x, y—Cartesian coordinates, , X, Y—Dimensionless coordinates, X=x/h, Y=y/h, a—Thermal diffusivity (k/rCp), m2/s, f—Solid volume fraction, kb—Boltzman constant, J/K, q1—Dimensionless temperature of hot flow, q2—Dimensionless temperature of cold flow, qm1—Dimensionless balk mean temperature of hot flow, qm3—Dimensionless balk mean temperature of cold flow, l—Heat transfer coefficient, W ? m–2? K–1, m—Dynamic viscosity, N ? s ? m-2, —Kinematic viscosity (m/r), m2/s, r—Density, kg ? m–3, s—Electrical conductivity, mS/cm. Subscripts C, 2—Cold, h, 1—Hot, eff—Effective, f—Pure fluid, nf—Nanofluid, out—Channel’s outlet, s—Nanoparticle.

    [1] Choi, SUS. Enhancing thermal conductivity of fluids with nanoparticles[J], 1995, 66: 99–105.

    [2] Eastman J A, Choi SUS, Li S, et al. Anomalously increased effective thermal conductivities of ethylene glycol-based nanofluids containing copper nanoparticles[J]., 2001, 78(6): 718–720.

    [3] Zhou D W. Heat transfer enhancement of copper nanofluid with acoustic cavitation[J]., 2004, 47: 3109–3117.

    [4] Moghadassi A R, MasoudHosseini S, Henneke D E. Effect of CuO nanoparticles in enhancing the thermal conductivities of mono ethylene glycol and paraffin fluids[J].2010, 49: 1900–1904.

    [5] Sahooli M, Sabbaghi S, Investigation of thermal properties of CuO nanoparticles on the ethylene glycol–water mixture[J]., 2013, 93: 254–257.

    [6] Izadi M, Behzadmehr A, Jalali-Vahida D. Numerical study of developing laminar forced convection of a nanofluid in an annulus[J]., 2009, 48(4): 2119–2129.

    [7] Bianco V, Chiacchio F, Manca O, et al. Numerical investigation of nanofluids forced convection in circular tubes[J]., 2009, 29: 3632–3642.

    [8] Santra A K, Sen S, Chakraborty N. Study of heat transfer due to laminarflow of copper–water nanofluid through two isothermally heated parallel plates[J]., 2009, 48: 391–400.

    [9] Manca O, Nardini S, Ricci D. A numerical study of nanofluid forced convection in ribbed channels[J].,2012, 37: 280–292.

    [10] Saha G, Paul CM. Numerical analysis of the heat transfer behaviour of water based Al2O3and TiO2nanofluids in a circular pipe under the turbulent flow condition[J]., 2014, 56: 96–108.

    [11] Peyghambarzadeh S M, Hashemabadi S H, SeifiJamnani M, et al. Improving the cooling performance of automobile radiator with Al2O3/water nanofluid[J].2011, 31: 1833–1838.

    [12] Hussein A M, Bakar R A, Kadirgama K. Study of forced convection nanofluid heat transfer in the automotive cooling system[J]., 2014, 2: 50–61.

    [13] Siegel R. Effect of magnetic field on forced convection heat transfer in a parallel plate channel[J]., 1958, 25: 415–420.

    [14] Alpher R A. Heat transfer in magneto hydrodynamic flow between parallel plates[J].,1961, 3(2): 108–112.

    [15] Back L H. Laminar heat transfer in electrically conducting fluids flowing in parallel plate channels[J]., 1968, 11: 1621–1636.

    [16] Lahjomri J, Oubarra A, Alemany A. Heat transfer by laminar Hartmann flow in thermal entrance region with a step change in wall temperatures: the Graetz problem extended[J]., 2002, 45(5): 1127–1148.

    [17] Lahjomri J, Zniber K, Oubarra A, et al. Heat transfer by laminar Hartmann’s flow in thermal entrance region with uniform wall heat flux: the Graetz problem extended[J].,2003, 44(1): 11–34.

    [18] Aminossadati S M, Raisi A, Ghasemi B. Effects of magnetic field on nanofluid forced convection in a partially heated microchannel., 2011, 46: 1373–1382.

    [19] Mahian O, Pop I, Sahin Z A, et al. Irreversibility analysis of a vertical annulus using TiO2/water nanofluid with MHD flow effects[J]., 2013, 64: 671–679.

    [20] SyamSundar L, Naik M T, Sharma K V, et al. Experimental investigation of forced convection heat transfer and friction factor in a tube with Fe3O4magnetic nanofluid[J]., 2011, 37: 65–71.

    [21] Ghofrani A, Dibaei M H, Hakim-Sima A, et al. Experimental investigation on laminar forced convection heat transfer of ferro-fluids under an alternating magnetic field[J]., 2013, 49: 193–200.

    [22] Malvandi A, Ganji D D. Magnetic field effect on nanoparticles migration and heat transfer of water/alumina nanofluid in a channel [J].,2014, 362: 172–179.

    [23] Raisi A, Ghasemi B, Aminossadati S M. A numerical study on the forced convection of laminar nanofluid in a microchannel with both slip and no-slip conditions[J].,, 2011, 59(2): 114–129.

    [24] Sheikholeslami M, Soleimani S, Gorji-Bandpy M, et al. Natural Convection of nanofluids in an enclosure between a circular and a sinusoidal cylinder in the presence of magnetic field[J]., 2012, 39: 1435–1443.

    [25] Brinkman H C. The viscosity of concentrated suspensions and solution[J]., 1952, 20: 571–581.

    [26] Patel H E, Sundararajan T, Pradeep T, et al. A micro-convection model for thermal conductivity of nanofluids[J]., 2005, 65(5): 863–869.

    [27] Patankar S V.[M]. New York: Hemisphere Publishing Corporation-Taylor and Francis Group, 1980.

    Biographical notes

    Afrasiab Raisi was born in Iran and currently is an associate professor atHe received his B. Eng. Degree in Mechanical Engineering fromin 1992, and the M.Sc. and Ph.D. degrees in Mechanical Engineering from, in 1995 and 2001, respectively. In 2001, he joined theas an assistant Professor, became an associate professor in 2015. His current research interests include forced and natural convection heat transfer, convective heat transfer in nanofluids and non-Newtonian fluids, melting and solidification and conduction microscale heat transfer.

    Tel: +98-3832324435;

    E-mail: afraisi@gmail.com; raisi@eng.sku.ac.ir

    Ahmad Qanbary was born in Iran in 1980 and currently is an engineer atHe received his M.Sc. degree in Mechanical Engineering (CFD branch) fromin 2013.

    E-mail: ahmadqanbary@yahoo.com

    Received August 17, 2015; revised May 30, 2016; accepted June 15, 2016

    ? Chinese Mechanical Engineering Society and Springer-Verlag Berlin Heidelberg 2016

    10.3901/CJME.2016.0615.072, available online at www.springerlink.com; www.cjmenet.com

    E-mail: raisi@eng.sku.ac.ir

    猜你喜歡
    節(jié)理繪制自動
    Art on coffee cups
    順傾節(jié)理邊坡開挖軟材料模型實驗設(shè)計與分析
    新疆阜康白楊河礦區(qū)古構(gòu)造應(yīng)力場特征
    自動捕盜機
    放學(xué)后
    童話世界(2018年17期)2018-07-30 01:52:02
    新疆阜康白楊河礦區(qū)構(gòu)造節(jié)理發(fā)育特征
    中國煤層氣(2018年3期)2018-07-27 11:41:42
    基于STM32的自動喂養(yǎng)機控制系統(tǒng)
    電子測試(2018年10期)2018-06-26 05:53:36
    關(guān)于自動駕駛
    汽車博覽(2016年9期)2016-10-18 13:05:41
    Stefan Greiner:我們?yōu)槭裁葱枰詣玉{駛?
    在轉(zhuǎn)變中繪制新藍圖
    成人免费观看视频高清| 亚洲avbb在线观看| 一级毛片女人18水好多| 亚洲中文字幕日韩| 在线永久观看黄色视频| 久久精品成人免费网站| 国产精品欧美亚洲77777| 久久精品亚洲精品国产色婷小说| av线在线观看网站| 黄色视频,在线免费观看| 日韩人妻精品一区2区三区| 两人在一起打扑克的视频| 国产日韩一区二区三区精品不卡| 又紧又爽又黄一区二区| 久久狼人影院| 色94色欧美一区二区| 国产精品久久久久久人妻精品电影 | 久久这里只有精品19| 最近最新中文字幕大全电影3 | 亚洲欧美激情在线| 人妻一区二区av| 亚洲一区二区三区欧美精品| 久久影院123| 老司机在亚洲福利影院| 国产精品久久久人人做人人爽| 国产免费福利视频在线观看| 欧美日韩亚洲国产一区二区在线观看 | 一级毛片电影观看| 一二三四在线观看免费中文在| 巨乳人妻的诱惑在线观看| 久久人妻福利社区极品人妻图片| 性少妇av在线| 女性被躁到高潮视频| 精品亚洲成a人片在线观看| 亚洲精品粉嫩美女一区| aaaaa片日本免费| 真人做人爱边吃奶动态| 精品亚洲乱码少妇综合久久| 黄频高清免费视频| 亚洲欧洲日产国产| 丰满饥渴人妻一区二区三| 国产一卡二卡三卡精品| 超色免费av| 亚洲精品国产色婷婷电影| 亚洲精品一卡2卡三卡4卡5卡| 国产精品麻豆人妻色哟哟久久| 成在线人永久免费视频| 精品欧美一区二区三区在线| av不卡在线播放| 国产亚洲精品第一综合不卡| 大型av网站在线播放| 欧美一级毛片孕妇| 国产在视频线精品| 久久精品熟女亚洲av麻豆精品| 亚洲国产精品一区二区三区在线| 美女扒开内裤让男人捅视频| 一级黄色大片毛片| 国产不卡av网站在线观看| 久久精品91无色码中文字幕| 久久久久精品国产欧美久久久| 欧美精品啪啪一区二区三区| 在线观看免费视频网站a站| 大型av网站在线播放| 欧美一级毛片孕妇| 中文亚洲av片在线观看爽 | 视频在线观看一区二区三区| 最近最新中文字幕大全电影3 | 日韩免费av在线播放| 国产一区二区三区综合在线观看| 国产一区二区在线观看av| 精品卡一卡二卡四卡免费| 97在线人人人人妻| 国产精品久久久久久人妻精品电影 | 日韩免费av在线播放| 五月开心婷婷网| 久久久精品国产亚洲av高清涩受| 亚洲综合色网址| 国产成+人综合+亚洲专区| 国产亚洲精品第一综合不卡| av超薄肉色丝袜交足视频| 高清在线国产一区| 婷婷成人精品国产| 色94色欧美一区二区| 性少妇av在线| 水蜜桃什么品种好| 99re6热这里在线精品视频| 最新的欧美精品一区二区| 法律面前人人平等表现在哪些方面| 国产成人免费无遮挡视频| 9热在线视频观看99| 老汉色∧v一级毛片| 男女无遮挡免费网站观看| 纵有疾风起免费观看全集完整版| 成人18禁高潮啪啪吃奶动态图| 麻豆乱淫一区二区| 啦啦啦 在线观看视频| 青草久久国产| 在线观看免费高清a一片| 18禁观看日本| 一夜夜www| 久久久久国产一级毛片高清牌| 国产精品国产高清国产av | a级毛片在线看网站| 亚洲精品国产色婷婷电影| 日韩中文字幕欧美一区二区| 高清视频免费观看一区二区| 老司机午夜十八禁免费视频| 曰老女人黄片| 亚洲精品美女久久久久99蜜臀| 久久久国产成人免费| 一区二区三区国产精品乱码| 免费看a级黄色片| 狠狠婷婷综合久久久久久88av| 成人三级做爰电影| 欧美亚洲日本最大视频资源| 日韩欧美免费精品| 高清av免费在线| 国产精品免费视频内射| 国产三级黄色录像| 热99re8久久精品国产| 香蕉久久夜色| 一区二区日韩欧美中文字幕| 99精品久久久久人妻精品| 丝袜美足系列| 日韩欧美一区视频在线观看| 精品熟女少妇八av免费久了| www日本在线高清视频| h视频一区二区三区| 久久久久国内视频| 国产精品电影一区二区三区 | 中文字幕另类日韩欧美亚洲嫩草| 成人影院久久| 成人国产一区最新在线观看| 母亲3免费完整高清在线观看| 欧美大码av| 成人av一区二区三区在线看| 狠狠狠狠99中文字幕| 欧美成人免费av一区二区三区 | 国产高清视频在线播放一区| 丝袜美腿诱惑在线| 国产精品欧美亚洲77777| 999久久久精品免费观看国产| 亚洲精品国产精品久久久不卡| 色婷婷av一区二区三区视频| 午夜两性在线视频| 香蕉丝袜av| 多毛熟女@视频| 999久久久精品免费观看国产| 久久精品人人爽人人爽视色| av一本久久久久| 日本黄色视频三级网站网址 | 亚洲人成伊人成综合网2020| 一级a爱视频在线免费观看| 亚洲成av片中文字幕在线观看| 热99re8久久精品国产| 在线十欧美十亚洲十日本专区| 80岁老熟妇乱子伦牲交| 三级毛片av免费| 老司机亚洲免费影院| 另类亚洲欧美激情| 色婷婷av一区二区三区视频| 成人特级黄色片久久久久久久 | 精品一区二区三区av网在线观看 | 两性夫妻黄色片| 亚洲av成人一区二区三| 欧美老熟妇乱子伦牲交| 国产成人免费观看mmmm| 久久久精品区二区三区| 久久人妻福利社区极品人妻图片| 9色porny在线观看| 一个人免费看片子| 国产色视频综合| 在线观看舔阴道视频| 久久亚洲精品不卡| 国产一区二区在线观看av| 成年版毛片免费区| 亚洲精品粉嫩美女一区| av一本久久久久| 午夜精品国产一区二区电影| 国产成人精品无人区| 黑人欧美特级aaaaaa片| 精品一区二区三卡| 99久久国产精品久久久| 亚洲自偷自拍图片 自拍| 丰满饥渴人妻一区二区三| 精品视频人人做人人爽| 久久久久久久精品吃奶| 一区二区三区精品91| 菩萨蛮人人尽说江南好唐韦庄| 免费久久久久久久精品成人欧美视频| 久久久久精品国产欧美久久久| av不卡在线播放| 在线观看人妻少妇| 老汉色av国产亚洲站长工具| 这个男人来自地球电影免费观看| 亚洲精品国产色婷婷电影| 中文亚洲av片在线观看爽 | 久久精品国产亚洲av香蕉五月 | 久久久久久久久久久久大奶| 国产97色在线日韩免费| 啦啦啦在线免费观看视频4| 久久久久久免费高清国产稀缺| 欧美变态另类bdsm刘玥| 亚洲精品中文字幕一二三四区 | 三上悠亚av全集在线观看| 国产伦理片在线播放av一区| 久久精品91无色码中文字幕| 国产人伦9x9x在线观看| 精品国产超薄肉色丝袜足j| 久久精品国产亚洲av香蕉五月 | 日本av免费视频播放| 亚洲 欧美一区二区三区| 亚洲欧美色中文字幕在线| 日本av免费视频播放| 老司机深夜福利视频在线观看| 亚洲男人天堂网一区| 十分钟在线观看高清视频www| 中文亚洲av片在线观看爽 | 欧美黄色淫秽网站| 成年版毛片免费区| 国产免费av片在线观看野外av| av超薄肉色丝袜交足视频| 久久久精品国产亚洲av高清涩受| 亚洲欧美日韩高清在线视频 | 亚洲五月婷婷丁香| 夫妻午夜视频| 黄色成人免费大全| 久久毛片免费看一区二区三区| 亚洲精品自拍成人| 久久久久精品国产欧美久久久| 天堂俺去俺来也www色官网| 精品一区二区三卡| 人成视频在线观看免费观看| 色综合欧美亚洲国产小说| 激情在线观看视频在线高清 | 极品教师在线免费播放| 精品国产乱子伦一区二区三区| 巨乳人妻的诱惑在线观看| 欧美日韩一级在线毛片| 下体分泌物呈黄色| 狠狠婷婷综合久久久久久88av| 50天的宝宝边吃奶边哭怎么回事| 久久毛片免费看一区二区三区| 亚洲成a人片在线一区二区| 激情视频va一区二区三区| 黄频高清免费视频| 91九色精品人成在线观看| kizo精华| 多毛熟女@视频| 十八禁高潮呻吟视频| 大香蕉久久成人网| 欧美激情久久久久久爽电影 | 欧美 亚洲 国产 日韩一| 亚洲精品美女久久av网站| √禁漫天堂资源中文www| av天堂在线播放| 久久性视频一级片| 久热爱精品视频在线9| 极品教师在线免费播放| 国产成人精品久久二区二区免费| 成年人免费黄色播放视频| 女人精品久久久久毛片| 欧美乱码精品一区二区三区| 午夜福利视频在线观看免费| 中文字幕另类日韩欧美亚洲嫩草| 大码成人一级视频| 成年人午夜在线观看视频| 精品人妻在线不人妻| 最新的欧美精品一区二区| 青青草视频在线视频观看| 亚洲全国av大片| 叶爱在线成人免费视频播放| 老司机午夜福利在线观看视频 | 久久久欧美国产精品| 欧美精品人与动牲交sv欧美| www.精华液| 国产一区二区三区视频了| 久久久国产精品麻豆| 99久久精品国产亚洲精品| 岛国在线观看网站| 欧美黑人精品巨大| 这个男人来自地球电影免费观看| 下体分泌物呈黄色| 国产福利在线免费观看视频| 91九色精品人成在线观看| 老汉色av国产亚洲站长工具| 天堂8中文在线网| 亚洲精品久久成人aⅴ小说| 性少妇av在线| 精品亚洲乱码少妇综合久久| 久久久精品国产亚洲av高清涩受| 天天躁夜夜躁狠狠躁躁| 在线观看免费午夜福利视频| 老汉色∧v一级毛片| 一边摸一边抽搐一进一小说 | 99re在线观看精品视频| 国产97色在线日韩免费| 精品亚洲成a人片在线观看| 99热国产这里只有精品6| 亚洲成人免费电影在线观看| 亚洲一区二区三区欧美精品| 一二三四在线观看免费中文在| 国产免费福利视频在线观看| av不卡在线播放| 久久国产精品大桥未久av| 精品国内亚洲2022精品成人 | av有码第一页| 丁香六月欧美| 欧美激情极品国产一区二区三区| 国产一区二区 视频在线| 久久亚洲精品不卡| 国产激情久久老熟女| 免费黄频网站在线观看国产| 国产区一区二久久| 久久久国产欧美日韩av| 成人国语在线视频| 免费黄频网站在线观看国产| 免费不卡黄色视频| 国产男靠女视频免费网站| 在线观看免费午夜福利视频| 狠狠精品人妻久久久久久综合| 久久精品成人免费网站| 国产视频一区二区在线看| 亚洲av欧美aⅴ国产| 久久精品亚洲熟妇少妇任你| 久久人妻熟女aⅴ| 日日夜夜操网爽| 久久中文看片网| 国产精品欧美亚洲77777| tube8黄色片| 亚洲成a人片在线一区二区| 欧美+亚洲+日韩+国产| 欧美日韩视频精品一区| 欧美成人午夜精品| 日韩欧美免费精品| 欧美久久黑人一区二区| 999精品在线视频| 精品人妻熟女毛片av久久网站| 99久久精品国产亚洲精品| 水蜜桃什么品种好| 国产一区二区激情短视频| 老司机福利观看| 久久久久久免费高清国产稀缺| 日韩欧美免费精品| 精品第一国产精品| 在线观看人妻少妇| 超色免费av| 亚洲午夜理论影院| 一进一出好大好爽视频| 国产精品秋霞免费鲁丝片| 国产精品一区二区精品视频观看| 国产午夜精品久久久久久| 99久久精品国产亚洲精品| 国产亚洲一区二区精品| 日韩欧美免费精品| 50天的宝宝边吃奶边哭怎么回事| 99久久人妻综合| 久久精品成人免费网站| 热re99久久国产66热| 一区二区日韩欧美中文字幕| 热99re8久久精品国产| 水蜜桃什么品种好| 丝袜人妻中文字幕| 精品少妇内射三级| 一进一出好大好爽视频| 性少妇av在线| 国产精品亚洲一级av第二区| 大陆偷拍与自拍| 在线播放国产精品三级| 97在线人人人人妻| 青草久久国产| 国产老妇伦熟女老妇高清| 最近最新中文字幕大全免费视频| 午夜激情av网站| 狠狠狠狠99中文字幕| 亚洲中文日韩欧美视频| 国产一卡二卡三卡精品| 欧美成人免费av一区二区三区 | 性少妇av在线| 一边摸一边抽搐一进一出视频| 19禁男女啪啪无遮挡网站| 午夜成年电影在线免费观看| 韩国精品一区二区三区| 日韩大码丰满熟妇| netflix在线观看网站| 国产不卡av网站在线观看| 老司机午夜十八禁免费视频| 亚洲三区欧美一区| 极品教师在线免费播放| 黄片播放在线免费| 国产真人三级小视频在线观看| av福利片在线| 国产精品98久久久久久宅男小说| 老司机影院毛片| 黑人欧美特级aaaaaa片| 一区二区日韩欧美中文字幕| 制服人妻中文乱码| 亚洲久久久国产精品| 看免费av毛片| videosex国产| 老司机深夜福利视频在线观看| 亚洲国产精品一区二区三区在线| 亚洲欧美精品综合一区二区三区| 国产精品偷伦视频观看了| 婷婷成人精品国产| 亚洲人成电影免费在线| 涩涩av久久男人的天堂| 午夜激情久久久久久久| 一进一出好大好爽视频| 看免费av毛片| 久久国产精品大桥未久av| 女人爽到高潮嗷嗷叫在线视频| 亚洲九九香蕉| 成年版毛片免费区| 成人手机av| 亚洲国产成人一精品久久久| 99久久精品国产亚洲精品| 少妇精品久久久久久久| 亚洲国产中文字幕在线视频| 精品乱码久久久久久99久播| av国产精品久久久久影院| 久久精品人人爽人人爽视色| 老司机深夜福利视频在线观看| 亚洲一码二码三码区别大吗| 亚洲精品成人av观看孕妇| 久久精品亚洲av国产电影网| 好男人电影高清在线观看| 夜夜夜夜夜久久久久| 丁香欧美五月| 中文字幕最新亚洲高清| 男女高潮啪啪啪动态图| 日韩 欧美 亚洲 中文字幕| 色尼玛亚洲综合影院| 精品乱码久久久久久99久播| 在线观看免费高清a一片| 叶爱在线成人免费视频播放| 精品福利观看| 欧美日韩视频精品一区| 极品少妇高潮喷水抽搐| 欧美另类亚洲清纯唯美| 久久av网站| 9色porny在线观看| 国产精品一区二区在线观看99| 夜夜骑夜夜射夜夜干| 69av精品久久久久久 | 国产国语露脸激情在线看| 女人高潮潮喷娇喘18禁视频| 久久久国产一区二区| 他把我摸到了高潮在线观看 | 中文字幕人妻丝袜制服| 黄色成人免费大全| 美女主播在线视频| 天天躁狠狠躁夜夜躁狠狠躁| 夫妻午夜视频| 欧美日本中文国产一区发布| 欧美乱码精品一区二区三区| 亚洲国产欧美一区二区综合| 免费在线观看影片大全网站| 最新在线观看一区二区三区| 精品人妻在线不人妻| 91老司机精品| 黄色视频,在线免费观看| 亚洲一区中文字幕在线| 国产精品成人在线| 欧美+亚洲+日韩+国产| 老司机午夜福利在线观看视频 | 日韩中文字幕视频在线看片| 国产精品久久久人人做人人爽| 黄色怎么调成土黄色| 日韩视频一区二区在线观看| 曰老女人黄片| 国产欧美日韩一区二区精品| 无遮挡黄片免费观看| 国产亚洲午夜精品一区二区久久| 天堂动漫精品| 精品久久久精品久久久| 1024视频免费在线观看| 在线亚洲精品国产二区图片欧美| 精品欧美一区二区三区在线| 免费在线观看完整版高清| 99国产精品99久久久久| 满18在线观看网站| 亚洲精品国产区一区二| 99久久精品国产亚洲精品| 精品一品国产午夜福利视频| 国产99久久九九免费精品| 日韩视频一区二区在线观看| 九色亚洲精品在线播放| 99精品久久久久人妻精品| 一个人免费看片子| 最新美女视频免费是黄的| 国产精品九九99| 亚洲va日本ⅴa欧美va伊人久久| a在线观看视频网站| 自线自在国产av| 777米奇影视久久| 欧美老熟妇乱子伦牲交| 精品国产国语对白av| 一级片免费观看大全| 午夜视频精品福利| 操美女的视频在线观看| 国产免费av片在线观看野外av| 免费在线观看日本一区| 久久中文字幕一级| 欧美精品亚洲一区二区| 欧美日韩福利视频一区二区| 丁香欧美五月| 黄色怎么调成土黄色| 中文字幕另类日韩欧美亚洲嫩草| 亚洲熟女毛片儿| 老司机午夜十八禁免费视频| 男女无遮挡免费网站观看| 国产在线精品亚洲第一网站| 欧美日韩亚洲综合一区二区三区_| 又紧又爽又黄一区二区| 亚洲人成伊人成综合网2020| 国产免费av片在线观看野外av| 国产1区2区3区精品| 交换朋友夫妻互换小说| 免费在线观看影片大全网站| 色综合欧美亚洲国产小说| 99re6热这里在线精品视频| 一区二区三区国产精品乱码| 日韩视频一区二区在线观看| av不卡在线播放| 美女视频免费永久观看网站| 一个人免费在线观看的高清视频| 亚洲午夜理论影院| 成人亚洲精品一区在线观看| 亚洲国产欧美网| 男女之事视频高清在线观看| 成年版毛片免费区| 高清毛片免费观看视频网站 | 欧美激情高清一区二区三区| 美女福利国产在线| 日韩欧美免费精品| 丰满迷人的少妇在线观看| 亚洲第一青青草原| 高潮久久久久久久久久久不卡| 国产在线视频一区二区| 精品少妇久久久久久888优播| 捣出白浆h1v1| 黄色毛片三级朝国网站| 黄色视频不卡| 丰满少妇做爰视频| 天堂中文最新版在线下载| 亚洲专区字幕在线| 精品福利永久在线观看| 午夜两性在线视频| av福利片在线| www.自偷自拍.com| aaaaa片日本免费| 日韩成人在线观看一区二区三区| 国产亚洲欧美精品永久| 另类亚洲欧美激情| 欧美国产精品一级二级三级| 日韩熟女老妇一区二区性免费视频| 天堂8中文在线网| 精品少妇黑人巨大在线播放| 啦啦啦视频在线资源免费观看| 在线天堂中文资源库| 久久精品国产综合久久久| 亚洲av第一区精品v没综合| 国产精品1区2区在线观看. | 国产欧美日韩精品亚洲av| 啦啦啦视频在线资源免费观看| 久久中文字幕人妻熟女| 色老头精品视频在线观看| 亚洲黑人精品在线| 在线观看免费视频日本深夜| 嫁个100分男人电影在线观看| 激情视频va一区二区三区| 日韩有码中文字幕| 国产高清视频在线播放一区| 性少妇av在线| avwww免费| 国产成人啪精品午夜网站| 国产精品久久久人人做人人爽| xxxhd国产人妻xxx| 国产av精品麻豆| 成年版毛片免费区| 男人操女人黄网站| 国产日韩欧美亚洲二区| 国产亚洲欧美精品永久| 国产精品免费一区二区三区在线 | 国产高清videossex| 99在线人妻在线中文字幕 | 精品一区二区三区视频在线观看免费 | 热re99久久国产66热| 最近最新中文字幕大全免费视频| 日韩一区二区三区影片| 久久久精品国产亚洲av高清涩受| 午夜福利欧美成人| 黄色怎么调成土黄色| 人人妻人人澡人人爽人人夜夜| 日韩免费高清中文字幕av| 久久人妻熟女aⅴ| 在线观看66精品国产| 麻豆国产av国片精品| 伊人久久大香线蕉亚洲五| 十分钟在线观看高清视频www| 久久精品aⅴ一区二区三区四区| 人人澡人人妻人| 亚洲av电影在线进入| 岛国在线观看网站| 桃红色精品国产亚洲av| 亚洲 国产 在线| 欧美一级毛片孕妇| 欧美黄色片欧美黄色片| 中文字幕高清在线视频| 中文字幕另类日韩欧美亚洲嫩草| 在线观看66精品国产| 在线观看免费午夜福利视频| tube8黄色片| 久久午夜综合久久蜜桃|