A spacetime DPG method for acoustic waves
|
|
- William Horn
- 5 years ago
- Views:
Transcription
1 A spacetime DPG method for acoustic waves Rainer Schneckenleitner JKU Linz January 30, 2018 ( JKU Linz ) Spacetime DPG method January 30, / 41
2 Overview 1 The transient wave problem 2 The broken weak formulation 3 Verification of the density condition 4 The ideal DPG method and a priori & a posteriori error estimates 5 Implementation of practical DPG and numerical results ( JKU Linz ) Spacetime DPG method January 30, / 41
3 Outline 1 The transient wave problem 2 The broken weak formulation 3 Verification of the density condition 4 The ideal DPG method and a priori & a posteriori error estimates 5 Implementation of practical DPG and numerical results ( JKU Linz ) Spacetime DPG method January 30, / 41
4 Problem formulation The wave equation Let Ω 0 R d be the spatial domain with boundary Ω 0 and let Ω = Ω 0 (0, T ) be the spacetime cylinder, with T > 0. Furthermore, let f L 2 (Ω) and g L 2 (Ω) d. Then the considered first order system for the wave equation is given by t q grad x µ = g, t µ div x q = f. Additionally, the wave equation is equipped with homogeneous initial and boundary conditions, i.e., µ t=0 = 0, q t=0 = 0, µ Ω0 (0,T ) = 0. ( JKU Linz ) Spacetime DPG method January 30, / 41
5 The formal wave operator Wave operator can be seen as first order distributional derivative operator A : L 2 (Ω) d+1 D (Ω) d+1 [ ] t u Au = q grad x u µ, t u µ div x u q where u L 2 (Ω) d+1 is split into u = [ uq u µ ] with u q L 2 (Ω) d and u µ L 2 (Ω). ( JKU Linz ) Spacetime DPG method January 30, / 41
6 The formal wave operator Theorem Let Ω be open. The space endowed with the norm is a Hilbert space. Proof. Blackboard. W (Ω) := {u L 2 (Ω) d+1 : Au L 2 (Ω) d+1 } = W u W = ( u 2 + Au 2 ) 1/2 ( JKU Linz ) Spacetime DPG method January 30, / 41
7 The formal wave operator The formal adjoint of A is the operator A and satisfies (Au, v) = (u, Av) for all u, v D(Ω) d+1. Furthermore, we introduce D : W W by Du, v W = (Au, v) + (u, Av) for all u, v W, where W denotes the dual space of W, and, W denotes the duality pairing in W. Assume u, v D(Ω) d+1, then Du, v W = u q (n t v q n x v µ ) + u µ (n t v µ n x v q )ds, Ω where n = (n x, n t ) is the outward unit normal to Ω R d+1 ( JKU Linz ) Spacetime DPG method January 30, / 41
8 The unbounded wave operator Definition of an unbounded operator, again denoted by A Domain dom(a) of A takes initial and boundary conditions into account We partition the spacetime boundary Ω into Moreover, we define Γ 0 = Ω 0 {0}, Γ T = Ω 0 {T }, Γ b = Ω 0 [0, T ], V = {v D(Ω) d+1 : v Γ0 = 0, v µ Γb = 0}, V = {v D(Ω) d+1 : v ΓT = 0, v µ Γb = 0}. ( JKU Linz ) Spacetime DPG method January 30, / 41
9 The unbounded wave operator Now we can define the unbounded operator with A : dom(a) L 2 (Ω) d+1 L 2 (Ω) d+1, [ ] t u Au = q grad x u µ, t u µ div x u q dom(a) = {u W : Du, v W = 0 for all v V }. D(Ω) d+1 dom(a) W = A is densly defined and has an adjoint A A equals A when applied to dom(a ) = {v L 2 (Ω) d+1 : l L 2 (Ω) d+1 s. t. (Au, v) = (u, l) for all u dom(a)} ( JKU Linz ) Spacetime DPG method January 30, / 41
10 The unbounded wave operator Throughout we denote dom(a) and dom(a ) endowed with the topology of W as V and V, respectively V and V are closed subsets of W = A, A are closed operators We have V = D(V ) := {w W : f, w W = 0 for all f D(V )} as well as the inclusions V V, V V. ( JKU Linz ) Spacetime DPG method January 30, / 41
11 Outline 1 The transient wave problem 2 The broken weak formulation 3 Verification of the density condition 4 The ideal DPG method and a priori & a posteriori error estimates 5 Implementation of practical DPG and numerical results ( JKU Linz ) Spacetime DPG method January 30, / 41
12 The broken weak formulation Partition of Ω into a mesh Ω h of open elements K with Ω = K Ω h K "Broken" analogue W h of W given by W h = {v L 2 (Ω) d+1 : A h v L 2 (Ω) d+1 }. Let A h be the wave operator applied element by element, i.e., (A h v) K = A(v K ), v W (K), K Ω h The operator D h : W h W h is defined by D h u, v h := D h u, v Wh = (A h u, v) + (u, A h v) for all u, v W h and, h denotes the duality pairing in W h ( JKU Linz ) Spacetime DPG method January 30, / 41
13 The broken weak formulation Let D h,v denote the restriction to V, i.e., and we define Q with is a complete space. q Q = D h,v = D h V Q := R(D h,v ) inf v D 1 h,v ({q}) v W ( JKU Linz ) Spacetime DPG method January 30, / 41
14 The broken weak formulation The bilinear form b((v, ρ), w) = (v, A h w) + ρ, w h on (L 2 (Ω) d+1 Q) W h leads to Broken weak formulation Let F W h. Find: u L2 (Ω) d+1 and λ Q such that b((u, λ), w) = F (w) for all w W h. (1) ( JKU Linz ) Spacetime DPG method January 30, / 41
15 The broken weak formulation Well-posedness of (1) is guaranteed, if V = D(V ), (2) A : V L 2 (Ω) d+1 is a bijection. (3) Theorem Suppose V is dense in V and V is dense in V. Then the conditions (2) and (3) are satisfied. Proof. Blackboard. ( JKU Linz ) Spacetime DPG method January 30, / 41
16 Outline 1 The transient wave problem 2 The broken weak formulation 3 Verification of the density condition 4 The ideal DPG method and a priori & a posteriori error estimates 5 Implementation of practical DPG and numerical results ( JKU Linz ) Spacetime DPG method January 30, / 41
17 Verification of the density condition Application of the density result for a hyperrectangle, i.e., for some a i > 0. Ω = Ω 0 (0, T ), Ω 0 = d (0, a i ), Theorem On the previously defined Ω, V is dense in V and V is dense in V. i=1 ( JKU Linz ) Spacetime DPG method January 30, / 41
18 Verification of the density condition Proof. The proof is divided into three steps. 1 Extension: Here, we extend a function in V using spatial reflections to a domain which has larger spatial extent than Ω. The operations R i, x = x 2x i e i, R i,+ x = x + 2(a i x i )e i perform reflections of x about x i = 0 and x i = a i for i = 1,..., d. The extended domains Q i are obtained recursively by Q 0 = Ω, Q i, = R 1 i, Q i 1, Q i,+ = R 1 i,+ Q i 1, Q i = Q i, Q i 1 Q i,+. The final extended domain is Q = Q d. ( JKU Linz ) Spacetime DPG method January 30, / 41
19 Verification of the density condition Then we need even and odd extensions of functions. Let G i,e, G i,o : L 2 (Q i 1 ) L 2 (Q i ) be defined by f (R i, x, t) if (x, t) Q i,, G i,e f (x, t) = f (R i,+ x, t) if (x, t) Q i,+, f (x, t) if (x, t) Q i 1. f (R i, x, t) if (x, t) Q i,, G i,o f (x, t) = f (R i,+ x, t) if (x, t) Q i,+, f (x, t) if (x, t) Q i 1. For vector valued functions v L 2 (Ω) d+1, we define G i v(x, t) = (G i,e v i )e i + j i (G i,o v j )e j. ( JKU Linz ) Spacetime DPG method January 30, / 41
20 Verification of the density condition Next, we define where E k = G k G k 1 G 1, E k = G k G k+1 G d, G i w(x, t) = (G i,ew i )e i + (G i,ow j )e j. j i with G G i,ow(x, t) = w(x, t) w(r 1 i, i,ew(x, t) = w(x, t) + w(r 1 i, x) w(r 1 i,+ x), x) + w(r 1 i,+ x). ( JKU Linz ) Spacetime DPG method January 30, / 41
21 Verification of the density condition It holds that (Ev, w) Q = (v, E w) for all v L 2 (Ω) d+1, w L 2 (Q) d+1. It can be proven that for any v V, AEv L 2 (Q) d+1, AEv coincides with EAv and Ev W (Q) ( JKU Linz ) Spacetime DPG method January 30, / 41
22 Verification of the density condition 2 Translation: In this step we translate up the previously obtained extension in time coordinate. Let v V and Ẽv be the extension of Ev by 0 to R d+1. For τ δ, δ > 0, the translation operator in time direction, i.e., (τ δ w)(x, t) = w(x, t δ) it holds that lim τ δg g L δ 0 2 (R d+1 ) = 0 for all g L2 (R d+1 ). With the restriction H δ from R d+1 to Q δ = d i=1 ( a i, 2a i ) ( δ, T + δ) it must be verified that AH δ τ δ Ẽv = H δ τ δ Ẽv. In particular we have H δ τ δ Ẽv W (Q δ ) whenever v V. ( JKU Linz ) Spacetime DPG method January 30, / 41
23 Verification of the density condition 3 Mollification: In this step we consider a v V and mollify the time-translated extension τ δ Ẽv. The used mollifier is given by where ρ 1 (x, t) = ρ ε (x, t) = ε (d+1) ρ 1 (ε 1 x, ε 1 t), { k exp( 1 1 x 2 t 2 ) if x 2 + t 2 < 1, 0 if x 2 + t 2 1 with k such that R d+1 ρ 1 = 1. To end this proof it suffices to show that is in V and v v ε Ω W ε 0 0. v ε = ρ ε τ δ Ẽv ( JKU Linz ) Spacetime DPG method January 30, / 41
24 Outline 1 The transient wave problem 2 The broken weak formulation 3 Verification of the density condition 4 The ideal DPG method and a priori & a posteriori error estimates 5 Implementation of practical DPG and numerical results ( JKU Linz ) Spacetime DPG method January 30, / 41
25 The ideal DPG method and a priori & a posteriori error estimates Approximation of the broken weak formulation by ideal DPG method Find: u h U h L 2 (Ω) d+1 and λ h Q h Q such that b((u h, λ h ), w h ) = F (w h ) for all w h T (U h Q h ), (4) where T : L 2 (Ω) d+1 Q W h is such that (T (v, ρ), w) h = b((v, ρ), w) for all w W h and any (v, ρ) L 2 (Ω) d+1 Q. The mixed formulation Find: ε h W h and (u h, λ h ) (U h Q h ) such that (ε h, w) h + b((u h, λ h ), w) = F (w) for all w W h, b((v, ρ), ε h ) = 0 is equivalent to formulation (4), see Seminar 08. for all (v, ρ) U h Q h ( JKU Linz ) Spacetime DPG method January 30, / 41
26 The ideal DPG method and a posteriori error estimates The expression η = ε h Wh = ε h 2 W (K) K Ω h 1/2 is an efficient and reliable a posteriori error estimator, see Seminar 08. ( JKU Linz ) Spacetime DPG method January 30, / 41
27 The ideal DPG method and a priori error estimates For a priori estimates we distinguish between Case A:: Ω h is a geometrically conforming mesh of (d+1)-simplices V h = {u V C(Ω) d+1 : u K P p+1 (K) d+1 for all K Ω h }, U h = {u L 2 (Ω) d+1 : u K P p (K) d+1 for all K Ω h }, Case B:: Ω h is a geometrically conforming mesh of hyperrectangles V h = {u V C(Ω) d+1 : u K Q p+1 (K) d+1 for all K Ω h }, U h = {u L 2 (Ω) d+1 : u K Q p (K) d+1 for all K Ω h }, where P p (K) and Q p (K) denote spaces of polynomials of total degree p, and of degree at most p in each variable, respectively. Q h = D h V h ( JKU Linz ) Spacetime DPG method January 30, / 41
28 The ideal DPG method and a priori error estimates Theorem Let u V H s+1 (Ω) d+1 and λ = D h u solve (1). Let U h and V h be as in one of the previously introduced cases depending on the mesh type, and Q h = D h V h. Then for (d 1)/2 < s p + 1. Proof. Blackboard. u u h + λ λ h Q ch s u H s+1 (Ω) d+1 ( JKU Linz ) Spacetime DPG method January 30, / 41
29 Outline 1 The transient wave problem 2 The broken weak formulation 3 Verification of the density condition 4 The ideal DPG method and a priori & a posteriori error estimates 5 Implementation of practical DPG and numerical results ( JKU Linz ) Spacetime DPG method January 30, / 41
30 Implementation of practical DPG and numerical results W h replaced by Y r h Case A:: Y r h = {w W h(ω) : w K P r (K) d+1 }, Case B:: Y r h = {w W h(ω) : w K Q r (K) d+1 }. Considered the mixed formulation Find: ε h Y r h and (u h, λ h ) (U h Q h ) such that (ε h, w) h + b((u h, λ h ), w) = F (w) for all w Y r h, b((v, ρ), ε h ) = 0 for all (v, ρ) U h Q h (5) r = p + d + 1 ( JKU Linz ) Spacetime DPG method January 30, / 41
31 Implementation of practical DPG and numerical results For the implementation of (5), λ h = D h z h for some z h V h and (ε h, w) h + b((u h, D h z h ), w) = F (w) for all w Y r h, b((v, D h r), ε h ) = 0 for all (v, r) U h V h (6) is considered. Decomposition of V h into V 0 h = {z V h : z K = 0 for all K Ω h } and remainder V 1 h = V h\v 0 h. b((v, D h V 0 h ), w) = 0, thus replace V h by V 1 h in (6). ( JKU Linz ) Spacetime DPG method January 30, / 41
32 Implementation of practical DPG and numerical results Yields the matrix equation [ A B B 0 ] [ e x ] = [ f 0 ], (7) where e and x are the vectors of coefficients in the basis expansion of ε h Y r h and (u h, z h ) U h V h, respectively, [A] kl = (y l, y k ) h, [B 0 ] ki = b((u i, 0), y k ), [B 1 ] = b((0, D h z j ), y k ) and B = [B 0, B 1 ]. r = p + d + 1 N (A) = N (B 0 ) = {0} B 1 may have a nontrivial kernel ( JKU Linz ) Spacetime DPG method January 30, / 41
33 Techniques to solve despite the null space Technique 1: Remaining orthogonal to null space in conjugate gradients Instead of (7) solve Schur complement system } B A {{ 1 B } x = B A 1 f =:C by means of CG ker C = ker B = ker B 1 Convergence if K n (C, r 0 ) remains l 2 orthogonal to ker C for all n x 0 = 0 = r 0 = B A 1 f R(B ) = (ker C) C n r 0 is orthogonal to ker C for all n 1 ( JKU Linz ) Spacetime DPG method January 30, / 41
34 Techniques to solve despite the null space Technique 2: Regularization of the linear system Rewriting B A 1 Bx = B A 1 f in block form yields [ ] B 0 A 1 B 0 B0 A 1 B 1 B1 A 1 B 0 B1 x = B A 1 f A 1 B 1 Solving the invertible system [ B 0 A 1 B 0 B0 A 1 B 1 B1 A 1 B 0 B1 A 1 B 1 + αm ] x = B A 1 f with the mass matrix M jl = (z l, z j ) and α > 0, e.g., α = ( JKU Linz ) Spacetime DPG method January 30, / 41
35 Convergence rates in two-dimensional spacetime Ω = (0, 1) 2, homogeneous boundary and initial conditions, the exact solution of the second order wave equation is given by φ(x, t) = sin(πx)sin 2 (πt) which results in a solution u = [ πcos(πx)sin 2 (πt) πsin(πx)sin(2πt) ] of the first order system g = 0, f = π 2 sin(πx)(2cos(2πt) + sin 2 (πt)) ( JKU Linz ) Spacetime DPG method January 30, / 41
36 Convergence rates in two-dimensional spacetime Figure: L 2 convergence rates for u h, taken from [J. Gopalakrishnan and P. Sepúlveda, (2017)] ( JKU Linz ) Spacetime DPG method January 30, / 41
37 Adaptivity Ω = (0, 1) 2, f = g = 0, homogeneous boundary conditions, i.e., µ = 0 and non-zero initial conditions with µ t=0 = φ 0, q t=0 = φ 0 φ 0 (x) = exp( 1000((x 0.5) 2 )) ( JKU Linz ) Spacetime DPG method January 30, / 41
38 Adaptivity Figure: Numerical pressure µ, adaptive refinement, setting p = 3, taken from [J. Gopalakrishnan and P. Sepúlveda, (2017)] ( JKU Linz ) Spacetime DPG method January 30, / 41
39 Convergence rates in three-dimensional spacetime Ω = (0, 1) 3, homogeneous boundary and initial conditions, the exact solution of the second order wave equation is given by which results in a solution φ(x, t) = sin(πx)sin(πy)t 2 u = of the first order system πcos(πx)sin(πy)t 2 πcos(πy)sin(2πx)t 2 2sin(πx)sin(πy)t g = 0, f = sin(πx)sin(πy)(2 + 2π 2 t 2 ) ( JKU Linz ) Spacetime DPG method January 30, / 41
40 Convergence rates in three-dimensional spacetime Figure: L 2 convergence rates for u h, taken from [J. Gopalakrishnan and P. Sepúlveda, (2017)] ( JKU Linz ) Spacetime DPG method January 30, / 41
41 References [1] L. Demkowicz, J. Gopalakrishnan, S. Nagaraj, and P. Sepúlveda. A spacetime DPG method for the Schrödinger equation. SIAM Journal on Numerical Analysis, 55(4): , [2] A. Ern, J.-L. Guermond, and G. Caplain. An intrinsic criterion for the bijectivity of Hilbert operators related to Friedrichs systems. Communications in Partial Differential Equations, 32: , [3] J. Gopalakrishnan and P. Sepúlveda. A spacetime DPG method for acoustic waves. arxiv: [math. NA], Sept ( JKU Linz ) Spacetime DPG method January 30, / 41
A spacetime DPG method for the wave equation in multiple dimensions
XXXX, 1 24 De Gruyter YYYY A spacetime DPG method for the wave equation in multiple dimensions Jay Gopalakrishnan and Paulina Sepúlveda Abstract. A spacetime discontinuous Petrov-Galerkin (DPG) method
More informationNumerical Solutions to Partial Differential Equations
Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University Nonconformity and the Consistency Error First Strang Lemma Abstract Error Estimate
More informationNumerical Solutions to Partial Differential Equations
Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University The Residual and Error of Finite Element Solutions Mixed BVP of Poisson Equation
More informationINF-SUP CONDITION FOR OPERATOR EQUATIONS
INF-SUP CONDITION FOR OPERATOR EQUATIONS LONG CHEN We study the well-posedness of the operator equation (1) T u = f. where T is a linear and bounded operator between two linear vector spaces. We give equivalent
More informationInequalities of Babuška-Aziz and Friedrichs-Velte for differential forms
Inequalities of Babuška-Aziz and Friedrichs-Velte for differential forms Martin Costabel Abstract. For sufficiently smooth bounded plane domains, the equivalence between the inequalities of Babuška Aziz
More informationNon-stationary Friedrichs systems
Department of Mathematics, University of Osijek BCAM, Bilbao, November 2013 Joint work with Marko Erceg 1 Stationary Friedrichs systems Classical theory Abstract theory 2 3 Motivation Stationary Friedrichs
More informationGoal. Robust A Posteriori Error Estimates for Stabilized Finite Element Discretizations of Non-Stationary Convection-Diffusion Problems.
Robust A Posteriori Error Estimates for Stabilized Finite Element s of Non-Stationary Convection-Diffusion Problems L. Tobiska and R. Verfürth Universität Magdeburg Ruhr-Universität Bochum www.ruhr-uni-bochum.de/num
More informationDiscontinuous Petrov-Galerkin Methods
Discontinuous Petrov-Galerkin Methods Friederike Hellwig 1st CENTRAL School on Analysis and Numerics for Partial Differential Equations, November 12, 2015 Motivation discontinuous Petrov-Galerkin (dpg)
More informationTraces and Duality Lemma
Traces and Duality Lemma Recall the duality lemma with H / ( ) := γ 0 (H ()) defined as the trace space of H () endowed with minimal extension norm; i.e., for w H / ( ) L ( ), w H / ( ) = min{ ŵ H () ŵ
More informationLECTURE 1: SOURCES OF ERRORS MATHEMATICAL TOOLS A PRIORI ERROR ESTIMATES. Sergey Korotov,
LECTURE 1: SOURCES OF ERRORS MATHEMATICAL TOOLS A PRIORI ERROR ESTIMATES Sergey Korotov, Institute of Mathematics Helsinki University of Technology, Finland Academy of Finland 1 Main Problem in Mathematical
More informationFinite-dimensional spaces. C n is the space of n-tuples x = (x 1,..., x n ) of complex numbers. It is a Hilbert space with the inner product
Chapter 4 Hilbert Spaces 4.1 Inner Product Spaces Inner Product Space. A complex vector space E is called an inner product space (or a pre-hilbert space, or a unitary space) if there is a mapping (, )
More information1 The Stokes System. ρ + (ρv) = ρ g(x), and the conservation of momentum has the form. ρ v (λ 1 + µ 1 ) ( v) µ 1 v + p = ρ f(x) in Ω.
1 The Stokes System The motion of a (possibly compressible) homogeneous fluid is described by its density ρ(x, t), pressure p(x, t) and velocity v(x, t). Assume that the fluid is barotropic, i.e., the
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 informationAN ALGEBRAIC APPROACH TO GENERALIZED MEASURES OF INFORMATION
AN ALGEBRAIC APPROACH TO GENERALIZED MEASURES OF INFORMATION Daniel Halpern-Leistner 6/20/08 Abstract. I propose an algebraic framework in which to study measures of information. One immediate consequence
More informationAnalysis of the practical DPG Method
Analysis of the practical DPG Method M. Mandlmayr Supervisor: o. Univ.-Prof. Dipl.-Ing. Dr. Ulrich LANGER January 9, 2018 M. Mandlmayr (JKU) DPG January 9, 2018 1 / 25 Overview 1 Ideal Testfunctions 2
More informationAdaptive methods for control problems with finite-dimensional control space
Adaptive methods for control problems with finite-dimensional control space Saheed Akindeinde and Daniel Wachsmuth Johann Radon Institute for Computational and Applied Mathematics (RICAM) Austrian Academy
More informationPreconditioned space-time boundary element methods for the heat equation
W I S S E N T E C H N I K L E I D E N S C H A F T Preconditioned space-time boundary element methods for the heat equation S. Dohr and O. Steinbach Institut für Numerische Mathematik Space-Time Methods
More informationFunctional Analysis Review
Outline 9.520: Statistical Learning Theory and Applications February 8, 2010 Outline 1 2 3 4 Vector Space Outline A vector space is a set V with binary operations +: V V V and : R V V such that for all
More informationBasic Concepts of Adaptive Finite Element Methods for Elliptic Boundary Value Problems
Basic Concepts of Adaptive Finite lement Methods for lliptic Boundary Value Problems Ronald H.W. Hoppe 1,2 1 Department of Mathematics, University of Houston 2 Institute of Mathematics, University of Augsburg
More informationYour first day at work MATH 806 (Fall 2015)
Your first day at work MATH 806 (Fall 2015) 1. Let X be a set (with no particular algebraic structure). A function d : X X R is called a metric on X (and then X is called a metric space) when d satisfies
More informationChapter 2 Finite Element Spaces for Linear Saddle Point Problems
Chapter 2 Finite Element Spaces for Linear Saddle Point Problems Remark 2.1. Motivation. This chapter deals with the first difficulty inherent to the incompressible Navier Stokes equations, see Remark
More informationSolutions of Selected Problems
1 Solutions of Selected Problems October 16, 2015 Chapter I 1.9 Consider the potential equation in the disk := {(x, y) R 2 ; x 2 +y 2 < 1}, with the boundary condition u(x) = g(x) r for x on the derivative
More informationASYMPTOTICALLY EXACT A POSTERIORI ESTIMATORS FOR THE POINTWISE GRADIENT ERROR ON EACH ELEMENT IN IRREGULAR MESHES. PART II: THE PIECEWISE LINEAR CASE
MATEMATICS OF COMPUTATION Volume 73, Number 246, Pages 517 523 S 0025-5718(0301570-9 Article electronically published on June 17, 2003 ASYMPTOTICALLY EXACT A POSTERIORI ESTIMATORS FOR TE POINTWISE GRADIENT
More informationTD 1: Hilbert Spaces and Applications
Université Paris-Dauphine Functional Analysis and PDEs Master MMD-MA 2017/2018 Generalities TD 1: Hilbert Spaces and Applications Exercise 1 (Generalized Parallelogram law). Let (H,, ) be a Hilbert space.
More informationErrata Applied Analysis
Errata Applied Analysis p. 9: line 2 from the bottom: 2 instead of 2. p. 10: Last sentence should read: The lim sup of a sequence whose terms are bounded from above is finite or, and the lim inf of a sequence
More informationA posteriori error estimation for elliptic problems
A posteriori error estimation for elliptic problems Praveen. C praveen@math.tifrbng.res.in Tata Institute of Fundamental Research Center for Applicable Mathematics Bangalore 560065 http://math.tifrbng.res.in
More informationSome SDEs with distributional drift Part I : General calculus. Flandoli, Franco; Russo, Francesco; Wolf, Jochen
Title Author(s) Some SDEs with distributional drift Part I : General calculus Flandoli, Franco; Russo, Francesco; Wolf, Jochen Citation Osaka Journal of Mathematics. 4() P.493-P.54 Issue Date 3-6 Text
More informationTitle: Localized self-adjointness of Schrödinger-type operators on Riemannian manifolds. Proposed running head: Schrödinger-type operators on
Title: Localized self-adjointness of Schrödinger-type operators on Riemannian manifolds. Proposed running head: Schrödinger-type operators on manifolds. Author: Ognjen Milatovic Department Address: Department
More informationOn the relationship of local projection stabilization to other stabilized methods for one-dimensional advection-diffusion equations
On the relationship of local projection stabilization to other stabilized methods for one-dimensional advection-diffusion equations Lutz Tobiska Institut für Analysis und Numerik Otto-von-Guericke-Universität
More informationWeak Galerkin Finite Element Scheme and Its Applications
Weak Galerkin Finite Element Scheme and Its Applications Ran Zhang Department of Mathematics Jilin University, China IMS, Singapore February 6, 2015 Talk Outline Motivation WG FEMs: Weak Operators + Stabilizer
More informationScientific Computing WS 2017/2018. Lecture 18. Jürgen Fuhrmann Lecture 18 Slide 1
Scientific Computing WS 2017/2018 Lecture 18 Jürgen Fuhrmann juergen.fuhrmann@wias-berlin.de Lecture 18 Slide 1 Lecture 18 Slide 2 Weak formulation of homogeneous Dirichlet problem Search u H0 1 (Ω) (here,
More informationThe Mortar Boundary Element Method
The Mortar Boundary Element Method A Thesis submitted for the degree of Doctor of Philosophy by Martin Healey School of Information Systems, Computing and Mathematics Brunel University March 2010 Abstract
More informationSpace-time Finite Element Methods for Parabolic Evolution Problems
Space-time Finite Element Methods for Parabolic Evolution Problems with Variable Coefficients Ulrich Langer, Martin Neumüller, Andreas Schafelner Johannes Kepler University, Linz Doctoral Program Computational
More informationPDEs in Image Processing, Tutorials
PDEs in Image Processing, Tutorials Markus Grasmair Vienna, Winter Term 2010 2011 Direct Methods Let X be a topological space and R: X R {+ } some functional. following definitions: The mapping R is lower
More informationLaplace s Equation. Chapter Mean Value Formulas
Chapter 1 Laplace s Equation Let be an open set in R n. A function u C 2 () is called harmonic in if it satisfies Laplace s equation n (1.1) u := D ii u = 0 in. i=1 A function u C 2 () is called subharmonic
More informationMixed Finite Element Methods. Douglas N. Arnold, University of Minnesota The 41st Woudschoten Conference 5 October 2016
Mixed Finite Element Methods Douglas N. Arnold, University of Minnesota The 41st Woudschoten Conference 5 October 2016 Linear elasticity Given the load f : Ω R n, find the displacement u : Ω R n and the
More informationContents: 1. Minimization. 2. The theorem of Lions-Stampacchia for variational inequalities. 3. Γ -Convergence. 4. Duality mapping.
Minimization Contents: 1. Minimization. 2. The theorem of Lions-Stampacchia for variational inequalities. 3. Γ -Convergence. 4. Duality mapping. 1 Minimization A Topological Result. Let S be a topological
More informationTHE FORM SUM AND THE FRIEDRICHS EXTENSION OF SCHRÖDINGER-TYPE OPERATORS ON RIEMANNIAN MANIFOLDS
THE FORM SUM AND THE FRIEDRICHS EXTENSION OF SCHRÖDINGER-TYPE OPERATORS ON RIEMANNIAN MANIFOLDS OGNJEN MILATOVIC Abstract. We consider H V = M +V, where (M, g) is a Riemannian manifold (not necessarily
More informationHigh order FEM vs. IgA by J.A.Evans, T.J.Hughes and A.Reali, 2014
High order FEM vs. IgA by J.A.Evans, T.J.Hughes and A.Reali, 2014 Rainer Schneckenleitner JKU Linz rainer.schneckenleitner@gmx.at Jan 20, 2017 Rainer Schneckenleitner (SS16) Isogeometric Analysis Jan 20,
More information2. Dual space is essential for the concept of gradient which, in turn, leads to the variational analysis of Lagrange multipliers.
Chapter 3 Duality in Banach Space Modern optimization theory largely centers around the interplay of a normed vector space and its corresponding dual. The notion of duality is important for the following
More informationTHE INVERSE FUNCTION THEOREM FOR LIPSCHITZ MAPS
THE INVERSE FUNCTION THEOREM FOR LIPSCHITZ MAPS RALPH HOWARD DEPARTMENT OF MATHEMATICS UNIVERSITY OF SOUTH CAROLINA COLUMBIA, S.C. 29208, USA HOWARD@MATH.SC.EDU Abstract. This is an edited version of a
More informationAnalysis Preliminary Exam Workshop: Hilbert Spaces
Analysis Preliminary Exam Workshop: Hilbert Spaces 1. Hilbert spaces A Hilbert space H is a complete real or complex inner product space. Consider complex Hilbert spaces for definiteness. If (, ) : H H
More informationMimetic Finite Difference methods
Mimetic Finite Difference methods An introduction Andrea Cangiani Università di Roma La Sapienza Seminario di Modellistica Differenziale Numerica 2 dicembre 2008 Andrea Cangiani (IAC CNR) mimetic finite
More informationYongdeok Kim and Seki Kim
J. Korean Math. Soc. 39 (00), No. 3, pp. 363 376 STABLE LOW ORDER NONCONFORMING QUADRILATERAL FINITE ELEMENTS FOR THE STOKES PROBLEM Yongdeok Kim and Seki Kim Abstract. Stability result is obtained for
More informationLocal pointwise a posteriori gradient error bounds for the Stokes equations. Stig Larsson. Heraklion, September 19, 2011 Joint work with A.
Local pointwise a posteriori gradient error bounds for the Stokes equations Stig Larsson Department of Mathematical Sciences Chalmers University of Technology and University of Gothenburg Heraklion, September
More informationTraces, extensions and co-normal derivatives for elliptic systems on Lipschitz domains
Traces, extensions and co-normal derivatives for elliptic systems on Lipschitz domains Sergey E. Mikhailov Brunel University West London, Department of Mathematics, Uxbridge, UB8 3PH, UK J. Math. Analysis
More informationFinite Element Clifford Algebra: A New Toolkit for Evolution Problems
Finite Element Clifford Algebra: A New Toolkit for Evolution Problems Andrew Gillette joint work with Michael Holst Department of Mathematics University of California, San Diego http://ccom.ucsd.edu/ agillette/
More informationICES REPORT The DPG Method for the Stokes Problem.
ICES REPORT 2-22 June 202 The DPG Method for the Stokes Problem. by N. Roberts, T. Bui-Thanh, and L. Demkowicz The Institute for Computational Engineering and Sciences The University of Texas at Austin
More informationA class of non-convex polytopes that admit no orthonormal basis of exponentials
A class of non-convex polytopes that admit no orthonormal basis of exponentials Mihail N. Kolountzakis and Michael Papadimitrakis 1 Abstract A conjecture of Fuglede states that a bounded measurable set
More informationUnbounded operators on Hilbert spaces
Chapter 1 Unbounded operators on Hilbert spaces Definition 1.1. Let H 1, H 2 be Hilbert spaces and T : dom(t ) H 2 be a densely defined linear operator, i.e. dom(t ) is a dense linear subspace of H 1.
More information1 Constrained Optimization
1 Constrained Optimization Let B be an M N matrix, a linear operator from the space IR N to IR M with adjoint B T. For (column vectors) x IR N and y IR M we have x B T y = Bx y. This vanishes for all y
More informationA Finite Element Method Using Singular Functions for Poisson Equations: Mixed Boundary Conditions
A Finite Element Method Using Singular Functions for Poisson Equations: Mixed Boundary Conditions Zhiqiang Cai Seokchan Kim Sangdong Kim Sooryun Kong Abstract In [7], we proposed a new finite element method
More informationNormality of adjointable module maps
MATHEMATICAL COMMUNICATIONS 187 Math. Commun. 17(2012), 187 193 Normality of adjointable module maps Kamran Sharifi 1, 1 Department of Mathematics, Shahrood University of Technology, P. O. Box 3619995161-316,
More informationFINITE ELEMENT APPROXIMATION OF STOKES-LIKE SYSTEMS WITH IMPLICIT CONSTITUTIVE RELATION
Proceedings of ALGORITMY pp. 9 3 FINITE ELEMENT APPROXIMATION OF STOKES-LIKE SYSTEMS WITH IMPLICIT CONSTITUTIVE RELATION JAN STEBEL Abstract. The paper deals with the numerical simulations of steady flows
More informationA brief introduction to finite element methods
CHAPTER A brief introduction to finite element methods 1. Two-point boundary value problem and the variational formulation 1.1. The model problem. Consider the two-point boundary value problem: Given a
More information2.3 Variational form of boundary value problems
2.3. VARIATIONAL FORM OF BOUNDARY VALUE PROBLEMS 21 2.3 Variational form of boundary value problems Let X be a separable Hilbert space with an inner product (, ) and norm. We identify X with its dual X.
More informationRecitation 1 (Sep. 15, 2017)
Lecture 1 8.321 Quantum Theory I, Fall 2017 1 Recitation 1 (Sep. 15, 2017) 1.1 Simultaneous Diagonalization In the last lecture, we discussed the situations in which two operators can be simultaneously
More informationFunctional Analysis Review
Functional Analysis Review Lorenzo Rosasco slides courtesy of Andre Wibisono 9.520: Statistical Learning Theory and Applications September 9, 2013 1 2 3 4 Vector Space A vector space is a set V with binary
More informationNumerical approximation of output functionals for Maxwell equations
Numerical approximation of output functionals for Maxwell equations Ferenc Izsák ELTE, Budapest University of Twente, Enschede 11 September 2004 MAXWELL EQUATIONS Assumption: electric field ( electromagnetic
More information1 Math 241A-B Homework Problem List for F2015 and W2016
1 Math 241A-B Homework Problem List for F2015 W2016 1.1 Homework 1. Due Wednesday, October 7, 2015 Notation 1.1 Let U be any set, g be a positive function on U, Y be a normed space. For any f : U Y let
More informationI teach myself... Hilbert spaces
I teach myself... Hilbert spaces by F.J.Sayas, for MATH 806 November 4, 2015 This document will be growing with the semester. Every in red is for you to justify. Even if we start with the basic definition
More informationMultiscale methods for time-harmonic acoustic and elastic wave propagation
Multiscale methods for time-harmonic acoustic and elastic wave propagation Dietmar Gallistl (joint work with D. Brown and D. Peterseim) Institut für Angewandte und Numerische Mathematik Karlsruher Institut
More informationNon-degeneracy of perturbed solutions of semilinear partial differential equations
Non-degeneracy of perturbed solutions of semilinear partial differential equations Robert Magnus, Olivier Moschetta Abstract The equation u + F(V (εx, u = 0 is considered in R n. For small ε > 0 it is
More informationCONVERGENCE THEORY. G. ALLAIRE CMAP, Ecole Polytechnique. 1. Maximum principle. 2. Oscillating test function. 3. Two-scale convergence
1 CONVERGENCE THEOR G. ALLAIRE CMAP, Ecole Polytechnique 1. Maximum principle 2. Oscillating test function 3. Two-scale convergence 4. Application to homogenization 5. General theory H-convergence) 6.
More informationA Posteriori Error Estimation Techniques for Finite Element Methods. Zhiqiang Cai Purdue University
A Posteriori Error Estimation Techniques for Finite Element Methods Zhiqiang Cai Purdue University Department of Mathematics, Purdue University Slide 1, March 16, 2017 Books Ainsworth & Oden, A posteriori
More informationWhere is matrix multiplication locally open?
Linear Algebra and its Applications 517 (2017) 167 176 Contents lists available at ScienceDirect Linear Algebra and its Applications www.elsevier.com/locate/laa Where is matrix multiplication locally open?
More informationNonlinear stabilization via a linear observability
via a linear observability Kaïs Ammari Department of Mathematics University of Monastir Joint work with Fathia Alabau-Boussouira Collocated feedback stabilization Outline 1 Introduction and main result
More informationL. Levaggi A. Tabacco WAVELETS ON THE INTERVAL AND RELATED TOPICS
Rend. Sem. Mat. Univ. Pol. Torino Vol. 57, 1999) L. Levaggi A. Tabacco WAVELETS ON THE INTERVAL AND RELATED TOPICS Abstract. We use an abstract framework to obtain a multilevel decomposition of a variety
More informationA variation on the infsup condition
Cinquième foire européenne aux éléments finis Centre International de Rencontres Mathématiques Marseille Luminy, 17-19 mai 2007 A variation on the infsup condition François Dubois CNAM Paris and University
More informationSpectral Continuity Properties of Graph Laplacians
Spectral Continuity Properties of Graph Laplacians David Jekel May 24, 2017 Overview Spectral invariants of the graph Laplacian depend continuously on the graph. We consider triples (G, x, T ), where G
More informationNumerische Mathematik
Numer. Math. (2012) 122:61 99 DOI 10.1007/s00211-012-0456-x Numerische Mathematik C 0 elements for generalized indefinite Maxwell equations Huoyuan Duan Ping Lin Roger C. E. Tan Received: 31 July 2010
More informationENERGY NORM A POSTERIORI ERROR ESTIMATES FOR MIXED FINITE ELEMENT METHODS
ENERGY NORM A POSTERIORI ERROR ESTIMATES FOR MIXED FINITE ELEMENT METHODS CARLO LOVADINA AND ROLF STENBERG Abstract The paper deals with the a-posteriori error analysis of mixed finite element methods
More information1 Continuity Classes C m (Ω)
0.1 Norms 0.1 Norms A norm on a linear space X is a function : X R with the properties: Positive Definite: x 0 x X (nonnegative) x = 0 x = 0 (strictly positive) λx = λ x x X, λ C(homogeneous) x + y x +
More informationarxiv: v1 [math.na] 29 Feb 2016
EFFECTIVE IMPLEMENTATION OF THE WEAK GALERKIN FINITE ELEMENT METHODS FOR THE BIHARMONIC EQUATION LIN MU, JUNPING WANG, AND XIU YE Abstract. arxiv:1602.08817v1 [math.na] 29 Feb 2016 The weak Galerkin (WG)
More information1. Introduction Since the pioneering work by Leray [3] in 1934, there have been several studies on solutions of Navier-Stokes equations
Math. Res. Lett. 13 (6, no. 3, 455 461 c International Press 6 NAVIER-STOKES EQUATIONS IN ARBITRARY DOMAINS : THE FUJITA-KATO SCHEME Sylvie Monniaux Abstract. Navier-Stokes equations are investigated in
More informationSplitting methods with boundary corrections
Splitting methods with boundary corrections Alexander Ostermann University of Innsbruck, Austria Joint work with Lukas Einkemmer Verona, April/May 2017 Strang s paper, SIAM J. Numer. Anal., 1968 S (5)
More informationMULTIGRID 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 informationDomain Decomposition Algorithms for an Indefinite Hypersingular Integral Equation in Three Dimensions
Domain Decomposition Algorithms for an Indefinite Hypersingular Integral Equation in Three Dimensions Ernst P. Stephan 1, Matthias Maischak 2, and Thanh Tran 3 1 Institut für Angewandte Mathematik, Leibniz
More informationSpectral element agglomerate AMGe
Spectral element agglomerate AMGe T. Chartier 1, R. Falgout 2, V. E. Henson 2, J. E. Jones 4, T. A. Manteuffel 3, S. F. McCormick 3, J. W. Ruge 3, and P. S. Vassilevski 2 1 Department of Mathematics, Davidson
More informationOn the distributional divergence of vector fields vanishing at infinity
Proceedings of the Royal Society of Edinburgh, 141A, 65 76, 2011 On the distributional divergence of vector fields vanishing at infinity Thierry De Pauw Institut de Recherches en Mathématiques et Physique,
More informationMultilevel Preconditioning of Graph-Laplacians: Polynomial Approximation of the Pivot Blocks Inverses
Multilevel Preconditioning of Graph-Laplacians: Polynomial Approximation of the Pivot Blocks Inverses P. Boyanova 1, I. Georgiev 34, S. Margenov, L. Zikatanov 5 1 Uppsala University, Box 337, 751 05 Uppsala,
More informationb i (x) u + c(x)u = f in Ω,
SIAM J. NUMER. ANAL. Vol. 39, No. 6, pp. 1938 1953 c 2002 Society for Industrial and Applied Mathematics SUBOPTIMAL AND OPTIMAL CONVERGENCE IN MIXED FINITE ELEMENT METHODS ALAN DEMLOW Abstract. An elliptic
More informationSPECTRAL PROPERTIES OF THE LAPLACIAN ON BOUNDED DOMAINS
SPECTRAL PROPERTIES OF THE LAPLACIAN ON BOUNDED DOMAINS TSOGTGEREL GANTUMUR Abstract. After establishing discrete spectra for a large class of elliptic operators, we present some fundamental spectral properties
More informationWELL-POSEDNESS FOR HYPERBOLIC PROBLEMS (0.2)
WELL-POSEDNESS FOR HYPERBOLIC PROBLEMS We will use the familiar Hilbert spaces H = L 2 (Ω) and V = H 1 (Ω). We consider the Cauchy problem (.1) c u = ( 2 t c )u = f L 2 ((, T ) Ω) on [, T ] Ω u() = u H
More informationMinimal periods of semilinear evolution equations with Lipschitz nonlinearity
Minimal periods of semilinear evolution equations with Lipschitz nonlinearity James C. Robinson a Alejandro Vidal-López b a Mathematics Institute, University of Warwick, Coventry, CV4 7AL, U.K. b Departamento
More informationConservative Control Systems Described by the Schrödinger Equation
Conservative Control Systems Described by the Schrödinger Equation Salah E. Rebiai Abstract An important subclass of well-posed linear systems is formed by the conservative systems. A conservative system
More informationExercises - Chapter 1 - Chapter 2 (Correction)
Université de Nice Sophia-Antipolis Master MathMods - Finite Elements - 28/29 Exercises - Chapter 1 - Chapter 2 Correction) Exercise 1. a) Let I =], l[, l R. Show that Cl) >, u C Ī) Cl) u H 1 I), u DĪ).
More informationControllability of the linear 1D wave equation with inner moving for
Controllability of the linear D wave equation with inner moving forces ARNAUD MÜNCH Université Blaise Pascal - Clermont-Ferrand - France Toulouse, May 7, 4 joint work with CARLOS CASTRO (Madrid) and NICOLAE
More informationFriedrich symmetric systems
viii CHAPTER 8 Friedrich symmetric systems In this chapter, we describe a theory due to Friedrich [13] for positive symmetric systems, which gives the existence and uniqueness of weak solutions of boundary
More informationKernel Method: Data Analysis with Positive Definite Kernels
Kernel Method: Data Analysis with Positive Definite Kernels 2. Positive Definite Kernel and Reproducing Kernel Hilbert Space Kenji Fukumizu The Institute of Statistical Mathematics. Graduate University
More informationMAT 578 FUNCTIONAL ANALYSIS EXERCISES
MAT 578 FUNCTIONAL ANALYSIS EXERCISES JOHN QUIGG Exercise 1. Prove that if A is bounded in a topological vector space, then for every neighborhood V of 0 there exists c > 0 such that tv A for all t > c.
More informationICES REPORT Global properties of DPG test spaces for convection-diffusion problems
ICES REPORT 13-05 February 2013 Global properties of DPG test spaces for convection-diffusion problems by Jesse Chan, Jay Gopalakrishnan, and Leszek Demkowicz The Institute for Computational Engineering
More informationICES REPORT Robust DPG Methods for Transient Convection-Diffusion
ICES REPORT 5- October 05 Robust DPG Methods for Transient Convection-Diffusion by Truman Ellis, Jesse Chan, and Leszek Demkowicz The Institute for Computational Engineering and Sciences The University
More informationOn invariant graph subspaces
On invariant graph subspaces Konstantin A. Makarov, Stephan Schmitz and Albrecht Seelmann Abstract. In this paper we discuss the problem of decomposition for unbounded 2 2 operator matrices by a pair of
More informationMicrolocal analysis and inverse problems Lecture 3 : Carleman estimates
Microlocal analysis and inverse problems ecture 3 : Carleman estimates David Dos Santos Ferreira AGA Université de Paris 13 Monday May 16 Instituto de Ciencias Matemáticas, Madrid David Dos Santos Ferreira
More informationMaxwell s equations in Carnot groups
Maxwell s equations in Carnot groups B. Franchi (U. Bologna) INDAM Meeting on Geometric Control and sub-riemannian Geometry Cortona, May 21-25, 2012 in honor of Andrey Agrachev s 60th birthday Researches
More informationA new approach for Kirchhoff-Love plates and shells
A new approach for Kirchhoff-Love plates and shells Walter Zulehner Institute of Computational Mathematics JKU Linz, Austria AANMPDE 10 October 2-6, 2017, Paleochora, Crete, Greece Walter Zulehner (JKU
More informationThe duality between the gradient and divergence operators on bounded Lipschitz domains
The duality between the gradient and divergence operators on bounded Lipschitz domains Mikael Kurula and Hans Zwart Department of Applied Mathematics The University of Twente October 15, 2012 Abstract
More informationApplied Analysis (APPM 5440): Final exam 1:30pm 4:00pm, Dec. 14, Closed books.
Applied Analysis APPM 44: Final exam 1:3pm 4:pm, Dec. 14, 29. Closed books. Problem 1: 2p Set I = [, 1]. Prove that there is a continuous function u on I such that 1 ux 1 x sin ut 2 dt = cosx, x I. Define
More informationSecond Order Elliptic PDE
Second Order Elliptic PDE T. Muthukumar tmk@iitk.ac.in December 16, 2014 Contents 1 A Quick Introduction to PDE 1 2 Classification of Second Order PDE 3 3 Linear Second Order Elliptic Operators 4 4 Periodic
More information