ISSN 975-333 Mapana J Sci,, 3(22), 25-232 https://doi.org/.2725/mjs.22.4 MHD and Thermal Dispersion-Radiation Effects on Non-Newtonian Fluid Saturated Non-Darcy Mixed Convective Flow with Melting Effect B D C N Prasad, * K Hemalatha, and J Siva Ram Prasad Abstract The problem of mixed convection along a vertical plate in a non-newtonian fluid saturated non-darcy porous medium in the presence of melting, thermal dispersionradiation effects for aiding and opposing external flows is analyzed. Similarity solution for the governing equations is obtained for the flow equations in steady state. The equations are numerically solved by using Runge-kutta fourth order method coupled with shooting technique. The effects of melting (M), thermal dispersion (D)- radiation (R), inertia (F) and mixed convection (Ra/Pe) on velocity and temperature distributions are examined in the presence of magnetic parameter MH. Keywords: Porous medium, Non-Newtonian Fluid, Melting, Thermal Dispersion, Radiation. Introduction Convection heat transfer in porous media in the presence of melting effect has received some attention in recent years. This stems from the fact that this topic has significant direct application * Department of Computer Applications, P.V.P. Institute of Technology, Vijayawada; bdcnprasad@gmail.com Department of Mathematics, V.R. Siddhartha Engineering College, Vijayawada; hemakphd@yahoo.co.in Department of Mathematics, V.R. Siddhartha Engineering College, Vijayawada; nanijsrp@yahoo.co.in Received: July 22, Reviewed: Aug. 22 25
B D C N Prasad, K Hemalatha, and J Siva Ram Prasad ISSN 975-333 in permafrost melting, frozen ground thawing, casting and bending processes as well as phase change metal. The study of melting effect is considered by many researchers in Newtonian fluids. Non- Newtonian power law fluids are so wide spread in industrial process and in the environment. The melting phenomena on free convection from a vertical front in a non- Newtonian fluid saturated porous matrix are analyzed by Poulikakos and Spatz. [] Nakayama and Koyama [2] studied the more general case of free convection over a non-isothermal body of arbitrary shape embedded in a porous medium. Rastogi and Poulikakos [3] examined the problem of double diffusive convection from a vertical plate in a porous medium saturated with a non-newtonian power law fluid. Shenoy [4] presented many interesting applications of non-newtonian power law fluids with yield stress on convective heat transport in fluid saturated porous media. Considering geothermal and oil reservoir engine-ering applications, Nakayama and Shenoy [5] studied a unified similarity transformation for Darcy and non-darcy forced, free and mixed convection heat transfer in non-newtonian inelastic fluid saturated porous media. Later Shenoy [6] studied non-darcy natural, forced and mixed convection heat transfer in non-newtonian power law fluid saturated porous media. Effect of melting and thermo-diffusion on natural convection heat mass transfer in a non-newtonian fluid saturated non-darcy porous medium was studied by R. R. Kairi and P. V. S. N. Murthy. [7] It is noted that the velocity, temperature and concentration profiles as well as the heat and mass transfer coefficients are significantly affected by the melting phenomena and thermaldiffusion in the medium. The non-linear behavior of non- Newtonian fluids in a porous matrix is quite different from that of Newtonian fluids in porous media. The prediction of heat or mass transfer characteristics for mixed or natural convection of non- Newtonian fluids in porous media is very important due to its practical engineering applications such as oil recovery and food processing. If the temperature of the surrounding fluid is rather high, radiation effects play an important role and this situation does exist in space technology. In such cases, one has to take into account the effect of thermal radiation and mass diffusion. On the other hand it is worth mentioning that heat transfer simultaneous 26
MHD and Thermal Dispersion-Radiation Effects Mapana J Sci,, 3(22) radiation and convection is very important in the context of space technology and processes involving high temperatures. Recently, A.Y. Bakier et al. [8] studied Group method analysis of melting effect on MHD mixed convection flow from a radiative vertical plate embedded in saturated porous medium for Newtonian fluids. He developed linear transformation group approach to simulate problem of hydro magnetic heat transfer by mixed convection along vertical plate in a liquid saturated porous medium in the presence of melting and thermal radiation effects for opposing external flow. He studied the effects of the pertinent parameters on the rate of the heat transfer in terms of the local Nusselt number at the solid-liquid interface. More recently Melting and radiation effects on mixed convection from a vertical surface embedded in a non-newtonian fluid saturated non-darcy porous medium for aiding and opposing external flows is analyzed by A. J. Chamka et al. [9] They obtained representative flow and heat transfer results for various combinations of physical parameters. The present paper is aimed at analyzing the effect of melting and thermal dispersion-radiation on steady mixed convective heat transfer from a vertical plate embedded in a non-newtonian power law fluid saturated non-darcy porous medium for aiding and opposing external flows. 2. Mathematical formulation A mixed convective heat transfer in a non-darcy porous medium saturated with a homogeneous non-newtonian fluid adjacent to a vertical plate, with a uniform wall temperature is considered. This plate constitutes the interface between the liquid phase and the solid phase during melting inside the porous matrix Figure. Schematic diagram of the problem 27
B D C N Prasad, K Hemalatha, and J Siva Ram Prasad ISSN 975-333 at steady state. The plate is at a constant temperature T m at which the material of the porous matrix melts. Fig. shows the co-ordinate and the flow model. The x-coordinate is measured along the plate and the y-coordinate normal to it. The solid phase is at temperature T < T m. A thin boundary layer exists close to the right of vertical plate and temperature changes smoothly through this layer from T m to T (T m < T ) which is the temperature of the fluid phase. Taking into account the effect of thermal dispersion the governing equations for steady non-darcy flow in a non-newtonian fluid saturated porous medium can be written as follows. The continuity equation is The momentum equation is + = () The energy equation is + + = u + v = α (2) (3) In the above equations, the term which represents the buoyancy forced effect on the flow field has ± signs. The plus sign indicates aiding buoyancy flow where as the negative sign stands for buoyancy opposed flow. Here u and v are the velocities along x and y directions respectively, n is the power-law fluid viscosity index, T is Temperature in the thermal boundary layer, K is Permeability, C is Forchheimer empirical constant, β is coefficient of thermal expansion, ν is Kinematics viscosity, ρ is Density, Cp is Specific heat at constant pressure, g is acceleration due to gravity, σ is electrical conductivity, B is magnetic field intensity and thermal diffusivity α=α m+α d, where α m is the molecular diffusivity and α d is the dispersion thermal diffusivity due to mechanical dispersion. As in the linear model proposed by Plumb[], the dispersion thermal diffusivity α d is proportional to the velocity component i.e α d = γud, where γ is the dispersion coefficient and d is the mean particle diameter. The radiative heat flux term q is written using the Rosseland approximation (Sparrow and Cess[], Raptis[2]) as 28
MHD and Thermal Dispersion-Radiation Effects Mapana J Sci,, 3(22) q = σ Where σ R is the Stefan - Boltzmann constant and a is the mean absorption coefficient. The physical boundary conditions for the present problem are (4) y=, T=T m, k = ρ[h + C (T T )]v (5) and y, T T, u = u (6) Where h sf and C s are latent heat of the solid and specific heat of the solid phases respectively and u is the assisting external flow velocity, k=αρc p is the effective thermal conductivity of the porous medium. The boundary condition (5) means that the temperature on the plate is constant and thermal flux of heat conduction to the melting surface is equal to the sum of the heat of melting and the heat required for raising the temperature of solid to its melting temperature T m. Introducing the stream function ψ with u =, and v = The continuity equation () will be satisfied and the equations (2) and (3) transform to ψ ψ + ψ + ν ψ = (α + γ ψ = β ν d) + σ ρ Introducing the similarity variables as ψ =f(η)(α m u x) /2, η = (7) (8), θ(η) =, the momentum equation (7) and energy equation (8) are reduced to and n(+mh) f f + 2Ff f = θ (9) ( + Df )θ + f + Df θ + R(θ + C ) θ + 3θ (θ + C ) = () where the prime symbol denotes the differentiation with respect to the similarity variable η and MH = is Magneto Hydro 29
B D C N Prasad, K Hemalatha, and J Siva Ram Prasad ISSN 975-333 Dynamic parameter, Ra x/pe x is the mixed convection parameter, Ra = ( ) is the local Rayleigh number, Pe = is the local peclet number. F = f (Pe ), f = is the nondarcian parameter. Pe d is the pore diameter dependent peclet number. D = is the dispersion parameter. C = is the temperature ratio, R = ( ) is the radiation parameter. Taking into consideration, the thermal dispersion effect together with melting, the boundary conditions (5) and (6) take the form η=,θ=,f()+{+df ()}2Mθ ()= () and η, θ=, f = (2) where M = ( ) is the melting parameter. ( ) 3. Solution Procedure The dimensionless equations (9) and () together with the boundary conditions () and (2) are solved numerically by means of the fourth order Runge-Kutta method coupled with double shooting technique. The solution thus obtained is matched with the given values of f ( ) and θ(). In addition the boundary condition η is approximated by η max = 8 which is found sufficiently large for the velocity and temperature to approach the relevant free stream properties. Numerical computations are carried out for F=,, ; D=,, ; Ra/Pe=, 2,5; M=,, 2; n =, 2.5, R=,, ; MH=,, 5, ; Cr =.. 4. Results and Discussion In order to get clear insight on the physics of the problem, a parametric study is performed and the obtained numerical results are displayed with the help of graphical illustrations. The parameters governing the physics of the present study are the melting (M), the mixed convection (Ra/Pe), the inertia (F), thermal dispersion (D), and Fluid viscosity index (n), radiation(r), 22
MHD and Thermal Dispersion-Radiation Effects Mapana J Sci,, 3(22) temperature ratio(cr). The numerical computations were carried out for the fixed value of buoyancy parameter Ra/Pe = for both the aiding and opposing external flows. The results of the parametric study are shown in figures 2-22. Effect of Melting Parameter on Velocity Profiles.3.2. Aiding External flow, Ra/Pe= M= M= M=2 M= M= M=2 f ' Opposing External flow, Ra/Pe= F=,n=,MH=,D=,R=,Cr=. 2 3 4 5 6 7 8 Figure 2 Figure 2 shows the effect of melting parameter M on velocity profiles for both aiding and opposing external flows. It is observed that different behavior exists between the velocity profiles in the presence of solid phase melting effect in the case of aiding and opposing external flows. Effect of Melting Parameter on Temperature Profile M= M= M=2.4.3 Aiding External flow, Ra/Pe=.2. F=,n=,MH=,D=,R=,Cr=. 2 3 4 5 6 7 8 Figure 3 22
B D C N Prasad, K Hemalatha, and J Siva Ram Prasad ISSN 975-333 The increase in the melting parameter M causes the increase of the velocity profiles for the case of aiding flow. Whereas the opposite behavior for the velocity profiles as M increases is found in the case of opposing flow case. Effect of Melting Parameter on Temperature Profile M= M= M=2.4.3 Opposing External flow, Ra/Pe=.2. F=,n=,MH=,D=,R=,Cr=. 2 3 4 5 6 7 8 Figure 4 Figure 3 and Figure 4 show the effect of the melting parameter M on temperature distributions for aiding and opposing external flows respectively. It is observed that as increasing the value of the melting parameter M, the temperature distributions decrease for both cases of aiding and opposing external flow conditions. Effect of Thermal Dispersion Perameter on Velocity Profiles.3.2. Aiding External flow, Ra/Pe= D= D= D= D= D= D= f ' Opposing External flow, Ra/Pe= F=,n=,M=2,MH=,R=,Cr=. 2 3 4 5 6 7 8 Figure 5 222
MHD and Thermal Dispersion-Radiation Effects Mapana J Sci,, 3(22) Effect of Thermal Dispersion Parameter on Temperature Profile D= D= D=.4.3 Aiding External flow, Ra/Pe=.2. F=,n=,M=2,MH=,R=,Cr=. 2 3 4 5 6 7 8 Figure 6 The effects of thermal dispersion parameter D on the velocity profiles for both cases of aiding and opposing external flow conditions are plotted in Figure 5. It is noted that for the case of aiding flow the velocity profiles increase with increase in the value of D. But this effect is found opposite in the case of opposing flow. Effect of Thermal Dispersion Parameter on Temperature Profile.4.3 D= D= D= Opposing External flow, Ra/Pe=.2. F=,n=,M=2,MH=,R=,Cr=. 2 3 4 5 6 7 8 Figure 7 The effects of thermal dispersion parameter D on temperature distributions for aiding and external flow conditions are plotted in Figure 6 and Figure 7 respectively. It is noted that increasing the values of D leads to decrease in the liquid temperature distributions in both cases. 223
B D C N Prasad, K Hemalatha, and J Siva Ram Prasad ISSN 975-333 Effect of Power Law Fluid Viscosity Index on Velocity Profile.4 F=,M=2, MH=,D=,R=,Cr=. n=.3 n=2.5 n= n=2.5.2 Aiding External flow, Ra/Pe=. f ' Opposing External flow, Ra/Pe= 2 3 4 5 6 7 8 Figure 8 Figure 8 shows the effects of power law fluid viscosity index n on velocity profiles for both aiding and opposing flow cases. It is observed that in aiding flow case as n value increases, velocity profiles decreases. Whereas opposite behavior obtained in velocity profiles in opposing flow case. Effect of Power Law Fluid Viscosity Index on Temperature Profile n= n=2.5.4.3 Aiding External flow, Ra/Pe=.2. F=,M=2,MH=,D=,R=,Cr=. 2 3 4 5 6 7 8 Figure 9 Figure 9 and Figure show the effects of power law fluid viscosity index n on temperature distributions for both aiding and opposing flow cases respectively. Significant effect is not found. In aiding flow case as n increases, the temperature distributions decrease. Whereas opposite result is found in opposing flow case. 224
MHD and Thermal Dispersion-Radiation Effects Mapana J Sci,, 3(22) Effect of Power Law Fluid Viscosity Index on Temperature Profile n= n=2.5.4.3 Opposing External flow, Ra/Pe=.2. F=,M=2,MH=,D=,R=,Cr=. 2 3 4 5 6 7 8 Figure Effect of Radiation Parameter on Velocity Profiles.3.2. n=,f=,m=2,mh=,d=,cr=. Aiding External Flow, Ra/Pe= R= R= R= R= R= data 5 f ' Opposing External Flow, Ra/Pe= 2 3 4 5 6 7 8 Figure Figure shows the effect of radiation parameter R on velocity profiles for both aiding and opposing external flows. The increase in the radiation parameter R causes the increase of the velocity profiles for aiding flow case. Whereas in the case of opposing flow the result is found opposite 225
B D C N Prasad, K Hemalatha, and J Siva Ram Prasad ISSN 975-333 Effect of Radiation Parameter on Temperature Profile R= R= R=.4.3 Aiding External Flow, Ra/Pe=.2. n=,f=,m=2,mh=,d=,cr=. 2 3 4 5 6 7 8 Figure 2.4.3.2 Effect of Radiation Parameter on Temperature Profile Opposing External Flow, Ra/Pe= R= R= R= n=,f=,m=2,mh=,d=,cr=.. 2 3 4 5 6 7 8 Figure 3 Figure 2 and Figure 3 show the effect of radiation parameter R on temperature distributions for aiding and opposing external flow conditions respectively. It is observed that same phenomena exist in both cases. As R increases the temperature distributions decrease for both cases. 226
MHD and Thermal Dispersion-Radiation Effects Mapana J Sci,, 3(22) Effect of MHD parameter on velocity profiles.3.2. Aiding External Flow, Ra/Pe= MH= MH=5 MH= MH= MH=5 MH= f ' Opposing External Flow, Ra/Pe= F=,M=2,n=,D=,R=,Cr=. 2 3 4 5 6 7 8 Dimensionless distance, Figure 4 Figure 4 shows the effect of Magneto Hydro Dynamic parameter MH on velocity profiles for both aiding and opposing flow cases. It is observed that in aiding flow case as MH value increases, velocity profiles decreases. Whereas opposite behavior obtained in velocity profiles in opposing flow case..4.3.2 Effect of MHD parameter on Temperature profile Aiding External Flow, Ra/Pe= F=,M=2,n=,D=,R=,Cr=. MH= MH=5 MH=. 2 3 4 5 6 7 8 Figure 5 227
B D C N Prasad, K Hemalatha, and J Siva Ram Prasad ISSN 975-333.4 Effect of MHD parameter on temperature profile Opposing External Flow, Ra/Pe= MH= MH=5 MH=.3.2. F=,M=2,n=,D=,R=,Cr=. 2 3 4 5 6 7 8 Figure 6 Figure 5 and Figure 6 show the effect of Magneto Hydro Dynamic parameter MH on temperature distributions for aiding and opposing external flow conditions respectively. In aiding flow case as MH increases, the temperature distributions decrease. Whereas opposite result is found in opposing flow case. Effect of Non-Darcy Porous Medium Parameter on Velocity Profiles 2.5 F= n=,m=2,mh=,d=,r=,cr=. F= F= 2 F= Aiding External flow, Ra/Pe= F= F=.5 f ' Opposing External flow, Ra/Pe= 2 3 4 5 6 7 8 Figure 7 Figure 7 shows the effect of Non-Darcy Porous Medium parameter F on velocity profiles for both aiding and opposing flow cases. It is observed that in aiding flow case as F value increases, velocity profiles decreases. Whereas opposite behavior obtained in velocity profiles in opposing flow case. 228
MHD and Thermal Dispersion-Radiation Effects Mapana J Sci,, 3(22) Figure 8 and Figure 9 show the effect of Non-Darcy Porous Medium parameter F on temperature distributions for aiding and opposing external flow conditions respectively. In aiding flow case as F increases, the temperature distributions decrease. Whereas opposite result is found in opposing flow case. Effect of Non-Darcy Porous Medium Parameter on Temperature Profile F= F= F=.4 Aiding External flow, Ra/Pe=.3.2. n=,m=2,mh=,d=,r=,cr=. 2 3 4 5 6 7 8 Figure 8 Effect of Non-Darcy Porous Medium Parameter on Temperature Profile F= F= F=.4.3 Opposing External flow, Ra/Pe=.2. n=,m=2,mh=,d=,r=,cr=. 2 3 4 5 6 7 8 Figure 9 229
B D C N Prasad, K Hemalatha, and J Siva Ram Prasad ISSN 975-333 Effect of Mxed Convection Parameter on Velocity Profiles.6.4.2 Aiding External Flow Ra/Pe= Ra/Pe=2 Ra/Pe=5 Ra/Pe= Ra/Pe=2 Ra/Pe=5 f ' Opposing External Flow.4.2 F=,n=,M=2,MH=,R=,D=,Cr=. 2 3 4 5 6 7 8 Figure 2 Figure 2 shows the effect of Mixed Convection parameter Ra/Pe on velocity profiles for both aiding and opposing flow cases. It is observed that in aiding flow case as Ra/Pe value increases, velocity profiles increases. Whereas opposite behavior obtained in velocity profiles in opposing flow case. Effect of Mixed Convection parameter on Temperature Profile Ra/Pe= Ra/Pe=2 Ra/Pe=5.4.3 Aiding External Flow.2. F=,n=,M=2,MH=,R=,D=,Cr=. 2 3 4 5 6 7 8 Figure 2 23
MHD and Thermal Dispersion-Radiation Effects Mapana J Sci,, 3(22) Effect of Mixed Convection Parameter on Temperature Profile Ra/Pe= Ra/Pe=2 Ra/Pe=5.4.3 Opposing External Flow.2 Figure 22 Figure 2 and Figure 22 show the effect of Mixed Convection parameter Ra/Pe on temperature distributions for aiding and opposing external flow conditions respectively. In aiding flow case as Ra/Pe increases, the temperature distributions increase. Whereas opposite result is found in opposing flow case. References. F=,n=,M=2,MH=,R=,D=,Cr=. 2 3 4 5 6 7 8 [] D Poulikakos and T L Spatz, Non-Newtonian natural convection at a melting front in a permeable solid matrix, Int. Commun. Heat and Mass Transfer, vol. 5, pp. 593-63, 988. [2] A Nakayama and H Koyama, Buoyancy induced flow of non- Newtonian fluids over a non isothermal body of arbitrary shape in a fluidsaturated porous medium, Appl. Scientific Research, vol. 48, pp. 55-7, 99. [3] S K Rostagi and D Poulikakos, Double-diffusion from a vertical surface in a porous region Saturated with a non-newtonian saturated with a non- Newtonian fluid, Int. J. Heat Mass Transfer, vol. 38, pp. 935-946, 995. [4] A V Shenoy, Non-Newtonian fluid heat transfer in a porous media, Advances in Heat transfer, vol. 24, pp. -9, 994. [5] A Nakayama and A V Shenoy, A unified similarity transformation for Darcy and non-darcy forced, free and mixed convection heat transfer in non-newtonian inelastic fluid saturated porous media, The Chemical Eng. J., vol. 5, pp. 33-45, 992. [6] A V Shenoy, Darcy-Forchheimer natural forced and mixed convection heat transfer in Non-Newtonian power law fluid saturated porous media, Transport in Porous Media, vol., pp. 29-24, 993. 23
B D C N Prasad, K Hemalatha, and J Siva Ram Prasad ISSN 975-333 [7] R R Kairi and P V S N Murty, Effect of melting and thermo-diffusion on natural convection heat mass transfer in a non-newtonian fluid saturated non-darcy porous medium, The Open Transport Phenomena J., vol., pp. 7-4, 29. [8] A Y Bakier et al., Group method analysis of melting effect on MHD mixed convection flow from radiate vertical plate embedded in a saturated porous media, Commun. Nonlinear Sci. and Numerical Simulation, vol. 4, pp. 26-27, 29. [9] A J Chamkha et al., Melting and radiation effects on mixed convection from a vertical surface embedded in a non-newtonian fluid saturated non-darcy porous medium for aiding and opposing external flows, Int. J. Phy. Sci., vol. 5, pp.22-224, 2. [] O A Plumb, The effect of thermal dispersion on heat transfer in packed bed boundary layers, Proc. st ASME/JSME Thermal Eng. Joint Conference, 2, pp. 7-2, 983. [] E M Sparrow and R D Cess, Radiation Heat Transfer, Washington: Hemisphere, 978. [2] A Rapits, Radiation and free convection flow through a porous medium, Int. Commun. Heat Mass Transfer, vol. 25, pp.289-295, 998. [3] W T Cheng and C H Lin, Melting effect on mixed convection heat transfer with aiding and opposing flows from the vertical plate in a liquid saturated porous medium, Int. J. Heat Mass Transfer, vol. 5, pp. 326-334, 27. 232