Appendix B. The Finite Difference Scheme
|
|
- Magnus Edwards
- 5 years ago
- Views:
Transcription
1 140 APPENDIXES Appendx B. The Fnte Dfference Scheme In ths appendx we present numercal technques whch are used to approxmate solutons of system A comprehensve treatment of theoretcal and mplementaton ssues of dscretzaton methods for advecton-dffusonreacton problems are gven n the monograph by Hundsdorfer and Verwer Hundsdorfer, Verwer 003. Very nterestng applcatons of these results for bomedcal problems are descrbed by Gersch and Chaplan Gersch, Chaplan 006, Gersch and Verwer Gersch, Verwer 00 see also references gven n these publcatons. Here we propose to use a lnearzed mplct backward Euler method for the approxmaton of the dffuson-reacton subproblems and the explct forward Euler method for soluton of the advecton subproblem. We have restrcted to the frst order methods due to ther robust stablty. Note that our man goal s to nvestgate the nfluence of a possble ll-posedness of the PDE system to the asymptotcal behavour of the soluton. Appendx B.1. The method of lnes. Dscretzaton n space We use the method of lnes MOL approach see, Gersch, Chaplan 006; Hundsdorfer, Verwer 003. At the frst step we approxmate the spatal dervatves n the PDE by applyng robust and accurate approxmatons targeted for specal physcal processes descrbed by dfferental equatons. Doman [0, 1] s covered by a dscrete unform grd ω h = {x : x = h, = 0,..., N 1}, x N = 1 wth the grd ponts x. On the semdscrete doman ω h k [0, T ] we defne functons U t = Ux, t, V t = V x, t, = 0,..., N 1, here U, V approxmate solutons ux, t, vx, t on the dscrete grd ω h at tme moment t. We also defne the forward and backward space fnte dfferences wth respect to x: x U = U +1 U h, x U = U U 1. h Usng the fnte volume approach, we approxmate the dffuson and reac-
2 APPENDIXES 141 ton terms n the PDE system 3.1 by the followng fnte dfference equatons: A DR1 U, = D x x U + γru 1 U, U p A DR V, U, = x x V + γ 1 + βu p V. B.57 The stencl of the dscrete scheme requres to use functons defned outsde of the grd ω h. We apply perodcty boundary condtons 3.3 to defne dscrete functons for < 0 or N: U = U N, U N 1+ = U 1 for > 0. B.58 The advecton term n equaton 3.1 depends on the varable velocty ax, t := χ 1 + αv v x, therefore the maxmum prncple s not vald for the respectve transport equaton. But problem 3.1 stll has non-negatve solutons, and ths property can be preserved on the dscrete level by applyng proper upwndng approxmatons. The dscrete spatal approxmaton of the velocty s computed by a + 1 t = χ 1 + αv + V +1 / xv. In the followng we consder the upwnd-based dscrete fluxes Gersch, Chaplan 006; Hundsdorfer, Verwer 003: [ F T U, a, + 1/ = a + 1 U + ψθ ] U +1 U, a+ 1 0, B.59 F T U, a, + 1/ = a + 1 wth the Koren lmter functon ψθ = max [ U+1 + ψ1/θ +1 ] U U +1, a+ 1 0, mn 1, θ, θ. The lmter ψ depends on the smoothness montor functon θ = U U 1 U +1 U. < 0,
3 14 APPENDIXES For ψ = 0 we get the standard frst-order upwnd flux F T UW U, a, + 1/ = max a + 1, 0 U + mn a + 1, 0 U +1. Let us denote the dscrete advecton operator as A T U, V, = 1 h FT U, a, + 1/ F T U, a, 1/. Then we obtan a nonlnear ODE system for the evaluaton of the approxmate sem-dscrete solutons du = A T U, V, + A DR1 U,, x ω h, B.60 dv = A DRV, U,. Appendx B.. Operator splttng methods In order to develop effcent solvers n tme for the obtaned large ODE systems we apply the splttng technques. They take nto account the dfferent nature of the dscrete operators defnng the advecton A T U, V, and the dffuson-reacton A DR1 U,, A DR V, U, terms. The system resolvng the sem-dscrete advecton process can be solved very effcently by usng explct solvers, whle the dffuson-reacton sem-dscrete system s stff and t requres an mplct treatment. Also we are nterested n preservng at the dscrete level the postvty and/or boundedness of the soluton, f such propertes hold for the dfferental ODE system. Frst we consder the symmetrcal splttng method also known as the Strang splttng Strang Gven approxmatons U n, V n at tme t n, we
4 APPENDIXES 143 compute solutons at t n+1 = t n + by the followng scheme: du = A T U, V n,, U t n = U n, t n t t n+ 1 = t n + /, B.61 du du dv = A DRV, U,, = A DR1 U,, U t n = U n+ 1, t n t t n+1, B.6 = A T U, V n,, U t n = U n+1, t n+ 1 t t n+1 V t n = V n, t n t t n+1. B.63 Here we also have spltted the gven ODE system nto two blocks wth respect to U and V functons. Lemma 1. Solutons of the splttng ODE problem B.61 B.63 are nonnegatve f U n 0 and V n 0 for all x ω h. Proof. The proof for the advecton subsystems follows from the constructon of the dscrete fluxes by usng the upwndng technque. The proof for the dffuson-reacton subsystems follows from the lemma n Gersch, Verwer 00 that the soluton of an ntal value problem for systems of ODEs dy = F t, Y t, t 0, Y 0 = Y 0 s non-negatve f and only f for all t and any vector V R m and all 1 m v = 0, v 0 for all = f t, V 0. It s easy to see that for the dffuson-reacton subsystems the dffuson parts of the matrces are dagonally domnant and all off-dagonal entres are nonpostve. For the reacton parts the requred estmates are also trvally satsfed. Lemma. If 0 U n C for all x ω h, then a soluton of the splttng ODE problem B.6 s also bounded U maxc, 1. Proof. Here we use the fact that U = 1 s a stable attractor of the reacton functon. Let U n = max U n. Then t follows from the defnton of A DR1 U, that n the worst case du > 0, f C < 1, du = 0, f C = 1,
5 144 APPENDIXES and The lemma s proved. du < 0, f C > 1. Appendx B.3. Numercal ntegraton of ODEs There are many numercal ntegraton methods for soluton of non-stff and stff ODEs. For detaled dscussons of these schemes we refer the reader to Gersch, Chaplan 006; Harer, Norset, Wanner 1993; Harer, Wanner 1996; Hundsdorfer, Verwer 003. Let ω be a unform tme grd ω = {t n : t n = n, n = 0,..., M, M = T f }, here s the tme step. For smplcty ths step sze s taken constant. In the followng, we consder numercal approxmatons U n, V n to the exact soluton values ux, t n, vx, t n at the grd ponts x, t n ω h ω. Remark 1. In Baronas, Šmkus 011, the explct forward Euler scheme s used to solve problem 3.1. Snce no detals are gven n Baronas, Šmkus 011 on approxmatons of spatal dervatves, we use dscrete operators ntroduced n prevous sectons and wrte the explct forward Euler scheme as: U n+1 V n+1 U n V n = A T U n, V n, + A DR1 U n,, = A DR V n, U n,. We note that ths scheme can be wrtten as a splttng algorthm: U n+ 1 U n U n+1 U n+ 1 V n+1 V n = A T U n, V n,, = A DR1 U n,, = A DR V n, U n,.
6 APPENDIXES 145 Thus the explct Euler scheme can be consdered as a specal case of splttng algorthms. Despte easy mplementaton and good parallel propertes of explct algorthms, the man drawback of the explct Euler method s that due to the condtonal stablty we must restrct the ntegraton step to Ch for stff dscrete dffuson-reacton subproblems. The Rosenbrock and mplct Runge-Kutta methods are successfully appled for ntegraton of a stff part of the splttng semdscrete-scheme,.e. dffuson reacton equatons B.6, B.63, see Gersch, Chaplan 006; Gersch, Verwer 00; Hundsdorfer, Verwer 003. Here we propose to use a lnearzed mplct backward Euler method for the approxmaton of the dffuson-reacton subproblems and the explct forward Euler method for soluton of the advecton subproblem. We have restrcted to the frst order methods due to ther robust stablty. Note that our man goal s to nvestgate the nfluence of a possble ll-posedness of the PDE system to the asymptotcal behavour of the soluton. We dscretze the semdscrete problem B.61 B.63 wth the fully dscrete scheme U n+ 1 3 U n 0.5 U n+ 3 U n+ 1 3 U n+1 U n V n+1 V n = A T U n, V n,, B.64 = D x x U n+ 3 + γru n U n+ 3, B.65 = A T U n+ 3, V n,, B.66 = A DR V n+1, U n+1,. B.67 We apply two splttngs of the advecton term, because then we use only half of the splttng step sze for the explct method. Ths doubles the stablty and postvty domans of the explct method see Gersch, Verwer 00. Next we prove that statements of Lemma 1 and hold also for solutons of the fully dscrete fnte dfference scheme B.64 B.67 Lemma 3. For a suffcently small tme step 0 solutons of the fnte dfference scheme B.64 B.67 are non-negatve f U n 0 and V n 0 for all x ω h. Proof. The proof for the advecton problems B.64 and B.66 follows from
7 146 APPENDIXES the constructon of the dscrete fluxes by usng the upwndng technque and selecton of a suffcently small tme step 0. The proof for the dffuson-reacton problems B.65 and B.67 follows from the maxmum prncple Samarsk 001. For example, consder problem B.65. We assume, that U n+ 3 = mn U n <N We wrte the dscrete equaton B.65 for U n γru n+ 1 3 n an explct form U n+ 3 = 1 + γr U n D h U n U n+ 3 1 U n+ 3 Snce U n and U n+ 3 ±1 U n+ 3 we get that U n Lemma 4. If 0 U n+ 1 3 C for all x ω h, then a soluton of the fnte dfference scheme B.65 s also bounded Proof. U n+ 3 maxc, 1, x ω h. The proof s based on the maxmum prncple and a specal form of the dscrete reacton term. Frst, we consder the case C 1. Let U n+ 3 = max U n+ 3. Then t follows from B.65 that 1+γrU n+ 1 3 U n+ 3 1+γr U n+ 1 3 = U n+ 3 Next we consder the case C > 1. Then we get that 1 n γru 1 + γru n U n+ 3 The lemma s proved. 1 n γru 1 + γru n+ 1 3 = 1 + U n γru n+ 1 3 < C.
Lecture 21: Numerical methods for pricing American type derivatives
Lecture 21: Numercal methods for prcng Amercan type dervatves Xaoguang Wang STAT 598W Aprl 10th, 2014 (STAT 598W) Lecture 21 1 / 26 Outlne 1 Fnte Dfference Method Explct Method Penalty Method (STAT 598W)
More informationNON-CENTRAL 7-POINT FORMULA IN THE METHOD OF LINES FOR PARABOLIC AND BURGERS' EQUATIONS
IJRRAS 8 (3 September 011 www.arpapress.com/volumes/vol8issue3/ijrras_8_3_08.pdf NON-CENTRAL 7-POINT FORMULA IN THE METHOD OF LINES FOR PARABOLIC AND BURGERS' EQUATIONS H.O. Bakodah Dept. of Mathematc
More informationAdditional Codes using Finite Difference Method. 1 HJB Equation for Consumption-Saving Problem Without Uncertainty
Addtonal Codes usng Fnte Dfference Method Benamn Moll 1 HJB Equaton for Consumpton-Savng Problem Wthout Uncertanty Before consderng the case wth stochastc ncome n http://www.prnceton.edu/~moll/ HACTproect/HACT_Numercal_Appendx.pdf,
More informationModule 3 LOSSY IMAGE COMPRESSION SYSTEMS. Version 2 ECE IIT, Kharagpur
Module 3 LOSSY IMAGE COMPRESSION SYSTEMS Verson ECE IIT, Kharagpur Lesson 6 Theory of Quantzaton Verson ECE IIT, Kharagpur Instructonal Objectves At the end of ths lesson, the students should be able to:
More informationCHAPTER 5 NUMERICAL EVALUATION OF DYNAMIC RESPONSE
CHAPTER 5 NUMERICAL EVALUATION OF DYNAMIC RESPONSE Analytcal soluton s usually not possble when exctaton vares arbtrarly wth tme or f the system s nonlnear. Such problems can be solved by numercal tmesteppng
More informationConvexity preserving interpolation by splines of arbitrary degree
Computer Scence Journal of Moldova, vol.18, no.1(52), 2010 Convexty preservng nterpolaton by splnes of arbtrary degree Igor Verlan Abstract In the present paper an algorthm of C 2 nterpolaton of dscrete
More informationNumerical Heat and Mass Transfer
Master degree n Mechancal Engneerng Numercal Heat and Mass Transfer 06-Fnte-Dfference Method (One-dmensonal, steady state heat conducton) Fausto Arpno f.arpno@uncas.t Introducton Why we use models and
More informationLecture 5.8 Flux Vector Splitting
Lecture 5.8 Flux Vector Splttng 1 Flux Vector Splttng The vector E n (5.7.) can be rewrtten as E = AU (5.8.1) (wth A as gven n (5.7.4) or (5.7.6) ) whenever, the equaton of state s of the separable form
More informationHigh resolution entropy stable scheme for shallow water equations
Internatonal Symposum on Computers & Informatcs (ISCI 05) Hgh resoluton entropy stable scheme for shallow water equatons Xaohan Cheng,a, Yufeng Ne,b, Department of Appled Mathematcs, Northwestern Polytechncal
More informationChapter 5. Solution of System of Linear Equations. Module No. 6. Solution of Inconsistent and Ill Conditioned Systems
Numercal Analyss by Dr. Anta Pal Assstant Professor Department of Mathematcs Natonal Insttute of Technology Durgapur Durgapur-713209 emal: anta.bue@gmal.com 1 . Chapter 5 Soluton of System of Lnear Equatons
More informationTransfer Functions. Convenient representation of a linear, dynamic model. A transfer function (TF) relates one input and one output: ( ) system
Transfer Functons Convenent representaton of a lnear, dynamc model. A transfer functon (TF) relates one nput and one output: x t X s y t system Y s The followng termnology s used: x y nput output forcng
More informationFormal solvers of the RT equation
Formal solvers of the RT equaton Formal RT solvers Runge- Kutta (reference solver) Pskunov N.: 979, Master Thess Long characterstcs (Feautrer scheme) Cannon C.J.: 970, ApJ 6, 55 Short characterstcs (Hermtan
More information2 Finite difference basics
Numersche Methoden 1, WS 11/12 B.J.P. Kaus 2 Fnte dfference bascs Consder the one- The bascs of the fnte dfference method are best understood wth an example. dmensonal transent heat conducton equaton T
More informationLecture 12: Discrete Laplacian
Lecture 12: Dscrete Laplacan Scrbe: Tanye Lu Our goal s to come up wth a dscrete verson of Laplacan operator for trangulated surfaces, so that we can use t n practce to solve related problems We are mostly
More information4DVAR, according to the name, is a four-dimensional variational method.
4D-Varatonal Data Assmlaton (4D-Var) 4DVAR, accordng to the name, s a four-dmensonal varatonal method. 4D-Var s actually a drect generalzaton of 3D-Var to handle observatons that are dstrbuted n tme. The
More informationCME 302: NUMERICAL LINEAR ALGEBRA FALL 2005/06 LECTURE 13
CME 30: NUMERICAL LINEAR ALGEBRA FALL 005/06 LECTURE 13 GENE H GOLUB 1 Iteratve Methods Very large problems (naturally sparse, from applcatons): teratve methods Structured matrces (even sometmes dense,
More informationHandout # 6 (MEEN 617) Numerical Integration to Find Time Response of SDOF mechanical system. and write EOM (1) as two first-order Eqs.
Handout # 6 (MEEN 67) Numercal Integraton to Fnd Tme Response of SDOF mechancal system State Space Method The EOM for a lnear system s M X + DX + K X = F() t () t = X = X X = X = V wth ntal condtons, at
More information6.3.4 Modified Euler s method of integration
6.3.4 Modfed Euler s method of ntegraton Before dscussng the applcaton of Euler s method for solvng the swng equatons, let us frst revew the basc Euler s method of numercal ntegraton. Let the general from
More information2.29 Numerical Fluid Mechanics Fall 2011 Lecture 12
REVIEW Lecture 11: 2.29 Numercal Flud Mechancs Fall 2011 Lecture 12 End of (Lnear) Algebrac Systems Gradent Methods Krylov Subspace Methods Precondtonng of Ax=b FINITE DIFFERENCES Classfcaton of Partal
More informationNUMERICAL DIFFERENTIATION
NUMERICAL DIFFERENTIATION 1 Introducton Dfferentaton s a method to compute the rate at whch a dependent output y changes wth respect to the change n the ndependent nput x. Ths rate of change s called the
More informationA Hybrid Variational Iteration Method for Blasius Equation
Avalable at http://pvamu.edu/aam Appl. Appl. Math. ISSN: 1932-9466 Vol. 10, Issue 1 (June 2015), pp. 223-229 Applcatons and Appled Mathematcs: An Internatonal Journal (AAM) A Hybrd Varatonal Iteraton Method
More informationModified Mass Matrices and Positivity Preservation for Hyperbolic and Parabolic PDEs
COMMUNICATIONS IN NUMERICAL METHODS IN ENGINEERING Commun. Numer. Meth. Engng 2000; 00:6 Prepared usng cnmauth.cls [Verson: 2000/03/22 v.0] Modfed Mass Matrces and Postvty Preservaton for Hyperbolc and
More informationErrors for Linear Systems
Errors for Lnear Systems When we solve a lnear system Ax b we often do not know A and b exactly, but have only approxmatons  and ˆb avalable. Then the best thng we can do s to solve ˆx ˆb exactly whch
More informationLectures - Week 4 Matrix norms, Conditioning, Vector Spaces, Linear Independence, Spanning sets and Basis, Null space and Range of a Matrix
Lectures - Week 4 Matrx norms, Condtonng, Vector Spaces, Lnear Independence, Spannng sets and Bass, Null space and Range of a Matrx Matrx Norms Now we turn to assocatng a number to each matrx. We could
More informationNumerical Transient Heat Conduction Experiment
Numercal ransent Heat Conducton Experment OBJECIVE 1. o demonstrate the basc prncples of conducton heat transfer.. o show how the thermal conductvty of a sold can be measured. 3. o demonstrate the use
More informationLab session: numerical simulations of sponateous polarization
Lab sesson: numercal smulatons of sponateous polarzaton Emerc Boun & Vncent Calvez CNRS, ENS Lyon, France CIMPA, Hammamet, March 2012 Spontaneous cell polarzaton: the 1D case The Hawkns-Voturez model for
More informationReport on Image warping
Report on Image warpng Xuan Ne, Dec. 20, 2004 Ths document summarzed the algorthms of our mage warpng soluton for further study, and there s a detaled descrpton about the mplementaton of these algorthms.
More informationA new Approach for Solving Linear Ordinary Differential Equations
, ISSN 974-57X (Onlne), ISSN 974-5718 (Prnt), Vol. ; Issue No. 1; Year 14, Copyrght 13-14 by CESER PUBLICATIONS A new Approach for Solvng Lnear Ordnary Dfferental Equatons Fawz Abdelwahd Department of
More informationParametric fractional imputation for missing data analysis. Jae Kwang Kim Survey Working Group Seminar March 29, 2010
Parametrc fractonal mputaton for mssng data analyss Jae Kwang Km Survey Workng Group Semnar March 29, 2010 1 Outlne Introducton Proposed method Fractonal mputaton Approxmaton Varance estmaton Multple mputaton
More informationpage 2 2 dscretzaton mantans ths stablty under a sutable restrcton on the tme step. SSP tme dscretzaton methods were frst developed by Shu n [20] and
page 1 A Survey of Strong Stablty Preservng Hgh Order Tme Dscretzatons Ch-Wang Shu Λ 1 Introducton Numercal soluton for ordnary dfferental equatons (ODEs) s an establshed research area. There are many
More informationInductance Calculation for Conductors of Arbitrary Shape
CRYO/02/028 Aprl 5, 2002 Inductance Calculaton for Conductors of Arbtrary Shape L. Bottura Dstrbuton: Internal Summary In ths note we descrbe a method for the numercal calculaton of nductances among conductors
More informationNew Method for Solving Poisson Equation. on Irregular Domains
Appled Mathematcal Scences Vol. 6 01 no. 8 369 380 New Method for Solvng Posson Equaton on Irregular Domans J. Izadan and N. Karamooz Department of Mathematcs Facult of Scences Mashhad BranchIslamc Azad
More informationThe Finite Element Method: A Short Introduction
Te Fnte Element Metod: A Sort ntroducton Wat s FEM? Te Fnte Element Metod (FEM) ntroduced by engneers n late 50 s and 60 s s a numercal tecnque for solvng problems wc are descrbed by Ordnary Dfferental
More informationPositivity-preserving time discretizations for production-destruction equations. with applications to non-equilibrium flows.
Postvty-preservng tme dscretzatons for producton-destructon equatons wth applcatons to non-equlbrum flows Juntao Huang and Ch-Wang Shu Abstract In ths paper, we construct a famly of modfed Patankar Runge-Kutta
More informationEPR Paradox and the Physical Meaning of an Experiment in Quantum Mechanics. Vesselin C. Noninski
EPR Paradox and the Physcal Meanng of an Experment n Quantum Mechancs Vesseln C Nonnsk vesselnnonnsk@verzonnet Abstract It s shown that there s one purely determnstc outcome when measurement s made on
More informationCanonical transformations
Canoncal transformatons November 23, 2014 Recall that we have defned a symplectc transformaton to be any lnear transformaton M A B leavng the symplectc form nvarant, Ω AB M A CM B DΩ CD Coordnate transformatons,
More informationMATH 5630: Discrete Time-Space Model Hung Phan, UMass Lowell March 1, 2018
MATH 5630: Dscrete Tme-Space Model Hung Phan, UMass Lowell March, 08 Newton s Law of Coolng Consder the coolng of a well strred coffee so that the temperature does not depend on space Newton s law of collng
More informationKey Words: Hamiltonian systems, canonical integrators, symplectic integrators, Runge-Kutta-Nyström methods.
CANONICAL RUNGE-KUTTA-NYSTRÖM METHODS OF ORDERS 5 AND 6 DANIEL I. OKUNBOR AND ROBERT D. SKEEL DEPARTMENT OF COMPUTER SCIENCE UNIVERSITY OF ILLINOIS AT URBANA-CHAMPAIGN 304 W. SPRINGFIELD AVE. URBANA, ILLINOIS
More informationChapter 12. Ordinary Differential Equation Boundary Value (BV) Problems
Chapter. Ordnar Dfferental Equaton Boundar Value (BV) Problems In ths chapter we wll learn how to solve ODE boundar value problem. BV ODE s usuall gven wth x beng the ndependent space varable. p( x) q(
More informationApplication of B-Spline to Numerical Solution of a System of Singularly Perturbed Problems
Mathematca Aeterna, Vol. 1, 011, no. 06, 405 415 Applcaton of B-Splne to Numercal Soluton of a System of Sngularly Perturbed Problems Yogesh Gupta Department of Mathematcs Unted College of Engneerng &
More informationMMA and GCMMA two methods for nonlinear optimization
MMA and GCMMA two methods for nonlnear optmzaton Krster Svanberg Optmzaton and Systems Theory, KTH, Stockholm, Sweden. krlle@math.kth.se Ths note descrbes the algorthms used n the author s 2007 mplementatons
More informationConsistency & Convergence
/9/007 CHE 374 Computatonal Methods n Engneerng Ordnary Dfferental Equatons Consstency, Convergence, Stablty, Stffness and Adaptve and Implct Methods ODE s n MATLAB, etc Consstency & Convergence Consstency
More informationDUE: WEDS FEB 21ST 2018
HOMEWORK # 1: FINITE DIFFERENCES IN ONE DIMENSION DUE: WEDS FEB 21ST 2018 1. Theory Beam bendng s a classcal engneerng analyss. The tradtonal soluton technque makes smplfyng assumptons such as a constant
More information2.29 Numerical Fluid Mechanics
REVIEW Lecture 10: Sprng 2015 Lecture 11 Classfcaton of Partal Dfferental Equatons PDEs) and eamples wth fnte dfference dscretzatons Parabolc PDEs Ellptc PDEs Hyperbolc PDEs Error Types and Dscretzaton
More informationNumerical Solution of Ordinary Differential Equations
Numercal Methods (CENG 00) CHAPTER-VI Numercal Soluton of Ordnar Dfferental Equatons 6 Introducton Dfferental equatons are equatons composed of an unknown functon and ts dervatves The followng are examples
More informationAA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS
AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS 1 / 74 AA214B: NUMERICAL METHODS FOR COMPRESSIBLE FLOWS The Fnte Dfference Method These sldes are based on the recommended textbook: Culbert B. Laney. Computatonal
More informationE91: Dynamics. E91: Dynamics. Numerical Integration & State Space Representation
E91: Dnamcs Numercal Integraton & State Space Representaton The Algorthm In steps of Δ f ( new ) f ( old ) df ( d old ) Δ Numercal Integraton of ODEs d d f() h Δ Intal value problem: Gven the ntal state
More informationNUMERICAL METHODS FOR FIRST ORDER ODEs
1723 - COMPUTATIONAL METHODS FOR FLOW IN POROUS MEDIA Sprng 2009 NUMERICAL METHODS FOR FIRST ORDER ODEs Lus Cueto-Felgueroso 1 PROBLEM STATEMENT Consder a system of frst order ordnary dfferental equatons
More informationCHAPTER 14 GENERAL PERTURBATION THEORY
CHAPTER 4 GENERAL PERTURBATION THEORY 4 Introducton A partcle n orbt around a pont mass or a sphercally symmetrc mass dstrbuton s movng n a gravtatonal potental of the form GM / r In ths potental t moves
More informationSalmon: Lectures on partial differential equations. Consider the general linear, second-order PDE in the form. ,x 2
Salmon: Lectures on partal dfferental equatons 5. Classfcaton of second-order equatons There are general methods for classfyng hgher-order partal dfferental equatons. One s very general (applyng even to
More informationPART 8. Partial Differential Equations PDEs
he Islamc Unverst of Gaza Facult of Engneerng Cvl Engneerng Department Numercal Analss ECIV 3306 PAR 8 Partal Dfferental Equatons PDEs Chapter 9; Fnte Dfference: Ellptc Equatons Assocate Prof. Mazen Abualtaef
More informationModelli Clamfim Equazioni differenziali 22 settembre 2016
CLAMFIM Bologna Modell 1 @ Clamfm Equazon dfferenzal 22 settembre 2016 professor Danele Rtell danele.rtell@unbo.t 1/22? Ordnary Dfferental Equatons A dfferental equaton s an equaton that defnes a relatonshp
More informationInexact Newton Methods for Inverse Eigenvalue Problems
Inexact Newton Methods for Inverse Egenvalue Problems Zheng-jan Ba Abstract In ths paper, we survey some of the latest development n usng nexact Newton-lke methods for solvng nverse egenvalue problems.
More informationModule 1 : The equation of continuity. Lecture 1: Equation of Continuity
1 Module 1 : The equaton of contnuty Lecture 1: Equaton of Contnuty 2 Advanced Heat and Mass Transfer: Modules 1. THE EQUATION OF CONTINUITY : Lectures 1-6 () () () (v) (v) Overall Mass Balance Momentum
More informationLecture 2: Numerical Methods for Differentiations and Integrations
Numercal Smulaton of Space Plasmas (I [AP-4036] Lecture 2 by Lng-Hsao Lyu March, 2018 Lecture 2: Numercal Methods for Dfferentatons and Integratons As we have dscussed n Lecture 1 that numercal smulaton
More informationChapter 4 The Wave Equation
Chapter 4 The Wave Equaton Another classcal example of a hyperbolc PDE s a wave equaton. The wave equaton s a second-order lnear hyperbolc PDE that descrbes the propagaton of a varety of waves, such as
More informationGeneralized Linear Methods
Generalzed Lnear Methods 1 Introducton In the Ensemble Methods the general dea s that usng a combnaton of several weak learner one could make a better learner. More formally, assume that we have a set
More informationENTROPY CONSERVATIVE SCHEMES AND ADAPTIVE MESH SELECTION FOR HYPERBOLIC CONSERVATION LAWS
Journal of Hyperbolc Dfferental Equatons Vol. 7, No. 3 () 383 44 c World Scentfc Publshng Company DOI:.4/S989677 ENTROPY CONSERVATIVE SCHEMES AND ADAPTIVE MESH SELECTION FOR HYPERBOLIC CONSERVATION LAWS
More informationNumerical Solution of two dimensional coupled viscous Burgers Equation using the Modified Cubic B-Spline Differential Quadrature Method
umercal Soluton of two dmensonal coupled vscous Burgers Equaton usng the odfed Cubc B-Splne Dfferental Quadrature ethod H. S. Shukla 1, ohammad Tamsr 1*, Vneet K. Srvastava, Ja Kumar 3 1 Department of
More informationEEE 241: Linear Systems
EEE : Lnear Systems Summary #: Backpropagaton BACKPROPAGATION The perceptron rule as well as the Wdrow Hoff learnng were desgned to tran sngle layer networks. They suffer from the same dsadvantage: they
More informationHidden Markov Models & The Multivariate Gaussian (10/26/04)
CS281A/Stat241A: Statstcal Learnng Theory Hdden Markov Models & The Multvarate Gaussan (10/26/04) Lecturer: Mchael I. Jordan Scrbes: Jonathan W. Hu 1 Hdden Markov Models As a bref revew, hdden Markov models
More informationSingular Value Decomposition: Theory and Applications
Sngular Value Decomposton: Theory and Applcatons Danel Khashab Sprng 2015 Last Update: March 2, 2015 1 Introducton A = UDV where columns of U and V are orthonormal and matrx D s dagonal wth postve real
More informationCSci 6974 and ECSE 6966 Math. Tech. for Vision, Graphics and Robotics Lecture 21, April 17, 2006 Estimating A Plane Homography
CSc 6974 and ECSE 6966 Math. Tech. for Vson, Graphcs and Robotcs Lecture 21, Aprl 17, 2006 Estmatng A Plane Homography Overvew We contnue wth a dscusson of the major ssues, usng estmaton of plane projectve
More informationLOW BIAS INTEGRATED PATH ESTIMATORS. James M. Calvin
Proceedngs of the 007 Wnter Smulaton Conference S G Henderson, B Bller, M-H Hseh, J Shortle, J D Tew, and R R Barton, eds LOW BIAS INTEGRATED PATH ESTIMATORS James M Calvn Department of Computer Scence
More informationNote 10. Modeling and Simulation of Dynamic Systems
Lecture Notes of ME 475: Introducton to Mechatroncs Note 0 Modelng and Smulaton of Dynamc Systems Department of Mechancal Engneerng, Unversty Of Saskatchewan, 57 Campus Drve, Saskatoon, SK S7N 5A9, Canada
More informationPerron Vectors of an Irreducible Nonnegative Interval Matrix
Perron Vectors of an Irreducble Nonnegatve Interval Matrx Jr Rohn August 4 2005 Abstract As s well known an rreducble nonnegatve matrx possesses a unquely determned Perron vector. As the man result of
More informationMarkov Chain Monte Carlo (MCMC), Gibbs Sampling, Metropolis Algorithms, and Simulated Annealing Bioinformatics Course Supplement
Markov Chan Monte Carlo MCMC, Gbbs Samplng, Metropols Algorthms, and Smulated Annealng 2001 Bonformatcs Course Supplement SNU Bontellgence Lab http://bsnuackr/ Outlne! Markov Chan Monte Carlo MCMC! Metropols-Hastngs
More informationSEMI-LAGRANGIAN SCHEMES FOR LINEAR AND FULLY NON-LINEAR DIFFUSION EQUATIONS
SEMI-LAGRANGIAN SCHEMES FOR LINEAR AND FULLY NON-LINEAR DIFFUSION EQUATIONS KRISTIAN DEBRABANT AND ESPEN R. JAKOBSEN Abstract. For lnear and fully non-lnear dffuson equatons of Bellman-Isaacs type, we
More informationFTCS Solution to the Heat Equation
FTCS Soluton to the Heat Equaton ME 448/548 Notes Gerald Recktenwald Portland State Unversty Department of Mechancal Engneerng gerry@pdx.edu ME 448/548: FTCS Soluton to the Heat Equaton Overvew 1. Use
More informationThe equation of motion of a dynamical system is given by a set of differential equations. That is (1)
Dynamcal Systems Many engneerng and natural systems are dynamcal systems. For example a pendulum s a dynamcal system. State l The state of the dynamcal system specfes t condtons. For a pendulum n the absence
More information[7] R.S. Varga, Matrix Iterative Analysis, Prentice-Hall, Englewood Clis, New Jersey, (1962).
[7] R.S. Varga, Matrx Iteratve Analyss, Prentce-Hall, Englewood ls, New Jersey, (962). [8] J. Zhang, Multgrd soluton of the convecton-duson equaton wth large Reynolds number, n Prelmnary Proceedngs of
More informationAerodynamics. Finite Wings Lifting line theory Glauert s method
α ( y) l Γ( y) r ( y) V c( y) β b 4 V Glauert s method b ( y) + r dy dγ y y dy Soluton procedure that transforms the lftng lne ntegro-dfferental equaton nto a system of algebrac equatons - Restrcted to
More informationVector Norms. Chapter 7 Iterative Techniques in Matrix Algebra. Cauchy-Bunyakovsky-Schwarz Inequality for Sums. Distances. Convergence.
Vector Norms Chapter 7 Iteratve Technques n Matrx Algebra Per-Olof Persson persson@berkeley.edu Department of Mathematcs Unversty of Calforna, Berkeley Math 128B Numercal Analyss Defnton A vector norm
More informationGlobal Optimization of Truss. Structure Design INFORMS J. N. Hooker. Tallys Yunes. Slide 1
Slde 1 Global Optmzaton of Truss Structure Desgn J. N. Hooker Tallys Yunes INFORMS 2010 Truss Structure Desgn Select sze of each bar (possbly zero) to support the load whle mnmzng weght. Bar szes are dscrete.
More informationGlobal Sensitivity. Tuesday 20 th February, 2018
Global Senstvty Tuesday 2 th February, 28 ) Local Senstvty Most senstvty analyses [] are based on local estmates of senstvty, typcally by expandng the response n a Taylor seres about some specfc values
More informationA quote of the week (or camel of the week): There is no expedience to which a man will not go to avoid the labor of thinking. Thomas A.
A quote of the week (or camel of the week): here s no expedence to whch a man wll not go to avod the labor of thnkng. homas A. Edson Hess law. Algorthm S Select a reacton, possbly contanng specfc compounds
More informationBallot Paths Avoiding Depth Zero Patterns
Ballot Paths Avodng Depth Zero Patterns Henrch Nederhausen and Shaun Sullvan Florda Atlantc Unversty, Boca Raton, Florda nederha@fauedu, ssull21@fauedu 1 Introducton In a paper by Sapounaks, Tasoulas,
More informationWeek3, Chapter 4. Position and Displacement. Motion in Two Dimensions. Instantaneous Velocity. Average Velocity
Week3, Chapter 4 Moton n Two Dmensons Lecture Quz A partcle confned to moton along the x axs moves wth constant acceleraton from x =.0 m to x = 8.0 m durng a 1-s tme nterval. The velocty of the partcle
More informationModelli Clamfim Equazione del Calore Lezione ottobre 2014
CLAMFIM Bologna Modell 1 @ Clamfm Equazone del Calore Lezone 17 15 ottobre 2014 professor Danele Rtell danele.rtell@unbo.t 1/24? Convoluton The convoluton of two functons g(t) and f(t) s the functon (g
More informationLecture 3: Shannon s Theorem
CSE 533: Error-Correctng Codes (Autumn 006 Lecture 3: Shannon s Theorem October 9, 006 Lecturer: Venkatesan Guruswam Scrbe: Wdad Machmouch 1 Communcaton Model The communcaton model we are usng conssts
More informationHongyi Miao, College of Science, Nanjing Forestry University, Nanjing ,China. (Received 20 June 2013, accepted 11 March 2014) I)ϕ (k)
ISSN 1749-3889 (prnt), 1749-3897 (onlne) Internatonal Journal of Nonlnear Scence Vol.17(2014) No.2,pp.188-192 Modfed Block Jacob-Davdson Method for Solvng Large Sparse Egenproblems Hongy Mao, College of
More informationAppendix for Causal Interaction in Factorial Experiments: Application to Conjoint Analysis
A Appendx for Causal Interacton n Factoral Experments: Applcaton to Conjont Analyss Mathematcal Appendx: Proofs of Theorems A. Lemmas Below, we descrbe all the lemmas, whch are used to prove the man theorems
More informationarxiv: v1 [math.ho] 18 May 2008
Recurrence Formulas for Fbonacc Sums Adlson J. V. Brandão, João L. Martns 2 arxv:0805.2707v [math.ho] 8 May 2008 Abstract. In ths artcle we present a new recurrence formula for a fnte sum nvolvng the Fbonacc
More informationProblem adapted reduced models based on Reaction-Diffusion Manifolds (REDIMs)
Problem adapted reduced models based on Reacton-Dffuson Manfolds (REDIMs) V Bykov, U Maas Thrty-Second Internatonal Symposum on ombuston, Montreal, anada, 3-8 August, 8 Problem Statement: Smulaton of reactng
More informationLecture 10 Support Vector Machines II
Lecture 10 Support Vector Machnes II 22 February 2016 Taylor B. Arnold Yale Statstcs STAT 365/665 1/28 Notes: Problem 3 s posted and due ths upcomng Frday There was an early bug n the fake-test data; fxed
More informationANSWERS. Problem 1. and the moment generating function (mgf) by. defined for any real t. Use this to show that E( U) var( U)
Econ 413 Exam 13 H ANSWERS Settet er nndelt 9 deloppgaver, A,B,C, som alle anbefales å telle lkt for å gøre det ltt lettere å stå. Svar er gtt . Unfortunately, there s a prntng error n the hnt of
More informationSome modelling aspects for the Matlab implementation of MMA
Some modellng aspects for the Matlab mplementaton of MMA Krster Svanberg krlle@math.kth.se Optmzaton and Systems Theory Department of Mathematcs KTH, SE 10044 Stockholm September 2004 1. Consdered optmzaton
More informationDifference Equations
Dfference Equatons c Jan Vrbk 1 Bascs Suppose a sequence of numbers, say a 0,a 1,a,a 3,... s defned by a certan general relatonshp between, say, three consecutve values of the sequence, e.g. a + +3a +1
More informationELASTIC WAVE PROPAGATION IN A CONTINUOUS MEDIUM
ELASTIC WAVE PROPAGATION IN A CONTINUOUS MEDIUM An elastc wave s a deformaton of the body that travels throughout the body n all drectons. We can examne the deformaton over a perod of tme by fxng our look
More informationDigital Signal Processing
Dgtal Sgnal Processng Dscrete-tme System Analyss Manar Mohasen Offce: F8 Emal: manar.subh@ut.ac.r School of IT Engneerng Revew of Precedent Class Contnuous Sgnal The value of the sgnal s avalable over
More information6) Derivatives, gradients and Hessian matrices
30C00300 Mathematcal Methods for Economsts (6 cr) 6) Dervatves, gradents and Hessan matrces Smon & Blume chapters: 14, 15 Sldes by: Tmo Kuosmanen 1 Outlne Defnton of dervatve functon Dervatve notatons
More informationChapter Newton s Method
Chapter 9. Newton s Method After readng ths chapter, you should be able to:. Understand how Newton s method s dfferent from the Golden Secton Search method. Understand how Newton s method works 3. Solve
More information8 Derivation of Network Rate Equations from Single- Cell Conductance Equations
Physcs 178/278 - Davd Klenfeld - Wnter 2015 8 Dervaton of Network Rate Equatons from Sngle- Cell Conductance Equatons We consder a network of many neurons, each of whch obeys a set of conductancebased,
More information1 GSW Iterative Techniques for y = Ax
1 for y = A I m gong to cheat here. here are a lot of teratve technques that can be used to solve the general case of a set of smultaneous equatons (wrtten n the matr form as y = A), but ths chapter sn
More informationSuppose that there s a measured wndow of data fff k () ; :::; ff k g of a sze w, measured dscretely wth varable dscretzaton step. It s convenent to pl
RECURSIVE SPLINE INTERPOLATION METHOD FOR REAL TIME ENGINE CONTROL APPLICATIONS A. Stotsky Volvo Car Corporaton Engne Desgn and Development Dept. 97542, HA1N, SE- 405 31 Gothenburg Sweden. Emal: astotsky@volvocars.com
More informationDiscretization. Consistency. Exact Solution Convergence x, t --> 0. Figure 5.1: Relation between consistency, stability, and convergence.
Chapter 5 Theory The numercal smulaton of PDE s requres careful consderatons of propertes of the approxmate soluton. A necessary condton for the scheme used to model a physcal problem s the consstency
More informationChp 3: Scalar Advection and Linear Hyperbolic Systems. By Prof. Dinshaw S. Balsara
Chp 3: Scalar Advecton and Lnear Hyperbolc Systems By Prof. Dnshaw S. Balsara 3.) Introducton We have seen the need for consstency and stablty n FDAs of PDEs. However, for the advecton equaton, whch s
More informationUsing T.O.M to Estimate Parameter of distributions that have not Single Exponential Family
IOSR Journal of Mathematcs IOSR-JM) ISSN: 2278-5728. Volume 3, Issue 3 Sep-Oct. 202), PP 44-48 www.osrjournals.org Usng T.O.M to Estmate Parameter of dstrbutons that have not Sngle Exponental Famly Jubran
More informationPhysics 181. Particle Systems
Physcs 181 Partcle Systems Overvew In these notes we dscuss the varables approprate to the descrpton of systems of partcles, ther defntons, ther relatons, and ther conservatons laws. We consder a system
More informationThe Feynman path integral
The Feynman path ntegral Aprl 3, 205 Hesenberg and Schrödnger pctures The Schrödnger wave functon places the tme dependence of a physcal system n the state, ψ, t, where the state s a vector n Hlbert space
More information