Preconditioning in H(div) and Applications
|
|
- Jade Jenkins
- 6 years ago
- Views:
Transcription
1 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 and multigrid tecniques for te system of linear algebraic equations wic arises from te finite element discretization of boundary value problems associated to te differential operator I grad div. Tese preconditioners are sown to be spectrally equivalent to te inverse of te operator and tus may be used to precondition iterative metods so tat any given error reduction may be acieved in a finite number of iterations independent of te mes discretization. We describe applications of tese results to te efficient solution of mixed and least squares finite element approximations of elliptic boundary value problems. 1.1 Introduction Tis paper summarizes te work of [AFW97], in wic we consider te solution of te system of linear algebraic equations wic arises from te finite element discretization of boundary value problems in two space dimensions for te differential operator I grad div. Te natural setting for te weak formulation of suc problems is te space: H(div) = { u L 2 (Ω) div u L 2 (Ω) }. Let (, ) denote te L 2 (Ω) inner product of bot scalar and vector-valued functions and J(u, v) := (u, v) + (div u, div v) 1 Department of Matematics, Penn State, University Park, PA dnamat.psu.edu 2 Dept. of Matematics, Rutgers University, New Brunswick, NJ falk@mat.rutgers.edu 3 Department of Informatics, University of Oslo, N-0316 Oslo, Norway. ragnar@ifi.uio.no 4 Te first autor was supported by NSF grants DMS and DMS and by te Institute for Matematics and its Applications. Te second autor was supported by NSF grant DMS Te tird autor was supported by Te Norwegian Researc Council under grants /431 and STP DD9 Proceedings Editor Petter Bjørstad, Magne Espedal and David Keyes c 1996 Jon Wiley & Sons Ltd.
2 2 D. N. ARNOLD, R. S. FALK, AND R. WINTHER denote te innerproduct on H(div). If f L 2 (Ω), te weak formulation is to find u H(div) suc tat for all v H(div), J(u, v) = (f, v). Tis corresponds to te boundary value problem (I grad div)u = f in Ω, div u = 0 on Ω. Note tat if u is a gradient, ten (I grad div)u = u+u, wile if u is a curl, ten (I grad div)u = u. A simple situation in wic te operator I grad div arises occurs in te computation of u = grad p, were p is te solution of te Diriclet problem p + p = g in Ω, p = 0 on Ω. Ten u H(div) satisfies J(u, v) = (g, div v) for all v H(div). Given a finite element subspace V of H(div), te natural finite element approximation sceme is: Find u V suc tat J(u, v ) = (f, v ) for all v V. We sall consider te case wen V consists of te Raviart Tomas space of index k 0, i.e., functions wic on eac triangle are of te form v(x, y) = p(x, y) + (x, y)q(x, y), p P k P k, q P k, (were P k denotes te polynomials of degree k) and for wic v n is continuous from triangle to triangle. Te goal is to find an efficient procedure for solving te discrete linear system corresponding to tis discretization, wic we write as J u = f. Denoting te eigenvalues of J by σ(j ), since te spectral condition number κ(j ) := max σ(j ) min σ(j ) of te operator J is O( 2 ), we will clearly need to precondition any standard iterative sceme if we want te number of iterations needed to acieve a given accuracy to be independent of. 1.2 Preconditioning in te abstract Let X L 2 be a finite dimensional normed vectorspace. We identify X and X as sets, but put te dual norm on te latter (dual wit respect to te L 2 inner product). Let A : X X be an L 2 -symmetric linear isomorpism. We suppose tat X is endowed wit an appropriate (energy) norm, i.e., we suppose tat A L(X,X ), A 1 L(X,X ) = O(1).
3 PRECONDITIONING IN H(DIV) AND APPLICATIONS 3 Given f X, we wis to solve A x = f by applying a standard iterative metod suc as CG or MINRES to te equation B A x = B f, were B : X X is an L 2 -symmetric, positive definite preconditioner. Our goal is to define B so tat te action of B is easily computable and κ(b A ) is bounded uniformly wit respect to. Since max σ(b A ) B A L(X,X ) A L(X,X ) B L(X,X ) and 1 min σ(b A ) (B A ) 1 L(X,X ) A 1 L(X,X ) B 1 L(X,X ) B is an effective preconditioner if B L(X,X ), B 1 L(X,X ) = O(1). In oter words, B is an effective preconditioner if it as te same mapping properties as A 1. Note tat te energy norm, and not te detailed structure of A, determine tese properties. Tus to solve te problem J u = f, we need to construct an efficiently computable operator K : V V for wic K L(V,V ), K 1 L(V,V ) = O(1). We will sow ow tis can be done using domain decomposition and multigrid tecniques. 1.3 Applications We are interested in te operator I grad div not for its own sake, but for its appearance in several important problems. Besides te example mentioned in te introduction, we will restrict our attention to two problems: te least squares formulation and te mixed formulation of second order scalar elliptic problems. Oter examples are discussed in [AFW97]. We first discuss te least squares variational principle. Consider te elliptic boundary value problem div(a grad p) = g in Ω, p = 0 on Ω, were te coefficient matrix A is assumed measurable, bounded, symmetric, and uniformly positive definite on Ω. Introducing u = A grad p leads to te first order system u A grad p = 0 in Ω, div u = g in Ω, p = 0 on Ω. Te least squares variational principle caracterizes te solution (u, p) as te minimizer of te functional v A grad q 2 + div v g 2
4 4 D. N. ARNOLD, R. S. FALK, AND R. WINTHER over te space H(div) H 1, were denotes te L 2 (Ω) norm and H 1 denotes te subspace of functions in H 1 (Ω) wic vanis on te boundary of Ω. Equivalently, we ave te weak formulation were B(u, p; v, q) = (g, div v) for all (v, q) H(div) H 1, B(u, p; v, q) = (u A grad p, v A grad q) + (div u, div v). To discretize te least squares formulation, we let X = V W be a finitedimensional subspace of H(div) H 1. Ten x := (u, p ) is te minimizer over X of v A grad q 2 + div v g 2, or in weak form, B(u, p ; v, q) = (g, div v) for all (v, q) X. Defining A : X X by (A x, y) = B(x, y) and f X by (f, (v, q)) = (g, div v), we may rewrite our problem as A x = f. Te key to te convergence teory for te least squares metod is te following teorem (cf. Pelivanov, Carey, Lazarov [PCL94] and Cai, Lazarov, Manteuffel, and McCormick [CLMM94]). Teorem 1.1 Te bilinear form B is an inner product on H(div) H 1 equivalent to te usual one. A direct consequence of te teorem is tat A : X X is symmetric, positive definite and satisfies A L(X,X ), A 1 L(X,X ) = O(1). Tus we need a preconditioner wit te opposite mapping properties. Since X = V W, we can coose a block diagonal preconditioner ( ) K 0 B =, 0 M were K is a good preconditioner in H(div), i.e., it maps like J 1 : V V, and M is a good preconditioner in H 1, i.e., it maps like 1 : W W. Hence we conclude tat a good preconditioner for te discrete least squares system is obtained using an H(div) preconditioner for te vector variable and a standard H 1 preconditioner for te scalar variable. We next consider a mixed variational formulation of tis boundary value problem. Te mixed variational principle caracterizes (u, p) as a saddle point of over H(div) L 2, or, in weak form, 1 2 (A 1 v, v) + (q, div v) (g, q), (A 1 u, v) + (p, div v) = 0 for all v H(div),
5 PRECONDITIONING IN H(DIV) AND APPLICATIONS 5 (div u, q) = (g, q) for all q L 2. Coosing X = V S H(div) L 2, we can define a discrete solution x = (u, p ) V S by restricting eiter te variational or weak formulation. Tis may be written A x = f, wit A : X X L 2 -symmetric but indefinite, since A as te form A = ( a b b t 0 ). Te convergence of tis metod depends on te coice of V ypoteses for te convergence analysis are te Brezzi conditions: and S. Te key (A 1 v, v) γ 1 v H(div) for all v V wit div v S, inf q S (q, div v) sup γ 2. v V q v H(div) Tese conditions are satisfied if V is te Raviart Tomas space of index k and S te space of (discontinuous) piecewise polynomials of degree k. Brezzi s teorem states tat if bot ypoteses are satisfied, ten A is an isomorpism and A 1 L(X,X ) may be bounded in terms of te γ i. We tus base our coice of B on te discrete version of te isomorpism ( ) A grad : H(div) L 2 H(div) L 2. div 0 We again use a simple block-diagonal preconditioner, wic tis time takes te form ( ) K 0 B =, 0 I were I is te identity on S and again K is a good preconditioner in H(div), i.e., it maps like J 1 : V V. We remark tat most oter work on preconditioning suc mixed metods uses te alternate isomorpism ( ) A grad : L 2 div 0 H 1 L 2 H 1, wic leads to a different (and less natural) coice of preconditioner. 1.4 An additive Scwarz preconditioner for J We let T H = {Ω n } N n=0, denote te coarse mes and T a refinement (te fine mes). We let {Ω n} N n=1 be an overlapping covering aligned wit te fine mes suc tat Ω n Ω n. We make te standard assumption of sufficient but bounded overlap. Let V n denote te Raviart Tomas space approximating H(div, Ω n) wit te boundary condition v n = 0 on Ω j \ Ω. Let V 0 denote te Raviart Tomas approximation to H(div, Ω) using te coarse mes.
6 6 D. N. ARNOLD, R. S. FALK, AND R. WINTHER Given f V, define u n V n by J(u n, v) = (f, v) for all v V n. Te additive Scwarz preconditioner is ten defined by K f := N n=0 u n. Our main result for tis domain decomposition preconditioner is te following teorem (cf. [AFW97] for te proof). Teorem 1.2 Tere exists a constant c independent of bot and H for wic κ(k J ) c. Following te teoretical framework of Dryja Widlund [DW90] or Xu [Xu92], a critical step of te proof is te following decomposition lemma. Lemma 1.1 For all v V, tere exist v n V n wit v = N n=0 v n and N v n 2 H(div) c v H(div). n=0 Te standard proof uses a partition of unity {θ n } N n=1 and takes v 0 V 0 a suitable approximation of v and v n = Π [θ n (v v 0 )] wit Π a suitable local projection into V. Te analysis leads to te following estimates. div v n c div[θ n (v v 0 )] c grad θ n L v v 0 + θ n L div(v v 0 ) ch 1 v v 0 + div(v v 0 ). In te standard elliptic case we bound te first term using v v 0 CH v 1. However it is not true tat v v 0 CH v H(div), so tis approac fails. We are able to get around tis problem by using a discrete Helmoltz decomposition, wic we now describe. Let V denote te Raviart Tomas space of index k, S te space of piecewise polynomials of degree k, and W te space of C 0 piecewise polynomials of degree k + 1. Ten we ave te following discrete Helmoltz decomposition. V = curl W grad S, were grad : S V is defined by (grad s, v) = (s, div v). Returning to te decomposition lemma, we write v = curl w + grad s and observe tat v 2 H(div) = curl w 2 + grad s 2 H(div). We ten decompose eac term separately. Since curl w H(div) w 1, we can use te standard decomposition lemma on w to write n n w = w j, w j 2 1 c w 2 1. j=0 j=0 Taking curls gives us te desired result on te curl w term. For v = grad s and v 0 = grad H s 0, were (s 0, v 0 ) is te mixed metod approximation to (s, v) in te space S 0 V 0, we can prove using standard results from te teory of mixed finite element approximations tat v v 0 CH v H(div),
7 PRECONDITIONING IN H(DIV) AND APPLICATIONS 7 and conclude te proof. Te key is tat altoug te above estimate does not old for all v V, it does old wen v = grad s. 1.5 V-cycle preconditioner We consider a nested sequence of meses, T 1, T 2,..., T N, and let V n be te Raviart- Tomas space of some fixed order subordinate to te mes T n. Tis gives a nested sequence of spaces V 1 V 2 V N = V and corresponding operators J n : V n V n. We also require smooters R n : V n V n wic we discuss below and te H(div)- projection operators P n : H(div) V n. Multigrid ten defines K n : V n V n recursively starting wit K 1 = J 1 1. We sall make use of te following multigrid convergence result. Teorem 1.3 Suppose tat for eac n = 1, 2,..., N te smooter R n is L 2 - symmetric and positive semi-definite and satisfies te conditions J([I R n J n ]v, v) 0 (R 1 n [I P n 1 ]v, [I P n 1 ]v) αj([i P n 1 ]v, [I P n 1 ]v). Ten tere exists a constant C independent of and N suc tat te eigenvalues of K J lie in te interval [1 δ, 1] were δ = C/(C + 2m), m denoting te number of smootings. For standard elliptic operators many smooters can be sown to satisfy te ypoteses, te simplest of wic is te scalar smooter. However, te proof for te scalar smooter and some oters fails in H(div) and te multigrid preconditioner constructed wit tese smooters is not effective. We sall consider an additive Scwarz smooter, defined in te following way. For eac vertex of te mes, consider te patc of elements containing tat vertex. Tese patces form an overlapping covering of Ω and so determine an additive Scwarz operator. We use tis operator as our smooter. Te verification of te first ypotesis is routine. Te standard proof of te second fails, but te difficulty can be surmounted by again using te discrete Helmoltz decomposition in a manner similar to tat used for te proof of domain decomposition. Te complete proof is given in [AFW97]. 1.6 Numerical Results First we made a numerical study of te condition number of J and te effect of preconditioning. In Table 1.1, te level m mes is a uniform triangulation of te unit square into 2 2m 1 triangles and as mes size = 1/2 m 1. Te space V is taken as te Raviart Tomas space of index 0 on tis mes. Te preconditioner K is te V-cycle multigrid preconditioner using one application of te standard additive Scwarz smooter wit te scaling factor taken to be 1/2. Te fift column of te table clearly displays te expected growt of te condition number of J as O( 2 ),
8 8 D. N. ARNOLD, R. S. FALK, AND R. WINTHER Table 1.1 Condition numbers for te operator J and for te preconditioned operator K J, and iterations counts to acieve an error reduction factor of level elements dim V κ(j ) κ(k J ) iterations / / / , / , /32 2,048 3,136 42, /64 8,192 12,416 8 and te sixt column te boundedness of te condition number of te preconditioned operator K J. As a second numerical study, we used te Raviart Tomas mixed metod to solve te factored Poisson equation u = grad p, div u = g in Ω, p = 0 on Ω, again on te unit square using te same sequence of meses as in te first example. We cose g = 2(x 2 +y 2 x y) so tat p = (x 2 x)(y 2 y). Te discrete solution (u, p ) belongs to te space V S, wit V te Raviart Tomas space described above and S te space of piecewise constant functions on te same mes. We solved te discrete equations bot wit a direct solver and by using te minimum residual metod preconditioned wit te block diagonal preconditioner aving as diagonal blocks K and te identity (as discussed previously). Full multigrid was used to initialize te minimum residual algoritm. Tat is, te computed solution at eac level was used as an initial guess at te next finer level, beginning wit te exact solution on te coarsest (two element) mes. In Table 1.2, we sow te condition number of te discrete operator A and of te preconditioned operator B A. Wile te former quantity grows linearly wit 1 (since tis is a first order system), te latter remains small. Table 1.2 Condition numbers for te indefinite operator A corresponding to te mixed system and for te preconditioned operator B A. level dim V dim S κ(a ) κ(b A ) / / / /
9 References [AFW97] Arnold D. N., Falk R. S., and Winter R. (1997) Preconditioning in H(div) and applications. Matematics of Computation to appear. [CLMM94] Cai Z., Lazarov R., Manteuffel T., and McCormick S. (1994) First-order system least squares for second-order partial differential equations: Part I. SIAM J. Numer. Anal. 31: [DW90] Dryja M. and Widlund O. B. (1990) Towards a unified teory of domain decomposition algoritms for elliptic problems. In Can T. F., Glowinski R., Périaux J., and Widlund O. B. (eds) Proc. Tird Int. Conf. on Domain Decomposition Mets., pages SIAM, Piladelpia. [PCL94] Pelivanov A. I., Carey G. F., and Lazarov R. D. (1994) Least-squares mixed finite elements for second-order elliptic problems. SIAM J. Numer. Anal. 31: [Xu92] Xu J. (1992) Iterative metods by space decomposition and subspace correction. SIAM Rev. 34: DD9 Proceedings Editor Petter Bjørstad, Magne Espedal and David Keyes c 1996 Jon Wiley & Sons Ltd.
MULTIGRID PRECONDITIONING IN H(div) ON NON-CONVEX POLYGONS* Dedicated to Professor Jim Douglas, Jr. on the occasion of his seventieth birthday.
MULTIGRID PRECONDITIONING IN H(div) ON NON-CONVEX POLYGONS* DOUGLAS N ARNOLD, RICHARD S FALK, and RAGNAR WINTHER Dedicated to Professor Jim Douglas, Jr on the occasion of his seventieth birthday Abstract
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 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 informationarxiv: v1 [math.na] 20 Jul 2009
STABILITY OF LAGRANGE ELEMENTS FOR THE MIXED LAPLACIAN DOUGLAS N. ARNOLD AND MARIE E. ROGNES arxiv:0907.3438v1 [mat.na] 20 Jul 2009 Abstract. Te stability properties of simple element coices for te mixed
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 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 informationA UNIFORM INF SUP CONDITION WITH APPLICATIONS TO PRECONDITIONING
A UNIFORM INF SUP CONDIION WIH APPLICAIONS O PRECONDIIONING KEN ANDRE MARDAL, JOACHIM SCHÖBERL, AND RAGNAR WINHER Abstract. A uniform inf sup condition related to a parameter dependent Stokes problem is
More informationA SPLITTING LEAST-SQUARES MIXED FINITE ELEMENT METHOD FOR ELLIPTIC OPTIMAL CONTROL PROBLEMS
INTERNATIONAL JOURNAL OF NUMERICAL ANALSIS AND MODELING Volume 3, Number 4, Pages 6 626 c 26 Institute for Scientific Computing and Information A SPLITTING LEAST-SQUARES MIED FINITE ELEMENT METHOD FOR
More informationA SADDLE POINT LEAST SQUARES APPROACH TO MIXED METHODS
A SADDLE POINT LEAST SQUARES APPROACH TO MIXED METHODS CONSTANTIN BACUTA AND KLAJDI QIRKO Abstract. We investigate new PDE discretization approaces for solving variational formulations wit different types
More informationSmoothed projections in finite element exterior calculus
Smooted projections in finite element exterior calculus Ragnar Winter CMA, University of Oslo Norway based on joint work wit: Douglas N. Arnold, Minnesota, Ricard S. Falk, Rutgers, and Snorre H. Cristiansen,
More informationc 2004 Society for Industrial and Applied Mathematics
SIAM J NUMER ANAL Vol 4, No, pp 86 84 c 004 Society for Industrial and Applied Matematics LEAST-SQUARES METHODS FOR LINEAR ELASTICITY ZHIQIANG CAI AND GERHARD STARKE Abstract Tis paper develops least-squares
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 informationFinite Element Methods for Linear Elasticity
Finite Element Metods for Linear Elasticity Ricard S. Falk Department of Matematics - Hill Center Rutgers, Te State University of New Jersey 110 Frelinguysen Rd., Piscataway, NJ 08854-8019 falk@mat.rutgers.edu
More informationFINITE ELEMENT EXTERIOR CALCULUS: FROM HODGE THEORY TO NUMERICAL STABILITY
BULLETIN (New Series) OF THE AMERICAN MATHEMATICAL SOCIETY Volume 47, Number 2, April 2010, Pages 281 354 S 0273-0979(10)01278-4 Article electronically publised on January 25, 2010 FINITE ELEMENT EXTERIOR
More informationNew Streamfunction Approach for Magnetohydrodynamics
New Streamfunction Approac for Magnetoydrodynamics Kab Seo Kang Brooaven National Laboratory, Computational Science Center, Building 63, Room, Upton NY 973, USA. sang@bnl.gov Summary. We apply te finite
More informationAN ANALYSIS OF THE EMBEDDED DISCONTINUOUS GALERKIN METHOD FOR SECOND ORDER ELLIPTIC PROBLEMS
AN ANALYSIS OF THE EMBEDDED DISCONTINUOUS GALERKIN METHOD FOR SECOND ORDER ELLIPTIC PROBLEMS BERNARDO COCKBURN, JOHNNY GUZMÁN, SEE-CHEW SOON, AND HENRYK K. STOLARSKI Abstract. Te embedded discontinuous
More informationAPPROXIMATION BY QUADRILATERAL FINITE ELEMENTS
MATHEMATICS OF COMPUTATION Volume 71, Number 239, Pages 909 922 S 0025-5718(02)01439-4 Article electronically publised on Marc 22, 2002 APPROXIMATION BY QUADRILATERAL FINITE ELEMENTS DOUGLAS N. ARNOLD,
More informationCrouzeix-Velte Decompositions and the Stokes Problem
Crouzeix-Velte Decompositions and te Stokes Problem PD Tesis Strauber Györgyi Eötvös Loránd University of Sciences, Insitute of Matematics, Matematical Doctoral Scool Director of te Doctoral Scool: Dr.
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 informationFigure 1. Figure. Anticipating te use of two-level domain decomposition preconditioners, we construct a triangulation of in te following way. Let be d
Lower Bounds for Two-Level Additive Scwarz Preconditioners wit Small Overlap * Susanne C. Brenner Department of Matematics University of Sout Carolina Columbia, SC 908 Summary. Lower bounds for te condition
More informationMixed Finite Element Methods for Incompressible Flow: Stationary Stokes Equations
Mixed Finite Element Metods for Incompressible Flow: Stationary Stoes Equations Ziqiang Cai, Carles Tong, 2 Panayot S. Vassilevsi, 2 Cunbo Wang Department of Matematics, Purdue University, West Lafayette,
More informationOverlapping domain decomposition methods for elliptic quasi-variational inequalities related to impulse control problem with mixed boundary conditions
Proc. Indian Acad. Sci. (Mat. Sci.) Vol. 121, No. 4, November 2011, pp. 481 493. c Indian Academy of Sciences Overlapping domain decomposition metods for elliptic quasi-variational inequalities related
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 informationDownloaded 11/15/17 to Redistribution subject to SIAM license or copyright; see
SIAM J. NUMER. ANAL. Vol. 55, No. 6, pp. 2787 2810 c 2017 Society for Industrial and Applied Matematics EDGE ELEMENT METHOD FOR OPTIMAL CONTROL OF STATIONARY MAXWELL SYSTEM WITH GAUSS LAW IRWIN YOUSEPT
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 information27. A Non-Overlapping Optimized Schwarz Method which Converges with Arbitrarily Weak Dependence on h
Fourteent International Conference on Domain Decomposition Metods Editors: Ismael Herrera, David E. Keyes, Olof B. Widlund, Robert Yates c 003 DDM.org 7. A Non-Overlapping Optimized Scwarz Metod wic Converges
More informationDOUGLAS N. ARNOLD, RICHARD S. FALK, AND R. WINTHER. Dedicated to Professor Ivo Babuška on the occasion of his seventieth birthday.
MATHEMATICS OF COMPUTATION Volume 66, Number 219, July 1997, Pages 957 984 S 0025-5718(97)00826-0 PRECONDITIONING IN H (div) AND APPLICATIONS DOUGLAS N ARNOLD, RICHARD S FALK, AND R WINTHER Dedicated to
More informationUniform estimate of the constant in the strengthened CBS inequality for anisotropic non-conforming FEM systems
Uniform estimate of te constant in te strengtened CBS inequality for anisotropic non-conforming FEM systems R. Blaeta S. Margenov M. Neytceva Version of November 0, 00 Abstract Preconditioners based on
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 informationNotes on Multigrid Methods
Notes on Multigrid Metods Qingai Zang April, 17 Motivation of multigrids. Te convergence rates of classical iterative metod depend on te grid spacing, or problem size. In contrast, convergence rates of
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 informationAN ANALYSIS OF NEW FINITE ELEMENT SPACES FOR MAXWELL S EQUATIONS
Journal of Matematical Sciences: Advances and Applications Volume 5 8 Pages -9 Available at ttp://scientificadvances.co.in DOI: ttp://d.doi.org/.864/jmsaa_7975 AN ANALYSIS OF NEW FINITE ELEMENT SPACES
More informationarxiv: v1 [math.na] 27 Jan 2014
L 2 -ERROR ESTIMATES FOR FINITE ELEMENT APPROXIMATIONS OF BOUNDARY FLUXES MATS G. LARSON AND ANDRÉ MASSING arxiv:1401.6994v1 [mat.na] 27 Jan 2014 Abstract. We prove quasi-optimal a priori error estimates
More informationBALANCING DOMAIN DECOMPOSITION FOR PROBLEMS WITH LARGE JUMPS IN COEFFICIENTS
MATHEMATICS OF COMPUTATION Volume 65, Number 216 October 1996, Pages 1387 1401 BALANCING DOMAIN DECOMPOSITION FOR PROBLEMS WITH LARGE JUMPS IN COEFFICIENTS JAN MANDEL AND MARIAN BREZINA Abstract. Te Balancing
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] 6 Dec 2010
MULTILEVEL PRECONDITIONERS FOR DISCONTINUOUS GALERKIN APPROXIMATIONS OF ELLIPTIC PROBLEMS WITH JUMP COEFFICIENTS BLANCA AYUSO DE DIOS, MICHAEL HOLST, YUNRONG ZHU, AND LUDMIL ZIKATANOV arxiv:1012.1287v1
More informationA LOCAL FOURIER ANALYSIS FRAMEWORK FOR FINITE-ELEMENT DISCRETIZATIONS OF SYSTEMS OF PDES
A LOCAL FOURIER ANALYSIS FRAMEWORK FOR FINITE-ELEMENT DISCRETIZATIONS OF SYSTEMS OF PDES SCOTT P. MACLACHLAN AND CORNELIS W. OOSTERLEE Abstract. Since teir popularization in te late 1970s and early 1980s,
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 informationDifferent Approaches to a Posteriori Error Analysis of the Discontinuous Galerkin Method
WDS'10 Proceedings of Contributed Papers, Part I, 151 156, 2010. ISBN 978-80-7378-139-2 MATFYZPRESS Different Approaces to a Posteriori Error Analysis of te Discontinuous Galerkin Metod I. Šebestová Carles
More informationON THE CONVERGENCE OF A DUAL-PRIMAL SUBSTRUCTURING METHOD. January 2000 Revised April 2000
ON THE CONVERGENCE OF A DUAL-PRIMAL SUBSTRUCTURING METHOD JAN MANDEL AND RADEK TEZAUR January 2000 Revised April 2000 Abstract In te Dual-Primal FETI metod, introduced by Farat et al [5], te domain is
More informationStokes Equation. Chapter Mathematical Formulations
Capter 5 Stokes Equation 5.1 Matematical Formulations We introduce some notations: Let v = (v 1,v 2 ) T or v = (v 1,v 2,v 3 ) T. Wen n = 2 we define curlv = v 2 x v 1 y, ) curlη =. ( η y, η x In fact,
More informationarxiv: v1 [math.na] 9 Sep 2015
arxiv:509.02595v [mat.na] 9 Sep 205 An Expandable Local and Parallel Two-Grid Finite Element Sceme Yanren ou, GuangZi Du Abstract An expandable local and parallel two-grid finite element sceme based on
More informationA SYMMETRIC NODAL CONSERVATIVE FINITE ELEMENT METHOD FOR THE DARCY EQUATION
A SYMMETRIC NODAL CONSERVATIVE FINITE ELEMENT METHOD FOR THE DARCY EQUATION GABRIEL R. BARRENECHEA, LEOPOLDO P. FRANCA 1 2, AND FRÉDÉRIC VALENTIN Abstract. Tis work introduces and analyzes novel stable
More informationFEM solution of the ψ-ω equations with explicit viscous diffusion 1
FEM solution of te ψ-ω equations wit explicit viscous diffusion J.-L. Guermond and L. Quartapelle 3 Abstract. Tis paper describes a variational formulation for solving te D time-dependent incompressible
More informationMATH745 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 informationSOME IMPLEMENTATIONS OF PROJECTION METHODS FOR NAVIER STOKES EQUATIONS
SOME IMPLEMENTATIONS OF PROJECTION METHODS FOR NAVIER STOKES EQUATIONS JEAN-LUC GUERMOND Abstract. Tis paper is concerned wit te implementation of spatially discrete versions of Corin Temam s projection
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 informationA Weak Galerkin Method with an Over-Relaxed Stabilization for Low Regularity Elliptic Problems
J Sci Comput (07 7:95 8 DOI 0.007/s095-06-096-4 A Weak Galerkin Metod wit an Over-Relaxed Stabilization for Low Regularity Elliptic Problems Lunji Song, Kaifang Liu San Zao Received: April 06 / Revised:
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 informationSome Applications of Fractional Step Runge-Kutta Methods
Some Applications of Fractional Step Runge-Kutta Metods JORGE, J.C., PORTERO, L. Dpto. de Matematica e Informatica Universidad Publica de Navarra Campus Arrosadia, s/n 3006 Pamplona Navarra SPAIN Abstract:
More informationParametric Spline Method for Solving Bratu s Problem
ISSN 749-3889 print, 749-3897 online International Journal of Nonlinear Science Vol4202 No,pp3-0 Parametric Spline Metod for Solving Bratu s Problem M Zarebnia, Z Sarvari 2,2 Department of Matematics,
More informationON THE CONSISTENCY OF THE COMBINATORIAL CODIFFERENTIAL
TRANSACTIONS OF THE AMERICAN MATHEMATICAL SOCIETY Volume 366, Number 10, October 2014, Pages 5487 5502 S 0002-9947(2014)06134-5 Article electronically publised on February 26, 2014 ON THE CONSISTENCY OF
More informationSMAI-JCM SMAI Journal of Computational Mathematics
SMAI-JCM SMAI Journal of Computational Matematics Compatible Maxwell solvers wit particles II: conforming and non-conforming 2D scemes wit a strong Faraday law Martin Campos Pinto & Eric Sonnendrücker
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 informationA SHORT INTRODUCTION TO BANACH LATTICES AND
CHAPTER A SHORT INTRODUCTION TO BANACH LATTICES AND POSITIVE OPERATORS In tis capter we give a brief introduction to Banac lattices and positive operators. Most results of tis capter can be found, e.g.,
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 informationParallel algorithm for convection diffusion system based on least squares procedure
Zang et al. SpringerPlus (26 5:69 DOI.86/s464-6-3333-8 RESEARCH Parallel algoritm for convection diffusion system based on least squares procedure Open Access Jiansong Zang *, Hui Guo 2, Hongfei Fu 3 and
More informationRobust approximation error estimates and multigrid solvers for isogeometric multi-patch discretizations
www.oeaw.ac.at Robust approximation error estimates and multigrid solvers for isogeometric multi-patc discretizations S. Takacs RICAM-Report 2017-32 www.ricam.oeaw.ac.at Robust approximation error estimates
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 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 informationHYBRIDIZED GLOBALLY DIVERGENCE-FREE LDG METHODS. PART I: THE STOKES PROBLEM
MATHEMATICS OF COMPUTATION Volume 75, Number 254, Pages 533 563 S 0025-5718(05)01804-1 Article electronically publised on December 16, 2005 HYBRIDIZED GLOBALLY DIVERGENCE-FREE LDG METHODS. PART I: THE
More informationJian-Guo Liu 1 and Chi-Wang Shu 2
Journal of Computational Pysics 60, 577 596 (000) doi:0.006/jcp.000.6475, available online at ttp://www.idealibrary.com on Jian-Guo Liu and Ci-Wang Su Institute for Pysical Science and Tecnology and Department
More informationarxiv: v1 [math.na] 19 Mar 2018
A primal discontinuous Galerkin metod wit static condensation on very general meses arxiv:1803.06846v1 [mat.na] 19 Mar 018 Alexei Lozinski Laboratoire de Matématiques de Besançon, CNRS UMR 663, Univ. Bourgogne
More informationCOMPUTATIONAL COMPARISON BETWEEN THE TAYLOR HOOD AND THE CONFORMING CROUZEIX RAVIART ELEMENT
Proceedings of ALGORITMY 5 pp. 369 379 COMPUTATIONAL COMPARISON BETWEEN E TAYLOR HOOD AND E CONFORMING OUZEIX RAVIART ELEMENT ROLF KRAHL AND EBERHARD BÄNSCH Abstract. Tis paper is concerned wit te computational
More informationGrad-div stabilization for the evolutionary Oseen problem with inf-sup stable finite elements
Noname manuscript No. will be inserted by te editor Grad-div stabilization for te evolutionary Oseen problem wit inf-sup stable finite elements Javier de Frutos Bosco García-Arcilla Volker Jon Julia Novo
More informationROBUST MULTISCALE ITERATIVE SOLVERS FOR NONLINEAR FLOWS IN HIGHLY HETEROGENEOUS MEDIA
ROBUST MULTISCALE ITERATIVE SOLVERS FOR NONLINEAR FLOWS IN HIGHLY HETEROGENEOUS MEDIA Y. EFENDIEV, J. GALVIS, S. KI KANG, AND R.D. LAZAROV Abstract. In tis paper, we study robust iterative solvers for
More informationComputing eigenvalues and eigenfunctions of Schrödinger equations using a model reduction approach
Computing eigenvalues and eigenfunctions of Scrödinger equations using a model reduction approac Suangping Li 1, Ziwen Zang 2 1 Program in Applied and Computational Matematics, Princeton University, New
More informationComputers and Mathematics with Applications. A nonlinear weighted least-squares finite element method for Stokes equations
Computers Matematics wit Applications 59 () 5 4 Contents lists available at ScienceDirect Computers Matematics wit Applications journal omepage: www.elsevier.com/locate/camwa A nonlinear weigted least-squares
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 the grad-div stabilization for the time-dependent Navier Stokes equations with inf-sup stable finite elements
arxiv:161.517v3 [mat.na] 2 May 217 Analysis of te grad-div stabilization for te time-dependent Navier Stokes equations wit inf-sup stable finite elements Javier de Frutos Bosco García-Arcilla Volker Jon
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 informationSTABILITY OF DISCRETE STOKES OPERATORS IN FRACTIONAL SOBOLEV SPACES
STABILITY OF DISCRETE STOKES OPERATORS IN FRACTIONAL SOBOLEV SPACES JEAN-LUC GUERMOND,, JOSEPH E PASCIAK Abstract Using a general approximation setting aving te generic properties of finite-elements, we
More informationAN EFFICIENT AND ROBUST METHOD FOR SIMULATING TWO-PHASE GEL DYNAMICS
AN EFFICIENT AND ROBUST METHOD FOR SIMULATING TWO-PHASE GEL DYNAMICS GRADY B. WRIGHT, ROBERT D. GUY, AND AARON L. FOGELSON Abstract. We develop a computational metod for simulating models of gel dynamics
More informationTHE INF-SUP STABILITY OF THE LOWEST ORDER TAYLOR-HOOD PAIR ON ANISOTROPIC MESHES arxiv: v1 [math.na] 21 Oct 2017
THE INF-SUP STABILITY OF THE LOWEST ORDER TAYLOR-HOOD PAIR ON ANISOTROPIC MESHES arxiv:1710.07857v1 [mat.na] 21 Oct 2017 GABRIEL R. BARRENECHEA AND ANDREAS WACHTEL Abstract. Uniform LBB conditions are
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 informationH(div) conforming and DG methods for incompressible Euler s equations
H(div) conforming and DG metods for incompressible Euler s equations Jonny Guzmán Filánder A. Sequeira Ci-Wang Su Abstract H(div) conforming and discontinuous Galerkin (DG) metods are designed for incompressible
More informationStability properties of a family of chock capturing methods for hyperbolic conservation laws
Proceedings of te 3rd IASME/WSEAS Int. Conf. on FLUID DYNAMICS & AERODYNAMICS, Corfu, Greece, August 0-, 005 (pp48-5) Stability properties of a family of cock capturing metods for yperbolic conservation
More informationFINITE ELEMENT DUAL SINGULAR FUNCTION METHODS FOR HELMHOLTZ AND HEAT EQUATIONS
J. KSIAM Vol.22, No.2, 101 113, 2018 ttp://dx.doi.org/10.12941/jksiam.2018.22.101 FINITE ELEMENT DUAL SINGULAR FUNCTION METHODS FOR HELMHOLTZ AND HEAT EQUATIONS DEOK-KYU JANG AND JAE-HONG PYO DEPARTMENT
More informationPARTITION OF UNITY FOR THE STOKES PROBLEM ON NONMATCHING GRIDS
PARTITION OF UNITY FOR THE STOES PROBLEM ON NONMATCHING GRIDS CONSTANTIN BACUTA AND JINCHAO XU Abstract. We consider the Stokes Problem on a plane polygonal domain Ω R 2. We propose a finite element method
More informationKey words. Sixth order problem, higher order partial differential equations, biharmonic problem, mixed finite elements, error estimates.
A MIXED FINITE ELEMENT METHOD FOR A SIXTH ORDER ELLIPTIC PROBLEM JÉRÔME DRONIOU, MUHAMMAD ILYAS, BISHNU P. LAMICHHANE, AND GLEN E. WHEELER Abstract. We consider a saddle point formulation for a sixt order
More informationLinearized Primal-Dual Methods for Linear Inverse Problems with Total Variation Regularization and Finite Element Discretization
Linearized Primal-Dual Metods for Linear Inverse Problems wit Total Variation Regularization and Finite Element Discretization WENYI TIAN XIAOMING YUAN September 2, 26 Abstract. Linear inverse problems
More informationDedicated to the 70th birthday of Professor Lin Qun
Journal of Computational Matematics, Vol.4, No.3, 6, 4 44. ACCELERATION METHODS OF NONLINEAR ITERATION FOR NONLINEAR PARABOLIC EQUATIONS Guang-wei Yuan Xu-deng Hang Laboratory of Computational Pysics,
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 informationA Least-Squares Finite Element Approximation for the Compressible Stokes Equations
A Least-Squares Finite Element Approximation for the Compressible Stokes Equations Zhiqiang Cai, 1 Xiu Ye 1 Department of Mathematics, Purdue University, 1395 Mathematical Science Building, West Lafayette,
More informationNumerical analysis of a free piston problem
MATHEMATICAL COMMUNICATIONS 573 Mat. Commun., Vol. 15, No. 2, pp. 573-585 (2010) Numerical analysis of a free piston problem Boris Mua 1 and Zvonimir Tutek 1, 1 Department of Matematics, University of
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 informationLIMITS AND DERIVATIVES CONDITIONS FOR THE EXISTENCE OF A LIMIT
LIMITS AND DERIVATIVES Te limit of a function is defined as te value of y tat te curve approaces, as x approaces a particular value. Te limit of f (x) as x approaces a is written as f (x) approaces, as
More informationEXPLICIT ERROR BOUNDS IN A CONFORMING FINITE ELEMENT METHOD
MATHEMATICS OF COMPUTATION Volume 68, Number 228, Pages 1379 1396 S 0025-5718(99)01093-5 Article electronically publised on February 24, 1999 EXPLICIT ERROR BOUNDS IN A CONFORMING FINITE ELEMENT METHOD
More informationMass Lumping for Constant Density Acoustics
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
More informationCLEMSON U N I V E R S I T Y
A Fractional Step θ-metod for Convection-Diffusion Equations Jon Crispell December, 006 Advisors: Dr. Lea Jenkins and Dr. Vincent Ervin Fractional Step θ-metod Outline Crispell,Ervin,Jenkins Motivation
More informationConvergence and Descent Properties for a Class of Multilevel Optimization Algorithms
Convergence and Descent Properties for a Class of Multilevel Optimization Algoritms Stepen G. Nas April 28, 2010 Abstract I present a multilevel optimization approac (termed MG/Opt) for te solution of
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 informationOptimization of stress modes by energy compatibility for 4-node hybrid quadrilaterals
INTERNATIONAL JOURNAL FOR NUMERIAL METHODS IN ENGINEERING Int. J. Numer. Met. Engng 2004; 59:293 33 (DOI: 0.002/nme.877) Optimization of stress modes by energy compatibility for 4-node ybrid quadrilaterals
More informationA domain decomposition convergence for elasticity equations
A domain decomposition convergence for elasticity equations A. Cakib a, A. Ellabib b, A. Nacaoui c, a Département de Matématiques Appliquées et Informatique FST de Beni-Mellal, Université Sultan My Slimane
More information4. The slope of the line 2x 7y = 8 is (a) 2/7 (b) 7/2 (c) 2 (d) 2/7 (e) None of these.
Mat 11. Test Form N Fall 016 Name. Instructions. Te first eleven problems are wort points eac. Te last six problems are wort 5 points eac. For te last six problems, you must use relevant metods of algebra
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 informationERROR CONTROL AND ADAPTIVITY FOR A VARIATIONAL MODEL PROBLEM DEFINED ON FUNCTIONS OF BOUNDED VARIATION
ERROR CONTROL AND ADAPTIVITY FOR A VARIATIONAL MODEL PROBLEM DEFINED ON FUNCTIONS OF BOUNDED VARIATION SÖREN BARTELS Abstract. We derive a fully computable, optimal a posteriori error estimate for te finite
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 informationLecture Note III: Least-Squares Method
Lecture Note III: Least-Squares Method Zhiqiang Cai October 4, 004 In this chapter, we shall present least-squares methods for second-order scalar partial differential equations, elastic equations of solids,
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 information