Numerical KDV equation by the Adomian decomposition method

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

On the Effective Region of Convergence of the Decomposition Series Solution

SUMMATION OF INFINITE SERIES REVISITED

Approximating Solutions for Ginzburg Landau Equation by HPM and ADM

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

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

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

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

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

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

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

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

Homotopy Analysis Method for Solving Fractional Sturm-Liouville Problems

UNIT 1: ANALYTICAL METHODS FOR ENGINEERS

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

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

1.225J J (ESD 205) Transportation Flow Systems

MODIFIED ADOMIAN DECOMPOSITION METHOD FOR SOLVING RICCATI DIFFERENTIAL EQUATIONS

Fixed Point Theorems for (, )-Uniformly Locally Generalized Contractions

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

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

Calculus BC 2015 Scoring Guidelines

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

A Novel Approach for Solving Burger s Equation

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

A Note on Integral Transforms and Differential Equations

VIM for Determining Unknown Source Parameter in Parabolic Equations

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

A NONLINEAR FOURTH-ORDER DIFFUSION-BASED MODEL FOR IMAGE DENOISING AND RESTORATION

6/10/2014. Definition. Time series Data. Time series Graph. Components of time series. Time series Seasonal. Time series Trend

Chapter 9 Autocorrelation

TAYLOR AND MACLAURIN SERIES

Chapter 11 Autocorrelation

Electrical Engineering Department Network Lab.

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

CS537. Numerical Analysis and Computing

Journal of Quality Measurement and Analysis JQMA 12(1-2) 2016, Jurnal Pengukuran Kualiti dan Analisis

Introduction to Earthquake Engineering Response Analysis

Extended Laguerre Polynomials

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

VARIOUS phenomena occurring in the applied sciences

STK4080/9080 Survival and event history analysis

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

Application of the Adomian Decomposition Method (ADM) and the SOME BLAISE ABBO (SBA) method to solving the diffusion-reaction equations

CS321. Numerical Analysis and Computing

(1) f ( Ω) Keywords: adjoint problem, a posteriori error estimation, global norm of error.

Comparison between Fourier and Corrected Fourier Series Methods

Four equations describe the dynamic solution to RBC model. Consumption-leisure efficiency condition. Consumption-investment efficiency condition

Single Degree of Freedom System Free Vibration

Inference of the Second Order Autoregressive. Model with Unit Roots

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

Single Degree of Freedom System Free Vibration

Mean Square Convergent Finite Difference Scheme for Stochastic Parabolic PDEs

Department of Mathematical and Statistical Sciences University of Alberta

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

Relations on the Apostol Type (p, q)-frobenius-euler Polynomials and Generalizations of the Srivastava-Pintér Addition Theorems

CHAPTER 2. Problem 2.1. Given: m k = k 1. Determine the weight of the table sec (b)

Stationarity and Unit Root tests

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

BEST LINEAR FORECASTS VS. BEST POSSIBLE FORECASTS

Approximately Quasi Inner Generalized Dynamics on Modules. { } t t R

Parametric Iteration Method for Solving Linear Optimal Control Problems

THE CONSERVATIVE DIFFERENCE SCHEME FOR THE GENERALIZED ROSENAU-KDV EQUATION

Linear System Theory

ME 321 Kinematics and Dynamics of Machines S. Lambert Winter 2002

Some Variants of Newton's Method with Fifth-Order and Fourth-Order Convergence for Solving Nonlinear Equations

Economics 8723 Macroeconomic Theory Problem Set 2 Professor Sanjay Chugh Spring 2017

If boundary values are necessary, they are called mixed initial-boundary value problems. Again, the simplest prototypes of these IV problems are:

B. Maddah INDE 504 Simulation 09/02/17

FRACTIONAL VARIATIONAL ITERATION METHOD FOR TIME-FRACTIONAL NON-LINEAR FUNCTIONAL PARTIAL DIFFERENTIAL EQUATION HAVING PROPORTIONAL DELAYS

Outline. simplest HMM (1) simple HMMs? simplest HMM (2) Parameter estimation for discrete hidden Markov models

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

Optimization of Rotating Machines Vibrations Limits by the Spring - Mass System Analysis

Sequences, Series, and All That

Lecture 15 First Properties of the Brownian Motion

METHOD OF THE EQUIVALENT BOUNDARY CONDITIONS IN THE UNSTEADY PROBLEM FOR ELASTIC DIFFUSION LAYER

LINEAR APPROXIMATION OF THE BASELINE RBC MODEL JANUARY 29, 2013

EXISTENCE THEORY OF RANDOM DIFFERENTIAL EQUATIONS D. S. Palimkar

BAYESIAN ESTIMATION METHOD FOR PARAMETER OF EPIDEMIC SIR REED-FROST MODEL. Puji Kurniawan M

A Study On (H, 1)(E, q) Product Summability Of Fourier Series And Its Conjugate Series

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

λ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

BE.430 Tutorial: Linear Operator Theory and Eigenfunction Expansion

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

Taylor Series Prediction of Time Series Data with Error Propagated by Artificial Neural Network

Curvilinear Motion: Normal and Tangential Components

Moment Generating Function

On Numerical Solution of Boundary Integral Equations of the Plane Elasticity Theory by Singular Integral Approximation Methods

Big O Notation for Time Complexity of Algorithms

Solutions to Problems 3, Level 4

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

Power Bus Decoupling Algorithm

LINEAR APPROXIMATION OF THE BASELINE RBC MODEL SEPTEMBER 17, 2013

A Note on Random k-sat for Moderately Growing k

Compact Finite Difference Schemes for Solving a Class of Weakly- Singular Partial Integro-differential Equations

Dynamic model of large amplitude vibration of a uniform cantilever beam carrying an intermediate lumped mass and rotary inertia

and the sum of its first n terms be denoted by. Convergence: An infinite series is said to be convergent if, a definite unique number., finite.

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

The Connection between the Basel Problem and a Special Integral

BIBECHANA A Multidisciplinary Journal of Science, Technology and Mathematics

Transcription:

America Joral o oder Physics ; () : -5 Pblished olie ay (hp://wwwsciecepblishiggropcom/j/ajmp) doi: 648/jajmp merical KDV eqaio by he Adomia decomposiio mehod Adi B Sedra Uiversié Ib Toail Faclé des Scieces Dépareme de Physiqe LHESI Kéira orocco Email address: msedra@icpi ( B Sedra) To cie his aricle: Adi B Sedra merical KDV Eqaio by he Adomia Decomposiio ehod America Joral o oder Physics Vol o pp -5 doi: 648/jajmp Absrac: Usig he Adomia decomposiio mehod (AD) we prese i his paper a merical approximaio o he solio o he oliear KDV eqaio The pricipal as cocers esseially he compaio o he Adomia polyomials or his ype o eqaio ad hereaer deermiig a sigiica crierio o esre he codiios or covergece o he mehod Keywords: KDV Eqaio merical Approach Adomia Decomposiio AD ehod or The KDV Eqaio Cosider he ollowig ormlaio o he KdV eqaio: + 6 + x x ( x ) x ( ) which ca be rewrie as ollows: () 6 ( ) ) ( x) () where / x represes he liear operaor o he eqaio; ad ( ) / is he o-liear cio Accordig o he Adomia decomposiio he solio is expressed as: ad he o-liear par by: wih: ) ) () A ( ) (4) d i A λ i por! dλ Ν (5) i λ By iegraig wih respec o ime ad sig he iiial codiios we have: So: + [ L 6 ( ) ] ) ( x) ds (6) ) ( x) L s) + 6 A ds (7) For he KdV eqaio he Adomia polyomials ca be expressed as ollows: A i i i This allows s o dedce ) amely: ) ( x) + ) [ + A ] wih he ollowig iiials codiios: ds (8) (9)

Adi e al: merical KdV eqaio by he Adomia decomposiio mehod ) sech We proceed i he ollowig o comp ) or [ ] approximae solio ~ ) ollowig ormla: x ad o be able o deermie he e ( x ) A () ad p he h order sig he ~ ) i ( x) i ) L i i s) + 6 A i ds i () Ths he approximae solio is : ~ ( x ) i ) merical esls i We will develop i he ollowig a compariso o he graphical preseaios bewee diere ieraios o he approximae solio o a give order ad he exac solio i a rage o space [ ] ad wo diere iervals o 5 5 as show below : ime [ ] ad [ ] () Fcio Graphic Preseaio Graphic Preseaio Based o graphical preseaios preseed above i shold be oed he ollowig observaios: The asympoic resls o he approximae solio rom abo 5 h order divergig a he growig ed o he ime

America Joral o oder Physics () : -5 parameer; Good agreeme bewee he exac solio ad he approximae solio i he redced ierval ime [ 5] To emphasize he secod remar we proceed o he graphical preseaio o wo cios: ( ) sech + 8 48 7 769 + 576 864 54888 4 ɶ ( ) () 6 8 ( )! The Taylor series expasio o ( ) ( ) which coverges or all have:! or : ( ) ( ) (5) (6) ad more we ( )!! ρ (7) We deie he cio: g ( ρ) ρ (8) where: The las preseaio clearly shows ha i he eighborhood o he redced ierval ime here are almos coicidece bewee he wo graphs ad ha beyod i he approximae solio diverges i a ssaied maer compared wih he exac solio which ca esimae he rae bewee he maximm vales o he wo cases o 4 imes We propose o calclae a coveie ime vale desigaed by τ rom which he approximae solio ~ ) begis o diverge Covergece Aalysis o he AD ehod or he KDV Eqaio We cosider ha : is a real aalyic cio i a disc B ( ) whose radis is ad L: veriies E ( ) C where C ad edow o he orm ( x) ( ) sp The ormlaio o Adomia mehod eables s o wrie: ( ) E( s) A ( ( s) ( s) ) ds (4) + Accordig o he Cachy esimaio here exiss e sch ha: Ths we have: ρ (9) ( ) ( g ) ( ) Usig he ormlaio o Adomia mehod we id: ( C ) () + () ( s) A ds C E A ds () Cg ( ) Cρ ( ) ( s) A ds C E A ds () C g ( ) g( ) C ρ ( ) + C + ( + ( ) ) + ρ ( ) (4)! ( + ) This allows s o dedce ha he Adomia series a are boded by he coverge series o powers based o

4 Adi e al: merical KdV eqaio by he Adomia decomposiio mehod ρ ( ) C or all [ τ ] τ i or case sdy C So we have : ( ) ( x) ) sech ad sig he Cachy esimae made earlier: ( )! We explici or he case : ad we ae: Le he solio ( ) ( )! x ( ) where (5) (6) (7) 6 (8) E be he operaor associaed wih he liear Cachy problem as ollows: sch as ( ) have: ( ) v E v Lv v( ) ad based o he above we ( ) C (9) E () Ths by ideiyig ad ollowig he expressio o he Forier rasorm o he solio ( x ) ˆ iξ ha: E( ) e v ca be dedced O he oher had sig he Parseval eqaliy we have: E + ( ) ( + ) ξ π π π + + ( + ξ ) ( + ξ ) Thereby veriyig he ieqaliy: wih : C i his case We have also : Eˆ( ) ˆ Eˆ( ) ˆ ˆ ( ) C () E (4) or by deiiio we ca wrie: + E( ) + ( C + ) ( C ) (5) + (6) The deermiaio o he hreshold covergece i ime eqal o τ he expressed by: C τ 6 (7) I he case o he KdV eqaio he liear operaor is: L ad as previosly meioed sig he Forier rasorm or he resolio o he liear problem we id: v + π iξ + iξx ) e ˆ ( ξ ) ) where: () ˆ + π iξx ( ξ ) e ( x) dx ξ () Aer calclaio we id ha:

America Joral o oder Physics () : -5 5 τ 5 Perormig he graphical compariso bewee he exac solio ad he approximaed oe a he TH order we id ha he deachme bewee he wo crves pracically sars rom he hreshold vale calclaed ad illsraed i he ollowig graphic: Throgh his wor i was possible o apply he AD mehod or he resolio o he KdV eqaio sig a polyomial o ime o TH order Solio ha has prove eecive i a give ime ierval Ths we have compleed or sdy wih a aalysis which allowed s o deermie he hreshold o covergece o he said solio eereces [] K Abbaoi ad Y Cherral Covergece o adomia s mehod applied o diereial eqaios Comp ah Appl (994) [] K Abbaoi ad Y Cherral Covergece o Adomia s mehod applied o oliear eqaios ah Comper odel (994) [] Himo K Abbaoi ad Y Cherral ew resls o covergece o Adomia s mehod Kyberees (999) [4] Y Zh Q Chag ad S W A ew algorihm or calclaig Adomia polyomials Appl ah Comp (5) 4 Coclsio [5] Ic O merical solio o parial diereial eqaio by he decomposiio mehod Kragjevac J ah (4)