Joural of Applied Fluid Mechaics, Vol. 9, No. 6, pp. 949-9, 6. Available olie at www.jafmolie.et, ISSN 7-7, EISSN 7-64. Variatioal Approach of Maragoi Mied Covectio Boudary Layer Flow with Pade Techique M. Gubes ad G. Oturac Departmet of Mathematics, Karamaoglu Mehmetbey Uiversity, Karama, 7, Turkey Departmet of Mathematics, Selcuk Uiversity, Koya, 4, Turkey Correspodig Author Email: mgubes@kmu.edu.tr (Received September 9, ; accepted February, 6) ABSTRACT I this paper, Variatioal Iteratio method with combiig Pade approimatio (Modified Variatioal Iteratio Method-MVIM) is performed to Maragoi covectio flow over the surface with buoyacy effects which is occurred gravity ad eteral pressure. After the appropriate trasformatio of equatios, we get the dimesioless form to solve umerically with modified variatioal iteratio method. We compare the our results with well-kow asymptotic epasio method used by Zhag Ya ad Zheg Liacu ad also compare with Fourth order Ruge Kutta solutio which are preseted i tables. Very efficiet ad accurate results are obtaied with preseted method. Keywords: Variatioal iteratio method; Maragoi covectio; Pade approimatio; Boudary layer flow. NOMENCLATURE f velocity fuctio g gravitatioal acceleratio Pr Pradtl umber S iterface T temperature T positive Temp. icremet alog the iterface compoet of velocity ue ( ) velocity outside boudary layer y compoet of velocity y Cartesia coordiates ormal to S Cartesia coordiate alog Marao mied covectio parameter similarity variable buoyacy forces parameter surface tesio temperature fuctio Lagrage multiplier stream fuctio. INTRODUCTION Maragoi covectio flow is stimulated by variatios of surface tesio throughout liquidliquid or liquid-gas surfaces ad it is importat i may fields of ature ad egieerig. Fudametal treatmet of Maragoi flow has bee aalyzed by (Gelles 978; Napolitao 98, 979, 978; Okao et al. 989). Some of umerical studies deped o the Maragoi covectio i various geometries have bee preseted by (Golia ad Viviai 98, 986), (Cristopher ad Wag ), (Pop et al. ), (Chamkha et al. 6), (Arafue ad Hirata 999), (Magyari ad Chamkha 8), (Aii et al. ), (Hamid ad Arifi 4), (Remeli et al. ), (Zhag ad Zheg ), (Che 7) ad (Al-Mudhaf ad Chamkha ). I real applicatios, mathematical problems are usually modeled by oliear ordiary ad partial differetial equatios such as physical ad egieerig applicatios. I geerally oliear models may ot be a eact solutios. Therefore, we try to fid approimate or umerical solutios of these models as see i refereces (Freidooimehr ad Rashidi ; Jhakal ). Oe of the most popular techique is variatioal iteratio method that is very powerful. I this cotet, we will cosider to eted the work studied by (Zhag ad Zheg 4) i order to fid aalytical solutios by usig variatioal iteratio method with combiig Pade approimatio called modified variatioal iteratio method (MVIM). Also, we will give compariso betwee MVIM ad
M. Gubes ad G. Oturac / JAFM, Vol. 9, No. 6, pp. 949-9, 6. asymptotic epasio method (Zhag ad Zheg 4). Additioally, ukow parameters of velocity gradiet f () a ad temperature gradiet () b which obtaied by MVIM is compared with Ruge-Kutta method as show i tables. The effects ad variatios of Pr, ad o velocity ad temperature profiles are preseted graphically. All the calculatios for solutios are provided by oly oe or two step iteratios. Thus, it is foud that the preset results are i very good agreemet with other kow result as preseted i Table, ad.. MATHEMATICAL DESCRIPTION OF PROBLEMS We will cosider two-dimesioal Maragoi boudary layer flow with buoyacy effects due to eteral pressure gradiet ad gravity. It's occurs alog a iterface S betwee two fluids as i Fig.. if y, u u ( ) u, T () e where u ad v are the velocity compoets correspodig to ad y ais respectively. ue ( ) is eteral velocity, T is fluid temperature, is Maragoi mied covectio parameter. Also, if, the buoyacy force is available ad if, the buoyacy force is ot available. Additioally, variatio of surface tesio as d T, mtt Tm. () dt By usig these facts ad boudary coditios, we ca write the trasform variables as 4 T y, T T ( ) ( y, ) 4 T f( ) u r, 4 T k T (4) Combiig ()-(4) ad cosiderig the literature (Golia ad Viviai 98, 986; Chamkha et al 6; Al-Mudhaf ad Chamkha ) we obtai the mai ordiary differetial system which is the reduced form of () as Fig.. Schematic of the problem. f f ff r Pr f f ad boudary coditios tur ito f (), f( ) r, f() k, () (6) Gravity g occurs throughout o iterface S, the surface tesio chages with temperature. Viscous dissipatio ad iterface tortuosity are egligible. Also the flow fields for two iterfacig fluids are idepedet (Golia ad Viviai 98, 986; Zhag ad Zheg 4). Cosiderig these iformatio, we ca write the goverig equatios for the Maragoi boudary layer with water based fluid as u v y u u du e u u v u T y e d y u T v T T y Pr y Also boudary coditios are give as u T if y,, v, T T y () (), ( ) Here, prime deotes the derivatives with respect to.. MODIFIED VARIATIONAL ITERATION METHOD (MVIM) Variatioal iteratio method (VIM) is oe of the powerful mathematical tool to solve various kids of liear ad oliear problems as show some of i ref. (He 7, 999, 997; He ad Hog 7). I order to basic defiitio of VIM, we cosider the followig geeral oliear problem (He 7, 999, 997; He ad Hog 7)?? EMBED Equatio.DSMT4 Lu( ) Ru( ) N u( ) g( ) (7) d m where L, m ad R are liear operators, d m N is a oliear operator ad g is give cotiuous fuctio. Accordig to the origially VIM, we costruct the correctio fuctioal as 9
M. Gubes ad G. Oturac / JAFM, Vol. 9, No. 6, pp. 949-9, 6. u ( ) u ( ) Lu ( s) Ru ( s) ( s) ds Nu ( s) g( s) (8) Here,, is a Lagrage multiplier (Iokuti et al. 978), u is cosidered as a restricted variatio, i.e. u. If we apply the variatio to correctio fuctioal (8) by usig variatioal aalysis, the we write dow as u ( ) u ( ) Lu ( s) Ru ( s) (, s) ds Nu ( s) g( s) Ku ( ) ( ) (, ) s u s ds g( s) (9) () From solutio of Euler-Lagrage problem show i (), we determie the Lagrage multiplier ad successive iteratios u ( ), are obtaied by usig Lagrage multiplier ad iitial approimatio u that satisfy, at least, the iitial ad boudary coditios with possible ukows. Cosequetly the eact solutio of (7) ca be obtaied by usig (He 7, 999, 997; He ad Hog 7) as u ( ) lim u( ) (). Pade Techic Some techiques eist to icrease the covergece of a give series. Amog them so-called Pade Techique is widely applied (Baker ad Morris 98). Supposed that a fuctio y ( ) is preseted by a power series c LM, Pade approimat is () L g g gl LM M O () hh hm which agree with () as far as possible. Here there are L idepedet umerator coefficiets ad M idepedet deomiator coefficiets, so makig LM ukow coefficiet i all. This is suggested that ormally LM, ought to fit the power series () amely L gg gl c M hh hm (4) L M O If the equatios equate with respect to L L LM,,...,, we write dow h c h c M LM M LM hc L hmclm hm clm hc L h c h c h c M L M L LM () If, c for cosistecy. Sice h, () become a set of M liear equatios for M ukow deomiator coefficiets ad also the umerator coefficiets g, g,..., g L follow immediately from (4) by equatig the coefficiets L M of,,,..., as g c, g cch, g c hc h c, mi LM, g c h c L L L (6) Thus, (6) ormally determie the Pade umerator ad deomiator. The LM, Pade approimat is costructed which agrees with c, through L M order. I order to fid the ifiite boudary coditios i (6) ad icrease covergece ad efficiecy of the series solutio (), we apply the Pade approimatio techic to (). Therefore, we combie the variatioal iteratio method ad pade techic so called modified variatioal iteratio method (MVIM). 4. SOLUTION PROCEDURE OF PROBLEMS Now, we will apply our proposed method MVIM to eqs. ()-(6) to obtai aalytic solutios. ad Let, assumed that f a b for the boudary coditios (6). By usig these cases, the iitial approimatios f ad which provided boudary coditios (6) are cosidered as 9
M. Gubes ad G. Oturac / JAFM, Vol. 9, No. 6, pp. 949-9, 6. f a ad b (7) where ab, are ukow coefficiets that will be obtaied by applyig boudary coditios (6). By usig the variatios theory (9)-(), Lagrage multipliers are foud as follow respectively ( ) f ad ( ) s s (8) Thus, Lagrage multipliers (8) put ito (8), the successive iteratio equatios are writte as f ( ) f ( ) f s( f s) ( s ) f ( s) f s r ds s ad s ( s ) f( s) s ds Pr f ss (9) () By applyig (7) to (9)-(), we obtai the solutios of () respect to boudary coditios (6) as follow 6 f a a 4 a b 6 4 Pr a a a a 4 6 Pr a 44 4 6 Prba 4 ( ) b Pra Pr abpr 6 ( a Pr Pr) 4 4 Pr( a ) Prb 8 () () (Pr a bpr b) Pr a b 4 Prb Pra Pr a 6 PrPr Pr a 6 7 4aPrb Prb 6 4 Pr a Pr a 4 For () ad (), iteratio process cotiues sufficietly (as see () ). I the series solutios () ad (), the ukow costat ab, which are deoted to velocity gradiet ad temperature gradiet respectively, are foud by applyig very efficiet approimatio called Pade Techique (-6) ad ifiite boudary coditios (6). For umerical values of, Pr ad (opposig buoyacy force), we fid ukow costats as f a.4746 b.886 () From () ad put, Pr ad ito (- ), the umerical solutio of () with MVIM is foud as follow respectively, f.4746.98.84 4.7788 6.4696.886.686.8446.76888 4.86.77499 6. RESULTS Ad DISCUSSION The dimesioless form () ad (6) of Maragoi boudary layer flow equatios () ad () are cosidered. These equatios have bee solved by modified variatioal iteratio method (MVIM) ad very efficiet ad accurate results are obtaied by MVIM. For both opposig ad favorable buoyacy effected Maragoi flow with various values of Pradtl umber (Pr) ad Maragoi mied covectio parameter, our results are compared with well kow asymptotic iteratio method (Zhag ad Zheg 4) ad fourth order Ruge-Kutta method i Table,,. 9
M. Gubes ad G. Oturac / JAFM, Vol. 9, No. 6, pp. 949-9, 6. Table Compariso of velocity ad temperature gradiets values for buoyacy force effects MVIM Ref (Zhag & Zheg) Numerical Pr f () () f () () f () ().4746 -.886.779 -.94868.46998947 -.78948746989.4667 -.976787667.98498-4.978.4946749794 -.974869986..6768 -.88497.8768 -.699.948 -.97689647..4846887 -.64944.4797 -.8848.4889987444 -.9988444994 MVIM Ref (Zhag & Zheg) Numerical Pr f () () f () () f () ().87 -.646677.768744 -.899.979974 -.676988967668.8667 -.94478.98999 4.67896.4779889 -.877974447.4766-4.94696767.68-7.999.67469494-4.78888688486.87 -.874.86 -.84987.64487497 -.644769.4874 -.877949.7 -.88.67964494 -.98786 Fourth order Ruge Kutta Table Compariso of Temperature ad Velocity gradiets for various values of Pr ad for MVIM Numerical Pr f () () f () ()..8646 -.799886.867897 -.7646446...8686 -.9666.7968766866 -.4896769868..787696 -.967.76994494 -.949848667648.6497 -.649446.94677 -.7467664944..44777 -.8886.498787896 -.8669849666.466 -.8987768.466786446 -.84849886 Fourth order Ruge Kutta Table Compariso of Temperature ad Velocity gradiets for various values of Pr ad for MVIM Numerical Pr f () () f () ()..96864 -.6867987.97687466 -.78797496...484466 -.7.489997666 -.489768697...8846 -.878469.647664987 -.64799488.46976 -.7.9879648 -.7449796446..478 -.877.96687746 -.6787849747799.87 -.874.64768 -.64467469889 Fourth order Ruge Kutta It is evidet from Table that our results better tha the results i ref. (Zhag ad Zheg 4) with compared umerical method. Also from Table ad, it tells us that our preseted method (MVIM) is efficiet ad powerful mathematical tool. Fig. - demostrate the variatios of velocity profiles for both opposig ad favorable buoyacy forces with differet values of Pradtl Number (Pr). Also, Fig. 4-7 show the variatios of temperature profiles as the same meaig of Fig. - with various values of Pr ad Maragoi mied covectio parameter respectively. 6. CONCLUSION I this paper, we cosider the oliear ordiary differetial equatios which correspods to Maragoi boudary layer flow with buoyacy effects. These equatios are solved by Modified variatioal iteratio method aalytically. The velocity gradiet f, temperature gradiet as well as the temperature ad velocity profiles are aalyzed ad compared for buoyacy 9
M. Gubes ad G. Oturac / JAFM, Vol. 9, No. 6, pp. 949-9, 6. opposed ad buoyacy favorable cases. Cosequetly, results show that MVIM is very powerful ad coveiet method for aalytical ad umerical solutios for oliear flow equatios. Fig.. Pradtl Number effects o temperature profile for eq. () with assistig buoyacy forces. Fig.. Pradtl Number effects o velocity profile for eq. () with opposig buoyacy forces. Fig. 6. Maragoi covectio parameter effects o temperature profile for eq. () with opposig buoyacy forces. Fig.. Pradtl Number effects o velocity profile for eq. () with assistig buoyacy forces Fig. 7. Maragoi covectio parameter effects o temperature profile for eq. () with assistig buoyacy forces. Fig. 4. Pradtl Number effects o temperature profile for eq. () with opposig buoyacy forces. ACKNOWLEDGEMENTS The authors would like to thak editor ad all referees for their cotributios. 94
M. Gubes ad G. Oturac / JAFM, Vol. 9, No. 6, pp. 949-9, 6. REFERENCES Aii, N. A., N. M. Arifi ad F. Ismail (). Radiatio effect o Maragoi covectio boudary layer flow of a aofluid. Mathematical Scieces, 6(). Al-Mudhaf, A. ad A. J. Chamkha (). Similarity solutios for MHD thermosolutal Maragoi covectio over a flat surface i the presece of heat geeratio or absorptio effects. Heat Mass Trasfer (4), -. Arafue, K. ad A. Hirata (999). Thermal ad solutal Maragoi covectio i I-Ga-Sb system. J. Crystal Growth (97), 8-87. Baker, G. A., P. (98). Graves-Morris, Padé Approimats, Ecyclopedia of Mathematics ad Its Applicatio, Parts I ad II, Addiso-Wesley Publishig Compay, New York. Chamkha, A. J., I. Pop I. ad H. S. Takhar (6). Maragoi mied covectio boudary layer flow. Meccaica (4), 9-. Chie Hsi, C. (7). Maragoi effects o forced covectio of power-law liquids i a thi film over a stretchig surface. Physics Letters A (7), -7. Christopher, D. M. ad B. X. Wag (). Similarity solutio for Maragoi covectio aroud a vapour bubble durig ucleatio ad growth. It. J. Heat Mass Trasfer (44), 799-8. Freidooimehr, N, ad M. M. Rashidi (). Dual Solutios for MHD Jeffery Hamel Nao-Fluid Flow i No-parallel Walls Usig Predictor Homotopy Aalysis Method. Joural of Applied Fluid Mechaics 8(4), 9-99. Gelles, S. H. (978). Microgravity studies i the liquid-phase immiscible system: alumiium idium. AIAA J. 6(), 4 48. Golia, C. ad A. Viviai (98). Maragoi buoyat boudary layers. L Aerotechica Missili e Spazio (64), 9. Golia, C. ad A. Viviai (986). No isobaric boudary layers related to Maragoi flows. Meccaica (), 4. Hamid, R. A. ad N. M. Arifi (4). The effect of wall suctio / ijectio a MHD Maragoi covectio boudary layer flow i aofluid. AIP Cof. Proceedigs 6, 86. He, J. H. (997). Variatioal iteratio method for delay di eretial equatios. Commu. Noliear Sci. Numer. Simulatio (4), -6. He, J. H. (999). Variatioal iteratio method-a kid of oliear aalytical techique: some eamples. Iteratioal Joural of No-Liear Mechaics (4), 699-78. He, J. H. (7). Variatioal iteratio method-some recet results ad ew iterpretatios. Joural of Computatioal ad Applied Mathematics - 7. He, J. H. ad W. Xu-Hog (7). Variatioal iteratio method: New developmet ad applicatios. Compt. Math. with Appl. (4), 88-894. Iokuti, M. ad et al. (978). Geeral use of the Lagrage multiplier i oliear mathematical physics. I: S. Nemat-Nasser (Ed.). Variatioal Method i the Mechaics of Solids, 6-6. Pergamo Press, Oford. Jhakal, A. K. (). MHD Boudary Layer Flow ear Stagatio Poit of Liear Stretchig Sheet with Variable Thermal Coductivity via He s Homotopy Perturbatio Method. Joural of Applied Fluid Mechaics 8(), 7-78. Magyari, E. ad A. J. Chamkha (8). Eact aalytical results for the thermosolutal MHD Maragoi boudary layers. It. J. Therm. Sci. (47), 848 87. Napolitao, L. G. (978). Microgravity fluid dyamics. I d Levitch coferece, Washigto. Napolitao, L. G. (979, Jue). Maragoi boudary layers. I Proceedigs of the rd Europea symposium o material sciece i space, Greoble, ESA SP-4. Napolitao, L. G. (98). Surface ad buoyacy drive free covectio. Acta Astroautica (9), 99. Okao, Y, M. Itoh ad A. Hirata (989). Natural ad Maragoi covectios i a twodimesioal rectagular ope boat. J. Chem. Eg Japa (), 7 8. Pop, I., A. Postelicu ad T. Grosa (). Thermosolutal Maragoi forced covectio boudary layers. Meccaica (6), -7. Remeli, A., N. M. Arifi, R. Nazar ad F. Ismail (). Maragoi mied covectio boudary layer flow i a aofluid. Malaysia Joural of Fudametal ad Applied Scieces 9(), 8-8. Zhag, Y. ad L. Zheg (). Aalysis of MHD thermosolutal Maragoi covectio with the heat geeratio ad a first-order chemical reactio. Chemical Egieerig Sciece (69), 449-4. Zhag, Y. ad L. Zheg (4). Similarity solutios of Maragoi covectio boudary layer flow with gravity ad eteral pressure. Chiese J. of Chem. Eg. (4), 6-69. 9