The solution of time harmonic wave equations using complete families of elementary solutions
|
|
- Duane Howard
- 6 years ago
- Views:
Transcription
1 The solution of time harmonic wave equations using complete families of elementary solutions Peter Monk Department of Mathematical Sciences University of Delaware Newark, DE USA Joint research with Tomi Huttunen and Teemu Luostari, Kuopio, Finland. Research supported in part by a grant from AFOSR
2 Linear time harmonic wave equations Goal: find numerical solutions of common time harmonic wave equations. Examples include: Helmholtz Equation: u + κ 2 n 2 u = 0.
3 Linear time harmonic wave equations Goal: find numerical solutions of common time harmonic wave equations. Examples include: Helmholtz Equation: Maxwell s equations: u + κ 2 n 2 u = 0. µ 1 u κ 2 ɛu = 0
4 Linear time harmonic wave equations Goal: find numerical solutions of common time harmonic wave equations. Examples include: Helmholtz Equation: Maxwell s equations: Navier s equation: u + κ 2 n 2 u = 0. µ 1 u κ 2 ɛu = 0 (λ + 2µ) ( u) µ ( u) + κ 2 ρu = 0 together with appropriate boundary conditions.
5 Examples of solutions Suppose the coefficients are constant: Helmholtz Equation: u + κ 2 n 2 u = 0, u = exp(iκnx d), d = 1.
6 Examples of solutions Suppose the coefficients are constant: Helmholtz Equation: u + κ 2 n 2 u = 0, u = exp(iκnx d), d = 1. Maxwell s equations: µ 1 u κ 2 ɛu = 0, u = p exp(iκ µɛx d), d = 1, p d.
7 Examples of solutions Suppose the coefficients are constant: Helmholtz Equation: u + κ 2 n 2 u = 0, u = exp(iκnx d), d = 1. Maxwell s equations: µ 1 u κ 2 ɛu = 0, u = p exp(iκ µɛx d), d = 1, p d. Navier s equation: (λ + 2µ) ( u) µ ( u) + κ 2 ρu = 0 u = p exp(iκ s x d) + q exp(iκ p x d) with p d, q d, κ p = ω ρ/(λ + 2µ), κ s = ω ρ/µ
8 Helmholtz equation Example: find numerical solutions of the Helmholtz equation: u + κ 2 n 2 u = 0.
9 Helmholtz Examples [Huttunen] 20 θ = x (mm) y (mm)
10 Model Acoustic Problem Given a bounded domain Ω R N, N = 2, 3, find u such that u + κ 2 n 2 u = 0 in Ω u iκu n = g on Γ := Ω. Assume the wave number κ is real. Important parameter: the wavelength (of a plane wave solution) is 2π/(κn), so large κ corresponds to short wavelength.
11 Standard Finite Elements Weak form: find u h S h H 1 (Ω) such that u h ξ κ 2 n 2 u h ξ dv iκ u h ξ da = Ω for all ξ S h. This gives rise to a matrix problem (S κ 2 M iκb) u = g Γ Ω gξ da The matrix is complex symmetric (not Hermitian) with eigenvalues having positive and negative real parts if k is large so The linear system is hard to solve efficiently. and the matrix has to be large...
12 High frequency: many DoF The discretization must become finer at high wave number κ Suppose n is constant. Melenk and Sauter 1 show, under several non-restrictive assumptions, that for a conforming hp finite element approximation of our model boundary value problem, a quasi-optimal error estimate holds with constant independent of κ if h p = δ and p = 1 + C ln(κ) nκ 1 Vienna University: ASC Report No. 31/2009, Wave-Number Explicit Convergence Analysis for Galerkin Discretizations of the Helmholtz Equation, J.M. Melenk and S. Sauter
13 High frequency: many DoF The discretization must become finer at high wave number κ Suppose n is constant. Melenk and Sauter 1 show, under several non-restrictive assumptions, that for a conforming hp finite element approximation of our model boundary value problem, a quasi-optimal error estimate holds with constant independent of κ if Note h p = δ and p = 1 + C ln(κ) nκ #DoF per wavelength = (2π/(κn))(p/h) = 2π/δ Best case: we need at least O(κ N ) degrees of freedom. 1 Vienna University: ASC Report No. 31/2009, Wave-Number Explicit Convergence Analysis for Galerkin Discretizations of the Helmholtz Equation, J.M. Melenk and S. Sauter
14 Other methods We will now look at a DG method using non-standard shape functions called the Ultra Weak Variational Formulation 2 (UWVF) For an h version analysis of plane wave DG methods in general (but not UWVF) see [ 3 ] below. For other methods using special basis functions see for example Trefftz methods, PUFEM 4, Discontinuous Enrichment 5, Integral equations 6,... 2 O. Cessenat and B. Després, Application of the ultra-weak variational formulation of elliptic PDEs to the 2-dimensional Helmholtz problem. SIAM J. Numer. Anal., 35:255 99, C. Gittelson, R. Hiptmair, and I. Perugia. Plane wave discontinuous Galerkin methods. ESAIM: Mathematical Modeling and Numerical Analysis, 43: , J.M. Melenk. On generalized finite element methods. PhD thesis, University of Maryland, College Park, MD, M. Amara, R. Djellouli, and C. Farhat. Convergence analysis of a discontinuous Galerkin method with plane waves and Lagrange multipliers for the solution of Helmholtz problems. SIAM J. Numer. Anal., 47: , F. Ecevit and F. Reitich. Analysis of multiple scattering iterations for high-frequency scattering problems. I: the two-dimensional case. Numer. Math., 114: , 2009.
15 Time harmonic symmetric hyperbolic system 7 Given a domain Ω R N, N = 2, 3, and temporal frequency ω, we seek u = u(x) C m (so U(x, t) = R(u(x) exp( iωt))) such that N iωau + A j u xj + Bu = 0 in Ω j=1 where all the matrices are real symmetric. Also A(x) is positive definite, B(x) positive semidefinite and we assume that the A j are constant. Boundary condition: (D N)u = g on Γ := Ω where D = N j=1 A jn j and n is the unit outward normal to Ω. N is chosen to enforce the desired boundary condition. 7 G.Gabard. Discontinuous Galerkin methods with plane waves for time-harmonic problems. J. Comput. Phys., 225(2): , 2007.
16 Example: Helmholtz s equation in R 2 We seek to approximate the function p that satisfies p + (ω/c) 2 n(x)p = 0 in Ω and ( ) c p c iω n + σp = Q p iω n σp + g on Γ Introduce v such that i(ω/c)v = p. Define k = ω/c. Then if u = (p, v 1, v 2 ) T, we see that u satisfies the hyperbolic system with B = 0 and A 1 = A = n(x)/c /c /c, A 2 =, ,
17 Boundary conditions for the 2D example The boundary condition matrix D = n 1 A 1 + n 2 A 2 0 n 1 n 2 D = n 1 0 0, n Let L ± = 1 2σ (±σ, n 1, n 2 ). Note that L ± u = (±σp n 1 v 1 n 2 v 2 )/ 2σ corresponds to an impedance boundary condition. We set D + = (L + ) T L + and D = (L ) T L so that D = D + + D On the boundary Γ we can choose N = D 2(D + QL (L + ) T ), Q 1.
18 The Mesh Approximate the domain Ω by a tetrahedral finite element mesh consisting of N h tetrahedra K m, m = 1,, N h of maximum diameter h. N iωau + A j u xj + Bu = 0 j=1 in Ω
19 Derivation of the UWVF Let K denote an element in the mesh with faces f K and outward normal n K. Let u K = u K. Multiply the equation by a smooth function φ K on K and integrate by parts: iω + K f K Au K φ K dv f K N j=1 D K u K φ K da = 0 ( A j u φ K x j + Bu φ K ) dv Here the overbar represents complex conjugation, and we use D K to denote D on K.
20 Derivation continued Splitting the flux D K = D K,+ + D K, and using the boundary condition and symmetry of the matrices we get N u K iωaφ K A j φ K x j + Bφ K dv K + K f K K Γ K j=1 ( (D K,+ u K ) φ K da + f ((D K + N)u K + g) φ K da = 0. f ) (D K, u K ) φ K da
21 Non standard basis Given Y K X K (L 2 ( K )) m we choose φ K to be a solution of the adjoint equation with impedance boundary condition iωaφ K N A j φ K x j + Bφ K = 0 on K j=1 We obtain (D K,+ u K ) φ K da + (D K, u K ) φ K da K K K Γ f K K ((D K + N)u K + g) φ K da = 0. f
22 We can write D K,± = ±L K,± (L K,± ) T and N = D K 2(D K, + QL K, (L K,+ ) T ). Using also L K, = L K,+ we get: (L K,+ ) T u K (L K,+ ) T φ K da K + ((L K,+ ) T u K ) (L K, ) T φ K da K f K K f Q((L K,+ ) T u K + g) (L K, ) T φ K da = 0. Γ K Set Y K = (L K,+ ) T φ K. Define the unknown impedance flux and operator F K : X K X K by X K = (L K,+ ) T u K and F K (Y K ) = (L K, ) T φ K.
23 The UWVF Summarizing, we see that X K X K (L 2 ( K )) m satisfies X K Y K da + X K F K (Y K ) da K K f K K f Q(X K + g) F K (Y K ) da = 0. K Γ for all Y K X K (L 2 ( K )) m and all tetraheda K in the mesh. This is the UWVF of Cessenat and Després before discretization Note that for the acoustic case m = 1, X K = L 2 ( K ), and for electromagnetism m = 3, while X K = L 2 t ( K ). Cessenat and Després prove existence of a unique solution for acoustics and electromagnetism.
24 Useful result Lemma (Isometry Lemma) For all Y K X K, F K (Y K ) L 2 ( K ) YK L 2 ( K ) with equality if B = 0. To see this: let Y K = (L K,+ ) T φ K. Then using the adjoint equation and integration by parts D K φ K φ K da = 2 φ K Bφ K dv. K Splitting D K = D K,+ + D K, and using the definition of F K completes the demonstration. K
25 The discrete UWVF Let Vp K denote a finite dimensional space of p := p(k ) solutions of the adjoint problem on an element K. Let { } Xp K = v X K : v = (L K,+ ) T φ K for some φ K Vp K The discrete problem is to seek X K h X K p such that K X K h YK da + K Γ K f K K f X K h Q(X K h + g) F K (Y K ) da = 0 for all Y K X K p and all elements K. F K (Y K ) da
26 Model Acoustic Problem Given a bounded domain Ω R N, N = 2, 3, find u such that Assume κ is real. u + κ 2 n 2 u = 0 in Ω u iκu n = g on Γ := Ω.
27 2D acoustic plane wave basis Let n K = n K assumed piecewise constant. Using p directions d K m = (cos(θ m ), sin(θ m )) T, θ m = 2πm/p we set { } = span exp(ikn K x d K m), 1 m p V K p A function Y K X K p can then be written Y K = 1 ( ) 1 p + I ym K exp(ikn K x d K 2 ik ν m). K m=1 Then F K (Y K ) is easy to compute since F K (Y K ) = 1 ( ) 1 p I ym K exp(ikn K x d K 2σ ik ν m). K m=1 X K h is expanded in the same way as YK.
28 Some comments on plane wave bases On straight tets/triangles, the UWVF integrals can be computed in a closed form, therefore assembly is rapid and we avoid the need for complicated quadratures (c.f. PUFEM). Ill-conditioning may occur when elements are small or at low frequencies since exp(ikx d) = 1 + ikx d +... Fourier-Bessel functions offer an attractive alternative.
29 Properties of the acoustic UWVF [Cessenat/Després, 2D] If p = 2µ + 1 X X h L 2 (Γ) Chµ 1/2 u C µ+1 (Ω) [Buffa-Monk, 2D] 8 Under the above conditions (with u h reconstructed from X h ), and if Ω is convex u u h L 2 (Ω) C(κ)hµ 1 u C µ+1 (Ω) 8 A. Buffa and P. Monk. Error estimates for the Ultra Weak Variational Formulation of the Helmholtz equation. ESAIM: Mathematical Modeling and Numerical Analysis, 42:925 40, P. Monk and D.Q. Wang. A least squares method for the Helmholtz equation. Comput. Meth. Appl. Mech. Eng., 175:121 36, 1999
30 Properties of the acoustic UWVF [Cessenat/Després, 2D] If p = 2µ + 1 X X h L 2 (Γ) Chµ 1/2 u C µ+1 (Ω) [Buffa-Monk, 2D] 8 Under the above conditions (with u h reconstructed from X h ), and if Ω is convex u u h L 2 (Ω) C(κ)hµ 1 u C µ+1 (Ω) Proof: We relate the UWVF back to an upwind DG method. Then usual DG results gives an estimate in DG norms. Using duality we can then estimate the L 2 norm, 9 loosing a half power of h along the way. 8 A. Buffa and P. Monk. Error estimates for the Ultra Weak Variational Formulation of the Helmholtz equation. ESAIM: Mathematical Modeling and Numerical Analysis, 42:925 40, P. Monk and D.Q. Wang. A least squares method for the Helmholtz equation. Comput. Meth. Appl. Mech. Eng., 175:121 36, 1999
31 Properties of the acoustic UWVF [Cessenat/Després, 2D] If p = 2µ + 1 X X h L 2 (Γ) Chµ 1/2 u C µ+1 (Ω) [Buffa-Monk, 2D] 8 Under the above conditions (with u h reconstructed from X h ), and if Ω is convex u u h L 2 (Ω) C(κ)hµ 1 u C µ+1 (Ω) Proof: We relate the UWVF back to an upwind DG method. Then usual DG results gives an estimate in DG norms. Using duality we can then estimate the L 2 norm, 9 loosing a half power of h along the way. Problem: we we usually fix a grid and use a p version. 8 A. Buffa and P. Monk. Error estimates for the Ultra Weak Variational Formulation of the Helmholtz equation. ESAIM: Mathematical Modeling and Numerical Analysis, 42:925 40, P. Monk and D.Q. Wang. A least squares method for the Helmholtz equation. Comput. Meth. Appl. Mech. Eng., 175:121 36, 1999
32 Much improved analysis [Hiptmair, Moiola, Perugia, 2009] 10 p-version error estimate (using κ-dependent norms): ( ) log(p) r 1/2 κ u u h L 2 (Ω) Chr 1 u r+1,κ,ω p 10 R. Hiptmair, A. Moiola, and I. Perugia, Plane wave discontinuous Galerkin methods for the 2D Helmholtz equation: analysis of the p-version
33 Much improved analysis [Hiptmair, Moiola, Perugia, 2009] 10 p-version error estimate (using κ-dependent norms): ( ) log(p) r 1/2 κ u u h L 2 (Ω) Chr 1 u r+1,κ,ω p [Cessenat, Despres] The discrete problem has the form (D C)x = b where D is Hermitian positive definite and the eigenvalues of D 1 C lie in the closure of the unit disk excluding 1. The condition number of D satisfies cond(d) C h 2p 2 10 R. Hiptmair, A. Moiola, and I. Perugia, Plane wave discontinuous Galerkin methods for the 2D Helmholtz equation: analysis of the p-version
34 Numerical results: 2D mesh refinement 11 We take u(x) = i 4 H(1) 0 (κ x x 0 ) with x 0 is outside the computational domain h = 0.5 h = T. Huttunen, P. Monk, and J.P. Kaipio. Computational aspects of the Ultra Weak Variational Formulation. J. Comput. Phys., 182:27 46, 2002.
35 Results for κ = 20, p = 15 (µ = 7) Mesh size h L 2 (Ω) Error (%) Order cond(d) Order Measured global convergence O(h 7.8 ). Cessenat and Després predict that the condition number of D will increase O(h 12 )
36 Results for κ = 40, p = 21 (µ = 10). Mesh size h L 2 (Ω) Error (%) Order cond(d) Order Measured global convergence O(h 10 ). Cessenat and Després predict that the condition number of D increases O(h 18 ).
37 Typical Electromagnetic Application [Tomi Huttunen] Example: Simulate microwave interaction with wood (part of a non-destructive testing project). A transmitting and receiving antenna are shown.
38 Maxwell s equations: Cavity problem E, H: Electric/magnetic field k : Wave-number ɛ r : Relative permittivity µ r : Relative permeability Ω: Bounded domain ikɛ r E H = 0 in Ω ikµ r H + E = 0 in Ω H n ηe T = Q(H n + ηe T ) 2ηg on Γ = Ω where η > 0 and Q 1.
39 Required continuity between elements Let E k = E Ωk and E j = E Ωj (similarly for H) then if E, H is to be a solution of the Maxwell s equations we need E k n K k = E j n K j and H k n K k = H j n K j on Σ j,k In the Ultra Weak Variational Formulation (UWVF) this is achieved by demanding that impedance data agree on the interfaces, so on Σ j,k E k n K k + η(h k n K k ) n K k = E j n K j + η(h j n K j ) n K j E k n K k η(h k n K k ) n K k = E j n K j η(h j n K j ) n K j where η > 0 is a parameter (function) on Σ j,k (e.g. η = R( µ r /ɛ r )).
40 12 13 Maxwell UWVF Let X k = E k n K k η(h k ) T Ωk and Y k = ξ k n K k η(ψ k ) T Ωk and let F k (Y k ) = ξ k n K k η(ψ k ) T Ωk then, for a tetrahedron surrounded by four other tetrahedra 1 Ω k η X ky k ds = 1 j Σ k,j η X jf k (Y k ) ds. Boundary faces are handled using the boundary condition in place of X j. 12 O. Cessenat. Application d une nouvelle formulation variationnelle aux équations d ondes harmoniques. Problèmes de Helmholtz 2D et de Maxwell 3D. PhD thesis, Université Paris IX Dauphine, T. Huttunen, M. Malinen, and P.B. Monk. Solving Maxwell s equations using the Ultra Weak Variational Formulation. J. Comput. Phys., 223:731 58, 2007.
41 The Discrete UWVF: plane waves For each element Ω k we choose p k directions d j, j = 1,..., p k on the unit sphere [Sloan et al.] and for each direction two polarizations q (1,2) j then define (choosing κ = ɛ r µ r ) X h k = 2 p k s=1 j=1 ( x j,s q s j n K k ηκ ) (d j q s j kµ ) T r exp(iκd j x) The test function is, for 1 r p k, 1 t 2 ( Yk h = q t r nk k ηκ ) (d r q t r kµ ) T exp(iκd r x) r Ωk In this case F k (Yk h ) is easy to compute: ( F k (Yk h ) = q t r nk k ηκ ) (d r q t r kµ ) T exp(iκd r x) r Ωk Ωk
42 Results for wood at 5GHz Using a Perfectly Matched Layer, and keeping a fixed mesh but varying the number of directions with the local wavelength. E x E y E x E y z (m) z (m) z (m) z (m) y (m) 70 db y (m) 70 db y (m) 70 db y (m) 70 db E z 0 E 0 E z 0 E 0 z (m) z (m) z (m) z (m) y (m) 70 db y (m) 70 db y (m) 70 db y (m) 70 db ɛ r = 3 ɛ r = 8
43 Scattering from a sphere Sphere of radius 0.25 inside cube [.5,.5] 3, κ = 100, ɛ = µ = 1 (λ = 0.06). PML width 0.1 (uses 3,474,770 degrees of freedom). Real E 1, UWVF Real E 2, UWVF Real E 2, UWVF
44 Iterative solution The UWVF linear system can be solved by simple iterative scheme. We use BiCGStab Bi CGStab Relative residual Iteration BiCGStab convergence for a problem having 3,474,770 degrees of freedom using a 24 processor cluster (2.8GHz P4, 48Gb memory total, 1000BaseT). Solution time is 451s using 25.3 GB memory (109 iterations).
45 Elastic UWVF: Navier s equation The elasticity problem: Navier s equation with boundary conditions µ u + (λ + µ) ( u) + ω 2 ρu = 0 T (n) (u) iσu = Q( T (n) (u) iσu) + g in Ω on Γ where the traction operator T (n) (u) is defined by T (n) (u) = 2µ u + λn u + µn u. n where n is an outward unit normal, g is the source term, ω is the angular frequency of the field, u is the time-harmonic displacement vector, λ and µ are the Lamé constants and ρ is the density of the medium.
46 The Elastic UWVF 14 [T. Luostari] We can derive the following UWVF of finding X k such that k k σ 1 X k ( T (n k ) (e k ) iσe k ) Ω k k Γ k Qσ 1 X k (T (n k ) (e k ) iσe k ) = k j k,j σ 1 X j (T (n k ) (e k ) iσe k ) Γ k σ 1 g (T (n k ) (e k ) iσe k ) (1) for all plane wave (S- or P-wave) test functions e k. Here X k = T (n k ) (u k ) iσu k on Ω k. So far, in 3D numerical simulations we use an ad hoc choice for coupling parameter (flux parameter) where I is the unit matrix. σ = ωρr{c P }I (2) 14 2D: T. Huttunen, P. Monk, F. Collino, and J.P. Kaipio. The Ultra Weak Variational Formulation for elastic wave problems. SIAM J. Sci. Comput., 25: , 2004.
47 Plane wave propagation in a unit cube The exact solution is of the form u = A 1 d exp(iκ P x d) + A 2 d SH exp(iκ S x d) + A 3 d SV exp(iκ S x d) where the wave numbers are κ P = ω/c P, κ S = ω/c S, the direction d [ ], d SH = d, d SV = d d and the amplitudes A 1 = A 2 = A 3 = 1. In addition, u P = 0 and u SH = u SV = 0. As a boundary condition we choose Q = 0.
48 Results for p-convergence 10 1 κ P = and κ S = p P =0.25p S p P =(1/3)p S p P =0.5p S p P =(2/3)p S p P =p S Relative error (%) max(cond(d k )) Number of basis functions Number of basis functions Results when κ P = , κ SH = κ SV = with different ratios between p P /p S and mesh size is fixed (24 tetrahedra in the unit cube).
49 Key questions We are starting to investigate: 1 Convergence to evanescent modes (e.g. Rayleigh waves) 2 Choice of number of directions 3 Choice of the coupling operator σ
50 Other developments Hybridization with finite element methods (acoustics) 15 PML 16 and coupling to integral equations 17 Fluid-solid problem in 2D P. Monk, J. Schöberl, and A. Sinwel. Hybridizing Raviart-Thomas elements for the Helmholtz equation. Electromagnetics, 30:149 76, T. Huttunen, J.P. Kaipio, and P. Monk. The perfectly matched layer for the ultra weak variational formulation of the 3D Helmholtz equation. Int. J. Numer. Meth. Eng., 61: , E. Darrigrand and P. Monk. Coupling of the Ultra-Weak Variational Formulation and an integral representation using a Fast Multipole Method in electromagnetism. J. Comput. Appl. Math., 204:400 7, T. Huttunen, J. Kaipio, and P. Monk. An ultra-weak method for acoustic fluid-solid interaction. J. Comput. Appl. Math., 213:166 85, 2008.
51 Conclusions The method is applicable to a variety of wave propagation problems using a common framework Excellent accuracy can be obtained if the complete family is well matched to the problem Robustness is an issue (particularly ill-conditioning)
Solving 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 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 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 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 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 informationTHE USE OF PLANE WAVES TO APPROXIMATE WAVE PROPAGATION IN ANISOTROPIC MEDIA * 1. Introduction
Journal of Computational Mathematics, Vol. 25, No. 3, 27, 35 367. THE USE OF PLANE WAVES TO APPROXIMATE WAVE PROPAGATION IN ANISOTROPIC MEDIA * Tomi Huttunen (Department of Physics, University of Kuopio,
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 informationA 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 Rome, 10th Apr 2017 Joint work with: Emmanuil
More informationULTRA-WEAK VARIATIONAL FORMULATION AND EFFICIENT INTEGRAL REPRESENTATION IN ELECTROMAGNETISM: A THOROUGH STUDY OF THE ALGORITHM COMPLEXITY
ESAIM: PROCEEDINGS, October 2007, Vol. 22, 190-197 Gabriel Caloz & Monique Dauge, Editors ULTRA-WEAK VARIATIONAL FORMULATION AND EFFICIENT INTEGRAL REPRESENTATION IN ELECTROMAGNETISM: A THOROUGH STUDY
More informationCombining the Ultra-Weak Variational Formulation and the Multilevel Fast Multipole Method
Combining the Ultra-Weak Variational Formulation and the Multilevel Fast Multipole Method Eric Darrigrand, Peter Monk To cite this version: Eric Darrigrand, Peter Monk. Combining the Ultra-Weak Variational
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 informationA 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 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 informationUltra-Weak Variational Formulation and Integral Representation using a Fast Multipole Method for the Equations of Electromagnetism
Ultra-Weak Variational Formulation and Integral Representation using a Fast Multipole Method for the Equations of Electromagnetism Eric Darrigrand and Peter Monk Université de Rennes 1 University of Delaware
More informationA Multigrid Method for Two Dimensional Maxwell Interface Problems
A Multigrid Method for Two Dimensional Maxwell Interface Problems Susanne C. Brenner Department of Mathematics and Center for Computation & Technology Louisiana State University USA JSA 2013 Outline A
More informationHybrid (DG) Methods for the Helmholtz Equation
Hybrid (DG) Methods for the Helmholtz Equation Joachim Schöberl Computational Mathematics in Engineering Institute for Analysis and Scientific Computing Vienna University of Technology Contributions by
More informationCoercivity of high-frequency scattering problems
Coercivity of high-frequency scattering problems Valery Smyshlyaev Department of Mathematics, University College London Joint work with: Euan Spence (Bath), Ilia Kamotski (UCL); Comm Pure Appl Math 2015.
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 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 informationAn inverse problem for eddy current equations
An inverse problem for eddy current equations Alberto Valli Department of Mathematics, University of Trento, Italy In collaboration with: Ana Alonso Rodríguez, University of Trento, Italy Jessika Camaño,
More informationANALYSIS OF A FINITE ELEMENT PML APPROXIMATION FOR THE THREE DIMENSIONAL TIME-HARMONIC MAXWELL PROBLEM
MATHEMATICS OF COMPUTATION Volume 77, Number 261, January 2008, Pages 1 10 S 0025-5718(07)02037-6 Article electronically published on September 18, 2007 ANALYSIS OF A FINITE ELEMENT PML APPROXIMATION FOR
More informationOptimized Schwarz Methods for Maxwell Equations with Discontinuous Coefficients
Optimized Schwarz Methods for Maxwell Equations with Discontinuous Coefficients Victorita Dolean, Martin Gander, Erwin Veneros To cite this version: Victorita Dolean, Martin Gander, Erwin Veneros. Optimized
More informationA Stable Spectral Difference Method for Triangles
A Stable Spectral Difference Method for Triangles Aravind Balan 1, Georg May 1, and Joachim Schöberl 2 1 AICES Graduate School, RWTH Aachen, Germany 2 Institute for Analysis and Scientific Computing, Vienna
More informationOptimization and Control Methods for Ultrasound Surgery
Procedings of COBEM 25 Copyright c 25 by ABCM November 6-, 25, Ouro Preto, MG Optimization and Control Methods for Ultrasound Surgery Matti Malinen Department of Applied Physics, University of Kuopio,
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 informationPREPRINT 2010:23. A nonconforming rotated Q 1 approximation on tetrahedra PETER HANSBO
PREPRINT 2010:23 A nonconforming rotated Q 1 approximation on tetrahedra PETER HANSBO Department of Mathematical Sciences Division of Mathematics CHALMERS UNIVERSITY OF TECHNOLOGY UNIVERSITY OF GOTHENBURG
More informationCONTINUED-FRACTION ABSORBING BOUNDARY CONDITIONS FOR THE WAVE EQUATION
Journal of Computational Acoustics, Vol. 8, No. 1 (2) 139 156 c IMACS CONTINUED-FRACTION ABSORBING BOUNDARY CONDITIONS FOR THE WAVE EQUATION MURTHY N. GUDDATI Department of Civil Engineering, North Carolina
More informationarxiv: v1 [math.na] 19 Nov 2018
A Bivariate Spline Solution to the Exterior Helmholtz Equation and its Applications Shelvean Kapita Ming-Jun Lai arxiv:1811.07833v1 [math.na] 19 Nov 2018 November 20, 2018 Abstract We explain how to use
More informationBETI for acoustic and electromagnetic scattering
BETI for acoustic and electromagnetic scattering O. Steinbach, M. Windisch Institut für Numerische Mathematik Technische Universität Graz Oberwolfach 18. Februar 2010 FWF-Project: Data-sparse Boundary
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 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 informationDG discretization of optimized Schwarz methods for Maxwell s equations
DG discretization of optimized Schwarz methods for Maxwell s equations Mohamed El Bouajaji, Victorita Dolean,3, Martin J. Gander 3, Stéphane Lanteri, and Ronan Perrussel 4 Introduction In the last decades,
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 WAVENUMBER INDEPENDENT BOUNDARY ELEMENT METHOD FOR AN ACOUSTIC SCATTERING PROBLEM
A WAVENUMBER INDEPENDENT BOUNDARY ELEMENT METHOD FOR AN ACOUSTIC SCATTERING PROBLEM S LANGDON AND S N CHANDLER-WILDE Abstract. In this paper we consider the impedance boundary value problem for the Helmholtz
More informationMéthodes de Galerkine Discontinues pour l imagerie sismique en domaine fréquentiel
Méthodes de Galerkine Discontinues pour l imagerie sismique en domaine fréquentiel Julien Diaz, Équipe Inria Magique 3D Séminaire de Mathématiques Appliquées, Laboratoire de Mathématiques Nicolas Oresme
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 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 informationHigh Frequency Scattering by Convex Polygons Stephen Langdon
Bath, October 28th 2005 1 High Frequency Scattering by Convex Polygons Stephen Langdon University of Reading, UK Joint work with: Simon Chandler-Wilde Steve Arden Funded by: Leverhulme Trust University
More informationTrefftz type method for 2D problems of electromagnetic scattering from inhomogeneous bodies.
Trefftz type method for 2D problems of electromagnetic scattering from inhomogeneous bodies. S. Yu. Reutsiy Magnetohydrodynamic Laboratory, P. O. Box 136, Mosovsi av.,199, 61037, Kharov, Uraine. e mail:
More informationAn adaptive fast multipole boundary element method for the Helmholtz equation
An adaptive fast multipole boundary element method for the Helmholtz equation Vincenzo Mallardo 1, Claudio Alessandri 1, Ferri M.H. Aliabadi 2 1 Department of Architecture, University of Ferrara, Italy
More informationClassical Scattering
Classical Scattering Daniele Colosi Mathematical Physics Seminar Daniele Colosi (IMATE) Classical Scattering 27.03.09 1 / 38 Contents 1 Generalities 2 Classical particle scattering Scattering cross sections
More informationNumerical approximation of output functionals for Maxwell equations
Numerical approximation of output functionals for Maxwell equations Ferenc Izsák ELTE, Budapest University of Twente, Enschede 11 September 2004 MAXWELL EQUATIONS Assumption: electric field ( electromagnetic
More informationThe Imaging of Anisotropic Media in Inverse Electromagnetic Scattering
The Imaging of Anisotropic Media in Inverse Electromagnetic Scattering Fioralba Cakoni Department of Mathematical Sciences University of Delaware Newark, DE 19716, USA email: cakoni@math.udel.edu Research
More informationAsymptotic Behavior of Waves in a Nonuniform Medium
Available at http://pvamuedu/aam Appl Appl Math ISSN: 1932-9466 Vol 12, Issue 1 June 217, pp 217 229 Applications Applied Mathematics: An International Journal AAM Asymptotic Behavior of Waves in a Nonuniform
More informationLecture Note III: Least-Squares Method
Lecture Note III: Least-Squares Method Zhiqiang Cai October 4, 004 In this chapter, we shall present least-squares methods for second-order scalar partial differential equations, elastic equations of solids,
More informationAdditive Schwarz method for scattering problems using the PML method at interfaces
Additive Schwarz method for scattering problems using the PML method at interfaces Achim Schädle 1 and Lin Zschiedrich 2 1 Zuse Institute, Takustr. 7, 14195 Berlin, Germany schaedle@zib.de 2 Zuse Institute,
More informationA SEAMLESS REDUCED BASIS ELEMENT METHOD FOR 2D MAXWELL S PROBLEM: AN INTRODUCTION
A SEAMLESS REDUCED BASIS ELEMENT METHOD FOR 2D MAXWELL S PROBLEM: AN INTRODUCTION YANLAI CHEN, JAN S. HESTHAVEN, AND YVON MADAY Abstract. We present a reduced basis element method (RBEM) for the time-harmonic
More informationTwo-level Domain decomposition preconditioning for the high-frequency time-harmonic Maxwell equations
Two-level Domain decomposition preconditioning for the high-frequency time-harmonic Maxwell equations Marcella Bonazzoli 2, Victorita Dolean 1,4, Ivan G. Graham 3, Euan A. Spence 3, Pierre-Henri Tournier
More informationScattering of electromagnetic waves by thin high contrast dielectrics II: asymptotics of the electric field and a method for inversion.
Scattering of electromagnetic waves by thin high contrast dielectrics II: asymptotics of the electric field and a method for inversion. David M. Ambrose Jay Gopalakrishnan Shari Moskow Scott Rome June
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 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 informationFull Wave Analysis of RF Signal Attenuation in a Lossy Rough Surface Cave Using a High Order Time Domain Vector Finite Element Method
Progress In Electromagnetics Research Symposium 2006, Cambridge, USA, March 26-29 425 Full Wave Analysis of RF Signal Attenuation in a Lossy Rough Surface Cave Using a High Order Time Domain Vector Finite
More informationA coupled BEM and FEM for the interior transmission problem
A coupled BEM and FEM for the interior transmission problem George C. Hsiao, Liwei Xu, Fengshan Liu, Jiguang Sun Abstract The interior transmission problem (ITP) is a boundary value problem arising in
More informationNedelec elements for computational electromagnetics
Nedelec elements for computational electromagnetics Per Jacobsson, June 5, 2007 Maxwell s equations E = jωµh () H = jωεe + J (2) (εe) = ρ (3) (µh) = 0 (4) J = jωρ (5) Only three of the equations are needed,
More informationThe Method of Fundamental Solutions applied to the numerical calculation of eigenfrequencies and eigenmodes for 3D simply connected domains
ECCOMAS Thematic Conference on Meshless Methods 2005 C34.1 The Method of Fundamental Solutions applied to the numerical calculation of eigenfrequencies and eigenmodes for 3D simply connected domains Carlos
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 informationA posteriori error estimates for the adaptivity technique for the Tikhonov functional and global convergence for a coefficient inverse problem
A posteriori error estimates for the adaptivity technique for the Tikhonov functional and global convergence for a coefficient inverse problem Larisa Beilina Michael V. Klibanov December 18, 29 Abstract
More informationComparison of a Finite Difference and a Mixed Finite Element Formulation of the Uniaxial Perfectly Matched Layer
Comparison of a Finite Difference and a Mixed Finite Element Formulation of the Uniaxial Perfectly Matched Layer V. A. Bokil a and M. W. Buksas b Center for Research in Scientific Computation a North Carolina
More informationRiesz bases of Floquet modes in semi-infinite periodic waveguides and implications
Riesz bases of Floquet modes in semi-infinite periodic waveguides and implications Thorsten Hohage joint work with Sofiane Soussi Institut für Numerische und Angewandte Mathematik Georg-August-Universität
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 informationAn eigenvalue method using multiple frequency data for inverse scattering problems
An eigenvalue method using multiple frequency data for inverse scattering problems Jiguang Sun Abstract Dirichlet and transmission eigenvalues have important applications in qualitative methods in inverse
More informationLocal Mesh Refinement with the PCD Method
Advances in Dynamical Systems and Applications ISSN 0973-5321, Volume 8, Number 1, pp. 125 136 (2013) http://campus.mst.edu/adsa Local Mesh Refinement with the PCD Method Ahmed Tahiri Université Med Premier
More informationBreast Cancer Detection by Scattering of Electromagnetic Waves
57 MSAS'2006 Breast Cancer Detection by Scattering of Electromagnetic Waves F. Seydou & T. Seppänen University of Oulu Department of Electrical and Information Engineering P.O. Box 3000, 9040 Oulu Finland
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 informationPDEs, part 1: Introduction and elliptic PDEs
PDEs, part 1: Introduction and elliptic PDEs Anna-Karin Tornberg Mathematical Models, Analysis and Simulation Fall semester, 2013 Partial di erential equations The solution depends on several variables,
More informationThe mathematics of scattering by unbounded, rough, inhomogeneous layers
The mathematics of scattering by unbounded, rough, inhomogeneous layers Simon N. Chandler-Wilde a Peter Monk b Martin Thomas a a Department of Mathematics, University of Reading, Whiteknights, PO Box 220
More informationEstimation of transmission eigenvalues and the index of refraction from Cauchy data
Estimation of transmission eigenvalues and the index of refraction from Cauchy data Jiguang Sun Abstract Recently the transmission eigenvalue problem has come to play an important role and received a lot
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 informationAMS subject classifications. Primary, 65N15, 65N30, 76D07; Secondary, 35B45, 35J50
A SIMPLE FINITE ELEMENT METHOD FOR THE STOKES EQUATIONS LIN MU AND XIU YE Abstract. The goal of this paper is to introduce a simple finite element method to solve the Stokes equations. This method is in
More informationGlobal Maxwellians over All Space and Their Relation to Conserved Quantites of Classical Kinetic Equations
Global Maxwellians over All Space and Their Relation to Conserved Quantites of Classical Kinetic Equations C. David Levermore Department of Mathematics and Institute for Physical Science and Technology
More information1 domain-based methods: convergence analysis for hp-fem. 2 discussion of some non-standard FEMs. 3 BEM for Helmholtz problems.
outline Helmholtz problems at large wavenumber k joint work with S. Sauter (Zürich) M. Löhndorf (Vienna) S. Langdon (Reading) 1 domain-based methods: convergence analysis for hp-fem 2 discussion of some
More informationarxiv: v1 [math.na] 11 Jul 2011
Multigrid Preconditioner for Nonconforming Discretization of Elliptic Problems with Jump Coefficients arxiv:07.260v [math.na] Jul 20 Blanca Ayuso De Dios, Michael Holst 2, Yunrong Zhu 2, and Ludmil Zikatanov
More informationADI iterations for. general elliptic problems. John Strain Mathematics Department UC Berkeley July 2013
ADI iterations for general elliptic problems John Strain Mathematics Department UC Berkeley July 2013 1 OVERVIEW Classical alternating direction implicit (ADI) iteration Essentially optimal in simple domains
More informationIntroduction to PML in time domain
Introduction to PML in time domain Alexander Thomann Introduction to PML in time domain - Alexander Thomann p.1 Overview 1 Introduction 2 PML in one dimension Classical absorbing layers One-dimensional
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 informationFirst-order overdetermined systems. for elliptic problems. John Strain Mathematics Department UC Berkeley July 2012
First-order overdetermined systems for elliptic problems John Strain Mathematics Department UC Berkeley July 2012 1 OVERVIEW Convert elliptic problems to first-order overdetermined form Control error via
More informationA Direct Method for reconstructing inclusions from Electrostatic Data
A Direct Method for reconstructing inclusions from Electrostatic Data Isaac Harris Texas A&M University, Department of Mathematics College Station, Texas 77843-3368 iharris@math.tamu.edu Joint work with:
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 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 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 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 informationc 1998 Society for Industrial and Applied Mathematics
SIAM J. NUMER. ANAL. Vol. 36, No. 1, pp. 251 274 c 1998 Society for Industrial and Applied Mathematics THE ADAPTIVE COMPUTATION OF FAR-FIELD PATTERNS BY A POSTERIORI ERROR ESTIMATION OF LINEAR FUNCTIONALS
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 informationWeierstrass Institute for Applied Analysis and Stochastics Direct and Inverse Elastic Scattering Problems for Diffraction Gratings
Weierstrass Institute for Applied Analysis and Stochastics Direct and Inverse Elastic Scattering Problems for Diffraction Gratings Johannes Elschner & Guanghui Hu Mohrenstrasse 39 10117 Berlin Germany
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 informationChapter 6. Finite Element Method. Literature: (tiny selection from an enormous number of publications)
Chapter 6 Finite Element Method Literature: (tiny selection from an enormous number of publications) K.J. Bathe, Finite Element procedures, 2nd edition, Pearson 2014 (1043 pages, comprehensive). Available
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 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 informationNew constructions of domain decomposition methods for systems of PDEs
New constructions of domain decomposition methods for systems of PDEs Nouvelles constructions de méthodes de décomposition de domaine pour des systèmes d équations aux dérivées partielles V. Dolean?? F.
More informationFinite volume method for nonlinear transmission problems
Finite volume method for nonlinear transmission problems Franck Boyer, Florence Hubert To cite this version: Franck Boyer, Florence Hubert. Finite volume method for nonlinear transmission problems. Proceedings
More informationINTERGRID OPERATORS FOR THE CELL CENTERED FINITE DIFFERENCE MULTIGRID ALGORITHM ON RECTANGULAR GRIDS. 1. Introduction
Trends in Mathematics Information Center for Mathematical Sciences Volume 9 Number 2 December 2006 Pages 0 INTERGRID OPERATORS FOR THE CELL CENTERED FINITE DIFFERENCE MULTIGRID ALGORITHM ON RECTANGULAR
More informationDG Methods for Aerodynamic Flows: Higher Order, Error Estimation and Adaptive Mesh Refinement
HONOM 2011 in Trento DG Methods for Aerodynamic Flows: Higher Order, Error Estimation and Adaptive Mesh Refinement Institute of Aerodynamics and Flow Technology DLR Braunschweig 11. April 2011 1 / 35 Research
More informationNitsche-type Mortaring for Maxwell s Equations
Nitsche-type Mortaring for Maxwell s Equations Joachim Schöberl Karl Hollaus, Daniel Feldengut Computational Mathematics in Engineering at the Institute for Analysis and Scientific Computing Vienna University
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 informationPARTITION OF UNITY FOR THE STOKES PROBLEM ON NONMATCHING GRIDS
PARTITION OF UNITY FOR THE STOES PROBLEM ON NONMATCHING GRIDS CONSTANTIN BACUTA AND JINCHAO XU Abstract. We consider the Stokes Problem on a plane polygonal domain Ω R 2. We propose a finite element method
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 informationElectromagnetism. Christopher R Prior. ASTeC Intense Beams Group Rutherford Appleton Laboratory
lectromagnetism Christopher R Prior Fellow and Tutor in Mathematics Trinity College, Oxford ASTeC Intense Beams Group Rutherford Appleton Laboratory Contents Review of Maxwell s equations and Lorentz Force
More informationA P4 BUBBLE ENRICHED P3 DIVERGENCE-FREE FINITE ELEMENT ON TRIANGULAR GRIDS
A P4 BUBBLE ENRICHED P3 DIVERGENCE-FREE FINITE ELEMENT ON TRIANGULAR GRIDS SHANGYOU ZHANG DEDICATED TO PROFESSOR PETER MONK ON THE OCCASION OF HIS 6TH BIRTHDAY Abstract. On triangular grids, the continuous
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 information