hp-dg timestepping for time-singular parabolic PDEs and VIs
|
|
- Caren Barrett
- 5 years ago
- Views:
Transcription
1 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
2 Motivation Well-posedness Application to FBM models dg Timestepping Examples Generalizations Conclusions Ch. Schwab September 27, 2011 p. 2
3 FBM Definition H (0, 1) Hurst parameter. A continuous Gaussian process W H on the probability space (Ω,A,F,P) with W H (0) = 0, W H (t) being a random variable with vanishing mean for all t 0 and E[W(t) H W H (s)] = 1 2 (t2h +s 2H t s 2H ), s,t 0, is called fractional Brownian motion (fbm). H = 1/2: Brownian motion. Ch. Schwab September 27, 2011 p. 3
4 FBM: Backward equations Consider the following problem: let the evolution of the underlying S under a pricing measure Q be given by ds(t) = rs(t)dt+σs(t)dw H (t), S(0) = S 0 > 0, 0 t T, where r,σ > 0 are constants. Then the price P(t) of a contingent claim at time t with payoff g(x) is given by P(t) = e r(t t) Ẽ Q [g(s T ) F H t ], Ẽ Q quasi-conditional expectation under the measure Q. Proof. See [Benth 03] or [Necula 02]. Ch. Schwab September 27, 2011 p. 4
5 FBM: Backward equations Under suitable integrability and smoothness assumptions holds the following analog of the Feynman-Kac theorem ([Benth 03], [Elliot et al. 03]). Theorem Let u(t, s) denote the solution of the backward equation t u+hσ 2 t γ s 2 ss u+rs s u ru = 0 }{{} A(t)u u(t,s) = g(s), s R + where γ = 2H 1, H (0,1) (so that γ ( 1,1)). Then P(t) = u(t,s(t)). Ch. Schwab September 27, 2011 p. 5
6 Ex. 1: Time-inhomogeneous Diffusion (t γ b,t γ σ,0), b R d,γ ( 1,1),σ R d d >0. A BS (t)ϕ(x) = t γ b ϕ(x)+t γ1 2 tr[σd2 ϕ(x)]. where A is defined by for γ ] 1,1[. t u t γ Au = f on J D, (1) A := γ +1 2 u(0) = u 0, (2) n j,k=1 a j,k (x), x j x k Ch. Schwab September 27, 2011 p. 6
7 Ex. 2: Time-inhomogeneous Jump-Diffusion (d = 1) (0,t δα 1 σ 2,k(t,z)) z t δ k(t,z) = t δα 1e M z 1+α, σ > 0, δ (0,1), M > 1, α (0,2) and t (0,T). A(t) of X(t) reads A(t) = A BS (t)+a J (t), ( A J (t)ϕ(x) = ϕ(x+z) ϕ(x) zϕ (x) ) k(t,z)dz, A BS (t)ϕ(x) = t δα 11 2 σ2 ϕ (x). Ch. Schwab September 27, 2011 p. 7
8 Ex. 2: Time-inhomogeneous Jump-Diffusion (d = 1) Weak formulation: Parametric bilinear forms a(, ), a J (, ) and a BS (, ) for ϕ(x),ψ(x) C 0 (Rd ) are given for 0 < t < T by Questions: a(t;ϕ,ψ) = a Diff (t;ϕ,ψ) +a J (t;ϕ,ψ), a J (t;ϕ,ψ) = (A J (t)ϕ,ψ), a Diff (t;ϕ,ψ) = (A BS (t)ϕ,ψ). meaning of (, )? H < 1/2 = γ = 2H 1 < 0: a Diff (t;ϕ,ϕ) as t 0. Ch. Schwab September 27, 2011 p. 8
9 Abstract forward problem For a bounded Lipschitz domain D R d, d 1 and a finite time interval I := (0, T), consider forward parabolic problem: find u such that t u t γ Au = f on I D, u(0) = g, A L(V,V ) nondegenerate, elliptic in V H H V. Ch. Schwab September 27, 2011 p. 9
10 Abstract forward problem Function Spaces: D R d bdd., Lipschitz, J = (0,T), T <. X := H 1 t γ/2 (J;V ) L 2 t γ/2 (J;V) = ( H 1 t γ/2 (J) V ) ( L 2 t γ/2 (J) V ), Y := L 2 t γ/2 (J;V) = L 2 t γ/2 (J) V, X (0 := {w X : w(0, ) = 0 in V }, X 0) := {w X : w(t, ) = 0 in V }, e.g. V := H 1 0(D), V = H 1 (D), L 2 t γ/2 (J) = C (J) L 2 t γ/2(j), H 1 t γ/2 (J) = C (J) H 1 t γ/2(j). Ch. Schwab September 27, 2011 p. 10
11 Well-posedness Possible approaches Semigroup approach: problems w. degeneracy in time. Lax-Milgram Lemma Works applied to space variable pointwise, but not uniformly in time Fails applied to space-time formulation inf sup or Babuska-Brezzi-conditions space-time variational formulation Ch. Schwab September 27, 2011 p. 11
12 BB-Conditions X, Y Hilbert spaces, B(, ) : X Y R bilinear form. BB-conditions: inf sup B(u, v) c stab > 0, (3) 0 u X 0 v Y u X v Y and 0 v Y : supb(u,v) > 0 (4) u X sup 0 u X,0 v Y B(u, v) u X v Y <. (5) Ch. Schwab September 27, 2011 p. 12
13 BB-Conditions Theorem Let B(, ) : X Y R satisfy (3)-(5). Then the problem: given f Y, find u X such that B(u,v) = f(v), v Y, (6) admits a unique variational solution u X with u X 1 c stab f Y. B(u,v) = T 0 (( u,v)+tγ a(u,v)) dt, a(u,v) = (Au,v), g = 0 (Homogeneous terminal conditions) Ch. Schwab September 27, 2011 p. 13
14 BB-Conditions For V := H0 1 (D) and L 2 t γ/2 (I) := C (0,1) L 2 t γ/2(i), u 2 L 2 t γ/2(i) := where X and Y given as B(u, v) inf sup 1 0 u X (0 0 v Y u X v Y 2 X := H 1 t γ/2 (I;V ) L 2 t γ/2 (I;V) = ( H 1 t γ/2 (I) V ) ( L 2 t γ/2 (I) V ), Y := L 2 t γ/2 (I;V) = L 2 t γ/2 (I) V, X (0 := {w X : w(0, ) = 0 in V }. I u 2 t γ dt. Ch. Schwab September 27, 2011 p. 14
15 Extensions Non-homogeneous terminal payoff functions. as natural BC s in (x, t) adjoint formulation B (u,v) = T for u Y, v X 0). 0 ( v,u)+t γ a(u,v) = f(v)+(v(0),g) The following embedding holds (P. Grisvard) where Λ = A 1/2 with X C 0 (I,D(Λ 1 2 γ 2 )), D(Λ 0 ) = V and D(Λ 1 ) = V. Ch. Schwab September 27, 2011 p. 15
16 Example: FBM (Reichmann 11) We consider the backward equation Transformations: t u+hσ 2 t 2H 1 S 2 SS u+rs S u ru = 0 u(t,s) = g(s). Time-reversal: t T τ Localization from R + to bounded domain /Weighted spaces Ch. Schwab September 27, 2011 p. 16
17 Localized forward formulation where B(v, w) = B(v,w) = f(v), T 0 (( v,w)+a(τ;v,w)) dτ a(τ;v,w) = σ2 2 (T τ)γ ( y v(τ,y), y w(τ,y)), f(v) = B( g, v) X := H 1 (T τ) γ/2 (I;V ) L 2 (T τ) γ/2 (I;V), Y := L 2 (T τ) γ/2 (I;V), V := H 1 0(D). Ch. Schwab September 27, 2011 p. 17
18 discontinuous Galerkin ( dg ) timestepping M = {I m } M m=1, partition of J = (0,T), r NM 0 dg orders. dg-timestepping: U V r (M;X) := {u : J V : u Im P rm (I m,v),m = 1 : M} such that for all v V r (M;V) B dg (U,v) = F dg (v), where M B dg (U,v) = (U,v) H dt+ I m + m=1 M m=1 I m t γ a(u,v)dt M ([U] m 1,v m 1 + ) H +(U 0 +,v+ 0 ) H m=2 F dg (v) = (u 0,v + 0 ) H + M m=1 I m (f(t),v) V Vdt Ch. Schwab September 27, 2011 p. 18
19 discontinuous Galerkin ( dg ) timestepping 1. Assume: Ex. 0 < α β < such that v V : α v 2 V V Av,v V β v 2 V then shall show that for a suitable dg Interpolant I holds ( u U L 2 C 1+ α ) u Iu t γ/2(i;v) β L 2 γ/2(i;v), t Ch. Schwab September 27, 2011 p. 19
20 discontinuous Galerkin ( dg ) timestepping 2. dg-interpolanti: Let I = (0,1). For u L 2 (I;V) continuous at t = 1 define Π r u P r (I,X), r 1, via the (Π r u u,q) H dt = 0, q P r 1 (I;V) and I Π r u(1 ) = u(1 ) V. Ch. Schwab September 27, 2011 p. 20
21 discontinuous Galerkin ( dg ) timestepping 3. Quasioptimality : u X (0 exact solution, U V r (M,V) time-semidiscrete dg solution. Let Iu V r (M,V) denote the dg interpolant of u be defined Then I m M : Iu Im = Π rm I m (u Im ). u U L 2 t γ/2(i;v) C for C > 0 independent of M and of r. ( 1+ α ) u Iu β L 2 γ/2(i;v), t Ch. Schwab September 27, 2011 p. 21
22 discontinuous Galerkin ( dg ) timestepping 4. scaled local dg interpolation error bounds Let I = (a,b), k = b a, r N 0 and u H s 0+1 for s 0 N 0. Then ex. C > 0 s.t. u Π r I u 2 L 2 (I,V) C Γ(r +1 s) max{1,r} 2 Γ(r +1+s) for any k > 0, r N 0, 0 s min{r,s 0 }, s real. ( ) k 2(s+1) u 2 H 2 s+1 (I,V), Ch. Schwab September 27, 2011 p. 22
23 discontinuous Galerkin ( dg ) timestepping 5. Solution regularity Assume u 0 H θ for 0 θ 1, g suff. regular A(t) t γ A, γ < 1 ( analytic case). Then u satisfies for 0 < a b < min(1,t) and for every k N u (k) (t) 2 H s t γ/2((a,b),v ) Cd 2(k+s) Γ(2k +3)a 2(k+s)+θ. Ch. Schwab September 27, 2011 p. 23
24 discontinuous Galerkin ( dg ) timestepping 6. Geometric Timesteps/ linear order vector: A geometric time partition M n,σ = {I m } n+1 m=1 with grading factor σ (0,1) and n+1 time steps I m, m = 1,2,...,n+1 is given by the nodes t 0 = 0, t m = Tσ n+1 m, 1 m n+1. Time step size: k m = I m = t m t m 1 k m = λt m 1, λ = 1 σ σ, 2 m M = n+1. Ch. Schwab September 27, 2011 p. 24
25 discontinuous Galerkin ( dg ) Timestepping 6. Geometric Timesteps/ linear order vector: A polynomial degree vector r = {r m } n+1 m=1 is called linear with slope ν 0 on the geometric partition M n,σ of (0,T) if r 1 = 0 and r m = νm for 2 m n+1. (Schötzau Diss. ETH 1999/ Schötzau and CS SINUM 2000). Ch. Schwab September 27, 2011 p. 25
26 discontinuous Galerkin ( dg ) timestepping 7. global dg interpolation error bounds There exist C,d > 0 such that I m M n,σ, 2 m n+1 and s m = α m r m with arbitrary α m (0,1), u Π rm I m u ( ) 2 L 2 Cσ(n m+2)θ γ (µd) 2αm(1 α m) 1 αm rm t γ/2(im,x) (1+α m ) 1+αm where µ = max{1,λ} and λ = 1 σ σ. C,d > 0 only depend on u 0 H θ and on γ ( 1,1). Ch. Schwab September 27, 2011 p. 26
27 discontinuous Galerkin ( dg ) Timestepping Theorem : [Reichmann & CS 2011] Consider the time-inhomogeneous forward problem on J = (0,1) with initial data u 0 H θ for some θ (0,1] and right hand side g. Discretize in time using dgfem on geometric partition M n,σ. Then for all degree vectors r = {r m } n+1 m=1 with slope ν ν 0 > 0 the semidiscrete dgfem solution U obtained in V r (M n,σ,v) converges exponentially w.r. to N, No. of time-dofs : u U L 2 t γ/2(j;v) C(σ,ν 0)exp( bn 1/2 ) Ch. Schwab September 27, 2011 p. 27
28 Examples Example 1: Exponential Convergence for fbm d = sigma=0.7 sigma=0.75 sigma=0.8 sigma=0.85 sigma= L error Number of time steps Figure: Convergence rates for different parameters σ and γ = 0.5. Ch. Schwab September 27, 2011 p. 28
29 Examples Example 2: Exponential Convergence time-inhomogeneous Lévy-like additive process 10 2 sigma=0.1 sigma=0.2 sigma=0.3 sigma=0.4 sigma= L error Number of time steps Figure: Convergence rates for different parameters σ and γ = 0. Ch. Schwab September 27, 2011 p. 29
30 Examples Example 2 (cont d): Exponential Convergence time-inhomogeneous Lévy-like additive process sigma=0.1 sigma=0.2 sigma=0.3 sigma=0.4 sigma=0.5 L error Number of time steps Figure: Convergence rates for different parameters σ and γ = 0.9. Ch. Schwab September 27, 2011 p. 30
31 Examples Example 3: algebraic convergence FBM, u 0 (x) = (x x 0 ) +, d = 1, dg(0) L2 finite element error Linfinity finite element error O(N 2 ) log error log number of mesh points Figure: Convergence rate in the L and L 2 norm at maturity of a plain vanilla European contract using dg(0) timestepping and linear C 0 -FEM in space. Hurst parameter H = 0.1, γ = 2H 1 = 0.8, graded mesh M (time grading parameter β = 5). Ch. Schwab September 27, 2011 p. 31
32 Examples Example 4: FBM PVI in X ( American Put ) H=0.1 H=0.3 H=0.5 H=0.7 H=0.9 Payoff 0.3 Option price Moneyness Figure: Prices of American puts in the FBM market model using dg(0) timestepping and linear FE in space: H [0.1,...,0.9] Ch. Schwab September 27, 2011 p. 32
33 Generalization γ-pseudohomogeneous, Gevrey time dependence: [ A BS (t)ϕ(x) = a(t) b ϕ(x)+ 1 ] 2 tr[σd2 ϕ(x)] where 0 < c 1 a(t)/t γ c 2 < in (0,T], and there exist d > 0 such that for some δ 1 k N : d k a(t) dt k (t) dk+1 (k!) δ t γ k, 0 < t T. Then u U L 2 t γ/2(j;v) C(q,ν 0)exp( bn 1/(1+δ) ) Ch. Schwab September 27, 2011 p. 33
34 Conclusions Well-posedness for time-degenerate parabolic problems hp-dg timestepping: exponential convergence in number N of space problems. (wavelet discretizations in (log-)price space: bounded condition number of matrices, ind. of dim.) h-dg: optimal conv. rates for time-graded M. time-degenerate PVIs: stable formulation, opt. rate for dg(0). Ch. Schwab September 27, 2011 p. 34
35 Conclusions Well-posedness for time-degenerate parabolic problems hp-dg timestepping: exponential convergence in number N of space problems. (wavelet discretizations in (log-)price space: bounded condition number of matrices, ind. of dim.) h-dg: optimal conv. rates for time-graded M. time-degenerate PVIs: stable formulation, opt. rate for dg(0). Thank you!!! Ch. Schwab September 27, 2011 p. 34
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 informationWELL POSEDNESS OF PROBLEMS I
Finite Element Method 85 WELL POSEDNESS OF PROBLEMS I Consider the following generic problem Lu = f, where L : X Y, u X, f Y and X, Y are two Banach spaces We say that the above problem is well-posed (according
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 informationFrom Completing the Squares and Orthogonal Projection to Finite Element Methods
From Completing the Squares and Orthogonal Projection to Finite Element Methods Mo MU Background In scientific computing, it is important to start with an appropriate model in order to design effective
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 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 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 informationVariational Formulations
Chapter 2 Variational Formulations In this chapter we will derive a variational (or weak) formulation of the elliptic boundary value problem (1.4). We will discuss all fundamental theoretical results that
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 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 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 informationELLIPTIC RECONSTRUCTION AND A POSTERIORI ERROR ESTIMATES FOR PARABOLIC PROBLEMS
ELLIPTIC RECONSTRUCTION AND A POSTERIORI ERROR ESTIMATES FOR PARABOLIC PROBLEMS CHARALAMBOS MAKRIDAKIS AND RICARDO H. NOCHETTO Abstract. It is known that the energy technique for a posteriori error analysis
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 informationPartial Differential Equations
Part II Partial Differential Equations Year 2015 2014 2013 2012 2011 2010 2009 2008 2007 2006 2005 2015 Paper 4, Section II 29E Partial Differential Equations 72 (a) Show that the Cauchy problem for u(x,
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 informationNumerische Mathematik
Numer. Math. (2003) 94: 195 202 Digital Object Identifier (DOI) 10.1007/s002110100308 Numerische Mathematik Some observations on Babuška and Brezzi theories Jinchao Xu, Ludmil Zikatanov Department of Mathematics,
More informationA Posteriori Error Estimation for hp-version Time-Stepping Methods for Parabolic Partial Differential Equations
Numerische Mathematik manuscript No. (will be inserted by the editor) A Posteriori Error Estimation for hp-version Time-Stepping Methods for Parabolic Partial Differential Equations Dominik Schötzau 1,
More informationGALERKIN TIME STEPPING METHODS FOR NONLINEAR PARABOLIC EQUATIONS
GALERKIN TIME STEPPING METHODS FOR NONLINEAR PARABOLIC EQUATIONS GEORGIOS AKRIVIS AND CHARALAMBOS MAKRIDAKIS Abstract. We consider discontinuous as well as continuous Galerkin methods for the time discretization
More informationDiscontinuous Galerkin Methods: Theory, Computation and Applications
Discontinuous Galerkin Methods: Theory, Computation and Applications Paola. Antonietti MOX, Dipartimento di Matematica Politecnico di Milano.MO. X MODELLISTICA E CALCOLO SCIENTIICO. MODELING AND SCIENTIIC
More informationPROBABILITY: LIMIT THEOREMS II, SPRING HOMEWORK PROBLEMS
PROBABILITY: LIMIT THEOREMS II, SPRING 218. HOMEWORK PROBLEMS PROF. YURI BAKHTIN Instructions. You are allowed to work on solutions in groups, but you are required to write up solutions on your own. Please
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 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 informationLecture 1. Finite difference and finite element methods. Partial differential equations (PDEs) Solving the heat equation numerically
Finite difference and finite element methods Lecture 1 Scope of the course Analysis and implementation of numerical methods for pricing options. Models: Black-Scholes, stochastic volatility, exponential
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 informationLecture Notes on PDEs
Lecture Notes on PDEs Alberto Bressan February 26, 2012 1 Elliptic equations Let IR n be a bounded open set Given measurable functions a ij, b i, c : IR, consider the linear, second order differential
More informationChapter 3 Conforming Finite Element Methods 3.1 Foundations Ritz-Galerkin Method
Chapter 3 Conforming Finite Element Methods 3.1 Foundations 3.1.1 Ritz-Galerkin Method Let V be a Hilbert space, a(, ) : V V lr a bounded, V-elliptic bilinear form and l : V lr a bounded linear functional.
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 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 informationMaximum-norm a posteriori estimates for discontinuous Galerkin methods
Maximum-norm a posteriori estimates for discontinuous Galerkin methods Emmanuil Georgoulis Department of Mathematics, University of Leicester, UK Based on joint work with Alan Demlow (Kentucky, USA) DG
More informationConvergence and optimality of an adaptive FEM for controlling L 2 errors
Convergence and optimality of an adaptive FEM for controlling L 2 errors Alan Demlow (University of Kentucky) joint work with Rob Stevenson (University of Amsterdam) Partially supported by NSF DMS-0713770.
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 informationSECOND ORDER TIME DISCONTINUOUS GALERKIN METHOD FOR NONLINEAR CONVECTION-DIFFUSION PROBLEMS
Proceedings of ALGORITMY 2009 pp. 1 10 SECOND ORDER TIME DISCONTINUOUS GALERKIN METHOD FOR NONLINEAR CONVECTION-DIFFUSION PROBLEMS MILOSLAV VLASÁK Abstract. We deal with a numerical solution of a scalar
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 informationSpace time finite and boundary element methods
Space time finite and boundary element methods Olaf Steinbach Institut für Numerische Mathematik, TU Graz http://www.numerik.math.tu-graz.ac.at based on joint work with M. Neumüller, H. Yang, M. Fleischhacker,
More informationPROBABILITY: LIMIT THEOREMS II, SPRING HOMEWORK PROBLEMS
PROBABILITY: LIMIT THEOREMS II, SPRING 15. HOMEWORK PROBLEMS PROF. YURI BAKHTIN Instructions. You are allowed to work on solutions in groups, but you are required to write up solutions on your own. Please
More informationMultigrid Methods for Saddle Point Problems
Multigrid Methods for Saddle Point Problems Susanne C. Brenner Department of Mathematics and Center for Computation & Technology Louisiana State University Advances in Mathematics of Finite Elements (In
More informationChapter 1 Foundations of Elliptic Boundary Value Problems 1.1 Euler equations of variational problems
Chapter 1 Foundations of Elliptic Boundary Value Problems 1.1 Euler equations of variational problems Elliptic boundary value problems often occur as the Euler equations of variational problems the latter
More informationDomain Perturbation for Linear and Semi-Linear Boundary Value Problems
CHAPTER 1 Domain Perturbation for Linear and Semi-Linear Boundary Value Problems Daniel Daners School of Mathematics and Statistics, The University of Sydney, NSW 2006, Australia E-mail: D.Daners@maths.usyd.edu.au
More informationUniformly Uniformly-ergodic Markov chains and BSDEs
Uniformly Uniformly-ergodic Markov chains and BSDEs Samuel N. Cohen Mathematical Institute, University of Oxford (Based on joint work with Ying Hu, Robert Elliott, Lukas Szpruch) Centre Henri Lebesgue,
More informationExponential multistep methods of Adams-type
Exponential multistep methods of Adams-type Marlis Hochbruck and Alexander Ostermann KARLSRUHE INSTITUTE OF TECHNOLOGY (KIT) 0 KIT University of the State of Baden-Wuerttemberg and National Laboratory
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 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 informationHigh order Cubature on Wiener space applied to derivative pricing and sensitivity estimations
1/25 High order applied to derivative pricing and sensitivity estimations University of Oxford 28 June 2010 2/25 Outline 1 2 3 4 3/25 Assume that a set of underlying processes satisfies the Stratonovich
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 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 informationProjected Surface Finite Elements for Elliptic Equations
Available at http://pvamu.edu/aam Appl. Appl. Math. IN: 1932-9466 Vol. 8, Issue 1 (June 2013), pp. 16 33 Applications and Applied Mathematics: An International Journal (AAM) Projected urface Finite Elements
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 informationThe Discontinuous Galerkin Finite Element Method
The Discontinuous Galerkin Finite Element Method Michael A. Saum msaum@math.utk.edu Department of Mathematics University of Tennessee, Knoxville The Discontinuous Galerkin Finite Element Method p.1/41
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 informationDyson series for the PDEs arising in Mathematical Finance I
for the PDEs arising in Mathematical Finance I 1 1 Penn State University Mathematical Finance and Probability Seminar, Rutgers, April 12, 2011 www.math.psu.edu/nistor/ This work was supported in part by
More informationSparse Tensor Galerkin Discretizations for First Order Transport Problems
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
More informationThe Azéma-Yor Embedding in Non-Singular Diffusions
Stochastic Process. Appl. Vol. 96, No. 2, 2001, 305-312 Research Report No. 406, 1999, Dept. Theoret. Statist. Aarhus The Azéma-Yor Embedding in Non-Singular Diffusions J. L. Pedersen and G. Peskir Let
More informationOptimal Left and Right Additive Schwarz Preconditioning for Minimal Residual Methods with Euclidean and Energy Norms
Optimal Left and Right Additive Schwarz Preconditioning for Minimal Residual Methods with Euclidean and Energy Norms Marcus Sarkis Worcester Polytechnic Inst., Mass. and IMPA, Rio de Janeiro and Daniel
More informationPART IV Spectral Methods
PART IV Spectral Methods Additional References: R. Peyret, Spectral methods for incompressible viscous flow, Springer (2002), B. Mercier, An introduction to the numerical analysis of spectral methods,
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 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 informationDiscontinuous Galerkin Time Discretization Methods for Parabolic Problems with Linear Constraints
J A N U A R Y 0 1 8 P R E P R N T 4 7 4 Discontinuous Galerkin Time Discretization Methods for Parabolic Problems with Linear Constraints gor Voulis * and Arnold Reusken nstitut für Geometrie und Praktische
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 informationR T (u H )v + (2.1) J S (u H )v v V, T (2.2) (2.3) H S J S (u H ) 2 L 2 (S). S T
2 R.H. NOCHETTO 2. Lecture 2. Adaptivity I: Design and Convergence of AFEM tarting with a conforming mesh T H, the adaptive procedure AFEM consists of loops of the form OLVE ETIMATE MARK REFINE to produce
More informationAn introduction to the mathematical theory of finite elements
Master in Seismic Engineering E.T.S.I. Industriales (U.P.M.) Discretization Methods in Engineering An introduction to the mathematical theory of finite elements Ignacio Romero ignacio.romero@upm.es October
More informationChapter 5 A priori error estimates for nonconforming finite element approximations 5.1 Strang s first lemma
Chapter 5 A priori error estimates for nonconforming finite element approximations 51 Strang s first lemma We consider the variational equation (51 a(u, v = l(v, v V H 1 (Ω, and assume that the conditions
More informationEULER MARUYAMA APPROXIMATION FOR SDES WITH JUMPS AND NON-LIPSCHITZ COEFFICIENTS
Qiao, H. Osaka J. Math. 51 (14), 47 66 EULER MARUYAMA APPROXIMATION FOR SDES WITH JUMPS AND NON-LIPSCHITZ COEFFICIENTS HUIJIE QIAO (Received May 6, 11, revised May 1, 1) Abstract In this paper we show
More informationA BIVARIATE SPLINE METHOD FOR SECOND ORDER ELLIPTIC EQUATIONS IN NON-DIVERGENCE FORM
A BIVARIAE SPLINE MEHOD FOR SECOND ORDER ELLIPIC EQUAIONS IN NON-DIVERGENCE FORM MING-JUN LAI AND CHUNMEI WANG Abstract. A bivariate spline method is developed to numerically solve second order elliptic
More informationIterative Methods for Linear Systems
Iterative Methods for Linear Systems 1. Introduction: Direct solvers versus iterative solvers In many applications we have to solve a linear system Ax = b with A R n n and b R n given. If n is large the
More informationSupraconvergence of a Non-Uniform Discretisation for an Elliptic Third-Kind Boundary-Value Problem with Mixed Derivatives
Supraconvergence of a Non-Uniform Discretisation for an Elliptic Third-Kind Boundary-Value Problem with Mixed Derivatives Etienne Emmrich Technische Universität Berlin, Institut für Mathematik, Straße
More informationhp-version Discontinuous Galerkin Finite Element Methods for Semilinear Parabolic Problems
Report no. 3/11 hp-version Discontinuous Galerkin Finite Element Methods for Semilinear Parabolic Problems Andris Lasis Endre Süli We consider the hp version interior penalty discontinuous Galerkin finite
More information10 The Finite Element Method for a Parabolic Problem
1 The Finite Element Method for a Parabolic Problem In this chapter we consider the approximation of solutions of the model heat equation in two space dimensions by means of Galerkin s method, using piecewise
More informationMesh Grading towards Singular Points Seminar : Elliptic Problems on Non-smooth Domain
Mesh Grading towards Singular Points Seminar : Elliptic Problems on Non-smooth Domain Stephen Edward Moore Johann Radon Institute for Computational and Applied Mathematics Austrian Academy of Sciences,
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 informationSpectral Properties of Elliptic Operators In previous work we have replaced the strong version of an elliptic boundary value problem
Spectral Properties of Elliptic Operators In previous work we have replaced the strong version of an elliptic boundary value problem L u x f x BC u x g x with the weak problem find u V such that B u,v
More informationSimple Examples on Rectangular Domains
84 Chapter 5 Simple Examples on Rectangular Domains In this chapter we consider simple elliptic boundary value problems in rectangular domains in R 2 or R 3 ; our prototype example is the Poisson equation
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 informationFINITE ELEMENT METHODS
FINITE ELEMENT METHODS Lecture notes arxiv:1709.08618v1 [math.na] 25 Sep 2017 Christian Clason September 25, 2017 christian.clason@uni-due.de https://udue.de/clason CONTENTS I BACKGROUND 1 overview of
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 informationDISCRETE MAXIMUM PRINCIPLES in THE FINITE ELEMENT SIMULATIONS
DISCRETE MAXIMUM PRINCIPLES in THE FINITE ELEMENT SIMULATIONS Sergey Korotov BCAM Basque Center for Applied Mathematics http://www.bcamath.org 1 The presentation is based on my collaboration with several
More informationState Space Representation of Gaussian Processes
State Space Representation of Gaussian Processes Simo Särkkä Department of Biomedical Engineering and Computational Science (BECS) Aalto University, Espoo, Finland June 12th, 2013 Simo Särkkä (Aalto University)
More information[2] (a) Develop and describe the piecewise linear Galerkin finite element approximation of,
269 C, Vese Practice problems [1] Write the differential equation u + u = f(x, y), (x, y) Ω u = 1 (x, y) Ω 1 n + u = x (x, y) Ω 2, Ω = {(x, y) x 2 + y 2 < 1}, Ω 1 = {(x, y) x 2 + y 2 = 1, x 0}, Ω 2 = {(x,
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 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 informationNonlinear Analysis 71 (2009) Contents lists available at ScienceDirect. Nonlinear Analysis. journal homepage:
Nonlinear Analysis 71 2009 2744 2752 Contents lists available at ScienceDirect Nonlinear Analysis journal homepage: www.elsevier.com/locate/na A nonlinear inequality and applications N.S. Hoang A.G. Ramm
More informationStochastic Volatility and Correction to the Heat Equation
Stochastic Volatility and Correction to the Heat Equation Jean-Pierre Fouque, George Papanicolaou and Ronnie Sircar Abstract. From a probabilist s point of view the Twentieth Century has been a century
More informationSparse Quadrature Algorithms for Bayesian Inverse Problems
Sparse Quadrature Algorithms for Bayesian Inverse Problems Claudia Schillings, Christoph Schwab Pro*Doc Retreat Disentis 2013 Numerical Analysis and Scientific Computing Disentis - 15 August, 2013 research
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 informationON ADDITIVE TIME-CHANGES OF FELLER PROCESSES. 1. Introduction
ON ADDITIVE TIME-CHANGES OF FELLER PROCESSES ALEKSANDAR MIJATOVIĆ AND MARTIJN PISTORIUS Abstract. In this note we generalise the Phillips theorem [1] on the subordination of Feller processes by Lévy subordinators
More informationNumerical Solutions to Partial Differential Equations
Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University Sobolev Embedding Theorems Embedding Operators and the Sobolev Embedding Theorem
More informationLAN property for sde s with additive fractional noise and continuous time observation
LAN property for sde s with additive fractional noise and continuous time observation Eulalia Nualart (Universitat Pompeu Fabra, Barcelona) joint work with Samy Tindel (Purdue University) Vlad s 6th birthday,
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 informationDynamical Reconstruction of Multi-Dimensional Image Sequences Using Optical Flows
Dynamical Reconstruction of Multi-Dimensional Image Sequences Using Optical Flows Dissertation zur Erlangung des Doktorgrades der Fakultät für Mathematik, Informatik und Naturwissenschaften der 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 informationANALYTIC SEMIGROUPS AND APPLICATIONS. 1. Introduction
ANALYTIC SEMIGROUPS AND APPLICATIONS KELLER VANDEBOGERT. Introduction Consider a Banach space X and let f : D X and u : G X, where D and G are real intervals. A is a bounded or unbounded linear operator
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 informationChapter Two: Numerical Methods for Elliptic PDEs. 1 Finite Difference Methods for Elliptic PDEs
Chapter Two: Numerical Methods for Elliptic PDEs Finite Difference Methods for Elliptic PDEs.. Finite difference scheme. We consider a simple example u := subject to Dirichlet boundary conditions ( ) u
More informationApproximation in Banach Spaces by Galerkin Methods
2 Approximation in Banach Spaces by Galerkin Methods In this chapter, we consider an abstract linear problem which serves as a generic model for engineering applications. Our first goal is to specify the
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 informationWavelet-in-time multigrid-in-space preconditioning of parabolic evolution equations
Wavelet-in-time multigrid-in-space preconditioning of parabolic evolution equations Roman Andreev To cite this version: Roman Andreev. Wavelet-in-time multigrid-in-space preconditioning of parabolic evolution
More informationELEMENTS OF PROBABILITY THEORY
ELEMENTS OF PROBABILITY THEORY Elements of Probability Theory A collection of subsets of a set Ω is called a σ algebra if it contains Ω and is closed under the operations of taking complements and countable
More informationDiscretization of SDEs: Euler Methods and Beyond
Discretization of SDEs: Euler Methods and Beyond 09-26-2006 / PRisMa 2006 Workshop Outline Introduction 1 Introduction Motivation Stochastic Differential Equations 2 The Time Discretization of SDEs Monte-Carlo
More information{σ x >t}p x. (σ x >t)=e at.
3.11. EXERCISES 121 3.11 Exercises Exercise 3.1 Consider the Ornstein Uhlenbeck process in example 3.1.7(B). Show that the defined process is a Markov process which converges in distribution to an N(0,σ
More informationWavelet-in-time multigrid-in-space preconditioning of parabolic evolution equations
www.oeaw.ac.at Wavelet-in-time multigrid-in-space preconditioning of parabolic evolution equations R. Andreev RICAM-Report 2014-27 www.ricam.oeaw.ac.at WAVELET-IN-TIME MULTIGRID-IN-SPACE PRECONDITIONING
More informationLECTURE 3: DISCRETE GRADIENT FLOWS
LECTURE 3: DISCRETE GRADIENT FLOWS Department of Mathematics and Institute for Physical Science and Technology University of Maryland, USA Tutorial: Numerical Methods for FBPs Free Boundary Problems and
More information