A space-time Trefftz method for the second order wave equation
|
|
- Isaac Preston
- 5 years ago
- Views:
Transcription
1 A space-time Trefftz method for the second order wave equation Lehel Banjai The Maxwell Institute for Mathematical Sciences Heriot-Watt University, Edinburgh Rome, 10th Apr 2017 Joint work with: Emmanuil Georgoulis (Leicester & Athens), Oluwaseun F Lijoka (HW) 1 / 36
2 Outline of the talk 1 Motivation 2 An interior-penalty space-time dg method 3 Damped wave equation 4 Numerical results 2 / 36
3 Outline 1 Motivation 2 An interior-penalty space-time dg method 3 Damped wave equation 4 Numerical results 3 / 36
4 Acoustic wave equation Find u(t) H0 1 (Ω), t [0, T ], s.t. (ü, v) L 2 (Ω) + (a u, v) L 2 (Ω) = 0 u(x, 0) = u 0 (x), u(x, 0) = v 0 (x), in Ω. for all v H 1 0 (Ω), 4 / 36
5 Acoustic wave equation Find u(t) H0 1 (Ω), t [0, T ], s.t. (ü, v) L 2 (Ω) + (a u, v) L 2 (Ω) = 0 u(x, 0) = u 0 (x), u(x, 0) = v 0 (x), in Ω. for all v H 1 0 (Ω), Initial data u 0 H0 1(Ω), v 0 L 2 (Ω). a(x) piecewise constant 0 < c a < a(x) < C a. Unique solution exists with u L 2 ([0, T ]; H0 1 (Ω)), u L 2 ([0, T ]; L 2 (Ω)), ü L 2 ([0, T ]; H 1 (Ω)) u C([0, T ]; H 1 0 (Ω)), u C([0, T ]; L 2 (Ω)). For smooth enough solution we have the transmission conditions u j = u k, a j n u j = a k n u k, on Ω j Ω k, u j = u Ωj, u k = u Ωk, where Ω j and Ω k are subsets of Ω with a a k in Ω k and a a j in Ω j. 4 / 36
6 How to discretize the wave equation? The usual approach: Construct a spatial mesh and a corresponding spatially discrete space: locally polynomial, continuous or discontinuous across the boundaries of the spatial elements (the spatial skeleton). Finite difference approximation in time. Solution computed by time-stepping. 5 / 36
7 How to discretize the wave equation? The usual approach: Construct a spatial mesh and a corresponding spatially discrete space: locally polynomial, continuous or discontinuous across the boundaries of the spatial elements (the spatial skeleton). Finite difference approximation in time. Solution computed by time-stepping. In this talk (Trefftz): Construct a space-time mesh and corresponding fully discrete space: In each space time element exact (polynomial or non-polynomial) exact solution of the wave equation. Discontinous across the space-time skeleton. Solve either by time-stepping or as a large system. 5 / 36
8 Frequency domain motivation Frequency domain Take cue from frequency domain û(x) where a j are directions, a j = 1. Motivation û ω 2 û = 0: k j=1 f j e iωx a j, For large ω, minimize the number of degrees of freedom per wavelength. 6 / 36
9 Frequency domain motivation Frequency domain Take cue from frequency domain û(x, ω) where a j are directions, a j = 1. k j=1 ˆf j (ω)e iωx a j, 6 / 36
10 Frequency domain motivation Frequency domain Take cue from frequency domain û(x, ω) where a j are directions, a j = 1. Time-domain Time-domain equivalent u(x, t) k j=1 k j=1 ˆf j (ω)e iωx a j, f j (t x a j ) 6 / 36
11 Frequency domain motivation Frequency domain Take cue from frequency domain û(x, ω) where a j are directions, a j = 1. Time-domain Time-domain equivalent u(x, t) k j=1 k k j=1 ˆf j (ω)e iωx a j, f j (t x a j ) p α j,l (t x a j ) l. j=1 l=0 6 / 36
12 (A bit of) Literature on Trefftz methods for waves Plenty of literature in the frequency domain O. Cessenat and B. Després, Application of an ultra weak variational formulation of elliptic PDEs to the two-dimensional Helmholtz equation, SIAM J. Numer. Anal., (1998) R. Hiptmair, A. Moiola, I. Perugia, A survey of Trefftz methods for the Helmholtz equation. Springer Lect. Notes Comput. Sci. Eng., 2016, pp / 36
13 (A bit of) Literature on Trefftz methods for waves Plenty of literature in the frequency domain O. Cessenat and B. Després, Application of an ultra weak variational formulation of elliptic PDEs to the two-dimensional Helmholtz equation, SIAM J. Numer. Anal., (1998) R. Hiptmair, A. Moiola, I. Perugia, A survey of Trefftz methods for the Helmholtz equation. Springer Lect. Notes Comput. Sci. Eng., 2016, pp Fewer in time-domain S. Petersen, C. Farhat, and R. Tezaur, A space-time discontinuous Galerkin method for the solution of the wave equation in the time domain, Internat. J. Numer. Methods Engrg. (2009) F. Kretzschmar, A. Moiola, I. Perugia, S. M. Schnepp, A priori error analysis of space-time Trefftz discontinuous Galerkin methods for wave problems, IMA J. Numer. Anal., 36(4) 2016, pp L. Banjai, E. Georgoulis, O Lijoka, A Trefftz polynomial space-time discontinuous Galerkin method for the second order wave equation, SIAM J. Numer. Anal (2017), pp / 36
14 Outline 1 Motivation 2 An interior-penalty space-time dg method 3 Damped wave equation 4 Numerical results 8 / 36
15 DG setting Time discretization 0 = t 0 < t 1 < < t N = T, I n = [t n, t n+1 ]; τ n = t n+1 t n. Spatial-mesh T n of Ω consisting of open simplices such that Ω = K Tn K. In each simplex K, a(x) is constant. Space-time slabs T n I n, h-space-time meshwidth. The skeleton of the space mesh denoted Γ n and ˆΓ n = Γ n 1 Γ n. Usual jump and average definitions (e = K + K Γ int ) and if e K + Ω, {u} e = 1 2 (u+ + u ), {v} e = 1 2 (v+ + v ), [u] e = u + n + + u n, [v] e = v + n + + v n, {v} e = v +, [u] e = u + n + Also u(t n ) = u(t + n ) u(t n ), u(t 0 ) = u(t + 0 ). 9 / 36
16 Local Trefftz spaces The space of piecewise polynomials on the time-space mesh denoted S h,p n = {u L 2 (Ω I n ) u K In P p (R d+1 ), K T n }, Let S h,p n,trefftz S n h,p with v(t, x) a v(t, x) = 0, The dimensions of the spaces S h,p n t I n, x K, for any v S h,p n,trefftz. and S h,p n,trefftz for spatial dimension d are 1D 2D 3D 1 poly 2 (p + 1)(p + 2) 1 6 (p + 1)(p + 2)(p + 3) 2D 1 4 (p + 4) Trefftz 2p + 1 (p + 1) (p + 1)(p + 2)(2p + 3) We expect the approximation properties of solutions of the wave equation to be the same for the two spaces of different dimension. 10 / 36
17 Constructing the polynomial spaces Choose directions ξ j : (t 1 l x ξ j ), l = 0, 1,..., p. a 11 / 36
18 Constructing the polynomial spaces Choose directions ξ j : (t 1 l x ξ j ), l = 0, 1,..., p. a Alternatively propagate polynomial initial condition: u(0) = x α k, u(0) = 0, and u(0) = 0, u(0) = x β k, with α k p and β k p 1, α k, β k multi-indices. 11 / 36
19 Constructing the polynomial spaces Choose directions ξ j : (t 1 l x ξ j ), l = 0, 1,..., p. a Alternatively propagate polynomial initial condition: u(0) = x α k, u(0) = 0, and u(0) = 0, u(0) = x β k, with α k p and β k p 1, α k, β k multi-indices. Important observation: trunctation of the Taylor expansion of exact solution is a polynomial solution of the wave equation. 11 / 36
20 The space on Ω [0, T ] is then defined as V h,p Trefftz = {u L2 (Ω [0, T ]) u Ω In S h,p n,trefftz, n = 0, 1..., N}. 12 / 36
21 The space on Ω [0, T ] is then defined as V h,p Trefftz = {u L2 (Ω [0, T ]) u Ω In S h,p n,trefftz, n = 0, 1..., N}. (Abuse of) Notation: u h V h,p Trefftz-discrete function on Ω [0, T ] u n S h,p n,trefftz, restriction of u on Ω I n. u ex the exact solution. 12 / 36
22 An interior penalty dg method We start with t n t n+1 [ Ω ü v + a u vdx Γ {a u} [ v] ds Γ [u] {a v} ds Testing with v = u gives where the energy is given by t n t n+1 + σ 0 Γ [u] [ v] ds ] dt = 0. d E(t, u)dt = 0, dt E(t, u) = 1 2 u(t) 2 Ω a u(t) 2 Ω σ 0 [u(t)] 2 Γ Γ {a u} [u] ds. Discrete inverse inequality in space and usual choice of penalty parameter gives E(t, u) / 36
23 Jumps in time Summing over n gives E(tN N 1 ) E(t+ 0 ) n=1 E(t n ) = 0. To give a sign to the extra terms (Hughes, Hulbert 88): 1 2 ( u(t n), u(t n )) L 2 (Ω) ( u(t n ), u(t + n )) L 2 (Ω) = 1 2 ( u(t n), u(t n ) ) L 2 (Ω). Do this for all the terms, including the stabilization. Obtain a dissipative method. 14 / 36
24 Space-time dg formulation N 1 a(u, v) = (ü, v) Ω In + ( u(t n ), v(t n + )) Ω n=0 + (a u, v) Ω In + ( a u(t n ), v(t n + )) Ω ({a u}, [ v]) Γn I n ( {a u(t n )}, [v(t + n )])ˆΓ n ([u], {a v}) Γn I n ( [u(t n )], {a v(t + n )})ˆΓn + (σ 0 [u], [ v]) Γn I n + (σ 0 [u(t n )], [v(t + n )])ˆΓ n and + (σ 1 [u], [v]) Γn I n + (σ 2 [a u], [a v]) Γn I n b init (v) = (v 0, v(t + 0 )) Ω + (a u 0, v(t + 0 )) Ω ({a u 0 }, [v(t + 0 )]) Γ0 ([u 0 ], {a v(t + 0 )}) Γ0 + (σ 0 [u 0 ], [v(t + 0 )]) Γ0. Find u h V h,p Trefftz (Ω [0, T ]) such that a(u h, v) = b init (v), v V h,p Trefftz (Ω [0, T ]). 15 / 36
25 Time-space dg as a time-stepping method a n (u, v) = (ü, v) Ω In + ( u(t + n ), v(t + n )) Ω + (a u, v) Ω In + (a u(t + n ), v(t + n )) Ω ({a u}, [ v]) Γn I n ({a u(t + n )}, [v(t + n )]) Γn ([u], {a v}) Γn I n ([u(t + n )], {a v(t + n )}) Γn + (σ 0 [u], [ v]) Γn I n + (σ 0 [u(t + n )], [v(t + n )]) Γn + (σ 1 [u], [v]) Γn I n + (σ 2 [a u], [a v]) Γn I n, b n (u, v) = ( u(t n ), v(t + n )) Ω + (a u(t n ), v(t + n )) Ω ({a u(t n )}, [v(t + n )]) Γn ([u(t n )], {a v(t + n )}) Γn 1 + (σ 0 [u(t n )], [v(t + n )])ˆΓ n, Find u n S h,p n,trefftz such that a n (u n, v) = b n (u n 1, v), v S h,p n,trefftz. 16 / 36
26 Consistency and stability Theorem The following statements hold: 1 The method is consistent for a sufficiently smooth solution u. 2 There exists a choice of σ 0 h 1, such that for any v S h,p n,trefftz and t I n the energy is bounded below as E(t, v) 1 2 v(t) 2 L 2 (Ω) a v(t) 2 L 2 (Ω). 3 Let u h V h,p Trefftz discrete solution. Then E(t N, uh ) E(t 1, u h ). 17 / 36
27 a(, ) 1/2 = - a norm on V h,p Trefftz a(w, w) = E h (tn, w) + E h(t 0 + N 1, w) + ( 1 2 ẇ(t n) 2 Ω a w(t n ) 2 Ω n=1 ( {a w(t n )}, [w(t n )] )ˆΓ n σ 0 [w(t n )] 2ˆΓn ) N 1 + n=0 ( σ 1 [w] 2 Γ n I n + σ 2 [a w] 2 Γ n I n ). 18 / 36
28 a(, ) 1/2 = - a norm on V h,p Trefftz a(w, w) = E h (tn, w) + E h(t 0 + N 1, w) + ( 1 2 ẇ(t n) 2 Ω a w(t n ) 2 Ω n=1 Theorem ( {a w(t n )}, [w(t n )] )ˆΓ n σ 0 [w(t n )] 2ˆΓn ) N 1 + n=0 ( σ 1 [w] 2 Γ n I n + σ 2 [a w] 2 Γ n I n ). a(v, v) 1/2 = v = 0 v = 0, for v V h,p Trefftz. Hence, the time-space dg method a(u, v) = b init (v), v V h,p Trefftz has a unique solution in V h,p Trefftz. The proof is by noticing that if v = 0 then v is a smooth solution of the wave equation, uniquely determined by the initial condition. 18 / 36
29 Convergence analysis If we prove continuity of a(, ) a(u, v) C u v, u cont. sol. + V h,p Trefftz, v V h,p Trefftz, 19 / 36
30 Convergence analysis If we prove continuity of a(, ) a(u, v) C u v, u cont. sol. + V h,p Trefftz, v V h,p Trefftz, we can use Galerkin orthogonality to show, for any v V h,p Trefftz u h v 2 = a(u h v, u h v) and hence we have quasi-optimality = a(u ex v, u h v) C u ex v u h v u h u ex inf u h v + v u ex v V h,p Trefftz inf v V h,p Trefftz v u ex + C v u ex 19 / 36
31 Integrating by parts a few times (this is how to implement the method) a(w, v) = Recall N 1 n=0 ( ({a ẇ}, [v]) Γn In (σ 0 [ẇ], [v]) Γn In ({ẇ}, [a v]) Γ int n I n + (σ 1 [w], [v]) Γn I n + (σ 2 [a w], [a v]) Γn I n ) N ( (ẇ(tn ), v(t n ) ) Ω + (a w(tn ), v(t n ) ) Ω n=1 ({a w(t n )}, [v(t n )] ) Γn ([w(t n )], {a v(t n )} ) Γn + (σ 0 [w(t n )], [v(t n )] ) Γn ). v 2 = E(t N, v) + E(t+ 0, v) + N 1 ( 1 2 v(t n) 2 Ω v(t n) 2 Ω n=1 + ( { v(t n )}, [v(t n )] )ˆΓn σ 0 [v(t n )] 2ˆΓ n ) N 1 + n=0 ( σ 1 [v] 2 Γ I n + σ 2 [ v] 2 Γ I n ), 20 / 36
32 Integrating by parts a few times (this is how to implement the method) a(w, v) = N 1 n=0 ( ({a ẇ}, [v]) Γn In (σ 0 [ẇ], [v]) Γn In ({ẇ}, [a v]) Γ int n I n + (σ 1 [w], [v]) Γn I n + (σ 2 [a w], [a v]) Γn I n ) N ( (ẇ(tn ), v(t n ) ) Ω + (a w(tn ), v(t n ) ) Ω n=1 Hence define, w 2 = 1 2 ({a w(t n )}, [v(t n )] ) Γn ([w(t n )], {a v(t n )} ) Γn + (σ 0 [w(t n )], [v(t n )] ) Γn ). N n=1 N 1 + n=0 ( ẇ(t n ) 2 Ω + a w(t n ) 2 Ω + σ 0 [w(t n )] 2 Γ n + σ 1/2 0 {a w(t ( σ 1 [w] 2 Γ n I n + σ 2 [a w] 2 Γ n I n + σ 1/2 2 {ẇ} 2 Γ int + σ 1/2 1 {a ẇ} 2 Γ n I n + σ 0 σ 1/2 1 [ẇ] 2 Γ n I n ). n I n 20 / 36
33 The choice of stabilization parameters and convergence Let τ n = t n+1 t n, h = diam(k), (x, t) K (t n, t n+1 ), K T n. diam(k)/ρ K c T, K T n, n = 0, 1,..., N 1, where ρ K is the radius of the inscribed circle of K. Assume space-time elements star-shaped with respect to a ball. Choice of parameters σ 0 = p 2 c T Ca 2 C inv (c a h) 1, σ 1 = C a p 3 (hτ n ) 1 σ2 = h(c a τ n ) 1. Theorem For sufficiently smooth solution and h τ u h u ex = O(h p 1/2 ). Proof uses truncated Taylor expansion. 21 / 36
34 Error estimate in mesh independent norm Using a Gronwall argument we can show for v S h,p n,trefftz v 2 Ω I n + a v 2 Ω I n τ n e C(t n+1 t n)/h n ( v(t n+1 ) 2 Ω + a v(t n+1) 2 Ω ), h n = min x Ω h(x, t), t I n. Let τ = max τ n and h = min h n. Then V 2 Ω (0,T ) + a V 2 Ω (0,T ) Cτe Cτ/ h V 2, V V h,p Trefftz. Proposition u h u ex 2 Ω [0,T ] + uh u ex 2 Ω [0,T ] C inf (τe Cτ/ h V u ex 2 V V h,p Trefftz + V u ex 2 Ω [0,T ] + a (V u ex ) 2 Ω [0,T ] ). Hence in mesh independent energy norm we expect error O(h p ). 22 / 36
35 Outline 1 Motivation 2 An interior-penalty space-time dg method 3 Damped wave equation 4 Numerical results 23 / 36
36 Wave equation with damping Damped wave equation: ü + α u u = 0. The extra term decreases the energy: d dt E(t) = α u 2, hence only a minor modification to the DG formulation needed. 24 / 36
37 Wave equation with damping Damped wave equation: ü + α u u = 0. The extra term decreases the energy: d dt E(t) = α u 2, hence only a minor modification to the DG formulation needed. However: truncations of the Taylor expansion are no longer solutions! 24 / 36
38 Wave equation with damping Damped wave equation: ü + α u u = 0. The extra term decreases the energy: d dt E(t) = α u 2, hence only a minor modification to the DG formulation needed. However: truncations of the Taylor expansion are no longer solutions! Grysa, Maciag, Adamczyk-Krasa 14 consider solutions of the form e αt p 1 (x, t) + p 2 (x, t) with p 1 and p 2 polynomial in x and t. 24 / 36
39 Solution formula in 1D Instead, use basis functions obtained by propagating polynomial initial data: u(x, 0) = u 0 (x) = x α j, u(x, 0) = v 0 (x) = 0, and with u 0 (x) = 0, v 0 (x) = x β j α j p β j p / 36
40 Solution formula in 1D Instead, use basis functions obtained by propagating polynomial initial data: u(x, 0) = u 0 (x) = x α j, u(x, 0) = v 0 (x) = 0, and with u 0 (x) = 0, v 0 (x) = x β j α j p β j p 1. In 1D solution is then given by the d Alambert-like formula u(x, t) = 1 2 [u 0(x t) + u 0 (x + t)] e αt/2 + α 4 e αt/ e αt/2 x+t x t x+t x t u 0 (s) {I 0 (ρ(s) α 2 ) + t ρ(s) I 1 (ρ(s) α )} 2 ds v 0 (s)i 0 (ρ(s) α 2 ) ds, ρ(s; x, t) = t 2 (x s) / 36
41 Solution formula ctd. Rearranging (the last term) te αt/2 [v 0(x + st) + v 0 (x st)] I 0 ( αt s 2 ) ds. For example for v 0 (x) = x 2 1 x 2 te αt/2 I 0 ( αt s 2 ) ds + t 3 e αt/2 0 s2 I 0 ( αt 2 1 s 2 ) ds / 36
42 Solution formula ctd. Rearranging (the last term) te αt/2 [v 0(x + st) + v 0 (x st)] I 0 ( αt s 2 ) ds. For example for v 0 (x) = x 2 1 x 2 te αt/2 I 0 ( αt s 2 ) ds + t 3 e αt/2 0 s2 I 0 ( αt 2 1 s 2 ) ds. Corresponding term in 3D 1 te αt/2 0 [v 0 (x + tsz) + Dv 0 (x + tsz)] zds z I 0 ( αt B(0,1) 2 1 s 2 ) ds 1 26 / 36
43 Solution formula ctd. Rearranging (the last term) te αt/2 [v 0(x + st) + v 0 (x st)] I 0 ( αt s 2 ) ds. For example for v 0 (x) = x 2 1 x 2 te αt/2 I 0 ( αt s 2 ) ds + t 3 e αt/2 0 s2 I 0 ( αt 2 1 s 2 ) ds. Corresponding term in 3D 1 te αt/2 0 [v 0 (x + tsz) + Dv 0 (x + tsz)] zds z I 0 ( αt B(0,1) 2 1 s 2 ) ds Hence need efficient representation of functions of the type 1 ϱ j (t) = 0 sj I 0 ( αt 2 1 s 2 ) ds 1 26 / 36
44 Damped wave equation: Efficient implementation Need efficient representation of functions of the type 1 ϱ j (t) = 0 sj I 0 ( αt 2 1 s 2 ) ds In Matlab: Chebfun works well. Note computations are done on a small space-time element, i.e. t (0, h) for a small h > 0. Can use truncated Taylor expansions of analytic functions I 0 and I 1 with efficiency increasing for decreasing h. We expect similar effect for any lower order terms. 27 / 36
45 Changes to the analysis Polynomial in space the same discrete inverse inequalities used + the extra term decreases energy choice of parameters, stability and quasi-optimality proof identical. Approximation properties of the discrete space (in 1D): Away from the boundary, in each space-time element K (t, t + ) project solution to K and neighbouring elements at time t to polynomials and propagate. At boundary, extend exact solution anti-symmetrically and again project and propagate. 28 / 36
46 Outline 1 Motivation 2 An interior-penalty space-time dg method 3 Damped wave equation 4 Numerical results 29 / 36
47 One dimensional setting Simple 1D setting: Ω = (0, 1), a 1. Initial data u 0 = e x 5/8 ( ) 2 δ, v 0 = 0, δ δ 0 = Interested in having few degrees of freedom for decreasing δ δ t = 1.5, δ = δ t = 1.5, δ = δ 0 / / 36
48 One dimensional setting Simple 1D setting: Ω = (0, 1), a 1. Initial data u 0 = e x 5/8 ( ) 2 δ, v 0 = 0, δ δ 0 = Energy of exact solution exact energy = 1 2 u x(x, 0) 2 Ω 2δ 1 y 2 e 2y 2 π dy = δ We compare with full polynomial space. Note for polynomial order p 2p + 1 Trefftz 1 2 (p + 1)(p + 2)full polynomial space. In all 1D experiments square space-time elements. 30 / 36
49 Error in dg-norm, δ = δ 0, T = 1/4: Error Error Error Trefftz poly Full poly p = 1 p = 2 p = 3 p = 4 p = h Trefftz p-convergence (fixed h in space and time): p = 1 p = 2 p = 3 p = 4 p = h p 31 / 36
50 Energy Energy Energy conservation For δ = δ 0 /4: Trefftz poly Full poly Energy, p=1 Energy, p=2 5 Energy, p=3 Energy, p=4 Exact Energy Time 10 Energy, p=1 Energy, p=2 5 Energy, p=3 Energy, p=4 Exact Energy Time 32 / 36
51 Error Relative error for decreasing δ error δ = ( δ 2 u(, T ) u h(, T ) 2 Ω + δ 2 u(, T ) u h(, T ) 2 Ω ) 1/ p=2, / = / 0 p=3, / = / 0 p=4,/ = / 0 p=5, / = / 0 p=2, / = / 0=2 p=3, / = / 0=2 p=4, / = / 0=2 p=5, / = / 0=2 p=2, / = / 0=4 p=3, / = / 0=4 p=4, / = / 0=4 p=5, / = / 0= h=/ 33 / 36
52 Error Error 2D experiment On square [0, 1] 2 with exact solution Energy of error at final time: u(x, y, t) = cos( 2πt) sin πx sin πy. error = ( 1 2 u(, T ) u h(, T ) 2 Ω u(, T ) u h(, T ) 2 Ω )1/2. Trefftz poly Full poly p = 1 p = 2 p = 3 p = Mesh-size p = 1 p = 2 p = 3 p = Mesh-size 34 / 36
53 Error Error Damped wave equation in 1D Error in dg norm. Exact solution on Ω = (0, 1): u(x, t) = e ( αt/2) sin(πx) cos π 2 α2 4 t + α 2 π 2 α 2 /4 sin π 2 α2 4 t. Trefftz poly Full poly p = 1 p = 2 p = 3 p = p = 1 p = 2 p = 3 p = Number of DOF Number of DOF 35 / 36
54 Conclusions A space-time interior penalty dg method for the acoustic wave equation in second order form: Can be considered either as implicit time-stepping method or a large space-time system. Allows Trefftz basis functions, polynomial in space. Fewer degrees of freedom than full polynomial method and requires integration only over the space-time skeleton. Stability and convergence analysis available. Also practical for damped wave equation. 36 / 36
55 Conclusions A space-time interior penalty dg method for the acoustic wave equation in second order form: To do: Can be considered either as implicit time-stepping method or a large space-time system. Allows Trefftz basis functions, polynomial in space. Fewer degrees of freedom than full polynomial method and requires integration only over the space-time skeleton. Stability and convergence analysis available. Also practical for damped wave equation. A posteriori error analysis Adaptivity in h, p, wave directions. Tent-pitching meshes allow quasi explicit time-stepping. p-analysis in higher dimension. Applications. 36 / 36
A space-time Trefftz method for the second order wave equation
A space-time Trefftz method for the second order wave equation Lehel Banjai The Maxwell Institute for Mathematical Sciences Heriot-Watt University, Edinburgh & Department of Mathematics, University of
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 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 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 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 informationApproximation by plane and circular waves
LMS EPSRC DURHAM SYMPOSIUM, 8 16TH JULY 2014 Building Bridges: Connections and Challenges in Modern Approaches to Numerical PDEs Approximation by plane and circular waves Andrea Moiola DEPARTMENT OF MATHEMATICS
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 informationA phase-based interior penalty discontinuous Galerkin method for the Helmholtz equation with spatially varying wavenumber
A phase-based interior penalty discontinuous Galerkin method for the Helmholtz equation with spatially varying wavenumber Chi Yeung Lam * and Chi-Wang Shu ** * Department of Mathematics, The Chinese University
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 informationA posteriori analysis of a discontinuous Galerkin scheme for a diffuse interface model
A posteriori analysis of a discontinuous Galerkin scheme for a diffuse interface model Jan Giesselmann joint work with Ch. Makridakis (Univ. of Sussex), T. Pryer (Univ. of Reading) 9th DFG-CNRS WORKSHOP
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 informationSuperconvergence of discontinuous Galerkin methods for 1-D linear hyperbolic equations with degenerate variable coefficients
Superconvergence of discontinuous Galerkin methods for -D linear hyperbolic equations with degenerate variable coefficients Waixiang Cao Chi-Wang Shu Zhimin Zhang Abstract In this paper, we study the superconvergence
More informationA local-structure-preserving local discontinuous Galerkin method for the Laplace equation
A local-structure-preserving local discontinuous Galerkin method for the Laplace equation Fengyan Li 1 and Chi-Wang Shu 2 Abstract In this paper, we present a local-structure-preserving local discontinuous
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 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 informationThe solution of time harmonic wave equations using complete families of elementary solutions
The solution of time harmonic wave equations using complete families of elementary solutions Peter Monk Department of Mathematical Sciences University of Delaware Newark, DE 19716 USA email: monk@math.udel.edu
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 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 informationarxiv: v1 [math.na] 19 Nov 2018
QUASINORMS IN SEMILINEAR ELLIPTIC PROBLEMS JAMES JACKAMAN AND TRISTAN PRYER arxiv:1811.07816v1 [math.na] 19 Nov 2018 (1.1) Abstract. In this note we examine the a priori and a posteriori analysis of discontinuous
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 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 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 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 informationDiscontinuous Galerkin methods for fractional diffusion problems
Discontinuous Galerkin methods for fractional diffusion problems Bill McLean Kassem Mustapha School of Maths and Stats, University of NSW KFUPM, Dhahran Leipzig, 7 October, 2010 Outline Sub-diffusion Equation
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 informationNumerical Methods for PDEs
Numerical Methods for PDEs Partial Differential Equations (Lecture 1, Week 1) Markus Schmuck Department of Mathematics and Maxwell Institute for Mathematical Sciences Heriot-Watt University, Edinburgh
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 informationPhD dissertation defense
Isogeometric mortar methods with applications in contact mechanics PhD dissertation defense Brivadis Ericka Supervisor: Annalisa Buffa Doctor of Philosophy in Computational Mechanics and Advanced Materials,
More informationplasmas Lise-Marie Imbert-Gérard, Bruno Després. July 27th 2011 Laboratoire J.-L. LIONS, Université Pierre et Marie Curie, Paris.
Lise-Marie Imbert-Gérard, Bruno Després. Laboratoire J.-L. LIONS, Université Pierre et Marie Curie, Paris. July 27th 2011 1 Physical and mathematical motivations 2 approximation of the solutions 3 4 Plan
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 informationGraded Mesh Refinement and Error Estimates of Higher Order for DGFE-solutions of Elliptic Boundary Value Problems in Polygons
Graded Mesh Refinement and Error Estimates of Higher Order for DGFE-solutions of Elliptic Boundary Value Problems in Polygons Anna-Margarete Sändig, Miloslav Feistauer University Stuttgart, IANS Journées
More informationSpatial and Modal Superconvergence of the Discontinuous Galerkin Method for Linear Equations
Spatial and Modal Superconvergence of the Discontinuous Galerkin Method for Linear Equations N. Chalmers and L. Krivodonova March 5, 014 Abstract We apply the discontinuous Galerkin finite element method
More informationA LAGRANGE MULTIPLIER METHOD FOR ELLIPTIC INTERFACE PROBLEMS USING NON-MATCHING MESHES
A LAGRANGE MULTIPLIER METHOD FOR ELLIPTIC INTERFACE PROBLEMS USING NON-MATCHING MESHES P. HANSBO Department of Applied Mechanics, Chalmers University of Technology, S-4 96 Göteborg, Sweden E-mail: hansbo@solid.chalmers.se
More informationAn a posteriori error indicator for discontinuous Galerkin discretizations of H(curl) elliptic partial differential equations
IMA Journal of Numerical Analysis 005) Page of 7 doi: 0.093/imanum/ An a posteriori error indicator for discontinuous Galerkin discretizations of Hcurl) elliptic partial differential equations PAUL HOUSTON
More informationSpace-time XFEM for two-phase mass transport
Space-time XFEM for two-phase mass transport Space-time XFEM for two-phase mass transport Christoph Lehrenfeld joint work with Arnold Reusken EFEF, Prague, June 5-6th 2015 Christoph Lehrenfeld EFEF, Prague,
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 informationFinite element approximation of the stochastic heat equation with additive noise
p. 1/32 Finite element approximation of the stochastic heat equation with additive noise Stig Larsson p. 2/32 Outline Stochastic heat equation with additive noise du u dt = dw, x D, t > u =, x D, t > u()
More informationRESIDUAL BASED ERROR ESTIMATES FOR THE SPACE-TIME DISCONTINUOUS GALERKIN METHOD APPLIED TO NONLINEAR HYPERBOLIC EQUATIONS
Proceedings of ALGORITMY 2016 pp. 113 124 RESIDUAL BASED ERROR ESTIMATES FOR THE SPACE-TIME DISCONTINUOUS GALERKIN METHOD APPLIED TO NONLINEAR HYPERBOLIC EQUATIONS VÍT DOLEJŠÍ AND FILIP ROSKOVEC Abstract.
More informationBUBBLE STABILIZED DISCONTINUOUS GALERKIN METHOD FOR PARABOLIC AND ELLIPTIC PROBLEMS
BUBBLE STABILIZED DISCONTINUOUS GALERKIN METHOD FOR PARABOLIC AND ELLIPTIC PROBLEMS ERIK BURMAN AND BENJAMIN STAMM Abstract. In this paper we give an analysis of a bubble stabilized discontinuous Galerkin
More informationLocal discontinuous Galerkin methods for elliptic problems
COMMUNICATIONS IN NUMERICAL METHODS IN ENGINEERING Commun. Numer. Meth. Engng 2002; 18:69 75 [Version: 2000/03/22 v1.0] Local discontinuous Galerkin methods for elliptic problems P. Castillo 1 B. Cockburn
More informationUne méthode de pénalisation par face pour l approximation des équations de Navier-Stokes à nombre de Reynolds élevé
Une méthode de pénalisation par face pour l approximation des équations de Navier-Stokes à nombre de Reynolds élevé CMCS/IACS Ecole Polytechnique Federale de Lausanne Erik.Burman@epfl.ch Méthodes Numériques
More informationDISCONTINUOUS GALERKIN FINITE ELEMENT METHOD FOR THE WAVE EQUATION. SIAM J. Numer. Anal., Vol. 44, pp , 2006
DISCONTINUOUS GALERIN INITE ELEMENT METHOD OR THE WAVE EQUATION MARCUS J. GROTE, ANNA SCHNEEBELI, AND DOMINI SCHÖTZAU SIAM J. Numer. Anal., Vol. 44, pp. 408-43, 006 Abstract. The symmetric interior penalty
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 informationAn introduction to Birkhoff normal form
An introduction to Birkhoff normal form Dario Bambusi Dipartimento di Matematica, Universitá di Milano via Saldini 50, 0133 Milano (Italy) 19.11.14 1 Introduction The aim of this note is to present an
More informationTong Sun Department of Mathematics and Statistics Bowling Green State University, Bowling Green, OH
Consistency & Numerical Smoothing Error Estimation An Alternative of the Lax-Richtmyer Theorem Tong Sun Department of Mathematics and Statistics Bowling Green State University, Bowling Green, OH 43403
More informationPIECEWISE LINEAR FINITE ELEMENT METHODS ARE NOT LOCALIZED
PIECEWISE LINEAR FINITE ELEMENT METHODS ARE NOT LOCALIZED ALAN DEMLOW Abstract. Recent results of Schatz show that standard Galerkin finite element methods employing piecewise polynomial elements of degree
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 informationQualifying Examination
Summer 24 Day. Monday, September 5, 24 You have three hours to complete this exam. Work all problems. Start each problem on a All problems are 2 points. Please email any electronic files in support of
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 informationHamburger Beiträge zur Angewandten Mathematik
Hamburger Beiträge zur Angewandten Mathematik Numerical analysis of a control and state constrained elliptic control problem with piecewise constant control approximations Klaus Deckelnick and Michael
More informationDiscontinuous Galerkin Method for interface problem of coupling different order elliptic equations
Discontinuous Galerkin Method for interface problem of coupling different order elliptic equations Igor Mozolevski, Endre Süli Federal University of Santa Catarina, Brazil Oxford University Computing Laboratory,
More informationFiltered scheme and error estimate for first order Hamilton-Jacobi equations
and error estimate for first order Hamilton-Jacobi equations Olivier Bokanowski 1 Maurizio Falcone 2 2 1 Laboratoire Jacques-Louis Lions, Université Paris-Diderot (Paris 7) 2 SAPIENZA - Università di Roma
More informationExistence and uniqueness of solutions for nonlinear ODEs
Chapter 4 Existence and uniqueness of solutions for nonlinear ODEs In this chapter we consider the existence and uniqueness of solutions for the initial value problem for general nonlinear ODEs. Recall
More informationDiscontinuous Galerkin and Finite Difference Methods for the Acoustic Equations with Smooth Coefficients. Mario Bencomo TRIP Review Meeting 2013
About Me Mario Bencomo Currently 2 nd year graduate student in CAAM department at Rice University. B.S. in Physics and Applied Mathematics (Dec. 2010). Undergraduate University: University of Texas at
More informationAn A Posteriori Error Estimate for Discontinuous Galerkin Methods
An A Posteriori Error Estimate for Discontinuous Galerkin Methods Mats G Larson mgl@math.chalmers.se Chalmers Finite Element Center Mats G Larson Chalmers Finite Element Center p.1 Outline We present an
More information1. Introduction. We consider the model problem that seeks an unknown function u = u(x) satisfying
A SIMPLE FINITE ELEMENT METHOD FOR LINEAR HYPERBOLIC PROBLEMS LIN MU AND XIU YE Abstract. In this paper, we introduce a simple finite element method for solving first order hyperbolic equations with easy
More informationTHE WAVE EQUATION. d = 1: D Alembert s formula We begin with the initial value problem in 1 space dimension { u = utt u xx = 0, in (0, ) R, (2)
THE WAVE EQUATION () The free wave equation takes the form u := ( t x )u = 0, u : R t R d x R In the literature, the operator := t x is called the D Alembertian on R +d. Later we shall also consider the
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 informationPREPRINT 2010:25. Fictitious domain finite element methods using cut elements: II. A stabilized Nitsche method ERIK BURMAN PETER HANSBO
PREPRINT 2010:25 Fictitious domain finite element methods using cut elements: II. A stabilized Nitsche method ERIK BURMAN PETER HANSBO Department of Mathematical Sciences Division of Mathematics CHALMERS
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 informationSolving Time-Harmonic Scattering Problems by the Ultra Weak Variational Formulation
Introduction Solving Time-Harmonic Scattering Problems by the Ultra Weak Variational Formulation Plane waves as basis functions Peter Monk 1 Tomi Huttunen 2 1 Department of Mathematical Sciences University
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 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 informationA Space-Time Boundary Element Method for the Wave 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 A Space-Time Boundary Element Method for the Wave Equation Marco Zank and Olaf Steinbach Institut für Numerische Mathematik Space-Time Methods for PDEs,
More informationOn a Discontinuous Galerkin Method for Surface PDEs
On a Discontinuous Galerkin Method for Surface PDEs Pravin Madhavan (joint work with Andreas Dedner and Bjo rn Stinner) Mathematics and Statistics Centre for Doctoral Training University of Warwick Applied
More informationOn Pressure Stabilization Method and Projection Method for Unsteady Navier-Stokes Equations 1
On Pressure Stabilization Method and Projection Method for Unsteady Navier-Stokes Equations 1 Jie Shen Department of Mathematics, Penn State University University Park, PA 1682 Abstract. We present some
More informationNUMERICAL STUDIES OF VARIATIONAL-TYPE TIME-DISCRETIZATION TECHNIQUES FOR TRANSIENT OSEEN PROBLEM
Proceedings of ALGORITMY 212 pp. 44 415 NUMERICAL STUDIES OF VARIATIONAL-TYPE TIME-DISCRETIZATION TECHNIQUES FOR TRANSIENT OSEEN PROBLEM NAVEED AHMED AND GUNAR MATTHIES Abstract. In this paper, we combine
More informationBUBBLE STABILIZED DISCONTINUOUS GALERKIN METHOD FOR STOKES PROBLEM
BUBBLE STABILIZED DISCONTINUOUS GALERKIN METHOD FOR STOKES PROBLEM ERIK BURMAN AND BENJAMIN STAMM Abstract. We propose a low order discontinuous Galerkin method for incompressible flows. Stability of the
More informationFinite difference method for heat equation
Finite difference method for heat equation Praveen. C praveen@math.tifrbng.res.in Tata Institute of Fundamental Research Center for Applicable Mathematics Bangalore 560065 http://math.tifrbng.res.in/~praveen
More information2tdt 1 y = t2 + C y = which implies C = 1 and the solution is y = 1
Lectures - Week 11 General First Order ODEs & Numerical Methods for IVPs In general, nonlinear problems are much more difficult to solve than linear ones. Unfortunately many phenomena exhibit nonlinear
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 informationKey words. Incompressible magnetohydrodynamics, mixed finite element methods, discontinuous Galerkin methods
A MIXED DG METHOD FOR LINEARIZED INCOMPRESSIBLE MAGNETOHYDRODYNAMICS PAUL HOUSTON, DOMINIK SCHÖTZAU, AND XIAOXI WEI Journal of Scientific Computing, vol. 40, pp. 8 34, 009 Abstract. We introduce and analyze
More informationApproximation of fluid-structure interaction problems with Lagrange multiplier
Approximation of fluid-structure interaction problems with Lagrange multiplier Daniele Boffi Dipartimento di Matematica F. Casorati, Università di Pavia http://www-dimat.unipv.it/boffi May 30, 2016 Outline
More informationOn an Approximation Result for Piecewise Polynomial Functions. O. Karakashian
BULLETIN OF THE GREE MATHEMATICAL SOCIETY Volume 57, 010 (1 7) On an Approximation Result for Piecewise Polynomial Functions O. arakashian Abstract We provide a new approach for proving approximation results
More informationNONLOCAL DIFFUSION EQUATIONS
NONLOCAL DIFFUSION EQUATIONS JULIO D. ROSSI (ALICANTE, SPAIN AND BUENOS AIRES, ARGENTINA) jrossi@dm.uba.ar http://mate.dm.uba.ar/ jrossi 2011 Non-local diffusion. The function J. Let J : R N R, nonnegative,
More informationSemi-Lagrangian Formulations for Linear Advection Equations and Applications to Kinetic Equations
Semi-Lagrangian Formulations for Linear Advection and Applications to Kinetic Department of Mathematical and Computer Science Colorado School of Mines joint work w/ Chi-Wang Shu Supported by NSF and AFOSR.
More informationFinite Difference and Finite Element Methods
Finite Difference and Finite Element Methods Georgy Gimel farb COMPSCI 369 Computational Science 1 / 39 1 Finite Differences Difference Equations 3 Finite Difference Methods: Euler FDMs 4 Finite Element
More informationAn Equal-order DG Method for the Incompressible Navier-Stokes Equations
An Equal-order DG Method for the Incompressible Navier-Stokes Equations Bernardo Cockburn Guido anschat Dominik Schötzau Journal of Scientific Computing, vol. 40, pp. 188 10, 009 Abstract We introduce
More informationScalar Conservation Laws and First Order Equations Introduction. Consider equations of the form. (1) u t + q(u) x =0, x R, t > 0.
Scalar Conservation Laws and First Order Equations Introduction. Consider equations of the form (1) u t + q(u) x =, x R, t >. In general, u = u(x, t) represents the density or the concentration of a physical
More informationSpace-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 informationC e n t r u m v o o r W i s k u n d e e n I n f o r m a t i c a
C e n t r u m v o o r W i s k u n d e e n I n f o r m a t i c a Modelling, Analysis and Simulation Modelling, Analysis and Simulation Bilinear forms for the recovery-based discontinuous Galerkin method
More informationOptimal Error Estimates for the hp Version Interior Penalty Discontinuous Galerkin Finite Element Method
Report no. 03/06 Optimal Error Estimates for the hp Version Interior Penalty Discontinuous Galerkin Finite Element Method Emmanuil H. Georgoulis and Endre Süli Oxford University Computing Laboratory Numerical
More informationFinite Differences for Differential Equations 28 PART II. Finite Difference Methods for Differential Equations
Finite Differences for Differential Equations 28 PART II Finite Difference Methods for Differential Equations Finite Differences for Differential Equations 29 BOUNDARY VALUE PROBLEMS (I) Solving a TWO
More informationIntegration of Vlasov-type equations
Alexander Ostermann University of Innsbruck, Austria Joint work with Lukas Einkemmer Verona, April/May 2017 Plasma the fourth state of matter 99% of the visible matter in the universe is made of plasma
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 informationApplied Math Qualifying Exam 11 October Instructions: Work 2 out of 3 problems in each of the 3 parts for a total of 6 problems.
Printed Name: Signature: Applied Math Qualifying Exam 11 October 2014 Instructions: Work 2 out of 3 problems in each of the 3 parts for a total of 6 problems. 2 Part 1 (1) Let Ω be an open subset of R
More informationarxiv: v1 [math.na] 21 Nov 2017
High Order Finite Difference Schemes for the Heat Equation Whose Convergence Rates are Higher Than Their Truncation Errors, A. Ditkowski arxiv:7.0796v [math.na] Nov 07 Abstract Typically when a semi-discrete
More informationStability of an abstract wave equation with delay and a Kelvin Voigt damping
Stability of an abstract wave equation with delay and a Kelvin Voigt damping University of Monastir/UPSAY/LMV-UVSQ Joint work with Serge Nicaise and Cristina Pignotti Outline 1 Problem The idea Stability
More informationRecovery-Based A Posteriori Error Estimation
Recovery-Based A Posteriori Error Estimation Zhiqiang Cai Purdue University Department of Mathematics, Purdue University Slide 1, March 2, 2011 Outline Introduction Diffusion Problems Higher Order Elements
More informationNumerical Analysis of Differential Equations Numerical Solution of Parabolic Equations
Numerical Analysis of Differential Equations 215 6 Numerical Solution of Parabolic Equations 6 Numerical Solution of Parabolic Equations TU Bergakademie Freiberg, SS 2012 Numerical Analysis of Differential
More informationA POSTERIORI ERROR ESTIMATES OF TRIANGULAR MIXED FINITE ELEMENT METHODS FOR QUADRATIC CONVECTION DIFFUSION OPTIMAL CONTROL PROBLEMS
A POSTERIORI ERROR ESTIMATES OF TRIANGULAR MIXED FINITE ELEMENT METHODS FOR QUADRATIC CONVECTION DIFFUSION OPTIMAL CONTROL PROBLEMS Z. LU Communicated by Gabriela Marinoschi In this paper, we discuss a
More informationA-posteriori error analysis of hp-version discontinuous Galerkin finite element methods for second-order quasilinear elliptic problems
Report no. NA-0/1 A-posteriori error analysis of hp-version discontinuous Galerkin finite element methods for second-order quasilinear elliptic problems Paul Houston 1, Endre Süli, and Thomas P. Wihler
More informationWEAK GALERKIN FINITE ELEMENT METHOD FOR SECOND ORDER PARABOLIC EQUATIONS
INERNAIONAL JOURNAL OF NUMERICAL ANALYSIS AND MODELING Volume 13, Number 4, Pages 525 544 c 216 Institute for Scientific Computing and Information WEAK GALERKIN FINIE ELEMEN MEHOD FOR SECOND ORDER PARABOLIC
More informationUltraconvergence of ZZ Patch Recovery at Mesh Symmetry Points
Ultraconvergence of ZZ Patch Recovery at Mesh Symmetry Points Zhimin Zhang and Runchang Lin Department of Mathematics, Wayne State University Abstract. The ultraconvergence property of the Zienkiewicz-Zhu
More informationPerformance comparison between hybridizable DG and classical DG methods for elastic waves simulation in harmonic domain
March 4-5, 2015 Performance comparison between hybridizable DG and classical DG methods for elastic waves simulation in harmonic domain M. Bonnasse-Gahot 1,2, H. Calandra 3, J. Diaz 1 and S. Lanteri 2
More informationA PENALTY FREE NON-SYMMETRIC NITSCHE TYPE METHOD FOR THE WEAK IMPOSITION OF BOUNDARY CONDITIONS
A PENALTY FREE NON-SYMMETRIC NITSCHE TYPE METHOD FOR THE WEAK IMPOSITION OF BOUNDARY CONDITIONS ERIK BURMAN Abstract. In this note we show that the non-symmetric version of the classical Nitsche s method
More informationLucio Demeio Dipartimento di Ingegneria Industriale e delle Scienze Matematiche
Scuola di Dottorato THE WAVE EQUATION Lucio Demeio Dipartimento di Ingegneria Industriale e delle Scienze Matematiche Lucio Demeio - DIISM wave equation 1 / 44 1 The Vibrating String Equation 2 Second
More informationAMS 529: Finite Element Methods: Fundamentals, Applications, and New Trends
AMS 529: Finite Element Methods: Fundamentals, Applications, and New Trends Lecture 25: Introduction to Discontinuous Galerkin Methods Xiangmin Jiao SUNY Stony Brook Xiangmin Jiao Finite Element Methods
More informationThe CG1-DG2 method for conservation laws
for conservation laws Melanie Bittl 1, Dmitri Kuzmin 1, Roland Becker 2 MoST 2014, Germany 1 Dortmund University of Technology, Germany, 2 University of Pau, France CG1-DG2 Method - Motivation hp-adaptivity
More information