Approximate Solutions for the Coupled Nonlinear. Equations Using the Homotopy Analysis Method

Similar documents
Mixture of a New Integral Transform and Homotopy Perturbation Method for Solving Nonlinear Partial Differential Equations

On Numerical Solutions of Two-Dimensional Boussinesq Equations by Using Adomian Decomposition and He's Homotopy Perturbation Method

The modified Exp-function method and its applications to the generalized K(n,n) and BBM equations with variable coefficients

Numerical KDV equation by the Adomian decomposition method

Homotopy Analysis Method for Solving Fractional Sturm-Liouville Problems

EXISTENCE THEORY OF RANDOM DIFFERENTIAL EQUATIONS D. S. Palimkar

VARIATIONAL ITERATION METHOD: A COMPUTATIONAL TOOL FOR SOLVING COUPLED SYSTEM OF NONLINEAR PARTIAL DIFFERENTIAL EQUATIONS

The Solution of the One Species Lotka-Volterra Equation Using Variational Iteration Method ABSTRACT INTRODUCTION

BE.430 Tutorial: Linear Operator Theory and Eigenfunction Expansion

Comparison between Fourier and Corrected Fourier Series Methods

Energy Density / Energy Flux / Total Energy in 1D. Key Mathematics: density, flux, and the continuity equation.

On The Geometrıc Interpretatıons of The Kleın-Gordon Equatıon And Solution of The Equation by Homotopy Perturbation Method

Application of Homotopy Analysis Method for Solving various types of Problems of Ordinary Differential Equations

Oscillator with Distributed Nonlinear Structure on a Segment of Lossy Transmission Line

New Applications of Adomian Decomposition Method. Emad A. Az-Zo'bi

Adomian Decomposition Method and its. Modification for Nonlinear. Abel's Integral Equation

Approximating Solutions for Ginzburg Landau Equation by HPM and ADM

The analysis of the method on the one variable function s limit Ke Wu

UNIT 1: ANALYTICAL METHODS FOR ENGINEERS

A NEW NUMERICAL TECHNIQUE FOR SOLVING FRACTIONAL PARTIAL DIFFERENTIAL EQUATIONS

F D D D D F. smoothed value of the data including Y t the most recent data.

Local Fractional Variational Iteration Method for Solving Nonlinear Partial Differential Equations within Local Fractional Operators

Outline. Review Homework Problem. Review Homework Problem II. Review Dimensionless Problem. Review Convection Problem

University of Mosul. From the SelectedWorks of Mohammed O. Al-Amr

Bernstein Direct Method for Solving. Variational Problems

Available online at J. Math. Comput. Sci. 4 (2014), No. 4, ISSN:

MODIFIED ADOMIAN DECOMPOSITION METHOD FOR SOLVING RICCATI DIFFERENTIAL EQUATIONS

VISCOSITY APPROXIMATION TO COMMON FIXED POINTS OF kn- LIPSCHITZIAN NONEXPANSIVE MAPPINGS IN BANACH SPACES

BEST LINEAR FORECASTS VS. BEST POSSIBLE FORECASTS

SOLVING OF THE FRACTIONAL NON-LINEAR AND LINEAR SCHRÖDINGER EQUATIONS BY HOMOTOPY PERTURBATION METHOD

REDUCED DIFFERENTIAL TRANSFORM METHOD FOR GENERALIZED KDV EQUATIONS. Yıldıray Keskin and Galip Oturanç

Reduced Differential Transform Method for Solving Klein Gordon Equations

Moment Generating Function

Research Article A Generalized Nonlinear Sum-Difference Inequality of Product Form

Exact solitary-wave Special Solutions for the Nonlinear Dispersive K(m,n) Equations by Means of the Homotopy Analysis Method

Ruled surfaces are one of the most important topics of differential geometry. The

PLASTIC BUCKLING OF SSSS THIN RECTANGULAR PLATES SUBJECTED TO UNIAXIAL COMPRESSION USING TAYLOR-MACLAURIN SHAPE FUNCTION

Mean Square Convergent Finite Difference Scheme for Stochastic Parabolic PDEs

Fresnel Dragging Explained

CHARACTERIZATIONS OF THE NON-UNIFORM IN TIME ISS PROPERTY AND APPLICATIONS

Curvilinear Motion: Normal and Tangential Components

APPROXIMATE SOLUTION OF FRACTIONAL DIFFERENTIAL EQUATIONS WITH UNCERTAINTY

Dynamic h-index: the Hirsch index in function of time

A Novel Approach for Solving Burger s Equation

Application of Homotopy Analysis Method for Solving various types of Problems of Partial Differential Equations

1. Solve by the method of undetermined coefficients and by the method of variation of parameters. (4)

1 Notes on Little s Law (l = λw)

An approximate solution for a generalized Hirota-Satsom coupled (Kdv) equation

Solitons in a system of three linearly coupled fiber gratings

TAKA KUSANO. laculty of Science Hrosh tlnlersty 1982) (n-l) + + Pn(t)x 0, (n-l) + + Pn(t)Y f(t,y), XR R are continuous functions.

6.2 The Moment-Curvature Equations

th m m m m central moment : E[( X X) ] ( X X) ( x X) f ( x)

1.225J J (ESD 205) Transportation Flow Systems

Online Supplement to Reactive Tabu Search in a Team-Learning Problem

Department of Mathematical and Statistical Sciences University of Alberta

Improvement Over General And Wider Class of Estimators Using Ranked Set Sampling

EXTERNALLY AND INTERNALLY POSITIVE TIME- VARYING LINEAR SYSTEMS

On the Effective Region of Convergence of the Decomposition Series Solution

FIXED FUZZY POINT THEOREMS IN FUZZY METRIC SPACE

Chapter 9 Autocorrelation

A Note on Integral Transforms and Differential Equations

S n. = n. Sum of first n terms of an A. P is

Chemistry 1B, Fall 2016 Topics 21-22

For demonstration of the concept of HAM, by considering general non-linear problem (1) Non-linear operator is N and v (t) [ ],

SUMMATION OF INFINITE SERIES REVISITED

Research Article The Analytical Solution of Some Fractional Ordinary Differential Equations by the Sumudu Transform Method

Extremal graph theory II: K t and K t,t

Optimum design of complementary transient experiments for estimating thermal properties

In this section we will study periodic signals in terms of their frequency f t is said to be periodic if (4.1)

HYPOTHESIS TESTING. four steps

Notes 03 largely plagiarized by %khc

Calculus Limits. Limit of a function.. 1. One-Sided Limits...1. Infinite limits 2. Vertical Asymptotes...3. Calculating Limits Using the Limit Laws.

ECE 570 Session 7 IC 752-E Computer Aided Engineering for Integrated Circuits. Transient analysis. Discuss time marching methods used in SPICE

ME 501A Seminar in Engineering Analysis Page 1

L-functions and Class Numbers

Chapter 11 Autocorrelation

Journal of Applied Science and Agriculture

CLOSED FORM EVALUATION OF RESTRICTED SUMS CONTAINING SQUARES OF FIBONOMIAL COEFFICIENTS

Turing-Computability of Solution of Hirota Equation Dianchen Lu1, a and Liming Fu1, b

Extended Laguerre Polynomials

Some Properties of Semi-E-Convex Function and Semi-E-Convex Programming*

NEWTON METHOD FOR DETERMINING THE OPTIMAL REPLENISHMENT POLICY FOR EPQ MODEL WITH PRESENT VALUE

ES 330 Electronics II Homework 03 (Fall 2017 Due Wednesday, September 20, 2017)

A note on deviation inequalities on {0, 1} n. by Julio Bernués*

Research Article Finite Difference and Sinc-Collocation Approximations to a Class of Fractional Diffusion-Wave Equations

A quadratic convergence method for the management equilibrium model

λiv Av = 0 or ( λi Av ) = 0. In order for a vector v to be an eigenvector, it must be in the kernel of λi

NONLINEAR SCHRÖDINGER EQUATION

Numerical Solution of Parabolic Volterra Integro-Differential Equations via Backward-Euler Scheme

MALAYSIAN JOURNAL OF MATHEMATICAL SCIENCES. Boundary Value Problem for the Higher Order Equation with Fractional Derivative

Common Fixed Point Theorem in Intuitionistic Fuzzy Metric Space via Compatible Mappings of Type (K)

Ideal Amplifier/Attenuator. Memoryless. where k is some real constant. Integrator. System with memory

Electrical Engineering Department Network Lab.

Recursive Identification of Hammerstein Systems with Polynomial Function Approximation

Section 8. Paraxial Raytracing

The Connection between the Basel Problem and a Special Integral

ODEs II, Supplement to Lectures 6 & 7: The Jordan Normal Form: Solving Autonomous, Homogeneous Linear Systems. April 2, 2003

Samuel Sindayigaya 1, Nyongesa L. Kennedy 2, Adu A.M. Wasike 3

Modified Decomposition Method for Solution of Fractional Partial Differential Equations of Two-Sided

E will be denoted by n

Transcription:

Applied Maheaical Scieces, Vol. 5,, o. 37, 89-86 Approxiae Solios for he Copled Noliear Eqaios Usig he Hooopy Aalysis Mehod Spig Qia a, b a Facly of Sciece, Jiags Uiersiy, Zhejiag, Jiags 3, Chia b Depare of Maheaics, Chagsh Isie of Techology, Chagsh, Jiags 55, Chia qsp@cslg.ed.c Li Wei a a Facly of Sciece, Jiags Uiersiy, Zhejiag, Jiags 3, Chia Absrac I his aricle, he hooopy aalysis ehod (HAM) is ipleeed o obai he approxiae solios of he oliear eolio eqaios i aheaical physics. The resls obaied by his ehod hae a good agreee wih oe obaied by oher ehods. I illsraes he alidiy ad he grea poeial of he hooopy aalysis ehod i solig parial differeial eqaios. Maheaics Sbjec Classificaio: 58B5 Keywords: Hooopy aalysis ehod; oliear eolio eqaio; approxiae solios. Irodcio I is well kow ha oliear dyaical syses arise i arios fields. A

8 Spig Qia ad Li Wei wealh of ehods hae bee deeloped o fid hese exac physically sigifica solios of a parial differeial eqaio hogh i is raher difficl. Soe of he os ipora ehods are Backld rasforaio [], he CK direc ehod [,3], ad he hoogeeos balace ehod [4]. The hooopy aalysis ehod (HAM) [5], is a powerfl ehod o sole o-liear probles. Based o hooopy of opology, he alidiy of he HAM is idepede of wheher or o here exis sall paraeers i he cosidered eqaio. Therefore, he HAM ca oercoe he foregoig resricios ad liiaios of perrbaio echiqes [6]. I rece years, his ehod has bee sccessflly eployed o sole ay ypes of o-liear probles [7,8,9]. I his paper, he hooopy aalysis ehod is sed o sole he syse of he approxiae eqaios for log waer waes [,], as follows: + /= x x xx { ( ) /= x wih he soliary wae solios [4] { xx [ ] = α [( α + β + γ) ] x, = α ah αx+ β + γ / / + α /+ c x, sec h x / /4 () () where αγ,,c are arbirary cosas ad β = cα + α /.. Aalysis of he ehod For he coeiece of he reader, we will firs prese a brief acco of he HAM. Le s cosider he followig differeial eqio : [ ( x ) ] Ν, = (3) where Ν is a oliear operaor, ( x, ) is a kow fcio, ad x ad deoe spaial ad eporal idepede ariables. For sipliciy, we igore all bodary or iiial codiios, which ca be reaed i he siilar way. By eas of geeralizig he radiioal hooopy ehod, we ca cosrc he so-called zero-order deforaio eqaio as: ( q) L[ φ( x, ; q) ( x, ) ] = q Η( x, ) Ν[ φ( x, ; q) ] h (4) where q [,] is he ebeddig paraeer, h is a o-zero axiliary paraeer,

Approxiae solios 8 Η( x, ) is a axiliary fcio, L is a axiliary liear operaor, a iiial gess of ( x, ) ad ( x,; q) whe q = ad q =, i holds φ ( x,; ) = ( x, ), φ ( x,; ) = ( x, ), x is φ is a kow fcio. Obiosly, respeciely. Ths, as q icreases fro o, he solio ( x,; q) φ aries fro he iiial gess ( x, ) o he solio ( x, ). Expadig φ ( x,; q) i Taylor series wih respec o q, we hae + = + (5) (,; ) (,) φ x q x q = where ( x ) ( xq,; ) φ, =! q q= (6) If he axiliary liear operaor, he iiial gess, he axiliary paraeer h, ad he axiliary fcio are properly chose, he series (5) coerges a q =, he we hae (, ) (, ) x = + x (7) + = which s be oe of solios of origial oliear eqaios. Accordig o he defiiio (6), he goerig eqaio ca be dedced fro he zero-order deforaio (4). Defie he ecor x x x x (, ) { (, ), (, ),, (, )} = L (8) Differeiaig eqaio (4) ies wih respec o he ebeddig paraeer q ad he seig q = ad fially diidig he by!, we hae he so-called h-order deforaio eqaio L [ ( x, ) χ ( x, )] = h H( x, ) R (9)

8 Spig Qia ad Li Wei where R ( ) = ( )! Ν [ φ ( xq,; )] q q= () ad {, χ =, > I order o assess he adaages ad accracy of HAM, we cosider he followig applicaio. () 3. Applicaio we shall deal wih he syse of he approxiae eqaios for log waer waes (). For sipliciy, α =, γ =, c = are sed i he aalyses. For applicaio of hooopy aalysis ehod, we choose he iiial codiios { (,) sec = x, = ah x + = x = h x ad he axiliary liear operaors ( x q) ( xq) φ,; ϕ,; L [ φ( x,; q) ] =, L [ ϕ( x,; q )] = () (3) wih he propery L [ C ] L [ C ] =, = where C ad C are cosa coefficies, φ ad ϕ are real fcios. Frherore, we defie he oliear operaors φ x,; q φ xq,; ϕ xq,; φ xq,; Ν [ φ( xq,; ), ϕ( xq,; )] = φ( xq,; ) + x x x (4), ϕ ϕ φ ϕ Ν [ φ( xq,; ), ϕ( xq,; )] = φ( xq,; ) ϕ( xq,; ) x x x x,; q xq,; xq,; xq,; where q [, ], φ ( x,; q) ad ( x,; q) ϕ are real fcios of x, ad q. Le (5)

Approxiae solios 83 h, h deoe he ozero axiliary paraeers. Usig he aboe defiiio, wih asspio H ( x, ) =, eqaios as follows, H x, =, we cosrc he zero-order deforaio ( q) L[ φ( xq, ; ) ( x, )] q N ( φ( xq, ; ), ϕ( xq, ; )) = h, (6) ( q) L[ ϕ( xq, ; ) ( x, )] q N ( φ( xq, ; ), ϕ( xq, ; )) = h (7) obiosly, whe q = ad q =, i is clear ha ( x, ; ) ( x, ), ( x, ;) ( x, ), ( x, ;) ( x, ), ( x, ;) ( x, ) φ = ϕ = φ = ϕ = (8) Boh of h ad h are properly chose so ha he ers ( x ) ( xq,; ) φ, =! q q= ad ( x ) ( xq,; ) ϕ, =! q q= (9) exis for ad he power series of q i he followig fors () = = (,; ) = (,) + (,), ϕ(,; ) = (,) + (,) φ x q x x q xq x x q are coerge a q =. So sig (9), we obai + +. () = = (, ) = + (, ), (, ) = + (, ) x x x x Accordig o he fdaeal heore of HAM, we hae he h-order deforaio eqaio [ χ ] χ [ ] L x, x, = h R,, L x, x, = h R,,() where i (3) x, x, x, x, R x x x x i,, = + i= x, x, x, x, R x x i i, (, ) (, = i i ) i= x i= x x (4)

84 Spig Qia ad Li Wei ad χ is defied by (). Now, he solio of he h-order deforaio eqaio () for becoes ( x, ) = χ (, ),, (, ) (, ), x + h L R x χ x L R = + h. (5) Noe ha he solios series () coai wo axiliary h ad h. For sipliciy, le = = h h h, he he approxiaios of ( x, ) ad (, ) x are oly depede o h. We wrie he differeial eqaios eed o calclae,,, L, ad 3,, 3, L, as follows x x x = h h + h h x x = h + h h h = ( h+ ) h + h h + x x x x = ( h+ ) h h + + + x x x x x = 3 ( + h) h + h h + + x x x x x = 3 ( + h ) h h + + + + + x x x x x x x M (6) i = ( h+ ) h + h h i x x i= x i i = ( h+ ) h h + i i x i= x i= x The ad i i,( i =,, 3, ) L copoes hae bee obaied sig he aple package. I is ipora o esre ha he series solios are coerge. We iesigae he iflece of he axiliary paraeer h o he coergece of he

Approxiae solios 85 series by ploig he so-called h -cres. The h cres of,,,,, ad x xx x a he poi ( x, ) = ( 3.5,5.5) are depiced i Fig. ad Fig.. The alid xx regio of h is a horizoal lie sege. I is obsered he alid regio for h is.5 < h <.5 as show i Fig. ad Fig.. Fig.. h -cres of ad is deriaies wih differe h i he case of ( x, ) = ( 3.5, 5.5) ad 8h-order approxiaio. Fig.. h -cres of ad is deriaies wih differe h i he case of ( x, ) = ( 3.5, 5.5) ad 8h-order approxiaio. I order o erify he alidiy of he HAM solio, he hree-diesioal plos of he absole error bewee he exac solios ad he solios series obaied by HAM for h =.8 are show i Fig. 3 ad Fig. 4. As show i hese figres, he behaior of he approxiae solios obaied by hooopy aalysis ehod agree well wih oe obaied by he exac solios (). Fig. 3 Absole error for he 8h-order approxiaio by HAM for (, ) Fig. 4 Absole error for he 8h-order approxiaio by HAM for (, ) x ad h =.8. x ad h =.8.

86 Spig Qia ad Li Wei 4. Coclsio I his paper, he hooopy aalysis ehod (HAM) is sed o obai he approxiae solios of he syse of he approxiae eqaios for log waer waes. The resls obaied by his ehod hae a good agreee wih oe obaied by oher ehods. The adaages of he HAM are illsraed. I is easy o see ha he HAM is a ery powerfl ad efficie echiqe i fidig aalyical solios of wide classes of oliear parial differeial eqaios. Refereces [] S.P. Qia, L.X. Tia, Nolocal Lie Bäckld syeries of he copled KdV syse, Physics Leer A. 364 (7), 35-38. [] S.Y. Lo, A oe o he ew siilariy redcios of he Bossiesq eqaio, Chiese Physics Leer A, 5 (99), 33-35. [3] P.J. Oler, Applicaio of Lie Grop o differeial Eqaios, New York,spriger, 993. [4] M.L. Wag, Y.B. Zho, Z.B. Li, Applicaio of a hoogeeos balace ehod o exac solios of oliear eqaios i aheaical physics, Physics Leer A, 6 (996), 67-75. [5] S.J. Liao, The proposed hooopy aalysis echiqes for he solio of oliear probles, Ph.D. disseraio, Shaghai Jiao Tog Uiersiy, Shaghai, 99 [i Eglish]. [6] J.H. He, Hooopy perrbaio ehod, Coper Mehods i Applied Mechaics Egieerig 78 (999), 57-6. [7] S.J. Liao, O he hooopy aalysis ehods for oliear probles, Applied Maheaics Copaio, 47 (4), 499-53. [8] S. Abbasbady, The applicaio of hooopy aalysis ehod o sole a geeralized Hiroa Sasa copled KdV eqaio, Physics Leer A, 36 (7), 478-483. [9] T. Haya, M. Sajid, O aalyic solio for hi fil flow of a forh grade flid dow a erical cylider, Physics Leer A,36 (7), 36-3. [] G.B. Whiha, Variaioal ehods ad applicaios o waer wae, Proceedigs of he Royal Sociey A, 99 (967), 6-5. [] L.T.F. Broer, Approxiae eqaios for log waer waes, Applied Scieific Research, 3 (975), 377-395. Receied: Noeber,