Modified Mass Matrices and Positivity Preservation for Hyperbolic and Parabolic PDEs
|
|
- Marian Bryant
- 5 years ago
- Views:
Transcription
1 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 Parabolc PDEs M. Berzns Λ Computatonal PDEs Unt, School of Computng, The Unversty of Leeds, Leeds LS2 9JT, UK. SUMMARY Modfcatons to the standard fnte element mass matrx are consdered wth the am of preservng the postvty of the dscrete soluton. The approach s used n connecton wth calculatng the ntal tme dervatve values for parabolc equatons and n connecton wth nonlnear Petrov-Galerkn schemes for hyperbolc equatons n one space dmenson. The extenson of the deas to unstructured meshes n two and three space dmensons s ndcated. Copyrght cfl 2000 John Wley & Sons, Ltd. KEY WORDS: Postvty Preservaton ; Mass Matrces, Fnte Element Methods. Introducton There are many stuatons n the numercal soluton of partal dfferental equatons n whch the computed soluton values should, on physcal grounds, reman non-negatve. One the smplest examples s that of the smple advecton equaton wth non-negatve ntal data whle other cases are those of concentratons of chemcal compounds n reactng flow calculatons. In the latter case preservng postvty s essental to avod the numercal calculaton becomng meanngless. Consder the soluton of the advecton equaton wth approprate ntal and boundary condton by usng the standard Galerkn method wth lnear bass (hat) functons φ (x) on a unformly spaced mesh x ; ;:::;N to get Z x+ x U t φ (x) dx Z x+ x U x φ (x) dx; ; :::;N; () where the approxmate soluton to ths p.d.e. as defned by U(x;t) N φ (x)u (t) where φ (x )δ. Evaluatng the ntegrals gves rse to the numercal scheme defned by Λ δx (U + U ) (2) Λ Correspondence to: M.Berzns, School of Computng, The Unversty of Leeds, Leeds LS2 9JT, UK. Receved 0 January 200 Copyrght cfl 2000 John Wley & Sons, Ltd. Revsed 22 March 200
2 2 M. BERZINS where δx s the unform mesh spacng n ths case and where du dt. Defnng the tme-dependent vector U by U [U ;:::;U N ] T allows ths system of equatons to be rewrtten n the form A (t) F(U(t)) (3) where the matrx A s referred to as the mass matrx. It s well-known that ths scheme s unsatsfactory n a very smlar way to that of lnear central dfference schemes, [8]. Many modfed Galerkn methods have been proposed to remedy ths stuaton. A survey of such methods s gven n [8] and ncludes Streamlne Upwnd Petrov-Galerkn (SUPG) methods [7] n whch the test functons are modfed to mprove the behavour of the method and Dscontnuous Galerkn (DG) methods [4, 7] n whch dscontnuous bass functons are used. There are many other approaches such as the modfed Petrov-Galerkn method of Cardle [3] n whch the test functon s modfed dfferently for the spatal and temporal terms. In ths case the numercal scheme that results s gven by + 6 ( β Λ ) δx (U + U α Λ )+ U 2U 2δx +U + (4) where β and α are the constants multplyng the Petrov-Galerkn addtonal polynomals n tme (cubc polynomal) and space (quadratc polynomal), see [3]. In the case of many of these methods t s clear that the magntude of unphyscal values s not as large as wth the standard Galerkn method and n the case of DG methods the mass matrx s the dentty matrx; ths makes t much easer to prove propertes such as postvty preservaton. The defnton used here for a postvty preservng scheme for the advecton equaton s one (see []) for whch the numercal soluton at tme t n+ may be wrtten n terms of the numercal soluton at tme t n n the form U (t n+ ) a U (t n ) where a ; and a 0 : (5) The key observaton wth regard to preservng postvty s due to Godunov [6] who proved that any scheme of better than frst order whch preserves postvty for the advecton equaton must be nonlnear. That s the coeffcents a n (5) above must depend on the numercal soluton to the p.d.e. Ths means that α and β n (4) must also depend on the soluton. In nvestgatng postvty preservng mass matrces and Galerkn fnte element methods for transent problems the startng pont wll be to rewrte the mass matrx as a postvty preservng matrx. Ths wll be appled to the soluton of a parabolc equaton. The same dea wll then be appled to hyperbolc equatons and lnked to the work of Cardle, [3] and to work on nonlnear fnte dfference schemes. Fnally the extenson of the approach to unstructured trangular and tetrahedral meshes wll be consdered. 2. Modfed Mass Matrces and the Intalsaton of Parabolc Equatons In tryng to solve parabolc equatons usng a Galerkn method-of-lnes approach Skeel and Berzns [9] showed that the ntal tme dervatves may have the wrong sgn. Ths s because the nverse of the mass matrx A may have have negatve entres. Suppose that the parabolc equaton s dscretsed n space to get a system of equatons of the form of (3) wth F 0. The ntal values of the tme dervatves are gven by solvng the equatons (3) for the ntal values of the tme dervatves (0). In computatonal experments tme dervatves wth the wrong sgn slow down the the ntegraton and Copyrght cfl 2000 John Wley & Sons, Ltd. Commun. Numer. Meth. Engng 2000; 00:6 Prepared usng cnmauth.cls
3 MODIFIED MASS MATRICES AND POSITIVITY PRESERVATION FOR PDES 3 may gve physcally msleadng soluton values. For these reasons Skeel and Berzns devsed a scheme that may be vewed as a lumped fnte element scheme n whch the mass matrx s replaced by the dentty matrx. The ssue of when a matrx may have an nverse consstng of postve entres s consdered n a large body of work on M matrces. See, for example, [2] who show that f A s a dagonally domnant M matrx wth negatve off dagonal entres then ts nverse A has only postve entres. The task s thus to modfy the mass matrx so that t has negatve off-dagonal entres. 2.. Dervaton of Modfed Mass Matrx For smplcty consder the case when lnear bass functons are used on a unform spatal mesh as n equatons (). The th row of the mass matrx s then gven by φ ;φ + + φ ;φ + + φ ;φ F : (6) Usng the dentty that on [x ;x + ] φ + φ + φ + gves: φ ; + + φ ;φ + + φ ;φ F : (7) + Defnng the rato s allows the th row of the mass matrx to be rewrtten as 2 3 φ ;φ φ ; φ ;φ 5 + F : (8) s Ths matrx s an M matrx f [:::] s postve (on a unform mesh ths requres 0 < s < ) as then the matrx s dagonally domnant wth negatve off-dagonal entres. In the case when s > on a unform mesh the th row s wrtten as φ ; + hφ ;φ + s φ ;φ F (9) whch agan s a row of an M matrx as [:::] s postve. In the case when s < 0 there appears h no alternatve but to dagonalse (lump) the matrx as φ ;φ + φ ;φ + φ ;φ +. The modfed matrx may also be wrtten as (wth approprate modfcatons f s 0) : (s + s ) h + 2s F δx : (0) In solvng the equatons (8) and (9) for the ntal values of the tme dervatves t s thus necessary to teratvely solve nonlnear equatons. Let m and s m be the values calculated at teraton m. The equatons solved n the case when s m > are then gven by m+ + γ m+ m In the case when 0 < s m < the teraton s defned by m+ + γ m+ + m F δx ; where γ F δx ; where γ» s m 6» s m 6 : () : (2) Copyrght cfl 2000 John Wley & Sons, Ltd. Commun. Numer. Meth. Engng 2000; 00:6 Prepared usng cnmauth.cls
4 4 M. BERZINS In order to llustrate that ths procedure produces ntal values of the tme dervatves wth the rght sgn the followng example s used. Skeel and Berzns [9] consder test examples such as the case when C the rght sde of equaton (3) s gven by F ß U t (x ;0) where C 0:fx < 0 and C :0 (Cx +:) otherwse. For a unform mesh of and 2 ponts across the nterval [;] Fgure shows that the.4 Mass matrx solve n.4 Mass matrx solve n Fgure. Mass Matrx Calculaton for Tme Dervatves - s true * s new, s orgnal method does result n tme dervatves of the correct sgn wthout the overshoots and undershoots of a standard Galerkn approach. The ssue of preservng postvty for the dffuson equaton has been consdered n very recent work by Farago and Horvath [5]. They show that for a d dmensonal problem usng the θ tme ntegraton method, t s necessary to restrct the choce of the parameter θ and the tmestep δt by 6θ d < δt δx 2 < 3d(θ ; where d ;2 Ther results also extends to three space dmensons. ) 3. Modfed Mass Matrces and Hyperbolc Equatons In general the tme dervatves may not have constant sgn and t s the non-negatvty of the soluton that must be consdered. In the case of the advecton equaton t s possble to rewrte equaton (4) as a smple explct method for hyperbolc equatons: AU (t n+ )AU (t n )+ δt F(U(t n )) (3) In usng the modfed mass matrx approach as part of a method for hyperbolc equatons t s nstructve to note other smlar approaches. The approxmaton of U x by a standard Galerkn method s dentcal to that of central fnte dfferences. A nonlnear central dfference method for hyperbolc equatons s gven by Swanson and Turkel [0] as : 2δx (U + U h )+ L 2δx + (U + U ) L (U U ) (4) where L ˆr +ˆr and ˆr U U 2 U U. The rght sde of ths may also be nterpreted as a nonlnear Petrov-Galerkn method n whch the test functon s ˆφ (x) where Λ dφ (x) ˆφ (x) φ (x)+δx L dx (5) Copyrght cfl 2000 John Wley & Sons, Ltd. Commun. Numer. Meth. Engng 2000; 00:6 Prepared usng cnmauth.cls
5 MODIFIED MASS MATRICES AND POSITIVITY PRESERVATION FOR PDES 5 where L Λ L f x < x < x + and L Λ L + otherwse. Although ths scheme s postvty preservng t s qute dffuse. A less dffuse scheme s one n whch L Λ s defned by L Λ V (r )r r +r +r f x + > x > x and L Λ V(r )r f x > x > x where V (r ) Routne manpulaton shows that ths defnton gves the well-known van Leer scheme, e.g []: " and r U + U U U. U (t) δx β (U (t) U (t)); where β V (r ) V (r ) + : (6) 2 2r An alternatve vew of ths scheme s thus as a nonlnear Petrov-Galerkn method. Requrng postvty for forward Euler tmesteppng requres the CFL type restrcton 0» β < δx δt. The nonlnear extenson of the type of Petrov-Galerkn method gven by (3) s thus gven by equaton (0) wth F defned by the rght sde of (6). An outlne proof that, when combned wth forward Euler tmesteppng, ths s postvty preservng follows from a modfed verson of (). (The proof for the case n (2) beng smlar). The teraton from () may be wrtten as: m+ (t) + γ» γ m (t)+ F δx ; m 0;;::: where γ s defned as n equaton () and so depends on s m and where F δx s defned by the rghthand sde of equaton (6). Hence ths equaton may also be wrtten as m+ (t) + γ " γ m (t) β δx (U (t) U (t)) # # ; m 0;;::: (7) The predcted values, 0 (t n+ ), are gven by equaton(6). An outlne of the approach used to defne soluton postvty n terms of equaton (5) may now be gven. The ntal guesses for the tme drvatves are gven by equaton (6): 0 (t h n+ ) U δx (t n ) U 2 (t n ) β (8) where β s also defned as n (6). Substtutng for 0 n (7) and applyng forward Euler tmesteppng gves: U (t n+ ) U (t δt n) δx( + γ) h β (U (t n ) U (t n )) + γβ (U (t n ) U 2 (t n )) where U (t n+ ) s the frst teraton estmate for U (t n+ ). Ths may be rewrtten as U (t n+ ) ( β E)U (t n)+(β γβ )EU (t n )+(γβ )EU 2 (t n ) (20) where E δt δx ( + γ) : In provng postvty of ths consder the worst case n equaton (9) and suppose that β (U (t n ) U (t n )) has a dfferent sgn to β (U (t n ) U 2 (t n )). It then follows that r s negatve and hence that β h + V (r ) 2 and β» V (r 2 ) 2r 2. Hence V (r ) β γβ + γ + γ V (r 2 ) 2 2r 2 (9) Copyrght cfl 2000 John Wley & Sons, Ltd. Commun. Numer. Meth. Engng 2000; 00:6 Prepared usng cnmauth.cls
6 6 M. BERZINS and n order to guarantee a postve soluton at the end of the frst teraton we need to mpose the condton γ < :0. Thus, from equatons () and (2), the teratve method s only appled f 7 < sm < 7 The same approach can then be used nductvely to prove postvty for second and subsequent teratons. Fgure 2 shows the numercal results obtaned when the method s appled to the advecton of a square pulse functon from x 0:2 to0:7 usng a spatal mesh of 5 equally spaced ponts. The fgure compares the soluton obtaned wth the van Leer method wth the new approach and shows the effect of usng a mass matrx s to gve a stll postve but more skewed profle than that obtaned from the van Leer method. Advecton cfl 50 Advecton cfl Solutons at t Solutons at t x 0. 0 x * Exact, + van Leer FV, new FE Fgure 2. Numercal soluton of advected square wave 4. Extenson to Trangular and Tetrahedral Meshes In the case of trangular mesh examples t s possble to use the same dea as n one space dmenson. In ths case, for example, the mass matrx for a mesh fragment consstng of three trangles wth node n common and a permeter consstng of nodes, k and l s gven by: (φ ;)+ φ ;φ + k φk ;φ + l φl ;φ (2) where φ ;φ 2 ha k + A l, φ ;φ 2 ha k + A l + A lk and where A k s the area of trangle wth nodes, and k. As the contrbuton of element k to the mass matrx s A k 2 h k : (22) Copyrght cfl 2000 John Wley & Sons, Ltd. Commun. Numer. Meth. Engng 2000; 00:6 Prepared usng cnmauth.cls
7 MODIFIED MASS MATRICES AND POSITIVITY PRESERVATION FOR PDES 7 The same deas as n one space dmenson may be used to rewrte 2 + k by usng the same approach as n Secton 2.. Consder the ; ;k trangle and let ( +k)2 be mdpont soluton value on k edge. The decomposton gven by h h 2 + k 2 + ( +k)2) k + 2 ( +k)2) (23) allows the terms on the rght sde of ths equaton to be vewed as second order approxmatons to second and frst space dervatves. Hence n dscardng these terms we ntroduce second order errors as n one space dmenson. The same dea extends to tetrahedral mesh examples. Consder a sngle tetrahedron wth ts four nodes labelled as ; ;k;l and assocated lnear bass functons φ ;φ ;φ k and φ l. The mass matrx ntegral assocated wth ths tetrahedron s h V kl (φ ;φ )+ φ ;φ + k φk ;φ + l φl ;φ : (24) Evaluatng the volume ntegral gves: V kl 20 h k + l whch may be rewrtten as three terms of the form V kl 20 appled as n one and two space dmensons. (25) h k and same deas 5. Summary In ths paper a novel approach to preservng postvty has been taken. The approach reles on usng a nonlnear form of the mass matrx n conuncton wth nonlnear Petrov-Galerkn type terms. The approach has applcatons n one,two and three space dmensonal cases, but further work s clearly needed to assess ts usefulness. REFERENCES. Berzns M, Ware JM. Postve Cell Centered Fnte Volume Dscretsaton Methods for Hyperbolc Equatons on Irregular Meshes, Appled Numercal Mathematcs 995; 6: Berman A, Plemmons RJ. Nonnegatve Matrces n the Mathematcal Scences. Acadedmc Press, Cardle JA. A modfcaton of the Petrov-Galerkn method for the transent convecton-dffuson equaton. Internatonal Journal for Numercal Methods n Engneerng 995; 38: Cockburn B, Karnadaks GE, Shu C-W. (eds). Dscontnuous Galerkn Methods, Theory Computatons and Applcatons. Lecture Notes n Computataonal Scence and Engneerng Sprnger Berln Hedelberg, 2000; Farago I, Horvath R. On the nonnegatvty conservaton of fnte element soluton of parabolc problems. In 3D Fnte Element Conference Proceedngs Krzek M., et al (eds). Gakkotosho Co., Tokyo Godunov SK. Fnte dfference method for the numercal computaton of dscontnuous solutons of the equatons of flud dynamcs. Math Sbornk 959; 47: Johnson C. Numercal Soluton of Partal Dfferental Equatons by the Fnte Element Method Cambrdge Unversty Press, Segal A. Fnte element methods for advecton-dffuson equatons. In Numercal Methods for Advecton Dffuson Problems. Notes on Numercal Flud Mechancs, Volume 45 Vreugdenhl CB, Koren B (eds). Veweg: Braunscweg/Wesbaden, 993; Skeel RD, Berzns M. A Method for the Spatal Dscretsaton of Parabolc Equatons, SIAM Journal on Scentfc Computng 990; (): Swanson RC, Turkel E. On central-dfference and upwnd schemes. Journal of Computatonal Physcs 992; 0: Copyrght cfl 2000 John Wley & Sons, Ltd. Commun. Numer. Meth. Engng 2000; 00:6 Prepared usng cnmauth.cls
Appendix 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 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 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 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 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 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 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 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 informationWorkshop: Approximating energies and wave functions Quantum aspects of physical chemistry
Workshop: Approxmatng energes and wave functons Quantum aspects of physcal chemstry http://quantum.bu.edu/pltl/6/6.pdf Last updated Thursday, November 7, 25 7:9:5-5: Copyrght 25 Dan Dll (dan@bu.edu) Department
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 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 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 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 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 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 informationBézier curves. Michael S. Floater. September 10, These notes provide an introduction to Bézier curves. i=0
Bézer curves Mchael S. Floater September 1, 215 These notes provde an ntroducton to Bézer curves. 1 Bernsten polynomals Recall that a real polynomal of a real varable x R, wth degree n, s a functon of
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 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 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 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 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 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 informationCollege of Computer & Information Science Fall 2009 Northeastern University 20 October 2009
College of Computer & Informaton Scence Fall 2009 Northeastern Unversty 20 October 2009 CS7880: Algorthmc Power Tools Scrbe: Jan Wen and Laura Poplawsk Lecture Outlne: Prmal-dual schema Network Desgn:
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 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 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 informationKernel Methods and SVMs Extension
Kernel Methods and SVMs Extenson The purpose of ths document s to revew materal covered n Machne Learnng 1 Supervsed Learnng regardng support vector machnes (SVMs). Ths document also provdes a general
More informationThe Finite Element Method
The Fnte Element Method GENERAL INTRODUCTION Read: Chapters 1 and 2 CONTENTS Engneerng and analyss Smulaton of a physcal process Examples mathematcal model development Approxmate solutons and methods of
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 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 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 informationLinear Approximation with Regularization and Moving Least Squares
Lnear Approxmaton wth Regularzaton and Movng Least Squares Igor Grešovn May 007 Revson 4.6 (Revson : March 004). 5 4 3 0.5 3 3.5 4 Contents: Lnear Fttng...4. Weghted Least Squares n Functon Approxmaton...
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 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 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 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 informationA New Refinement of Jacobi Method for Solution of Linear System Equations AX=b
Int J Contemp Math Scences, Vol 3, 28, no 17, 819-827 A New Refnement of Jacob Method for Soluton of Lnear System Equatons AX=b F Naem Dafchah Department of Mathematcs, Faculty of Scences Unversty of Gulan,
More informationSolving Fractional Nonlinear Fredholm Integro-differential Equations via Hybrid of Rationalized Haar Functions
ISSN 746-7659 England UK Journal of Informaton and Computng Scence Vol. 9 No. 3 4 pp. 69-8 Solvng Fractonal Nonlnear Fredholm Integro-dfferental Equatons va Hybrd of Ratonalzed Haar Functons Yadollah Ordokhan
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 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 informationU.C. Berkeley CS294: Beyond Worst-Case Analysis Luca Trevisan September 5, 2017
U.C. Berkeley CS94: Beyond Worst-Case Analyss Handout 4s Luca Trevsan September 5, 07 Summary of Lecture 4 In whch we ntroduce semdefnte programmng and apply t to Max Cut. Semdefnte Programmng Recall that
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 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 informationTime-Varying Systems and Computations Lecture 6
Tme-Varyng Systems and Computatons Lecture 6 Klaus Depold 14. Januar 2014 The Kalman Flter The Kalman estmaton flter attempts to estmate the actual state of an unknown dscrete dynamcal system, gven nosy
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 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 informationBezier curves. Michael S. Floater. August 25, These notes provide an introduction to Bezier curves. i=0
Bezer curves Mchael S. Floater August 25, 211 These notes provde an ntroducton to Bezer curves. 1 Bernsten polynomals Recall that a real polynomal of a real varable x R, wth degree n, s a functon of the
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 informationThe Exact Formulation of the Inverse of the Tridiagonal Matrix for Solving the 1D Poisson Equation with the Finite Difference Method
Journal of Electromagnetc Analyss and Applcatons, 04, 6, 0-08 Publshed Onlne September 04 n ScRes. http://www.scrp.org/journal/jemaa http://dx.do.org/0.46/jemaa.04.6000 The Exact Formulaton of the Inverse
More informationBeyond Zudilin s Conjectured q-analog of Schmidt s problem
Beyond Zudln s Conectured q-analog of Schmdt s problem Thotsaporn Ae Thanatpanonda thotsaporn@gmalcom Mathematcs Subect Classfcaton: 11B65 33B99 Abstract Usng the methodology of (rgorous expermental mathematcs
More informationNormally, in one phase reservoir simulation we would deal with one of the following fluid systems:
TPG4160 Reservor Smulaton 2017 page 1 of 9 ONE-DIMENSIONAL, ONE-PHASE RESERVOIR SIMULATION Flud systems The term sngle phase apples to any system wth only one phase present n the reservor In some cases
More informationStructure and Drive Paul A. Jensen Copyright July 20, 2003
Structure and Drve Paul A. Jensen Copyrght July 20, 2003 A system s made up of several operatons wth flow passng between them. The structure of the system descrbes the flow paths from nputs to outputs.
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 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 informationA MODIFIED METHOD FOR SOLVING SYSTEM OF NONLINEAR EQUATIONS
Journal of Mathematcs and Statstcs 9 (1): 4-8, 1 ISSN 1549-644 1 Scence Publcatons do:1.844/jmssp.1.4.8 Publshed Onlne 9 (1) 1 (http://www.thescpub.com/jmss.toc) A MODIFIED METHOD FOR SOLVING SYSTEM OF
More informationAPPENDIX A Some Linear Algebra
APPENDIX A Some Lnear Algebra The collecton of m, n matrces A.1 Matrces a 1,1,..., a 1,n A = a m,1,..., a m,n wth real elements a,j s denoted by R m,n. If n = 1 then A s called a column vector. Smlarly,
More informationIntroduction to Vapor/Liquid Equilibrium, part 2. Raoult s Law:
CE304, Sprng 2004 Lecture 4 Introducton to Vapor/Lqud Equlbrum, part 2 Raoult s Law: The smplest model that allows us do VLE calculatons s obtaned when we assume that the vapor phase s an deal gas, and
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 informationNumerical methods for hyperbolic balance laws with discontinuous flux functions and applications in radiotherapy
Numercal methods for hyperbolc balance laws wth dscontnuous flux functons and applcatons n radotherapy Numersche Verfahren für hyperbolsche Blanzglechungen mt unstetger Flussfunkton und Anwendungen n der
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 informationCubic Trigonometric B-Spline Applied to Linear Two-Point Boundary Value Problems of Order Two
World Academy of Scence Engneerng and echnology Internatonal Journal of Mathematcal and omputatonal Scences Vol: No:0 00 ubc rgonometrc B-Splne Appled to Lnear wo-pont Boundary Value Problems of Order
More informationVARIATION OF CONSTANT SUM CONSTRAINT FOR INTEGER MODEL WITH NON UNIFORM VARIABLES
VARIATION OF CONSTANT SUM CONSTRAINT FOR INTEGER MODEL WITH NON UNIFORM VARIABLES BÂRZĂ, Slvu Faculty of Mathematcs-Informatcs Spru Haret Unversty barza_slvu@yahoo.com Abstract Ths paper wants to contnue
More informationLinear Algebra and its Applications
Lnear Algebra and ts Applcatons 4 (00) 5 56 Contents lsts avalable at ScenceDrect Lnear Algebra and ts Applcatons journal homepage: wwwelsevercom/locate/laa Notes on Hlbert and Cauchy matrces Mroslav Fedler
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 informationOn the Interval Zoro Symmetric Single-step Procedure for Simultaneous Finding of Polynomial Zeros
Appled Mathematcal Scences, Vol. 5, 2011, no. 75, 3693-3706 On the Interval Zoro Symmetrc Sngle-step Procedure for Smultaneous Fndng of Polynomal Zeros S. F. M. Rusl, M. Mons, M. A. Hassan and W. J. Leong
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 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 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 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 informationCHAPTER 4d. ROOTS OF EQUATIONS
CHAPTER 4d. ROOTS OF EQUATIONS A. J. Clark School o Engneerng Department o Cvl and Envronmental Engneerng by Dr. Ibrahm A. Assakka Sprng 00 ENCE 03 - Computaton Methods n Cvl Engneerng II Department o
More informationYong Joon Ryang. 1. Introduction Consider the multicommodity transportation problem with convex quadratic cost function. 1 2 (x x0 ) T Q(x x 0 )
Kangweon-Kyungk Math. Jour. 4 1996), No. 1, pp. 7 16 AN ITERATIVE ROW-ACTION METHOD FOR MULTICOMMODITY TRANSPORTATION PROBLEMS Yong Joon Ryang Abstract. The optmzaton problems wth quadratc constrants often
More informationA Cartesian-grid integrated-rbf method for viscoelastic flows
Home Search Collectons Journals About Contact us My IOPscence A Cartesan-grd ntegrated-rbf method for vscoelastc flows Ths artcle has been downloaded from IOPscence. Please scroll down to see the full
More informationDeriving the X-Z Identity from Auxiliary Space Method
Dervng the X-Z Identty from Auxlary Space Method Long Chen Department of Mathematcs, Unversty of Calforna at Irvne, Irvne, CA 92697 chenlong@math.uc.edu 1 Iteratve Methods In ths paper we dscuss teratve
More information1 Matrix representations of canonical matrices
1 Matrx representatons of canoncal matrces 2-d rotaton around the orgn: ( ) cos θ sn θ R 0 = sn θ cos θ 3-d rotaton around the x-axs: R x = 1 0 0 0 cos θ sn θ 0 sn θ cos θ 3-d rotaton around the y-axs:
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 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 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 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 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 informationThe Multiple Classical Linear Regression Model (CLRM): Specification and Assumptions. 1. Introduction
ECONOMICS 5* -- NOTE (Summary) ECON 5* -- NOTE The Multple Classcal Lnear Regresson Model (CLRM): Specfcaton and Assumptons. Introducton CLRM stands for the Classcal Lnear Regresson Model. The CLRM s also
More informationAPPROXIMATE PRICES OF BASKET AND ASIAN OPTIONS DUPONT OLIVIER. Premia 14
APPROXIMAE PRICES OF BASKE AND ASIAN OPIONS DUPON OLIVIER Prema 14 Contents Introducton 1 1. Framewor 1 1.1. Baset optons 1.. Asan optons. Computng the prce 3. Lower bound 3.1. Closed formula for the prce
More informationLeast squares cubic splines without B-splines S.K. Lucas
Least squares cubc splnes wthout B-splnes S.K. Lucas School of Mathematcs and Statstcs, Unversty of South Australa, Mawson Lakes SA 595 e-mal: stephen.lucas@unsa.edu.au Submtted to the Gazette of the Australan
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 informationFinite Element Modelling of truss/cable structures
Pet Schreurs Endhoven Unversty of echnology Department of Mechancal Engneerng Materals echnology November 3, 214 Fnte Element Modellng of truss/cable structures 1 Fnte Element Analyss of prestressed structures
More informationCOMPUTATIONAL METHODS AND ALGORITHMS Vol. II - Finite Element Method - Jacques-Hervé SAIAC
COMPUTATIONAL METHODS AND ALGORITHMS Vol. II - Fnte Element Method - Jacques-Hervé SAIAC FINITE ELEMENT METHOD Jacques-Hervé SAIAC Départment de Mathématques, Conservatore Natonal des Arts et Méters, Pars,
More informationSTUDY ON TWO PHASE FLOW IN MICRO CHANNEL BASED ON EXPERI- MENTS AND NUMERICAL EXAMINATIONS
Blucher Mechancal Engneerng Proceedngs May 0, vol., num. www.proceedngs.blucher.com.br/evento/0wccm STUDY ON TWO PHASE FLOW IN MICRO CHANNEL BASED ON EXPERI- MENTS AND NUMERICAL EXAMINATIONS Takahko Kurahash,
More informationThe lower and upper bounds on Perron root of nonnegative irreducible matrices
Journal of Computatonal Appled Mathematcs 217 (2008) 259 267 wwwelsevercom/locate/cam The lower upper bounds on Perron root of nonnegatve rreducble matrces Guang-Xn Huang a,, Feng Yn b,keguo a a College
More informationMATH 241B FUNCTIONAL ANALYSIS - NOTES EXAMPLES OF C ALGEBRAS
MATH 241B FUNCTIONAL ANALYSIS - NOTES EXAMPLES OF C ALGEBRAS These are nformal notes whch cover some of the materal whch s not n the course book. The man purpose s to gve a number of nontrval examples
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 informationCIVL 8/7117 Chapter 10 - Isoparametric Formulation 42/56
CIVL 8/77 Chapter 0 - Isoparametrc Formulaton 4/56 Newton-Cotes Example Usng the Newton-Cotes method wth = ntervals (n = 3 samplng ponts), evaluate the ntegrals: x x cos dx 3 x x dx 3 x x dx 4.3333333
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 informationMaejo International Journal of Science and Technology
Maejo Int. J. Sc. Technol. () - Full Paper Maejo Internatonal Journal of Scence and Technology ISSN - Avalable onlne at www.mjst.mju.ac.th Fourth-order method for sngularly perturbed sngular boundary value
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 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 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 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 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 informationResearch Article Cubic B-Spline Collocation Method for One-Dimensional Heat and Advection-Diffusion Equations
Appled Mathematcs Volume 22, Artcle ID 4587, 8 pages do:.55/22/4587 Research Artcle Cubc B-Splne Collocaton Method for One-Dmensonal Heat and Advecton-Dffuson Equatons Joan Goh, Ahmad Abd. Majd, and Ahmad
More information= = = (a) Use the MATLAB command rref to solve the system. (b) Let A be the coefficient matrix and B be the right-hand side of the system.
Chapter Matlab Exercses Chapter Matlab Exercses. Consder the lnear system of Example n Secton.. x x x y z y y z (a) Use the MATLAB command rref to solve the system. (b) Let A be the coeffcent matrx and
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 information