THE CONSERVATIVE DIFFERENCE SCHEME FOR THE GENERALIZED ROSENAU-KDV EQUATION

Similar documents
Numerical Methods for Finding Multiple Solutions of a Superlinear Problem

Partial Differential Equations

Interactions of Soliton Waves for a Generalized Discrete KdV Equation

EXISTENCE OF CONCENTRATED WAVES FOR THE DIFFERENT TYPE OF NONLINEARITY. V. Eremenko and N. Manaenkova

LINEARIZATION OF NONLINEAR EQUATIONS By Dominick Andrisani. dg x. ( ) ( ) dx

Three-Step Iterative Methods with Sixth-Order Convergence for Solving Nonlinear Equations

MAT2400 Assignment 2 - Solutions

LECTURE 13 SPURIOUS REGRESSION, TESTING FOR UNIT ROOT = C (1) C (1) 0! ! uv! 2 v. t=1 X2 t

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.

too many conditions to check!!

Introduction. Question: Why do we need new forms of parametric curves? Answer: Those parametric curves discussed are not very geometric.

Approximation of the Likelihood Ratio Statistics in Competing Risks Model Under Informative Random Censorship From Both Sides

Open problem in orthogonal polynomials

In this document, if A:

Application of Digital Filters

Fluids Lecture 17 Notes

Stability of Solution for Nonlinear Singular Systems with Delay

Approximate Solutions of Set-Valued Stochastic Differential Equations

A collocation method for singular integral equations with cosecant kernel via Semi-trigonometric interpolation

MATHEMATICS I COMMON TO ALL BRANCHES

Research Article Two Expanding Integrable Models of the Geng-Cao Hierarchy

Definition 2.1 (The Derivative) (page 54) is a function. The derivative of a function f with respect to x, represented by. f ', is defined by

Sequences A sequence of numbers is a function whose domain is the positive integers. We can see that the sequence

The z Transform. The Discrete LTI System Response to a Complex Exponential

Research Article A Note on Ergodicity of Systems with the Asymptotic Average Shadowing Property

AN OPEN-PLUS-CLOSED-LOOP APPROACH TO SYNCHRONIZATION OF CHAOTIC AND HYPERCHAOTIC MAPS

62. Power series Definition 16. (Power series) Given a sequence {c n }, the series. c n x n = c 0 + c 1 x + c 2 x 2 + c 3 x 3 +

Comparison Study of Series Approximation. and Convergence between Chebyshev. and Legendre Series

The Comparison Adomian Decomposition Method and Differential Quadrature Method for Solving Some Nonlinear Partial Diferential Equations

The Modification of BCC Model Using Facet Analysis

Gain scheduled observer state feedback controller for rational LPV systems

Alternating Series. 1 n 0 2 n n THEOREM 9.14 Alternating Series Test Let a n > 0. The alternating series. 1 n a n.

PAijpam.eu ON TENSOR PRODUCT DECOMPOSITION

Uniform Strict Practical Stability Criteria for Impulsive Functional Differential Equations

Definition 4.2. (a) A sequence {x n } in a Banach space X is a basis for X if. unique scalars a n (x) such that x = n. a n (x) x n. (4.

NUMBER OF SPANNING TREES OF NEW JOIN GRAPHS

IJITE Vol.2 Issue-11, (November 2014) ISSN: Impact Factor

Estimation of Backward Perturbation Bounds For Linear Least Squares Problem

which are generalizations of Ceva s theorem on the triangle

An almost sure invariance principle for trimmed sums of random vectors

Some q-analogues of Fibonacci, Lucas and Chebyshev polynomials with nice moments

Numerical Solution of the First-Order Hyperbolic Partial Differential Equation with Point-Wise Advance

Linear Differential Equations of Higher Order Basic Theory: Initial-Value Problems d y d y dy

A sequence of numbers is a function whose domain is the positive integers. We can see that the sequence

A Note on Generalization of Semi Clean Rings

5.6 Absolute Convergence and The Ratio and Root Tests

INFINITE SEQUENCES AND SERIES

NUMERICAL METHOD FOR SINGULARLY PERTURBED DELAY PARABOLIC PARTIAL DIFFERENTIAL EQUATIONS

MONOTONE DIFFERENCE SCHEMES STABILIZED BY DISCRETE MOLLIFICATION FOR STRONGLY DEGENERATE PARABOLIC EQUATIONS

DECOMPOSITION METHOD FOR SOLVING A SYSTEM OF THIRD-ORDER BOUNDARY VALUE PROBLEMS. Park Road, Islamabad, Pakistan

Newton Homotopy Solution for Nonlinear Equations Using Maple14. Abstract

Chapter Vectors

Sequences, Mathematical Induction, and Recursion. CSE 2353 Discrete Computational Structures Spring 2018

wavelet collocation method for solving integro-differential equation.

6a Time change b Quadratic variation c Planar Brownian motion d Conformal local martingales e Hints to exercises...

Solution of Differential Equation from the Transform Technique

Castiel, Supernatural, Season 6, Episode 18

ECE-S352 Introduction to Digital Signal Processing Lecture 3A Direct Solution of Difference Equations

Finite Difference Analysis of 2-Dimensional Acoustic Wave with a Signal Function

MAT1026 Calculus II Basic Convergence Tests for Series

On a class of convergent sequences defined by integrals 1

Product measures, Tonelli s and Fubini s theorems For use in MAT3400/4400, autumn 2014 Nadia S. Larsen. Version of 13 October 2014.

6. Cox Regression Models. (Part I)

WELL-BALANCED SCHEMES FOR CONSERVATION LAWS WITH SOURCE TERMS BASED ON A LOCAL DISCONTINUOUS FLUX FORMULATION

Applied Mathematics Letters. Asymptotic distribution of two-protected nodes in random binary search trees

Research Article Some E-J Generalized Hausdorff Matrices Not of Type M

Numerical Astrophysics: hydrodynamics

Math 2784 (or 2794W) University of Connecticut

Discrete-Time Systems, LTI Systems, and Discrete-Time Convolution

Riesz-Fischer Sequences and Lower Frame Bounds

Decoupling Zeros of Positive Discrete-Time Linear Systems*

1.3 Convergence Theorems of Fourier Series. k k k k. N N k 1. With this in mind, we state (without proof) the convergence of Fourier series.

Unit 6: Sequences and Series

Similarity Solutions to Unsteady Pseudoplastic. Flow Near a Moving Wall

On forward improvement iteration for stopping problems

A convergence result for the Kuramoto model with all-to-all coupling

REGULARIZATION OF CERTAIN DIVERGENT SERIES OF POLYNOMIALS

Number of Spanning Trees of Circulant Graphs C 6n and their Applications

Chapter 4. Fourier Series

Universal source coding for complementary delivery

CHAPTER 5 SOME MINIMAX AND SADDLE POINT THEOREMS

A 2nTH ORDER LINEAR DIFFERENCE EQUATION

On Arithmetic Means of Sequences Generated by a Periodic Function

Finite Difference Approximation for First- Order Hyperbolic Partial Differential Equation Arising in Neuronal Variability with Shifts

Periodic solutions for a class of second-order Hamiltonian systems of prescribed energy

Physics 116A Solutions to Homework Set #1 Winter Boas, problem Use equation 1.8 to find a fraction describing

Applying Differential Transformation Method to. the One-Dimensional Planar Bratu Problem

4. Linear Classification. Kai Yu

A new iterative algorithm for reconstructing a signal from its dyadic wavelet transform modulus maxima

Zeros of Polynomials

A remark on p-summing norms of operators

Chapter 6 Infinite Series

The Choquet Integral with Respect to Fuzzy-Valued Set Functions

Physics 324, Fall Dirac Notation. These notes were produced by David Kaplan for Phys. 324 in Autumn 2001.

Section 5.5. Infinite Series: The Ratio Test

Subject: Differential Equations & Mathematical Modeling -III. Lesson: Power series solutions of Differential Equations. about ordinary points

SEQUENCES AND SERIES

μ are complex parameters. Other

NEW FAST CONVERGENT SEQUENCES OF EULER-MASCHERONI TYPE

The Perturbation Bound for the Perron Vector of a Transition Probability Tensor

Transcription:

Zho, J., et al.: The Coservative Differece Scheme for the Geeralized THERMAL SCIENCE, Year 6, Vol., Sppl. 3, pp. S93-S9 S93 THE CONSERVATIVE DIFFERENCE SCHEME FOR THE GENERALIZED ROSENAU-KDV EQUATION by J ZHOU a*, Maobo ZHENG b, ad Re-Xi JIANG c a School of Mathematics ad Statistics, Yagtze Normal Uiversity, Chogqig, Chia b Chegd Techological Uiversity, Chegd, Chia c Eglish Teachig ad Research Grop, Forteeth Middle School of Flig, Chogqig, Chia Itrodctio Origial scietific paper DOI:.98/TSCI6S393Z I this paper, merical soltios for the geeralized Rosea-KdV eqatio are cosidered via the eergy ad mometm coservative o-liear implicit fiite differece scheme. Uiqe eistece of the coservative properties of the soltios for the differece scheme is show. Nmerical reslts demostrate that the scheme is efficiet ad reliable. Key words: geeralized Rosea-KdV eqatio, differece scheme, properties, merical eperimet The well-kow Korteweg-de Vries (KdV) eqatio [-4]: t + + = () has bee sed to describe the wave propagatio ad spread iteractio. I this paper, we cosider the followig geeralized Rosea-KdV eqatio: p t t + + + + ( ) = () where p is a iteger. Whe p eq. () is called as sal Rosea-KdV eqatio. I [, ], the solitary soltios for the geeralized Rosea-KdV eqatio with sal solitary asatze method were discssed ad the two ivariats for the geeralized Rosea- KdV eqatio were give. I [], the two types of solito soltio, i. e., solitary wave soltio ad siglar solito, were researched. Frthermore, they also sed pertrbatio theory ad semi-variatio priciple to stdy the pertrbed geeralized Rosea-KdV eqatio aalytically. I [3], asatze method was applied to obtai the topological solito soltio or shock soltio of this eqatio. Moreover, three methods, that is, asatze method, G'/G-epasio method as well as the ep-fctio method were applied to etract a few more soltios to this eqatio i [4]. Bt the merical method to the iitial-bodary vale problem of geeralized Rosea-KdV eqatio has ot bee stdied till ow. I [5, 6], two coservative differece schemes for the geeralized Rosea-KdV eqatio were proposed. Bt their schemes ca oly preserve oe coservative law. * Correspodig athor; e-mail: flzzklm@6.com

S94 Zho, J., et al.: The Coservative Differece Scheme for the Geeralized THERMAL SCIENCE, Year 6, Vol., Sppl. 3, pp. S93-S9 I this paper, we propose a coservative o-liear Crak-Nicolso-implicit differece scheme for the eqatio. The stdies show that the scheme does ot eed to select aother scheme to help iitial comptatio sch as the average liear scheme i [5]. It shold be oted that the merical simlatios show that the scheme preserves two coservative ivariats, which is better tha those reslts i [5, 6]. Hece, i this paper, we propose a coservative two-level o-liear implicit fiite differece scheme for the geeralized Rosea-KdV eq. () with the bodary coditios: ad iitial coditio: [, ]: X (,) t = X (,) t ( X,) t = ( X,) t l r l r ( X,) t = ( X,) t t [, T] (3) l r (,) = ( ) (4) The iitial-bodary vale problem presets the followig coservative properties Xr Xr (5) X X Mt ( ) = d = d = M() l l Xr (6) L L X Et ( ) = ( + )d = + = E() l Whe Xl, Xr, the iitial-bodary vale problem, eqs. ()-(4), ad the Cachy problem, eq. (), are cosistet. Coservative implicit differece scheme I this sectio, we first give some otatio sed i this paper, ad propose the coservative differece scheme for the problem of eqs. ()-(4). As sal, deote = X l + h, t = τ, J, N, where h = (X r X l )/J ad τ are the iform the spatial ad temporal step size, respectively. Let ( h, τ ), h J J+ Z = { = ( ) = = = J + }. Throghot this paper, we deote C as a geeral costat idepedet of h ad τ. Defie the differece operators, ier prodct ad orms are: t ( ) ( ) + + + + h + t ( )ˆ t ( ) h + t, J = v = h v, p p i p i ( ) ( )ˆ =,, + h + + h = ma J I view of ( ) = /( + p) Σ ( ) [7], we ca costrct the followig coservative implicit fiite differece scheme for the problems, eqs. ()-(4): i ( ) ( ) p + / + / + / + / t ˆ ˆ t p ( ) ( ) ( p i ) ( ) + + + + = (7) + ˆ

Zho, J., et al.: The Coservative Differece Scheme for the Geeralized THERMAL SCIENCE, Year 6, Vol., Sppl. 3, pp. S93-S9 S95 = ( ) J (8) = J ( ) ˆ = ( J) ˆ ( ) = ( J ) = (9), h Lemma.. [8] For ay two mesh fctios v Z, oe has: v, = v,, ˆ, v= v, ˆ, v, =, v () The we have: J,, =, = () Frthermore, if ( ) = ( ) the: = () To prove the eistece of soltio for scheme, eqs. (7)-(9), the followig Browder fied poit theorem shold be itrodced. For the proof, see [9]. Lemma.. (Browder fied poit theorem) Let H be a fiite dimesioal ier prodct space. Sppose that g: H H is cotios ad there eists a α sch that * * * g ( ), >, H, = α. The there eists H sch that g ( ) = ad α. Theorem.3. There eists Z h satisfyig the differece scheme (7)-(9). Proof: For N, we assme that,,, Z h satisfy the differece scheme (.)-(.3). Net we prove that there eists + satisfyig eqs. (7)-(9). Defie a operator g o Z as follows: h p τ i p i τ ˆ τ + + p (3) + ˆ gv ( ) = v + v + v + v + ( v ) ( v ) = By comptig the ier prodct of eq, (3) with v, we get: Therefore: vˆ, v vˆ, v p i p i v ( v ), ˆ v = + p gv ( ), v = v, v + v, v v v + v v v v v v + + + v + + v v + It is obvios that gv ( ), v, for all v Z h with v = + +. It follows * from Lemma. that there eists v Z h sch that g(v * ) =. Let + = v *, ad it ca be proved that + is the soltio of the scheme (7)-(9). (4)

S96 Zho, J., et al.: The Coservative Differece Scheme for the Geeralized THERMAL SCIENCE, Year 6, Vol., Sppl. 3, pp. S93-S9 H [ Xl, Xr ], the the soltio of the iitial- Lemma.4. [6] Sppose that bodary vale problem ()-(4) satisfies: L C, L C, C (5) Theorem.5. Sppose that H [ Xl, Xr ], the the schemes (7)-(9) are coservative for discrete mometm ad eergy, that is: = J = = = = M h M M (6) = + = = = E E E (7) Proof: Mltiplyig eq. (7) with h ad smmig p for from to J, from the bodary coditio i eq. (9), ad lemma., we get: where J + ( ) = (8) = h Therefore, eq. (6) is easily gotte from eq. (8). By comptig the ier prodct of eq. (7) with + (i. e. + + ), we have: h + + + J + + + + + ( + ) = τ ˆ ˆ + + ( + ) ( ) + + P = τ i p i p + + P = p + ˆ By the defiitio of ( ) t, it follows from the first term of eq. (9) that: J + ( + ) = = τ ( + ) () τ From eq. (), it follows from the secod ad the third term of eq. (9) that: J + + J + + + = + = = = ˆ J + + J + + = + = () = = (9)

Zho, J., et al.: The Coservative Differece Scheme for the Geeralized THERMAL SCIENCE, Year 6, Vol., Sppl. 3, pp. S93-S9 S97 Similarly: = () J + + = ˆ With the help of the bodary coditio i eqs. (9) ad (), it follows from the forth term of eq. (9) that: J + + + [( ) ( ) ] = = τ τ (3) I view of eq. (), it follows from the last term of eq. (9) that: P, i p i + J p h + + + = p + = ˆ h = p + h = p + i+ p i J p + + = p i i+ J p + + + = ˆ Let i = p ( i+ ). Obviosly, if i =, the i = p. If i = p, the i =. It follows from eq. (4) that: i p i + J h + + + + = = p + = i = p ˆ ˆ (4) P, P, (5) Therefore: P, + = (6) By the previos reslts of eqs. ()-(3) ad eq. (6), we have: + + + = The, by the defiitio of E, eq. (7) holds, which implies that the differece scheme is coservative for eergy. I order to prove the bodedess ad the coservative law of the merical soltios, we lead ito the followig lemma [8]. Lemma.6 (Discrete Sobolev s ieqality) There eist two costat C ad C sch that: (7) C C + (8)

S98 Zho, J., et al.: The Coservative Differece Scheme for the Geeralized THERMAL SCIENCE, Year 6, Vol., Sppl. 3, pp. S93-S9 which yields: Theorem.7. Sppose H [ Xl, Xr ], the the soltios C Proof: It follows from eq. (7) that:, C C (,, N) C, By Lemma. ad Schwartz ieqality, we get: + C From Lemma.6, we have C (,, N). Nmerical eperimets of eqs. (7)-(9) satisfy: C (9) I this sectio, we preset some merical eperimets to verify theoretical reslts obtaied i previos sectios. We ow cosider two cases: p = 3 ad p = 5, respectively. Whe p = 3, the solito soltio is: 5 + 4 t (, ) = 5 + 3 4sech (5 4) t 4 4 + ad the iitial coditio is: (3) (3) 5 + 4 (,) = 5 + 3 4sech (3) 4 4 Whe p = 5, the solito soltio is: t (, ) 4 4 ( 5 34)sech 5 34 = + + (5 34) t 5 3 + ad the iitial coditio is: (,) = ( 5 + 34) 5 sech 5 + 34 (34) 3 4 4 First, we simlate the wave graph of the merical soltio of the implicit o-liear scheme eqs. (7)-(9). The compariso of merical soltio (, t ) betwee differet time step ad space step at varios times are give i fig. whe p = 3. Similarly, we ca get the almost same figres for differet time step ad space step at differet times whe p = 5, respectively. (33) Figre. Wave graph of (, t) at varios times whe p = 3 ad τ = h =.5

Zho, J., et al.: The Coservative Differece Scheme for the Geeralized THERMAL SCIENCE, Year 6, Vol., Sppl. 3, pp. S93-S9 S99 Meawhile, we also list the coservatio ivariats M ad E at differet time i tabs. ad whe p = 3. Similarly, we ca get the coservatio ivariats M ad E at differet time whe p = 5. Table. The mometm of differet time i differet time step ad space step whe p = 3 (h, τ) T = s T = s T = 3 s T = 4 s (/4, /4) 4.8989794854547 4.89897948495597 4.8989796576935 4.898977554546477 (/8, /8) 4.898979485364376 4.898979485795 4.8989795335864 4.8989789393 (/6, /6) 4.89897948537968 4.8989794857347 4.89897949535347 4.8989793555465 (/3, /3) 4.898979484544 4.898979485494 4.8989794888958 4.8989793985586 Whe p = 5. M = 7.9364389 ad E = 3.793878834 at differet time. These reslts also verify that the proposed scheme is coservative for two qatities M ad E. Table. The eergy of differet time i differet time step ad space step whe p = 3 (h, τ) T = s T = s T = 3 s T = 4 s (/4, /4).6858993335.685899339.68589933636.68589933979 (/8, /8).6854385683.6854385575.68543854444.68543854348 (/6, /6).685466653.685466935.6854669785.685466948 (/3, /3).685474939659.685474956936.68547495648.68547495885 Coclsio I this paper, we costrcted o-liear-implicit fiite differece scheme for the geeral Rosea-KdV eqatio ad ivestigated some properties of its merical soltio. We proved that the o-liear scheme preserved the discrete mass ad eergy coservatio, respectively. The proposed scheme is the coditioally stable ad secod-order covergece by the discrete eergy method. The reslts show that the scheme is reliable ad efficiet. Refereces [] Ami, E., Solitary Wave Soltios for Geeralized Rosea-KdV Eqatio, Commicatios i Theoretical Physics, 55 (), 3, pp. 396-398 [] Polia, R., et al., Pertrbatio of Dispersive Shallow Water Waves, Ocea Egieerig, 63 (3), May, pp. -7 [3] Saha, A., Topological -Solito Soltios for the Geeralized Rosea-Kdv Eqatio, Fdametal Joral of Mathematical Physics, (),, pp. 9-3 [4] Ghodrat, E., et al., Topological Solitos ad Other Soltios of the Rosea-KdV Eqatio with Power Law No-Liearity, Romaia Joral of Physics, 58 (3), -, pp. - [5] Zheg, M. B., et al., A Average Liear Differece Scheme for the Geeralized Rosea-KdV Eqatio, Joral of Applied Mathematics, 4 (4), ID 793 [6] Lo, Y., et al., Coservative Differece Scheme for Geeralized Rosea-KdV Eqatio, Advaces i Mathematical Physics, 4 (4), ID 98698 [7] Zo, J. M., et al., A New Coservative Differece Scheme for the Geeral Rosea-RLW Eqatio, Bodary Vale Problems, (), ID 566

S9 Zho, J., et al.: The Coservative Differece Scheme for the Geeralized THERMAL SCIENCE, Year 6, Vol., Sppl. 3, pp. S93-S9 [8] Zho, Y., Applicatios of Discrete Fctioal Aalysis to the Fiite Differece Method, Iteratioal Academic Pblishers, Beiig, Chia, 99 [9] Browder, F. E., Eistece ad Uiqeess Theorems for Soltios of No-Liear Bodary Vale Problems, Proceedigs, Symposia i Applied Mathematics, New York, USA, 7 (965), 6, pp. 4-49 Paper sbmitted: December 5, 5 Paper revised: Febrary 5, 6 Paper accepted: March 5, 6