Mass Lumping for Constant Density Acoustics
|
|
- Albert Cain
- 5 years ago
- Views:
Transcription
1 Lumping 1 Mass Lumping for Constant Density Acoustics William W. Symes ABSTRACT Mass lumping provides an avenue for efficient time-stepping of time-dependent problems wit conforming finite element spatial discretization. Typical justifications for mass lumping use quadrature error estimates wic do not old for nonsmoot coefficients. In tis paper, I sow tat te mass-lumped semidiscrete system for te constant-density acoustic wave equation wit Q 1 elements exibits optimal order convergence even wen te coefficient bulk modulus is merely bounded and measurable, provided tat te rigt-and side possesses some smootness in time. INTRODUCTION Mass lumping is te process of replacing te mass matrix, occurring in a finite element formulation of a time-dependent or eigenvalue problem, wit a diagonal matrix. Suc diagonal replacement is essential in time-domain finite element approximation of wave propagation: time steps must be fine enoug to resolve te oscillations caracteristic of wave problems, and solution of a nontrivial linear system at every time step is proibitively expensive for large-scale D or 3D simulations te reader can take tis expense as a definition of large-scale. Coen 001 describes many instances of tis construction for spectral C 0 element approximation of acoustic, elastic, and electromagnetic waves for example. According to standard teory for conforming spectral FEM Coen, 001, lumping is possible witout loss of accuracy troug use of a Gauss-Lobatto quadrature rule wit nodes coinciding wit tose of te nodal conforming basis of te metod. Tis approac to mass-lumping relies for its justification on quadrature error estimates suc as tat presented by?, problem 4.1.5: one obtains an error of order O k in te mass integral over a regular familly of elements of diameter for sufficiently regular trial solutions, provided tat te quadrature rule is exact for polynomials of degree k 1, and te multiplier mass density or oter coefficient belongs to te space W k,. For lowest order k = 1, energy errors of order follow so long as te quadrature rule is exact for linear functions and te coefficient in te mass integral is of class W 1,. Tis justification fails wen te coefficients factors in te mass matrix integrand lack sufficient regularity, for instance exibit discontinuities along piecewise smoot interfaces. Te purpose of tis note is to sketc a different argument for lowest order masslumping, wic olds for coefficients tat are merely bounded and measurable, and wic
2 Symes preserves convergence wit optimal order O in energy. Te essential ingredient is sufficient smootness of te solution to yield optimal order convergence in energy of te Galerkin approximation. Tese two caracteristics are somewat in tension: lack of regularity in te coefficients limits te regularity of te solution. In tis note I study te simplest special case, te constant-density acoustic wave equation, for wic te required degree of regularity is easily establised for solutions wic are smoot in time. Smootness in time is a natural feature of some problems, suc as simulation of seismograms, in wic wavefields ave fixed temporal bandwidt. Time regularity induces just enoug spatial regularity in te constant-density case to give bot optimal order convergence in energy and te same order of convergence for te mass-lumped approximation. A by-product of te argument is te observation tat wit additional coefficient regularity a global LIpsitz bound, for example, te error in te lumped-matrix approximation is te same, to leading order, as in te consistent approximation. In tis report, I demonstrate tis stricter approximation property only for 1D problems. Te constant density acoustic wave equation is a very special case. Even for te variable density wave equation, coefficient discontinuities generally destroy te regularity underlying te optimal order error estimate, even for lowest-order elements, along wit te corresponding convergence rate for te mass-lumped system. It is owever possible to coose special elements for wic optimal order convergence is restored. For isolated discontinuities along interfaces, tis approac is known as immersed finite elements IFE Li and Ito, 006. Owadi as recently sown ow to construct suitable elements for very general self-adjoint elliptic second order problems, vastly generalizing te IFE metod Owadi and Zang, 007; Owadi and Zang, 008. It seems possible to base a mass-lumping strategy for more general wave problems on tese concepts. Tis report begins wit a review of te properties of weak solutions of te acoustic wave equation. In te next two sections I establis general properties of te Galerkin approximation, and derive an error estimate for for mass lumping wit Q 1 elements on a regular mes. Te more refined estimate for Lipsitz coefficients in 1D is te topic of te next section. I end wit a brief discussion of prospects for generalization. An Appendix reviews te standard error estimate for Galerkin metods applied to te wave equation. WEAK SOLUTIONS Te constant density acoustic wave equation relates te fields ux, t, ηx, and fx, t troug η u t u = f. 1 In te acoustic setting, u is excess pressure; η is related to te sound velocity c, bulk modulus κ, and density ρ by η = c = ρ/κ. Te rigt-and side f is eiter te divergence of a specific body force or te time rate of cange of a constitutive law defect. I assume tat u and f are causal, tat is, vanis for large negative times, and tat te wave equation 1 olds in te domain Ω R, Ω R d. Te field u vanises on te
3 Lumping 3 boundary Ω. Te weak form of 1 is: dt ηu, d dt φ + au, φ = dt f, φ, φ C 0R, H 1 0Ω, in wic, is te inner product in L Ω and a, is te usual energy form: av, w = v w 3 Weak finite energy causal solutions of 1 are known to exist and be uniquely determined by teir data in oter words, f, under very weak ypoteses on η and f: log η sould be bounded and measurable, and f sould be locally square integrable as a function of time wit values in L Ω. A standard argument using energy estimates and duality Lions, 197; Stolk, 000 establises te existence of a causal weak solution u C 1 R, H0Ω 1 of, satisfying for any α > 0 e[u]t 1 η u u t, t t t Ω + aut, ut 1 dτ e αt τ fτ. 4 4α Te symbol on te RHS of te preceding equation denotes te norm in L Ω: g = g, g. In tis generality, very little can be said about te convergence of finite element approximations, except tat it takes place: te cited references establis existence of weak solutions precisely by demonstrating convergence of a generic Galerkin approximation. Tis report exploits te improved approximation possible wen f is smooter in t. For instance, if f Hloc k R, L Ω, ten te same energy estimates sow tat te weak solution possesses additional regularity: u H k+1 loc R, L Ω. For k = 1, te strong form of te equation olds and it follows tat u C 0 R, H Ω. Te solution is tus smoot enoug for to establis te usual optimal order error estimates for P 1 or Q 1 Galerkin approximations. SEMIDISCRETE GALERKIN METHOD Coose finite dimensional subspaces V H 1 0Ω, define Galerkin approximation u C R, V by d dt ηu, φ + au, φ = f, φ, φ V, 5 plus appropriate initial conditions. Enumerate a basis {v j : 1 j N } of V, write N ux, t = Uj tv j x. 6 j=1
4 4 Symes Define mass, stiffness, and load matrices M, K, F by M i,j = ηv i, v j, K i,j = av i, v j, F j t = f, t, v j. 7 Ten te Galerkin system 5 is equivalent to te system of ODEs M d U dt + K U = F 8 I presume tat te family of subspaces V H 1 0Ω as te standard approximation properties: tere exists C 0, 0 < C C independent of so tat d C I M d C I, and a similar inequality olds for te matrix M 1 in wic η is replaces by 1; a similar elementwise bound olds: for any i, j, 0 M i,j C d ; for all g H Ω H 1 0Ω, tere is v V for wic g v H 1 0 Ω C g H Ω; for all v V, and v H Ω C 1 v H 1 0 Ω, v H 1 0 Ω C 1 v L Ω. Q 1 elements on a regular rectangular grid, for example, ave tese properties REF. Note tat te first of tese properties implies tat u = i U i v i C d U T U u L Ω C d U T U. 9 Also du e[u ] = 1 T M du + U T K U E[U ]. 10 dt dt An argument similar to tat used to establis 4 yields E[u ]t 1 4α dτ e αt τ F τ T M 1 F τ. 11 As remarked at te end of te last section, if f C 1 R, L Ω and f = 0 for t << 0, ten te solution u is smoot enoug tat te usual optimal order error estimate in energy olds: e[u u ]t C f t. 1 See te appendix for a derivation. If f as yet more square-integrable derivatives, corresponding estimates old for iger time derivatives of u u.
5 Lumping 5 MASS LUMPING FOR RECTANGULAR Q 1 ELEMENTS In tis section I suppose tat Ω is a rectangle, wic after translation could be taken to be Ω = [0, l 1 ] [0, l ]... [0, l d ]. Te usual reflection construction makes H 1 0Ω isomorpic to te odd subspace of H 1 perr d, were subscript per signifies multiple periodicity wit multiperiod l 1,..., l d. I will use tis identification witout comment in wat follows. Also note tat te odd multiperiodic extenstion of H 1 0Ω H Ω is a subspace of H perr d. Let = 1,..., d be te vector of cell side lengts, and V H 1 0Ω te corresponding space of Q 1 elements. For a multiindex α = α 1,...α d, let v α be te element of V satisfying v α α 1 1,..., α d d = 1; v α β 1 1,..., β d d = 0, α β N d. Set = {β Z d : β i 1, i = 1,..., d}. Ten M U α = β M α,α+βu α+β = β M α,α+βu α+β U α + M α U α, in wic M α = β M α,α+β is te diagonal α, α entry of te usual row-sum lumped mass matrix. Let Ũ be te solution of te mass-lumped Galerkin system M d Ũ dt + K Ũ = F. 13 For any solution of suc a system, an energy inequality similar to 11 olds: for any α > 0, du Ẽ[u ]t 1 T du M + U T K U t dt dt 1 dτ e αt τ F τ T M 1 F τ. 14 4α A little algebra sows tat te difference between Ũ and te Galerkin trajectory U satisfies [ ] M d dt U Ũ + K U Ũ = d U M α+β α,α+β d Uα. 15 α dt dt β According to 14, Ẽ[U Ũ ] 1 dτ e αt τ 4α
6 6 Symes [ β M, +β d U +β dt ] T [ d U τ M 1 M dt, +β β C d β dτ d U +β dt d U dt d U +β dt ] d U τ. dt, 16 in wic = max 1,..., d and C from ere on will stand for a constant depending on te coefficients and te geometry of te problem but not on. Define te unitary translation operator T β Since v α+β = Tβ v α, d Uα+β dt α d U α dt T β ux = ux + diagβ. v α = α d U α dt : L perr d L perr d by v α β d U α dt v α = T β I u t. Using te bounds on mass matrices described in te last section, obtain d U +β d U d U C d α+β d U α v dt dt dt dt α = C d T β I u t α L Ω L Ω. 17 Combining 16 and 17 gives Ẽ[U Ũ ] C β dτ T β I u t τ L Ω. 18 Suppose tat in addition to te previous ypoteses, f Hloc 3 R, L Ω. Ten 4 applied to te u/ t gives u u t t t C dτ f t τ, wence Ẽ[U Ũ ] C β dτ L Ω T β I u t τ L Ω L Ω + C dτ f t τ As noted earlier, under tese ypoteses on f, u/ t C 0 R, H Ω, and even C 0 R, HperR d, wit u t t t C dτ 3 f t τ 3. H Ω L Ω L Ω
7 Lumping 7 A straigtforward argument based on te trace teorem sows tat T β I u t t L Ω C dτ 3 f t τ 3 L Ω. Combining te last few inequalities, obtain e[ũ u] CẼ[U Ũ ] + e[u u] C dτ 3 f t τ 3 L Ω. 19 Tat is, te error in te lumped mass Galerkin approximation is te same order as te error in te consistent mass approximation. MASS LUMPING FOR P 1 ELEMENTS FOR DIMENSION 1 If te coefficient η as a bit more regularity, ten a more precise estimate tan 19 is possible. In tis section I describe te case d = 1; it seems likely tat similar estimates old for d > 1. In tis section, d = 1, Ω = [0, 1]. Define te element mes x i = i, i = 0,...N + 1 = 1. 0 Coose for V te piecewise linear functions wit nodes {x i }, tat is, te P1 or equally well Q1 element space. Te usual nodal basis is given by v j x = x x j 1, x j 1 x x j x j+1 x, x j x x j+1 0, else Ten K i,i = /, K i,i±1 = 1/, and K i,j = 0 for i j > 1, wit appropriate modifications at te endpoints. Te usual prescription for mass lumping is to replace M by te diagonal matrix M of its row sums. Note tat M is also tridiagonal. Tus it is always possible to write 1 M i, = M i, + ak i, + bd i,. D is te centered difference matrix, te tridiagonal matrix wit A little algebra sows tat D i,i±1 = ±1, D i,i = 0, i =,...N 1. a = M i,i+1 + M i,i 1 P +, b = M i,i+1 M i,i 1 P.
8 8 Symes Te lumped version of te Galerkin ODE system is M d Ũ dt + K Ũ = F. 3 A little more algebra sows tat te difference between its solution Ũ and te Galerkin trajectory U satisfies Since M d dt U Ũ + K U Ũ = P + K d U + P dt D d U dt. 4 M i,i = 1 te first of te following two inequalities follows: were So te lumped energy form 0 ηv i, η min M i,i η max i = 1,..., N, 5 η min I M η max I 6 η min = inf 0 x 1 ηx, η max = sup 0 x 1 ηx. Ẽ U = 1 du T dt M du dt + U T KU is positive definite, as is te consistent energy E U = 1 du T M du dt dt + U T K U. Note tat u = i U i v i E U = e[u], 7 and for suitable positive constants C and C. C E [U] Ẽ [U] C E [U] 8 Te standard energy estimate for causal solutions of system 8 see Appendix is E U t 1 e αt τ F τ T M 1 F τ 9 4α for any α > 0. A similar estimate olds for solutions of te lumped system; in particular, Ẽ U Ũ 1 [ t e αt τ P + K d U T α dt M P 1 + K d U dt
9 Lumping 9 + P D d U T dt M P ] 1 D d U 30 dt To estimate te first term on te RHS of 30, note first tat te equation satisfied by te 3rd time derivative of U, implies via 9 tat d 4 U dt 4 M d5 U dt 5, M d4 U 1 dt 4 4α On te oter and, te second time derivative satisfies M d4 U dt 4 + K d3 U dt 3 = d3 F dt 3, 31 d dτe αt τ 3 F dt, M 1 d3 F 3 dt K d U dt = d F dt, 33 wence d 4 U dt, M d4 U = K d U d F, M 1 K d U d F 4 dt 4 dt dt dt dt 1 K d U dt, M 1 K d U d F dt dt, M 1 d F dt From tis point on, I will use C to denote a constant depending on upper and lower bounds for η and on an upper bound for te time t, but independent of. It follows from te matrix bounds 5 tat P ± is bounded by C in any weigted l norm wit diagonal weigt matrix. Tus P + K d U dt t, M 1 P + K d U dt t C 4 K d U dt t, M 1 K d U dt t C 4 K d U dt t, M 1 K d U dt t d C 4 F dt tt, M 1 d F d 3 dt t F + dt, M 1 d3 F dt 3 To bound te second term in 30, note tat except in te first and last rows, K satisfies K = D + D = D D +, D ± U k = ± 1 U k±1 U k, 34
10 10 Symes wereas Also note tat D = D+ + D. D ± T = D. Denote by Γ te operator norm of P in te weigted l norm wit weigt matrix M 1. According to 5, Γ = O in general. Ten P D d U T M P 1 D d U dt dt Γ D + + D d U T dt M D D d U dt [ Γ D + d U T D + d U + D d U T D ] d U η min dt dt dt dt [ d Γ U T D D + d U d U T + D ] + D d U dt dt dt dt η min 4Γ d U T d U K 8Γ d U E η min dt dt η min dt d C Γ T F d M 1 F dt dt Combine 30, 35, and 36 to yield Ẽ [U Ũ ] C 4 + Γ Define ũ = i Ũ i v i ; ten 3 i=0 e[u ũ ] = E [U Ũ ] CẼ [U Ũ ] C 4 + Γ On te oter and, te rigt and side is C 4 + Γ i=0. d i T F d M 1 i F dt i 3 i=0 3 i f t i i=0. 36 dt i. 37 d i T F d M 1 i F dt i dt i 38 We distinguis two cases. First, suppose tat we know only tat log η is bounded and measurable. Ten Γ = O, and tis is sarp if η as discontinuities. Accordingly, te Galerkin error bound 1 combines wit 38 and 39 to yield 3 e[u ũ ] C i f t i
11 Lumping 11 Second, suppose tat η is Lipsitz continuous. Ten Γ = O 3, and we get instead e[u ũ ] = e[u u ] + O 4, 41 tat is, lumping te mass matrix does not cange te leading beaviour of te error at all. REFERENCES Ciarlet, P., 00, Te finite element metod for elliptic problems: SIAM. Coen, G. C., 001, Higer order numerical metods for transient wave equations: Springer. Li, Z. and K. Ito, 006, Te immersed interface metod: Numerical solutions of pdes involving interfaces and irregular domains: SIAM. Lions, J., 197, Nonomogeneous boundary value problems and applications: Springer Verlag. Owadi, H. and L. Zang, 007, Metric based upscaling: Communications in Pure and Applied Matematics, 60, Owadi, H. and L. Zang, 008, Homogenization of te acoustic wave equation wit a continuum of scales: Computer Metods in Applied Mecanics and Engineering, , Stolk, C., 000, On te modeling and inversion of seismic data: PD tesis, Universiteit Utrect..
MATH745 Fall MATH745 Fall
MATH745 Fall 5 MATH745 Fall 5 INTRODUCTION WELCOME TO MATH 745 TOPICS IN NUMERICAL ANALYSIS Instructor: Dr Bartosz Protas Department of Matematics & Statistics Email: bprotas@mcmasterca Office HH 36, Ext
More informationA h u h = f h. 4.1 The CoarseGrid SystemandtheResidual Equation
Capter Grid Transfer Remark. Contents of tis capter. Consider a grid wit grid size and te corresponding linear system of equations A u = f. Te summary given in Section 3. leads to te idea tat tere migt
More informationConsider a function f we ll specify which assumptions we need to make about it in a minute. Let us reformulate the integral. 1 f(x) dx.
Capter 2 Integrals as sums and derivatives as differences We now switc to te simplest metods for integrating or differentiating a function from its function samples. A careful study of Taylor expansions
More informationPoisson Equation in Sobolev Spaces
Poisson Equation in Sobolev Spaces OcMountain Dayligt Time. 6, 011 Today we discuss te Poisson equation in Sobolev spaces. It s existence, uniqueness, and regularity. Weak Solution. u = f in, u = g on
More informationOrder of Accuracy. ũ h u Ch p, (1)
Order of Accuracy 1 Terminology We consider a numerical approximation of an exact value u. Te approximation depends on a small parameter, wic can be for instance te grid size or time step in a numerical
More information5 Ordinary Differential Equations: Finite Difference Methods for Boundary Problems
5 Ordinary Differential Equations: Finite Difference Metods for Boundary Problems Read sections 10.1, 10.2, 10.4 Review questions 10.1 10.4, 10.8 10.9, 10.13 5.1 Introduction In te previous capters we
More informationTHE STURM-LIOUVILLE-TRANSFORMATION FOR THE SOLUTION OF VECTOR PARTIAL DIFFERENTIAL EQUATIONS. L. Trautmann, R. Rabenstein
Worksop on Transforms and Filter Banks (WTFB),Brandenburg, Germany, Marc 999 THE STURM-LIOUVILLE-TRANSFORMATION FOR THE SOLUTION OF VECTOR PARTIAL DIFFERENTIAL EQUATIONS L. Trautmann, R. Rabenstein Lerstul
More information1 The concept of limits (p.217 p.229, p.242 p.249, p.255 p.256) 1.1 Limits Consider the function determined by the formula 3. x since at this point
MA00 Capter 6 Calculus and Basic Linear Algebra I Limits, Continuity and Differentiability Te concept of its (p.7 p.9, p.4 p.49, p.55 p.56). Limits Consider te function determined by te formula f Note
More informationPolynomial Interpolation
Capter 4 Polynomial Interpolation In tis capter, we consider te important problem of approximatinga function fx, wose values at a set of distinct points x, x, x,, x n are known, by a polynomial P x suc
More information232 Calculus and Structures
3 Calculus and Structures CHAPTER 17 JUSTIFICATION OF THE AREA AND SLOPE METHODS FOR EVALUATING BEAMS Calculus and Structures 33 Copyrigt Capter 17 JUSTIFICATION OF THE AREA AND SLOPE METHODS 17.1 THE
More informationlecture 26: Richardson extrapolation
43 lecture 26: Ricardson extrapolation 35 Ricardson extrapolation, Romberg integration Trougout numerical analysis, one encounters procedures tat apply some simple approximation (eg, linear interpolation)
More informationA Hybrid Mixed Discontinuous Galerkin Finite Element Method for Convection-Diffusion Problems
A Hybrid Mixed Discontinuous Galerkin Finite Element Metod for Convection-Diffusion Problems Herbert Egger Joacim Scöberl We propose and analyse a new finite element metod for convection diffusion problems
More informationA = h w (1) Error Analysis Physics 141
Introduction In all brances of pysical science and engineering one deals constantly wit numbers wic results more or less directly from experimental observations. Experimental observations always ave inaccuracies.
More information1. Introduction. Consider a semilinear parabolic equation in the form
A POSTERIORI ERROR ESTIMATION FOR PARABOLIC PROBLEMS USING ELLIPTIC RECONSTRUCTIONS. I: BACKWARD-EULER AND CRANK-NICOLSON METHODS NATALIA KOPTEVA AND TORSTEN LINSS Abstract. A semilinear second-order parabolic
More informationMA455 Manifolds Solutions 1 May 2008
MA455 Manifolds Solutions 1 May 2008 1. (i) Given real numbers a < b, find a diffeomorpism (a, b) R. Solution: For example first map (a, b) to (0, π/2) and ten map (0, π/2) diffeomorpically to R using
More informationNumerical Experiments Using MATLAB: Superconvergence of Nonconforming Finite Element Approximation for Second-Order Elliptic Problems
Applied Matematics, 06, 7, 74-8 ttp://wwwscirporg/journal/am ISSN Online: 5-7393 ISSN Print: 5-7385 Numerical Experiments Using MATLAB: Superconvergence of Nonconforming Finite Element Approximation for
More informationDiscretization of Multipole Sources in a Finite Difference. Setting for Wave Propagation Problems
Discretization of Multipole Sources in a Finite Difference Setting for Wave Propagation Problems Mario J. Bencomo 1 and William Symes 2 1 Institute for Computational and Experimental Researc in Matematics,
More informationSuperconvergence of energy-conserving discontinuous Galerkin methods for. linear hyperbolic equations. Abstract
Superconvergence of energy-conserving discontinuous Galerkin metods for linear yperbolic equations Yong Liu, Ci-Wang Su and Mengping Zang 3 Abstract In tis paper, we study superconvergence properties of
More informationPolynomial Interpolation
Capter 4 Polynomial Interpolation In tis capter, we consider te important problem of approximating a function f(x, wose values at a set of distinct points x, x, x 2,,x n are known, by a polynomial P (x
More informationPreconditioning in H(div) and Applications
1 Preconditioning in H(div) and Applications Douglas N. Arnold 1, Ricard S. Falk 2 and Ragnar Winter 3 4 Abstract. Summarizing te work of [AFW97], we sow ow to construct preconditioners using domain decomposition
More informationChapter 5 FINITE DIFFERENCE METHOD (FDM)
MEE7 Computer Modeling Tecniques in Engineering Capter 5 FINITE DIFFERENCE METHOD (FDM) 5. Introduction to FDM Te finite difference tecniques are based upon approximations wic permit replacing differential
More informationHow to Find the Derivative of a Function: Calculus 1
Introduction How to Find te Derivative of a Function: Calculus 1 Calculus is not an easy matematics course Te fact tat you ave enrolled in suc a difficult subject indicates tat you are interested in te
More informationQuasiperiodic phenomena in the Van der Pol - Mathieu equation
Quasiperiodic penomena in te Van der Pol - Matieu equation F. Veerman and F. Verulst Department of Matematics, Utrect University P.O. Box 80.010, 3508 TA Utrect Te Neterlands April 8, 009 Abstract Te Van
More informationNUMERICAL DIFFERENTIATION. James T. Smith San Francisco State University. In calculus classes, you compute derivatives algebraically: for example,
NUMERICAL DIFFERENTIATION James T Smit San Francisco State University In calculus classes, you compute derivatives algebraically: for example, f( x) = x + x f ( x) = x x Tis tecnique requires your knowing
More information3 Parabolic Differential Equations
3 Parabolic Differential Equations 3.1 Classical solutions We consider existence and uniqueness results for initial-boundary value problems for te linear eat equation in Q := Ω (, T ), were Ω is a bounded
More informationarxiv: v2 [math.na] 5 Jul 2017
Trace Finite Element Metods for PDEs on Surfaces Maxim A. Olsanskii and Arnold Reusken arxiv:1612.00054v2 [mat.na] 5 Jul 2017 Abstract In tis paper we consider a class of unfitted finite element metods
More informationDecay of solutions of wave equations with memory
Proceedings of te 14t International Conference on Computational and Matematical Metods in Science and Engineering, CMMSE 14 3 7July, 14. Decay of solutions of wave equations wit memory J. A. Ferreira 1,
More informationSome Error Estimates for the Finite Volume Element Method for a Parabolic Problem
Computational Metods in Applied Matematics Vol. 13 (213), No. 3, pp. 251 279 c 213 Institute of Matematics, NAS of Belarus Doi: 1.1515/cmam-212-6 Some Error Estimates for te Finite Volume Element Metod
More informationA Finite Element Primer
A Finite Element Primer David J. Silvester Scool of Matematics, University of Mancester d.silvester@mancester.ac.uk. Version.3 updated 4 October Contents A Model Diffusion Problem.................... x.
More information1 Calculus. 1.1 Gradients and the Derivative. Q f(x+h) f(x)
Calculus. Gradients and te Derivative Q f(x+) δy P T δx R f(x) 0 x x+ Let P (x, f(x)) and Q(x+, f(x+)) denote two points on te curve of te function y = f(x) and let R denote te point of intersection of
More informationGeneric maximum nullity of a graph
Generic maximum nullity of a grap Leslie Hogben Bryan Sader Marc 5, 2008 Abstract For a grap G of order n, te maximum nullity of G is defined to be te largest possible nullity over all real symmetric n
More informationNumerical Differentiation
Numerical Differentiation Finite Difference Formulas for te first derivative (Using Taylor Expansion tecnique) (section 8.3.) Suppose tat f() = g() is a function of te variable, and tat as 0 te function
More informationEffect of Numerical Integration on Meshless Methods
Effect of Numerical Integration on Mesless Metods Ivo Babuška Uday Banerjee Jon E. Osborn Qingui Zang Abstract In tis paper, we present te effect of numerical integration on mesless metods wit sape functions
More informationDifferentiation in higher dimensions
Capter 2 Differentiation in iger dimensions 2.1 Te Total Derivative Recall tat if f : R R is a 1-variable function, and a R, we say tat f is differentiable at x = a if and only if te ratio f(a+) f(a) tends
More information2.8 The Derivative as a Function
.8 Te Derivative as a Function Typically, we can find te derivative of a function f at many points of its domain: Definition. Suppose tat f is a function wic is differentiable at every point of an open
More informationSolving Continuous Linear Least-Squares Problems by Iterated Projection
Solving Continuous Linear Least-Squares Problems by Iterated Projection by Ral Juengling Department o Computer Science, Portland State University PO Box 75 Portland, OR 977 USA Email: juenglin@cs.pdx.edu
More informationMIXED DISCONTINUOUS GALERKIN APPROXIMATION OF THE MAXWELL OPERATOR. SIAM J. Numer. Anal., Vol. 42 (2004), pp
MIXED DISCONTINUOUS GALERIN APPROXIMATION OF THE MAXWELL OPERATOR PAUL HOUSTON, ILARIA PERUGIA, AND DOMINI SCHÖTZAU SIAM J. Numer. Anal., Vol. 4 (004), pp. 434 459 Abstract. We introduce and analyze a
More informationClick here to see an animation of the derivative
Differentiation Massoud Malek Derivative Te concept of derivative is at te core of Calculus; It is a very powerful tool for understanding te beavior of matematical functions. It allows us to optimize functions,
More informationVolume 29, Issue 3. Existence of competitive equilibrium in economies with multi-member households
Volume 29, Issue 3 Existence of competitive equilibrium in economies wit multi-member ouseolds Noriisa Sato Graduate Scool of Economics, Waseda University Abstract Tis paper focuses on te existence of
More informationFunction Composition and Chain Rules
Function Composition and s James K. Peterson Department of Biological Sciences and Department of Matematical Sciences Clemson University Marc 8, 2017 Outline 1 Function Composition and Continuity 2 Function
More informationInf sup testing of upwind methods
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING Int. J. Numer. Met. Engng 000; 48:745 760 Inf sup testing of upwind metods Klaus-Jurgen Bate 1; ;, Dena Hendriana 1, Franco Brezzi and Giancarlo
More informationA Mixed-Hybrid-Discontinuous Galerkin Finite Element Method for Convection-Diffusion Problems
A Mixed-Hybrid-Discontinuous Galerkin Finite Element Metod for Convection-Diffusion Problems Herbert Egger Joacim Scöberl We propose and analyse a new finite element metod for convection diffusion problems
More informationarxiv: v1 [math.na] 7 Mar 2019
Local Fourier analysis for mixed finite-element metods for te Stokes equations Yunui He a,, Scott P. MacLaclan a a Department of Matematics and Statistics, Memorial University of Newfoundland, St. Jon
More informationParameter Fitted Scheme for Singularly Perturbed Delay Differential Equations
International Journal of Applied Science and Engineering 2013. 11, 4: 361-373 Parameter Fitted Sceme for Singularly Perturbed Delay Differential Equations Awoke Andargiea* and Y. N. Reddyb a b Department
More informationFinite Difference Method
Capter 8 Finite Difference Metod 81 2nd order linear pde in two variables General 2nd order linear pde in two variables is given in te following form: L[u] = Au xx +2Bu xy +Cu yy +Du x +Eu y +Fu = G According
More informationERROR BOUNDS FOR THE METHODS OF GLIMM, GODUNOV AND LEVEQUE BRADLEY J. LUCIER*
EO BOUNDS FO THE METHODS OF GLIMM, GODUNOV AND LEVEQUE BADLEY J. LUCIE* Abstract. Te expected error in L ) attimet for Glimm s sceme wen applied to a scalar conservation law is bounded by + 2 ) ) /2 T
More informationError estimates for a semi-implicit fully discrete finite element scheme for the mean curvature flow of graphs
Interfaces and Free Boundaries 2, 2000 34 359 Error estimates for a semi-implicit fully discrete finite element sceme for te mean curvature flow of graps KLAUS DECKELNICK Scool of Matematical Sciences,
More informationMVT and Rolle s Theorem
AP Calculus CHAPTER 4 WORKSHEET APPLICATIONS OF DIFFERENTIATION MVT and Rolle s Teorem Name Seat # Date UNLESS INDICATED, DO NOT USE YOUR CALCULATOR FOR ANY OF THESE QUESTIONS In problems 1 and, state
More informationarxiv: v1 [math.na] 28 Apr 2017
THE SCOTT-VOGELIUS FINITE ELEMENTS REVISITED JOHNNY GUZMÁN AND L RIDGWAY SCOTT arxiv:170500020v1 [matna] 28 Apr 2017 Abstract We prove tat te Scott-Vogelius finite elements are inf-sup stable on sape-regular
More informationJournal of Computational and Applied Mathematics
Journal of Computational and Applied Matematics 94 (6) 75 96 Contents lists available at ScienceDirect Journal of Computational and Applied Matematics journal omepage: www.elsevier.com/locate/cam Smootness-Increasing
More informationA First-Order System Approach for Diffusion Equation. I. Second-Order Residual-Distribution Schemes
A First-Order System Approac for Diffusion Equation. I. Second-Order Residual-Distribution Scemes Hiroaki Nisikawa W. M. Keck Foundation Laboratory for Computational Fluid Dynamics, Department of Aerospace
More informationSECTION 3.2: DERIVATIVE FUNCTIONS and DIFFERENTIABILITY
(Section 3.2: Derivative Functions and Differentiability) 3.2.1 SECTION 3.2: DERIVATIVE FUNCTIONS and DIFFERENTIABILITY LEARNING OBJECTIVES Know, understand, and apply te Limit Definition of te Derivative
More informationFinite Difference Methods Assignments
Finite Difference Metods Assignments Anders Söberg and Aay Saxena, Micael Tuné, and Maria Westermarck Revised: Jarmo Rantakokko June 6, 1999 Teknisk databeandling Assignment 1: A one-dimensional eat equation
More information3.4 Worksheet: Proof of the Chain Rule NAME
Mat 1170 3.4 Workseet: Proof of te Cain Rule NAME Te Cain Rule So far we are able to differentiate all types of functions. For example: polynomials, rational, root, and trigonometric functions. We are
More informationarxiv: v1 [math.na] 12 Mar 2018
ON PRESSURE ESTIMATES FOR THE NAVIER-STOKES EQUATIONS J A FIORDILINO arxiv:180304366v1 [matna 12 Mar 2018 Abstract Tis paper presents a simple, general tecnique to prove finite element metod (FEM) pressure
More informationPart VIII, Chapter 39. Fluctuation-based stabilization Model problem
Part VIII, Capter 39 Fluctuation-based stabilization Tis capter presents a unified analysis of recent stabilization tecniques for te standard Galerkin approximation of first-order PDEs using H 1 - conforming
More informationOn convergence of the immersed boundary method for elliptic interface problems
On convergence of te immersed boundary metod for elliptic interface problems Zilin Li January 26, 2012 Abstract Peskin s Immersed Boundary (IB) metod is one of te most popular numerical metods for many
More informationAnalysis of A Continuous Finite Element Method for H(curl, div)-elliptic Interface Problem
Analysis of A Continuous inite Element Metod for Hcurl, div)-elliptic Interface Problem Huoyuan Duan, Ping Lin, and Roger C. E. Tan Abstract In tis paper, we develop a continuous finite element metod for
More informationarxiv: v1 [math.na] 17 Jul 2014
Div First-Order System LL* FOSLL* for Second-Order Elliptic Partial Differential Equations Ziqiang Cai Rob Falgout Sun Zang arxiv:1407.4558v1 [mat.na] 17 Jul 2014 February 13, 2018 Abstract. Te first-order
More informationVariational Localizations of the Dual Weighted Residual Estimator
Publised in Journal for Computational and Applied Matematics, pp. 192-208, 2015 Variational Localizations of te Dual Weigted Residual Estimator Tomas Ricter Tomas Wick Te dual weigted residual metod (DWR)
More informationPOLYNOMIAL AND SPLINE ESTIMATORS OF THE DISTRIBUTION FUNCTION WITH PRESCRIBED ACCURACY
APPLICATIONES MATHEMATICAE 36, (29), pp. 2 Zbigniew Ciesielski (Sopot) Ryszard Zieliński (Warszawa) POLYNOMIAL AND SPLINE ESTIMATORS OF THE DISTRIBUTION FUNCTION WITH PRESCRIBED ACCURACY Abstract. Dvoretzky
More informationdu(l) 5 dl = U(0) = 1 and 1.) Substitute for U an unspecified trial function into governing equation, i.e. dx + = =
Consider an ODE of te form: Finite Element Metod du fu g d + wit te following Boundary Conditions: U(0) and du(l) 5 dl.) Substitute for U an unspecified trial function into governing equation, i.e. ^ U
More informationPrecalculus Test 2 Practice Questions Page 1. Note: You can expect other types of questions on the test than the ones presented here!
Precalculus Test 2 Practice Questions Page Note: You can expect oter types of questions on te test tan te ones presented ere! Questions Example. Find te vertex of te quadratic f(x) = 4x 2 x. Example 2.
More informationMath 1210 Midterm 1 January 31st, 2014
Mat 110 Midterm 1 January 1st, 01 Tis exam consists of sections, A and B. Section A is conceptual, wereas section B is more computational. Te value of every question is indicated at te beginning of it.
More informationExam 1 Review Solutions
Exam Review Solutions Please also review te old quizzes, and be sure tat you understand te omework problems. General notes: () Always give an algebraic reason for your answer (graps are not sufficient),
More informationContinuity and Differentiability of the Trigonometric Functions
[Te basis for te following work will be te definition of te trigonometric functions as ratios of te sides of a triangle inscribed in a circle; in particular, te sine of an angle will be defined to be te
More information, meant to remind us of the definition of f (x) as the limit of difference quotients: = lim
Mat 132 Differentiation Formulas Stewart 2.3 So far, we ave seen ow various real-world problems rate of cange and geometric problems tangent lines lead to derivatives. In tis section, we will see ow to
More informationLEAST-SQUARES FINITE ELEMENT APPROXIMATIONS TO SOLUTIONS OF INTERFACE PROBLEMS
SIAM J. NUMER. ANAL. c 998 Society for Industrial Applied Matematics Vol. 35, No., pp. 393 405, February 998 00 LEAST-SQUARES FINITE ELEMENT APPROXIMATIONS TO SOLUTIONS OF INTERFACE PROBLEMS YANZHAO CAO
More informationREVIEW LAB ANSWER KEY
REVIEW LAB ANSWER KEY. Witout using SN, find te derivative of eac of te following (you do not need to simplify your answers): a. f x 3x 3 5x x 6 f x 3 3x 5 x 0 b. g x 4 x x x notice te trick ere! x x g
More information1. Introduction. We consider the model problem: seeking an unknown function u satisfying
A DISCONTINUOUS LEAST-SQUARES FINITE ELEMENT METHOD FOR SECOND ORDER ELLIPTIC EQUATIONS XIU YE AND SHANGYOU ZHANG Abstract In tis paper, a discontinuous least-squares (DLS) finite element metod is introduced
More informationContinuity and Differentiability Worksheet
Continuity and Differentiability Workseet (Be sure tat you can also do te grapical eercises from te tet- Tese were not included below! Typical problems are like problems -3, p. 6; -3, p. 7; 33-34, p. 7;
More informationMath 102 TEST CHAPTERS 3 & 4 Solutions & Comments Fall 2006
Mat 102 TEST CHAPTERS 3 & 4 Solutions & Comments Fall 2006 f(x+) f(x) 10 1. For f(x) = x 2 + 2x 5, find ))))))))) and simplify completely. NOTE: **f(x+) is NOT f(x)+! f(x+) f(x) (x+) 2 + 2(x+) 5 ( x 2
More informationCHAPTER 4. Elliptic PDEs
CHAPTER 4 Elliptic PDEs One of te main advantages of extending te class of solutions of a PDE from classical solutions wit continuous derivatives to weak solutions wit weak derivatives is tat it is easier
More informationlim 1 lim 4 Precalculus Notes: Unit 10 Concepts of Calculus
Syllabus Objectives: 1.1 Te student will understand and apply te concept of te limit of a function at given values of te domain. 1. Te student will find te limit of a function at given values of te domain.
More informationPre-Calculus Review Preemptive Strike
Pre-Calculus Review Preemptive Strike Attaced are some notes and one assignment wit tree parts. Tese are due on te day tat we start te pre-calculus review. I strongly suggest reading troug te notes torougly
More informationChapter 2 Limits and Continuity
4 Section. Capter Limits and Continuity Section. Rates of Cange and Limits (pp. 6) Quick Review.. f () ( ) () 4 0. f () 4( ) 4. f () sin sin 0 4. f (). 4 4 4 6. c c c 7. 8. c d d c d d c d c 9. 8 ( )(
More informationThe Convergence of a Central-Difference Discretization of Rudin-Osher-Fatemi Model for Image Denoising
Te Convergence of a Central-Difference Discretization of Rudin-Oser-Fatemi Model for Image Denoising Ming-Jun Lai 1, Bradley Lucier 2, and Jingyue Wang 3 1 University of Georgia, Atens GA 30602, USA mjlai@mat.uga.edu
More informationMATH1151 Calculus Test S1 v2a
MATH5 Calculus Test 8 S va January 8, 5 Tese solutions were written and typed up by Brendan Trin Please be etical wit tis resource It is for te use of MatSOC members, so do not repost it on oter forums
More informationERROR BOUNDS FOR FINITE-DIFFERENCE METHODS FOR RUDIN OSHER FATEMI IMAGE SMOOTHING
ERROR BOUNDS FOR FINITE-DIFFERENCE METHODS FOR RUDIN OSHER FATEMI IMAGE SMOOTHING JINGYUE WANG AND BRADLEY J. LUCIER Abstract. We bound te difference between te solution to te continuous Rudin Oser Fatemi
More informationBob Brown Math 251 Calculus 1 Chapter 3, Section 1 Completed 1 CCBC Dundalk
Bob Brown Mat 251 Calculus 1 Capter 3, Section 1 Completed 1 Te Tangent Line Problem Te idea of a tangent line first arises in geometry in te context of a circle. But before we jump into a discussion of
More information1. Consider the trigonometric function f(t) whose graph is shown below. Write down a possible formula for f(t).
. Consider te trigonometric function f(t) wose grap is sown below. Write down a possible formula for f(t). Tis function appears to be an odd, periodic function tat as been sifted upwards, so we will use
More informationMANY scientific and engineering problems can be
A Domain Decomposition Metod using Elliptical Arc Artificial Boundary for Exterior Problems Yajun Cen, and Qikui Du Abstract In tis paper, a Diriclet-Neumann alternating metod using elliptical arc artificial
More informationChapter 4: Numerical Methods for Common Mathematical Problems
1 Capter 4: Numerical Metods for Common Matematical Problems Interpolation Problem: Suppose we ave data defined at a discrete set of points (x i, y i ), i = 0, 1,..., N. Often it is useful to ave a smoot
More informationPolynomial Functions. Linear Functions. Precalculus: Linear and Quadratic Functions
Concepts: definition of polynomial functions, linear functions tree representations), transformation of y = x to get y = mx + b, quadratic functions axis of symmetry, vertex, x-intercepts), transformations
More informationOSCILLATION OF SOLUTIONS TO NON-LINEAR DIFFERENCE EQUATIONS WITH SEVERAL ADVANCED ARGUMENTS. Sandra Pinelas and Julio G. Dix
Opuscula Mat. 37, no. 6 (2017), 887 898 ttp://dx.doi.org/10.7494/opmat.2017.37.6.887 Opuscula Matematica OSCILLATION OF SOLUTIONS TO NON-LINEAR DIFFERENCE EQUATIONS WITH SEVERAL ADVANCED ARGUMENTS Sandra
More informationThe Priestley-Chao Estimator
Te Priestley-Cao Estimator In tis section we will consider te Pristley-Cao estimator of te unknown regression function. It is assumed tat we ave a sample of observations (Y i, x i ), i = 1,..., n wic are
More informationContinuity. Example 1
Continuity MATH 1003 Calculus and Linear Algebra (Lecture 13.5) Maoseng Xiong Department of Matematics, HKUST A function f : (a, b) R is continuous at a point c (a, b) if 1. x c f (x) exists, 2. f (c)
More informationSECTION 1.10: DIFFERENCE QUOTIENTS LEARNING OBJECTIVES
(Section.0: Difference Quotients).0. SECTION.0: DIFFERENCE QUOTIENTS LEARNING OBJECTIVES Define average rate of cange (and average velocity) algebraically and grapically. Be able to identify, construct,
More informationThe derivative function
Roberto s Notes on Differential Calculus Capter : Definition of derivative Section Te derivative function Wat you need to know already: f is at a point on its grap and ow to compute it. Wat te derivative
More informationIntroduction to Derivatives
Introduction to Derivatives 5-Minute Review: Instantaneous Rates and Tangent Slope Recall te analogy tat we developed earlier First we saw tat te secant slope of te line troug te two points (a, f (a))
More information1.72, Groundwater Hydrology Prof. Charles Harvey Lecture Packet #9: Numerical Modeling of Groundwater Flow
1.7, Groundwater Hydrology Prof. Carles Harvey Lecture Packet #9: Numerical Modeling of Groundwater Flow Simulation: Te prediction of quantities of interest (dependent variables) based upon an equation
More informationSolution. Solution. f (x) = (cos x)2 cos(2x) 2 sin(2x) 2 cos x ( sin x) (cos x) 4. f (π/4) = ( 2/2) ( 2/2) ( 2/2) ( 2/2) 4.
December 09, 20 Calculus PracticeTest s Name: (4 points) Find te absolute extrema of f(x) = x 3 0 on te interval [0, 4] Te derivative of f(x) is f (x) = 3x 2, wic is zero only at x = 0 Tus we only need
More informationMath 34A Practice Final Solutions Fall 2007
Mat 34A Practice Final Solutions Fall 007 Problem Find te derivatives of te following functions:. f(x) = 3x + e 3x. f(x) = x + x 3. f(x) = (x + a) 4. Is te function 3t 4t t 3 increasing or decreasing wen
More information2.11 That s So Derivative
2.11 Tat s So Derivative Introduction to Differential Calculus Just as one defines instantaneous velocity in terms of average velocity, we now define te instantaneous rate of cange of a function at a point
More informationSection 2.7 Derivatives and Rates of Change Part II Section 2.8 The Derivative as a Function. at the point a, to be. = at time t = a is
Mat 180 www.timetodare.com Section.7 Derivatives and Rates of Cange Part II Section.8 Te Derivative as a Function Derivatives ( ) In te previous section we defined te slope of te tangent to a curve wit
More informationCS522 - Partial Di erential Equations
CS5 - Partial Di erential Equations Tibor Jánosi April 5, 5 Numerical Di erentiation In principle, di erentiation is a simple operation. Indeed, given a function speci ed as a closed-form formula, its
More informationFINITE ELEMENT EXTERIOR CALCULUS FOR PARABOLIC EVOLUTION PROBLEMS ON RIEMANNIAN HYPERSURFACES
FINITE ELEMENT EXTERIOR CALCULUS FOR PARABOLIC EVOLUTION PROBLEMS ON RIEMANNIAN HYPERSURFACES MICHAEL HOLST AND CHRIS TIEE ABSTRACT. Over te last ten years, te Finite Element Exterior Calculus (FEEC) as
More informationERROR ESTIMATES FOR THE DISCONTINUOUS GALERKIN METHODS FOR PARABOLIC EQUATIONS. 1. Introduction. We consider the parabolic PDE of the form,
ERROR ESTIMATES FOR THE DISCONTINUOUS GALERKIN METHODS FOR PARABOLIC EQUATIONS K. CHRYSAFINOS AND NOEL. J. WALKINGTON Abstract. We analyze te classical discontinuous Galerkin metod for a general parabolic
More informationLecture 15. Interpolation II. 2 Piecewise polynomial interpolation Hermite splines
Lecture 5 Interpolation II Introduction In te previous lecture we focused primarily on polynomial interpolation of a set of n points. A difficulty we observed is tat wen n is large, our polynomial as to
More informationA method of Lagrange Galerkin of second order in time. Une méthode de Lagrange Galerkin d ordre deux en temps
A metod of Lagrange Galerkin of second order in time Une métode de Lagrange Galerkin d ordre deux en temps Jocelyn Étienne a a DAMTP, University of Cambridge, Wilberforce Road, Cambridge CB3 0WA, Great-Britain.
More information