Analytical Solution of A Differential Equation that Predicts the Weather Condition by Lorenz Equations Using Homotopy Perturbation Method

Similar documents
A Study of the Solutions of the Lotka Volterra. Prey Predator System Using Perturbation. Technique

UNSTEADY HEAT TRANSFER

Fourier Series and Parseval s Relation Çağatay Candan Dec. 22, 2013

16.512, Rocket Propulsion Prof. Manuel Martinez-Sanchez Lecture 3: Ideal Nozzle Fluid Mechanics

The Procedure Abstraction Part II: Symbol Tables and Activation Records

Revisiting what you have learned in Advanced Mathematical Analysis

Section 2: The Z-Transform

UNSTEADY STATE HEAT CONDUCTION

3.4 Repeated Roots; Reduction of Order

A modified hyperbolic secant distribution

Relation between Fourier Series and Transform

Inventory Management Model with Quadratic Demand, Variable Holding Cost with Salvage value

The model proposed by Vasicek in 1977 is a yield-based one-factor equilibrium model given by the dynamic

LINEAR 2 nd ORDER DIFFERENTIAL EQUATIONS WITH CONSTANT COEFFICIENTS

Single Correct Type. cos z + k, then the value of k equals. dx = 2 dz. (a) 1 (b) 0 (c)1 (d) 2 (code-v2t3paq10) l (c) ( l ) x.

5.1-The Initial-Value Problems For Ordinary Differential Equations

Math 266, Practice Midterm Exam 2

VIBRATION ANALYSIS OF CURVED SINGLE-WALLED CARBON NANOTUBES EMBEDDED IN AN ELASTIC MEDIUM BASED ON NONLOCAL ELASTICITY

Exact Solutions for Some Nonlinear Partial Differential Equations in Mathematical Physics

International Journal on Recent and Innovation Trends in Computing and Communication ISSN: Volume: 5 Issue:

EEE 303: Signals and Linear Systems

Engine Thrust. From momentum conservation

Reliability Analysis of a Bridge and Parallel Series Networks with Critical and Non- Critical Human Errors: A Block Diagram Approach.

Systems of First Order Linear Differential Equations

1 Finite Automata and Regular Expressions

Decline Curves. Exponential decline (constant fractional decline) Harmonic decline, and Hyperbolic decline.

Systems of First Order Linear Differential Equations

Laplace Transform. National Chiao Tung University Chun-Jen Tsai 10/19/2011

(A) 1 (B) 1 + (sin 1) (C) 1 (sin 1) (D) (sin 1) 1 (C) and g be the inverse of f. Then the value of g'(0) is. (C) a. dx (a > 0) is

Midterm. Answer Key. 1. Give a short explanation of the following terms.

Chapter 3. The Fourier Series

Lecture 4: Laplace Transforms

Inventory Model with Quadratic Demand under the Two Warehouse Management System

The Laplace Transform

A THREE COMPARTMENT MATHEMATICAL MODEL OF LIVER

A Study on the Nature of an Additive Outlier in ARMA(1,1) Models

Incremental DFT Based Search Algorithm for Similar Sequence

Fourier. Continuous time. Review. with period T, x t. Inverse Fourier F Transform. x t. Transform. j t

Jonathan Turner Exam 2-10/28/03

Boyce/DiPrima 9 th ed, Ch 7.8: Repeated Eigenvalues

Charging of capacitor through inductor and resistor

Equations and Boundary Value Problems

Week 06 Discussion Suppose a discrete random variable X has the following probability distribution: f ( 0 ) = 8

An Indian Journal FULL PAPER. Trade Science Inc. A stage-structured model of a single-species with density-dependent and birth pulses ABSTRACT

Speaker Identification using Spectrograms of Varying Frame Sizes

Ch 1.2: Solutions of Some Differential Equations

Stability and Optimal Harvesting of Modified Leslie-Gower Predator-Prey Model

Inverse Fourier Transform. Properties of Continuous time Fourier Transform. Review. Linearity. Reading Assignment Oppenheim Sec pp.289.

Lecture contents. Bloch theorem k-vector Brillouin zone Almost free-electron model Bands Effective mass Holes. NNSE 508 EM Lecture #9

Advanced Microeconomics II. Lijun Pan Nagoya University

CSE 245: Computer Aided Circuit Simulation and Verification

Boyce/DiPrima 9 th ed, Ch 2.1: Linear Equations; Method of Integrating Factors

INTEGRALS. Exercise 1. Let f : [a, b] R be bounded, and let P and Q be partitions of [a, b]. Prove that if P Q then U(P ) U(Q) and L(P ) L(Q).

FREE VIBRATION AND BENDING ANALYSES OF CANTILEVER MICROTUBULES BASED ON NONLOCAL CONTINUUM MODEL

Elementary Differential Equations and Boundary Value Problems

Chapter 5 The Laplace Transform. x(t) input y(t) output Dynamic System

Canadian Journal of Physics. Kantowski-Sachs modified holographic Ricci dark energy model in Saez-Ballester theory of gravitation

Chahrazed L Journal of Scientific and Engineering Research, 2018, 5(4): and

Chapter 4 Multifield Surface Bone Remodeling

An Optimal Ordering Policy for Inventory Model with. Non-Instantaneous Deteriorating Items and. Stock-Dependent Demand

A MATHEMATICAL MODEL FOR NATURAL COOLING OF A CUP OF TEA

A Production Inventory Model for Different Classes of Demands with Constant Production Rate Considering the Product s Shelf-Life Finite

CS 541 Algorithms and Programs. Exam 2 Solutions. Jonathan Turner 11/8/01

UNSTEADY FLOW OF A FLUID PARTICLE SUSPENSION BETWEEN TWO PARALLEL PLATES SUDDENLY SET IN MOTION WITH SAME SPEED

HIGHER ORDER DIFFERENTIAL EQUATIONS

On the Derivatives of Bessel and Modified Bessel Functions with Respect to the Order and the Argument

Advanced Engineering Mathematics, K.A. Stroud, Dexter J. Booth Engineering Mathematics, H.K. Dass Higher Engineering Mathematics, Dr. B.S.

Advanced Queueing Theory. M/G/1 Queueing Systems

2 T. or T. DSP First, 2/e. This Lecture: Lecture 7C Fourier Series Examples: Appendix C, Section C-2 Various Fourier Series

FOURIER ANALYSIS Signals and System Analysis

Lecture 11 Waves in Periodic Potentials Today: Questions you should be able to address after today s lecture:

Major: All Engineering Majors. Authors: Autar Kaw, Luke Snyder

Physics 160 Lecture 3. R. Johnson April 6, 2015

Control Systems. Modelling Physical Systems. Assoc.Prof. Haluk Görgün. Gears DC Motors. Lecture #5. Control Systems. 10 March 2013

More on FT. Lecture 10 4CT.5 3CT.3-5,7,8. BME 333 Biomedical Signals and Systems - J.Schesser

K x,y f x dx is called the integral transform of f(x). The function

ELECTRIC VELOCITY SERVO REGULATION

Lecture 2: Current in RC circuit D.K.Pandey

Bicomplex Version of Laplace Transform

Numerical Simulation for the 2-D Heat Equation with Derivative Boundary Conditions

Behaviors and Attitudes

PHA Second Exam. Fall On my honor, I have neither given nor received unauthorized aid in doing this assignment.

A Simple Method for Determining the Manoeuvring Indices K and T from Zigzag Trial Data

AR(1) Process. The first-order autoregressive process, AR(1) is. where e t is WN(0, σ 2 )

² Metres. Jack & Bore. Wesley Brooks Memorial Conservation Area (Fairy Lake) Directional Drilling** East Holland River. Tom Taylor Trail.

Lecture 21 : Graphene Bandstructure

On the Speed of Heat Wave. Mihály Makai

5. An object moving along an x-coordinate axis with its scale measured in meters has a velocity of 6t

Quality Improvement of Unbalanced Three-phase Voltages Rectification

Section 4.3 Logarithmic Functions

ANSWERS TO EVEN NUMBERED EXERCISES IN CHAPTER 2

Math 3301 Homework Set 6 Solutions 10 Points. = +. The guess for the particular P ( ) ( ) ( ) ( ) ( ) ( ) ( ) cos 2 t : 4D= 2

Stability of time-varying linear system

On the Existence and uniqueness for solution of system Fractional Differential Equations

Contraction Mapping Principle Approach to Differential Equations

Lecture 1: Numerical Integration The Trapezoidal and Simpson s Rule

1 Introduction to Modulo 7 Arithmetic

( ) ( ) + = ( ) + ( )

Available at Vol. 1, No. 1 (2006) pp

Library Support. Netlist Conditioning. Observe Point Assessment. Vector Generation/Simulation. Vector Compression. Vector Writing

Transcription:

Globl Journl of Pur nd Applid Mhmis. ISSN 0973-768 Volum 3, Numbr 207, pp. 8065-8074 Rsrh Indi Publiions hp://www.ripubliion.om Anlyil Soluion of A Diffrnil Equion h Prdis h Whr Condiion by Lornz Equions Using Homoopy Prurbion Mhod S. Muhukumr, C. Thngpndi 2, S. Mhlkshmi 2, M. Vrmuni 2 Rgisrr,Tmil Univrsiy,Thnjvur, Tmilndu, Indi. 2 Dprmn of Mhmis, Th Mdur Collg, Mduri, Tmilndu, Indi. *Corrsponding uhor C : Thngpndi, Absr Th Lornz quion hs md qulifying hos possibl whih hs inspird mny mhmiins o rsrh nd sudy hos [2]. Chos hory is h brnh of mhmis fousd on h bhviour of dynmil sysms h r highly snsiiv o iniil sysms. Choi bhviour xiss in mny nurl sysms suh s whr nd lim. Th drminisi nur of h sysm dos no mk bhviour prdibl. This bhviour is known s drminisi hos or simply. Approxim nlyil soluion of Lornz quion is obind by Homoopy prurbion mhod HPM.Furhrmor, in his work numril simulion of h problm is lso rpord using Silb/Mlb progrm. Agrmns bwn nlyil nd numril rsuls r nod. Th nlyil rsul rpord in his work is usful o undrsnd h bhvior of h sysm. Kywords: Lornz quions; Chos; Homoopy prurbion mhod HPM; Mhmil modling; S vribls.. INTRODUCTION Th Lornz sysm is sysm of ordinry diffrnil quions firs sudid by Edwrd Lornz. In priulr, h Lornz ror is s of hoi soluions of h Lornz sysm whih whn plod rsmbl burfly or figur igh. Smll diffrns in iniil ondiions yild widly divrging ouoms for suh dynmil Corrsponding uhor:c. Thngpndimhs@gmil.om

8066 S. Muhukumr l sysms populrly rfrrd o s h burfly ff- rndring long -rm prdiion of h bhvior impossibl in gnrl. Th ph h ld o Lornz o hs quions bgn wih n ffor o find simpl modl problm whih h mhods usd for sisil whr forsing would fil. Th Lornz quions r lso onnd o ohr physil phnomnon []. Th Lornz quion hs md qulifying hos possibl whih hs inspird mny mhmiins o rsrh nd sudy hos [2]. Chos hory is h brnh of mhmis fousd on h bhvior of dynmil sysms h r highly snsiiv o iniil sysms. Choi bhvior xiss in mny nurl sysms suh s whr nd lim. Th drminisi nur of h sysm dos no mk bhvior prdibl. This bhvior is known s drminisi hos or simply hos. This bhvior n b sudid hrough nlysis of hoi mhmil modl. Chos hory hs ppliions in svrl disiplins suh s morology, soiology, physis, ology, onomis, biology. In hoi sysms, h unriny in fors inrss xponnilly wih h lpsd im. Th modl h inrodud [3] n b hough of gross simplifiion of on fur of mosphr nmly h fluid moion drivn by h hrml buoyny known s onvion. Th modl dsribs h onvion moion of fluid in smll, idlizd Ryligh Bnrd ll. Curry [4] hs shown h if h mod runion is no don, bu insd suffiin mods r rind o giv numril onvrgn, h hos dispprs. On h ohr hnd MLuglin nd Mrin [5] hs showd h hos is obind for hr dimnsion vrsion. Th xprimnl sysm dsribd by Lornz quions is h Riik dynmo homopolr gnror wih h oupu fd bk hrough induors nd rsisors o h oil gnring h mgni fild [6]. Th oupld irui nd roion quions n b rdud o h Lornz form nd xprimns [7]. In Mrh 963, Lornz wro h h wnd o inrodu, ordinry diffrnil quions whos soluion xmpl h simpls xmpl of drminisi non priodi flow nd fini mpliud onvion. In his ppr h xmins h work of Brry Slzmn nd John Ryligh whil inorporing svrl physil phnomn [8,9]. Lornz usd hr vribls o onsru simpl modl bsd on h 2 dimnsionl rprsnion of rh s mosphr. Th purpos of his ppr is o driv h pproxim xprssions of s vribls of Lornz quions using Homoopy prurbion mhod for ll vlus of prmrs. 2. Mhmil modling nd nlysis L us onsidr h diffrnil quion s follows dx x y d

Anlyil Soluion of A Diffrnil Equion h Prdis h Whr. 8067 dy bx y xz 2 d dz d z xy 3 whr, b, nd r prmrs. Th iniil ondiions r A =0; x=, y=, z= 4 3. Anlyil soluion of Lornz quion using Homoopy Prurbion mhod Nonlinr sysm of quions plys n imporn rol in physis, hmisry nd biology. Consruing of priulr x soluion for hs quions rmins n imporn problm. In h ps mny uhors minly hd pid nion o sudy h soluion of nonlinr quions by using vrious mhods. Th Homoopy prurbion mhod hs bn workd ou ovr numbr of yrs by numrous uhors. Th Homoopy prurbion mhod HPM ws proposd by H nd ws sussfully pplid o uonomous ordinry diffrnil quions o nonlinr polyrysllin solids nd ohr filds. In his mhod h soluion produr is vry simpl nd only fw irions ld o high ur soluions whih r vlid for whol soluion domin. By solving h Eqns. - 4, using Homoopy Prurbion mhod w obin h nlyil soluions of Lornz quion s follows: x y - b b - 2 2 2 5 b 6

8068 S. Muhukumr l - 2 2 2 2 2 2 b z 7 4. DISCUSSION Eqns. 5-7 r h simpl nlyil xprssions of s vribl for ll vlus of prmrs, b nd. Figur, rprsns s vribl x vrsus im for fixd vlus of b=.2 nd =0.2. Figur 2, rprsns s vribl y vrsus im for fixd vlu of b=.5. Figur 3, rprsns s vribl z vrsus im for fixd vlus of =2 nd b = 0.0. 5. CONCLUSION Approxim nlyil soluions of h Lornz quions r prsnd using Homoopy Prurbion mhod. A simpl nd nw mhod of siming h s vribls r drivd. This soluion produr n b sily xndd o ll kinds of non-linr diffrnil quions wih vrious omplx boundry ondiions in nzym subsr rion diffusion prosss. Th xprssions providd in his work r usful o undrsnd h bhvior of h sysm. Rfrns [] Glik, J. Chos: Mking Nw Sin. Pnguin Books, Nw York, NY, 987. [2] S.H. Srogz, Nonlinr Dynmis nd hos. Addison Wsly, Rding, A, 994. [3] E. Lornz, Trnsions of h Nw York Admy of Sins, 409 963. [4] J.H. Curry, Commun. Mh. Phys. 60, 93 978. [5] J.B. MLughlin nd P.C. Mrin, Phys. Rv. A2, 86 975. [6] E.A. Jkson, Cmbridg Univrsiy Prss, Vol.2, hp. 7 99. [7] K.A. Robbins, Mh. Pro. Cmb. Phil. So. 82, 309 977. [8] Brdly, Lrry, Chos nd Frls. 200. [9] Viswnh, Diwkr. Th Frl Propry of h Lornz Aror. Physi D: Nonlinr Phnomn, Volum 90, Issus 2, Mrh 2004.

Anlyil Soluion of A Diffrnil Equion h Prdis h Whr. 8069 Appndix A: Bsi onps of h Homoopy Prurbion mhod To xplin his mhod l us onsidr h following funion Aw fr 0; r A. Wih h boundry ondiions of w Bw, 0 ; r A.2 n Whr A, B, f r nd Г r gnrl diffrnil opror, boundry opror, known nlyi funion nd h boundry of h domin Ω rspivly. Gnrlly spking h opror A n b dividd ino linr pr L nd nonlinr pr N. Eqn. A. n b wrin s Lw Nw fr 0 A.3 W onsru homoopy zr,p : Ω x [0,] R whih sisfis Hz,p p[ L w L z] p[ A z f r] 0, p [0,],r 0 or Hz,p Lw0 L z pl w0 p[ N z f r] 0 A.4 whr p ϵ [0,] is n mbdding prmr, whil w0 is n iniil pproximion of Eqn.A.,whih sisfis h boundry ondiions. Obviously, from Eqn.A.4 w will hv A.5 Hz,0 Lz Lw0 0 Hz, A z f r A.6 Th hnging pross of p from zro o uniy is jus h of zr, p from w0 o wr.in Topology, his is lld dformion, whil L z L w nd A z f r r lld 0 Homoopy. Aording o h HPM, w n firs us mbdding prmr p s smll prmr, nd ssum h h soluions of Eqns. A.3 nd A.4 n b wrin s powr sris in p. 2 z z pz p z... A.7 0 2 sing p= rsuls in h pproxim soluion of Eqn.A. w lim p z z0 z z2... A.8 Th ombinion of h Prurbion mhod nd h Homoopy mhod is lld h HPM, whih limins h drwbks of h rdiionl Prurbion mhods whil kping ll is dvngs.

8070 S. Muhukumr l Appndix B: Anlyil soluions for h dimnsionlss onnrions In his ppndix, w indi how Eqns. 5 7 in his ppr hs bn drivd. To find h soluions of Eqn., w onsru Homoopy s follows, from Eqn.A., dx dx p x p x y 0 B. d d Th iniil pproximions r s follows 2 x x0 px p x... B.2 2 0 dx0 P : x0 0 d B.3 dx P : x y0 0 d B.4 dx P 2 2 : x2 y 0 d B.5 Solving h bov Eqns. B. o B.5, w g x 0 B.6 x B.7 b x2 Adding h bov Eqns. B.6 o B.8, w g Eqn. 5 Similrly w n obin Eqns. 6 nd 7. Appndix C funion p4num opions= ods 'RlTol',-6,'Ss','on'; %iniil ondiions Xo= [; ; ]; spn = [0,]; i [,X] = od45@tsfunion,spn,xo,opions; o B.8

Anlyil Soluion of A Diffrnil Equion h Prdis h Whr. 807 figur hold on %plo, X:,,'-' plo, X:,2,'-' %plo, X:,3,'-' lgnd'x','x2','x3' ylbl'x' xlbl'' rurn funion [dx_d]= TsFunion,x =2;b=0.0;=; dx_d = -*x+*x2; dx_d2 =b*x-x2-x*x3; dx_d3 =-*x3+x*x2; dx_d = dx_d'; rurn Figur : Plo of S vribl x vrsus im vrious vlus of h prmrs. Solid lins rprsn numril soluions whrs h dod lin rprsns nlyil soluions.

8072 S. Muhukumr l Figur 2: Plo of S vribls x vrsus im vrious vlus of h prmrs. Solid lins rprsn numril soluions whrs h dod lin rprsns nlyil soluions. Figur 3: Plo of S vribls y vrsus im vrious vlus of h prmrs nd. Solid lins rprsn numril soluions whrs h dod lin rprsns nlyil soluions.

Anlyil Soluion of A Diffrnil Equion h Prdis h Whr. 8073 Figur 4: Plo of S vribls y vrsus im vrious vlus of h prmrs nd. Solid lins rprsn numril soluions whrs h dod lin rprsns nlyil soluions. Figur 5: Plo of S vribls z vrsus im vrious vlus of h prmrs nd b.

8074 S. Muhukumr l Solid lins rprsn numril soluions whrs h dod lin rprsns nlyil soluions. Figur 6: Plo of S vribls z vrsus im vrious vlus of h prmrs nd b. Solid lins rprsn numril soluions whrs h dod lin rprsns nlyil soluions.