A CLASS OF DISCONTINUOUS PETROV-GALERKIN METHODS. PART I: THE TRANSPORT EQUATION
|
|
- Gabriel Perry
- 6 years ago
- Views:
Transcription
1 A CLASS OF DISCONTINUOUS PETROV-GALERKIN METHODS. PART I: THE TRANSPORT EQUATION L. DEMKOWICZ AND J. GOPALAKRISHNAN Abstract. Considering a simple model transport problem, we present a new finite element metod. Wile te new metod fits in te class of discontinuous Galerkin (DG) metods, it differs from standard DG and streamline diffusion metods, in tat it uses a space of discontinuous trial functions tailored for stability. Te new metod, unlike te older approaces, yields optimal estimates for te primal variable in bot te element size and polynomial degree p, and outperforms te standard upwind DG metod.. Introduction We introduce a new Petrov-Galerkin metod for advective problems. Wile it belongs in te class of discontinuous Galerkin (DG) metods, unlike te standard upwind DG metod, we are able to prove optimal and p error estimates in te L 2 -norm for our discrete solution on general meses (were, as usual, is te mes size and p is te polynomial degree). Te metod includes a separate outflux approximation on element interfaces and a space of non-standard test functions designed for stability. Te boundary value problem tat is te subject of tis paper is posed on a polyedral domain Ω. Given f and g, we need to find a finite element approximation to te solution u of β u = f on Ω, (a) u = g on in Ω. (b) We only consider te case of constant β in tis paper (but extensions are possible, as mentioned in Section 5). Te inflow boundary in Ω appearing in (), is defined, letting n denote te unit outward normal, by in Ω = { x Ω : β n( x) < 0}, (2) i.e., in Ω denotes te global inflow boundary. Wile non-finite-element numerical tecniques can be designed for tis problem (e.g. te metod of caracteristics), we aim for finite elements because of its versatility in andling complicated domains as well as certain regular and singular perturbations of te above problem. A regular perturbation of () is β u + α( x)u = f. (3) A singular perturbation of () is obtained by te addition of a small viscosity term wit second derivatives. Tis is arder to analyze. Witin te domain of finite element metods for (), tere are two broad categories (see [] for a review). One is te very popular Te work of te first autor was supported by DOE troug Predictive Engineering Science (PECOS) Center at ICES (PI: Bob Moser), and by a researc contract wit Boeing. Te second autor was supported in part by te National Science Foundation under grants , , and an Oden fellowsip at ICES.
2 2 L. DEMKOWICZ AND J. GOPALAKRISHNAN streamline diffusion metod [3] and its descendants. Te oter category is composed of DG metods. Since our contribution fits in te latter, we sall now review previous works in tis category in detail. Te well known first papers proposing and analyzing te original DG metod for () are [4, 5, 8]. To distinguis tis metod from our DG metod, we will call te original DG metod te upwind DG metod and denote it by UDG, wile we call ours te discontinuous Petrov-Galerkin metod and denote it by DPG. It is proved in [5] tat if u is te UDG approximation, ten (for a fixed p) it satisfies u u L 2 (Ω) O( s ) for some s p + dictated by te regularity of te exact solution. Tis result was improved by [4] werein it was sown tat te rate of convergence is in fact O( s /2 ). In bot cases convergence wit respect to p was not studied. Even if we set aside te p-convergence issue, notice tat bot te results are suboptimal in, as te best approximation error of te finite element space is O( s ). For some special classes of meses owever, many autors ave observed (and proved) te optimal rate of convergence of te UDG metod [8, 9] wit respect to. Noneteless, on general meses, te suboptimal rate of convergence cannot be improved, as sown by a numerical example in [7] using a particular quasiuniform mes and a smoot exact solution. To express te sentiment of many, we quote from [8] tat te mecanisms tat induce te loss of /2 in te order of convergence of te L 2 -norm of te error are not very well known yet. An p analysis of te UDG sceme was first provided in [3]. Tey considered te regular perturbation (3) under te assumption tat 0 < c 0 α( x) x Ω (4) Because of tis assumption, tey are able to control te L 2 (Ω)-norm of te solution. Tey also introduced a stabilization parameter into te original upwind DG metod. A few years later, te paper [2] extended te results of [3] in several directions, providing a unified teory for an p version of te streamline diffusion metod, as well as te upwind DG metod. Teir analysis did not assume (4), rater tey let te advection vector β depend on x and assumed 0 < c 0 2 β( x) + α( x) x Ω, (5) as a consequence of wic tey ave stability in c 0 u L 2 (Ω). (Note tat for te case we intend to study (), te rigt and side of (5) evaluates to zero, so (5) does not old.) Bot papers analyzed te stabilized version of te upwind DG metod and bot relied on te proper coice of te stabilization parameter. Wile te results of [2] are optimal in p, tose of [3] are suboptimal in p. In all tese works, te L 2 (Ω)-rate of convergence of te error wit respect to remained suboptimal by /2. In contrast, our results do not exibit tis suboptimality, nor do we add any stabilization parameter. We believe our metod is te first in te finite element (not just DG) family of metods for te transport equation wic as provably optimal convergence rates on very general meses. Te design of our metod is guided by a generalization of Céa lemma due to Babuška [, 5]. We only need a simple version of te result, wic we now describe using te following notations (all our spaces are over R): Let X, X X, and V be Banac spaces and let a (, ) be a bilinear form on X V. Suppose te exact solution U X satisfies and te
3 DPG METHOD 3 discrete solution U X satisfies a (U U, v ) = 0 for all v V. (6) If te bilinear form is continuous in te sense tat tere is a C > 0 suc tat a (w, v ) C w X v V for all w X, v V (7) and also te inf-sup condition, i.e., tere is a C 2 > 0 suc tat C 2 w X sup v V a (w, v ) v V for all w X, (8) ten, as is well known, te following teorem can be formulated: Teorem.. Under te above setting, we ave te following error estimate: U U X ( + C ) inf U w X. C 2 w X Proof. Te argument is simple and standard: U U X = U w X + w U X U w X + C 2 U w X + C 2 a (w U, v ) sup by (8), v V v V a (w U, v ) sup by (6), v V v V U w X + C C 2 w U X by (7). Tis finises te proof. Many refinements and improvements of suc a teorem are known. But our purpose in going troug te above simple argument is to clearly sow tat te test space need not ave approximation properties. Hence, in designing Petrov-Galerkin metods, wile we must coose trial spaces wit good approximation properties, we may design test spaces solely to obtain good stability properties. Tis will be our guiding principle in designing our metod. In fact, te test spaces we propose sortly can ave discontinuities inside te mes elements. Many researcers ave put te above principle to good use. In fact, even te abbreviation we use for our new metod DPG metod, as been previously used [4, 6] for oter metods. Te teme in tese works is te searc for stable test spaces using bubbles or oter polynomials. Our test space functions, in contrast, need not be polynomial on an element, and indeed, need not even be continuous. We are also not te first to consider suc functions wit discontinuities witin a finite element. Suc elements are routinely used in X-FEM and similar metods [2] for difficult simulations like crack propagation. However, we use discontinuities solely for stability purposes, and solely in test spaces. Our trial spaces, being standard polynomial spaces, possess provably good approximation properties. Our metod also introduces a new flux unknown on te element interfaces. Tis is in line wit te recent developments on ybridized DG (HDG) metods [9]. HDG metods tat extend te ideas in [9] to te case of convection can be found in recent works [0, 6]. Tese metods are constructed by defining a independent flux variable on te element
4 4 L. DEMKOWICZ AND J. GOPALAKRISHNAN interfaces wic can solved for first, after wic te internal variables can be locally solved for. Wile tis can be tougt of as akin to static condensation, additional advantages can be exploited, suc as easy stabilization [6] using a penalty parameter. However, p- independent stability for suc metods as not been proved yet. Wile we borrow te idea of letting te fluxes be independent variables in te design of our metod, our metod does not ave stabilization parameters, and as p-independent stability. We organize our presentation suc tat a spectral version of te metod is first exibited (in Section 2). Details regarding te new space of test functions and te stability estimates for te metod on a single element are presented in tat section. Section 3 ten presents te composite metod on a triangular mes. Optimal L 2 error estimates are proved in Teorem 3.2 tere. We conclude in Section 5 opining on important future directions. Proofs of a few tecnical estimates are gatered in Appendix A. 2. Te spectral metod on one element We start by considering te one-element case to fix te ideas and study te element spaces. In oter words, we let Ω be an N-simplex (N = 2, 3 are te cases of most practical importance) wic we denote by K trougout tis section. We now describe and analyze te DPG metod on te single element K. Define te outflow and inflow boundaries of K by out K = { x K : β n > 0}, in K = { x K : β n < 0}. Here, and trougout te paper, we use n to generically denote te unit outward normal on te boundary of a polygonal domain (te domain will vary at different occurrences of n, but will always be clear from te context, e.g., above it is K, wile in (2) it was Ω). Note tat we omitted te β n = 0 case in (9). Tis, and any concerns regarding computational inaccuracies in determination of strict inequalities of (9), will be allayed in Petrov-Galerkin metod using a new test space. Before we describe te metod, we need to introduce a new test space. For tis, we need a coordinate system aligned wit te flow see Figure. Let e (),..., e (N), denote te standard unit vectors and let e β denote te unit vector in te e β -direction. Completing { e β } to form an ortonormal basis { e β, e (2) (N) β,... e β } for R N, we write te new coordinates as η i, i.e., x = x e () + + x N e (N) = η e β + η 2 e (2) β + + η N e (N) β. Any function defined on out K can be extended into K in suc a way tat te extension is constant along te β-direction, i.e., te extension is independent of η. We call tis te extension from outflow boundary and denote it by E out. We can tink of te outflow surface as te grap of a function η = F out (η 2,..., η N ). Ten, a function on te outflow surface can be expressed (after eliminating te η variable) as φ(η 2,..., η N ). Tese coordinates make te expression for E out φ trivial, namely E out φ(η, η 2,..., η N ) = φ(η 2,..., η N ). Let us consider te question of computing a polynomial approximation of te transport solution u satisfying () wit Ω = K. Te starting point in deriving te DPG metod is (9a) (9b) (0)
5 DPG METHOD 5 η = (u, β v) K + β n u, v K 2.. Petrov-Galerkin metod using a new test space Note tat if out = (u, β v) K + φ, v outk + β n g, v in K () V (K) consists of wit φ = β n u on Before out K, we tedescribe outflux variable. te metod, we need to introduce Te metod finds an approximation to te pair (u, φ) using te spaces non-polynomial sp a new test space. For tis, we need a coordinate system M l ( out K) = {µ : µ F P l (F ), for all faces F witin K. aligned wit te flow. oflet K contained e (),...,e in (N) out, K}, denote te standard unit vectors V l (K) = and η P l let (K) e + β (2a) Te spectral DP Edenote out (M l+ te ( out unit K)). vector in te (2b) approximations (u Here, for any domain e β -direction. D, we use Completing P l (D) to denote {e β } to teform spacean of polynomials ortonormal of degree at most l, and fying basis {e β,e (2) β η,...e(n) β } for R N P l (K) = {η p l : p l, we P l write (K)}. te new coordinates as η i,i.e., (u p, β v Note tat if out K consists of a single face of K, ten V l (K) consists of polynomials, but in general, V l (K) is a non-polynomial space, as it can ave lines of discontinuity witin K. Te spectral DPG metod, motivated x = x e () by (), defines + + x N e (N) approximations (u p, φ p+ ) P p (K) M p+ ( out K), satisfying (0) for all v V = η e β η 2 e (2) β + + η N e (N) p (K). (u p, β v) K + φ p+, v outk β. = (f, v) K (3) β n g, v in K, a((u p,φ p+ ),v) Any function defined on out can be extended into K for all v V p (K). For compact notation, we write in suc a way tat te extension is constant along te β- b(v) a( (u p, φ p+ ), v) def = (u p, β direction, i.e., te extension is v) K + φ p+, v outk independent of η. We call X p (K) b(v) def = (f, v) K β n g, v in K, tis te extension from outflow boundary and denote it by X p (K) def = P p (K) M p+ ( out K), so tat (3) can E out. so tat (3) can be written as te problem of finding (u p, φ p+ ) X p (K) suc tat (u p,φ p+ ) X p (K Let us consider te question of computing a polynomial a( (u p, φ p+ ), v) = b(v) v V p (K). approximation of te transport solution u satisfying () a((u p,φ p wit Ω = K. Te starting point in deriving te DPG Tis is a Petrovtrial spaces are metod is to multiply (a) by a test function v and inte- di using te spaces η 2 β K M ( out K)= V (K) = Figure. Streamline coordinates regarding computational inaccuracies in determination of to multiply (a) strict by a inequalities test functionof v and (9), integrate will be allayed by parts: in 2.3. (f, v) = ( β u, v) K Here, for any dom of polynomials of η P (
6 6 L. DEMKOWICZ AND J. GOPALAKRISHNAN Tis is a Petrov-Galerkin metod because te test and trial spaces are different. Te solvability of tis metod will be clear sortly, as soon as we study te new finite element space. Proposition 2.. Te component spaces η P p (K) and E out (M p+ ( out K)) are linearly independent, so tat V p (K) can be written as te direct sum V p (K) = η P p (K) E out (M p+ ( out K)). (4) Proof. Suppose tere is an r P p (K) and µ M p+ ( out K) suc tat η r = E out µ, ten, letting η denote te distributional derivative wit respect to te variable η, we ave η (η r) = η (E out µ) = 0. Integrating from te η = 0 plane (extending te polynomial r to all R 3 if necessary) we find tat t r(t, η 2,..., η N ) = t 0 η (η r(η,..., η N )) ds = 0, for any t, so r 0, wic in turn implies tat µ = 0 as well. Tus te sum in (4) is indeed a direct sum. Te degrees of freedom of V p (K) can be inferred from te next lemma. Lemma 2.. Given any µ in M p+ ( out K) and any w in P p (K), tere is a unique v in V p (K) satisfying β v = w, (5a) v outk = µ. (5b) Proof. We construct te required v as a sum of two functions v and v 2, defined below. First set η v (η,..., η N ) = 0 β w(s, η 2,..., η N ) ds, were β denotes te lengt of te vector β. Note tat v 2, wen restricted to eac face of out K is a polynomial of degree at most p +. Now let v 2 = E out (µ v outk). It is easy to verify tat v = v + v 2 satisfies bot te equalities of (5). Te uniqueness of v is also easy to see: If v satisfies (5) wit w = 0 and µ = 0, ten integrating into K from te outflow boundary along te β-direction, we find tat v 0. One can now define a new finite element formally in te sense of Ciarlet [7] as follows. Let Σ denote te set of linear functionals l q, l η defined by l q (v) = ( β v, q) K, for all q P p (K), l η (v) = v, η F, for all η P p+ (F ), for all faces F out K. Ten, wit Σ as te set of degrees of freedom (d.o.f), te geometry-space-d.o.f triple (K, V p (K), Σ) defines a unisolvent finite element because of Lemma 2.. Altoug V p (K)
7 DPG METHOD 7 may contain discontinuous functions, implementation of tis finite element can proceed using standard finite element tecnology using te above d.o.fs. Proposition 2.2. Tere is a unique solution to te spectral DPG metod (3). Proof. Te system (3) is square. Tis follows because te space of trial functions as dimension dim(x p (K)) = dim(p p (K)) + dim(m p+ ( out K)) = dim(η P p (K)) + dim(e out (M p+ ( out K))) wic equals te dimension of te test space V p (K), due to Proposition 2.. Hence it suffices to sow tat if te rigt and side of (3) vanises, te only possible solution is trivial. But tis is immediate from Lemma 2., by wic we may coose v suc tat β v = u p and v outk = φ p+. Let us now outline an error analysis. We want to apply Teorem.. Te continuity of te form a(, ) is evident from a( (u p, φ p+ ), v) = (u p, β v) K + φ p+, v Ko u p K β v K + /2 K φ p+ outk /2 K v outk (6) ( u p 2 K + K φ p+ 2 outk) /2 v K were we ave cosen te test space norm ( v K = β /2 v 2 K + K v 2 outk). Te inf-sup condition also olds, because applying Lemma 2., we can find a v suc tat β v = u p and v outk = K φ p+, so a( (u p, φ p+ ), v) sup v V p(k) v K = u p 2 K + K φ p+ 2 outk v K ( up 2 K + K φ p+ 2 outk) /2 v K v K. Tus applying Teorem., we obtain te following error estimate for te spectral metod: u u p 2 K + K φ φ p+ 2 K ( 2 inf u wp 2 K + K φ µ p+ K) 2. (w p,µ p+ ) X p(k) Here K = diam(k). Reviewing te above argument, te introduction of factors K and K in (6) may seem arbitrary, and indeed it is. Wile we cose tese factors above for simplicity, later we will need to introduce different mes dependent factors to obtain stability in more appropriate norms. Altoug we used te inf-sup condition to motivate te design of our metod, te way we used it above does not provide us wit te best error estimate possible. It is possible to obtain better error estimates, and indeed to sow tat te error for eac variable is decoupled. In fact, bot variables possess superconvergence in te spectral case, as we see next.
8 8 L. DEMKOWICZ AND J. GOPALAKRISHNAN Teorem 2.. Let Π W and Π M denote te L 2 -ortogonal projections into P p (K) and M p+ ( out K), respectively. Ten φ p+ Π M φ = 0, u p Π W u = 0. Proof. Subtracting () from (3), we find tat (u p Π M u, β v) K + φ p+ Π M φ, v outk = 0 for all v in V p (K). By Lemma 2., we can coose v in V p (K) so tat v outk = E out (φ p+ Π M φ) and β v = Π M u u p. Hence te identities of te teorem follow Local stability estimates. Altoug our error analysis of te spectral metod is already completed (by Teorem 2.), we want to obtain a few refined stability estimates for te solutions. Tis will be needed later, wen we analyze te multi-element version of te metod. To tis end, we introduce te following local solution operators: Q p : L 2 ( in K) M p+ ( out K), U p : L 2 ( in K) P p (K). Given φ in L 2 ( in K), te local solutions Q p φ and U p φ are defined by te spectral metod (3) wit f = 0, element by element, as follows: (U p φ, β v) K + Q p φ, v outk = φ, v in K (7) for all v V p (K). Also define solution operators for element sources by Q f p : L 2 (K) M p+ ( out K), U f p : L 2 (K) P p (K) (U f pz, β v) K + Q f pz, v outk = (z, v) K, (8) for all v V p (K) and any z in L 2 (K). Te next teorem is te main result of tis subsection. Let β,k denote te lengt of te longest line segment, in te β-direction, contained in K, wit one endpoint in in K and te oter in out K (as marked in Fig. 0(b)). We denote by n F a unit normal to a mes face F. For any subcollection F of te set of faces of K were β n 0, we define te following norms: Here te sums run over all mes faces F contained in F. ψ 2 β,f = F F β n F ψ 2 F, (9) v 2 β,f = F F β n F v 2 F. (20) Teorem 2.2. Te local solution operators ave te following stability estimates: Q p φ 2 β, outk φ 2 β, ink, U p φ 2 K β,k β φ 2 β, ink, (2a) (2b)
9 DPG METHOD 9 8 L. DEMKOWICZ AND J. GOPALAKRISHNAN outk ( β n 0ere) K outk β 0 β Figure 2.. AAoutflowedgedegeneratingas degenerating as β β 00 and 2.3. Remarks on te robustness of te outflux variable. Notice tat on faces were β n =0,weavenotyetdefinedanyoutfluxapproximation. Q f Sinceteexactoutflux pz 2 (defined in (26)) vanises on suc faces, we β,k β, outk set β te z 2 K, (22a) discrete outflux φ to zero on tose faces. Computationally, te identification Uof f pz inflow 2 K 2 β,k and outflow faces can only be as accurate as te round-off errors permit, i.e., te strict inequalities β 2 z 2 K, (22b) (9)canonlybeimplementedup for to round-off. all φ L 2 ( Tis begs te question: On faces F were β in K) and z in L 2 (K). n F is close to zero, is φ p too sensitive to weter we set F to be an inflow or outflow face? AWe proof nowof argue te teorem tat it is is not. givenfor in Appendix clarity, consider A. te situation of Figure (wic is entirely typical). Te solution φ p of (3) can be written as φ p = Q p g + Q f pf. Ten,byte 2.3. estimates Remarks of Teorem on te2.2, robustness we ave of te outflux variable. Notice tat on faces were β n = 0, we ave not yet defined any outflux approximation. Since te exact outflux (defined in (26)) vanises on suc faces, we set te discrete outflux φ φ p β, outk g β, ink + /2 to zero on tose β,k faces. β f K. Computationally, te identification of inflow and outflow /2 faces can only be as accurate as te round-off errors permit, Tis continues to old as β i.e., te strict inequalities (9) can only be implemented up to round-off. Tis begs te question: n approaces On zero faceson F an were outflow β n face F is(as close in Fig. to zero, ). Note is φ p tat too sensitive te rigtto and weter side remains we set Fbounded to during tis limiting process. But since te left and side norm contains te factor β be an inflow or outflow face? We now argue tat it is not. For n /2 clarity, on te consider outflowteface situation were β of Figure n approaces 2 (wic zero, is entirely te solution typical). φ p must Teapproac solution φzero p of tere. (3) can be written as φ p = Q p g + Q f pf. Ten, by te estimates It is tis of robustness Teorem 2.2, in we φ p ave tat led us to introduce it as te new outflux variable to be solved for. Te coice of suc a variable is not traditional. In ybridized metods for elliptic problems [5] one usually sets te trace of te primal variable u as a new unknown. φ In tis spirit, one could contemplate p β, outk a metod g β, ink were + /2 β,k β one f finds K. /2 u p and λ p satisfying Tis continues to (uold p, β as β v) n K approaces + β n λ p,v zero outk on=(f,v) an outflow K face β n (as g, v in in Fig. K, 2). Note tat (23) te rigt and side remains bounded during tis limiting process. But since te left and side instead normof contains (3). Suc te metods factor ave β n appeared /2 on te in te outflow literature face were testandard β n approaces UDGmetod zero, te cansolution be recast φ p into must ybrid approac formzero wittere. upwind traces set as a new unknown λ,seee.g., te Itmetod is tis robustness in [6] for te in purely φ p tat convective led us tocase introduce and teitrecent as teworkof[]).however,an new outflux variable to be equation solved suc for. as Te (23) coice can lead of suc to matrices a variable close is not to singularity traditional. inin situations ybridized likemetods in Figurefor. elliptic Indeed, problems te λ in [9] (23) onecannot usuallybe sets uniquely te trace defined of te were primal β n variable =0. u as a new unknown.
10 0 L. DEMKOWICZ AND J. GOPALAKRISHNAN In tis spirit, one could contemplate a metod were one finds u p and λ p satisfying (u p, β v) K + β n λ p, v outk = (f, v) K β n g, v in K, (23) instead of (3). Suc metods ave appeared in te literature te standard UDG metod can be recast into ybrid form wit te upwind traces set as a new unknown λ, see e.g., te metod in [0] for te purely convective case). However, an equation suc as (23) can lead to matrices close to singularity in situations like in Figure 2. Indeed, te λ in (23) cannot be uniquely defined were β n = Approximations using a triangular mes In tis section, we move to te multi-element case, restricting ourselves to te case of two space dimensions. We consider te composite DPG metod on a mes of te domain Ω, and a general polynomial degree p. Since te degree is uniformly set to p on all elements tis is te p-version of te metod. We assume tat te domain Ω is mesed by a geometrically conforming mes of triangles, denoted by T, and set = max{ K : K T }. Define te (discontinuous) finite element spaces W = {w : w K P p (K), for all mes elements K T }, M = { µ : µ F P p+ (F ), for all mes edges F not on in Ω}, V = { v : v K V p (K), for all mes elements K T }. Here V p (K) is te new finite element we introduced in Section 2. (24a) (24b) (24c) 3.. Definition of te composite metod. We put togeter te spectral metod on eac element to get te following composite metod on te mes T : Find (u, φ ) def = W M satisfying X for all v V, were a ( (u, φ ), v ) = (f, v ) Ω β n g, v in Ω, (25) a ( (u, φ ), v ) = ] [ φ, v outk φ, v ink\ inω (u, β v ) K. K T We call φ te outflux approximation, because on comparison wit (), we expect it to approximate φ def = βu n + (26) on every mes edge F, were n + is te unit normal on F wit its direction cosen so tat β n + > 0. On edges were β n = 0, we set φ = 0. Equation () also sows tat tis φ togeter wit te exact solution u, satisfies a ( (u, φ), v) = (f, v) Ω β n g, v in Ω, (27)
11 DPG METHOD for all v in V = {ν L 2 (Ω) : ν K V (K) for all mes elements K}, were V (K) is as in (55). In oter words, te DPG metod is consistent. In order to simplify notation and write φ, v in K for φ, v in K\ in Ω, we identify te space M wit its extension by zero to in Ω, i.e., we identify M and {µ : µ F P p+ (F ) for all mes edge F and µ in Ω = 0}, to be te same (cf. te definition in (24b)) Mes dependent norms. We collect te definitions of various mes dependent norms ere for ready reference. First, we extend te definition of te norms ψ 2 and v β,f β,f to more general collections F of edges from te entire mes (now not necessarily edges of one triangle) simply by letting te sums in (9) and (20) runs over all suc edges in F. As before, only edges wit β n 0 are picked. Using tese extended definitions, we now define more norms on a subcollection of mes elements R, as follows: v 2,β,R = β,k v 2 β,, (28) outk K R ψ 2 =, β,r β,k ψ 2 β K R ψ 2 = ψ 2, β,r β,k β K R, outk, (29), outk, (30) ψ 2, β,r = ψ Q pψ 2, β,r. (3) Tese norm notations ave mnemonic subscripts of, /, β, or /β etc., wic indicates a factor of β,k, / β,k, β n, or / β n (resp.) in te sums over mes objects defining te norms. Wit tese notations, te norm on te trial space is defined by (u, φ) 2 = K T ( u U p φ 2 K + ψ Q p ψ 2, β,k Te norm on te test space is defined by v 2 = K T β v 2 K + v 2,β,K. It is obvious from te construction of V p (K) tat te above is a norm on V. To sow tat (u, φ) is a norm on te test space, we first need to recall te following. Lemma 3.. In any triangulation T of Ω, tere is at least one triangle K T suc tat in K in Ω. Proof. See [5]. Proposition 3.. Te functional (u, φ) is a norm on te space X = L 2 (Ω) M, were M = {µ : µ F L 2 (F ) for all mes edges F and µ in Ω = 0}. Proof. Suppose (u, φ) = 0. Ten, in particular φ Q p φ β, outk = 0 (32) ).
12 2 L. DEMKOWICZ AND J. GOPALAKRISHNAN for all elements K. By Lemma 3., tere is an element K wit in K in Ω. On tis element, φ in K = 0, so Q p φ outk = 0. Hence (32) implies tat φ outk = 0. Proceeding inductively over all elements we find tat φ = 0. Tis in turn implies tat u U p φ = u = Preliminary error analysis. It is possible to get stability and error estimates for te DPG metod in te mes dependent norm immediately. Solvability of te metod follows as a particular consequence. Teorem 3.. Te bilinear form a (, ) satisfies te inf-sup and continuity conditions: olds for all u W and φ M, and a ( (u, φ ), v ) sup (u, φ ), v V v a ( (u, φ), v ) (u, φ) v for all u L 2 (Ω), φ M, and v V. If (u, φ) and (u, φ ) are te exact and discrete solutions, ten (u u, φ φ ) 2 2 inf (u w, φ µ ) 2. (w,µ ) W M Proof. We first rewrite (25) using te definition of te local solution operators in (7) as follows: a ( (u, φ ), v ) = K T [ (u, β v ) K + φ, v outk] φ, v in K = K T (u U p φ, β v ) K + φ Q p φ, v outk. (33) Let us prove te inf-sup condition. By Lemma 2., we can coose v element by element, suc tat β v K = u U p (φ in K ), Ten ( φ v outk = Q p (φ in )) K outk β,k β. n a ( (u, φ ), v ) = K T u U p φ 2 K + φ Q p φ 2, β,k = (u, φ ) v. Tis proves te inf-sup condition. Te proof of te continuity inequality also follows from (33) after applying te Caucy- Scwarz inequality. Te error estimate ten immediately follows from Teorem..
13 DPG METHOD 3 2 L. DEMKOWICZ AND J. GOPALAKRISHNAN S 3 S 2 S S 0 S 2 S S 0 S S 0 S 0 Figure Figure 2. Amesbeingsplitintolayers(flowisalongtepositivey-direction). 3. A split into layers is along te positive y- direction). By Lemma 3., S 0 is not empty. Next, setting T (0) = T,werecursivelydefine,forall integers l>0, te sets Te above analysis is in te most natural norms suggested by te metod itself. In fact, we ave sown above T (l) tat = T (l ) \ S l S l = {K a T (l) : in K in Ω (l) }, ( (u, φ ), v ) sup = (u, φ ) were Ω (l) is te domain formed by te union of all triangles of T (l), v.repeatedapplication V v of Lemma 3. sows tat S l is nonempty for every l, unlesst (l) is empty. So tis process exausts for all allu elements W, of φ te Mmes,. Noneteless, at some finite convergence value ofrates l, say inl standard = n. Te norms domain (like L 2 ) def are not obvious Λ l = Ω\ Ω (l) from te teorem. Te remainder of tis section is devoted to obtaining is mesed by S 0 S S l.itsoutflowboundaryisdenotedby convergence rates in L 2 -like norms. Γ l = out Λ l. All mes 3.4. Aedges, Poincaré excluding inequality. tose on in As Ω, aand first excluding step to tose obtainwere stability β n =0,arecontained and error estimates in in Γ l Lfor 2 -norm, some lwe prove n. Wittesenotationsandobservations,weimposeamildassumption a Poincaré inequality. Notice tat te norm in (3) is applied to te on te difference mes and ψ prove Q te Poincaré inequality. outk p(ψ in ). Tis is te difference between te outflow values of ψ and its K Assumption inflow values 3.. Define transported te layer ontowidt out K. bydiscrete analogues of te Poincaré inequality bound L 2 -like norms using sum of squares of differences (under some discrete analogue of a zero boundary condition). Terefore, d l =max{ it is peraps β,k : K not S surprising l }. tat suc a discrete Poincaré We assume type estimate tat our can meses also be arefound suc for tate norm in (3) (recalling tat any ψ in M satisfies a zero boundary condition on te inflow n boundary). To establis te inequality, we first dneed l Lto Ω define layers of mes elements marcing from in Ω. Te first layer is definedl=0 by were L Ω is a fixed constant depending on Ω, but independent of te meses and element sizes. We also assume tat tere iss a 0 fixed = {Kconstant T : C in K suc in tat Ω}. if K and K are any two
14 4 L. DEMKOWICZ AND J. GOPALAKRISHNAN By Lemma 3., S 0 is not empty. Next, setting T (0) integers l > 0, te sets Λ l T (l) = T (l ) \ S l S l = {K T (l) : in K in Ω (l) }, = T, we recursively define, for all were Ω (l) is te domain formed by te union of all triangles of T (l). Repeated application of Lemma 3. sows tat S l is nonempty for every l, unless T (l) is empty. So tis process exausts all elements of te mes, at some finite value of l, say l = n. Te domain def = Ω \ Ω (l) is mesed by S 0 S S l. Its outflow boundary is denoted by Γ l = out Λ l. All mes edges, excluding tose on in Ω, and excluding tose were β n = 0, are contained in Γ l for some l n. Wit tese notations and observations, we impose a mild assumption on te mes and prove te Poincaré inequality. Assumption 3.. Define te layer widt by d l = max{ β,k : K S l }. We assume tat our meses are suc tat n d l L Ω l=0 were L Ω is a fixed constant depending on Ω, but independent of te meses and element sizes. We also assume tat tere is a fixed constant C suc tat if K and K are any two neigboring elements, ten eiter te inequality olds or it olds after excanging K and K. β,k C β,k Above and in te remainder, te letter C, wit or witout subscripts, denotes generic constants independent of elements K and polynomial degree p. Teir value at different occurrences may vary. Assumption 3. permits quite general meses, including anisotropic refinements (e.g., to capture sock waves). Loosely speaking, te first part assumes tat elements witin eac layer ave about te same lengt in te flow direction. Te second part assumes tat element lengts in te flow direction for neigboring elements are comparable. Quasiuniform meses obviously satisfy te assumption. Teorem 3.2. If Assumption 3. olds, ten for all ψ M, we ave ψ 2, β,t C 0L 2 Ω ψ 2, β,t. (were C 0 is independent of T and p, and te norms are as in (29) and (3)). Proof. Te first step is a local inequality. On any triangle K, ψ β, outk ψ Q pψ β, outk + Q pψ β, outk ψ Q p ψ β, outk + ψ β, ink,
15 DPG METHOD 5 were te last inequality was due to Teorem 2.2 (2). Consequently, ψ 2 β, outk ψ 2 β, ink ( ψ β, outk + ψ ) (34) β, ink ψ Qp ψ β, outk olds on every mes element K. Te next step is to sum over R l = S 0 S S l. Ten, we ave ) K R l ( ψ 2 β, outk ψ 2 β, ink K R l ( ψ β, outk + ψ β, ink) ψ Qp ψ β, outk. Wen rewriting te left and side as a sum over te mes edges, we observe tat all contributions due to edges interior to Λ l cancel out. Furtermore, since ψ M, it vanises on te inflow boundary in Λ l in Ω. Terefore only contributions on te outflow boundary of Λ l, namely Γ l, remains, i.e., ψ 2 (35) β,γ l ( ψ β, outk + ψ ) β, ink ψ Qp ψ, outk. β K R l We apply Caucy-Scwarz inequality to te rigt and side of (35) to get ( ψ 2 ( 2 β,γ l β,k ψ 2 + ) ) /2 β, outk ψ 2 β, ink K R l ( ) /2 ψ Q p ψ 2 β,k β, outk K R l By Assumption 3., te lengts β,k for neigboring elements are comparable, so tat wic implies K R l β,k ψ 2 β, ink C K R l β,k ψ 2 β, outk ( ) /2 ψ 2 C β,γ l β,k ψ 2 β, outk K R l ( ) /2 ψ Q p ψ 2 β,k β, outk K R l l=0 = C ψ, β,r ψ l, β,r. l Te final step involves multiplying by te layer widt d l, and summing over l, n n d l ψ 2 d β,γ l l C ψ, β,r ψ l, β,r l l=0 CL Ω ψ, β,t ψ, β,t, (36)
16 6 L. DEMKOWICZ AND J. GOPALAKRISHNAN were we increased te norms over R l to norms over te wole mes T. By te definition of d l, we know tat so we ave β,k d l, for all K S l, ψ 2, β,t = n l=0 β,k ψ 2 β, outk K S l n d l ψ 2 β, outk l=0 K S l L Ω C ψ, β,t ψ, β,t, were te last inequality is due to (36). Tis proves te result L 2 -error estimates. Now we prove an optimal error estimate for u in te standard L 2 (Ω)-norm, as well as error estimates for te outflux approximation φ in L 2 -like norms on mes edges. Te estimates are wit constants independent of te polynomial degree p and mes T. Teorem 3.3. Suppose (u, φ) is te exact solution, (u, φ ) is te discrete solution. Ten we ave te following error estimates: If Assumption 3. olds, letting we ave φ φ 2 3 inf φ ψ, β,k p+ 2 ψ p+ M β,k β ( ) ɛ(φ) def = inf φ ψ 2 + C ψ M, 0L 2 β,t Ω φ ψ 2,, β,t, K. (37) φ φ 2 2 ɛ(φ) and (38), β,t u u 2 Ω 4 β ɛ(φ) + 2 inf u w 2 w W Ω. (39) Proof. Te first step of tis proof consists of identifying te equations satisfied by te following discrete error functions: ε u = Π W u u, ε φ = Π Mφ φ,
17 DPG METHOD 7 were Π W and Π M are L 2 -ortogonal projections into W and M, respectively. Subtracting te exact and discrete equations, 0 = K T ( (u u, β v ) K ) + φ φ, v outk φ φ, v in K = K T ( (Π W u u, β v ) K + Π M φ φ, v outk φ φ, v in K were te introduction of Π W is possible as β maps V p (K) into P p (K). Notice tat in te last term we may not replace φ by Π M φ as te inflow traces of v need not be polynomial. But we may add and subtract Π M φ witin it. Ten, we can express te identity using te discrete error functions as follows: or equivalently (see (33)), a ( (ε u, ε φ ), v ) = K T φ Π M φ, v in K, ), on every mes element K. In oter words, (ε u U p ε φ, β v ) K + ε φ Q pε φ, v outk = φ Π M φ, v in K, ε u U p ε φ = U p(φ Π M φ), ε φ Q pε φ = Q p(φ Π M φ). We sall use tese identities to bound te errors. Te second step of tis proof involves bounding ε φ. By (4b), (40) (4a) (4b) ε φ Q pε φ β, outk = Q p(φ Π M φ) β, outk φ Π M φ β, ink, (42) by te bound (2) on Q p ( ) of Teorem 2.2. Terefore, φ φ, β,k φ Π M φ, β,k + εφ, β,k = /2 β,k (φ Π Mφ) Q p (φ Π M φ) β, outk + /2 β,k εφ Q pε φ β, outk /2 β,k φ Π Mφ β, outk + 2 /2 β,k φ Π Mφ β, ink. By te inequality of aritmetic and geometric means, tis implies φ φ 2, 3,K β,k φ Π Mφ 2, β β, outk ink
18 8 L. DEMKOWICZ AND J. GOPALAKRISHNAN wic gives te first error estimate of te teorem. Te tird step again involves a bound on ε φ, but now in a weaker norm. We use te Poincaré inequality of Teorem 3.2: ε φ 2 C, 0L 2 β,t Ω ε φ 2, β,t C 0 L 2 Ω φ Π M φ 2 β, in K, (43) β,k K T were we ave also used (42) in te last step. Splitting φ φ = (φ Π M φ) + ε φ as before, and estimating, we prove te second estimate (38) of te teorem. Now we come to te final step, were we obtain te error estimate for u, namely (39), as follows: ε u K U p ε φ K ε u U p ε φ K by triangle inequality = U p (φ Π M φ) K by (4a) /2 β,k β /2 φ Π Mφ β, ink by Teorem 2.2. Using Teorem 2.2 again, we ave ε u K /2 β,k ( ε φ β /2 β, ink + φ Π ) Mφ β, ink. Summing over all elements and using (43), we obtain ε u 2 Ω 2 ( ε φ β 2 +, β,t β,k φ Π M φ 2 ) β, in K K T 2 β K T ( C0 L 2 Ω β,k φ Π M φ 2 β, in K + β,k φ Π M φ 2 β, in K ), from wic te estimate (39) follows Interior fluxes by postprocessing. Te advective flux is q = βu. A natural approximation to tis is q = βu, were u is te computed solution. However q is not conservative. It is possible to easily adapt an idea of [8] to generate a conservative flux by a simple postprocessing sceme. Our postprocessed flux is denoted by q. Its restriction to eac element K lies in te Raviart-Tomas space P p+ (K) + xp p+ (K) and is defined by ( q, r) K = ( β u, r) K, r P p (K), (44a) q n, µ F = ϕ, µ n F, µ P p+ (F ), (44b)
19 DPG METHOD 9 were (44b) is imposed on all edges F of K, φ n, on out K, φ n, on in K \ in Ω, ϕ = β g, on in K in Ω, 0, on faces were β n = 0, and, as before, n denotes te unit outward normal of K. Since (44) uses te well known degrees of freedom of te Raviart-Tomas space, it is easy to see tat (44) uniquely defines a q in H(div, Ω). Teorem 3.4. Write u β def = β u = q and u β, def = q. Ten, q is conservative, and u β, Π p+ u β = 0 (45) were Π p+ is te L 2 (K)-ortogonal projection on P p+ (K). Proof. First of all, observe tat P p+ (K) V p (K). Indeed, if z p+ is any polynomial in P p+ (K), ten finding a v in V p (K) as in Lemma 2. suc tat β v = β z p+, (46a) v outk = z p+ outk, (46b) we find from (46a) tat z p+ v must be a function tat is constant along te β-direction. But te same function vanises on out K by (46b), so it must vanis everywere. Hence z p+ v V p (K). As a consequence, we ave, for all z in P p+ (K), ( q, z) K + q n, z K = (u, β z) K + φ, z outk φ, z in K + β n g, z in K in Ω by (44), = (f, z) K = ( β u, z) K by (25). Wit q = βu, we can rewrite tis as ( q, z) K = ( q, z), for all z P p+ (K), (47) wic proves (45). To sow tat q is conservative, let D be any subdomain formed by a union of some or all of te mes elements. Ten (47) applied wit z = implies q n ds = q n ds, D were we ave also used te H(div, Ω)-conformity of q. Tis sows tat te net outward fluxes of q and q coincide on any suc subdomain D. D 4. Numerical studies We present a few small scale numerical studies to illustrate and confirm te teory. For more realistic and larger size simulations, we will need full p adaptivity, wic will be presented elsewere.
20 20 L. DEMKOWICZ AND J. GOPALAKRISHNAN 4.. Implementation aspects. In comparison wit te UDG metod, we do ave more equations to solve. However, te reason our metod remains competitive is tat we are able to solve for te fluxes (φ ) first and locally recover u afterwards. Tis results in a dimensional reduction tat is quite attractive for ig p. Indeed, wit respect to p, te number of degrees of freedom for φ is O(p), wile te corresponding number for u is O(p 2 ). Tis dimensional reduction is quite analogous to ybridized metods [9], were te so-called Lagrange multipliers are solved for first, and te remaining variables are recovered locally, element by element. To state te matrix form of te metod, we need bases for W, M, and V. For M, we set a basis consisting of functions supported only on one edge, wile te bases for te oter spaces consist of functions supported on only one triangle. Let Φ and U denote te vector of coefficients of φ and u in teir basis expansion, respectively. To sow ow one solves for fluxes first, observe tat te global finite element space V can be split into te outflow extensions E out (M ) plus a linearly independent remainder R. We enumerate te global basis for V suc tat te basis functions in E out (M ) come first. Wit v set to suc test functions, (25) reduces to an equation involving only φ : [ ] φ, v outk φ, v in K\ in Ω (48) K T = (f, v ) Ω β n g, v in Ω for all v E out (M ). As a consequence, te stiffness matrix of (25) is block triangular: ( ( ) A 0 Φ = F (49) B C) U were A, B and C are sparse matrix blocks, and F is te corresponding load vector. Tis clearly sows tat we can solve for Φ first and te result can ten be used to compute u (or U) locally, element by element. Additionally, note tat in te decomposition V = E out (M ) R, all te functions in V aving discontinuities inside elements lie in te E out (M )-component. Terefore, wen making te stiffness matrix of te bilinear form a ( (u, φ ), v ), tere is no need to integrate basis functions wit discontinuities witin te element. Indeed, as seen in (48), te term involving integration witin elements, namely (u, β v ), vanises for v in E out (M ), so tere is no need to compute tis integral. Finally, suppose we enumerate te basis for M in suc a way tat functions on Γ l appear before tose on Γ l (see te definition of te layers in 3.4). Ten te matrix A in (49) is a square triangular matrix. In tis case, Φ can be found fast using backsubstitution. Since tis is te only globally coupled system to solve, te remaining (local) computations for u can be locally optimized (or parallelized) A one-dimensional example. We begin wit a simple one-dimensional example to sow te enanced accuracy and stability of te DPG metod compared to te DG metod. Te model problem on Ω = (0, ) tat we sall consider is u = f, u(0) = u 0, wit data u 0 R and f L 2 (0, ) set so tat te exact solution is u(x) = arctan(α(x x 0 ))
21 DPG METHOD DPG DG L 2 error in u p Figure 4. Convergence for te one-dimensional example were x 0 = and α = 00. We use te analogue of our spectral DPG metod (of Section 2) in one dimension, i.e., K = (0, ) and te test space is V p (K) = P p+ (K). Figure 4 sows te error in u vs. p in a semilog plot. Wile bot metods exibit exponential convergence, te error for te DPG metod is one order less tan tat for te DG metod. Te gap between te two convergence curves grows wit increasing coefficient α tat controls te steepness of te solution. Te different beavior of te two metods is furter illustrated in Figures 5(a) and 5(b) tat sow te graps of te approximate solutions for some values of p. Bot metods deliver te exact value of te flux at x =. Contrary to te DG metod, owever, te DPG solution u witin te element is disconnected from value at and delivers te best approximation error in L 2 -norm. Tis manifests itself wit a less oscillatory beavior of te solution. Our numerics also confirmed Teorem 2., i.e., wen we compared te DPG solution wit te L 2 -projection of te exact solution, we found te difference to be zero (of te order of round-off errors) and p convergence to a solution wit discontinuity. Following [2], we consider a square domain Ω = (, ) 2 werein an advection problem is set up so tat β = (, 9/0) and te exact solution is sin(π(x + ) 2 /4 sin(π(y 9x/0)/2) for x, 9x/0x < y <, u(x, y) = e 5(x2 +(y 9x/0) 2 ) for x, y < 9x/0. Tis u as a discontinuity along te line y = 9x/0. As in [2], we first mes Ω by an n n quadrilateral mes aligned wit te line of discontinuity (see [2, Fig. 0] for details). Since our metod is based on triangular meses ([2] only considered quadrilateral meses), we furter split eac of te quadrilaterals into two triangles by connecting teir diagonals of positive slope. Anoter difference wit te experiment in [2] is tat wile [2] needed to include a stabilizing lower order reaction term for teoretical reasons, we do not need to,
22 22 L. DEMKOWICZ AND J. GOPALAKRISHNAN 0 Te spectral DPG solutions 0.5 u(x).5 2 Exact solution 2.5 p= p=3 p= x 0 (a) Te DPG solution Te spectral DG solutions 0.5 u(x).5 2 Exact solution 2.5 p= p=3 p= x (b) Te DG solution Figure 5. Te enanced stability of DPG approximations, compared to te DG ones, in te case of a model transport problem wit a layer at te rigt end ence we solve te purely advective problem. We report te errors in u and φ for te cases n = 5, 9, 7, 33 and p =,..., 5 in Figure 6. It is clear from Figure 6 tat te DPG metod converged under and p refinements as if te solution were infinitely regular (see Figure 6). We observed a similar beavior for te DG metod as well. Suc observations demonstrate te advantage of using discontinuous finite elements wit meses aligned wit te solution discontinuity. For te standard DG metod, tis observation was made in [2], were DG is compared wit te streamline diffusion metod. Wile te streamline diffusion solution can also be improved by aligning meses wit sock lines [2], its convergence rate remains limited by te regularity of te
23 DPG METHOD 23 (a) Log-log plots of u u for various and p values (b) Log-log plots of flux errors φ φ β,e Figure 6. Te p convergence of te DPG metod (log-log plots) for te Houston-Scwab-Süli example [2] solution, as observed in [2]. Since our purpose is a comparison of DG and DPG, we note ere tat te errors we observed for te DG metod (not sown in te figure), as in 4.2, remained iger tan tat of te DPG metod. Figure 6 suggests tat exponential rates of convergence can be obtained wit te DPG metod, even for discontinuous solutions, once an adaptive strategy to align te meses wit te socks is implemented. Te data tat was used to plot Figure 6 also sows te rate of te convergence. Figure 6(a) indicates tat u u converges at te rate p+ (for p =... 5). Tis is in accordance wit Teorem 3.3. On te oter and, Figure 6(b) sows tat te fluxes converge suc tat φ φ, β,t goes to zero at te rate of p+2.5, wic is one order more tan wat we were able to prove in Teorem 3.3.
24 24 L. DEMKOWICZ AND J. GOPALAKRISHNAN (a) -refinement: DPG converges at a iger rate (b) p-refinement: DPG gives lower error Figure 7. Te DPG metod performs better tan te DG metod under bot and p refinements for te Peterson example [7]. Bot plots are in a log-log scale Te Peterson example. Next, we consider te well-known Peterson example [7]. Peterson s mes of te unit square is a specific quasiuniform mes of te type appearing in Figure 8, obtained from an n n partition of te unit square, but wit additional vertical lines troug some mes vertices suc tat te new lines divide te unit square into, say m, vertical strips (see [7] for details). On suc meses, we apply te UDG and DPG metod to te problem u/ y = 0, u(x, 0) = sin 6x, x (0, ). Te results in Figure 7 sow tat te DPG metod delivers better results tan te standard DG metod. Figure 7(a) is obtained wit p = and four progressively finer meses: first wit n = 6, m = 3, second
25 DPG METHOD 25 (a) Te DPG solution (b) Te DG solution Figure 8. DPG solutions exibit less crosswind diffusion tan DG solutions. In tis example, te flow is upward. wit n = 2, m = 6, tird wit n = 24, m = 8, and fourt wit n = 48, m = 6. Tis allows us to study te case of -refinement. On te oter and, Figure 7(b) is obtained by fixing te mes wit n = 6, m = 3, and increasing te degree p from 0 troug 5. Te results clearly demonstrate te well-known 0.5-order suboptimality of te DG metod in te -refinement case, but more interestingly, also sow tat te DPG metod converges at te optimal rate. We also present te solutions u obtained wit bot DG and DPG metods wen te inflow data is u(x, 0) = x 2 in Fig. 8. We used piecewise constant approximations (p = 0). Observe tat as te flow proceeds upwards, te values diffuse orizontally for DG case, wile tere is little indication of suc crosswind diffusion for te DPG metod wit te optimal test functions. 5. Concluding remarks 5.. Summary. We proved optimal error estimates for u in and p. Using te known best approximation estimates [20] for a quasiuniform mes of messize, we obtain as a corollary of Teorem 3.3 tat C u u Ω s p u s H s (Ω) + ( )s2 /2 β /2 p u s H 2 /2 s 2 (Ω) (50) wit some C independent of and p and any 0 s p + and 0 s 2 p + 2. If te solution is smoot, coosing s = p+ and s 2 = p+2, we obtain te optimal -convergence rate of O( p+ ). Te p-convergence rate is also optimal. Note owever tat tere may be room for improvement in te tecniques used in te current proof. In Teorem 3.3 and (50), note tat te regularity on u tat we need to obtain te optimal convergence rate is iger tan expected. Furtermore, for te fluxes, we proved a rate of convergence tat is suboptimal by one order, but numerical experiments
26 26 L. DEMKOWICZ AND J. GOPALAKRISHNAN DPG METHOD 23 Inflow faces witout data Foreground points Background points Foreground edges Edges covered by a face Figure 8. Top view of a 3D mes. To te 3D mes, keep te background 9. points Toponview te plane of aof 3Dte mes. paper To andobtain raise te te foreground 3D mes, points keep onete back- Figure ground unit from points te paper. on tete plane inflow of te direction paper is emanating and raise from te foreground te paper. points one unit from te paper. Te inflow direction is emanating from te paper. to formulate a generalization of our metod by combining wit temixedmetod. Tis is because te space for u in our metod and in te mixed metod is te same (one can indicate proceed, tat forte instance, fluxesalong do converge te same lines at optimal as in [6]). order. However, A sarper provinganalysis p estimates sedding wit ligt on tese p independent issues can constants be valuable. for suc a metod does not appear to be an easy task. Te case of variable transport direction β can be andled by projecting β into an appropriate finite element space. In particular, a low order metod canbeimmediatelywritten more accurate solutions tan In every numerical experiment we performed, witout exception, te DPG metod gave down if we approximate β te standard DG metod. Numerical experiments by its lowest-order Raviart-Tomas interpolant Π RT β. confirm Ten, te wenever teoretically β =0,teapproximantΠ proved optimal rate of RT β convergence is a piecewise constant for te vector primal field variable wit continuous course, interelement te full normal potential components. of te Hence, metod by simply is not modifying indicated our test by space an estimate V p (K) to like (50). u. Of Indeed, align one witcan Π RT β acieve witin eac exponential element K, rate weobtainalow-ordergeneralizationoftedpg of convergence by constructing a variable-degree analogue metodwere to te case eacofelement variable β Kwit is endowed minimal modifications. wit its ownto trial obtain space truly of ig-order degree p K. Te fluxes metods mustfor bevariable cosenβ,ourtestspacemusttakeintoaccounttevariableadvection. to be of degree one plus te maximum of te degrees from adjacent elements. Tis and A full oter study generalizations of te variable becomedegree clearer, version, once we view togeter te metod wit an presented p-adaptive ere sceme wic as aautomatically particular caseselects of a more orgeneral p refinement paradigminfor different constructing locations, scemes, is postponed not only forto a future te transport problem or its singular perturbation, but indeed for a muc larger class work. of problems. Te new paradigm is to find a space of optimal test functions tailored for any particular problem. It turns out tat te framework of DG metods offers a unique opportunity to locally compute test functions tat approac optimality in te sense of 5.2. reproducing Te tree-dimensional te exact stability properties case. Te of te entire boundary analysis value of problem Sectionat 2 te olds discrete in any dimension, level. including Tis approac tree. and Teitsole connection reasonto wete restricted space of test ourselves functions topresented te two-dimensional ere will case in Section be clarified 3 is intat a sequel. Lemma 3. does not old for general tetraedral meses. Tis fact does not seem to be well known, so we indicate a counterexample in Figure 9. In te mes sown terein, tere is no Appendix tetraedron A. Proof wose ofinflow Teorem boundary 2.2 is wolly contained in te global inflow boundary of te domain. Tis section is devoted to a proof of Teorem 2.2. Before we embark on tis, we need Intosituations develop several like intermediate in Figure 9, results. we may Recall coose tat Eto out appropriately denotes te extension bisect from te problematic te tetraedra outflow boundary. so tat Now, Lemma let us3. alsodoes defineold te extension at eacfrom stage. inflowalternately, boundary, denotedby wile generating tetraedral E in as follows. meses Given byaadvancing function µ front in K,tefunctionE mesing algoritms, in µ on Kwe coincides can ensure wit µ tat on elements at eac in K and front is independent are suc tat of η Lemma (i.e., it is 3. constant is satisfied. along te β-direction). Notice tat Lemma 3. is required not only for te analysis of our metod, but also for practical implementation. Indeed, te process of splitting te mes into layers S l (see 3.4) yields an enumeration of flux degrees of freedom tat makes te matrix for φ s triangular, and ence suitable for fast solution by backsubstitution. Terefore any extra effort to obtain meses of te required type is likely to pay off in te final analysis. (Note tat Lemma 3. is needed to obtain triangular matrices even wen using te standard UDG metod.) Tese and related algoritmic issues will be studied in te future.
A Class of Discontinuous Petrov Galerkin Methods. Part I: The Transport Equation
Portland State University PDXScolar Matematics and Statistics Faculty Publications and Presentations Fariborz Masee Department of Matematics and Statistics 200 A Class of Discontinuous Petrov Galerkin
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 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 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 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 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 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 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 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 informationCONSTRUCTIVELY WELL-POSED APPROXIMATION METHODS WITH UNITY INF SUP AND CONTINUITY CONSTANTS FOR PARTIAL DIFFERENTIAL EQUATIONS
MATHEMATICS OF COMPUTATION Volume 82, Number 284, October 203, Pages 923 952 S 0025-578(203)02697-X Article electronically publised on April 23, 203 CONSTRUCTIVELY WELL-POSED APPROXIMATION METHODS WITH
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 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 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 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 informationArbitrary order exactly divergence-free central discontinuous Galerkin methods for ideal MHD equations
Arbitrary order exactly divergence-free central discontinuous Galerkin metods for ideal MHD equations Fengyan Li, Liwei Xu Department of Matematical Sciences, Rensselaer Polytecnic Institute, Troy, NY
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 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 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 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 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 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 informationCopyright c 2008 Kevin Long
Lecture 4 Numerical solution of initial value problems Te metods you ve learned so far ave obtained closed-form solutions to initial value problems. A closedform solution is an explicit algebriac formula
More informationNew DPG techniques for designing numerical schemes
New DPG techniques for designing numerical schemes Jay Gopalakrishnan University of Florida Collaborator: Leszek Demkowicz October 2009 Massachusetts Institute of Technology, Boston Thanks: NSF Jay Gopalakrishnan
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 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 informationThe Laplace equation, cylindrically or spherically symmetric case
Numerisce Metoden II, 7 4, und Übungen, 7 5 Course Notes, Summer Term 7 Some material and exercises Te Laplace equation, cylindrically or sperically symmetric case Electric and gravitational potential,
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 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 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 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 informationA MULTILEVEL PRECONDITIONER FOR THE INTERIOR PENALTY DISCONTINUOUS GALERKIN METHOD
A MULTILEVEL PRECONDITIONER FOR THE INTERIOR PENALTY DISCONTINUOUS GALERKIN METHOD KOLJA BRIX, MARTIN CAMPOS PINTO, AND WOLFGANG DAHMEN Abstract. In tis article we present a multilevel preconditioner for
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 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] 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 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 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 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 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 informationDiscontinuous Galerkin Methods for Relativistic Vlasov-Maxwell System
Discontinuous Galerkin Metods for Relativistic Vlasov-Maxwell System He Yang and Fengyan Li December 1, 16 Abstract e relativistic Vlasov-Maxwell (RVM) system is a kinetic model tat describes te dynamics
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 informationAdaptive Finite Element Method
39 Capter 3 Adaptive inite Element Metod 31 Introduction As already pointed out in capter 2, singularities occur in interface problems Wen discretizing te problem (221) wit inite Elements, te singularities
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 informationarxiv: v1 [math.na] 20 Nov 2018
An HDG Metod for Tangential Boundary Control of Stokes Equations I: Hig Regularity Wei Gong Weiwei Hu Mariano Mateos Jon R. Singler Yangwen Zang arxiv:1811.08522v1 [mat.na] 20 Nov 2018 November 22, 2018
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 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 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 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 informationAnalytic Functions. Differentiable Functions of a Complex Variable
Analytic Functions Differentiable Functions of a Complex Variable In tis capter, we sall generalize te ideas for polynomials power series of a complex variable we developed in te previous capter to general
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 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 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 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 informationHybridization and Postprocessing Techniques for Mixed Eigenfunctions
Portland State University PDXScolar Matematics and Statistics Faculty Publications and Presentations Fariborz Masee Department of Matematics and Statistics 2010 Hybridization and Postprocessing Tecniques
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 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 informationOn convexity of polynomial paths and generalized majorizations
On convexity of polynomial pats and generalized majorizations Marija Dodig Centro de Estruturas Lineares e Combinatórias, CELC, Universidade de Lisboa, Av. Prof. Gama Pinto 2, 1649-003 Lisboa, Portugal
More informationHOMEWORK HELP 2 FOR MATH 151
HOMEWORK HELP 2 FOR MATH 151 Here we go; te second round of omework elp. If tere are oters you would like to see, let me know! 2.4, 43 and 44 At wat points are te functions f(x) and g(x) = xf(x)continuous,
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 informationFourier Type Super Convergence Study on DDGIC and Symmetric DDG Methods
DOI 0.007/s095-07-048- Fourier Type Super Convergence Study on DDGIC and Symmetric DDG Metods Mengping Zang Jue Yan Received: 7 December 06 / Revised: 7 April 07 / Accepted: April 07 Springer Science+Business
More informationA trace finite element method for a class of coupled bulk-interface transport problems
Numerical Analysis and Scientific Computing Preprint Seria A trace finite element metod for a class of coupled bulk-interface transport problems S. Gross M.A. Olsanskii A. Reusken Preprint #28 Department
More informationEXTENSION OF A POSTPROCESSING TECHNIQUE FOR THE DISCONTINUOUS GALERKIN METHOD FOR HYPERBOLIC EQUATIONS WITH APPLICATION TO AN AEROACOUSTIC PROBLEM
SIAM J. SCI. COMPUT. Vol. 26, No. 3, pp. 821 843 c 2005 Society for Industrial and Applied Matematics ETENSION OF A POSTPROCESSING TECHNIQUE FOR THE DISCONTINUOUS GALERKIN METHOD FOR HYPERBOLIC EQUATIONS
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 informationCombining functions: algebraic methods
Combining functions: algebraic metods Functions can be added, subtracted, multiplied, divided, and raised to a power, just like numbers or algebra expressions. If f(x) = x 2 and g(x) = x + 2, clearly f(x)
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 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 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 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 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 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 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 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 information5.1 We will begin this section with the definition of a rational expression. We
Basic Properties and Reducing to Lowest Terms 5.1 We will begin tis section wit te definition of a rational epression. We will ten state te two basic properties associated wit rational epressions and go
More informationch (for some fixed positive number c) reaching c
GSTF Journal of Matematics Statistics and Operations Researc (JMSOR) Vol. No. September 05 DOI 0.60/s4086-05-000-z Nonlinear Piecewise-defined Difference Equations wit Reciprocal and Cubic Terms Ramadan
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 informationarxiv: v1 [math.na] 9 Mar 2018
A simple embedded discrete fracture-matrix model for a coupled flow and transport problem in porous media Lars H. Odsæter a,, Trond Kvamsdal a, Mats G. Larson b a Department of Matematical Sciences, NTNU
More informationOne-Sided Position-Dependent Smoothness-Increasing Accuracy-Conserving (SIAC) Filtering Over Uniform and Non-uniform Meshes
DOI 10.1007/s10915-014-9946-6 One-Sided Position-Dependent Smootness-Increasing Accuracy-Conserving (SIAC) Filtering Over Uniform and Non-uniform Meses JenniferK.Ryan Xiaozou Li Robert M. Kirby Kees Vuik
More informationUniversity Mathematics 2
University Matematics 2 1 Differentiability In tis section, we discuss te differentiability of functions. Definition 1.1 Differentiable function). Let f) be a function. We say tat f is differentiable at
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 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 informationExcursions in Computing Science: Week v Milli-micro-nano-..math Part II
Excursions in Computing Science: Week v Milli-micro-nano-..mat Part II T. H. Merrett McGill University, Montreal, Canada June, 5 I. Prefatory Notes. Cube root of 8. Almost every calculator as a square-root
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 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 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 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 informationALGEBRA AND TRIGONOMETRY REVIEW by Dr TEBOU, FIU. A. Fundamental identities Throughout this section, a and b denotes arbitrary real numbers.
ALGEBRA AND TRIGONOMETRY REVIEW by Dr TEBOU, FIU A. Fundamental identities Trougout tis section, a and b denotes arbitrary real numbers. i) Square of a sum: (a+b) =a +ab+b ii) Square of a difference: (a-b)
More informationSymmetry Labeling of Molecular Energies
Capter 7. Symmetry Labeling of Molecular Energies Notes: Most of te material presented in tis capter is taken from Bunker and Jensen 1998, Cap. 6, and Bunker and Jensen 2005, Cap. 7. 7.1 Hamiltonian Symmetry
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 informationAdvancements In Finite Element Methods For Newtonian And Non-Newtonian Flows
Clemson University TigerPrints All Dissertations Dissertations 8-3 Advancements In Finite Element Metods For Newtonian And Non-Newtonian Flows Keit Galvin Clemson University, kjgalvi@clemson.edu Follow
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 informationFlavius Guiaş. X(t + h) = X(t) + F (X(s)) ds.
Numerical solvers for large systems of ordinary differential equations based on te stocastic direct simulation metod improved by te and Runge Kutta principles Flavius Guiaş Abstract We present a numerical
More informationLecture 10: Carnot theorem
ecture 0: Carnot teorem Feb 7, 005 Equivalence of Kelvin and Clausius formulations ast time we learned tat te Second aw can be formulated in two ways. e Kelvin formulation: No process is possible wose
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 informationExercises for numerical differentiation. Øyvind Ryan
Exercises for numerical differentiation Øyvind Ryan February 25, 2013 1. Mark eac of te following statements as true or false. a. Wen we use te approximation f (a) (f (a +) f (a))/ on a computer, we can
More informationHigher order unfitted isoparametric space-time FEM on moving domains
Higer order unfitted isoparametric space-time FEM on moving domains Master s tesis by: Janosc Preuß Supervisor: Jun.-Prof. Dr. Cristop Lerenfeld Second Assessor: Prof. Dr. Gert Lube prepared at te Institute
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 informationNONNEGATIVITY OF EXACT AND NUMERICAL SOLUTIONS OF SOME CHEMOTACTIC MODELS
NONNEGATIVITY OF EXACT AND NUMERICAL SOLUTIONS OF SOME CHEMOTACTIC MODELS PATRICK DE LEENHEER, JAY GOPALAKRISHNAN, AND ERICA ZUHR Abstract. We investigate nonnegativity of exact and numerical solutions
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 informationA Local Projection Stabilization/Continuous Galerkin Petrov Method for Incompressible Flow Problems
A Local Projection Stabilization/Continuous Galerkin Petrov Metod for Incompressible Flow Problems Naveed Amed, Volker Jon, Gunar Matties, Julia Novo Marc 22, 218 Abstract A local projection stabilization
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 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 information