• <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)變中繪制新藍圖
    亚州av有码| 中出人妻视频一区二区| 精品一区二区三区视频在线| 久久午夜福利片| 午夜爱爱视频在线播放| 卡戴珊不雅视频在线播放| 日韩强制内射视频| 波多野结衣巨乳人妻| 久久国内精品自在自线图片| 国产精品久久久久久亚洲av鲁大| 在线天堂最新版资源| 免费黄网站久久成人精品| 一级毛片我不卡| 校园人妻丝袜中文字幕| 狂野欧美激情性xxxx在线观看| 国产成年人精品一区二区| 在线观看av片永久免费下载| 99热这里只有是精品在线观看| 久久韩国三级中文字幕| 欧美日韩乱码在线| 淫秽高清视频在线观看| 美女黄网站色视频| 你懂的网址亚洲精品在线观看 | 哪里可以看免费的av片| 久久久精品欧美日韩精品| 免费观看人在逋| 在线观看午夜福利视频| 亚洲欧美清纯卡通| 久久人人精品亚洲av| 国产 一区 欧美 日韩| 国产黄色视频一区二区在线观看 | 国产毛片a区久久久久| 高清毛片免费看| 麻豆国产av国片精品| 男人的好看免费观看在线视频| 欧美成人a在线观看| 日本免费a在线| 亚洲国产精品成人综合色| 成人午夜高清在线视频| 久久精品国产清高在天天线| 国产精品国产高清国产av| 中文资源天堂在线| 99久久人妻综合| 午夜免费激情av| 国产成人a∨麻豆精品| 午夜福利视频1000在线观看| 亚洲成人久久性| 国产免费男女视频| 十八禁国产超污无遮挡网站| 国产大屁股一区二区在线视频| 国产精品麻豆人妻色哟哟久久 | av黄色大香蕉| 欧美xxxx性猛交bbbb| 春色校园在线视频观看| 91av网一区二区| 色综合色国产| 18+在线观看网站| 亚洲最大成人中文| 内射极品少妇av片p| 婷婷精品国产亚洲av| 欧美高清成人免费视频www| 高清毛片免费看| 欧美成人精品欧美一级黄| 天美传媒精品一区二区| 亚洲欧美日韩卡通动漫| 亚洲在久久综合| 99在线视频只有这里精品首页| 国产成人影院久久av| 久久九九热精品免费| 久久精品国产亚洲网站| 禁无遮挡网站| 少妇猛男粗大的猛烈进出视频 | 亚洲一区高清亚洲精品| 日韩 亚洲 欧美在线| 一进一出抽搐动态| 99久久精品国产国产毛片| 能在线免费观看的黄片| 给我免费播放毛片高清在线观看| 性插视频无遮挡在线免费观看| 69av精品久久久久久| 国产精品久久视频播放| 禁无遮挡网站| av在线播放精品| 卡戴珊不雅视频在线播放| 亚洲欧美清纯卡通| 最近的中文字幕免费完整| 国产一级毛片在线| 欧美另类亚洲清纯唯美| 国产av一区在线观看免费| 桃色一区二区三区在线观看| 久久久久性生活片| 国产一级毛片在线| 欧美色视频一区免费| 麻豆成人午夜福利视频| 少妇人妻精品综合一区二区 | 少妇被粗大猛烈的视频| 女同久久另类99精品国产91| 你懂的网址亚洲精品在线观看 | 老司机福利观看| 性欧美人与动物交配| 亚洲欧美清纯卡通| 国产精品一及| 男女做爰动态图高潮gif福利片| 变态另类丝袜制服| 麻豆国产97在线/欧美| 一级毛片电影观看 | 亚洲18禁久久av| av国产免费在线观看| 变态另类丝袜制服| 欧美zozozo另类| 国产成人a∨麻豆精品| 日韩欧美精品免费久久| 亚洲国产色片| 少妇熟女欧美另类| 国产精品永久免费网站| 欧美极品一区二区三区四区| 亚洲欧美日韩高清在线视频| 97超碰精品成人国产| 人人妻人人看人人澡| 91精品国产九色| 大又大粗又爽又黄少妇毛片口| 在线天堂最新版资源| 久久精品人妻少妇| 久久精品国产自在天天线| 久久久久久久久大av| 日韩人妻高清精品专区| 国产爱豆传媒在线观看| 18禁黄网站禁片免费观看直播| 国产伦精品一区二区三区视频9| 国产一区亚洲一区在线观看| 搡老妇女老女人老熟妇| 国产亚洲91精品色在线| 国产av不卡久久| 噜噜噜噜噜久久久久久91| 好男人视频免费观看在线| 韩国av在线不卡| 精品99又大又爽又粗少妇毛片| 最近最新中文字幕大全电影3| 亚洲不卡免费看| 欧美成人免费av一区二区三区| 日本三级黄在线观看| 最近最新中文字幕大全电影3| 国产探花在线观看一区二区| 一本精品99久久精品77| 女的被弄到高潮叫床怎么办| 中文字幕制服av| 欧美一级a爱片免费观看看| 99热这里只有是精品50| 国产激情偷乱视频一区二区| 久久久久免费精品人妻一区二区| 亚洲精华国产精华液的使用体验 | 国产伦一二天堂av在线观看| 国产成人一区二区在线| 一个人看的www免费观看视频| 国产精品.久久久| 不卡一级毛片| 久久久午夜欧美精品| 99热只有精品国产| 有码 亚洲区| 性插视频无遮挡在线免费观看| 亚洲精品久久国产高清桃花| 国产人妻一区二区三区在| 亚洲久久久久久中文字幕| 午夜久久久久精精品| 亚洲国产日韩欧美精品在线观看| 国产伦一二天堂av在线观看| 日本熟妇午夜| 婷婷亚洲欧美| 不卡一级毛片| 久久久成人免费电影| 91av网一区二区| 国产乱人视频| 97人妻精品一区二区三区麻豆| 99热网站在线观看| 床上黄色一级片| 国产蜜桃级精品一区二区三区| 国国产精品蜜臀av免费| 女人十人毛片免费观看3o分钟| 国产老妇伦熟女老妇高清| 日本黄色片子视频| 精品99又大又爽又粗少妇毛片| 直男gayav资源| 晚上一个人看的免费电影| 国产精品伦人一区二区| 日韩成人av中文字幕在线观看| 国产精品美女特级片免费视频播放器| 久久人人爽人人片av| 色视频www国产| 亚洲人与动物交配视频| 欧美色欧美亚洲另类二区| 人妻少妇偷人精品九色| 亚洲图色成人| 99热全是精品| 乱系列少妇在线播放| 深爱激情五月婷婷| 岛国在线免费视频观看| 欧美一区二区亚洲| 久久精品国产亚洲网站| 亚洲欧美日韩高清在线视频| 亚洲在久久综合| 在线观看av片永久免费下载| 国产一区二区三区av在线 | 久久久久久大精品| 亚洲欧美日韩高清专用| 如何舔出高潮| 亚洲七黄色美女视频| 少妇猛男粗大的猛烈进出视频 | 午夜福利在线观看吧| av女优亚洲男人天堂| 国产乱人偷精品视频| 国产精品人妻久久久久久| 亚洲激情五月婷婷啪啪| 亚洲精品久久久久久婷婷小说 | 乱码一卡2卡4卡精品| 99久久九九国产精品国产免费| ponron亚洲| 高清日韩中文字幕在线| 精品久久国产蜜桃| 91久久精品国产一区二区成人| 精品久久久久久久末码| 午夜精品在线福利| 国产又黄又爽又无遮挡在线| 日韩在线高清观看一区二区三区| 97超碰精品成人国产| 人妻制服诱惑在线中文字幕| 久久99精品国语久久久| 精品久久国产蜜桃| 午夜老司机福利剧场| 日本撒尿小便嘘嘘汇集6| 精品日产1卡2卡| 91精品一卡2卡3卡4卡| 夫妻性生交免费视频一级片| 日韩精品有码人妻一区| 噜噜噜噜噜久久久久久91| 国产真实乱freesex| 精品久久久久久久久久免费视频| 亚洲国产日韩欧美精品在线观看| 亚洲欧美精品自产自拍| 99久久中文字幕三级久久日本| videossex国产| 久久人妻av系列| 亚洲人成网站在线播放欧美日韩| 成人国产麻豆网| 一级毛片久久久久久久久女| 国产真实乱freesex| 中文字幕人妻熟人妻熟丝袜美| 午夜精品国产一区二区电影 | 麻豆国产97在线/欧美| 一区二区三区免费毛片| 欧美潮喷喷水| 99视频精品全部免费 在线| 欧美精品一区二区大全| 日韩精品有码人妻一区| 久久精品综合一区二区三区| 国产黄片视频在线免费观看| 亚洲av熟女| 亚洲国产欧洲综合997久久,| 午夜老司机福利剧场| 九九爱精品视频在线观看| 国产精品1区2区在线观看.| 久久亚洲精品不卡| 亚洲高清免费不卡视频| 内地一区二区视频在线| 欧美高清成人免费视频www| 亚洲精品自拍成人| 国产亚洲91精品色在线| 一级毛片电影观看 | 三级国产精品欧美在线观看| 精品午夜福利在线看| 毛片一级片免费看久久久久| 热99re8久久精品国产| 亚洲精品乱码久久久v下载方式| 蜜桃亚洲精品一区二区三区| 久久人人爽人人片av| 亚洲在线自拍视频| 三级国产精品欧美在线观看| 亚洲欧美成人综合另类久久久 | 欧美一区二区精品小视频在线| 久久久久久伊人网av| 午夜老司机福利剧场| 亚洲不卡免费看| h日本视频在线播放| 国产三级在线视频| АⅤ资源中文在线天堂| 成人漫画全彩无遮挡| 少妇熟女欧美另类| 色5月婷婷丁香| 日韩欧美精品v在线| 国产成人福利小说| 26uuu在线亚洲综合色| 久久婷婷人人爽人人干人人爱| 岛国在线免费视频观看| 亚洲久久久久久中文字幕| 一个人看的www免费观看视频| 国产精品久久久久久亚洲av鲁大| 91aial.com中文字幕在线观看| 最近中文字幕高清免费大全6| 日韩欧美国产在线观看| 亚洲精品亚洲一区二区| 欧美zozozo另类| 精品国内亚洲2022精品成人| 久久精品国产自在天天线| 春色校园在线视频观看| 成人综合一区亚洲| 国产白丝娇喘喷水9色精品| 久久中文看片网| 美女高潮的动态| 成人高潮视频无遮挡免费网站| 亚洲18禁久久av| 国产精品伦人一区二区| 国产精品国产高清国产av| 国产精品三级大全| 亚洲精品乱码久久久v下载方式| 欧美一区二区精品小视频在线| 淫秽高清视频在线观看| 性欧美人与动物交配| 不卡视频在线观看欧美| 男女下面进入的视频免费午夜| 国产成人a∨麻豆精品| 成人二区视频| 国产精品野战在线观看| 人妻夜夜爽99麻豆av| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 舔av片在线| 成人午夜高清在线视频| 如何舔出高潮| 精品少妇黑人巨大在线播放 | 亚洲婷婷狠狠爱综合网| 欧美色欧美亚洲另类二区| 国产在线精品亚洲第一网站| 波野结衣二区三区在线| 97超视频在线观看视频| 1024手机看黄色片| 一级毛片aaaaaa免费看小| 搞女人的毛片| 国产三级在线视频| 国产一区二区在线av高清观看| 亚洲aⅴ乱码一区二区在线播放| 好男人在线观看高清免费视频| 亚洲在线自拍视频| 日韩一本色道免费dvd| 一个人免费在线观看电影| 亚洲av不卡在线观看| 一个人看视频在线观看www免费| 村上凉子中文字幕在线| av卡一久久| 蜜桃亚洲精品一区二区三区| 一级黄色大片毛片| 97超视频在线观看视频| 久99久视频精品免费| 色尼玛亚洲综合影院| 三级男女做爰猛烈吃奶摸视频| 欧美性猛交黑人性爽| 最近2019中文字幕mv第一页| 九九爱精品视频在线观看| 小说图片视频综合网站| 久久久久久伊人网av| 简卡轻食公司| 麻豆一二三区av精品| 国产精品福利在线免费观看| 亚洲av一区综合| 99久久人妻综合| 12—13女人毛片做爰片一| 亚洲美女视频黄频| 国产高清有码在线观看视频| 乱码一卡2卡4卡精品| 久久精品91蜜桃| 悠悠久久av| 99在线人妻在线中文字幕| 久久人人爽人人片av| 亚洲欧美日韩东京热| 亚洲欧洲国产日韩| 99久久无色码亚洲精品果冻| 啦啦啦观看免费观看视频高清| 三级男女做爰猛烈吃奶摸视频| 亚洲自偷自拍三级| 99热这里只有是精品在线观看| 丝袜喷水一区| 日韩人妻高清精品专区| а√天堂www在线а√下载| 在线播放国产精品三级| 白带黄色成豆腐渣| 精华霜和精华液先用哪个| 国产精华一区二区三区| 欧美一级a爱片免费观看看| 国内少妇人妻偷人精品xxx网站| 激情 狠狠 欧美| 久久久国产成人精品二区| 亚洲丝袜综合中文字幕| 黄片无遮挡物在线观看| 亚洲无线在线观看| 国产美女午夜福利| 国产大屁股一区二区在线视频| 精品不卡国产一区二区三区| 精品久久久久久久久亚洲| 亚洲精品国产av成人精品| 男人狂女人下面高潮的视频| 在线天堂最新版资源| 亚洲内射少妇av| 午夜激情欧美在线| 中文资源天堂在线| 亚洲精品自拍成人| 国产精品爽爽va在线观看网站| 国产一级毛片在线| 你懂的网址亚洲精品在线观看 | 麻豆乱淫一区二区| 十八禁国产超污无遮挡网站| 久久精品国产自在天天线| 99热这里只有是精品在线观看| 欧美三级亚洲精品| 爱豆传媒免费全集在线观看| 老师上课跳d突然被开到最大视频| 亚洲精品日韩av片在线观看| 欧美xxxx性猛交bbbb| 在线观看免费视频日本深夜| 一本久久精品| 国产熟女欧美一区二区| av在线天堂中文字幕| 中文字幕av在线有码专区| 日韩强制内射视频| 亚洲欧美日韩东京热| 性插视频无遮挡在线免费观看| 熟女电影av网| 亚洲人与动物交配视频| 亚洲精品日韩av片在线观看| 尤物成人国产欧美一区二区三区| 五月玫瑰六月丁香| 中文字幕av在线有码专区| 在线观看美女被高潮喷水网站| 黄色配什么色好看| 欧美xxxx黑人xx丫x性爽| 国产一级毛片七仙女欲春2| 日韩 亚洲 欧美在线| 亚洲欧美精品专区久久| 亚洲人成网站在线播| 国产私拍福利视频在线观看| 麻豆一二三区av精品| 我要搜黄色片| 小说图片视频综合网站| 熟女人妻精品中文字幕| 三级男女做爰猛烈吃奶摸视频| av卡一久久| 中文字幕av在线有码专区| 久久中文看片网| 一卡2卡三卡四卡精品乱码亚洲| 欧美+亚洲+日韩+国产| 亚洲国产色片| 成年女人看的毛片在线观看| 丝袜美腿在线中文| 搞女人的毛片| 一级二级三级毛片免费看| 亚洲av中文字字幕乱码综合| 久久久久久大精品| 中国美白少妇内射xxxbb| 一级毛片电影观看 | 超碰av人人做人人爽久久| 99久久人妻综合| 精品久久久久久久久亚洲| 欧美3d第一页| 欧美又色又爽又黄视频| 能在线免费观看的黄片| 人妻夜夜爽99麻豆av| 日本欧美国产在线视频| 九九在线视频观看精品| 岛国在线免费视频观看| 熟妇人妻久久中文字幕3abv| 青春草亚洲视频在线观看| 免费av观看视频| 91在线精品国自产拍蜜月| 亚洲色图av天堂| 女人十人毛片免费观看3o分钟| 欧美+日韩+精品| 一级毛片电影观看 | 国产片特级美女逼逼视频| 国语自产精品视频在线第100页| 我要搜黄色片| 欧美日韩国产亚洲二区| 狂野欧美激情性xxxx在线观看| 99久久精品一区二区三区| 国产国拍精品亚洲av在线观看| 成人美女网站在线观看视频| 男女做爰动态图高潮gif福利片| 国产成人精品一,二区 | 高清午夜精品一区二区三区 | 国产三级在线视频| 中文字幕人妻熟人妻熟丝袜美| 亚州av有码| 国产欧美日韩精品一区二区| 亚洲av熟女| 精品一区二区三区人妻视频| 色视频www国产| 久久久久久伊人网av| 欧美变态另类bdsm刘玥| 免费人成视频x8x8入口观看| 一进一出抽搐gif免费好疼| 国产精品久久久久久久久免| 丰满人妻一区二区三区视频av| 欧美最黄视频在线播放免费| a级毛片a级免费在线| 51国产日韩欧美| 男人舔女人下体高潮全视频| 国产探花在线观看一区二区| 最近中文字幕高清免费大全6| 国产白丝娇喘喷水9色精品| 亚洲欧美日韩高清专用| 亚洲成人中文字幕在线播放| 国产成人精品一,二区 | 国产极品精品免费视频能看的| 非洲黑人性xxxx精品又粗又长| 亚洲欧洲国产日韩| 麻豆一二三区av精品| 观看美女的网站| 亚洲国产精品合色在线| 免费人成视频x8x8入口观看| 国产高清不卡午夜福利| 一本一本综合久久| 亚洲激情五月婷婷啪啪| 久久草成人影院| 一区福利在线观看| 中文在线观看免费www的网站| 精品午夜福利在线看| 日韩视频在线欧美| 长腿黑丝高跟| 国产精品电影一区二区三区| 国国产精品蜜臀av免费| 亚洲av中文字字幕乱码综合| 国产精品爽爽va在线观看网站| 免费av观看视频| 久久综合国产亚洲精品| 不卡视频在线观看欧美| 国产精品嫩草影院av在线观看| 亚洲精品日韩在线中文字幕 | a级毛片免费高清观看在线播放| 岛国毛片在线播放| 国产成人aa在线观看| 热99在线观看视频| 18+在线观看网站| 日韩欧美精品免费久久| 亚洲欧美精品专区久久| 日本一本二区三区精品| 69人妻影院| 一进一出抽搐gif免费好疼| 99热这里只有精品一区| 久久人人爽人人爽人人片va| 少妇高潮的动态图| 美女 人体艺术 gogo| 免费av毛片视频| 亚洲无线观看免费| 在现免费观看毛片| 少妇丰满av| 亚洲人成网站在线观看播放| av视频在线观看入口| 22中文网久久字幕| 99久久成人亚洲精品观看| a级毛片a级免费在线| 国产亚洲精品久久久久久毛片| 免费在线观看成人毛片| 免费不卡的大黄色大毛片视频在线观看 | 91久久精品国产一区二区成人| 日韩一本色道免费dvd| 精品久久久久久久久久久久久| 18禁在线播放成人免费| 色哟哟哟哟哟哟| a级毛色黄片| 国产三级中文精品| 丝袜美腿在线中文| 青春草国产在线视频 | 国产精品av视频在线免费观看| 老女人水多毛片| 舔av片在线| 久久久久免费精品人妻一区二区| 久久午夜福利片| 一进一出抽搐gif免费好疼| 国产精品伦人一区二区| 亚洲成人中文字幕在线播放| 校园春色视频在线观看| 久久精品夜夜夜夜夜久久蜜豆| 久久人妻av系列| 久久久久久伊人网av| 久久久精品大字幕| 国产视频首页在线观看| 国内精品一区二区在线观看| 插阴视频在线观看视频| 在线观看av片永久免费下载| 久久久国产成人精品二区| 男女做爰动态图高潮gif福利片| 亚洲欧美日韩高清在线视频| 国产精品爽爽va在线观看网站| 亚洲性久久影院| 狠狠狠狠99中文字幕| 日本-黄色视频高清免费观看| 午夜精品国产一区二区电影 | 亚洲人与动物交配视频| 久久亚洲精品不卡| 国产午夜精品论理片| 99久国产av精品国产电影| 日本五十路高清| 国产精品久久久久久久久免| 国产男人的电影天堂91| 国产精品一区二区在线观看99 | 成人午夜高清在线视频| 男人舔奶头视频| 国产又黄又爽又无遮挡在线| 欧美激情在线99| 99国产极品粉嫩在线观看| 成人综合一区亚洲| 国产成人a∨麻豆精品| 91午夜精品亚洲一区二区三区| 男女那种视频在线观看| 99久久无色码亚洲精品果冻| 久久99蜜桃精品久久| 成人一区二区视频在线观看| 夜夜爽天天搞| 国产精品无大码| 少妇的逼水好多|