High Frequency Scattering by Convex Polygons Stephen Langdon
|
|
- Jasper Dean
- 5 years ago
- Views:
Transcription
1 Bath, October 28th 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 of Bath, October 28th 2005
2 Bath, October 28th Contents Background - high frequency scattering problems Scattering by a convex polygon Understanding solution behaviour Approximating the solution efficiently Galerkin method Collocation method Numerical Results
3 Bath, October 28th Background - scattering problems
4 Bath, October 28th Scattering theory Concerned with the effect obstacles or inhomogeneities have on an incident wave. Direct problem: Compute scattered field when incident field and scattering obstacle are known. Inverse problem: Compute shape / properties of scatterer when scattered field is known.
5 Bath, October 28th Applications Prospecting for hydrocarbons Searching for archaeological remains Scanning for tumours Modelling road / rail / aircraft noise Defence applications Identifying objects via radar / sonar
6 Bath, October 28th Mathematical formulation for time harmonic waves Propagation of sound in a homogeneous media is mathematically described by the wave equation ũ(x, t) 1 c 2 2 ũ(x, t) = 0, (1) t2 where c denotes the speed of sound, and ũ is the induced pressure. Considering only time harmonic waves of the form ũ(x, t) = u(x)e iωt, where ω denotes frequency, the wave equation (1) then reduces to the Helmholtz equation: u + k 2 u = 0, where k = ω/c is the wavenumber.
7 Bath, October 28th Analytical solution possible only for very simple obstacles. In general we need to use a numerical method. Finite Element Method Need to cover whole domain in a mesh Leads to large sparse systems Boundary Element Method Only need to cover boundary in a mesh Leads to smaller full systems For each approach - difficulties when k is large.
8 Bath, October 28th Problem: Find u such that FEM for a model 1D problem u + k 2 u = 0, 0 < x < 1 u(0) = 1 u (1) iku(1) = 0 Approximate by piecewise linear FEM on a uniform grid of size h. THEOREM (Ihlenburg and Babuška) Let u h denote the FE solution. Then for h sufficiently small This is VERY BAD when k is large. u u h H 1 (0,1) C(hk + h 2 k 3 ).
9 Bath, October 28th First difficulty - local errors To accurately model a wave potential the rule of thumb is 10 elements per wave length. This is expensive for high wavenumbers. Figure 1: Influence of the number of elements per wavelength on the accuracy of numerical wave modelling
10 Bath, October 28th Second difficulty - pollution errors Additional error caused by innacurate modelling of wavelength - can build up across numerical model. Figure 2: Pollution error in numerical wave modelling
11 Bath, October 28th Scattering by convex polygons u i, incident wave u + k 2 u = 0 obstacle u = 0 D Γ x 2 x 1
12 Bath, October 28th u i, incident wave u + k 2 u = 0 obstacle u = 0 D Γ x 2 Green s representation theorem: u(x) = u i (x) Φ(x, y) u (y)ds(y), x D, n Γ x 1 where Φ(x, y) := i 4 H(1) 0 (k x y ).
13 Bath, October 28th u i, incident wave u + k 2 u = 0 obstacle u = 0 D Γ x 2 x 1 From Green s representation theorem (Burton & Miller 1971): ( ) 1 u Φ(x, y) u 2 n (x) + + iηφ(x, y) (y)ds(y) = f(x), x Γ, n(x) n where Γ f(x) := ui n (x) + iηui (x).
14 Bath, October 28th u i, incident wave u + k 2 u = 0 2D polygon u = 0 D Γ x 2 x 1 From Green s representation theorem: ( ) 1 u Φ(x, y) u 2 n (x) + + iηφ(x, y) (y)ds(y) = f(x), x Γ. n(x) n Γ Theorem (follows from Burton & Miller 1971, Selepov 1969) If η R, η 0, then this integral equation is uniquely solvable in L 2 (Γ).
15 Bath, October 28th u n (x) + Γ ( Φ(x, y) n(x) ) u + iηφ(x, y) (y)ds(y) = f(x), x Γ. n Conventional BEM: Apply a Galerkin method, approximating u/ n by a piecewise polynomial of degree P, leading to a linear system to solve with N degrees of freedom.
16 Bath, October 28th u n (x) + Γ ( Φ(x, y) n(x) ) u + iηφ(x, y) (y)ds(y) = f(x), x Γ. n Conventional BEM: Apply a Galerkin method, approximating u/ n by a piecewise polynomial of degree P, leading to a linear system to solve with N degrees of freedom. Problem: N of order of kl, where L is linear dimension, so cost is O(N 2 ) to compute full matrix and apply iterative solver.
17 Bath, October 28th u n (x) + Γ ( Φ(x, y) n(x) ) u + iηφ(x, y) (y)ds(y) = f(x), x Γ. n Conventional BEM: Apply a Galerkin method, approximating u/ n by a piecewise polynomial of degree P, leading to a linear system to solve with N degrees of freedom. Problem: N of order of kl, where L is linear dimension, so cost is O(N 2 ) to compute full matrix and apply iterative solver... or close to O(N) if a fast multipole method (e.g. Amini & Profit 2003, Darve 2004) is used. This is fantastic but still infeasible as kl.
18 Bath, October 28th u n (x) + Alternative: Γ ( Φ(x, y) n(x) ) u + iηφ(x, y) (y)ds(y) = f(x), x Γ. n Reduce N by using new basis functions, e.g. (i) approximate u/ n by taking a large number of plane waves and multiplying these by conventional piecewise polynomial basis functions (Perrey-Debain et al. 2003, 2004). This is very successful (in 2D, 3D, for acoustic/elastic waves and Neumann/impedance b.c.s), reducing number of degrees of freedom per wavelength from e.g to close to 2. However N still increases proportional to kl.
19 Bath, October 28th u n (x) + Alternative: Γ ( Φ(x, y) n(x) ) u + iηφ(x, y) (y)ds(y) = f(x), x Γ. n Reduce N by using new basis functions, e.g. (ii) for convex scatterers, remove some of the oscillation by factoring out the oscillation of the incident wave, i.e. writing u ui (y) = (y) F (y) n n and approximating F by a conventional BEM (e.g. Abboud et al. 1994, Darrigrand 2002, Bruno et al 2004, Domínguez, Graham and Smyshlyaev 2005).
20 Bath, October 28th Alternative: Reduce N by using new basis functions, e.g. (ii) for convex scatterers, remove some of the oscillation by factoring out the oscillation of the incident wave, i.e. writing u ui (y) = (y) F (y) n n ( ) and approximating F by a conventional BEM. For smooth obstacles this works well: equation ( ) holds with F (y) 2 on the illuminated side (physical optics) and F (y) 0 in the shadow zone.
21 Bath, October 28th (ii) for convex scatterers, remove some of the oscillation by factoring out the oscillation of the incident wave, i.e. writing u ui (y) = (y) F (y) n n ( ) and approximating F by a conventional BEM. Not very effective for non-smooth scatterers.
22 Bath, October 28th Understanding solution behaviour Let G(x, y) := Φ(x, y) Φ(x, y ) be the Dirichlet Green function for the left half-plane Ω. By Green s representation theorem, u(x) = u i (x) + u r G(x, y) (x) + u(y)ds(y), x Ω. n(y) Ω\Γ
23 Bath, October 28th Understanding solution behaviour In the left half-plane Ω, u(x) = u i (x) + u r (x) + u (x) = 2 ui n n (x)+2 Ω\Γ Ω\Γ G(x, y) n(y) u(y)ds(y) 2 Φ(x, y) n(x) n(y) u(y)ds(y), x γ = Ω Γ.
24 Bath, October 28th Explicitly, where s is distance along γ, and φ(s) and ψ(s) are k 1 u/ n and u, at distance s along γ, φ(s) = P.O. + i 2 [ e iks v + (s) + e iks v (s) ] where v + (s) := k 0 and F (z) := e iz H (1) 1 (z)/z F ( k(s s 0 ) ) e iks 0 ψ(s 0 )ds 0.
25 Bath, October 28th φ(s) = P.O. + i 2 [ e iks v + (s) + e iks v (s) ] where v + (s) := k 0 F ( k(s s 0 ) ) e iks 0 ψ(s 0 )ds 0. Now F (z) := e iz H (1) 1 (z)/z which is non-oscillatory, in that F (n) (z) = O(z 3/2 n ) as z.
26 Bath, October 28th φ(s) = P.O. + i 2 [ e iks v + (s) + e iks v (s) ] where v + (s) := k 0 F ( k(s s 0 ) ) e iks 0 ψ(s 0 )ds 0. Now F (z) := e iz H (1) 1 (z)/z which is non-oscillatory, in that F (n) (z) = O(z 3/2 n ) as z. v (n) + (s) = O(k n (ks) 1/2 n ) as ks.
27 Bath, October 28th where φ(s) = P.O. + i [ e iks v + (s) + e iks v (s) ] 2 ( k n v (n) + (s) = O (ks) 1/2 n) as ks and (by separation of variables local to the corner), k n v (n) + (s) = O ( (ks) α n) as ks 0, where α < 1/2 depends on the corner angle.
28 Bath, October 28th φ(s) = P.O. + i 2 [ e iks v + (s) + e iks v (s) ] where k n v (n) + (s) = O ( (ks) 1/2 n) as ks O ((ks) α n ) as ks 0, where α < 1/2 depends on the corner angle. Thus approximate φ(s) P.O. + i 2 [ e iks V + (s) + e iks V (s) ], where V + and V are piecewise polynomials on graded meshes.
29 Bath, October 28th Thus approximate φ(s) P.O. + i 2 [ e iks V + (s) + e iks V (s) ], where V + and V are piecewise polynomials on graded meshes. Figure 3: Scattering by a square
30 Bath, October 28th Thus approximate φ(s) P.O. + i 2 [ e iks V + (s) + e iks V (s) ], where V + and V are piecewise polynomials on graded meshes. Figure 4: Scattering by a square
31 Bath, October 28th Thus approximate φ(s) P.O. + i 2 [ e iks V + (s) + e iks V (s) ], where V + and V are piecewise polynomials on graded meshes. Figure 5: Scattering by a square
32 Bath, October 28th Approximation error Theorem: If V + is the best L 2 approximation from the approximation space, then where k 1/2 v + V + 2 C p n 1/2 (1 + log 1/2 (kl)) N p+1, N degrees of freedom p = polynomial degree L = max side length n = number of sides of polygon
33 Bath, October 28th Boundary integral equation method Integral equation in parametric form ϕ(s) + Kϕ(s) = F (s), where ϕ(s) := 1 u (x(s)) P.O.. k n Theorem. The operator (I + K) : L 2 (Γ) L 2 (Γ) is bijective with bounded inverse (I + K) 1 2 C, so that the integral equation has exactly one solution.
34 Bath, October 28th Boundary integral equation method Integral equation in parametric form ϕ(s) + Kϕ(s) = F (s), where ϕ(s) := 1 u (x(s)) P.O.. k n Difficulty 1 The operator (I + K) : L 2 (Γ) L 2 (Γ) is bijective with bounded inverse (I + K) 1 2 C(k), where the dependence of C(k) on k is not clear.
35 Bath, October 28th Approximation space: seek ϕ N (s) = M v j ρ j (s) V N, j=1 where ρ j (s) := e ±iks piecewise polynomial supported on graded mesh. Question: how do we compute v j?
36 Bath, October 28th Galerkin method To solve seek ϕ NG V N such that ϕ(s) + Kϕ(s) = F (s), (I + P NG K)ϕ NG = P NG F, where P NG is the orthogonal projection onto the approximation space. Equivalently (ϕ NG, ρ) + (Kϕ NG, ρ) = (F, ρ), ρ V N, M v j [(ρ j, ρ m ) + (Kρ j, ρ m )] = (F, ρ m ). j=1 If ρ j, ρ m supported on same side of polygon, integrals not oscillatory.
37 Bath, October 28th Galerkin method Theorem. For N N, the operator (I + P NG K) : L 2 (Γ) V N is bijective with bounded inverse (I + P NG K) 1 2 C s.
38 Bath, October 28th Galerkin method Difficulty 2. For N N (k), the operator (I + P NG K) : L 2 (Γ) V N is bijective with bounded inverse (I + P NG K) 1 2 C s (k), where the dependence of N (k) and C s (k) on k is not clear.
39 Bath, October 28th Collocation method To solve seek ϕ NC V N such that ϕ(s) + Kϕ(s) = F (s), (I + P NC K)ϕ NC = P NC F, where P NC is the interpolatory projection onto the approximation space. Equivalently ϕ NC (s m ) + Kϕ NC (s m ) = F (s m ), m = 1,..., M, M v j [ρ j (s m ) + Kρ j (s m )] = F (s m ). j=1 If ρ j supported on same side of polygon as s m, integrals not oscillatory.
40 Bath, October 28th Collocation method We have not shown that (I + P NC K) : L 2 (Γ) V N is bijective with bounded inverse.
41 Bath, October 28th Galerkin vs. Collocation: error analysis Theorem There exists a constant C p > 0, independent of k, such that for N N k 1/2 ϕ ϕ NG 2 C p C s sup u(x) n1/2 (1 + log 1/2 (kl/n)) x D N p+1, k 1/2 u(x) u NG (x) C p C s sup u(x) n1/2 (1 + log 1/2 (kl/n)) x D N p+1. Stability and convergence not proven for collocation scheme.
42 Bath, October 28th Galerkin vs. Collocation: conditioning Galerkin: mass matrix M G := [(ρ j, ρ m )] has condm (1 + σ)/(1 σ), where min(y + j σ max, y m) max(y + j 1, y m 1 ) (y + j y+ j 1 )(y m ym 1 ) < 1, and if side lengths and angles are equal we can prove σ < ( ) 1/2N log k 1. kl Collocation: difficulty with choice of collocation points, M C := [ρ j (s m )] may be ill conditioned.
43 Bath, October 28th Galerkin vs. Collocation: implementation Galerkin: need to evaluate numerically many integrals of form a [ ] d H (1) 0 (k s 2 + t 2 ) + ith(1) 1 (k s 2 + t 2 ) e ik(σ jt σ m s) dt ds. s2 + t 2 b c Collocation: need to evaluate numerically many integrals of form [ ] b H (1) 0 (k s 2 m + t 2 ) + ith(1) 1 (k s 2 m + t 2 ) e ikσjt dt. s 2 m + t 2 a Collocation method easier to implement
44 Bath, October 28th Efficient evaluation of highly oscillatory integrals Consider Conventional wisdom... b a g(x)e ikf(x) dx. oscillatory integrals are difficult. oscillatory nature of solution increases errors. If we use interpolatory quadrature, this is certainly true - we need at least N quadrature points per wavelength, N 1, to achieve a reasonable level of accuracy.
45 Bath, October 28th Figure 6: Interpolation of e ikx2, k = 10.
46 Bath, October 28th Figure 7: Interpolation of e ikx2, k = 10.
47 Bath, October 28th Figure 8: Interpolation of e ikx2, k = 10.
48 Bath, October 28th Figure 9: Interpolation of e ikx2, k = 10.
49 Bath, October 28th Figure 10: Interpolation of e ikx2, k = 100.
50 Bath, October 28th Figure 11: Interpolation of e ikx2, k = 100.
51 Bath, October 28th Figure 12: Interpolation of e ikx2, k = 100.
52 Bath, October 28th Figure 13: Interpolation of e ikx2, k = 100.
53 Bath, October 28th Some alternative approaches 1. Filon quadrature (Iserles, N orsett)
54 Bath, October 28th Some alternative approaches 1. Filon quadrature (Iserles, Nörsett)... rather than interpolating the whole function g(x)e ikf(x), we just interpolate the non-oscillatory part, i.e. b a g(x)e ikf(x) dx b a [ gj χ j (x)] e ikf(x) dx, = b g j χ j (x)e ikf(x) dx. a
55 Bath, October 28th Some alternative approaches 1. Filon quadrature (Iserles, Nörsett)... rather than interpolating the whole function g(x)e ikf(x), we just interpolate the non-oscillatory part, i.e. b a g(x)e ikf(x) dx g j b a χ j (x)e ikf(x) dx. If hk 1, error O(h p+1 ). If hk 1, error O(h p+1 /(hk) 2 ). So performance actually improves as k increases.
56 Bath, October 28th Some alternative approaches 1. Filon quadrature (Iserles, Nörsett)... rather than interpolating the whole function g(x)e ikf(x), we just interpolate the non-oscillatory part, i.e. b a g(x)e ikf(x) dx g j b a χ j (x)e ikf(x) dx. If hk 1, error O(h p+1 ). If hk 1, error O(h p+1 /(hk) 2 ). So performance actually improves as k increases. However, we need to know the moments, hence for scattering problems application of Filon quadrature is difficult.
57 Bath, October 28th Some alternative approaches 1. Filon quadrature (Iserles, Nörsett) 2. Steepest descents (Huybrechs, Vanderwalle)
58 Bath, October 28th Some alternative approaches 1. Filon quadrature (Iserles, Nörsett) 2. Steepest descents (Huybrechs, Vanderwalle) 3. Asymptotic methods (Bruno, Reitich et. al., Iserles, Nörsett)
59 Bath, October 28th Asymptotic methods - the good news Oscillation leads to cancellation, and we can take advantage of this. Provided f (x) 0, b a g(x)e ikf(x) dx = 1 k [boundary terms] + i k b a g 1 (x)e ikf(x) dx, = 1 k [boundary terms] + 1 [boundary terms] k2 1 b k 2 g 2 (x)e ikf(x) dx, = n j=1 a boundary terms k n + O ( ) 1 k n+1.
60 Bath, October 28th Asymptotic methods - the bad news b a g(x)e ikf(x) dx = n j=1 boundary terms k n + ( ) n i b g n e ikf(x) dx, k a where g 0 = g, g n = ( gn 1 f ) =... = O ( ) 2n 1 f. As an asymptotic estimate - no problems. As an error estimate, care is needed, especially if f is small.
61 Bath, October 28th scattering by a square, k = 5 scattering by a square, k = 10 Numerical results
62 Bath, October 28th Numerical results (scattering by a square) Solution minus P.O. approximation;
63 Bath, October 28th Numerical results (scattering by a square) Solution minus P.O. approximation;
64 Bath, October 28th Numerical results (scattering by a square) Solution minus P.O. approximation;
65 Bath, October 28th Numerical results (scattering by a square) Solution minus P.O. approximation;
66 Bath, October 28th Numerical results (scattering by a square) Solution minus P.O. approximation;
67 Bath, October 28th Numerical results (scattering by a square) Solution minus P.O. approximation;
68 Bath, October 28th Numerical results (scattering by a square) Solution minus P.O. approximation;
69 Bath, October 28th Numerical results (scattering by a square) Exact solution minus P.O. approximation, k = 5;
70 Bath, October 28th Numerical results (scattering by a square) Exact solution minus P.O. approximation, k = 10;
71 Bath, October 28th Numerical results (scattering by a square) Exact solution minus P.O. approximation, k = 20;
72 Bath, October 28th Numerical results (scattering by a square) Exact solution minus P.O. approximation, k = 40;
73 Bath, October 28th Table 1: Relative errors, k = 10 k N dof ϕ ϕ NG 2 ϕ 2 ϕ ϕ NC 2 ϕ
74 Bath, October 28th Table 2: Relative errors, k = 160 k N dof ϕ ϕ NG 2 ϕ 2 ϕ ϕ NC 2 ϕ
75 Bath, October 28th k M N ϕ ϕ N 2 ϕ ϕ N 2 / ϕ 2 COND Table 3: Errors and relative L 2 errors, various k, N = 64
76 Bath, October 28th k N u N u 256 u 256 ( π, 3π) u N u 256 u 256 (3π, 3π) u N u 256 u 256 (3π, π) Table 4: Relative errors, for u N (x)
77 Bath, October 28th k N u N u 256 u 256 ( π, 3π) u N u 256 u 256 (3π, 3π) u N u 256 u 256 (3π, π) Table 5: Relative errors, for u N (x)
78 Bath, October 28th What we actually are computing... The difference between the exact solution and a leading order approximation; Figure 14: square, k = 5
79 Bath, October 28th What we actually are computing... The difference between the exact solution and a leading order approximation; Figure 15: square, k = 10
80 Bath, October 28th What we actually are computing... The difference between the exact solution and a leading order approximation; Figure 16: square, k = 20
81 Bath, October 28th What we actually are computing... The difference between the exact solution and a leading order approximation; Figure 17: square, k = 40
82 Bath, October 28th Summary and Conclusions Using Green s representation theorem in a half-plane we can understand behaviour of the field on the boundary and its derivatives for scattering by a convex polygon (extends to convex polyhedron in 3D) For a convex polygon, design of an optimal graded mesh for piecewise polynomial approximation is then straightforward The number of degrees of freedom need only grow logarithmically with the wavenumber to maintain a fixed accuracy Ongoing considerations Galerkin vs. Collocation - stability and convergence analysis Better schemes for evaluating oscillatory integrals hp ideas
High Frequency Scattering by Convex Polygons Stephen Langdon
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 of Maryland, September 2005
More informationA GALERKIN BOUNDARY ELEMENT METHOD FOR HIGH FREQUENCY SCATTERING BY CONVEX POLYGONS
A GALERKIN BOUNDARY ELEMENT METHOD FOR HIGH FREQUENCY SCATTERING BY CONVEX POLYGONS S N CHANDLER-WILDE AND S LANGDON Abstract. In this paper we consider the problem of time-harmonic acoustic scattering
More informationError analysis and fast solvers for high-frequency scattering problems
Error analysis and fast solvers for high-frequency scattering problems I.G. Graham (University of Bath) Woudschoten October 2014 High freq. problem for the Helmholtz equation Given an object Ω R d, with
More informationA Galerkin boundary element method for high frequency scattering by convex polygons
A Galerkin boundary element method for high frequency scattering by convex polygons Article Published Version Chandler Wilde, S. N. and Langdon, S. (2007) A Galerkin boundary element method for high frequency
More informationHigh Frequency Boundary Element Methods for Scattering by Convex Polygons
High Frequency Boundary Element Methods for Scattering by Convex Polygons J. D. Coleman August 1, 006 1 Acknowledgements I would like to acknowledge the EPSRC for financially supporting this work. I would
More informationA GALERKIN BOUNDARY ELEMENT METHOD FOR HIGH FREQUENCY SCATTERING BY CONVEX POLYGONS
A GALERKIN BOUNDARY ELEMENT METHOD FOR HIGH FREQUENCY SCATTERING BY CONVEX POLYGONS READING UNIVERSITY NUMERICAL ANALYSIS REPORT 7/06 S N CHANDLER-WILDE AND S LANGDON Abstract. In this paper we consider
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 informationWhen is the error in the h BEM for solving the Helmholtz equation bounded independently of k?
BIT manuscript No. (will be inserted by the editor) When is the error in the h BEM for solving the Helmholtz equation bounded independently of k? I. G. Graham M. Löhndorf J. M. Melenk E. A. Spence Received:
More informationAcoustic scattering : high frequency boundary element methods and unified transform methods
Acoustic scattering : high frequency boundary element methods and unified transform methods Book or Report Section Accepted Version Chandler Wilde, S.N. and Langdon, S. (2015) Acoustic scattering : high
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 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 informationEfficient evaluation of highly oscillatory acoustic scattering surface integrals
Efficient evaluation of highly oscillatory acoustic scattering surface integrals M. Ganesh a, S. Langdon b,,1, I. H. Sloan c a Department of Mathematical and Computer Sciences, Colorado School of Mines,
More informationarxiv: v1 [math.na] 11 Mar 2019 A high frequency boundary element method for scattering by a class of multiple obstacles
(2019) Page 1 of 39 arxiv:1903.04449v1 [math.na] 11 Mar 2019 A high frequency boundary element method for scattering by a class of multiple obstacles ANDREW GIBBS, DEPT. OF COMPUTER SCIENCE, K.U. LEUVEN,
More informationDepartment of Mathematics and Statistics
School of Mathematical and Physical Sciences Department of Mathematics and Statistics Preprint MPS-11-18 5 November 11 A High Frequency hp Boundary Element Method for Scattering by Convex Polygons by D.P.
More informationBoundary integral methods in high frequency scattering
Boundary integral methods in high frequency scattering Simon N. Chandler-Wilde Ivan G. Graham Abstract In this article we review recent progress on the design, analysis and implementation of numerical-asymptotic
More informationA High Frequency Boundary Element Method for Scattering by a Class of Nonconvex Obstacles
School of Mathematical and Physical Sciences Department of Mathematics and Statistics Preprint MPS-2012-04 28 May 2012 A High Frequency Boundary Element Method for Scattering by a Class of Nonconvex Obstacles
More informationA high frequency boundary element method for scattering by a class of nonconvex obstacles
A high frequency boundary element method for scattering by a class of nonconvex obstacles Article Accepted Version Chandler Wilde, S. N., Hewett, D. P., Langdon, S. and Twigger, A. (2015) A high frequency
More informationShadow boundary effects in hybrid numerical-asymptotic methods for high frequency scattering
Shadow boundary effects in hybrid numerical-asymptotic methods for high frequency scattering arxiv:1409.3776v [math.na] 3 Jun 015 D. P. Hewett Department of Mathematics and Statistics, University of Reading,
More informationElectromagnetic wave propagation. ELEC 041-Modeling and design of electromagnetic systems
Electromagnetic wave propagation ELEC 041-Modeling and design of electromagnetic systems EM wave propagation In general, open problems with a computation domain extending (in theory) to infinity not bounded
More informationHigh Frequency Solution Behaviour in Wave Scattering Problems Simon Chandler-Wilde University of Reading, UK
High Frequency Solution Behaviour in Wave Scattering Problems Simon Chandler-Wilde University of Reading, UK www.reading.ac.uk/~sms03snc With: Steve Langdon, Andrea Moiola (Reading), Ivan Graham, Euan
More informationNumerical methods for high frequency scattering by multiple obstacles
Numerical methods for high frequency scattering by multiple obstacles Andrew Gibbs University of Reading A thesis submitted for the degree of Doctor of Philosophy March 2017 Declaration I confirm that
More informationFast Hybrid Algorithms for High Frequency Scattering
Fast Hybrid Algorithms for High Frequency Scattering Björn Engquist, Khoa Tran, and Lexing Ying Citation: AIP Conference Proceedings 1106, 3 (2009); doi: 10.1063/1.3117111 View online: http://dx.doi.org/10.1063/1.3117111
More informationInverse Obstacle Scattering
, Göttingen AIP 2011, Pre-Conference Workshop Texas A&M University, May 2011 Scattering theory Scattering theory is concerned with the effects that obstacles and inhomogenities have on the propagation
More informationA new 9-point sixth-order accurate compact finite difference method for the Helmholtz equation
A new 9-point sixth-order accurate compact finite difference method for the Helmholtz equation Majid Nabavi, M. H. Kamran Siddiqui, Javad Dargahi Department of Mechanical and Industrial Engineering, Concordia
More informationA high frequency $hp$ boundary element method for scattering by convex polygons
A high frequency $hp$ boundary element method for scattering by convex polygons Article Published Version Hewett, D. P., Langdon, S. and Melenk, J. M. (213) A high frequency $hp$ boundary element method
More informationA Survey of Computational High Frequency Wave Propagation II. Olof Runborg NADA, KTH
A Survey of Computational High Frequency Wave Propagation II Olof Runborg NADA, KTH High Frequency Wave Propagation CSCAMM, September 19-22, 2005 Numerical methods Direct methods Wave equation (time domain)
More informationConvergence and optimality of an adaptive FEM for controlling L 2 errors
Convergence and optimality of an adaptive FEM for controlling L 2 errors Alan Demlow (University of Kentucky) joint work with Rob Stevenson (University of Amsterdam) Partially supported by NSF DMS-0713770.
More 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 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 informationA frequency-independent boundary element method for scattering by two-dimensional screens and apertures
A frequency-independent boundary element method for scattering by two-dimensional screens and apertures arxiv:1401.2786v2 [math.na] 11 Aug 2014 D. P. Hewett, S. Langdon and S. N. Chandler-Wilde Department
More informationFrequency Domain Integral Equations for Acoustic and Electromagnetic Scattering Problems
Frequency Domain Integral Equations for Acoustic and Electromagnetic Scattering Problems Christophe Geuzaine, U of Liege Fernando Reitich, U of Minnesota Catalin Turc, UNC Charlotte Outline Integral equations
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 informationA Numerical Approach To The Steepest Descent Method
A Numerical Approach To The Steepest Descent Method K.U. Leuven, Division of Numerical and Applied Mathematics Isaac Newton Institute, Cambridge, Feb 15, 27 Outline 1 One-dimensional oscillatory integrals
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 information1. How do error estimators for the Galerkin FEM depend on &, ft? 2. Can the classical Galerkin approach be improved towards a linear rule?
Reliability of finite element methods for the numerical computation of waves F. Ihlenburg*, I. Babuska*, S. Sauted "Institute for Physical Science and Technology, University of Maryland at College Park,
More informationON THE INFLUENCE OF THE WAVENUMBER ON COMPRESSION IN A WAVELET BOUNDARY ELEMENT METHOD FOR THE HELMHOLTZ EQUATION
INTERNATIONAL JOURNAL OF NUMERICAL ANALYSIS AND MODELING Volume 4, Number 1, Pages 48 62 c 2007 Institute for Scientific Computing and Information ON THE INFLUENCE OF THE WAVENUMBER ON COMPRESSION IN A
More informationPolynomial Approximation: The Fourier System
Polynomial Approximation: The Fourier System Charles B. I. Chilaka CASA Seminar 17th October, 2007 Outline 1 Introduction and problem formulation 2 The continuous Fourier expansion 3 The discrete Fourier
More informationThe Factorization Method for Inverse Scattering Problems Part I
The Factorization Method for Inverse Scattering Problems Part I Andreas Kirsch Madrid 2011 Department of Mathematics KIT University of the State of Baden-Württemberg and National Large-scale Research Center
More informationFinite and Boundary Element Methods in Acoustics
Finite and Boundary Element Methods in Acoustics W. Kreuzer, Z. Chen, H. Waubke Austrian Academy of Sciences, Acoustics Research Institute ARI meets NuHAG Kreuzer, Chen, Waubke (ARI) FEM-BEM-FMM ARI meets
More informationA far-field based T-matrix method for three dimensional acoustic scattering
ANZIAM J. 50 (CTAC2008) pp.c121 C136, 2008 C121 A far-field based T-matrix method for three dimensional acoustic scattering M. Ganesh 1 S. C. Hawkins 2 (Received 14 August 2008; revised 4 October 2008)
More informationNumerical Solutions of Laplacian Problems over L-Shaped Domains and Calculations of the Generalized Stress Intensity Factors
WCCM V Fifth World Congress on Computational Mechanics July 7-2, 2002, Vienna, Austria Eds.: H.A. Mang, F.G. Rammerstorfer, J. Eberhardsteiner Numerical Solutions of Laplacian Problems over L-Shaped Domains
More informationInverse scattering problem from an impedance obstacle
Inverse Inverse scattering problem from an impedance obstacle Department of Mathematics, NCKU 5 th Workshop on Boundary Element Methods, Integral Equations and Related Topics in Taiwan NSYSU, October 4,
More informationarxiv: v1 [physics.comp-ph] 28 Jan 2008
A new method for studying the vibration of non-homogeneous membranes arxiv:0801.4369v1 [physics.comp-ph] 28 Jan 2008 Abstract Paolo Amore Facultad de Ciencias, Universidad de Colima, Bernal Díaz del Castillo
More informationFast Multipole BEM for Structural Acoustics Simulation
Fast Boundary Element Methods in Industrial Applications Fast Multipole BEM for Structural Acoustics Simulation Matthias Fischer and Lothar Gaul Institut A für Mechanik, Universität Stuttgart, Germany
More informationA GWBEM method for high frequency acoustic scattering
A GWBEM method for high frequency acoustic scattering E. Perrey-Debain, H. Bériot, M. Ben Tahar and C. Vayssade Université de Technologie de Compiègne, Lab. Roberval UMR 653, Dept. Acoustique, BP 60319,
More informationModelling in photonic crystal structures
Modelling in photonic crystal structures Kersten Schmidt MATHEON Nachwuchsgruppe Multiscale Modelling and Scientific Computing with PDEs in collaboration with Dirk Klindworth (MATHEON, TU Berlin) Holger
More 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 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 informationAnton ARNOLD. with N. Ben Abdallah (Toulouse), J. Geier (Vienna), C. Negulescu (Marseille) TU Vienna Institute for Analysis and Scientific Computing
TECHNISCHE UNIVERSITÄT WIEN Asymptotically correct finite difference schemes for highly oscillatory ODEs Anton ARNOLD with N. Ben Abdallah (Toulouse, J. Geier (Vienna, C. Negulescu (Marseille TU Vienna
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 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 informationWeighted Residual Methods
Weighted Residual Methods Introductory Course on Multiphysics Modelling TOMASZ G. ZIELIŃSKI bluebox.ippt.pan.pl/ tzielins/ Table of Contents Problem definition. oundary-value Problem..................
More informationScientific Computing I
Scientific Computing I Module 8: An Introduction to Finite Element Methods Tobias Neckel Winter 2013/2014 Module 8: An Introduction to Finite Element Methods, Winter 2013/2014 1 Part I: Introduction to
More informationOn domain decomposition preconditioners for finite element approximations of the Helmholtz equation using absorption
On domain decomposition preconditioners for finite element approximations of the Helmholtz equation using absorption Ivan Graham and Euan Spence (Bath, UK) Collaborations with: Paul Childs (Emerson Roxar,
More informationShifted Laplace and related preconditioning for the Helmholtz equation
Shifted Laplace and related preconditioning for the Helmholtz equation Ivan Graham and Euan Spence (Bath, UK) Collaborations with: Paul Childs (Schlumberger Gould Research), Martin Gander (Geneva) Douglas
More informationFundamental Solutions and Green s functions. Simulation Methods in Acoustics
Fundamental Solutions and Green s functions Simulation Methods in Acoustics Definitions Fundamental solution The solution F (x, x 0 ) of the linear PDE L {F (x, x 0 )} = δ(x x 0 ) x R d Is called the fundamental
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 informationFEM-FEM and FEM-BEM Coupling within the Dune Computational Software Environment
FEM-FEM and FEM-BEM Coupling within the Dune Computational Software Environment Alastair J. Radcliffe Andreas Dedner Timo Betcke Warwick University, Coventry University College of London (UCL) U.K. Radcliffe
More informationImproved near-wall accuracy for solutions of the Helmholtz equation using the boundary element method
Center for Turbulence Research Annual Research Briefs 2006 313 Improved near-wall accuracy for solutions of the Helmholtz equation using the boundary element method By Y. Khalighi AND D. J. Bodony 1. Motivation
More informationSimple Examples on Rectangular Domains
84 Chapter 5 Simple Examples on Rectangular Domains In this chapter we consider simple elliptic boundary value problems in rectangular domains in R 2 or R 3 ; our prototype example is the Poisson equation
More informationAMS 529: Finite Element Methods: Fundamentals, Applications, and New Trends
AMS 529: Finite Element Methods: Fundamentals, Applications, and New Trends Lecture 3: Finite Elements in 2-D Xiangmin Jiao SUNY Stony Brook Xiangmin Jiao Finite Element Methods 1 / 18 Outline 1 Boundary
More 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 informationCalculation of Highly Oscillatory Integrals by Quadrature Method
Calculation of Highly Oscillatory Integrals by Quadrature Methods Krishna Thapa Texas A&M University, College Station, TX May 18, 211 Outline Purpose Why Filon Type Quadrature Methods Our Expectations
More informationDiscrete Projection Methods for Integral Equations
SUB Gttttingen 7 208 427 244 98 A 5141 Discrete Projection Methods for Integral Equations M.A. Golberg & C.S. Chen TM Computational Mechanics Publications Southampton UK and Boston USA Contents Sources
More informationFast methods for integral equations
Fast methods for integral equations Matthieu Aussal Research engineer Centre de Mathématiques Appliquées de l École Polytechnique Route de Saclay - 91128 Palaiseau CEDEX France Journée de rentrée du CMAP
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 informationAdaptive Collocation with Kernel Density Estimation
Examples of with Kernel Density Estimation Howard C. Elman Department of Computer Science University of Maryland at College Park Christopher W. Miller Applied Mathematics and Scientific Computing Program
More informationarxiv: v3 [math.ap] 4 Jan 2017
Recovery of an embedded obstacle and its surrounding medium by formally-determined scattering data Hongyu Liu 1 and Xiaodong Liu arxiv:1610.05836v3 [math.ap] 4 Jan 017 1 Department of Mathematics, Hong
More informationSubdiffusion in a nonconvex polygon
Subdiffusion in a nonconvex polygon Kim Ngan Le and William McLean The University of New South Wales Bishnu Lamichhane University of Newcastle Monash Workshop on Numerical PDEs, February 2016 Outline Time-fractional
More informationA DIRECT IMAGING ALGORITHM FOR EXTENDED TARGETS
A DIRECT IMAGING ALGORITHM FOR EXTENDED TARGETS SONGMING HOU, KNUT SOLNA, AND HONGKAI ZHAO Abstract. We present a direct imaging algorithm for extended targets. The algorithm is based on a physical factorization
More informationNon-Conforming Finite Element Methods for Nonmatching Grids in Three Dimensions
Non-Conforming Finite Element Methods for Nonmatching Grids in Three Dimensions Wayne McGee and Padmanabhan Seshaiyer Texas Tech University, Mathematics and Statistics (padhu@math.ttu.edu) Summary. In
More informationAnalysis of thin acoustic barriers over an
Analysis of thin acoustic barriers over an impedance plane using a dual BEM formulation L.A. de Lacerda*\ L.C. Wrobel*, WJ. Mansur^ "Wessex Institute of Technology, Ashurst, Southampton, UK ^COPPE/URRJ,
More informationfor compression of Boundary Integral Operators. Steven Paul Nixon B.Sc.
Theory and Applications of the Multiwavelets for compression of Boundary Integral Operators. Steven Paul Nixon B.Sc. Institute for Materials Research School of Computing, Science & Engineering, University
More informationA frequency-independent boundary element method for scattering by two-dimensional screens and apertures
IMA Journal of Numerical Analysis Advance Access published October 16, 2014 IMA Journal of Numerical Analysis (2014) Page 1 of 31 doi:10.1093/imanum/dru043 A frequency-independent boundary element method
More informationMaximum norm estimates for energy-corrected finite element method
Maximum norm estimates for energy-corrected finite element method Piotr Swierczynski 1 and Barbara Wohlmuth 1 Technical University of Munich, Institute for Numerical Mathematics, piotr.swierczynski@ma.tum.de,
More informationFastMMLib: a generic Fast Multipole Method library. Eric DARRIGRAND. IRMAR Université de Rennes 1
FastMMLib: a generic Fast Multipole Method library Eric DARRIGRAND joint work with Yvon LAFRANCHE IRMAR Université de Rennes 1 Outline Introduction to FMM Historically SVD principle A first example of
More informationScientific Computing WS 2018/2019. Lecture 15. Jürgen Fuhrmann Lecture 15 Slide 1
Scientific Computing WS 2018/2019 Lecture 15 Jürgen Fuhrmann juergen.fuhrmann@wias-berlin.de Lecture 15 Slide 1 Lecture 15 Slide 2 Problems with strong formulation Writing the PDE with divergence and gradient
More 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 informationFinite Element Analysis of Acoustic Scattering
Frank Ihlenburg Finite Element Analysis of Acoustic Scattering With 88 Illustrations Springer Contents Preface vii 1 The Governing Equations of Time-Harmonic Wave Propagation, 1 1.1 Acoustic Waves 1 1.1.1
More informationPART IV Spectral Methods
PART IV Spectral Methods Additional References: R. Peyret, Spectral methods for incompressible viscous flow, Springer (2002), B. Mercier, An introduction to the numerical analysis of spectral methods,
More informationThe Closed Form Reproducing Polynomial Particle Shape Functions for Meshfree Particle Methods
The Closed Form Reproducing Polynomial Particle Shape Functions for Meshfree Particle Methods by Hae-Soo Oh Department of Mathematics, University of North Carolina at Charlotte, Charlotte, NC 28223 June
More informationQBX and the DPIE for the Maxwell Equations
University of Illinois @ Urbana-Champaign Fall 2017 - CS 598 APK The Goal For this project, the goals were to implement a model for the Python package pytential for tackling the Maxwell Equations for perfect
More informationAn asymptotic model for computing the scattering of waves by small heterogeneities in the time domain
An asymptotic model for computing the scattering of waves by small heterogeneities in the time domain Simon Marmorat EPI POems INRIA Saclay / CEA List Joint work with Patrick Joly (EPI POems INRIA/ENSTA/CNRS)
More informationarxiv: v2 [math.na] 13 Jan 2018
A Galerkin BEM for high-frequency scattering problems based on frequency dependent changes of variables arxiv:609.06v [math.na] 3 Jan 08 Fatih Ecevit Boğaziçi University, Department of Mathematics, Bebek
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 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 informationInverse obstacle scattering problems using multifrequency measurements
Inverse obstacle scattering problems using multifrequency measurements Nguyen Trung Thành Inverse Problems Group, RICAM Joint work with Mourad Sini *** Workshop 3 - RICAM special semester 2011 Nov 21-25
More informationThe Pseudospectral Method
The Pseudospectral Method Heiner Igel Department of Earth and Environmental Sciences Ludwig-Maximilians-University Munich Heiner Igel Computational Seismology 1 / 43 Outline 1 Introduction Motivation History
More informationFinite Elements. Colin Cotter. January 15, Colin Cotter FEM
Finite Elements January 15, 2018 Why Can solve PDEs on complicated domains. Have flexibility to increase order of accuracy and match the numerics to the physics. has an elegant mathematical formulation
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 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 informationPROJECTION METHODS FOR DYNAMIC MODELS
PROJECTION METHODS FOR DYNAMIC MODELS Kenneth L. Judd Hoover Institution and NBER June 28, 2006 Functional Problems Many problems involve solving for some unknown function Dynamic programming Consumption
More informationVon Neumann Analysis of Jacobi and Gauss-Seidel Iterations
Von Neumann Analysis of Jacobi and Gauss-Seidel Iterations We consider the FDA to the 1D Poisson equation on a grid x covering [0,1] with uniform spacing h, h (u +1 u + u 1 ) = f whose exact solution (to
More informationRobust exponential convergence of hp-fem for singularly perturbed systems of reaction-diffusion equations
Robust exponential convergence of hp-fem for singularly perturbed systems of reaction-diffusion equations Christos Xenophontos Department of Mathematics and Statistics University of Cyprus joint work with
More informationKirchhoff, Fresnel, Fraunhofer, Born approximation and more
Kirchhoff, Fresnel, Fraunhofer, Born approximation and more Oberseminar, May 2008 Maxwell equations Or: X-ray wave fields X-rays are electromagnetic waves with wave length from 10 nm to 1 pm, i.e., 10
More informationComputer simulation of multiscale problems
Progress in the SSF project CutFEM, Geometry, and Optimal design Computer simulation of multiscale problems Axel Målqvist and Daniel Elfverson University of Gothenburg and Uppsala University Umeå 2015-05-20
More informationThe method of lines (MOL) for the diffusion equation
Chapter 1 The method of lines (MOL) for the diffusion equation The method of lines refers to an approximation of one or more partial differential equations with ordinary differential equations in just
More informationFast and accurate methods for the discretization of singular integral operators given on surfaces
Fast and accurate methods for the discretization of singular integral operators given on surfaces James Bremer University of California, Davis March 15, 2018 This is joint work with Zydrunas Gimbutas (NIST
More informationImplementation of Sparse Wavelet-Galerkin FEM for Stochastic PDEs
Implementation of Sparse Wavelet-Galerkin FEM for Stochastic PDEs Roman Andreev ETH ZÜRICH / 29 JAN 29 TOC of the Talk Motivation & Set-Up Model Problem Stochastic Galerkin FEM Conclusions & Outlook Motivation
More informationFEM Convergence for PDEs with Point Sources in 2-D and 3-D
FEM Convergence for PDEs with Point Sources in 2-D and 3-D Kourosh M. Kalayeh 1, Jonathan S. Graf 2 Matthias K. Gobbert 2 1 Department of Mechanical Engineering 2 Department of Mathematics and Statistics
More information