A robust computational method for the Schrödinger equation cross sections using an MG-Krylov scheme
|
|
- Denis Hodge
- 5 years ago
- Views:
Transcription
1 A robust computational method for the Schrödinger equation cross sections using an MG-Krylov scheme International Conference On Preconditioning Techniques For Scientific And Industrial Applications June 2015, Eindhoven, The Netherlands University of Antwerp, Applied Mathematics S. Cools and W. Vanroose Contact:
2 Overview [1] Motivation [2] The Helmholtz and Schrödinger equations [3] The classical far field map calculation [4] The complex contour approach for the far field map [5] Schrödinger cross sections: the MG-CCCS preconditioner [6] Conclusions & discussion
3 Motivation Simulation of quantum mechanical molecular break-up reactions. Electron wave character CERN Large Hadron Collider Practical use: determine material composition and structure.
4 Schrödinger equation Full 3D break-up problem description two particle system (d = 2) Schrödinger equation (6-dimensional) (H E) ψ(r 1, r 2) = ζ(r 1, r 2). with outgoing wave ψ i.f.o. r 1, r 2 R 3 and Hamiltonian H = 1 2m r 1 1 2m r 2 + V 1(r 1) + V 2(r 2) + V 12(r 1, r 2). Right-hand side source term ζ(r 1, r 2 ) models incoming electron (electron scattering), or [Rescigno Baertschy Isaacs McCurdy 1999] dipole operator (photo-ionization). [Vanroose Martin Rescigno McCurdy 2005]
5 Partial wave expansion Spherical coordinate representation: r 1 = (ρ 1, Θ 1 ), r 2 = (ρ 2, Θ 2 ) l 1 l 2 ψ(r 1, r 2) = ψ l1 m 1,l 2 m 2 (ρ 1, ρ 2)Y l1 m 1 (Θ 1)Y l2 m 2 (Θ 2), l 1 =0 l 2 =0 m 2 = l 2 m 1 = l 1 leads to a coupled system for the radial wave components ψ l1 m 1,l 2 m 2 (ρ 1, ρ 2) ψ = ψ l 1 m 1 (ρ1, ρ2),l 2 m 2. Coupled partial wave system 2m 1 l 1,l 2 + V l1 m 1 l 2 m 2 ;l 1 m 1 l 2 m 2 E V l1 m 1 l 2 m 2 ;l 1 m 1 l 2 m 2 V l m 1 1 l 2 m 2 ;l 1 m 1 l 2 m 2 2m 1 l 1,l + V l 2 1 m 1 l 2 m 2 ;l 1 m 1 l 2 m E 2. Note: system decouples for V 1, V 2 and V 12 spherically symmetric..... ψ = ζ.
6 Relation to Helmholtz equation Diagonal blocks assume form of 2D driven Schrödinger equation ( 12 + V (x) E ) u(x) = g(x), x = (x, y) R 2. Helmholtz equation (d-dimensional) (equivalent formulation) ( k 2 (x) ) u(x) = f (x), x R d, with k 2 (x) = 2(E V (x)) and f (x) = 2g(x). Properties. potential V (x) and hence wavenumber k(x) is analytic, experimental observations (cross sections) are far field maps. [McCurdy Baertschy Rescigno 2004]
7 Helmholtz equation Representation of the physics behind a wave scattering at an object χ defined on a compact area O located within a domain Ω R d. Scattered wave solution u sc (x) satisfies inhomogeneous Helmholtz ( k 2 (x) ) u sc (x) = f (x), x Ω R d, with f (x) = k 2 0 χ(x)u in(x). Aim: calculate far field amplitude map F (α) u in O α
8 Far field map Analytic solution on whole R d using Green s function: u(x ) = G(x, x ) k0 2 χ(x) [u Ω }{{} in (x) + u sc (x)] dx, x R d. }{{} Green s function scattered wave [Colton Kress 1998] Calculate u in any point x R d outside the numerical domain Ω, using only the information inside the numerical domain. Computation: Split the far field integral into a sum I 1 + I 2, with I 1 = G(x, x )χ(x)u in (x)dx and I 2 = G(x, x )χ(x)u sc (x)dx Ω Ω }{{}}{{} all factors known explicitly requires u sc(x) for x Ω
9 State-of-the-art Helmholtz solvers. Preconditioned Krylov methods Solve M 1 Au = M 1 f using Krylov methods, with M = ϱk(x) 2 (preconditioner) where A = k(x) 2 and Mu = f easily solvable iteratively. ϱ = 1: original Helmholtz operator [von Helmholtz 19th century] ϱ = 0: Laplacian [Bayliss Goldstein Turkel 1983] ϱ < 0: shifted Laplacian or screened Poisson operator [Laird 2001] ϱ C: complex shifted Laplacian (CSL): ϱ = α + βi [Erlangga Vuik Oosterlee 2004]
10 Far field map Complex contour approach. I{z} For u and χ analytical the far field integral I 2 = G(x, x )χ(x)u sc (x)dx Ω Z 1 Z 2 γ R{z} can be calculated over a complex contour Z = Z 1 + Z 2, rather than over the real domain Ω, i.e. I 2 = G(x, x )χ(z)u sc (z)dz + Z 1 G(x, x )χ(z)u sc (z)dz. Z 2
11 Far field map Complex contour approach. I{z} For u and χ analytical the far field integral I 2 = G(x, x )χ(x)u sc (x)dx Ω Z 1 Z 2 γ R{z} can be calculated over a complex contour Z = Z 1 + Z 2, rather than over the real domain Ω, i.e. I 2 = G(x, x )χ(z)u sc (z)dz + G(x, x )χ(z)u sc (z)dz. Z } 1 Z {{} 2 requires u sc(z) for z Z 1
12 Helmholtz on complex contour Complex shifted Laplacian (CSL) system with shift parameter β R ( (1 + iβ)k 2 (x) ) u(x) = f (x) is efficiently solvable using multigrid. [Erlangga Oosterlee Vuik 2004] FD discretization: ( 1h 2 L (1 + iβ)k2 ) u h = f h, with Laplacian stencil matrix L. Division by (1 + iβ) yields ( ) 1 (1 + iβ)h L 2 k2 u h = f h 1 + iβ, the original Helmholtz system with complex h = 1 + iβ h = ρe iγ h. [Reps Vanroose bin Zubair 2010]
13 Contour approach - 2D validation Object of interest 5-point FD stencil, n x n y = I{z} Z1 γ Z2 R{z} Real domain with ECS χ(x) (θ ECS = 45 ) Complex contour χ(z) (γ = 14.6 )
14 Contour approach - 2D validation Scattered wave solution 5-point FD stencil, n x n y = I{z} Z1 γ Z2 R{z} Real domain with ECS u(x) LU factorization Complex contour u(z) V(1,1) cycles (tol res = 10 6 )
15 Contour approach - 2D validation Far field amplitude map 5-point FD stencil, n x n y = I{z} Z1 γ Z2 R{z} Real domain with ECS F (α) LU factorization Complex contour F (α) V(1,1) cycles (tol res = 10 6 )
16 Contour approach - 3D validation 3D damped Helmholtz solver (γ = 10 ) n x n y n z k 0 = 1/4 k 0 = 1/2 k 0 = 1 k 0 = 2 k 0 = 4 10 (0.79s.) 9 (4.65s.) 9 (44.2s.) 9 (352s.) 9 (2778s.) (0.92s.) 10 (4.96s.) 10 (48.3s.) 10 (390s.) 9 (2797s.) (0.62s.) 13 (6.59s.) 11 (54.6s.) 10 (387s.) 10 (3079s.) (0.28s.) 8 (4.24s.) 13 (63.9s.) 11 (428s.) 10 (3006s.) (0.20s.) 2 (1.35s.) 7 (36.1s.) 13 (503s.) 11 (3306s.) GMRES(3)-smoothed V(1,1) cycles (tol res = 10 6 ) k 0h = 0.625
17 Contour approach - 3D validation 3D damped Helmholtz solver (γ = 10, k 0 = 1) n x n y n z CPU time 0.20 s s s s. 462 s. r 2 3.3e-5 7.9e-5 2.7e-5 1.1e-5 4.6e s. 1.0e-5 GMRES(3)-smoothed FMG(1,1) cycle [Vasseur 2012] Serial implementation, Intel Core i7-2720qm 2.20GHz CPU, 6MB Cache, 8GB RAM.
18 Schrödinger cross sections Single ionization amplitude = single ionization probability s n(e) = φ kn (x)φ n(y) [g(x, y) V 12(x, y)u(x, y) ] dx dy. Ω }{{}}{{} Green s function scattered wave Double ionization cross section = double ionization probability d k1,k 2 (E) = φ k1 (x)φ k2 (y) [g(x, y) V 12(x, y)u(x, y) ] dx dy. Ω }{{}}{{} Green s function scattered wave Single ionization Double ionization
19 Schrödinger on complex contour 2D driven Schrödinger equation ( 12 + V1(x) + V2(y) + V12(x, y) E ) u(x, y) = f (x, y), with x, y 0, V i potentials, E R energy of the system. Multigrid convergence factor damped Schrödinger on full complex grid with γ 8.5 GMRES(3)-smoothed multigrid V(1,1) cycles Observation: Poor convergence for 1<E <0.
20 Schrödinger spectral analysis Single ionization Double ionization E < 0 E = 0 E > imag 4 imag 4 imag real real real bound states with near-zero e.v. s destroy multigrid convergence
21 Coupled channel approximation Bound states are characterized by eigenstates of 1D Hamiltonians H 1φ n(x) = λ nφ n(x), H 2ϕ n(y) = µ nϕ n(y), with λ n < 0 and µ n < 0, hence for M n x, L n y approximate M L u(x, y) A m(y)φ m(x) + B l (x)ϕ l (y). m=1 l=1 [Heller Reinhardt 1973] [McCarthy Stelbovics 1983] Coupled channel correction step (CCCS). M L u (k+1) (x, y) = u (k) (x, y) + em(y)φ A m(x) + el B (x)ϕ l (y). m=1 l=1 Determine coefficients as solution of M+L 1D Schrödinger systems.
22 Numerical results: convergence Convergence rate of MG-CCCS as stand-alone solver/preconditioner Exponential model problem Temkin-Poet model problem V1 (x) = 4.5 exp( x 2 ) V1 (x) = 1/x V12 (x, y ) = 2 exp( 0.1(x + y )2 ) V12 (x, y ) = 1/ max(x, y ) Ω = [0, 20]2, nx ny = Ω = [0, 100]2, nx ny =
23 Numerical results: cross sections 2D Temkin Poet model problem Potentials: V1 (x) = 1/x, V2 (y ) = 1/y and V12 (x, y ) = 1/ max(x, y ) Discretization: Ω = [0, 108]2 with nx ny = spectral element grid Solver: LU with θecs = 30 (real) vs. MG-CCCS with γ = 9 (complex) Single differential cross section Total cross section E =1 E [0, 2]
24 In this work we presented... Conclusions Proof-of-concept for complex contour approach proposed in [Cools Reps Vanroose 2013] and application to impact scattering. Coupled Channel Correction Step (CCCS) after each MG V-cycle accounts for presence of localized bound states. Fast and robust method for the computation of the ionization cross sections for electron-impact models (for any energy E). MG-CCCS validated on 2D Temkin-Poet model problem, convergence as solver/preconditioner shows O(N) scalability. Outlook Generalization to 3D Schrödinger partial wave systems. Analysis of bound states and influence of complex rotation γ for general discretizations.
25 References Y.A. Erlangga, C.W. Oosterlee, and C. Vuik, A novel multigrid based preconditioner for heterogeneous Helmholtz problems, SIAM Journal on Scientific Computing 27(4), pp , W. Vanroose, D.A. Horner, F. Martin, T.N. Rescigno and C.W. McCurdy. Double photoionization of aligned molecular hydrogen. Physical Review A, 74(5), pp , O.G. Ernst and M.J. Gander, Why it is difficult to solve Helmholtz problems with classical iterative methods, Numerical Analysis of Multiscale Problems, pp , S. Cools, B. Reps and W. Vanroose, An efficient multigrid method calculation of the far field map for Helmholtz problems, SIAM Journal on Scientific Computing 36(3), pp. B367-B395, S. Cools and W. Vanroose, A fast and robust computational method for the ionization cross sections of the driven Schrödinger equation using an O(N) multigrid-based scheme, under review, arxiv: , 2015.
arxiv: v5 [math.na] 19 Mar 2014
AN EFFICIENT MULTIGRID CALCULATION OF THE FAR FIELD MAP FOR HELMHOLTZ AND SCHRÖDINGER EQUATIONS SIEGFRIED COOLS, BRAM REPS, AND WIM VANROOSE arxiv:1211.4461v5 [math.na] 19 Mar 2014 Abstract. In this paper
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 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 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 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 informationON A ROBUST ITERATIVE METHOD FOR HETEROGENEOUS HELMHOLTZ PROBLEMS FOR GEOPHYSICS APPLICATIONS
INTERNATIONAL JOURNAL OF NUMERICAL ANALYSIS AND MODELING Volume 2, Supp, Pages 197 28 c 25 Institute for Scientific Computing and Information ON A ROBUST ITERATIVE METHOD FOR HETEROGENEOUS HELMHOLTZ PROBLEMS
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 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 informationarxiv: v2 [math.na] 29 Oct 2011
arxiv:115.347v2 [math.na] 29 Oct 211 Analyzing the wave number dependency of the convergence rate of a multigrid preconditioned Krylov method for the Helmholtz equation with an absorbing layer B. Reps
More informationT-matrix calculations for the electron-impact ionization of hydrogen in the Temkin-Poet model
T-matrix calculations for the electron-impact ionization of hydrogen in the Temkin-Poet model M. S. Pindzola, D. Mitnik, and F. Robicheaux Department of Physics, Auburn University, Auburn, Alabama 36849
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 informationarxiv: v1 [physics.atom-ph] 10 Jun 2016
Potential splitting approach to e-h and e-he + scattering with zero total angular momentum E. Yarevsky 1, S. L. Yakovlev 1, and N. Elander 2 1 St Petersburg State University, 7/9 Universitetskaya nab.,
More informationA MULTIGRID-BASED SHIFTED-LAPLACIAN PRECONDITIONER FOR A FOURTH-ORDER HELMHOLTZ DISCRETIZATION.
A MULTIGRID-BASED SHIFTED-LAPLACIAN PRECONDITIONER FOR A FOURTH-ORDER HELMHOLTZ DISCRETIZATION. N.UMETANI, S.P.MACLACHLAN, C.W. OOSTERLEE Abstract. In this paper, an iterative solution method for a fourth-order
More informationarxiv: v2 [math.na] 5 Jan 2015
A SADDLE POINT NUMERICAL METHOD FOR HELMHOLT EQUATIONS RUSSELL B. RICHINS arxiv:1301.4435v2 [math.na] 5 Jan 2015 Abstract. In a previous work, the author and D.C. Dobson proposed a numerical method for
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 informationCONVERGENCE BOUNDS FOR PRECONDITIONED GMRES USING ELEMENT-BY-ELEMENT ESTIMATES OF THE FIELD OF VALUES
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 CONVERGENCE BOUNDS FOR PRECONDITIONED GMRES USING ELEMENT-BY-ELEMENT
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 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 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 information(This is a sample cover image for this issue. The actual cover is not yet available at this time.)
(This is a sample cover image for this issue. The actual cover is not yet available at this time.) This article appeared in a journal published by Elsevier. The attached copy is furnished to the author
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 informationThe solution of the discretized incompressible Navier-Stokes equations with iterative methods
The solution of the discretized incompressible Navier-Stokes equations with iterative methods Report 93-54 C. Vuik Technische Universiteit Delft Delft University of Technology Faculteit der Technische
More informationComparison of multigrid and incomplete LU shifted-laplace preconditioners for the inhomogeneous Helmholtz equation
Comparison of multigrid and incomplete LU shifted-laplace preconditioners for the inhomogeneous Helmholtz equation Y.A. Erlangga, C. Vuik, C.W. Oosterlee Delft Institute of Applied Mathematics, Delft University
More informationMultigrid based preconditioners for the numerical solution of two-dimensional heterogeneous problems in geophysics
Technical Report RAL-TR-2007-002 Multigrid based preconditioners for the numerical solution of two-dimensional heterogeneous problems in geophysics I. S. Duff, S. Gratton, X. Pinel, and X. Vasseur January
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 informationScattering theory with the Coulomb potential
Journal of Physics: Conference Series Scattering theory with the Coulomb potential To cite this article: A S Kadyrov et al 2009 J. Phys.: Conf. Ser. 194 012017 View the article online for updates and enhancements.
More informationTime-dependent close-coupling calculations of the triple-differential cross section for electron-impact ionization of hydrogen
PHYSICAL REVIEW A, VOLUME 65, 471 Time-dependent close-coupling calculations of the triple-differential cross section for electron-impact ionization of hydrogen J. Colgan, M. S. Pindzola, and F. J. Robicheaux
More informationFast solvers for steady incompressible flow
ICFD 25 p.1/21 Fast solvers for steady incompressible flow Andy Wathen Oxford University wathen@comlab.ox.ac.uk http://web.comlab.ox.ac.uk/~wathen/ Joint work with: Howard Elman (University of Maryland,
More informationProblem 1: Spin 1 2. particles (10 points)
Problem 1: Spin 1 particles 1 points 1 Consider a system made up of spin 1/ particles. If one measures the spin of the particles, one can only measure spin up or spin down. The general spin state of a
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 informationExponential integrators for semilinear parabolic problems
Exponential integrators for semilinear parabolic problems Marlis Hochbruck Heinrich-Heine University Düsseldorf Germany Innsbruck, October 2004 p. Outline Exponential integrators general class of methods
More informationEigenmodes for coupled harmonic vibrations. Algebraic Method for Harmonic Oscillator.
PHYS208 spring 2008 Eigenmodes for coupled harmonic vibrations. Algebraic Method for Harmonic Oscillator. 07.02.2008 Adapted from the text Light - Atom Interaction PHYS261 autumn 2007 Go to list of topics
More information1 Assignment 1: Nonlinear dynamics (due September
Assignment : Nonlinear dynamics (due September 4, 28). Consider the ordinary differential equation du/dt = cos(u). Sketch the equilibria and indicate by arrows the increase or decrease of the solutions.
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 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 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 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 informationVarious problems in electron-atom collision theory
Various problems in electron-atom collision theory I. Bray, K. Bartschat, D. V. Fursa, A. S. Kadyrov and A. T. Stelbovics Centre for Atomic, Molecular, and Surface Physics, Murdoch University, Perth 615
More informationChapter 5 Fast Multipole Methods
Computational Electromagnetics; Chapter 1 1 Chapter 5 Fast Multipole Methods 5.1 Near-field and far-field expansions Like the panel clustering, the Fast Multipole Method (FMM) is a technique for the fast
More informationLaplace s Equation. Chapter Mean Value Formulas
Chapter 1 Laplace s Equation Let be an open set in R n. A function u C 2 () is called harmonic in if it satisfies Laplace s equation n (1.1) u := D ii u = 0 in. i=1 A function u C 2 () is called subharmonic
More informationMultipole-Based Preconditioners for Sparse Linear Systems.
Multipole-Based Preconditioners for Sparse Linear Systems. Ananth Grama Purdue University. Supported by the National Science Foundation. Overview Summary of Contributions Generalized Stokes Problem Solenoidal
More informationColumbia University Department of Physics QUALIFYING EXAMINATION
Columbia University Department of Physics QUALIFYING EXAMINATION Wednesday, January 13, 2016 3:10PM to 5:10PM Modern Physics Section 4. Relativity and Applied Quantum Mechanics Two hours are permitted
More informationThe implementation of the Helmholtz problem on a Maxeler machine
Delft University of Technology Faculty of Electrical Engineering, Mathematics and Computer Science Delft Institute of Applied Mathematics The implementation of the Helmholtz problem on a Maxeler machine
More informationM.A. Botchev. September 5, 2014
Rome-Moscow school of Matrix Methods and Applied Linear Algebra 2014 A short introduction to Krylov subspaces for linear systems, matrix functions and inexact Newton methods. Plan and exercises. M.A. Botchev
More informationRBF-FD Approximation to Solve Poisson Equation in 3D
RBF-FD Approximation to Solve Poisson Equation in 3D Jagadeeswaran.R March 14, 2014 1 / 28 Overview Problem Setup Generalized finite difference method. Uses numerical differentiations generated by Gaussian
More informationSparse Matrix Techniques for MCAO
Sparse Matrix Techniques for MCAO Luc Gilles lgilles@mtu.edu Michigan Technological University, ECE Department Brent Ellerbroek bellerbroek@gemini.edu Gemini Observatory Curt Vogel vogel@math.montana.edu
More informationEdwin van der Weide and Magnus Svärd. I. Background information for the SBP-SAT scheme
Edwin van der Weide and Magnus Svärd I. Background information for the SBP-SAT scheme As is well-known, stability of a numerical scheme is a key property for a robust and accurate numerical solution. Proving
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 informationChemistry 532 Problem Set 7 Spring 2012 Solutions
Chemistry 53 Problem Set 7 Spring 01 Solutions 1. The study of the time-independent Schrödinger equation for a one-dimensional particle subject to the potential function leads to the differential equation
More informationFast numerical methods for solving linear PDEs
Fast numerical methods for solving linear PDEs P.G. Martinsson, The University of Colorado at Boulder Acknowledgements: Some of the work presented is joint work with Vladimir Rokhlin and Mark Tygert at
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 informationIterative Methods for Sparse Linear Systems
Iterative Methods for Sparse Linear Systems Luca Bergamaschi e-mail: berga@dmsa.unipd.it - http://www.dmsa.unipd.it/ berga Department of Mathematical Methods and Models for Scientific Applications University
More informationAn Angular Multigrid Acceleration Method for S n Equations with Highly Forward-Peaked Scattering
An Angular Multigrid Acceleration Method for S n Equations with Highly Forward-Peaked Scattering Bruno Turcksin & Jean C. Ragusa & Jim E. Morel Texas A&M University, Dept. of Nuclear Engineering Bruno
More informationMotivation: Sparse matrices and numerical PDE's
Lecture 20: Numerical Linear Algebra #4 Iterative methods and Eigenproblems Outline 1) Motivation: beyond LU for Ax=b A little PDE's and sparse matrices A) Temperature Equation B) Poisson Equation 2) Splitting
More informationIterative Methods and Multigrid
Iterative Methods and Multigrid Part 3: Preconditioning 2 Eric de Sturler Preconditioning The general idea behind preconditioning is that convergence of some method for the linear system Ax = b can be
More information1. Thomas-Fermi method
1. Thomas-Fermi method We consider a system of N electrons in a stationary state, that would obey the stationary Schrödinger equation: h i m + 1 v(r i,r j ) Ψ(r 1,...,r N ) = E i Ψ(r 1,...,r N ). (1.1)
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 informationOpen-source finite element solver for domain decomposition problems
1/29 Open-source finite element solver for domain decomposition problems C. Geuzaine 1, X. Antoine 2,3, D. Colignon 1, M. El Bouajaji 3,2 and B. Thierry 4 1 - University of Liège, Belgium 2 - University
More informationMultigrid Method for 2D Helmholtz Equation using Higher Order Finite Difference Scheme Accelerated by Krylov Subspace
201, TextRoad Publication ISSN: 2090-27 Journal of Applied Environmental and Biological Sciences www.textroad.com Multigrid Method for 2D Helmholtz Equation using Higher Order Finite Difference Scheme
More informationAlgebraic Multigrid Preconditioners for Computing Stationary Distributions of Markov Processes
Algebraic Multigrid Preconditioners for Computing Stationary Distributions of Markov Processes Elena Virnik, TU BERLIN Algebraic Multigrid Preconditioners for Computing Stationary Distributions of Markov
More informationA Central Compact-Reconstruction WENO Method for Hyperbolic Conservation Laws
A Central Compact-Reconstruction WENO Method for Hyperbolic Conservation Laws Kilian Cooley 1 Prof. James Baeder 2 1 Department of Mathematics, University of Maryland - College Park 2 Department of Aerospace
More informationIncomplete Cholesky preconditioners that exploit the low-rank property
anapov@ulb.ac.be ; http://homepages.ulb.ac.be/ anapov/ 1 / 35 Incomplete Cholesky preconditioners that exploit the low-rank property (theory and practice) Artem Napov Service de Métrologie Nucléaire, Université
More informationAsymptotics of generalized eigenfunctions on manifold with Euclidean and/or hyperbolic ends
Asymptotics of generalized eigenfunctions on manifold with Euclidean and/or hyperbolic ends Kenichi ITO (University of Tokyo) joint work with Erik SKIBSTED (Aarhus University) 3 July 2018 Example: Free
More informationIterative methods for positive definite linear systems with a complex shift
Iterative methods for positive definite linear systems with a complex shift William McLean, University of New South Wales Vidar Thomée, Chalmers University November 4, 2011 Outline 1. Numerical solution
More 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 informationMATH325 - QUANTUM MECHANICS - SOLUTION SHEET 11
MATH35 - QUANTUM MECHANICS - SOLUTION SHEET. The Hamiltonian for a particle of mass m moving in three dimensions under the influence of a three-dimensional harmonic oscillator potential is Ĥ = h m + mω
More informationEfficient Solvers for Stochastic Finite Element Saddle Point Problems
Efficient Solvers for Stochastic Finite Element Saddle Point Problems Catherine E. Powell c.powell@manchester.ac.uk School of Mathematics University of Manchester, UK Efficient Solvers for Stochastic Finite
More informationUncertainty in Molecular Photoionization!
Uncertainty in Molecular Photoionization! Robert R. Lucchese! Department of Chemistry! Texas A&M University Collaborators:! At Texas A&M: R. Carey, J. Lopez, J. Jose! At ISMO, Orsay, France: D. Dowek and
More informationPoisson Solver, Pseudopotentials, Atomic Forces in the BigDFT code
CECAM Tutorial on Wavelets in DFT, CECAM - LYON,, in the BigDFT code Kernel Luigi Genovese L_Sim - CEA Grenoble 28 November 2007 Outline, Kernel 1 The with Interpolating Scaling Functions in DFT for Interpolating
More informationHigh Performance Nonlinear Solvers
What is a nonlinear system? High Performance Nonlinear Solvers Michael McCourt Division Argonne National Laboratory IIT Meshfree Seminar September 19, 2011 Every nonlinear system of equations can be described
More informationQuantum wires, orthogonal polynomials and Diophantine approximation
Quantum wires, orthogonal polynomials and Diophantine approximation Introduction Quantum Mechanics (QM) is a linear theory Main idea behind Quantum Information (QI): use the superposition principle of
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 informationBETI for acoustic and electromagnetic scattering
BETI for acoustic and electromagnetic scattering O. Steinbach, M. Windisch Institut für Numerische Mathematik Technische Universität Graz Oberwolfach 18. Februar 2010 FWF-Project: Data-sparse Boundary
More 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 informationLECTURE 5 APPLICATIONS OF BDIE METHOD: ACOUSTIC SCATTERING BY INHOMOGENEOUS ANISOTROPIC OBSTACLES DAVID NATROSHVILI
LECTURE 5 APPLICATIONS OF BDIE METHOD: ACOUSTIC SCATTERING BY INHOMOGENEOUS ANISOTROPIC OBSTACLES DAVID NATROSHVILI Georgian Technical University Tbilisi, GEORGIA 0-0 1. Formulation of the corresponding
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 informationChapter 7: Bounded Operators in Hilbert Spaces
Chapter 7: Bounded Operators in Hilbert Spaces I-Liang Chern Department of Applied Mathematics National Chiao Tung University and Department of Mathematics National Taiwan University Fall, 2013 1 / 84
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 informationChapter. 5 Bound States: Simple Case
Announcement Course webpage http://highenergy.phys.ttu.edu/~slee/2402/ Textbook PHYS-2402 Lecture 12 HW3 (due 3/2) 13, 15, 20, 31, 36, 41, 48, 53, 63, 66 ***** Exam: 3/12 Ch.2, 3, 4, 5 Feb. 26, 2015 Physics
More informationMultilevel preconditioning techniques with applications Yousef Saad Department of Computer Science and Engineering University of Minnesota
Multilevel preconditioning techniques with applications Yousef Saad Department of Computer Science and Engineering University of Minnesota Syracuse University, Sep. 20, 2013 First: Joint work with Ruipeng
More informationScattering theory of current-induced forces. Reinhold Egger Institut für Theoretische Physik, Univ. Düsseldorf
Scattering theory of current-induced forces Reinhold Egger Institut für Theoretische Physik, Univ. Düsseldorf Overview Current-induced forces in mesoscopic systems: In molecule/dot with slow mechanical
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 informationQuantum Physics III (8.06) Spring 2007 FINAL EXAMINATION Monday May 21, 9:00 am You have 3 hours.
Quantum Physics III (8.06) Spring 2007 FINAL EXAMINATION Monday May 21, 9:00 am You have 3 hours. There are 10 problems, totalling 180 points. Do all problems. Answer all problems in the white books provided.
More informationSplitting Iteration Methods for Positive Definite Linear Systems
Splitting Iteration Methods for Positive Definite Linear Systems Zhong-Zhi Bai a State Key Lab. of Sci./Engrg. Computing Inst. of Comput. Math. & Sci./Engrg. Computing Academy of Mathematics and System
More informationAsymptotic behaviour of the heat equation in twisted waveguides
Asymptotic behaviour of the heat equation in twisted waveguides Gabriela Malenová Faculty of Nuclear Sciences and Physical Engineering, CTU, Prague Nuclear Physics Institute, AS ČR, Řež Graphs and Spectra,
More informationCourant Institute of Mathematical Sciences, New York University, New York, NY 10012,
EFFICIENT VARIABLE-COEFFICIENT FINITE-VOLUME STOKES SOLVERS MINGCHAO CAI, ANDY NONAKA, JOHN B. BELL, BOYCE E. GRIFFITH, AND ALEKSANDAR DONEV Abstract. We investigate several robust preconditioners for
More informationMath 115 ( ) Yum-Tong Siu 1. Derivation of the Poisson Kernel by Fourier Series and Convolution
Math 5 (006-007 Yum-Tong Siu. Derivation of the Poisson Kernel by Fourier Series and Convolution We are going to give a second derivation of the Poisson kernel by using Fourier series and convolution.
More informationExtended Foldy Lax Approximation on Multiple Scattering
Mathematical Modelling and Analysis Publisher: Taylor&Francis and VGTU Volume 9 Number, February 04, 85 98 http://www.tandfonline.com/tmma http://dx.doi.org/0.3846/3969.04.893454 Print ISSN: 39-69 c Vilnius
More informationLecture 5 Scattering theory, Born Approximation. SS2011: Introduction to Nuclear and Particle Physics, Part 2 2
Lecture 5 Scattering theory, Born Approximation SS2011: Introduction to Nuclear and Particle Physics, Part 2 2 1 Scattering amplitude We are going to show here that we can obtain the differential cross
More informationMASSACHUSETTS INSTITUTE OF TECHNOLOGY 5.76 Modern Topics in Physical Chemistry Spring, Problem Set #2
Reading Assignment: Bernath Chapter 5 MASSACHUSETTS INSTITUTE O TECHNOLOGY 5.76 Modern Topics in Physical Chemistry Spring 994 Problem Set # The following handouts also contain useful information: C &
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 informationarxiv: v1 [quant-ph] 12 Mar 2015
Quasi Sturmian Basis in Two-Electron Continuum Problems A. S. Zaytsev 1, L. U. Ancarani, and S. A. Zaytsev 1 1 Pacific National University, Khabarovsk, 680035, Russia and Equipe TMS, SRSMC, UMR CNRS 7565,
More informationOptimal Control of the Schrödinger Equation on Many-Core Architectures
Optimal Control of the Schrödinger Equation on Many-Core Architectures Manfred Liebmann Institute for Mathematics and Scientific Computing University of Graz manfred.liebmann@uni-graz.at July 15, 2014
More informationMultigrid and Domain Decomposition Methods for Electrostatics Problems
Multigrid and Domain Decomposition Methods for Electrostatics Problems Michael Holst and Faisal Saied Abstract. We consider multigrid and domain decomposition methods for the numerical solution of electrostatics
More informationPreconditioning techniques for highly indefinite linear systems Yousef Saad Department of Computer Science and Engineering University of Minnesota
Preconditioning techniques for highly indefinite linear systems Yousef Saad Department of Computer Science and Engineering University of Minnesota Institut C. Jordan, Lyon, 25/03/08 Introduction: Linear
More informationElectron-helium scattering within the S-wave model
PHYSICAL REVIEW A, VOLUME 65, 032701 Electron-helium scattering within the S-wave model Chris Plottke, 1 Igor Bray, 1 Dmitry V. Fursa, 2 and Andris T. Stelbovics 1 1 Centre for Atomic, Molecular and Surface
More informationRadial Basis Functions generated Finite Differences (RBF-FD) for Solving High-Dimensional PDEs in Finance
1 / 27 Radial Basis Functions generated Finite Differences (RBF-FD) for Solving High-Dimensional PDEs in Finance S. Milovanović, L. von Sydow Uppsala University Department of Information Technology Division
More informationTheory of optically thin emission line spectroscopy
Theory of optically thin emission line spectroscopy 1 Important definitions In general the spectrum of a source consists of a continuum and several line components. Processes which give raise to the continuous
More informationApplied Nuclear Physics (Fall 2006) Lecture 2 (9/11/06) Schrödinger Wave Equation
22.101 Applied Nuclear Physics (Fall 2006) Lecture 2 (9/11/06) Schrödinger Wave Equation References -- R. M. Eisberg, Fundamentals of Modern Physics (Wiley & Sons, New York, 1961). R. L. Liboff, Introductory
More informationFEM and Sparse Linear System Solving
FEM & sparse system solving, Lecture 7, Nov 3, 2017 1/46 Lecture 7, Nov 3, 2015: Introduction to Iterative Solvers: Stationary Methods http://people.inf.ethz.ch/arbenz/fem16 Peter Arbenz Computer Science
More information