SYNChronised numerical methods
|
|
- Arthur George
- 6 years ago
- Views:
Transcription
1 1 SYNChronsed numercal methods By Alan S Dawes AWE, Computatonal Physcs Group, Aldermaston, Berkshre, RG7 4PR, UK Abstract The dscretzaton of partal dfferental equatons wll produce naccuraces, and n partcular symmetry errors. Tradtonally, symmetry s ftted nto the numercal method[7] based on three approaches. The frst s to algn the mesh wth the underlyng symmetry and adapt the numercal method accordngly[4, 14]. A second approach s to use the natural coordnate frame and to use a coordnate transformaton from physcal to computatonal space. For ths type of strategy the partcular symmetry s contaned wthn the Jacoban for the transformaton and s embedded n a Geometrc Conservaton Law[17]. Fnally, a thrd approach s to modfy the numercal method to reduce the asymmetrc errors, nherent wthn the truncaton terms, by usng more nformaton from the wder computatonal mesh stencl[13]. Partal dfferental equatons reman unchanged when certan transformatons are performed; these symmetres leave the equatons nvarant after a gven transformaton. Le group theory has been shown to be a general way of determnng these symmetres[15, 3, 16, 12] and has been used to fnd analytc solutons to many partal dfferental equatons, and n partcular, hydrodynamc related equatons[5, 2]. However, n general the dscretzed equatons do not preserve these symmetres. By extendng Le groups to the dscrete system[8, 9, 6, 10, 11] we wll develop gudng prncples for SYNChronsed numercal methods where the SYmmetres of the partal dfferental equatons are Numercally Captured. By usng ths framework we wll show, through the concept of numercal symmetry breakng, why some numercal technques are more sutable than others for a gven applcaton. c Brtsh Crown Copyrght MOD/23 1. Fundamentals Le group theory s the marrage between co-ordnate geometry and group theory. Transformatons form a group and are contnuous based upon the varaton of a group parameter. Le group theory s based upon transformatons of both the ndependent and dependent varables. For m ndependent varables (x 1,...,x m ) and n dependent varables (u 1,...,u n ) these defne a m+n dmensonal phase space n whch all ponts are treated as ndependent. They are transformed n a contnuous manner based on the group parameter and the arbtrarly defned pontwse functons thus, x k = F k (x 1,...,x m,u 1,...,u n ;), ũ j = F j (x 1,...,x m,u 1,...,u n ;), (1.1) for 1 k m and m+1 j n+m and where the tlde represents the transformed varable. It s assumed that all of the functons are contnuous and contnuously dfferentable (C ) n all of the varables and. Note, there s an dentfy value = o that leaves the co-ordnates unchanged. Whence, x k = F k (x 1,...,x m,u 1,...,u n ; o ), u j = F j (x 1,...,x m,u 1,...,u n ; o ). (1.2) Mathematcs n Defence 23
2 SYNChronsed numercal methods 2 If we consder an nfntesmal transformatons from o o + then they can be wrtten as, x k = x k + ξ k (x 1,,u 1, ) +O( 2 ), ũ j = u j + η k (x 1,,u 1, ) +O( 2 ), (1.3) where ξ k and η j are respectvely the nfntesmal generators for the k-th ndependent and j-th dependent varables. The generator for the group can now be wrtten, X = ξ ξ m + η η n. (1.4) x 1 x m u 1 u n The computatonal scentst wll be famlar wth the use of transformatons mappng from physcal to computatonal space[17]. However, for these types of transformatons they are geometrc n nature. Furthermore, they assume that both the tme t and dependent varables (u 1,...,u n ) reman unchanged after the transformaton. When compared to the transformatons Eq.(1.1) they represent a subset of the Le group transformatons. For the general partal dfferental equaton, G(x k,u,u,m,u,mn, )=0, where u,m = u x m, u,mn = 2 u x m x n, (1.5) the transformed PDE s related to the orgnal PDE by the transformaton thus, G( x k,ũ,ũ, m,ũ, mn, ;)=G(x k,u,u,m,u,mn, ; o )+ (X G)+O( 2 ), (1.6) where 0 s the dentfy value that leaves the PDE unchanged, the tlde symbol represents the transformed PDE and the generator X has been prolonged[15] to nclude transformatons of the partal dervatves found wthn G( ) = 0; see next secton. For nvarance of the PDE X G( )=H( )G( ) for the arbtrary functonh( ) provded G( )=0. Whence, G( x k,ũ,ũ, m,ũ, mn, ;)=G(x k,u,u,m,u,mn, ; o ). (1.7) 2. Extenson to Numercal Methods For a gven numercal method (NM) expressed as, G NM (x n,tn,τn,h n +,hn,un 1,un,un +1,un+1, )=0, (2.1) t can be wrtten after the Le group transformaton as, G NM ( ;)=G NM ( ; o )+(X NM G NM ) +O( 2 ), (2.2) wherex NM s the group operator whch has been prolonged to nclude the transformatons of the numercal method. Ths operator must now nclude terms that represent the dscrete dependent varables and ndependent quanttes (mesh postons) wth respect to. For example, Fgure 1 llustrates one way that the ponts(x n,tn,un ) can vary after a transformaton. The group operator X NM wll leave the PDE nvarant and for numercal nvarance, X NM G NM ( )=H NM ( )G NM ( ), (2.3) for the arbtrary dscrete functonh NM (x n,tn,τn,h n +,hn,un 1,un,un +1,un+1, ) snce G NM ( )= 0. Whence, G NM ( ;)=G NM ( ; o ). (2.4) When ths s satsfed the numercal method wll be SYNChronsed wth the PDEs. Alternatvely, a numercal method satsfyng ths wll be known as a SYNC scheme; SYmmetres of the partal dfferental equatons are Numercally Captured. Mathematcs n Defence 23
3 INVARIANT TRUNCATION ERROR 3 t ( x n+1, t n+1,ũ n+1 ) (x n+1,t n+1,u n+1 ) (x+1 n+1,tn+1 +1,un+1 +1 ) ( x (x n,tn,un ) +1 n+1, tn+1 +1,ũn+1 +1 ) (x n ,tn +1,un +1 ) x ( x n, tn,ũn ) ( x n +1, tn +1,ũn +1 ) 3. Invarant Truncaton Error u FIGURE 1. Relatonshp between dscrete Le co-ordnate ponts. Usng Taylor s theorem a numercal method can be shown to be equal to a modfed equaton (ME) n whch the orgnal homogeneous PDE now equals a truncaton error (TE) source term thus, G ME ( )=G( )=TE( ). (3.1) By applyng the prolonged group operator to the modfed equaton t wll be nvarant, lke the PDE, when the truncaton error s also nvarant thus, X NM G ME = δ(te) δ = 0. (3.2) Ths s the condton for a SYNChronsed numercal method. 4. Numercal Symmetry Breakng The equatons of hydrodynamcs can be wrtten n conservaton law form as, u t + f(u) = 0, (4.1) x where the dependent varable u represents some conserved quantty, such as mass. Ths equaton can be wrtten numercally n explct Euleran conservaton law form usng a fxed mesh (please refer to the books by Hrsch[1] for general descrptons of types of numercal methods) as, G NM ( ; o )=u n+1 u n + τn h ( F(un p,...,un +q ;h,τn ) F(u n 1 p,...,un 1+q ;h,τn ))=0, (4.2) Mathematcs n Defence 23
4 SYNChronsed numercal methods 4 where h s the grd spacng (assumed fxed and unform), τ n the tme-step, and F( ) are the numercal flux functons. For smplcty we shall use the notaton F+1/2 n = F(un p,...,un +q ;h,τn ). The ntegers p and q defne the sze of the mesh stencl p+q+1 for the flux functons. The condton that ths wll be a SYNC scheme sx NM G NM =H NM G NM. That s, δu n+1 ( δ δun 1 δ + δτ n ) ( h δ τn δh δf h 2 (F n δ +1/2 F 1/2 n )+ τn n ) +1/2 δfn 1/2 =H NM G NM, (4.3) h δ δ We consder ths framework for the non-lnear vscous Burgers equaton, u t + u u x = u ν 2 x 2, (4.4) where ν s a parameter representng vscosty. The group structure for ths equaton s gven by Hydon, and conssts of fve groups representng scale nvarance, projectve nvarance, Gallean nvarance, temporal and spatal translaton. Wthout loss of generalty we shall consder a three pont stencl (p = 0 and q = 1) and by prolongng the fve group operators yelds, XNM 1 = xn x n + 2t n t u n n u n + h h + 2τn τ u n n 1 u n u n +1 1 u n +1 u n+1 XNM 2 = xn tn +(t x n ) 2 +(x n t n n un tn ) + ht u n n h +(tn +t n+1 )τ n + τ n u n+1 XNM 3 = tn x n + u n + u n + 1 u n +[1 τ n S + t (u x )] +1 u n+1 (4.5) X 4 NM = t n X 5 NM = x n where S t + s a shft operator. We do not provde the whole structure forxnm 2 due to ts complexty. From XNM 4 G NM andxnm 5 G NM t s mmedately obvous thath NM = 0 and that G NM cannot be a functon of x n and t n. ForXNM 1 G NM t can be shown (usng the method of characterstcs) that the flux can be wrtten n terms of dscrete nvarants, J 4 =G(J 1,J 2,J 3 ), (4.6) whereg s an arbtrary flux functon wth the dscrete nvarants gven by, J 1 = un +1 u n, J 2 = h2 τ n, J 3 = u n h, J 4 = τ n F n +1/2. (4.7) In practce the values of the dscrete nvarants wll be restrcted to satsfy numercal stablty and monotoncty. It can be shown that the nvarant flux functon Eq.(4.6) does not satsfy the Gallean and Projectve groups. Therefore, a numercal scheme of the type defned by Eq.(4.2) wll break the symmetry nherent n the PDE. The lack of nvarance les n the ntroducton of gradents wthn some of the prolonged group operators. The Gallean operator ncludes the gradent term [1 τ n S t + (u x )] and represents the relatve moton between the fxed and movng frames and s llustrated n Fgure 2. At t = t n both frames are concdent, but after a short perod τ n the Gallean frame has moved relatve to the fxed frame. Wthn the fxed frame x n = x n+1 and wthn the movng frame x n = x n+1, but the relatve postons between the dfferent frames changes. The value of ũ n+1 moves along a trajectory relatve to the fxed frame and ts value can be thought of as beng nterpolated from the values wthn the fxed frame. Mathematcs n Defence 23
5 CONCLUSIONS 5 x n+1 τ n x n+1 t n+1 = t n + τ n x n h xn x n + h t = t n FIGURE 2. Relatonshp between fxed (black) and movng frames (red). Ths s the source of the gradent term. Note, ths type of analyss can be appled to the projectve operator. A numercal scheme based upon a fxed stencl can break certan symmetres. A scheme that s SYNChronsed wth all of the PDEs symmetres s one where the mesh moves to ts natural postons. In ths paper we wll not consder how ths can be done due to space lmtatons. However, the reader should consder numercal methods based on dscrete nvarants[19, 10, 11] and those based on Movng Frames[18, 8]. 5. Conclusons Usng Le group theory we have presented a concse mathematcal defnton for the nvarance of partal dfferental equatons under a gven transformaton. In ths paper ths framework was extended to ts representatve numercal method and t was shown that t can reman nvarant only f ts truncaton error remans nvarant. Numercal schemes of ths type are SYNChronsed wth ther respectve PDEs symmetres. In ths paper we have shown how numercal nvarance errors are due to numercal symmetry breakng. A scheme that s SYNChronsed s one where the mesh moves to ts natural postons. REFERENCES [1] C. Hrsch. Numercal Computaton of Internal and External Flows, Volumes 1 and 2. John Wley & Sons. [2] S. Atzen and J. Meyer-Ter-Vehn. The Physcs of Inertal Fuson. Clarendon Press, Oxford, [3] G.W. Bluman and J.D. Cole. Smlarty Methods for Dfferental Equatons. Sprnger-Verlag, New York, [4] E. Caramana and P. Whalen. Numercal Preservaton of Symmetry Propertes of Contnuum Problems. J. Comp. Phys. 141, pages , [5] S.V. Coggeshall and R.A. Axford. Le Group Invarance Propertes of Radaton Hydrodynamcs Equatons and Ther Assocated Smlarty Solutons. J. Phys. Fluds, pages , [6] A.S. Dawes. Invarant Numercal Methods. ICFD Conference on Numercal Methods n Fluds, Readng Unversty, UK, March, [7] A.S. Dawes. Sources of Cartesan Mesh Induced Asymmetres Based Upon the Lagrangan + Remap Method. Numercal Methods for Mult-Materal Flud Flows, Praque, CZ, September 2007, [8] A.S. Dawes. Invarant Numercal Methods. Int. J. Num. Meth. Fluds, 56, pages , [9] A.S. Dawes. Invarant Numercal Methods. Symmetry Preservng Methods n CFD and Beyond, ICFD Workshop, Unversty of Readng, Mathematcs n Defence 23
6 SYNChronsed numercal methods 6 [10] V.A. Dorodntsyn. Fnte dfference methods entrely nhertng the symmetry of the orgnal equatons. n Modern Group Analyss: Advanced analytcal and computatonal methods n mathematcal physcs, edted by N. Ibragmov (CRC Press, New York), [11] V.A. Dorodntsyn. Symmetry of Fnte Dfference Equatons. n CRC Handbook of Le Group Analyss of Dfferental Equatons, Volume 1: Symmetres, Exact Solutons and Conservaton Laws, edted by N. Ibragmov (CRC Press, New York), [12] P.E. Hydon. Symmetry Methods for Dfferental Equatons. Cambrdge Unversty Press, [13] A. Kumar. Isotropc Fnte-Dfference. J. Comp. Phys. 2, pages , [14] Margoln and M. Shashkov. Usng a Curvlnear Grd to Construct Symmetry Preservng Dscretsatons for Lagrangan Gas Dynamcs. J. Comp. Phys. 149, pages , [15] P. Olver. Applcaton of Le Groups to Dfferental Equatons. Sprnger, [16] H. Stephan. Dfferental Equatons. Ther soluton usng symmetres. Cambrdge Unversty Press, [17] P.D. Thomas and C.K. Lombard. Geometrc Conservaton Law and ts Applcatons to Flow Computatons on Movng Grds. AIAA Journal, Vol. 17, No. 10, pages , [18] M.Chhay et al. Comparson of some Le-symmetry-based ntegrators. J. Comp. Phys., Volume 230 Issue 5, pages , 21. [19] Bakrova and Dorodntsyn and Kozlov. Symmetry preservng dscrete schemes for the heat transfer equatons. J. Phys. A: Math. Gen., 30, pages , 1997 Mathematcs n Defence 23
χ x B E (c) Figure 2.1.1: (a) a material particle in a body, (b) a place in space, (c) a configuration of the body
Secton.. Moton.. The Materal Body and Moton hyscal materals n the real world are modeled usng an abstract mathematcal entty called a body. Ths body conssts of an nfnte number of materal partcles. Shown
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 informationKinematics of Fluids. Lecture 16. (Refer the text book CONTINUUM MECHANICS by GEORGE E. MASE, Schaum s Outlines) 17/02/2017
17/0/017 Lecture 16 (Refer the text boo CONTINUUM MECHANICS by GEORGE E. MASE, Schaum s Outlnes) Knematcs of Fluds Last class, we started dscussng about the nematcs of fluds. Recall the Lagrangan and Euleran
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 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 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 informationMathematical Preparations
1 Introducton Mathematcal Preparatons The theory of relatvty was developed to explan experments whch studed the propagaton of electromagnetc radaton n movng coordnate systems. Wthn expermental error the
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 informationAppendix B. The Finite Difference Scheme
140 APPENDIXES Appendx B. The Fnte Dfference Scheme In ths appendx we present numercal technques whch are used to approxmate solutons of system 3.1 3.3. A comprehensve treatment of theoretcal and mplementaton
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 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 informationGroup Analysis of Ordinary Differential Equations of the Order n>2
Symmetry n Nonlnear Mathematcal Physcs 997, V., 64 7. Group Analyss of Ordnary Dfferental Equatons of the Order n> L.M. BERKOVICH and S.Y. POPOV Samara State Unversty, 4430, Samara, Russa E-mal: berk@nfo.ssu.samara.ru
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 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 informationPhysics 5153 Classical Mechanics. Principle of Virtual Work-1
P. Guterrez 1 Introducton Physcs 5153 Classcal Mechancs Prncple of Vrtual Work The frst varatonal prncple we encounter n mechancs s the prncple of vrtual work. It establshes the equlbrum condton of a mechancal
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 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 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 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 informationA Solution of the Harry-Dym Equation Using Lattice-Boltzmannn and a Solitary Wave Methods
Appled Mathematcal Scences, Vol. 11, 2017, no. 52, 2579-2586 HIKARI Ltd, www.m-hkar.com https://do.org/10.12988/ams.2017.79280 A Soluton of the Harry-Dym Equaton Usng Lattce-Boltzmannn and a Soltary Wave
More informationSymmetric Lie Groups and Conservation Laws in Physics
Symmetrc Le Groups and Conservaton Laws n Physcs Audrey Kvam May 1, 1 Abstract Ths paper eamnes how conservaton laws n physcs can be found from analyzng the symmetrc Le groups of certan physcal systems.
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 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 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 informationLecture 13 APPROXIMATION OF SECOMD ORDER DERIVATIVES
COMPUTATIONAL FLUID DYNAMICS: FDM: Appromaton of Second Order Dervatves Lecture APPROXIMATION OF SECOMD ORDER DERIVATIVES. APPROXIMATION OF SECOND ORDER DERIVATIVES Second order dervatves appear n dffusve
More informationPhysics 5153 Classical Mechanics. D Alembert s Principle and The Lagrangian-1
P. Guterrez Physcs 5153 Classcal Mechancs D Alembert s Prncple and The Lagrangan 1 Introducton The prncple of vrtual work provdes a method of solvng problems of statc equlbrum wthout havng to consder the
More informationn α j x j = 0 j=1 has a nontrivial solution. Here A is the n k matrix whose jth column is the vector for all t j=0
MODULE 2 Topcs: Lnear ndependence, bass and dmenson We have seen that f n a set of vectors one vector s a lnear combnaton of the remanng vectors n the set then the span of the set s unchanged f that vector
More informationComparative Studies of Law of Conservation of Energy. and Law Clusters of Conservation of Generalized Energy
Comparatve Studes of Law of Conservaton of Energy and Law Clusters of Conservaton of Generalzed Energy No.3 of Comparatve Physcs Seres Papers Fu Yuhua (CNOOC Research Insttute, E-mal:fuyh1945@sna.com)
More informationON A DETERMINATION OF THE INITIAL FUNCTIONS FROM THE OBSERVED VALUES OF THE BOUNDARY FUNCTIONS FOR THE SECOND-ORDER HYPERBOLIC EQUATION
Advanced Mathematcal Models & Applcatons Vol.3, No.3, 2018, pp.215-222 ON A DETERMINATION OF THE INITIAL FUNCTIONS FROM THE OBSERVED VALUES OF THE BOUNDARY FUNCTIONS FOR THE SECOND-ORDER HYPERBOLIC EUATION
More informationMechanics Physics 151
Mechancs Physcs 5 Lecture 0 Canoncal Transformatons (Chapter 9) What We Dd Last Tme Hamlton s Prncple n the Hamltonan formalsm Dervaton was smple δi δ p H(, p, t) = 0 Adonal end-pont constrants δ t ( )
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 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 informationInner Product. Euclidean Space. Orthonormal Basis. Orthogonal
Inner Product Defnton 1 () A Eucldean space s a fnte-dmensonal vector space over the reals R, wth an nner product,. Defnton 2 (Inner Product) An nner product, on a real vector space X s a symmetrc, blnear,
More informationLecture 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 informationLinear Regression Analysis: Terminology and Notation
ECON 35* -- Secton : Basc Concepts of Regresson Analyss (Page ) Lnear Regresson Analyss: Termnology and Notaton Consder the generc verson of the smple (two-varable) lnear regresson model. It s represented
More informationA Solution of Porous Media Equation
Internatonal Mathematcal Forum, Vol. 11, 016, no. 15, 71-733 HIKARI Ltd, www.m-hkar.com http://dx.do.org/10.1988/mf.016.6669 A Soluton of Porous Meda Equaton F. Fonseca Unversdad Naconal de Colomba Grupo
More informationProfessor Terje Haukaas University of British Columbia, Vancouver The Q4 Element
Professor Terje Haukaas Unversty of Brtsh Columba, ancouver www.nrsk.ubc.ca The Q Element Ths document consders fnte elements that carry load only n ther plane. These elements are sometmes referred to
More informationAn Algorithm to Solve the Inverse Kinematics Problem of a Robotic Manipulator Based on Rotation Vectors
An Algorthm to Solve the Inverse Knematcs Problem of a Robotc Manpulator Based on Rotaton Vectors Mohamad Z. Al-az*, Mazn Z. Othman**, and Baker B. Al-Bahr* *AL-Nahran Unversty, Computer Eng. Dep., Baghdad,
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 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 informationIntegrals and Invariants of Euler-Lagrange Equations
Lecture 16 Integrals and Invarants of Euler-Lagrange Equatons ME 256 at the Indan Insttute of Scence, Bengaluru Varatonal Methods and Structural Optmzaton G. K. Ananthasuresh Professor, Mechancal Engneerng,
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 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 informationModule 2. Random Processes. Version 2 ECE IIT, Kharagpur
Module Random Processes Lesson 6 Functons of Random Varables After readng ths lesson, ou wll learn about cdf of functon of a random varable. Formula for determnng the pdf of a random varable. Let, X be
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 informationOne-sided finite-difference approximations suitable for use with Richardson extrapolation
Journal of Computatonal Physcs 219 (2006) 13 20 Short note One-sded fnte-dfference approxmatons sutable for use wth Rchardson extrapolaton Kumar Rahul, S.N. Bhattacharyya * Department of Mechancal Engneerng,
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 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 informationPHYS 705: Classical Mechanics. Calculus of Variations II
1 PHYS 705: Classcal Mechancs Calculus of Varatons II 2 Calculus of Varatons: Generalzaton (no constrant yet) Suppose now that F depends on several dependent varables : We need to fnd such that has a statonary
More informationarxiv: v2 [math.na] 6 Dec 2016
Symmetry-Preservng Numercal Schemes arxv:1608.02557v2 [math.na] 6 Dec 2016 Alexander Bhlo Francs Valquette Department of Mathematcs and Statstcs Department of Mathematcs Memoral Unversty SUNY New Paltz
More informationA particle in a state of uniform motion remain in that state of motion unless acted upon by external force.
The fundamental prncples of classcal mechancs were lad down by Galleo and Newton n the 16th and 17th centures. In 1686, Newton wrote the Prncpa where he gave us three laws of moton, one law of gravty,
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 informationGeorgia Tech PHYS 6124 Mathematical Methods of Physics I
Georga Tech PHYS 624 Mathematcal Methods of Physcs I Instructor: Predrag Cvtanovć Fall semester 202 Homework Set #7 due October 30 202 == show all your work for maxmum credt == put labels ttle legends
More informationResearch Article Green s Theorem for Sign Data
Internatonal Scholarly Research Network ISRN Appled Mathematcs Volume 2012, Artcle ID 539359, 10 pages do:10.5402/2012/539359 Research Artcle Green s Theorem for Sgn Data Lous M. Houston The Unversty of
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 information(Online First)A Lattice Boltzmann Scheme for Diffusion Equation in Spherical Coordinate
Internatonal Journal of Mathematcs and Systems Scence (018) Volume 1 do:10.494/jmss.v1.815 (Onlne Frst)A Lattce Boltzmann Scheme for Dffuson Equaton n Sphercal Coordnate Debabrata Datta 1 *, T K Pal 1
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 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 information9 Characteristic classes
THEODORE VORONOV DIFFERENTIAL GEOMETRY. Sprng 2009 [under constructon] 9 Characterstc classes 9.1 The frst Chern class of a lne bundle Consder a complex vector bundle E B of rank p. We shall construct
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 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 informationPHYS 705: Classical Mechanics. Canonical Transformation II
1 PHYS 705: Classcal Mechancs Canoncal Transformaton II Example: Harmonc Oscllator f ( x) x m 0 x U( x) x mx x LT U m Defne or L p p mx x x m mx x H px L px p m p x m m H p 1 x m p m 1 m H x p m x m m
More informationAffine transformations and convexity
Affne transformatons and convexty The purpose of ths document s to prove some basc propertes of affne transformatons nvolvng convex sets. Here are a few onlne references for background nformaton: http://math.ucr.edu/
More informationQuantum Particle Motion in Physical Space
Adv. Studes Theor. Phys., Vol. 8, 014, no. 1, 7-34 HIKARI Ltd, www.-hkar.co http://dx.do.org/10.1988/astp.014.311136 Quantu Partcle Moton n Physcal Space A. Yu. Saarn Dept. of Physcs, Saara State Techncal
More informationLecture 3: Probability Distributions
Lecture 3: Probablty Dstrbutons Random Varables Let us begn by defnng a sample space as a set of outcomes from an experment. We denote ths by S. A random varable s a functon whch maps outcomes nto the
More informationON THE JACOBIAN CONJECTURE
v v v Far East Journal of Mathematcal Scences (FJMS) 17 Pushpa Publshng House, Allahabad, Inda http://www.pphm.com http://dx.do.org/1.17654/ms1111565 Volume 11, Number 11, 17, Pages 565-574 ISSN: 97-871
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 informationLagrangian Field Theory
Lagrangan Feld Theory Adam Lott PHY 391 Aprl 6, 017 1 Introducton Ths paper s a summary of Chapter of Mandl and Shaw s Quantum Feld Theory [1]. The frst thng to do s to fx the notaton. For the most part,
More informationPolynomial Regression Models
LINEAR REGRESSION ANALYSIS MODULE XII Lecture - 6 Polynomal Regresson Models Dr. Shalabh Department of Mathematcs and Statstcs Indan Insttute of Technology Kanpur Test of sgnfcance To test the sgnfcance
More informationBOUNDEDNESS OF THE RIESZ TRANSFORM WITH MATRIX A 2 WEIGHTS
BOUNDEDNESS OF THE IESZ TANSFOM WITH MATIX A WEIGHTS Introducton Let L = L ( n, be the functon space wth norm (ˆ f L = f(x C dx d < For a d d matrx valued functon W : wth W (x postve sem-defnte for all
More informationThermal-Fluids I. Chapter 18 Transient heat conduction. Dr. Primal Fernando Ph: (850)
hermal-fluds I Chapter 18 ransent heat conducton Dr. Prmal Fernando prmal@eng.fsu.edu Ph: (850) 410-6323 1 ransent heat conducton In general, he temperature of a body vares wth tme as well as poston. In
More informationFoundations of Arithmetic
Foundatons of Arthmetc Notaton We shall denote the sum and product of numbers n the usual notaton as a 2 + a 2 + a 3 + + a = a, a 1 a 2 a 3 a = a The notaton a b means a dvdes b,.e. ac = b where c s an
More informationCHAPTER 6. LAGRANGE S EQUATIONS (Analytical Mechanics)
CHAPTER 6 LAGRANGE S EQUATIONS (Analytcal Mechancs) 1 Ex. 1: Consder a partcle movng on a fxed horzontal surface. r P Let, be the poston and F be the total force on the partcle. The FBD s: -mgk F 1 x O
More informationTensor Smooth Length for SPH Modelling of High Speed Impact
Tensor Smooth Length for SPH Modellng of Hgh Speed Impact Roman Cherepanov and Alexander Gerasmov Insttute of Appled mathematcs and mechancs, Tomsk State Unversty 634050, Lenna av. 36, Tomsk, Russa RCherepanov82@gmal.com,Ger@npmm.tsu.ru
More informationMore metrics on cartesian products
More metrcs on cartesan products If (X, d ) are metrc spaces for 1 n, then n Secton II4 of the lecture notes we defned three metrcs on X whose underlyng topologes are the product topology The purpose of
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 informationModule 3: Element Properties Lecture 1: Natural Coordinates
Module 3: Element Propertes Lecture : Natural Coordnates Natural coordnate system s bascally a local coordnate system whch allows the specfcaton of a pont wthn the element by a set of dmensonless numbers
More information3.1 Expectation of Functions of Several Random Variables. )' be a k-dimensional discrete or continuous random vector, with joint PMF p (, E X E X1 E X
Statstcs 1: Probablty Theory II 37 3 EPECTATION OF SEVERAL RANDOM VARIABLES As n Probablty Theory I, the nterest n most stuatons les not on the actual dstrbuton of a random vector, but rather on a number
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 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 informationSupplementary Notes for Chapter 9 Mixture Thermodynamics
Supplementary Notes for Chapter 9 Mxture Thermodynamcs Key ponts Nne major topcs of Chapter 9 are revewed below: 1. Notaton and operatonal equatons for mxtures 2. PVTN EOSs for mxtures 3. General effects
More informationA PROCEDURE FOR SIMULATING THE NONLINEAR CONDUCTION HEAT TRANSFER IN A BODY WITH TEMPERATURE DEPENDENT THERMAL CONDUCTIVITY.
Proceedngs of the th Brazlan Congress of Thermal Scences and Engneerng -- ENCIT 006 Braz. Soc. of Mechancal Scences and Engneerng -- ABCM, Curtba, Brazl,- Dec. 5-8, 006 A PROCEDURE FOR SIMULATING THE NONLINEAR
More informationThe Order Relation and Trace Inequalities for. Hermitian Operators
Internatonal Mathematcal Forum, Vol 3, 08, no, 507-57 HIKARI Ltd, wwwm-hkarcom https://doorg/0988/mf088055 The Order Relaton and Trace Inequaltes for Hermtan Operators Y Huang School of Informaton Scence
More informationCOMPOSITE BEAM WITH WEAK SHEAR CONNECTION SUBJECTED TO THERMAL LOAD
COMPOSITE BEAM WITH WEAK SHEAR CONNECTION SUBJECTED TO THERMAL LOAD Ákos Jósef Lengyel, István Ecsed Assstant Lecturer, Professor of Mechancs, Insttute of Appled Mechancs, Unversty of Mskolc, Mskolc-Egyetemváros,
More informationInternational Journal of Pure and Applied Sciences and Technology
Int. J. Pure Appl. Sc. Technol., 4() (03), pp. 5-30 Internatonal Journal of Pure and Appled Scences and Technology ISSN 9-607 Avalable onlne at www.jopaasat.n Research Paper Schrödnger State Space Matrx
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 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 informationComputational Astrophysics
Computatonal Astrophyscs Solvng for Gravty Alexander Knebe, Unversdad Autonoma de Madrd Computatonal Astrophyscs Solvng for Gravty the equatons full set of equatons collsonless matter (e.g. dark matter
More informationcoordinates. Then, the position vectors are described by
Revewng, what we have dscussed so far: Generalzed coordnates Any number of varables (say, n) suffcent to specfy the confguraton of the system at each nstant to tme (need not be the mnmum number). In general,
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 informationMechanics Physics 151
Mechancs Physcs 151 Lecture 3 Lagrange s Equatons (Goldsten Chapter 1) Hamlton s Prncple (Chapter 2) What We Dd Last Tme! Dscussed mult-partcle systems! Internal and external forces! Laws of acton and
More informationAdvanced Quantum Mechanics
Advanced Quantum Mechancs Rajdeep Sensarma! sensarma@theory.tfr.res.n ecture #9 QM of Relatvstc Partcles Recap of ast Class Scalar Felds and orentz nvarant actons Complex Scalar Feld and Charge conjugaton
More informationIntegrals and Invariants of
Lecture 16 Integrals and Invarants of Euler Lagrange Equatons NPTEL Course Varatonal Methods and Structural Optmzaton G. K. Ananthasuresh Professor, Mechancal Engneerng, Indan Insttute of Scence, Banagalore
More informationOpen Systems: Chemical Potential and Partial Molar Quantities Chemical Potential
Open Systems: Chemcal Potental and Partal Molar Quanttes Chemcal Potental For closed systems, we have derved the followng relatonshps: du = TdS pdv dh = TdS + Vdp da = SdT pdv dg = VdP SdT For open systems,
More informationClassical Field Theory
Classcal Feld Theory Before we embark on quantzng an nteractng theory, we wll take a dverson nto classcal feld theory and classcal perturbaton theory and see how far we can get. The reader s expected to
More informationCinChE Problem-Solving Strategy Chapter 4 Development of a Mathematical Model. formulation. procedure
nhe roblem-solvng Strategy hapter 4 Transformaton rocess onceptual Model formulaton procedure Mathematcal Model The mathematcal model s an abstracton that represents the engneerng phenomena occurrng n
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 informationSolving Nonlinear Differential Equations by a Neural Network Method
Solvng Nonlnear Dfferental Equatons by a Neural Network Method Luce P. Aarts and Peter Van der Veer Delft Unversty of Technology, Faculty of Cvlengneerng and Geoscences, Secton of Cvlengneerng Informatcs,
More informationSo far: simple (planar) geometries
Physcs 06 ecture 5 Torque and Angular Momentum as Vectors SJ 7thEd.: Chap. to 3 Rotatonal quanttes as vectors Cross product Torque epressed as a vector Angular momentum defned Angular momentum as a vector
More informationPower law and dimension of the maximum value for belief distribution with the max Deng entropy
Power law and dmenson of the maxmum value for belef dstrbuton wth the max Deng entropy Bngy Kang a, a College of Informaton Engneerng, Northwest A&F Unversty, Yanglng, Shaanx, 712100, Chna. Abstract Deng
More information