{'fraction_non_alphanumeric': 0.029024812463935373, 'fraction_numerical': 0.012290825158684363, 'mean_word_length': 5.0344707520891365, 'pattern_counts': {'":': 0, '<': 0, '<?xml version=': 0, '>': 1, 'https://': 0, 'lorem ipsum': 0, 'www.': 0, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 0, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
Natural Convection of Third Grade Non-Newtonian Fluid Flow in a Changeable Semi-permeable Medium Considering the Effects of Variable Heat Source and Magnetic Field

Natural Convection of Third Grade Non-Newtonian Fluid Flow in a Changeable Semi-permeable Medium Considering the Effects of Variable Heat Source and Magnetic Field. ABSTRACT The main purpose of this paper is to conduct momentum and energy analysis on non-Newtonian fluid flow in a semi-permeable medium between two infinite vertical flat plates considering the influences of medium porosity, heat source and magnetic field. In this study, porosity, heat source and magnetic field terms are considered to be variable. Initially, partial differential equations of momentum and energy are derived and turned into ordinary differential equations utilizing a similarity solution. Afterwards, a system of differential equations is solved by Least Square Method (LSM) and subsequently reliable functions are proposed for temperature and velocity distributions. To assess the precision of this method, the equations are also solved by a numerical method and an analytical approach known as Galerkin Method (GM). Regarding the comparisons, it can be implied that when the values of magnetic field, heat source and porosity are constant along the channel width, both LSM and GM methods are efficient and appropriate agreement is observed between the results. However, for the case in which these parameters linearly vary along the channel width, Galerkin approach shows more accurate results. Ultimately, the influences of Hartman number, porosity and heat source parameters on the velocity and temperature distributions are discussed.


INTRODUCTION
Natural convection is a mechanism of heat transport in which no external source influences the fluid motion. In fact, density difference due to the gravitational field is the major causality of this type of heat transport. Due to its potential applications in industry and engineering, this mechanism has attracted a great deal of attentions among researchers. In many cases, accurate anticipation of the behavior of fluid flows can remove harmful effects such as mechanical noises and vibrations from industrial systems (Jianu and Rosen, 2017). A very conventional industrial application of natural convection is free cooling of electronic devices. This heat transfer mechanism can exist in small-scale devices such as computer chips to large-scale physical systems such as geothermal systems, heat exchangers, chemical catalytic reactors, fibre and granular insulation, packed beds, and petroleum reservoirs. In the last couple of decades, natural convection heat transfer for Newtonian and non-Newtonian fluids streaming between two infinite parallel vertical plates has been studied by different researchers. In these systems, heat is transferred by natural convection mechanism from vertical plates to fluids. Analysis of this mechanism is significant especially when the moving fluid is minimally affected by forced convection (McCabe et al., 1993). Natural convection heat transfer between vertical flat plates was analysed for several different non-Newtonian fluids by Bruce and Na (Bruce and Na, 1967). In another investigation, they considered and analysed several natural convection problems in laminar regimes (Shenoy and Mashelkar, 1982). Rajagopal (Rajagopal and Na, 1985) conducted a thermodynamic analysis on different constitutive functions in the area of fluid mechanics. It should be pointed out that in fluid mechanics problems, to yield reliable results, inherent non-linear behavior of systems must be considered. However, solving nonlinear equations is often very difficult. To overcome the complexities, various methods such as instance Homotopy Perturbation Method (HPM) (Nayfeh, 1973;Ganji et al., 2010;Rafei and Ganji, 2006;He, 2006;Soltanian et al., 2010;Esmaeilpour and Ganji, 2007;Ganji et al., 2007;Ganji et al., 2006), Variational Iteration Method (VIM) (Wang and He, 2007), Homotopy Analysis Method (HAM) (Hassan and Rashidi, 2013;Shahzad and Ramzan, 2012), Adomian Decomposition Method (ADM) (Abbasbandy, 2007), Least Square Method (Shi and Ren, 2010), Galerkin Method (Hessari, 2014) and Collocation Method (CM) (Petroudi et al., 2014) have been introduced by different researchers. It is worth mentioning that the aforementioned methods offer accurate successive approximations for solution of ODEs. Etbaeitabari and Domairy (Etbaeitabari et al., 2013) employed VIM to evaluate the natural convection of a non-Newtonian fluid streaming between two infinite vertical plates (Etbaeitabari et al., 2013). Ziabakhsh and Domairy (Ziabakhsh and Domairry, 2009) analytically investigated the natural convection of a non-Newtonian fluid between two infinite parallel vertical flat plates considering the nature of the fluid. Hatami et al. (Hatami et al., 2014) numerically and analytically performed a heat transfer and flow analysis on a non-Newtonian third grade nanofluid flow in porous medium considering magnetic field.  analytically solve the problem of nanofluid hydrothermal behavior in presence of magnetic field and Brownian motion of nanoparticles using Differential Transformation Method.  used differential transformation method to develop a theoretical solution for the case of laminar steady magnetohydrodymanics nanofluid streaming in a semi-permeable channel. Azimi et al. (Azimi et al., 2014) analytically studied the heat transfer of an unsteady nanofluid flow between two moving parallel plates applying Galerkin Optimal Homotopy Asymptotic Method. In another paper, Azimi et al. (Azimi and Riazi, 2016) proposed an analytical solution for the problem of a steady nanofluid streaming through a convergent-divergent channel utilizing Reconstruction of Variational Iteration method. Shahriari et al. (Shahriari et al., 2018) analytically investigated the effect of viscous dissipation on temperature distribution of a two-dimensional unsteady nanofluid streaming between two moving parallel plates utilizing Akbari-Ganji Method. The main purpose of the present paper, for the first time, is to carry out a momentum and energy analysis on non-Newtonian fluid flow between two infinite vertical flat plates in a changeable semi-permeable medium considering variable heat source and magnetic field terms. Initially, partial differential equations of momentum and energy are derived and turned into ordinary differential equations using a similarity solution. Afterwards, a system of differential equations is solved by LSM and subsequently reliable functions are presented for temperature and velocity distributions. In order to evaluate the reliability of this method, the governing equations are also solved and compared by a numerical method and an analytical approach known as GM. Subsequently, the influences of Hartman number, porosity and heat source parameters on the velocity and temperature distributions are presented. Finally, variations of velocity with change types of porosity and magnetic field are studied. Figure 1 illustrates the schematic of the problem analysed in the current paper. As can be seen in this figure, the system includes two vertical flat plates (the distance between the plates is considered to be 2b) and a non-Newtonian fluid streaming between the plates. In this investigation, a steady, incompressible, non-Newtonian fluid flow between two parallel vertical plates in a semi-permeable medium with variable heat source and magnetic field terms is considered. The schematic of the problem is shown in Figure 1. The temperatures of the walls located at = + and = − are constant and assumed to be 2 and 1 ( 1 > 2 ), respectively. This temperature difference results in upward motion of fluid near the wall at = − and downward motion of fluid near the wall at = + .

GOVERNING EQUATIONS
Generally, a viscous fluid is governed by continuity and Navier-Stokes equations. Since the fluid is considered to be incompressible, the conservation equations of momentum and energy are stated as follows (Rajagopal and Na, 1985;Hatami et al., 2014;: where and 0 are porosity of medium and heat source parameter, respectively. Stress in a third-grade non-Newtonian fluid ( * ) is presented in the following (Hatami et al., 2014): where 3 , 1 and 2 are temperature-dependent material modules which. In Eq. 3, − represents the spherical stress related to the restraint of incompressibility and is velocity gradient. The definitions of the kinematical tensors, 1 and 2 , can be found in (Hatami et al., 2014). Variable magnetic field, porosity of medium and heat source are defined as follows: By using the following similarity variables, the governing momentum and energy equations can be defined in simpler forms (Rajagopal and Na, 1985): By substituting the above parameters into the Navier-Stokes and energy equations, following equation are obtained: where , , , , and are Prandtl number, Eckert number, dimensionless non-Newtonian viscosity, dimensionless porosity parameter, Hartman number and dimensionless heat source parameter, respectively. These parameters are defined in the following: In this study, for solving the coupled and non-linear Eqs. (6) and (7), least square and Galerkin methods are used.

THEORETICAL METHODS
Least square and Galerkin methods are analytical methods operating based on Weighted Residuals Methods. In these methods, a differential operator represented by is applied on the function of to generate function .
where function is approximated by � and is a linear combination of a set of functions.
Since cannot be calculated by applying operator on the approximated function , an error will be appeared in the calculations which can be defined as follows: Minimizing the average residual of Eq. 13 is the main purpose of these three procedures.
In LSM, can be calculated by the following correlation: In GM, is defined as follows: It is notable that the number of weight functions are equal to the number of unknown parameters of the approximated function which must satisfy the considered initial and boundary conditions.

Comparison of Analytical and Numerical Results
As previously expressed, LSM is used for solving the non-linear governing equations. Analysis and solution of the governing equations are implemented in Matlab. By applying the boundary conditions, following approximated functions are considered for velocity and temperature distributions, respectively. It is worth mentioning that these functions are calculated by using try-and-error technique.
By using the velocity and temperature approximated functions and applying the LSM, unknown coefficients in Eqs. 17 and 18 will be obtained. For validation purpose, the values of non-dimensionalized numbers including , , , , and are taken to be unity and constant. By doing some algebraic computations, following equations are achieved for velocity and temperature distributions: Results calculated by the current numerical and analytical solutions for flow velocity and temperature are indicated in Figures 2 and 3, respectively. As can be observed in these figures, there is an appropriate agreement between the numerical and analytical results. Hence, it can be concluded that for this problem LSM is an efficient technique for solving the momentum and energy equations.
Tables 1 and 2 list the values of velocity and temperature calculated by using the numerical and analytical approaches which includes LSM and GM for two different cases. As presented in Table 1, for the first case, δ = Pr = Ec = M = P = Ha = 1. In this case, porosity and magnetic field parameters are constant along the channel width. The average error related to LSM for non-dimensionalized functions of velocity and temperature are found to be 0.034% and 0.16% in this case, respectively. The average error related to GM for non-dimensionalized functions of velocity and temperature are found to be 0.035% and 0.17%, respectively. By comparing the values of errors, it is implied that precision and efficiency of both of the methods are considerable high. As listed in Table 2, calculations are conducted for the case in which δ = Pr = Ec = M = P = Ha = 1. For this case, porosity and magnetic field parameters are proportional to 1 − and linearly vary along the channel width. In this case the average error related to LSM for velocity and temperature functions are found to be 0.5% and 3.36%, and for GM the values of average error are found to be 0.088% and 0.42%, respectively. Results show that when porosity, heat source and magnetic field parameters are constant along the channel width, both of the method are efficient. Nevertheless, for the case in which the parameters vary linearly, Galerkin method possesses better performance.

Effects of Porosity on Velocity and Temperature Distributions
Figures 4 and 5 describe the influences of dimensionless porosity parameter on the velocity and temperature when the values of the other non-dimensionalized parameters including δ, Pr, Ec, Ha and M are unity. According to these figures, it can be seen that porosity parameter which appears in the form of (1 − ) linearly varies. With regard to Figure 4, it can be observed that with an increase in the value of porosity parameter, velocity reduces. In other words, rising this parameter results in an increment in porosity and a reduction in permeability leading to a greater resistance against the fluid flow and subsequently a decline in the value of velocity. The effect of channel width on the temperature distribution for several values of porosity parameter is delineated in Figure  5. As can be seen in Figure 5, variation of temperature with channel width for different values of porosity parameters are very close to each other. Therefore, it can be concluded that porosity has negligible impact on the flow temperature.  (1 − ), and ( + 1) in the momentum equation, respectively. As can be seen in Figure 6, when porosity parameter is constant, velocity distribution is symmetric along the channel width. In descending order of porosity, velocity of the flow is zero before reaching the centerline of the channel and also on the left-hand side of the centerline. In this case, the value of resistance force due to the flow impermeability on the left-hand side of the centerline is greater on the left-hand side of the centerline. In this regard, velocity of the flow will be zero in this position. In ascending order of porosity, velocity of the flow is zero on the right-hand side of the centerline.

Effects of Heat Source Term on Velocity and Temperature Distributions
In Figures 7 and 8, the influences of non-dimensionalized heat source parameter on velocity and temperature distributions when the other parameters including δ, Ha, P, Pr and Ec are unity are indicated. Regarding Figure 7, with an increment in the value of heat source parameter, velocity of flow reduces along the channel width. The reason is that any change in temperature due to existence of heat source will lead to remarkable variations in buoyancy force and subsequently the flow velocity. With respect to Figure 8, it can be seen that an enhancement in the value of heat source parameter results in a reduction in temperature along the channel width. According to Figure 8, it must be pointed out that by adding heat source term in energy equation, temperature varies nonlinearly. With increasing the value of heat source parameter, curvature of temperature distribution increases. In these figures, all the considered dimensionless parameters except the Hartman number are assumed to be unity. In this case, magnetic field parameter is appeared in the form of 2 (1 − ) in momentum equation. As shown in Figure 9, with an increment in the value of Hartman number, velocity in each point of the channel width declines. In fact, increasing the Hartman number rises the Lorentz force which is related to the magnetic field leading to an increment in resistance against the flow and subsequently reduction in the value of velocity. In Figure 10, the effect of Hartman number on temperature distribution is shown. According to this figure, Hartman number has insignificance impact on temperature. In other words, the effect of magnetic field on temperature is negligible. Figure 11 illustrates the variation of velocity with position for several change types of It should be noted that in this figure all the considered non-dimensionalized parameters are taken to be unity. In descending order, constant value and ascending order of magnetic field, this parameter is appeared in forms of 2 (1 − ), 2 and 2 ( + 1), respectively, in the momentum equation. As can be seen in Figure 11, when magnetic field parameter is constant, velocity distribution is symmetric along the channel width. In descending order of magnetic field, velocity of the flow is zero before reaching the centerline of the channel as well as the left-hand side of the centerline. In this case, the value of resistance force due to the magnetic field on the left-hand side of the centerline is further. In this regard, velocity of the flow will be zero in this position. In ascending order of magnetic field, velocity of the flow is zero on the right-hand side of the centerline because of lower resistance in this region.