Model reduction of wave propagation
|
|
- Hugo Dean
- 5 years ago
- Views:
Transcription
1 Model reduction of wave propagation via phase-preconditioned rational Krylov subspaces Delft University of Technology and Schlumberger V. Druskin, R. Remis, M. Zaslavsky, Jörn Zimmerling November 8, 2017 Jörn Zimmerling (TU Delft) Model reduction of wave propagation November 8, / 37
2 Motivation Second order wave equation with wave operator A Receiver Source PML 5500 Au s 2 u = δ(x x S ) Assume N grid steps in every spatial direction Scaling of surface seismic in 3D: # Grid points O(N 3 ) # Sources O(N 2 ) # Frequencies O(N) Overall O(N 6 )ψ(n 3 ) x-direction [m] y-direction [m] Wavespeed [m/s] (a) Section of wave speed profile of the Marmousi model. Jörn Zimmerling (TU Delft) Model reduction of wave propagation November 8, / 37
3 Goal of this work Simulate and compress large scale wave fields in modern high performance computing environment (parallel CPU and GPU environment) Use projection based model order reduction to Approximate transfer function reduce # of frequencies needed to solve reduce # of sources to be considered reduce # number of grid points needed Jörn Zimmerling (TU Delft) Model reduction of wave propagation November 8, / 37
4 Introduction Simulating and compressing large scale wave fields Au [l] s 2 u [l] = δ(x x [l] S ), (1) With the wave operator given by A ν 2, Laplace frequency s We consider a Multiple-Input Multiple-Output problem Define fields U = [u [1], u [2],..., u [Ns] ] and sources B = [ δ(x x [1] S ), δ(x x[2] S ),..., δ(x x[ns] S )] We are interested in the transfer function (Receivers and Sources coincide) F(s) = B H U(s)dx (2) Open Domains Jörn Zimmerling (TU Delft) Model reduction of wave propagation November 8, / 37
5 Problem Formulation After finite difference discretization with PML (A(s) s 2 I)U = ˆB Step sizes inside the PML h i = α i + β i s Frequency dependent A(s) caused by absorbing boundary Q(s)U = B with Q(s) C N N Q(s) propterties sparse complex symmetric (reciprocity) Schwarz reflection principle Q(s) = Q(s) passive (nonlinear NR 1 Re < 0) 1 NR:W {A(s)} = { s C : x H A(s)x = 0 x C k \0 } Jörn Zimmerling (TU Delft) Model reduction of wave propagation November 8, / 37
6 Problem Formulation Transfer function from sources to receivers F(B, s) = B H Q(s) 1 B Reduced order modeling of transfer function over frequency range F m (B, s) = V m B H (VmQ(s)V H m ) 1 VmB H with V m C N m valid in a range of s, and easy to store Motivation FD grid overdiscretized w.r.t. Nyquist approximation F(B, s) to noise level PML introduces losses limited I/O map Jörn Zimmerling (TU Delft) Model reduction of wave propagation November 8, / 37
7 Outline 1 Problem formulation 2 Model order reduction Rational Krylov subspaces 3 Phase-Preconditioning 4 Numerical Experiments 5 Conclusions Jörn Zimmerling (TU Delft) Model reduction of wave propagation November 8, / 37
8 Structure preserving rational Krylov subspaces Preserve: symmetry (w.r.t. matrix, frequency), passivity Block rational Krylov subspace approach κ = s 1,..., s m K m (κ) = span { Q(s 1 ) 1 B,..., Q(s m ) 1 B } K 2m Re = span {Km (κ), K m (κ)} Let V m be a (real) basis for KRe m then with reduced order model (via Galerkin condition) we approximate R m (s) = V H mq(s)v m F m (s) = (V H mb) H R m (s) 1 V H mb, Jörn Zimmerling (TU Delft) Model reduction of wave propagation November 8, / 37
9 Structure preserving rational Krylov subspaces Numerical Range: W {R m (s)} W {Q(s)} Proof: x H mr m (s)x m = (V m x m ) H Q(s)(V m x m ) R m (s) is passive F m (s) is Hermite interpolant of F(s) on κ κ Q(κ) 1 B KRe 2m + uniqueness of Galerkin Schwarz reflection and symmetry hold aswell Jörn Zimmerling (TU Delft) Model reduction of wave propagation November 8, / 37
10 RKS for a Resonant cavity 6 5 FDFD Response RKS Response m=60 RKS Response m=20 Response [a.u] x-direction Receiver Source PML Normalized Frequency y-direction RKS has excellent convergence if singular Hankel values of system decay fast (few contributing eigenvectors) F m (s) is a [2m 1/2m] rational function Jörn Zimmerling (TU Delft) Model reduction of wave propagation November 8, / 37
11 Problem of RKS in Geophysics - Nyquist Limit 0.2 Response [a.u.] Normalized Frequency Long travel times δ(t T) F exp( st) Nyquist sampling of ω = π T F(s) = B H Q(s) 1 Bdx is oscillatory Jörn Zimmerling (TU Delft) Model reduction of wave propagation November 8, / 37
12 Filion Quadrature Filion quadrature deals with oscillatory integral F(s) = exp(st) f(t)dt quadrature requires s t 1 F(s) t n a n exp(s n t) f(n t) Filion quadrature makes a n function of s t F(s) t n a n (s t) exp(s n t) f(n t) Make projection basis s dependent Frequency dependence from asymptotic s i (WKB) Jörn Zimmerling (TU Delft) Model reduction of wave propagation November 8, / 37
13 Phase-Preconditioning I - 1D We can overcome the Nyquist demand by splitting the wavefield into oscillatory and smooth part u(s j ) = exp( s j T eik )c out (s j ) + exp(s j T eik )c in (s j ). (3) Oscillatory phase term obtained from high frequency asymptotics Eikonal equation T [l] eik 2 = 1 ν 2 Amplitudes c out/in are smooth Motivated by Filon quadrature Handle oscillatory part analytically Quadrature with smooth amplitudes Note: Splitting not unique Jörn Zimmerling (TU Delft) Model reduction of wave propagation November 8, / 37
14 Phase-Preconditioning II Projection on frequency dependent Reduced Order Basis K 2m EIK (κ, s) = span{ exp( st eik) c out (s 1 ),..., exp( st eik ) c out (s m ), exp( st eik ) c in (s s ),..., exp( st eik ) c in (s m )} Preserve Schwartz reflection principle K 4m EIK;R (κ, s) = span { KEIK 2m (κ, s), K2m EIK (κ, s)} (4) equivalent to changing Operator Coefficients from Galerkin condition m m u m (s) = α i (s) exp( st eik ) c out (s i )+ β i (s) exp(st eik ) c in (s i )+... i=1 i=1 Jörn Zimmerling (TU Delft) Model reduction of wave propagation November 8, / 37
15 Phase-Preconditioning III Non-uniqueness of splitting resolved by one-way WEQ c out (s j ) = ν ( ) sj exp( s j T eik ) 2s j ν u(s j) x x S u(s j), (5) c in (s j ) = ν ( ) sj exp( s j T eik ) 2s j ν u(s j) + x x S u(s j). (6) Effects of Phase preconditioning on Number of Interpolation points Size of the computational Grid MIMO problems Computational Scheme Jörn Zimmerling (TU Delft) Model reduction of wave propagation November 8, / 37
16 Projection on frequency dependent space Let V m;eik (s) be a real basis of KEIK;R 4m (κ, s) The reduced order model is given by large inner products on GPU This preserves symmetry Schwarz reflection principle passivity Interpolation R m;eik (s) = V H m;eik (s)q(s)v m;eik(s). Jörn Zimmerling (TU Delft) Model reduction of wave propagation November 8, / 37
17 Number of interpolation points needed Double interpolation of transfer function still holds F m (s) = F m (s) and d ds F m(s) = d F(s) with s κ κ. (7) ds Number of interpolation point needed dependent on complexity medium, not latest arrival Proposition: Let a 1D problem have l homogenous layers. Then there exist m l + 1 non-coinciding interpolation points, such that the solution u m;eik (s) = u. Jörn Zimmerling (TU Delft) Model reduction of wave propagation November 8, / 37
18 Illustration of Proposition Wavespeed Source L1 L2 L3 L4 L5 L6 L7 Real part field Source L1 L2 L3 L4 L5 L6 L7 Imag(C) outgooing Source L1 L2 L3 L4 L5 L6 L7 Imag(C) incoming 2 0 Source L1 L2 L3 L4 L5 L6 L7 Amplitudes are constants in layers + left and right of source Basis is complete after l + 1 iterations Jörn Zimmerling (TU Delft) Model reduction of wave propagation November 8, / 37
19 Phase-Preconditioning higher dimensions/mimo Split with dimension specific function u(s j ) [l] = g(s j T [l] eik )c out(s j ) [l] + g( s j T [l] eik )c in(s j ) [l], (8) g(x) obtained from WKB approximation One way wave equations along T eik used for decomposition In 2D is we use g(x) = K 0 (x) for outgoing Multiple T [l] direction c in (s j ) = eik s j T sign(im (s j ))iπ for multiple sources [l] account for multiple ] [K 1 (s j T) u(s j ) + K 0 (s j T) ν2 T u(s j ) s j Jörn Zimmerling (TU Delft) Model reduction of wave propagation November 8, / 37
20 Size of the computational Grid Receiver Source PML x-direction [m] wavespeed [m/s] y-direction [m] (b) Section of the wave speed profile of the smoothed Marmousi model. (c) Real part of the wavefield u [4]. (d) Real part of the amplitude c [4] out. Jörn Zimmerling (TU Delft) Model reduction of wave propagation November 8, / 37
21 Numerical Experiments - I Configurations (Neumann boundary condition on top) Layered medium Travel time dominated 5 Sources and 5 Receivers x 4m Comp. Size 829x480 points Size 3160 m x 1920m range c m/s Range Quadrature 0-40 Hz x-direction [m] Receiver Source PML y-direction [m] wavespeed [m/s] (e) Section of the wave speed profile of the smoothed Marmousi model. Jörn Zimmerling (TU Delft) Model reduction of wave propagation November 8, / 37
22 Numerical Experiments response [a.u.] Full Response RKS m=20 PPRKS m=20 Interpolation points Source 1 to Receiver 5 PPRKS clearly outperforms RKS normalized frequency (f) Real part of the frequency-domain transfer function Jörn Zimmerling (TU Delft) Model reduction of wave propagation November 8, / 37
23 Numerical Experiments I Time-domain convergence of the RKS and the PPRKS 1 0 RKS PPRKS 10 log of error number of interpolation frequencies m (g) Jörn Zimmerling (TU Delft) Model reduction of wave propagation November 8, / 37
24 Computational Grid Amplitudes c in/out are much smoother than the wavefield ROM can extrapolate to high frequencies oscillatory part is handled analytically Two-grid approach: Amplitudes can be computed on coarse grid U c = Q course (s) 1 B c Interpolate amplitudes to fine grid Projection of operator and evaluation are performed on fine grid F c;m (s) = B H { [V c;m (s)] H Q fine (s)v c;m (s) } 1 B Solution gets gauged to the fine grid (no interpolation anymore) Jörn Zimmerling (TU Delft) Model reduction of wave propagation November 8, / 37
25 Phase-Preconditioning SVD Amplitudes are smooth in space and can become redundant Reduce amplitudes via SVD of [c out c in ] c j SVD Amplitudes have no source information 10log normalized singular values RKS Contracted amplitude basis index singular value (h) Singular values of normalized (c out c in ) # singular values larger than 0.01 versus # sources with m = 40 N src [c out, c in ] u m N src Reduction of sources Jörn Zimmerling (TU Delft) Model reduction of wave propagation November 8, / 37
26 Phase-Preconditioning SVD Assume s ir then we obtain [ u m [l] N src M SVD a [l] rj (s) = α [l] rj r=1 j=1 ] T [ ] g(st [r] eik )cj SVD g( st [r] eik )cj SVD (9) where M SVD 2mN src. Coefficients from Galerkin condition c j SVD is no longer source dependent Jörn Zimmerling (TU Delft) Model reduction of wave propagation November 8, / 37
27 Computational Scheme - Coarse Grid - CPU ROM Construction phase Initialize Simulation Compute T [l] eik Solve coarse problem single shot/frequency Q coarse (κ i )u [l] (κ i ) = b [l] Embarrassingly Parallel Compute SVD of c [l] out and c [l] in Jörn Zimmerling (TU Delft) Model reduction of wave propagation November 8, / 37
28 Computational Scheme - Fine Grid - GPU Compute SVD of c [l] out and c [l] in Evaluate ROM single Frequency s e R m;eik = V m;eik (s e ) H Q fine V m;eik (s e ) F c,m = B H s V m;eik (s e )R 1 m;eik V m;eik(s e ) H B s ROM Evaluation Phase Compute inverse Fourier Transform ˆF m;c (t) = F 1 F m;c (iω) Jörn Zimmerling (TU Delft) Model reduction of wave propagation November 8, / 37
29 Numerical Experiments - Grid Coarsening Coarsen the computation mesh by factor 16 (4 in each direction) Interpolation points until 5.5 ppw (m = 40, N src = 12, M = 100) FD error averaged over all traces RKS m=40 stepsize=0.5 uniform shifts FDFD m=500 stepsize=0.6 PP-RKS m=40 stepsize=4.0 shifts used in PP-RKS normalized frequency (i) Relative error err average x-direction [m] y-direction [m] (j) Configuration Speed [m/s] Figure: Smooth Marmousi test configuration with grid coarsening. Jörn Zimmerling (TU Delft) Model reduction of wave propagation November 8, / 37
30 Numerical Experiments - Grid Coarsening Projection of fine operator gauges the ROM Direct evaluation of coarse operator not accurate FD error averaged over all traces Direct evaluation of 500 frequency points PPRKS with m=40 shifts points per wavelength (k) err average ROM;coarse versus erraverage FD;coarse Figure: Smooth Marmousi test configuration with grid coarsening. Jörn Zimmerling (TU Delft) Model reduction of wave propagation November 8, / 37
31 Numerical Experiments - Grid Coarsening Ricker Wavelet with cut-off frequency of 2.7 ppw on coarse grid 10-6 response [a.u.] Comparison ROM normalized time (l) Time domain trace from the left most source to the right most receiver after m = 40 interpolation points. Jörn Zimmerling (TU Delft) Model reduction of wave propagation November 8, / 37
32 Numerical Experiment III Reiceiver Source PML Resonant Borehole in smooth Geology Resonant behavior causes long runtimes 6 Surface- and 8 BH source-receiver pairs x-direction [m] wavespeed [m/s] y-direction [m] 1500 (m) Simulated configuration. Jörn Zimmerling (TU Delft) Model reduction of wave propagation November 8, / 37
33 Numerical Experiments III x-direction [m] Reiceiver Source PML Speed [m/s] x-direction [m] Reiceiver Source PML Speed [m/s] u [l] (s j ) = g(s j T [l] eik )c[l] out;eik (s j) + g( s j T [l] eik )c[l] in;eik (s j), u [l] (s j ) = g(s j T [l] eik;cm )c[l] out;cm (s j) + g( s j T [l] eik;cm )c[l] in;cm (s j) y-direction [m] (n) Isosurfaces T eik y-direction [m] (o) Isosurf. T eik;cm. m = 40, N src = 14, M SVD = 30 c in/out;eik/cm Jörn Zimmerling (TU Delft) Model reduction of wave propagation November 8, / 37
34 Numerical Experiments III Ricker Wavelet with cut-off frequency of 2.7 ppw on coarse grid response [a.u.] Comparison ROM normalized time (p) Time-domain trace of the coinciding source receiver pair number 1 after m = 40 interpolation points, together with the comparison solution. Jörn Zimmerling (TU Delft) Model reduction of wave propagation November 8, / 37
35 Numerical Experiments III Ricker Wavelet with cut-off frequency of 2.7 ppw on coarse grid response [a.u.] normalized time (q) Time-domain trace from source number 7 inside the borehole to the rightmost surface receiver number 14 after m = 40 interpolation points. Jörn Zimmerling (TU Delft) Model reduction of wave propagation November 8, / 37
36 Conclusions All three challenges (Grid size, Nr of Sources, Nr of interpolation points) can be reduced with phase preconditioning Projection on frequency dependent basis allows ROM beyond the Nyquist limit Can be used for other oscillatory PDEs that have asymptotic solutions Work shifted from solvers to inner products Significantly compressed the ROM into coarse amplitudes Jörn Zimmerling (TU Delft) Model reduction of wave propagation November 8, / 37
37 Paper V. Druskin, R. Remis, M. Zaslavsky and J. Zimmerling, Compressing Large-Scale Wave Propagation Models via Phase-Preconditioned Rational Krylov Subspaces, arxiv: Thanks 2 2 STW (project 14222, Good Vibrations) and Schlumberger Doll-Research Jörn Zimmerling (TU Delft) Model reduction of wave propagation November 8, / 37
38 Numerical Experiments IIII Coarsen the computation mesh by factor 16 (4 in each direction) Interpolation points until 5.5 ppw (m = 40, N src = 12, M = 150) FD error averaged over all traces RKS m=40 step size=1.0 FDFD m=500 step size=4.0 FDFD m=500 step size=1.2 PP-RKS m=40 step size= normalized frequency (g) Relative error err average x-direction [m] Receiver Source PML y-direction [m] (h) Configuration Figure: Marmousi test configuration with grid coarsening Wavespeed [m/s] Jörn Zimmerling (TU Delft) Model reduction of wave propagation November 8, / 4
39 Numerical Experiments IIII Ricker Wavelet with cut-off frequency of 2.7 ppw on coarse grid response [a.u.] Comparison ROM normalized time (a) Time domain trace from the left most source to the right most receiver after m = 40 interpolation points. Jörn Zimmerling (TU Delft) Model reduction of wave propagation November 8, / 4
40 Computational Complexity Cost of the basis computation and evaluation of the ROM Basis Computation on CPU 3, Evaluation GPU 4 Basis Computation comparison Computation Time Block solve fine grid Q fine (s i ) 1 B 10.3s Single solve fine grid Q fine (s i ) 1 b 4.1s Block solve coarse grid Q coarse (s i ) 1 B 0.6s Single solve coarse grid Q coarse (s i ) 1 b 0.2s Evaluation Step Computation Time Scaling Computing phase functions exp(iωt eik ) s N src N f Hadamard Products exp(iωt eik ) c SVD s M SVD N src N f Galerkin inner product V m;eik (s e ) H Q fine V m;eik (s e ) 1.752s N f MSVD 2 N2 src 3 Solved using UMFPACK v on a 4-Core Intel i CPU@3.40 GHz with parallel BLAS level-3 routines 4 Double precision python implementation on an Nvidia GTX 1080 Ti Jörn Zimmerling (TU Delft) Model reduction of wave propagation November 8, / 4
41 Dispersion Correction At 5.5 ppw with a second order scheme we dispersion analytical travel time does not correspond to numerical use ν [l] in decomposition to cancel highest s 2 term ( exp 2sT [l] eik ) k i=1 ( ) D xi exp st [l] eik 2 = s2 (10) ν [l]2 Jörn Zimmerling (TU Delft) Model reduction of wave propagation November 8, / 4
Efficient Wavefield Simulators Based on Krylov Model-Order Reduction Techniques
1 Efficient Wavefield Simulators Based on Krylov Model-Order Reduction Techniques From Resonators to Open Domains Rob Remis Delft University of Technology November 3, 2017 ICERM Brown University 2 Thanks
More informationNonlinear seismic imaging via reduced order model backprojection
Nonlinear seismic imaging via reduced order model backprojection Alexander V. Mamonov, Vladimir Druskin 2 and Mikhail Zaslavsky 2 University of Houston, 2 Schlumberger-Doll Research Center Mamonov, Druskin,
More informationNonlinear acoustic imaging via reduced order model backprojection
Nonlinear acoustic imaging via reduced order model backprojection Alexander V. Mamonov 1, Vladimir Druskin 2, Andrew Thaler 3 and Mikhail Zaslavsky 2 1 University of Houston, 2 Schlumberger-Doll Research
More informationSeismic imaging and multiple removal via model order reduction
Seismic imaging and multiple removal via model order reduction Alexander V. Mamonov 1, Liliana Borcea 2, Vladimir Druskin 3, and Mikhail Zaslavsky 3 1 University of Houston, 2 University of Michigan Ann
More informationarxiv: v1 [math.na] 1 Apr 2015
Nonlinear seismic imaging via reduced order model backprojection Alexander V. Mamonov, University of Houston; Vladimir Druskin and Mikhail Zaslavsky, Schlumberger arxiv:1504.00094v1 [math.na] 1 Apr 2015
More informationSolving Symmetric Indefinite Systems with Symmetric Positive Definite Preconditioners
Solving Symmetric Indefinite Systems with Symmetric Positive Definite Preconditioners Eugene Vecharynski 1 Andrew Knyazev 2 1 Department of Computer Science and Engineering University of Minnesota 2 Department
More informationIterative Helmholtz Solvers
Iterative Helmholtz Solvers Scalability of deflation-based Helmholtz solvers Delft University of Technology Vandana Dwarka March 28th, 2018 Vandana Dwarka (TU Delft) 15th Copper Mountain Conference 2018
More informationAdaptive algebraic multigrid methods in lattice computations
Adaptive algebraic multigrid methods in lattice computations Karsten Kahl Bergische Universität Wuppertal January 8, 2009 Acknowledgements Matthias Bolten, University of Wuppertal Achi Brandt, Weizmann
More informationMultigrid absolute value preconditioning
Multigrid absolute value preconditioning Eugene Vecharynski 1 Andrew Knyazev 2 (speaker) 1 Department of Computer Science and Engineering University of Minnesota 2 Department of Mathematical and Statistical
More informationPreconditioned Locally Minimal Residual Method for Computing Interior Eigenpairs of Symmetric Operators
Preconditioned Locally Minimal Residual Method for Computing Interior Eigenpairs of Symmetric Operators Eugene Vecharynski 1 Andrew Knyazev 2 1 Department of Computer Science and Engineering University
More informationA decade of fast and robust Helmholtz solvers
A decade of fast and robust Helmholtz solvers Werkgemeenschap Scientific Computing Spring meeting Kees Vuik May 11th, 212 1 Delft University of Technology Contents Introduction Preconditioning (22-28)
More informationSpectral analysis of complex shifted-laplace preconditioners for the Helmholtz equation
Spectral analysis of complex shifted-laplace preconditioners for the Helmholtz equation C. Vuik, Y.A. Erlangga, M.B. van Gijzen, and C.W. Oosterlee Delft Institute of Applied Mathematics c.vuik@tudelft.nl
More informationMulti-Domain Approaches for the Solution of High-Frequency Time-Harmonic Propagation Problems
Académie universitaire Wallonie Europe Université de Liège Faculté des Sciences Appliquées Collège de doctorat en Électricité, électronique et informatique Multi-Domain Approaches for the Solution of High-Frequency
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 informationOn the choice of abstract projection vectors for second level preconditioners
On the choice of abstract projection vectors for second level preconditioners C. Vuik 1, J.M. Tang 1, and R. Nabben 2 1 Delft University of Technology 2 Technische Universität Berlin Institut für Mathematik
More informationAlgebraic Multigrid as Solvers and as Preconditioner
Ò Algebraic Multigrid as Solvers and as Preconditioner Domenico Lahaye domenico.lahaye@cs.kuleuven.ac.be http://www.cs.kuleuven.ac.be/ domenico/ Department of Computer Science Katholieke Universiteit Leuven
More informationA PRECONDITIONER FOR THE HELMHOLTZ EQUATION WITH PERFECTLY MATCHED LAYER
European Conference on Computational Fluid Dynamics ECCOMAS CFD 2006 P. Wesseling, E. Oñate and J. Périaux (Eds) c TU Delft, The Netherlands, 2006 A PRECONDITIONER FOR THE HELMHOLTZ EQUATION WITH PERFECTLY
More informationA Jacobi Davidson Method with a Multigrid Solver for the Hermitian Wilson-Dirac Operator
A Jacobi Davidson Method with a Multigrid Solver for the Hermitian Wilson-Dirac Operator Artur Strebel Bergische Universität Wuppertal August 3, 2016 Joint Work This project is joint work with: Gunnar
More informationOverview Avoiding Cycle Skipping: Model Extension MSWI: Space-time Extension Numerical Examples
Overview Avoiding Cycle Skipping: Model Extension MSWI: Space-time Extension Numerical Examples Guanghui Huang Education University of Chinese Academy of Sciences, Beijing, China Ph.D. in Computational
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 informationOn complex shifted Laplace preconditioners for the vector Helmholtz equation
On complex shifted Laplace preconditioners for the vector Helmholtz equation C. Vuik, Y.A. Erlangga, M.B. van Gijzen, C.W. Oosterlee, D. van der Heul Delft Institute of Applied Mathematics c.vuik@tudelft.nl
More informationAN HELMHOLTZ ITERATIVE SOLVER FOR THE THREE-DIMENSIONAL SEISMIC IMAGING PROBLEMS?
European Conference on Computational Fluid Dynamics ECCOMAS CFD 2006 P. Wesseling, E. Oñate, J. Périaux (Eds) c TU Delft, Delft The Netherlands, 2006 AN HELMHOLTZ ITERATIVE SOLVER FOR THE THREE-DIMENSIONAL
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 informationarxiv: v1 [math.na] 6 Nov 2017
Efficient boundary corrected Strang splitting Lukas Einkemmer Martina Moccaldi Alexander Ostermann arxiv:1711.02193v1 [math.na] 6 Nov 2017 Version of November 6, 2017 Abstract Strang splitting is a well
More informationApproximate Methods for Time-Reversal Processing of Large Seismic Reflection Data Sets. University of California, LLNL
147th Annual Meeting: Acoustical Society of America New York City, May 25, 2004 Approximate Methods for Time-Reversal Processing of Large Seismic Reflection Data Sets Speaker: James G. Berryman University
More informationA PML absorbing boundary condition for 2D viscoacoustic wave equation in time domain: modeling and imaging
A PML absorbing boundary condition for 2D viscoacoustic wave equation in time domain: modeling and imaging Ali Fathalian and Kris Innanen Department of Geoscience, University of Calgary Summary The constant-q
More informationAcoustic Wave Equation
Acoustic Wave Equation Sjoerd de Ridder (most of the slides) & Biondo Biondi January 16 th 2011 Table of Topics Basic Acoustic Equations Wave Equation Finite Differences Finite Difference Solution Pseudospectral
More informationStabilization and Acceleration of Algebraic Multigrid Method
Stabilization and Acceleration of Algebraic Multigrid Method Recursive Projection Algorithm A. Jemcov J.P. Maruszewski Fluent Inc. October 24, 2006 Outline 1 Need for Algorithm Stabilization and Acceleration
More informationSolving PDEs with Multigrid Methods p.1
Solving PDEs with Multigrid Methods Scott MacLachlan maclachl@colorado.edu Department of Applied Mathematics, University of Colorado at Boulder Solving PDEs with Multigrid Methods p.1 Support and Collaboration
More informationA fast method for the solution of the Helmholtz equation
A fast method for the solution of the Helmholtz equation Eldad Haber and Scott MacLachlan September 15, 2010 Abstract In this paper, we consider the numerical solution of the Helmholtz equation, arising
More informationPseudospectral and High-Order Time-Domain Forward Solvers
Pseudospectral and High-Order Time-Domain Forward Solvers Qing H. Liu G. Zhao, T. Xiao, Y. Zeng Department of Electrical and Computer Engineering Duke University DARPA/ARO MURI Review, August 15, 2003
More informationFourier transforms, Generalised functions and Greens functions
Fourier transforms, Generalised functions and Greens functions T. Johnson 2015-01-23 Electromagnetic Processes In Dispersive Media, Lecture 2 - T. Johnson 1 Motivation A big part of this course concerns
More informationAn iterative solver enhanced with extrapolation for steady-state high-frequency Maxwell problems
An iterative solver enhanced with extrapolation for steady-state high-frequency Maxwell problems K. Hertel 1,2, S. Yan 1,2, C. Pflaum 1,2, R. Mittra 3 kai.hertel@fau.de 1 Lehrstuhl für Systemsimulation
More informationWe G Model Reduction Approaches for Solution of Wave Equations for Multiple Frequencies
We G15 5 Moel Reuction Approaches for Solution of Wave Equations for Multiple Frequencies M.Y. Zaslavsky (Schlumberger-Doll Research Center), R.F. Remis* (Delft University) & V.L. Druskin (Schlumberger-Doll
More informationKasetsart University Workshop. Multigrid methods: An introduction
Kasetsart University Workshop Multigrid methods: An introduction Dr. Anand Pardhanani Mathematics Department Earlham College Richmond, Indiana USA pardhan@earlham.edu A copy of these slides is available
More informationDistributed Sensor Imaging
Distributed Sensor Imaging George Papanicolaou Department of Mathematics - Stanford University URL: http://georgep.stanford.edu In collaboration with: Gregoire Derveaux, INRIA, France Chrysoula Tsogka,
More informationMarchenko redatuming: advantages and limitations in complex media Summary Incomplete time reversal and one-sided focusing Introduction
Marchenko redatuming: advantages and limitations in complex media Ivan Vasconcelos*, Schlumberger Gould Research, Dirk-Jan van Manen, ETH Zurich, Matteo Ravasi, University of Edinburgh, Kees Wapenaar and
More information1. Fast Iterative Solvers of SLE
1. Fast Iterative Solvers of crucial drawback of solvers discussed so far: they become slower if we discretize more accurate! now: look for possible remedies relaxation: explicit application of the multigrid
More informationAn improved two-grid preconditioner for the solution of three-dimensional Helmholtz problems in heterogeneous media
An improved two-grid preconditioner for the solution of three-dimensional Helmholtz problems in heterogeneous media Henri Calandra, Serge Gratton, Xavier Pinel and Xavier Vasseur Technical Report TR/PA/12/2
More informationComputational Linear Algebra
Computational Linear Algebra PD Dr. rer. nat. habil. Ralf-Peter Mundani Computation in Engineering / BGU Scientific Computing in Computer Science / INF Winter Term 2018/19 Part 4: Iterative Methods PD
More informationAdaptive Multigrid for QCD. Lattice University of Regensburg
Lattice 2007 University of Regensburg Michael Clark Boston University with J. Brannick, R. Brower, J. Osborn and C. Rebbi -1- Lattice 2007, University of Regensburg Talk Outline Introduction to Multigrid
More informationAn Efficient Two-Level Preconditioner for Multi-Frequency Wave Propagation Problems
An Efficient Two-Level Preconditioner for Multi-Frequency Wave Propagation Problems M. Baumann, and M.B. van Gijzen Email: M.M.Baumann@tudelft.nl Delft Institute of Applied Mathematics Delft University
More informationComputational Electromagnetics Definitions, applications and research
Computational Electromagnetics Definitions, applications and research Luis E. Tobón Pontificia Universidad Javeriana Seminario de investigación Departamento de Electrónica y Ciencias de la Computación
More informationMultigrid Methods and their application in CFD
Multigrid Methods and their application in CFD Michael Wurst TU München 16.06.2009 1 Multigrid Methods Definition Multigrid (MG) methods in numerical analysis are a group of algorithms for solving differential
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 informationMigration with Implicit Solvers for the Time-harmonic Helmholtz
Migration with Implicit Solvers for the Time-harmonic Helmholtz Yogi A. Erlangga, Felix J. Herrmann Seismic Laboratory for Imaging and Modeling, The University of British Columbia {yerlangga,fherrmann}@eos.ubc.ca
More informationq is large, to ensure accuracy of the spatial discretization 1 In practice, q = 6 is a good choice.
Time-stepping beyond CFL: a locally one-dimensional scheme for acoustic wave propagation Leonardo Zepeda-Núñez 1*), Russell J. Hewett 1, Minghua Michel Rao 2, and Laurent Demanet 1 1 Dept. of Mathematics
More informationKrylov Subspace Methods for Nonlinear Model Reduction
MAX PLANCK INSTITUT Conference in honour of Nancy Nichols 70th birthday Reading, 2 3 July 2012 Krylov Subspace Methods for Nonlinear Model Reduction Peter Benner and Tobias Breiten Max Planck Institute
More informationNonlinear Iterative Solution of the Neutron Transport Equation
Nonlinear Iterative Solution of the Neutron Transport Equation Emiliano Masiello Commissariat à l Energie Atomique de Saclay /DANS//SERMA/LTSD emiliano.masiello@cea.fr 1/37 Outline - motivations and framework
More informationSPARSE SOLVERS POISSON EQUATION. Margreet Nool. November 9, 2015 FOR THE. CWI, Multiscale Dynamics
SPARSE SOLVERS FOR THE POISSON EQUATION Margreet Nool CWI, Multiscale Dynamics November 9, 2015 OUTLINE OF THIS TALK 1 FISHPACK, LAPACK, PARDISO 2 SYSTEM OVERVIEW OF CARTESIUS 3 POISSON EQUATION 4 SOLVERS
More informationFDTD for 1D wave equation. Equation: 2 H Notations: o o. discretization. ( t) ( x) i i i i i
FDTD for 1D wave equation Equation: 2 H = t 2 c2 2 H x 2 Notations: o t = nδδ, x = iδx o n H nδδ, iδx = H i o n E nδδ, iδx = E i discretization H 2H + H H 2H + H n+ 1 n n 1 n n n i i i 2 i+ 1 i i 1 = c
More informationA Fast Fourier transform based direct solver for the Helmholtz problem. AANMPDE 2017 Palaiochora, Crete, Oct 2, 2017
A Fast Fourier transform based direct solver for the Helmholtz problem Jari Toivanen Monika Wolfmayr AANMPDE 2017 Palaiochora, Crete, Oct 2, 2017 Content 1 Model problem 2 Discretization 3 Fast solver
More informationAMS526: Numerical Analysis I (Numerical Linear Algebra)
AMS526: Numerical Analysis I (Numerical Linear Algebra) Lecture 24: Preconditioning and Multigrid Solver Xiangmin Jiao SUNY Stony Brook Xiangmin Jiao Numerical Analysis I 1 / 5 Preconditioning Motivation:
More informationGEOPHYSICAL INVERSE THEORY AND REGULARIZATION PROBLEMS
Methods in Geochemistry and Geophysics, 36 GEOPHYSICAL INVERSE THEORY AND REGULARIZATION PROBLEMS Michael S. ZHDANOV University of Utah Salt Lake City UTAH, U.S.A. 2OO2 ELSEVIER Amsterdam - Boston - London
More information- Part 4 - Multicore and Manycore Technology: Chances and Challenges. Vincent Heuveline
- Part 4 - Multicore and Manycore Technology: Chances and Challenges Vincent Heuveline 1 Numerical Simulation of Tropical Cyclones Goal oriented adaptivity for tropical cyclones ~10⁴km ~1500km ~100km 2
More informationAttenuation compensation in least-squares reverse time migration using the visco-acoustic wave equation
Attenuation compensation in least-squares reverse time migration using the visco-acoustic wave equation Gaurav Dutta, Kai Lu, Xin Wang and Gerard T. Schuster, King Abdullah University of Science and Technology
More informationData-to-Born transform for inversion and imaging with waves
Data-to-Born transform for inversion and imaging with waves Alexander V. Mamonov 1, Liliana Borcea 2, Vladimir Druskin 3, and Mikhail Zaslavsky 3 1 University of Houston, 2 University of Michigan Ann Arbor,
More informationKarhunen-Loève Approximation of Random Fields Using Hierarchical Matrix Techniques
Institut für Numerische Mathematik und Optimierung Karhunen-Loève Approximation of Random Fields Using Hierarchical Matrix Techniques Oliver Ernst Computational Methods with Applications Harrachov, CR,
More informationA Comparison of Adaptive Algebraic Multigrid and Lüscher s Inexact Deflation
A Comparison of Adaptive Algebraic Multigrid and Lüscher s Inexact Deflation Andreas Frommer, Karsten Kahl, Stefan Krieg, Björn Leder and Matthias Rottmann Bergische Universität Wuppertal April 11, 2013
More informationAn Efficient Low Memory Implicit DG Algorithm for Time Dependent Problems
An Efficient Low Memory Implicit DG Algorithm for Time Dependent Problems P.-O. Persson and J. Peraire Massachusetts Institute of Technology 2006 AIAA Aerospace Sciences Meeting, Reno, Nevada January 9,
More informationProjected Schur Complement Method for Solving Non-Symmetric Saddle-Point Systems (Arising from Fictitious Domain Approaches)
Projected Schur Complement Method for Solving Non-Symmetric Saddle-Point Systems (Arising from Fictitious Domain Approaches) Jaroslav Haslinger, Charles University, Prague Tomáš Kozubek, VŠB TU Ostrava
More informationSome Geometric and Algebraic Aspects of Domain Decomposition Methods
Some Geometric and Algebraic Aspects of Domain Decomposition Methods D.S.Butyugin 1, Y.L.Gurieva 1, V.P.Ilin 1,2, and D.V.Perevozkin 1 Abstract Some geometric and algebraic aspects of various domain decomposition
More informationMultilevel low-rank approximation preconditioners Yousef Saad Department of Computer Science and Engineering University of Minnesota
Multilevel low-rank approximation preconditioners Yousef Saad Department of Computer Science and Engineering University of Minnesota SIAM CSE Boston - March 1, 2013 First: Joint work with Ruipeng Li Work
More informationSolving the Stochastic Steady-State Diffusion Problem Using Multigrid
Solving the Stochastic Steady-State Diffusion Problem Using Multigrid Tengfei Su Applied Mathematics and Scientific Computing Advisor: Howard Elman Department of Computer Science Sept. 29, 2015 Tengfei
More informationAspects of Multigrid
Aspects of Multigrid Kees Oosterlee 1,2 1 Delft University of Technology, Delft. 2 CWI, Center for Mathematics and Computer Science, Amsterdam, SIAM Chapter Workshop Day, May 30th 2018 C.W.Oosterlee (CWI)
More informationNested Krylov methods for shifted linear systems
Nested Krylov methods for shifted linear systems M. Baumann, and M. B. van Gizen Email: M.M.Baumann@tudelft.nl Delft Institute of Applied Mathematics Delft University of Technology Delft, The Netherlands
More informationExponentially Convergent Sparse Discretizations and Application to Near Surface Geophysics
Exponentially Convergent Sparse Discretizations and Application to Near Surface Geophysics Murthy N. Guddati North Carolina State University November 9, 017 Outline Part 1: Impedance Preserving Discretization
More informationAcoustics Analysis of Speaker ANSYS, Inc. November 28, 2014
Acoustics Analysis of Speaker 1 Introduction ANSYS 14.0 offers many enhancements in the area of acoustics. In this presentation, an example speaker analysis will be shown to highlight some of the acoustics
More informationConstrained Minimization and Multigrid
Constrained Minimization and Multigrid C. Gräser (FU Berlin), R. Kornhuber (FU Berlin), and O. Sander (FU Berlin) Workshop on PDE Constrained Optimization Hamburg, March 27-29, 2008 Matheon Outline Successive
More informationStability of Krylov Subspace Spectral Methods
Stability of Krylov Subspace Spectral Methods James V. Lambers Department of Energy Resources Engineering Stanford University includes joint work with Patrick Guidotti and Knut Sølna, UC Irvine Margot
More informationPreface to the Second Edition. Preface to the First Edition
n page v Preface to the Second Edition Preface to the First Edition xiii xvii 1 Background in Linear Algebra 1 1.1 Matrices................................. 1 1.2 Square Matrices and Eigenvalues....................
More informationExponential multistep methods of Adams-type
Exponential multistep methods of Adams-type Marlis Hochbruck and Alexander Ostermann KARLSRUHE INSTITUTE OF TECHNOLOGY (KIT) 0 KIT University of the State of Baden-Wuerttemberg and National Laboratory
More informationRobust solution of Poisson-like problems with aggregation-based AMG
Robust solution of Poisson-like problems with aggregation-based AMG Yvan Notay Université Libre de Bruxelles Service de Métrologie Nucléaire Paris, January 26, 215 Supported by the Belgian FNRS http://homepages.ulb.ac.be/
More informationPreconditioned space-time boundary element methods for the heat equation
W I S S E N T E C H N I K L E I D E N S C H A F T Preconditioned space-time boundary element methods for the heat equation S. Dohr and O. Steinbach Institut für Numerische Mathematik Space-Time Methods
More 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 informationPartial Differential Equations and the Finite Element Method
Partial Differential Equations and the Finite Element Method Pavel Solin The University of Texas at El Paso Academy of Sciences ofthe Czech Republic iwiley- INTERSCIENCE A JOHN WILEY & SONS, INC, PUBLICATION
More informationIndefinite and physics-based preconditioning
Indefinite and physics-based preconditioning Jed Brown VAW, ETH Zürich 2009-01-29 Newton iteration Standard form of a nonlinear system F (u) 0 Iteration Solve: Update: J(ũ)u F (ũ) ũ + ũ + u Example (p-bratu)
More informationLowrank RTM for converted wave imaging
Lowrank RM for converted wave imaging Lorenzo Casasanta * (CGG), Zhiguang Xue (U Austin), Sam Gray (CGG) Summary his paper is an attempt to fill the technology gap existing between pure P- and PS-wave
More informationDomain Decomposition Preconditioners for Spectral Nédélec Elements in Two and Three Dimensions
Domain Decomposition Preconditioners for Spectral Nédélec Elements in Two and Three Dimensions Bernhard Hientzsch Courant Institute of Mathematical Sciences, New York University, 51 Mercer Street, New
More informationAcceleration of Time Integration
Acceleration of Time Integration edited version, with extra images removed Rick Archibald, John Drake, Kate Evans, Doug Kothe, Trey White, Pat Worley Research sponsored by the Laboratory Directed Research
More informationOverview of sparse system identification
Overview of sparse system identification J.-Ch. Loiseau 1 & Others 2, 3 1 Laboratoire DynFluid, Arts et Métiers ParisTech, France 2 LIMSI, Université d Orsay CNRS, France 3 University of Washington, Seattle,
More informationImplicit Solution of Viscous Aerodynamic Flows using the Discontinuous Galerkin Method
Implicit Solution of Viscous Aerodynamic Flows using the Discontinuous Galerkin Method Per-Olof Persson and Jaime Peraire Massachusetts Institute of Technology 7th World Congress on Computational Mechanics
More informationThe failure of the sounding assumption in electroseismic investigations
The failure of the sounding assumption in electroseismic investigations F.D. Fourie, J.F. Botha Institute for Groundwater Studies, University of the Free State, PO Box 339, Bloemfontein, 9300, South Africa
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 informationPath integrals and the classical approximation 1 D. E. Soper 2 University of Oregon 14 November 2011
Path integrals and the classical approximation D. E. Soper University of Oregon 4 November 0 I offer here some background for Sections.5 and.6 of J. J. Sakurai, Modern Quantum Mechanics. Introduction There
More informationOUTLINE ffl CFD: elliptic pde's! Ax = b ffl Basic iterative methods ffl Krylov subspace methods ffl Preconditioning techniques: Iterative methods ILU
Preconditioning Techniques for Solving Large Sparse Linear Systems Arnold Reusken Institut für Geometrie und Praktische Mathematik RWTH-Aachen OUTLINE ffl CFD: elliptic pde's! Ax = b ffl Basic iterative
More informationFrequency-domain ray series for viscoelastic waves with a non-symmetric stiffness matrix
Frequency-domain ray series for viscoelastic waves with a non-symmetric stiffness matrix Ludě Klimeš Department of Geophysics, Faculty of Mathematics Physics, Charles University, Ke Karlovu 3, 121 16 Praha
More informationElastic wave-equation migration for laterally varying isotropic and HTI media. Richard A. Bale and Gary F. Margrave
Elastic wave-equation migration for laterally varying isotropic and HTI media Richard A. Bale and Gary F. Margrave a Outline Introduction Theory Elastic wavefield extrapolation Extension to laterally heterogeneous
More informationA reduced basis method and ROM-based optimization for batch chromatography
MAX PLANCK INSTITUTE MedRed, 2013 Magdeburg, Dec 11-13 A reduced basis method and ROM-based optimization for batch chromatography Yongjin Zhang joint work with Peter Benner, Lihong Feng and Suzhou Li Max
More informationScalable Domain Decomposition Preconditioners For Heterogeneous Elliptic Problems
Scalable Domain Decomposition Preconditioners For Heterogeneous Elliptic Problems Pierre Jolivet, F. Hecht, F. Nataf, C. Prud homme Laboratoire Jacques-Louis Lions Laboratoire Jean Kuntzmann INRIA Rocquencourt
More informationParameterized Partial Differential Equations and the Proper Orthogonal D
Parameterized Partial Differential Equations and the Proper Orthogonal Decomposition Stanford University February 04, 2014 Outline Parameterized PDEs The steady case Dimensionality reduction Proper orthogonal
More informationThe mathematics that models wavefield physics in engineering applications - A voyage through the landscape of fundamentals by Adrianus T.
The mathematics that models wavefield physics in engineering applications - A voyage through the landscape of fundamentals by Adrianus T. de Hoop Delft University of Technology Laboratory of Electromagnetic
More informationIndex. C 2 ( ), 447 C k [a,b], 37 C0 ( ), 618 ( ), 447 CD 2 CN 2
Index advection equation, 29 in three dimensions, 446 advection-diffusion equation, 31 aluminum, 200 angle between two vectors, 58 area integral, 439 automatic step control, 119 back substitution, 604
More informationMultilevel Low-Rank Preconditioners Yousef Saad Department of Computer Science and Engineering University of Minnesota. Modelling 2014 June 2,2014
Multilevel Low-Rank Preconditioners Yousef Saad Department of Computer Science and Engineering University of Minnesota Modelling 24 June 2,24 Dedicated to Owe Axelsson at the occasion of his 8th birthday
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 informationFinite Element Solver for Flux-Source Equations
Finite Element Solver for Flux-Source Equations Weston B. Lowrie A thesis submitted in partial fulfillment of the requirements for the degree of Master of Science in Aeronautics Astronautics University
More informationAMS526: Numerical Analysis I (Numerical Linear Algebra for Computational and Data Sciences)
AMS526: Numerical Analysis I (Numerical Linear Algebra for Computational and Data Sciences) Lecture 19: Computing the SVD; Sparse Linear Systems Xiangmin Jiao Stony Brook University Xiangmin Jiao Numerical
More informationSeismic inverse scattering by reverse time migration
Seismic inverse scattering by reverse time migration Chris Stolk 1, Tim Op t Root 2, Maarten de Hoop 3 1 University of Amsterdam 2 University of Twente 3 Purdue University MSRI Inverse Problems Seminar,
More informationLecture 8: Fast Linear Solvers (Part 7)
Lecture 8: Fast Linear Solvers (Part 7) 1 Modified Gram-Schmidt Process with Reorthogonalization Test Reorthogonalization If Av k 2 + δ v k+1 2 = Av k 2 to working precision. δ = 10 3 2 Householder Arnoldi
More informationFourier Coefficients of Far Fields and the Convex Scattering Support John Sylvester University of Washington
Fourier Coefficients of Far Fields and the Convex Scattering Support John Sylvester University of Washington Collaborators: Steve Kusiak Roland Potthast Jim Berryman Gratefully acknowledge research support
More information