Sparse Tensor Galerkin Discretizations for First Order Transport Problems
|
|
- Octavia Gilbert
- 5 years ago
- Views:
Transcription
1 Sparse Tensor Galerkin Discretizations for First Order Transport Problems Ch. Schwab R. Hiptmair, E. Fonn, K. Grella, G. Widmer ETH Zürich, Seminar for Applied Mathematics IMA WS Novel Discretization Methods Oct 31-Nov Acknowledgement: Swiss Nat l Science Foundation and ERC W. Dahmen and G. Kutyniok Priority Research Programme 1324 German Research Foundation (DFG).
2 Ch. Schwab (SAM, ETH Zurich) Sparse Tensor Approximation Novel Discretization Methods 2 / 38 Outline 1 Problem 2 Variational Formulations 3 Galerkin discretization Without boundary conditions With boundary conditions 4 Numerical experiments 5 Tensor Product Wavelet Galerkin Approach (Widmer etal. JCP 2008) Inflow Boundary: Example Preconditioning Subspace Correction Preconditioner Experiment I 6 Adaptive Approximation Adaptivity: Rationale Numerical Experiment II Numerical Experiment III 7 Conclusion
3 Ch. Schwab (SAM, ETH Zurich) Sparse Tensor Approximation Novel Discretization Methods 3 / 38 The Radiative Transfer Problem Stationary monochromatic radiative transfer problem: (s x + κ(x)) u(x, s) = κ(x)i b (x), (x, s) D S d S u(x, s) = 0, x Γ (s) := D {s : s n(x) < 0}. x D R d D bounded, Lipschitz, s S d S 5-dimensional phase space (d D, d S ) = (3, 2) D R d D bounded: zero inflow BCs Model for general, parametric transport problems (UQ for Neutron Transport, v x ψ(x, v; t)...) Numerical approaches: MCM, discrete ordinates, P N approximation Complexity h d D D h d S S, slow convergence due to curse of dimensionality Sparse tensorization: Wavelets in S d S ( Widmer et al. (2008) ) Current: Sparse Adaptive Tensor Approximation with 1 spectral discretization in S d S ( this talk) 2 hierarchic, directional representation systems in D ( W. Dahmen s talk)
4 Variational formulation cg unstable = a) Stabilization or b) Petrov Galerkin (Talk of W. Dahmen) Stabilized LSQ formulation (Manteuffel, Ressel, Starke: SINUM(2000)): Find u V 0 such that with a(u, v) = l(v) v V 0 a(u, v) := (ɛ(s x u + κu), s x v + κv) L2 (D S), ɛ(x) = 1 κ(x) l(v) := (ɛκi b, s x v + κv) L2 (D S d S ) V := { u L 2 (D S d S ) : s x u L 2 (D S d S ) } V 0 := {u V : u = 0 on D, s n(x) < 0} a(, ) is continuous and coercive on V 0 V 0 Energy Norm u A := a(u, u) Ch. Schwab (SAM, ETH Zurich) Sparse Tensor Approximation Novel Discretization Methods 4 / 38
5 Ch. Schwab (SAM, ETH Zurich) Sparse Tensor Approximation Novel Discretization Methods 5 / 38 Discretization Cartesian product domain D S d S span V with product basis functions from tensor product (TP) space H 1 (D) L 2 (S d S ) Inflow Boundary s-dependent: TP structure broken Hierarchical function spaces allow for incremental resolution refinement: On D: Piecewise polynomials (here p = 1) on dyadically refined mesh TD l V l D := S p,1 (D, T l D) On S d S : Spherical harmonics S (d S) n,m V l D = V l 1 D W l D up to order l H 1 (D) V l S := P d S l = span{s (d S) n,m : n = 0,..., l; m = 1,..., m n,ds } L 2 (S d S ) V l S = V l 1 S W l S Finite dimensional Full Tensor Product (FTP) space: V L = VD L VS N = 0 l D L 0 l S N W l D D W l S S H 1 (D) L 2 (S d S )
6 Ch. Schwab (SAM, ETH Zurich) Sparse Tensor Approximation Novel Discretization Methods 6 / 38 Structure of FTP space Illustration with VD L=3 = S (1,1) ([0, 1] [0, 1], TD 3 N=3 ) and VS = P 1 3 W l S S W 0 D W 0 S W 0 D W 1 S W 0 D W 2 S W 0 D W 3 S W 1 D W 0 S W 1 D W 1 S W 1 D W 2 S W 1 D W 3 S W 2 D W 0 S W 2 D W 1 S W 2 D W 2 S W 2 D W 3 S W 3 D W 0 S W 3 D W 1 S W 3 D W 2 S W 3 D W 3 S W l D D
7 Ch. Schwab (SAM, ETH Zurich) Sparse Tensor Approximation Novel Discretization Methods 7 / 38 Structure of Sparse Tensor Product (STP) space Illustration with VD L=3 = S (1,1) ([0, 1] [0, 1], TD 3 N=3 ) and VS = P 1 3 W l S S W 0 D W 0 S W 0 D W 1 S W 0 D W 2 S W 0 D W 3 S W 1 D W 0 S W 1 D W 1 S W 1 D W 2 S W 1 D W 3 S W 2 D W 0 S W 2 D W 1 S W 2 D W 2 S W 2 D W 3 S W 3 D W 0 S W 3 D W 1 S W 3 D W 2 S W 3 D W 3 S W l D D
8 Ch. Schwab (SAM, ETH Zurich) Sparse Tensor Approximation Novel Discretization Methods 8 / 38 Sparse Tensor Product (STP) Space FTP Finite Element FTP space: Finite Element STP space: V L = ˆV L := 0 l D L 0 l S N 0 f (l D,l S ) L with cutoff function f : [0, L] [0, N] R. Dimensionality: W l D D W l D D W l S S W l S S, dim(v L ) N D N S dim( ˆV L ) N D log N S + N S log N D for f (l D, l S ) = l D + l S L N
9 Ch. Schwab (SAM, ETH Zurich) Sparse Tensor Approximation Novel Discretization Methods 9 / 38 Approximation properties Error estimates for approximation u L of u H s+1,t (D S), s [0, p], t N 0 : Full tensor case u L V L : u u L H 1,0 (D S) max{n t/d S S, N s/d D } u (H 1,t H 1+s,0 )(D S) Sparse tensor case u L ˆV L with f (l D, l S ) = l D + L log 2 (l S +1) log 2 (N+1) : u u L H 1,0 (D S) log 2 ( N 1/d D ) max{n t/d S S, N s/d D } u H 1+s,t (D S) Convergence rate (almost) retained while number of dofs reduced. To equilibrate contributions: N h s/t
10 Product basis functions satisfying boundary conditions Zero inflow conditions: u(x, s) = 0 for x D, s n(x) < 0 Physical basis functions and corresponding outflow regions: x 0 0 x x x x π π π 2 s s s π 3π 2 π 3π 2 π x 1 0 3π 2 Partition of sphere into angular regions S q, mesh T S = {S q } n S q=1 : S 2 π 2 S 1 π 0 S 3 S 4 3π 2 Ch. Schwab (SAM, ETH Zurich) Sparse Tensor Approximation Novel Discretization Methods 10 / 38
11 Ch. Schwab (SAM, ETH Zurich) Sparse Tensor Approximation Novel Discretization Methods 11 / 38 Tensor product space with boundary conditions To maintain TP structure and strongly satisfy b.c.s add restricted angular functions to VS L: Introduce space P d S,T S N of angular functions restricted to elements of mesh T S, e. g. (d S = 1) Compressed spherical harmonics Legendre polynomials on S q. Split physical spaces into boundary and non-boundary subspaces: V l D = V l D V l D,0 Tensorize V 0 -conforming combinations only: V0 L := VD,0 L P d S N V D L P d S,T S N Sparse version: ˆV L 0 := V L D,0 ˆ P d S N V L D ˆ P d S,T S N
12 Ch. Schwab (SAM, ETH Zurich) Sparse Tensor Approximation Novel Discretization Methods 12 / 38 Structure of TP spaces with boundary conditions Active dofs for L = 3, N = 8 for FTP and STP space (here d D + d S = 2 + 1): Physical dofs FTP l D S S 1 S 2 S 3 S Angular dofs Physical dofs STP Angular dofs 1853 dofs 275 dofs
13 Ch. Schwab (SAM, ETH Zurich) Sparse Tensor Approximation Novel Discretization Methods 13 / 38 Solution of discretized variational problem Benefits of TP structure: Separable quadrature for stiffness matrices over D and S System of equations is matrix equation: B n UC n = L Matrix-vector multiplication faster CG to solve system (currently no preconditioning). To evaluate solution: Look at Relative error of u in X = H 1,0, A, L 2 -norm on D S: err X = u u L X / u X n
14 Ch. Schwab (SAM, ETH Zurich) Sparse Tensor Approximation Novel Discretization Methods 14 / 38 Numerical experiments D = [0, 1] [0, 1], S = S 1, κ(x) = 1. Experiment 1 Solution prescribed: u(x, ϕ) = U(x)Y (ϕ) Experiment 2 U(x) = ( 1(x 1 1) 2 + 1)( 4(x 2 1/2) 2 + 1)x 1 x 2 ( 4(x 2 1)) 2 π ϕ + 1 if 0 ϕ < π 2 2 Y (ϕ) = π (ϕ 3π 2 ) if 3π 2 < ϕ < 2π 0 else U(x) as in Exp. 1 ( ) 2 sin(8ϕ) Y (ϕ) = 8ϕ if 0 ϕ < π 2 3π 2 < ϕ < 2π 0 else
15 Ch. Schwab (SAM, ETH Zurich) Sparse Tensor Approximation Novel Discretization Methods 15 / 38 Convergence rates Experiment Error u (testside2), FT/ST, N=32, L refined 10 0 Error u (testside2), FT/ST, L=5, N refined u u N / u H10 error FT A error FT L2 error FT H10 error ST A error ST L2 error ST # Dofs p=2 p=1/2 p=1 u u N / u In this setup: error limited by physical contribution p= # Dofs In physical space convergence rate 1/2 for H 1,0 -error as expected In angular space resolution refinement not effective ST method achieves better error per employed dofs ratio p=1/2 p=1
16 Ch. Schwab (SAM, ETH Zurich) Sparse Tensor Approximation Novel Discretization Methods 16 / 38 Convergence rates Experiment Error u (testside1), FT/ST, N=32, L refined 10 0 Error u (testside1), FT/ST, L=5, N refined u u N / u H10 error FT A error FT L2 error FT H10 error ST A error ST L2 error ST # Dofs p=2 p=1/2 p=1 u u N / u Here error limited by angular contribution Physical refinement not effective H10 error FT A error FT L2 error FT H10 error ST A error ST L2 error ST p= Angular convergence rate between 1/2 and 1 for H 1,0 -error # Dofs p=1/2 ST LSQ FEM achieves better error per dofs ratio for angular refinement p=1
17 Ch. Schwab (SAM, ETH Zurich) Sparse Tensor Approximation Novel Discretization Methods 17 / 38 Tensor Product Wavelet Galerkin D R 2 disc, dimd = 2, d S = 3 Meshes T 3 D, T 3 S 2 :
18 Ch. Schwab (SAM, ETH Zurich) Sparse Tensor Approximation Novel Discretization Methods 18 / 38 Multiscale Decompositions Nested meshes T l D nested spaces V l D := S0 1 (T l D ) TS l nested spaces VS l := S 1 0 (TS l) Complement ( detail, hierarchical surplus ) Spaces V l D = W l D V l 1 D V l S = W l S V l 1 S
19 Ch. Schwab (SAM, ETH Zurich) Sparse Tensor Approximation Novel Discretization Methods 19 / 38 Multiscale Decompositions In D R d : + In S d 1 :
20 Ch. Schwab (SAM, ETH Zurich) Sparse Tensor Approximation Novel Discretization Methods 20 / 38 Multiscale Decompositions Full Tensor Galerkin trial/test space V L = VD L VS L = WD i W j S 0 i,j L Sparse tensor product space: V L := WD i W j S dim V L = O(ND L log NS L + NS L log ND) L 0 i+j L
21 Ch. Schwab (SAM, ETH Zurich) Sparse Tensor Approximation Novel Discretization Methods 21 / 38 Sparse Tensor Product Space V L := WD i W j S 0 i+j L
22 Ch. Schwab (SAM, ETH Zurich) Sparse Tensor Approximation Novel Discretization Methods 22 / 38 Graded ST Space Problem: Homogeneous Essential BC s in Γ (s) destroy TP structure Solution: 1 Enforce Dirichlet BC s weakly via Multipliers 2 Enforce Dirichlet BC s strongly: local isotropic mesh refinement along Γ (s), s S 1. Graded ST space
23 Ch. Schwab (SAM, ETH Zurich) Sparse Tensor Approximation Novel Discretization Methods 23 / 38 φ Inflow Boundary: Example y 4 3 x 2 γ(x) φ D Tangent in x to D y x
24 Ch. Schwab (SAM, ETH Zurich) Sparse Tensor Approximation Novel Discretization Methods 24 / 38 Graded Sparse Tensor Product Space D S1 gamma Augment V L 0 := V 0 V L by tensor product nodal basis functions (in neighborhood of γ) multiscale tensor product is dictionary, frame (no Riesz basis anymore) dd Moderate increase of #d.o.f. by a logarithmic factor L. dim V 0 L = O(L 2 ( M L ) O(LML ) ) Approximation rate (proved for n = 2, d = 1, 2): C > 0 : L : û L V L 0 : u û L S C log h L h L u H 2,1 (D S 2 ).
25 Ch. Schwab (SAM, ETH Zurich) Sparse Tensor Approximation Novel Discretization Methods 25 / 38 Preconditioning Multiplicative subspace correction preconditioner
26 Ch. Schwab (SAM, ETH Zurich) Sparse Tensor Approximation Novel Discretization Methods 25 / 38 Preconditioning Multiplicative subspace correction preconditioner transformation in standard basis
27 Ch. Schwab (SAM, ETH Zurich) Sparse Tensor Approximation Novel Discretization Methods 25 / 38 Preconditioning Multiplicative subspace correction preconditioner transformation in standard basis solution of decoupled, low dimensional transport problems (AMG)
28 Ch. Schwab (SAM, ETH Zurich) Sparse Tensor Approximation Novel Discretization Methods 25 / 38 Preconditioning Multiplicative subspace correction preconditioner transformation in standard basis solution of decoupled, low dimensional transport problems (AMG) Effort d.o.f. (Overhead of AMG ignored)
29 Ch. Schwab (SAM, ETH Zurich) Sparse Tensor Approximation Novel Discretization Methods 26 / 38 Preconditioning: Numerical Experiment Regular ST space Absorption coefficient: κ(x) = 1 Blackbody radiation: I b (x) = e 4 x 2
30 Ch. Schwab (SAM, ETH Zurich) Sparse Tensor Approximation Novel Discretization Methods 27 / 38 Numerical Experiment I Throughout: 2.5 D setting, D = {x R 2 x < 1}, n = d = 2 smooth solution u(x, s) due to large (constant) absorption 10 Emission f(r) distance to origin r Blackbody intensity f ( x ) Heat flux q(x)
31 Numerical Experiment I Efficacy of ST a priori DOF selection criterion: 1.18e+00 1 Level in Solid Angle e e 03 Size of multiscale contributions Ch. Schwab (SAM, ETH Zurich) 1.13e 01 l1l 2 l 3 Level in Physical Space l 1.00e largest multiscale contributions of solution in full tensor product space Sparse Tensor Approximation Novel Discretization Methods 28 / 38
32 Ch. Schwab (SAM, ETH Zurich) Sparse Tensor Approximation Novel Discretization Methods 29 / 38 Numerical Experiment I Numerical solutions: Incident Radiation G(x 1,0) Level 0 Full Level 1 Full Level 1 Sparse Level 2 Full Level 2 Sparse Level 3 Full Level 3 Sparse x 1 Incident radiation Net emission (div q)
33 Ch. Schwab (SAM, ETH Zurich) Sparse Tensor Approximation Novel Discretization Methods 30 / 38 Numerical Experiment I Convergence (T l D, T l S by uniform regular refinement): Relative error of the incident radiation in the L 2 (D)-norm Relative error of the heat flux in the L 2 (D)-norm
34 Ch. Schwab (SAM, ETH Zurich) Sparse Tensor Approximation Novel Discretization Methods 31 / 38 Adaptivity: Rationale Propagation of singularities No elliptic lifting theorems for radiative transfer BVP!
35 Ch. Schwab (SAM, ETH Zurich) Sparse Tensor Approximation Novel Discretization Methods 31 / 38 Adaptivity: Rationale Propagation of singularities No elliptic lifting theorems for radiative transfer BVP! Example (dimd = 2, d S emitting region = 1): small compact
36 Ch. Schwab (SAM, ETH Zurich) Sparse Tensor Approximation Novel Discretization Methods 31 / 38 Adaptivity: Rationale Propagation of singularities No elliptic lifting theorems for radiative transfer BVP! Example (dimd = 2, d S emitting region = 1): small compact Note: more smoothness with (isotropic) scattering
37 Ch. Schwab (SAM, ETH Zurich) Sparse Tensor Approximation Novel Discretization Methods 32 / 38 Numerical Experiment II highly non-uniform absorption/source u non-smooth flux Hea
38 Ch. Schwab (SAM, ETH Zurich) Sparse Tensor Approximation Novel Discretization Methods 33 / 38 Numerical Experiment II (Limited) efficacy of a priori selection criterion: Size of multiscale contributions largest multiscale contributions of solution in full tensor product space
39 Ch. Schwab (SAM, ETH Zurich) Sparse Tensor Approximation Novel Discretization Methods 34 / 38 Numerical Experiment II Numerical solutions: Incident radiation Net emission
40 Ch. Schwab (SAM, ETH Zurich) Sparse Tensor Approximation Novel Discretization Methods 35 / 38 Numerical Experiment II Convergence (T l D, T l S by uniform regular refinement): Relative error of the incident radiation in the L 2 (D)-norm Relative error of the heat flux in the L 2 (D)-norm
41 Ch. Schwab (SAM, ETH Zurich) Sparse Tensor Approximation Novel Discretization Methods 36 / 38 Numerical Experiment III { 5, x 0.08 Absorption: κ(x, y) = 3, otherwise. { e 20 x x 0 2 e Source: I b (x) = 5, x x , otherwise. x 0 = ( ) Threshold ν = 0.04
42 Ch. Schwab (SAM, ETH Zurich) Sparse Tensor Approximation Novel Discretization Methods 37 / 38 Conclusion STP combined with spectral and MR discretization of S d S essential b.c.s on Γ (s) strongly satisfied & TP structure maintained STP method results in discretization with far fewer dofs than FTP Theoretical error estimates without b.c.s shown to be (almost) identical for FTP and STP Error estimates with b.c.s numerically confirm theoretical estimates Adaptive mesh refinements in D and in S d S Essential BC s and Directional Adaptivity: stable frame discretizations of transport problems in bounded domains D
43 Ch. Schwab (SAM, ETH Zurich) Sparse Tensor Approximation Novel Discretization Methods 38 / 38 References Konstantin Grella, Christoph Schwab. Sparse tensor spherical harmonics approximation in radiative transfer. Technical report, SAM, ETH Zürich, reports G. Widmer, R. Hiptmair, and Ch. Schwab. Sparse adaptive finite elements for radiative transfer. Journal of Computational Physics, 227: , Gisela Widmer. Sparse Finite Elements for Radiative Transfer. PhD thesis, ETH Zürich, No
Space-time sparse discretization of linear parabolic equations
Space-time sparse discretization of linear parabolic equations Roman Andreev August 2, 200 Seminar for Applied Mathematics, ETH Zürich, Switzerland Support by SNF Grant No. PDFMP2-27034/ Part of PhD thesis
More informationImplementation of Sparse Wavelet-Galerkin FEM for Stochastic PDEs
Implementation of Sparse Wavelet-Galerkin FEM for Stochastic PDEs Roman Andreev ETH ZÜRICH / 29 JAN 29 TOC of the Talk Motivation & Set-Up Model Problem Stochastic Galerkin FEM Conclusions & Outlook Motivation
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 informationTrefftz-discontinuous Galerkin methods for time-harmonic wave problems
Trefftz-discontinuous Galerkin methods for time-harmonic wave problems Ilaria Perugia Dipartimento di Matematica - Università di Pavia (Italy) http://www-dimat.unipv.it/perugia Joint work with Ralf Hiptmair,
More informationTrefftz-Discontinuous Galerkin Methods for Acoustic Scattering - Recent Advances
Trefftz-Discontinuous Galerkin Methods for Acoustic Scattering - Recent Advances Ilaria Perugia Dipartimento di Matematica Università di Pavia (Italy) In collaboration with Ralf Hiptmair, Christoph Schwab,
More informationScientific Computing WS 2018/2019. Lecture 15. Jürgen Fuhrmann Lecture 15 Slide 1
Scientific Computing WS 2018/2019 Lecture 15 Jürgen Fuhrmann juergen.fuhrmann@wias-berlin.de Lecture 15 Slide 1 Lecture 15 Slide 2 Problems with strong formulation Writing the PDE with divergence and gradient
More informationFrom the Boundary Element Domain Decomposition Methods to Local Trefftz Finite Element Methods on Polyhedral Meshes
From the Boundary Element Domain Decomposition Methods to Local Trefftz Finite Element Methods on Polyhedral Meshes Dylan Copeland 1, Ulrich Langer 2, and David Pusch 3 1 Institute of Computational Mathematics,
More informationWAVELETS, SHEARLETS AND GEOMETRIC FRAMES: PART II
WAVELETS, SHEARLETS AND GEOMETRIC FRAMES: PART II Philipp Grohs 1 and Axel Obermeier 2 October 22, 2014 1 ETH Zürich 2 ETH Zürich, supported by SNF grant 146356 OUTLINE 3. Curvelets, shearlets and parabolic
More informationLECTURE # 0 BASIC NOTATIONS AND CONCEPTS IN THE THEORY OF PARTIAL DIFFERENTIAL EQUATIONS (PDES)
LECTURE # 0 BASIC NOTATIONS AND CONCEPTS IN THE THEORY OF PARTIAL DIFFERENTIAL EQUATIONS (PDES) RAYTCHO LAZAROV 1 Notations and Basic Functional Spaces Scalar function in R d, d 1 will be denoted by u,
More informationIsogeometric Analysis:
Isogeometric Analysis: some approximation estimates for NURBS L. Beirao da Veiga, A. Buffa, Judith Rivas, G. Sangalli Euskadi-Kyushu 2011 Workshop on Applied Mathematics BCAM, March t0th, 2011 Outline
More informationLecture 8: Boundary Integral Equations
CBMS Conference on Fast Direct Solvers Dartmouth College June 23 June 27, 2014 Lecture 8: Boundary Integral Equations Gunnar Martinsson The University of Colorado at Boulder Research support by: Consider
More informationFrom the Boundary Element DDM to local Trefftz Finite Element Methods on Polyhedral Meshes
www.oeaw.ac.at From the Boundary Element DDM to local Trefftz Finite Element Methods on Polyhedral Meshes D. Copeland, U. Langer, D. Pusch RICAM-Report 2008-10 www.ricam.oeaw.ac.at From the Boundary Element
More informationLecture 9 Approximations of Laplace s Equation, Finite Element Method. Mathématiques appliquées (MATH0504-1) B. Dewals, C.
Lecture 9 Approximations of Laplace s Equation, Finite Element Method Mathématiques appliquées (MATH54-1) B. Dewals, C. Geuzaine V1.2 23/11/218 1 Learning objectives of this lecture Apply the finite difference
More informationSchwarz Preconditioner for the Stochastic Finite Element Method
Schwarz Preconditioner for the Stochastic Finite Element Method Waad Subber 1 and Sébastien Loisel 2 Preprint submitted to DD22 conference 1 Introduction The intrusive polynomial chaos approach for uncertainty
More informationThe Helmholtz Equation
The Helmholtz Equation Seminar BEM on Wave Scattering Rene Rühr ETH Zürich October 28, 2010 Outline Steklov-Poincare Operator Helmholtz Equation: From the Wave equation to Radiation condition Uniqueness
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 informationDeep Neural Networks and Partial Differential Equations: Approximation Theory and Structural Properties. Philipp Christian Petersen
Deep Neural Networks and Partial Differential Equations: Approximation Theory and Structural Properties Philipp Christian Petersen Joint work Joint work with: Helmut Bölcskei (ETH Zürich) Philipp Grohs
More informationAdaptive Petrov-Galerkin methods for first order transport equations. Author(s): Dahmen, Wolfgang; Huang, Chunyan; Schwab, Christoph; Welper, Gerrit
Research Collection Report Adaptive Petrov-Galerkin methods for first order transport equations Author(s): Dahmen, Wolfgang; Huang, Chunyan; Schwab, Christoph; Welper, Gerrit Publication Date: 2011 Permanent
More informationComputation of operators in wavelet coordinates
Computation of operators in wavelet coordinates Tsogtgerel Gantumur and Rob Stevenson Department of Mathematics Utrecht University Tsogtgerel Gantumur - Computation of operators in wavelet coordinates
More information1 Discretizing BVP with Finite Element Methods.
1 Discretizing BVP with Finite Element Methods In this section, we will discuss a process for solving boundary value problems numerically, the Finite Element Method (FEM) We note that such method is a
More informationOverlapping Schwarz preconditioners for Fekete spectral elements
Overlapping Schwarz preconditioners for Fekete spectral elements R. Pasquetti 1, L. F. Pavarino 2, F. Rapetti 1, and E. Zampieri 2 1 Laboratoire J.-A. Dieudonné, CNRS & Université de Nice et Sophia-Antipolis,
More informationA note on accurate and efficient higher order Galerkin time stepping schemes for the nonstationary Stokes equations
A note on accurate and efficient higher order Galerkin time stepping schemes for the nonstationary Stokes equations S. Hussain, F. Schieweck, S. Turek Abstract In this note, we extend our recent work for
More informationBasic Principles of Weak Galerkin Finite Element Methods for PDEs
Basic Principles of Weak Galerkin Finite Element Methods for PDEs Junping Wang Computational Mathematics Division of Mathematical Sciences National Science Foundation Arlington, VA 22230 Polytopal Element
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 informationIntegral Representation Formula, Boundary Integral Operators and Calderón projection
Integral Representation Formula, Boundary Integral Operators and Calderón projection Seminar BEM on Wave Scattering Franziska Weber ETH Zürich October 22, 2010 Outline Integral Representation Formula Newton
More informationHierarchical Parallel Solution of Stochastic Systems
Hierarchical Parallel Solution of Stochastic Systems Second M.I.T. Conference on Computational Fluid and Solid Mechanics Contents: Simple Model of Stochastic Flow Stochastic Galerkin Scheme Resulting Equations
More informationRobust exponential convergence of hp-fem for singularly perturbed systems of reaction-diffusion equations
Robust exponential convergence of hp-fem for singularly perturbed systems of reaction-diffusion equations Christos Xenophontos Department of Mathematics and Statistics University of Cyprus joint work with
More informationModelling in photonic crystal structures
Modelling in photonic crystal structures Kersten Schmidt MATHEON Nachwuchsgruppe Multiscale Modelling and Scientific Computing with PDEs in collaboration with Dirk Klindworth (MATHEON, TU Berlin) Holger
More informationWeight-adjusted DG methods for elastic wave propagation in arbitrary heterogeneous media
Weight-adjusted DG methods for elastic wave propagation in arbitrary heterogeneous media Jesse Chan Department of Computational and Applied Math, Rice University ICOSAHOM 2018 July 12, 2018 Chan (CAAM)
More informationFast and accurate methods for the discretization of singular integral operators given on surfaces
Fast and accurate methods for the discretization of singular integral operators given on surfaces James Bremer University of California, Davis March 15, 2018 This is joint work with Zydrunas Gimbutas (NIST
More informationNumerical methods for PDEs FEM convergence, error estimates, piecewise polynomials
Platzhalter für Bild, Bild auf Titelfolie hinter das Logo einsetzen Numerical methods for PDEs FEM convergence, error estimates, piecewise polynomials Dr. Noemi Friedman Contents of the course Fundamentals
More informationTrefftz-Discontinuous Galerkin Methods for the Time-Harmonic Maxwell Equations
Trefftz-Discontinuous Galerkin Methods for the Time-Harmonic Maxwell Equations Ilaria Perugia Dipartimento di Matematica - Università di Pavia (Italy) http://www-dimat.unipv.it/perugia Joint work with
More informationFinite Elements. Colin Cotter. February 22, Colin Cotter FEM
Finite Elements February 22, 2019 In the previous sections, we introduced the concept of finite element spaces, which contain certain functions defined on a domain. Finite element spaces are examples of
More informationFast Structured Spectral Methods
Spectral methods HSS structures Fast algorithms Conclusion Fast Structured Spectral Methods Yingwei Wang Department of Mathematics, Purdue University Joint work with Prof Jie Shen and Prof Jianlin Xia
More informationAlgorithms for Scientific Computing
Algorithms for Scientific Computing Hierarchical Methods and Sparse Grids d-dimensional Hierarchical Basis Michael Bader Technical University of Munich Summer 208 Intermezzo/ Big Picture : Archimedes Quadrature
More informationAdaptive algebraic multigrid methods in lattice computations
Adaptive algebraic multigrid methods in lattice computations Karsten Kahl Bergische Universität Wuppertal January 8, 2009 Acknowledgements Matthias Bolten, University of Wuppertal Achi Brandt, Weizmann
More informationOptimal multilevel preconditioning of strongly anisotropic problems.part II: non-conforming FEM. p. 1/36
Optimal multilevel preconditioning of strongly anisotropic problems. Part II: non-conforming FEM. Svetozar Margenov margenov@parallel.bas.bg Institute for Parallel Processing, Bulgarian Academy of Sciences,
More informationTensor Sparsity and Near-Minimal Rank Approximation for High-Dimensional PDEs
Tensor Sparsity and Near-Minimal Rank Approximation for High-Dimensional PDEs Wolfgang Dahmen, RWTH Aachen Collaborators: Markus Bachmayr, Ron DeVore, Lars Grasedyck, Endre Süli Paris, Oct.11, 2013 W.
More informationFast evaluation of boundary integral operators arising from an eddy current problem. MPI MIS Leipzig
1 Fast evaluation of boundary integral operators arising from an eddy current problem Steffen Börm MPI MIS Leipzig Jörg Ostrowski ETH Zürich Induction heating 2 Induction heating: Time harmonic eddy current
More informationStabilization and Acceleration of Algebraic Multigrid Method
Stabilization and Acceleration of Algebraic Multigrid Method Recursive Projection Algorithm A. Jemcov J.P. Maruszewski Fluent Inc. October 24, 2006 Outline 1 Need for Algorithm Stabilization and Acceleration
More informationNewton-Multigrid Least-Squares FEM for S-V-P Formulation of the Navier-Stokes Equations
Newton-Multigrid Least-Squares FEM for S-V-P Formulation of the Navier-Stokes Equations A. Ouazzi, M. Nickaeen, S. Turek, and M. Waseem Institut für Angewandte Mathematik, LSIII, TU Dortmund, Vogelpothsweg
More informationMultilevel stochastic collocations with dimensionality reduction
Multilevel stochastic collocations with dimensionality reduction Ionut Farcas TUM, Chair of Scientific Computing in Computer Science (I5) 27.01.2017 Outline 1 Motivation 2 Theoretical background Uncertainty
More informationAdaptive Finite Element Methods Lecture 1: A Posteriori Error Estimation
Adaptive Finite Element Methods Lecture 1: A Posteriori Error Estimation Department of Mathematics and Institute for Physical Science and Technology University of Maryland, USA www.math.umd.edu/ rhn 7th
More informationfor compression of Boundary Integral Operators. Steven Paul Nixon B.Sc.
Theory and Applications of the Multiwavelets for compression of Boundary Integral Operators. Steven Paul Nixon B.Sc. Institute for Materials Research School of Computing, Science & Engineering, University
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 informationTwo-scale Dirichlet-Neumann preconditioners for boundary refinements
Two-scale Dirichlet-Neumann preconditioners for boundary refinements Patrice Hauret 1 and Patrick Le Tallec 2 1 Graduate Aeronautical Laboratories, MS 25-45, California Institute of Technology Pasadena,
More informationTrefftz-Discontinuous Galerkin Methods for Maxwell s Equations
Trefftz-Discontinuous Galerkin Methods for Maxwell s Equations Ilaria Perugia Faculty of Mathematics, University of Vienna Vienna P D E Ralf Hiptmair (ETH Zürich), Andrea Moiola (University of Reading)
More informationKrylov subspace methods for linear systems with tensor product structure
Krylov subspace methods for linear systems with tensor product structure Christine Tobler Seminar for Applied Mathematics, ETH Zürich 19. August 2009 Outline 1 Introduction 2 Basic Algorithm 3 Convergence
More informationSpline Element Method for Partial Differential Equations
for Partial Differential Equations Department of Mathematical Sciences Northern Illinois University 2009 Multivariate Splines Summer School, Summer 2009 Outline 1 Why multivariate splines for PDEs? Motivation
More informationDomain Decomposition Preconditioners for Spectral Nédélec Elements in Two and Three Dimensions
Domain Decomposition Preconditioners for Spectral Nédélec Elements in Two and Three Dimensions Bernhard Hientzsch Courant Institute of Mathematical Sciences, New York University, 51 Mercer Street, New
More informationRobust Domain Decomposition Preconditioners for Abstract Symmetric Positive Definite Bilinear Forms
www.oeaw.ac.at Robust Domain Decomposition Preconditioners for Abstract Symmetric Positive Definite Bilinear Forms Y. Efendiev, J. Galvis, R. Lazarov, J. Willems RICAM-Report 2011-05 www.ricam.oeaw.ac.at
More informationStandard Finite Elements and Weighted Regularization
Standard Finite Elements and Weighted Regularization A Rehabilitation Martin COSTABEL & Monique DAUGE Institut de Recherche MAthématique de Rennes http://www.maths.univ-rennes1.fr/~dauge Slides of the
More informationMaximum norm estimates for energy-corrected finite element method
Maximum norm estimates for energy-corrected finite element method Piotr Swierczynski 1 and Barbara Wohlmuth 1 Technical University of Munich, Institute for Numerical Mathematics, piotr.swierczynski@ma.tum.de,
More informationVector and scalar penalty-projection methods
Numerical Flow Models for Controlled Fusion - April 2007 Vector and scalar penalty-projection methods for incompressible and variable density flows Philippe Angot Université de Provence, LATP - Marseille
More informationNumerical Solution I
Numerical Solution I Stationary Flow R. Kornhuber (FU Berlin) Summerschool Modelling of mass and energy transport in porous media with practical applications October 8-12, 2018 Schedule Classical Solutions
More informationA Robust Preconditioner for the Hessian System in Elliptic Optimal Control Problems
A Robust Preconditioner for the Hessian System in Elliptic Optimal Control Problems Etereldes Gonçalves 1, Tarek P. Mathew 1, Markus Sarkis 1,2, and Christian E. Schaerer 1 1 Instituto de Matemática Pura
More informationMatrix assembly by low rank tensor approximation
Matrix assembly by low rank tensor approximation Felix Scholz 13.02.2017 References Angelos Mantzaflaris, Bert Juettler, Boris Khoromskij, and Ulrich Langer. Matrix generation in isogeometric analysis
More informationA High-Order Galerkin Solver for the Poisson Problem on the Surface of the Cubed Sphere
A High-Order Galerkin Solver for the Poisson Problem on the Surface of the Cubed Sphere Michael Levy University of Colorado at Boulder Department of Applied Mathematics August 10, 2007 Outline 1 Background
More informationDiscontinuous Galerkin Methods
Discontinuous Galerkin Methods Joachim Schöberl May 20, 206 Discontinuous Galerkin (DG) methods approximate the solution with piecewise functions (polynomials), which are discontinuous across element interfaces.
More informationComputer simulation of multiscale problems
Progress in the SSF project CutFEM, Geometry, and Optimal design Computer simulation of multiscale problems Axel Målqvist and Daniel Elfverson University of Gothenburg and Uppsala University Umeå 2015-05-20
More informationChapter 12. Partial di erential equations Di erential operators in R n. The gradient and Jacobian. Divergence and rotation
Chapter 12 Partial di erential equations 12.1 Di erential operators in R n The gradient and Jacobian We recall the definition of the gradient of a scalar function f : R n! R, as @f grad f = rf =,..., @f
More informationOn the interplay between discretization and preconditioning of Krylov subspace methods
On the interplay between discretization and preconditioning of Krylov subspace methods Josef Málek and Zdeněk Strakoš Nečas Center for Mathematical Modeling Charles University in Prague and Czech Academy
More informationKarhunen-Loève Approximation of Random Fields Using Hierarchical Matrix Techniques
Institut für Numerische Mathematik und Optimierung Karhunen-Loève Approximation of Random Fields Using Hierarchical Matrix Techniques Oliver Ernst Computational Methods with Applications Harrachov, CR,
More information1. Let a(x) > 0, and assume that u and u h are the solutions of the Dirichlet problem:
Mathematics Chalmers & GU TMA37/MMG800: Partial Differential Equations, 011 08 4; kl 8.30-13.30. Telephone: Ida Säfström: 0703-088304 Calculators, formula notes and other subject related material are not
More informationFinal Ph.D. Progress Report. Integration of hp-adaptivity with a Two Grid Solver: Applications to Electromagnetics. David Pardo
Final Ph.D. Progress Report Integration of hp-adaptivity with a Two Grid Solver: Applications to Electromagnetics. David Pardo Dissertation Committee: I. Babuska, L. Demkowicz, C. Torres-Verdin, R. Van
More informationA Finite Element Method for the Surface Stokes Problem
J A N U A R Y 2 0 1 8 P R E P R I N T 4 7 5 A Finite Element Method for the Surface Stokes Problem Maxim A. Olshanskii *, Annalisa Quaini, Arnold Reusken and Vladimir Yushutin Institut für Geometrie und
More informationADAPTIVE ANISOTROPIC PETROV-GALERKIN METHODS FOR FIRST ORDER TRANSPORT EQUATIONS
ADAPTIVE ANISOTROPIC PETROV-GALERKIN METHODS FOR FIRST ORDER TRANSPORT EQUATIONS W. DAHMEN, G. KUTYNIOK, W.-Q LIM, C. SCHWAB, AND G. WELPER Abstract. This paper builds on recent developments of adaptive
More informationINTRODUCTION TO FINITE ELEMENT METHODS
INTRODUCTION TO FINITE ELEMENT METHODS LONG CHEN Finite element methods are based on the variational formulation of partial differential equations which only need to compute the gradient of a function.
More informationDual-Primal Isogeometric Tearing and Interconnecting Solvers for Continuous and Discontinuous Galerkin IgA Equations
Dual-Primal Isogeometric Tearing and Interconnecting Solvers for Continuous and Discontinuous Galerkin IgA Equations Christoph Hofer and Ulrich Langer Doctoral Program Computational Mathematics Numerical
More informationFast Sparse Spectral Methods for Higher Dimensional PDEs
Fast Sparse Spectral Methods for Higher Dimensional PDEs Jie Shen Purdue University Collaborators: Li-Lian Wang, Haijun Yu and Alexander Alekseenko Research supported by AFOSR and NSF ICERM workshop, June
More informationFast Algorithms for the Computation of Oscillatory Integrals
Fast Algorithms for the Computation of Oscillatory Integrals Emmanuel Candès California Institute of Technology EPSRC Symposium Capstone Conference Warwick Mathematics Institute, July 2009 Collaborators
More informationProjected Schur Complement Method for Solving Non-Symmetric Saddle-Point Systems (Arising from Fictitious Domain Approaches)
Projected Schur Complement Method for Solving Non-Symmetric Saddle-Point Systems (Arising from Fictitious Domain Approaches) Jaroslav Haslinger, Charles University, Prague Tomáš Kozubek, VŠB TU Ostrava
More informationSolving PDEs with Multigrid Methods p.1
Solving PDEs with Multigrid Methods Scott MacLachlan maclachl@colorado.edu Department of Applied Mathematics, University of Colorado at Boulder Solving PDEs with Multigrid Methods p.1 Support and Collaboration
More informationHybridized Discontinuous Galerkin Methods
Hybridized Discontinuous Galerkin Methods Theory and Christian Waluga Joint work with Herbert Egger (Uni Graz) 1st DUNE User Meeting, Stuttgart Christian Waluga (AICES) HDG Methods October 6-8, 2010 1
More informationhp-dg timestepping for time-singular parabolic PDEs and VIs
hp-dg timestepping for time-singular parabolic PDEs and VIs Ch. Schwab ( SAM, ETH Zürich ) (joint with O. Reichmann) ACMAC / Heraklion, Crete September 26-28, 2011 September 27, 2011 Motivation Well-posedness
More information1. Introduction. The Stokes problem seeks unknown functions u and p satisfying
A DISCRETE DIVERGENCE FREE WEAK GALERKIN FINITE ELEMENT METHOD FOR THE STOKES EQUATIONS LIN MU, JUNPING WANG, AND XIU YE Abstract. A discrete divergence free weak Galerkin finite element method is developed
More informationSpherical Wave Approach (SWA) for the electromagnetic modelling of 3D GPR scenarios. Nicola Tedeschi
Spherical Wave Approach (SWA) for the electromagnetic modelling of 3D GPR scenarios Nicola Tedeschi Outline Solutions of the vector wave equation Spherical wave approach Spherical reflected and transmitted
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 informationIterative methods for positive definite linear systems with a complex shift
Iterative methods for positive definite linear systems with a complex shift William McLean, University of New South Wales Vidar Thomée, Chalmers University November 4, 2011 Outline 1. Numerical solution
More informationAdaptive Finite Element Methods Lecture Notes Winter Term 2017/18. R. Verfürth. Fakultät für Mathematik, Ruhr-Universität Bochum
Adaptive Finite Element Methods Lecture Notes Winter Term 2017/18 R. Verfürth Fakultät für Mathematik, Ruhr-Universität Bochum Contents Chapter I. Introduction 7 I.1. Motivation 7 I.2. Sobolev and finite
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 informationNumerical Analysis of Differential Equations Numerical Solution of Elliptic Boundary Value
Numerical Analysis of Differential Equations 188 5 Numerical Solution of Elliptic Boundary Value Problems 5 Numerical Solution of Elliptic Boundary Value Problems TU Bergakademie Freiberg, SS 2012 Numerical
More informationTrefftz-DG solution to the Helmholtz equation involving integral equations
Trefftz-DG solution to the Helmholtz equation involving integral equations H. Barucq, A. Bendali, M. Fares, V. Mattesi, S. Tordeux Magique 3D Inria Bordeaux Sud Ouest LMA UMR CNRS 5142 INSA of Toulouse
More informationSparse Tensor Methods for PDEs with Stochastic Data II: Convergence Rates of gpcfem
Sparse Tensor Methods for PDEs with Stochastic Data II: Convergence Rates of gpcfem Christoph Schwab Seminar for Applied Mathematics ETH Zürich Albert Cohen (Paris, France) & Ron DeVore (Texas A& M) R.
More informationarxiv: v2 [math.na] 8 Sep 2017
arxiv:1704.06339v [math.na] 8 Sep 017 A Monte Carlo approach to computing stiffness matrices arising in polynomial chaos approximations Juan Galvis O. Andrés Cuervo September 3, 018 Abstract We use a Monte
More informationEfficient Augmented Lagrangian-type Preconditioning for the Oseen Problem using Grad-Div Stabilization
Efficient Augmented Lagrangian-type Preconditioning for the Oseen Problem using Grad-Div Stabilization Timo Heister, Texas A&M University 2013-02-28 SIAM CSE 2 Setting Stationary, incompressible flow problems
More informationThe All-floating BETI Method: Numerical Results
The All-floating BETI Method: Numerical Results Günther Of Institute of Computational Mathematics, Graz University of Technology, Steyrergasse 30, A-8010 Graz, Austria, of@tugraz.at Summary. The all-floating
More informationAMS 529: Finite Element Methods: Fundamentals, Applications, and New Trends
AMS 529: Finite Element Methods: Fundamentals, Applications, and New Trends Lecture 3: Finite Elements in 2-D Xiangmin Jiao SUNY Stony Brook Xiangmin Jiao Finite Element Methods 1 / 18 Outline 1 Boundary
More informationIsogeometric mortaring
Isogeometric mortaring E. Brivadis, A. Buffa, B. Wohlmuth, L. Wunderlich IMATI E. Magenes - Pavia Technical University of Munich A. Buffa (IMATI-CNR Italy) IGA mortaring 1 / 29 1 Introduction Splines Approximation
More informationNonlinear, Transient Conduction Heat Transfer Using A Discontinuous Galerkin Hierarchical Finite Element Method
Nonlinear, Transient Conduction Heat Transfer Using A Discontinuous Galerkin Hierarchical Finite Element Method by Jerome Charles Sanders B.S. in Physics, May 2002 The College of New Jersey A Thesis submitted
More informationBoundary Integral Equations on the Sphere with Radial Basis Functions: Error Analysis
Boundary Integral Equations on the Sphere with Radial Basis Functions: Error Analysis T. Tran Q. T. Le Gia I. H. Sloan E. P. Stephan Abstract Radial basis functions are used to define approximate solutions
More informationContents as of 12/8/2017. Preface. 1. Overview...1
Contents as of 12/8/2017 Preface 1. Overview...1 1.1 Introduction...1 1.2 Finite element data...1 1.3 Matrix notation...3 1.4 Matrix partitions...8 1.5 Special finite element matrix notations...9 1.6 Finite
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 informationOUTLINE ffl CFD: elliptic pde's! Ax = b ffl Basic iterative methods ffl Krylov subspace methods ffl Preconditioning techniques: Iterative methods ILU
Preconditioning Techniques for Solving Large Sparse Linear Systems Arnold Reusken Institut für Geometrie und Praktische Mathematik RWTH-Aachen OUTLINE ffl CFD: elliptic pde's! Ax = b ffl Basic iterative
More informationChapter 1: The Finite Element Method
Chapter 1: The Finite Element Method Michael Hanke Read: Strang, p 428 436 A Model Problem Mathematical Models, Analysis and Simulation, Part Applications: u = fx), < x < 1 u) = u1) = D) axial deformation
More informationError analysis and fast solvers for high-frequency scattering problems
Error analysis and fast solvers for high-frequency scattering problems I.G. Graham (University of Bath) Woudschoten October 2014 High freq. problem for the Helmholtz equation Given an object Ω R d, with
More informationInstitut de Recherche MAthématique de Rennes
LMS Durham Symposium: Computational methods for wave propagation in direct scattering. - July, Durham, UK The hp version of the Weighted Regularization Method for Maxwell Equations Martin COSTABEL & Monique
More informationHigh order Galerkin appoximations for parametric second order elliptic partial differential equations
Research Collection Report High order Galerkin appoximations for parametric second order elliptic partial differential equations Author(s): Nistor, Victor; Schwab, Christoph Publication Date: 2012 Permanent
More informationA very short introduction to the Finite Element Method
A very short introduction to the Finite Element Method Till Mathis Wagner Technical University of Munich JASS 2004, St Petersburg May 4, 2004 1 Introduction This is a short introduction to the finite element
More informationStochastic multiscale modeling of subsurface and surface flows. Part III: Multiscale mortar finite elements for coupled Stokes-Darcy flows
Stochastic multiscale modeling of subsurface and surface flows. Part III: Multiscale mortar finite elements for coupled Stokes-Darcy flows Ivan otov Department of Mathematics, University of Pittsburgh
More information