On the polarizability and capacitance of the cube
|
|
- Dominick Jared Stevenson
- 5 years ago
- Views:
Transcription
1 On the polarizability and capacitance of the cube Johan Helsing Lund University Double Seminar at CSC, May 10, 2012
2 Acknowledgement The work presented has in part been carried out in cooperation with or supported by: Alexandru Aleman Anders Karlsson Ross McPhedran Graeme Milton Rikard Ojala Karl-Mikael Perfekt Daniel Sjöberg the Swedish Research Council
3 A dielectric inclusion in a dielectric medium e ǫ 2 ǫ 1 V E S ν r
4 The electrostatic boundary value problem lim U(r) = e r U(r) continuous ǫ 2 ν r U int (r) = ǫ 1 ν r U ext (r) U(r) = 0 U(r) = 0
5 The polarizability α is given by α = (ǫ 2 ǫ 1 ) (e U(r))dr. V
6 The polarizability α is given by α = (ǫ 2 ǫ 1 ) When does α exist? V (e U(r))dr.
7 The polarizability α is given by α = (ǫ 2 ǫ 1 ) V (e U(r))dr. When does α exist? For what ǫ 1, ǫ 2 does the cube have an α?
8 The polarizability α is given by α = (ǫ 2 ǫ 1 ) V (e U(r))dr. When does α exist? For what ǫ 1, ǫ 2 does the cube have an α? Is it difficult to compute α?
9 The polarizability α is given by α = (ǫ 2 ǫ 1 ) V (e U(r))dr. When does α exist? For what ǫ 1, ǫ 2 does the cube have an α? Is it difficult to compute α? Does a small rounding of a corner correspond to a small relative perturbation of input data?
10 The polarizability α is given by α = (ǫ 2 ǫ 1 ) V (e U(r))dr. When does α exist? For what ǫ 1, ǫ 2 does the cube have an α? Is it difficult to compute α? Does a small rounding of a corner correspond to a small relative perturbation of input data? Is it important to compute α?
11 The polarizability α is given by α = (ǫ 2 ǫ 1 ) V (e U(r))dr. When does α exist? For what ǫ 1, ǫ 2 does the cube have an α? Is it difficult to compute α? Does a small rounding of a corner correspond to a small relative perturbation of input data? Is it important to compute α? How high accuracy should one demand of a fast α-solver?
12 The polarizability α is given by α = (ǫ 2 ǫ 1 ) V (e U(r))dr. When does α exist? For what ǫ 1, ǫ 2 does the cube have an α? Is it difficult to compute α? Does a small rounding of a corner correspond to a small relative perturbation of input data? Is it important to compute α? How high accuracy should one demand of a fast α-solver? Are there many failed attempts to compute α for non-smooth S in the literature?
13 Equivalent integral equation formulation Ansatz which solves the PDE and takes care of r : U(r) = e r + G(r,r )ρ(r )dσ r, where ρ(r) is an unknown layer density. S
14 Equivalent integral equation formulation Ansatz which solves the PDE and takes care of r : U(r) = e r + G(r,r )ρ(r )dσ r, S where ρ(r) is an unknown layer density. Insertion of U(r) in the boundary condition at S: (ǫ 2 + ǫ 1 ) (ǫ 2 ǫ 1 ) ρ(r) + 2 G(r,r )ρ(r )dσ r = 2(e ν r ), r S. ν r S
15 Equivalent integral equation formulation Ansatz which solves the PDE and takes care of r : U(r) = e r + G(r,r )ρ(r )dσ r, S where ρ(r) is an unknown layer density. Insertion of U(r) in the boundary condition at S: (ǫ 2 + ǫ 1 ) (ǫ 2 ǫ 1 ) ρ(r) + 2 G(r,r )ρ(r )dσ r = 2(e ν r ), r S. ν r Abbreviated notation: S ( z + K) ρ(r) = g(r). (1)
16 Equivalent integral equation formulation Ansatz which solves the PDE and takes care of r : U(r) = e r + G(r,r )ρ(r )dσ r, S where ρ(r) is an unknown layer density. Insertion of U(r) in the boundary condition at S: (ǫ 2 + ǫ 1 ) (ǫ 2 ǫ 1 ) ρ(r) + 2 G(r,r )ρ(r )dσ r = 2(e ν r ), r S. ν r Abbreviated notation: The polarizability from ρ(r): α = ǫ 1 S ( z + K) ρ(r) = g(r). (1) S ρ(r)(e r) dσ r. (2)
17 When does have a (unique) solution? ( z + K) ρ(r) = g(r)
18 When does ( z + K) ρ(r) = g(r) have a (unique) solution? Well, when z is not an eigenvalue of K.
19 When does ( z + K) ρ(r) = g(r) have a (unique) solution? Well, when z is not an eigenvalue of K. What does the spectrum of K look like?
20 When does ( z + K) ρ(r) = g(r) have a (unique) solution? Well, when z is not an eigenvalue of K. What does the spectrum of K look like? If S is smooth, then K is compact and has a discrete spectrum accumulating at zero. All eigenvalues are real and less than one in modulus.
21 When does ( z + K) ρ(r) = g(r) have a (unique) solution? Well, when z is not an eigenvalue of K. What does the spectrum of K look like? If S is smooth, then K is compact and has a discrete spectrum accumulating at zero. All eigenvalues are real and less than one in modulus. The electrostatic equation is uniquely solvable for z / [ 1,1], and often otherwise as well.
22 When does ( z + K) ρ(r) = g(r) have a (unique) solution? Well, when z is not an eigenvalue of K. What does the spectrum of K look like? If S is smooth, then K is compact and has a discrete spectrum accumulating at zero. All eigenvalues are real and less than one in modulus. The electrostatic equation is uniquely solvable for z / [ 1,1], and often otherwise as well. Eigenvalues z correspond to PLASMONS in α(z).
23 The Superellipse The superellipse r 1 k + r 2 k = 1 is a circle for k = 2 and numerically a square for k = Superellipse at k=2, 8, 32,
24 What does the spectrum of K look like for a superellipse?
25 What does the spectrum of K look like for a superellipse? superellipse r 1 k + r 2 k = eigenvalues z i dark plasmons bright plasmons k
26 A zoom at low k reveals a fascinating pattern superellipse r k + r k = dark plasmons bright plasmons positive eigenvalues z i k Are there any phase transitions or critical points here?
27 What does the polarizability α(z), z [ 1,1], look like for a superellipse (almost square)?
28 What does the polarizability α(z), z [ 1,1], look like for a superellipse (almost square)? superellipse r 1 k + r 2 k =1 with k= α(x) x A myriad of PLASMONS.
29 Can we guess the spectrum of K for a square by studying the superellipse in the limit k? superellipse r 1 k + r 2 k = eigenvalues z i dark plasmons bright plasmons k
30 Can we guess the spectrum of K for a square by studying the superellipse in the limit k? superellipse r 1 k + r 2 k = eigenvalues z i dark plasmons bright plasmons k The answer depends on what function space we consider.
31 What does the limit polarizability α + (x) = lim y 0 + α(x + iy) look like for a square? Does it exist at all? and if so, in what function spaces lie the corresponding limit solutions ρ + (r) and U + (r)?
32 What does the limit polarizability α + (x) = lim y 0 + α(x + iy) look like for a square? Does it exist at all? and if so, in what function spaces lie the corresponding limit solutions ρ + (r) and U + (r)? The square: polarizability α + (x) real{α + (x)} 8 imag{α + (x)} x Yes, α + (x) exists, but do we see any PLASMONS?
33 Comparison: superellipse r 1 k + r 2 k =1 with k=10 12 The square: polarizability α(x) x α + (x) real{α + (x)} 8 imag{α + (x)} x
34 Theory for Lipschitz sets The polarizability α + (x) exists. The polarizability α(z) admits an integral representation dµ(x) α(z) = x z = µ (x)dx, σ µ x z R where σ µ = {x : µ (x) > 0} and the measure dµ(x) relate to α + (x) as µ (x) = I{α + (x)}/π, x R. The spectral measure is of great interest in theoretical material science. One can show that for z L, L : σ µ L [ 1,1], the electrostatic equation has no finite-energy solution U(r). For all other z there are solutions U(r). This has consequences.
35 What about the cube?
36 What about the cube? The cube: polarizability real{α + (x)} imag{α + (x)} 10 α + (x) x No one has succeeded in computing this (in its entirety) before. The reason being, perhaps, the behavior of U + (r) and ρ + (r).
37 Some cube records Table: Reference values, estimated relative errors, and best previous estimates of the polarizability α + (x) and the capacitance C of the cube. present reference values relerr previous results relerr α + ( 1) α + ( 0.6) i 10 8 α + (0.25) i 10 7 α + (1) C
38 10 0 The cube: convergence Estimated relative error in α + (x) and C α + ( 1) α + ( 0.6) α + (0.25) α + (1) C Number of discretization points n Figure: Convergence of α + (x) and the capacitance C of a unit cube. The values of x correspond to: cube with infinite permittivity (x = 1), resonance in corners (x = 0.6), resonance along edges (x = 0.25), and cube with zero permittivity (x = 1).
39 Table: Numerical progress for the capacitance of the unit cube. value year Author Reitan Higgins Greenspan Silverman Cochran Brown Goto Shi Yoshida Douglas Zho Hubbard Given Hubbard Douglas Read Mansfield Douglas Garboczi Bai Lonngren Hwang Mascagni Hwang Mascagni Wintle Mascagni Simonov Read Ong Lim Lazić Štefančić Abraham Mukhopadhyay Majumdar Hwang Mascagni Won Bontzios Dimopoulos Hatzopoulos Helsing Perfekt
40 Numerical method The computation of α + (x) relies on the solution ρ + (r) to (I + λk) ρ + (r) = λg(r), but how is it possible to obtain accurate numerical approximations to ρ + (r), which barely lies in W 1,?
41 Numerical method The computation of α + (x) relies on the solution ρ + (r) to (I + λk) ρ + (r) = λg(r), but how is it possible to obtain accurate numerical approximations to ρ + (r), which barely lies in W 1,? Well, thanks to K = K + K, ρ(r) = (I + λk )ρ(r), ( I + λk (I + λk ) 1) ρ(r) = λg(r), which gives a solution ρ + (r) in L.
42 Discretization Discretization of (I + λk) ρ(r) = λg(r) gives (I coa + λk coa R) ρ coa = λg coa, (3) where the small block diagonal matrix R = P T W (I fin + λk fin ) 1 P is obtained via a fast recursion.
43 Polarizability Superellipse and Square The cube Numerical method Details To the movies γ 1 Figure: A coarse mesh with eight quadrature panels on a closed surface. A fine mesh with 14 panels is created by refinement in a direction toward the corner γ nz = nz = Figure: Kfin = K fin + K fin nz =
44 Details (I + λk)ρ(r) = λg(r). Assumption: K is compact some distance away from the corners γ j and g(r) is piecewise smooth. Let K(r,r ) be the kernel of K. Split K(r,r ) into two functions K(r,r ) = K (r,r ) + K (r,r ), where K (r,r ) is zero except for when r and r simultaneously lie close to the same γ j. In this latter case K (r,r ) is zero. The kernel split corresponds to an operator split K = K + K, where K is compact. Discretization on a fine mesh gives (I fin + λk fin + λk fin ) ρ fin = λg fin.
45 Details The transform ρ fin = (I fin + λk fin ) 1 ρ fin leads to ( I fin + λk fin (I fin + λk fin ) 1) ρ fin = λg fin. Note: K fin (I fin + λk fin ) 1 is, in practice, a compact operator ρ fin is piecewise smooth neither g fin, ρ fin, nor K fin needs a fine mesh for resolution.
46 Details Let P be a prolongation operator from the coarse mesh to the fine mesh and let Q be a restriction operator in the opposite direction Then QP = I. g fin = Pg coa ρ fin = P ρ coa K fin = PK coa PT W. Compare reduced SVD. Here P W = W fin PWcoa 1 and W has quadrature weights on the diagonal. Taken together ( ) I coa + λk coa PT W (I fin + λk fin ) 1 P ρ coa = λg coa. (4)
47 Details With (4) becomes R = P T W (I fin + λk fin ) 1 P (I coa + λk coa R) ρ coa = λg coa. We only need the fine mesh for the construction of the small block diagonal matrix R. One block per corner. The blocks of R can be obtained via a fast recursion, where step i inverts and compresses contributions to R involving the outermost panels on level i of a locally n-ply refined mesh. Newton s method can accelerate the recursion in situations that demand extreme refinement (very many recursion steps). This corresponds to inverting a giant full matrix in sublinear time. Boundary singularities do no longer pose any practical problems.
48 The fast recursion R i = P T Wbc ( F{R 1 (i 1) } + I b + λk b) 1 Pbc, i = 1,...,n. part of coarse mesh i = 1 i = 2 refined n = 4 i = 3 Figure: Recursion on a refined mesh surrounding a corner.
49 Matlab Kmat=Kinit(zloc,wloc,nzloc,96); MAT=eye(96)+lambda*Kmat; starl=[17:80]; R=inv(MAT(starL,starL)); myerr=1; while myerr>eps Rold=R; MAT(starL,starL)=inv(R); R=Pwbc *inv(mat)*pbc; myerr=norm(r-rold, fro )/norm(r, fro ); end
50 An extreme example 20 years too late? Figure: A unit cell of a random checkerboard with a million squares. The permittivity ratio is ǫ 2 /ǫ 1 = 10 6 and the effective permittivity is obtained with a relative error of 10 9.
51 References J. Helsing and K.-M. Perfekt (2012) On the polarizability and capacitance of the cube, arxiv: v2. J. Helsing, R.C. McPhedran, and G.W. Milton (2011) Spectral super-resolution in metamaterial composites, New J. Phys., 13(11), J. Helsing (2011) The effective conductivity of arrays of squares: large random unit cells and extreme contrast ratios, J. Comput. Phys., 230(20), J. Helsing (2011) A fast and stable solver for singular integral equations on piecewise smooth curves, SIAM J. Sci. Comput., 33(1), J. Helsing (2009) Integral equation methods for elliptic problems with boundary conditions of mixed type, J. Comput. Phys., 228(23), pp J. Helsing and R. Ojala (2008) Corner singularities for elliptic problems: integral equations, graded meshes, quadrature, and compressed inverse preconditioning, J. Comput. Phys., 227(20),
52 Staggered array of square cylinders at p = Staggered array of square cylinders at p = effective permittivity ε eff (σ) real{εeff} imag{εeff} permittivity ratio σ ε eff (σ)ε eff (1/σ ) permittivity ratio σ Figure: Effective permittivity (polarizability) and error estimates.
Solving integral equations on non-smooth domains
Solving integral equations on non-smooth domains Johan Helsing Lund University Talk at the IMA Hot Topics workshop, August 2-5, 2010 σ yy A B C σ yy Acknowledgement The work presented in this talk has
More informationExplicit kernel-split panel-based Nyström schemes for planar or axisymmetric Helmholtz problems
z Explicit kernel-split panel-based Nyström schemes for planar or axisymmetric Helmholtz problems Johan Helsing Lund University Talk at Integral equation methods: fast algorithms and applications, Banff,
More informationThe effective conductivity of arrays of squares: large random unit cells and extreme contrast ratios
The effective conductivity of arrays of squares: large random unit cells and extreme contrast ratios Johan Helsing 1 Numerical Analysis, Centre for Mathematical Sciences, Lund University, Box 118, SE-221
More informationSolving integral equations on piecewise smooth boundaries using the RCIP method: a tutorial
Solving integral equations on piecewise smooth boundaries using the RCIP method: a tutorial Helsing, Johan Published in: Abstract and Applied Analysis DOI: 10.1155/2013/938167 Published: 2013-01-01 Link
More informationOn the polarizability and capacitance of the cube
On the polarizability and capacitance of the cube Johan Helsing, Karl-Mikael Perfekt Centre for Mathematical ciences Lund University, Box 118, E-221 00 Lund, weden Abstract An efficient integral equation
More informationSpectral super-resolution in metamaterial composites
Spectral super-resolution in metamaterial composites J. Helsing 1, R. C. McPhedran 2 and G. W. Milton 3 1 Numerical Analysis, Centre for Mathematical Sciences, Lund University, Box 118, SE-221 00 Lund,
More informationThe effective conductivity of random checkerboards
The effective conductivity of random checkerboards Johan Helsing 1 Numerical Analysis, Centre for Mathematical Sciences, Lund University, Box 118, SE-221 00 LUND, Sweden Abstract An algorithm is presented
More informationA robust and accurate solver of Laplace s equation with general boundary conditions on general domains in the plane. 1.
A robust and accurate solver of Laplace s equation with general boundary conditions on general domains in the plane Rikard Ojala Numerical Analysis, Centre for Mathematical Sciences, Lund University, Box
More informationAnalysis of the anomalous localized resonance
Hyundae Lee(Inha University, Korea) Joint work with Habib Ammari, Giulio Ciraolo, Hyeonbae Kang, Graeme Milton. UCI, June 22, 2012 A conference on inverse problems in honor of Gunther Uhlmann Outline Introduction
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 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 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 informationWeighted Regularization of Maxwell Equations Computations in Curvilinear Polygons
Weighted Regularization of Maxwell Equations Computations in Curvilinear Polygons Martin Costabel, Monique Dauge, Daniel Martin and Gregory Vial IRMAR, Université de Rennes, Campus de Beaulieu, Rennes,
More informationAn explicit kernel-split panel-based Nyström scheme for integral equations on axially symmetric surfaces
An explicit kernel-split panel-based Nyström scheme for integral equations on axially symmetric surfaces Helsing, Johan; Karlsson, Anders Published in: Journal of Computational Physics DOI: 10.1016/j.jcp.014.04.053
More informationConvergence Order Studies for Elliptic Test Problems with COMSOL Multiphysics
Convergence Order Studies for Elliptic Test Problems with COMSOL Multiphysics Shiming Yang and Matthias K. Gobbert Abstract. The convergence order of finite elements is related to the polynomial order
More informationarxiv: v2 [math.ap] 29 Sep 2017
The spectra of harmonic layer potential operators on domains with rotationally symmetric conical points Johan Helsing Lund University, Box 8, SE-22 Lund, Sweden Karl-Mikael Perfekt Department of Mathematics
More informationDetermination of normalized electric eigenfields in microwave cavities with sharp edges
Determination of normalied electric eigenfields in microwave cavities with sharp edges Johan Helsing and Anders Karlsson October 6, 5 Abstract The magnetic field integral equation for aially symmetric
More informationAn Introduction to Algebraic Multigrid (AMG) Algorithms Derrick Cerwinsky and Craig C. Douglas 1/84
An Introduction to Algebraic Multigrid (AMG) Algorithms Derrick Cerwinsky and Craig C. Douglas 1/84 Introduction Almost all numerical methods for solving PDEs will at some point be reduced to solving A
More informationThe Distribution Function
The Distribution Function As we have seen before the distribution function (or phase-space density) f( x, v, t) d 3 x d 3 v gives a full description of the state of any collisionless system. Here f( x,
More informationThe Helmholtz Equation
The Helmholtz Equation Seminar BEM on Wave Scattering Rene Rühr ETH Zürich October 28, 2010 Outline Steklov-Poincare Operator Helmholtz Equation: From the Wave equation to Radiation condition Uniqueness
More informationA Recursive Trust-Region Method for Non-Convex Constrained Minimization
A Recursive Trust-Region Method for Non-Convex Constrained Minimization Christian Groß 1 and Rolf Krause 1 Institute for Numerical Simulation, University of Bonn. {gross,krause}@ins.uni-bonn.de 1 Introduction
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 informationA SHORT NOTE COMPARING MULTIGRID AND DOMAIN DECOMPOSITION FOR PROTEIN MODELING EQUATIONS
A SHORT NOTE COMPARING MULTIGRID AND DOMAIN DECOMPOSITION FOR PROTEIN MODELING EQUATIONS MICHAEL HOLST AND FAISAL SAIED Abstract. We consider multigrid and domain decomposition methods for the numerical
More informationMathematics Qualifying Exam Study Material
Mathematics Qualifying Exam Study Material The candidate is expected to have a thorough understanding of engineering mathematics topics. These topics are listed below for clarification. Not all instructors
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 informationSummer course: Plasmonic resonances
Summer course: Plasmonic resonances Faouzi TRIKI Grenoble-Alpes University, France Journées d été des mathématiciens tunisiens à l étranger Mediterranean Institute For The Mathematical Sciences July 20,
More informationLinear Algebra Massoud Malek
CSUEB Linear Algebra Massoud Malek Inner Product and Normed Space In all that follows, the n n identity matrix is denoted by I n, the n n zero matrix by Z n, and the zero vector by θ n An inner product
More informationNonlinear stability of time-periodic viscous shocks. Margaret Beck Brown University
Nonlinear stability of time-periodic viscous shocks Margaret Beck Brown University Motivation Time-periodic patterns in reaction-diffusion systems: t x Experiment: chemical reaction chlorite-iodite-malonic-acid
More informationarxiv: v2 [math.ap] 9 Sep 2014
Localized and complete resonance in plasmonic structures Hoai-Minh Nguyen and Loc Hoang Nguyen April 23, 2018 arxiv:1310.3633v2 [math.ap] 9 Sep 2014 Abstract This paper studies a possible connection between
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 informationPresentation of Dark Matter as violation of superposition principle: is it quantum non-locality of energy? Abstract
Presentation of Dark Matter as violation of superposition principle: is it quantum non-locality of energy? Dmitri Martila (eestidima@gmail.com) Former Tartu University Lääne 9-51, Tartu 50605, Estonia
More informationComplex Analysis Math 185A, Winter 2010 Final: Solutions
Complex Analysis Math 85A, Winter 200 Final: Solutions. [25 pts] The Jacobian of two real-valued functions u(x, y), v(x, y) of (x, y) is defined by the determinant (u, v) J = (x, y) = u x u y v x v y.
More informationNumerical Programming I (for CSE)
Technische Universität München WT 1/13 Fakultät für Mathematik Prof. Dr. M. Mehl B. Gatzhammer January 1, 13 Numerical Programming I (for CSE) Tutorial 1: Iterative Methods 1) Relaxation Methods a) Let
More informationNumerical Methods I Non-Square and Sparse Linear Systems
Numerical Methods I Non-Square and Sparse Linear Systems Aleksandar Donev Courant Institute, NYU 1 donev@courant.nyu.edu 1 MATH-GA 2011.003 / CSCI-GA 2945.003, Fall 2014 September 25th, 2014 A. Donev (Courant
More informationAN ITERATION. In part as motivation, we consider an iteration method for solving a system of linear equations which has the form x Ax = b
AN ITERATION In part as motivation, we consider an iteration method for solving a system of linear equations which has the form x Ax = b In this, A is an n n matrix and b R n.systemsof this form arise
More informationIII. Iterative Monte Carlo Methods for Linear Equations
III. Iterative Monte Carlo Methods for Linear Equations In general, Monte Carlo numerical algorithms may be divided into two classes direct algorithms and iterative algorithms. The direct algorithms provide
More informationA Balancing Algorithm for Mortar Methods
A Balancing Algorithm for Mortar Methods Dan Stefanica Baruch College, City University of New York, NY 11, USA Dan Stefanica@baruch.cuny.edu Summary. The balancing methods are hybrid nonoverlapping Schwarz
More informationOverview. 1. What range of ε eff, µ eff parameter space is accessible to simple metamaterial geometries? ``
MURI-Transformational Electromagnetics Innovative use of Metamaterials in Confining, Controlling, and Radiating Intense Microwave Pulses University of New Mexico August 21, 2012 Engineering Dispersive
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 informationA COLLOCATION METHOD FOR SOLVING THE EXTERIOR NEUMANN PROBLEM
TUDIA UNIV. BABEŞ BOLYAI, MATHEMATICA, Volume XLVIII, Number 3, eptember 2003 A COLLOCATION METHOD FOR OLVING THE EXTERIOR NEUMANN PROBLEM ANDA MICULA Dedicated to Professor Gheorghe Micula at his 60 th
More informationLecture 4. Diffusing photons and superradiance in cold gases
Lecture 4 Diffusing photons and superradiance in cold gases Model of disorder-elastic mean free path and group velocity. Dicke states- Super- and sub-radiance. Scattering properties of Dicke states. Multiple
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 informationA Broadband Flexible Metamaterial Absorber Based on Double Resonance
Progress In Electromagnetics Research Letters, Vol. 46, 73 78, 2014 A Broadband Flexible Metamaterial Absorber Based on Double Resonance ong-min Lee* Abstract We present a broadband microwave metamaterial
More informationIntegration, differentiation, and root finding. Phys 420/580 Lecture 7
Integration, differentiation, and root finding Phys 420/580 Lecture 7 Numerical integration Compute an approximation to the definite integral I = b Find area under the curve in the interval Trapezoid Rule:
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 informationProblems of Corner Singularities
Stuttgart 98 Problems of Corner Singularities Monique DAUGE Institut de Recherche MAthématique de Rennes Problems of Corner Singularities 1 Vertex and edge singularities For a polyhedral domain Ω and a
More informationIntroduction to Multigrid Methods Sampling Theory and Elements of Multigrid
Introduction to Multigrid Methods Sampling Theory and Elements of Multigrid Gustaf Söderlind Numerical Analysis, Lund University Contents V3.15 What is a multigrid method? Sampling theory Digital filters
More informationLine Search Methods for Unconstrained Optimisation
Line Search Methods for Unconstrained Optimisation Lecture 8, Numerical Linear Algebra and Optimisation Oxford University Computing Laboratory, MT 2007 Dr Raphael Hauser (hauser@comlab.ox.ac.uk) The Generic
More informationTowards Classification of Separable Pauli Equations
Proceedings of Institute of Mathematics of NAS of Ukraine 2002, Vol. 43, Part 2, 768 773 Towards Classification of Separable Pauli Equations Alexander ZHALIJ Institute of Mathematics of NAS of Ukraine,
More informationRegularizations of Singular Integral Operators (joint work with C. Liaw)
1 Outline Regularizations of Singular Integral Operators (joint work with C. Liaw) Sergei Treil Department of Mathematics Brown University April 4, 2014 2 Outline 1 Examples of Calderón Zygmund operators
More informationPreliminary Examination in Numerical Analysis
Department of Applied Mathematics Preliminary Examination in Numerical Analysis August 7, 06, 0 am pm. Submit solutions to four (and no more) of the following six problems. Show all your work, and justify
More informationChapter 5 HIGH ACCURACY CUBIC SPLINE APPROXIMATION FOR TWO DIMENSIONAL QUASI-LINEAR ELLIPTIC BOUNDARY VALUE PROBLEMS
Chapter 5 HIGH ACCURACY CUBIC SPLINE APPROXIMATION FOR TWO DIMENSIONAL QUASI-LINEAR ELLIPTIC BOUNDARY VALUE PROBLEMS 5.1 Introduction When a physical system depends on more than one variable a general
More informationPart III: Impurities in Luttinger liquids
Functional RG for interacting fermions... Part III: Impurities in Luttinger liquids 1. Luttinger liquids 2. Impurity effects 3. Microscopic model 4. Flow equations 5. Results S. Andergassen, T. Enss (Stuttgart)
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 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 informationANALOGY BETWEEN ELECTROSTATIC FIELD AND HEAT TRANSFER SIMPLE MODELS
ANALOGY BETWEEN ELECTROSTATIC FIELD AND HEAT TRANSFER SIMPLE MODELS PROF. RNDR ING. MILOSLAV KOŠEK, CSC. ING. JAN VODOLAN. Abstract: The perfect analogy allows solving of heat transfer problems by a lot
More informationSOLVING ILL-POSED LINEAR SYSTEMS WITH GMRES AND A SINGULAR PRECONDITIONER
SOLVING ILL-POSED LINEAR SYSTEMS WITH GMRES AND A SINGULAR PRECONDITIONER LARS ELDÉN AND VALERIA SIMONCINI Abstract. Almost singular linear systems arise in discrete ill-posed problems. Either because
More informationFaIMS: A fast algorithm for the inverse medium problem with multiple frequencies and multiple sources for the scalar Helmholtz equation
FaIMS: A fast algorithm for the inverse medium problem with multiple frequencies and multiple sources for the scalar Helmholtz equation S. Chaillat G. Biros Abstract We consider the inverse medium problem,
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 informationarxiv: v1 [math.na] 11 Jul 2011
Multigrid Preconditioner for Nonconforming Discretization of Elliptic Problems with Jump Coefficients arxiv:07.260v [math.na] Jul 20 Blanca Ayuso De Dios, Michael Holst 2, Yunrong Zhu 2, and Ludmil Zikatanov
More informationNumerical Analysis Preliminary Exam 10 am to 1 pm, August 20, 2018
Numerical Analysis Preliminary Exam 1 am to 1 pm, August 2, 218 Instructions. You have three hours to complete this exam. Submit solutions to four (and no more) of the following six problems. Please start
More informationStress intensity factors for a crack in front of an inclusion
Stress intensity factors for a crack in front of an inclusion Johan Helsing Department of Solid Mechanics and NADA, Royal Institute of Technology, SE-100 44 Stockholm, Sweden Email: helsing@nada.kth.se,
More informationFEM Convergence for PDEs with Point Sources in 2-D and 3-D
FEM Convergence for PDEs with Point Sources in -D and 3-D Kourosh M. Kalayeh 1, Jonathan S. Graf, and Matthias K. Gobbert 1 Department of Mechanical Engineering, University of Maryland, Baltimore County
More informationarxiv: v1 [math.ap] 11 Jan 2014
THE UNIFIED TRANSFORM FOR THE MODIFIED HELMHOLTZ EQUATION IN THE EXTERIOR OF A SQUARE A. S. FOKAS AND J. LENELLS arxiv:4.252v [math.ap] Jan 24 Abstract. The Unified Transform provides a novel method for
More informationSolutions to Tutorial 7 (Week 8)
The University of Sydney School of Mathematics and Statistics Solutions to Tutorial 7 (Week 8) MATH2962: Real and Complex Analysis (Advanced) Semester 1, 2017 Web Page: http://www.maths.usyd.edu.au/u/ug/im/math2962/
More information7.3 The Jacobi and Gauss-Siedel Iterative Techniques. Problem: To solve Ax = b for A R n n. Methodology: Iteratively approximate solution x. No GEPP.
7.3 The Jacobi and Gauss-Siedel Iterative Techniques Problem: To solve Ax = b for A R n n. Methodology: Iteratively approximate solution x. No GEPP. 7.3 The Jacobi and Gauss-Siedel Iterative Techniques
More information1. How much charge is stored in a capacitor, whose capacitance C = 2µF, connected to a 12V battery?
IMP 113: 2 nd test (Union College: Spring 2010) Instructions: 1. Read all directions. 2. In keeping with the Union College policy on academic honesty, you should neither accept nor provide unauthorized
More informationStandard Finite Elements and Weighted Regularization
Standard Finite Elements and Weighted Regularization A Rehabilitation Martin COSTABEL & Monique DAUGE Institut de Recherche MAthématique de Rennes http://www.maths.univ-rennes1.fr/~dauge Slides of the
More informationNonlinear equations. Norms for R n. Convergence orders for iterative methods
Nonlinear equations Norms for R n Assume that X is a vector space. A norm is a mapping X R with x such that for all x, y X, α R x = = x = αx = α x x + y x + y We define the following norms on the vector
More informationTHE solution of the absolute value equation (AVE) of
The nonlinear HSS-like iterative method for absolute value equations Mu-Zheng Zhu Member, IAENG, and Ya-E Qi arxiv:1403.7013v4 [math.na] 2 Jan 2018 Abstract Salkuyeh proposed the Picard-HSS iteration method
More informationBindel, Fall 2011 Intro to Scientific Computing (CS 3220) Week 3: Wednesday, Jan 9
Problem du jour Week 3: Wednesday, Jan 9 1. As a function of matrix dimension, what is the asymptotic complexity of computing a determinant using the Laplace expansion (cofactor expansion) that you probably
More informationAn Iteratively Regularized Projection Method for Nonlinear Ill-posed Problems
Int. J. Contemp. Math. Sciences, Vol. 5, 2010, no. 52, 2547-2565 An Iteratively Regularized Projection Method for Nonlinear Ill-posed Problems Santhosh George Department of Mathematical and Computational
More informationFast and accurate numerical solution to an elastostatic problem involving ten thousand randomly oriented cracks
Fast and accurate numerical solution to an elastostatic problem involving ten thousand randomly oriented cracks Johan Helsing (helsing@nada.kth.se) Department of Solid Mechanics and NADA, Royal Institute
More informationLECTURE-13 : GENERALIZED CAUCHY S THEOREM
LECTURE-3 : GENERALIZED CAUCHY S THEOREM VED V. DATAR The aim of this lecture to prove a general form of Cauchy s theorem applicable to multiply connected domains. We end with computations of some real
More informationLECTURE NOTES ELEMENTARY NUMERICAL METHODS. Eusebius Doedel
LECTURE NOTES on ELEMENTARY NUMERICAL METHODS Eusebius Doedel TABLE OF CONTENTS Vector and Matrix Norms 1 Banach Lemma 20 The Numerical Solution of Linear Systems 25 Gauss Elimination 25 Operation Count
More informationScientific Computing WS 2018/2019. Lecture 9. Jürgen Fuhrmann Lecture 9 Slide 1
Scientific Computing WS 2018/2019 Lecture 9 Jürgen Fuhrmann juergen.fuhrmann@wias-berlin.de Lecture 9 Slide 1 Lecture 9 Slide 2 Simple iteration with preconditioning Idea: Aû = b iterative scheme û = û
More informationEfficiency-based h- and hp-refinement strategies for finite element methods
NUMERICAL LINEAR ALGEBRA WITH APPLICATIONS Numer. Linear Algebra Appl. 7; : 5 [Version: /9/8 v.] Efficiency-based h- and hp-refinement strategies for finite element methods H. De Sterck, Thomas A. Manteuffel,
More informationSpace-time sparse discretization of linear parabolic equations
Space-time sparse discretization of linear parabolic equations Roman Andreev August 2, 200 Seminar for Applied Mathematics, ETH Zürich, Switzerland Support by SNF Grant No. PDFMP2-27034/ Part of PhD thesis
More informationEigenvalues of Robin Laplacians on infinite sectors and application to polygons
Eigenvalues of Robin Laplacians on infinite sectors and application to polygons Magda Khalile (joint work with Konstantin Pankrashkin) Université Paris Sud /25 Robin Laplacians on infinite sectors 1 /
More informationAn efficient multigrid solver based on aggregation
An efficient multigrid solver based on aggregation Yvan Notay Université Libre de Bruxelles Service de Métrologie Nucléaire Graz, July 4, 2012 Co-worker: Artem Napov Supported by the Belgian FNRS http://homepages.ulb.ac.be/
More informationCS 542G: The Poisson Problem, Finite Differences
CS 542G: The Poisson Problem, Finite Differences Robert Bridson November 10, 2008 1 The Poisson Problem At the end last time, we noticed that the gravitational potential has a zero Laplacian except at
More informationFractional Spectral and Spectral Element Methods
Fractional Calculus, Probability and Non-local Operators: Applications and Recent Developments Nov. 6th - 8th 2013, BCAM, Bilbao, Spain Fractional Spectral and Spectral Element Methods (Based on PhD thesis
More informationCLASS XII WB SET A PHYSICS
PHYSICS 1. Two cylindrical straight and very long non magnetic conductors A and B, insulated from each other, carry a current I in the positive and the negative z-direction respectively. The direction
More informationElectrodynamics I Midterm - Part A - Closed Book KSU 2005/10/17 Electro Dynamic
Electrodynamics I Midterm - Part A - Closed Book KSU 5//7 Name Electro Dynamic. () Write Gauss Law in differential form. E( r) =ρ( r)/ɛ, or D = ρ, E= electricfield,ρ=volume charge density, ɛ =permittivity
More informationElliptic Problems / Multigrid. PHY 604: Computational Methods for Physics and Astrophysics II
Elliptic Problems / Multigrid Summary of Hyperbolic PDEs We looked at a simple linear and a nonlinear scalar hyperbolic PDE There is a speed associated with the change of the solution Explicit methods
More informationModeling and Simulation of Artificial Core-Shell Based Nanodielectrics for Electrostatic Capacitors Applications
Modeling and Simulation of Artificial Core-Shell Based Nanodielectrics for Electrostatic Capacitors Applications -D. Musuwathi Ekanath 1, N. Badi* 1,2, A. Bensaoula 2 1 Center for Advanced Materials, University
More informationConservation of mass. Continuum Mechanics. Conservation of Momentum. Cauchy s Fundamental Postulate. # f body
Continuum Mechanics We ll stick with the Lagrangian viewpoint for now Let s look at a deformable object World space: points x in the object as we see it Object space (or rest pose): points p in some reference
More informationNotation. 0,1,2,, 1 with addition and multiplication modulo
Notation Q,, The set of all natural numbers 1,2,3, The set of all integers The set of all rational numbers The set of all real numbers The group of permutations of distinct symbols 0,1,2,,1 with addition
More informationAdaptation of the Lanczos Algorithm for the Solution of Buckling Eigenvalue Problems
Original Article Adaptation of the Lanczos Algorithm for the Solution of Buckling Eigenvalue Problems Abstract An adaptation of the conventional Lanczos algorithm is proposed to solve the general symmetric
More informationPhysics 202 Midterm 1 Practice Exam
Physics 202 Midterm 1 Practice Exam Name: Student ID: Section: TA (please circle): Daniel Crow Scott Douglas Yutao Gong Taylor Klaus Aaron Levine Andrew Loveridge Jason Milhone Hojin Yoo Instructions:
More information3. Numerical integration
3. Numerical integration... 3. One-dimensional quadratures... 3. Two- and three-dimensional quadratures... 3.3 Exact Integrals for Straight Sided Triangles... 5 3.4 Reduced and Selected Integration...
More informationInverses of Square Matrices
Inverses of Square Matrices A. Havens Department of Mathematics University of Massachusetts, Amherst February 23-26, 2018 Outline 1 Defining Inverses Inverses for Products and Functions Defining Inverse
More informationOptimal multilevel preconditioning of strongly anisotropic problems.part II: non-conforming FEM. p. 1/36
Optimal multilevel preconditioning of strongly anisotropic problems. Part II: non-conforming FEM. Svetozar Margenov margenov@parallel.bas.bg Institute for Parallel Processing, Bulgarian Academy of Sciences,
More informationFDTD Analysis on Optical Confinement Structure with Electromagnetic Metamaterial
Memoirs of the Faculty of Engineering, Okayama University, Vol. 44, pp. 1-6, January 21 FDTD Analysis on Optical Confinement Structure with Electromagnetic Metamaterial Shinji NAGAI, Ryosuke UMEDA, Kenji
More informationThe purpose of this lecture is to present a few applications of conformal mappings in problems which arise in physics and engineering.
Lecture 16 Applications of Conformal Mapping MATH-GA 451.001 Complex Variables The purpose of this lecture is to present a few applications of conformal mappings in problems which arise in physics and
More informationShort note on compact operators - Monday 24 th March, Sylvester Eriksson-Bique
Short note on compact operators - Monday 24 th March, 2014 Sylvester Eriksson-Bique 1 Introduction In this note I will give a short outline about the structure theory of compact operators. I restrict attention
More informationE. YABLONOVITCH photonic crystals by using level set methods
Appl. Phys. B 81, 235 244 (2005) DOI: 10.1007/s00340-005-1877-3 Applied Physics B Lasers and Optics C.Y. KAO 1, Maximizing band gaps in two-dimensional S. OSHER 2 E. YABLONOVITCH photonic crystals by using
More informationLAB 5. INTRODUCTION Finite Difference Method
LAB 5 In previous two computer labs, you have seen how the analytical techniques for solving electrostatic problems can be approximately solved using numerical methods. For example, you have approximated
More informationOverlapping Schwarz preconditioners for Fekete spectral elements
Overlapping Schwarz preconditioners for Fekete spectral elements R. Pasquetti 1, L. F. Pavarino 2, F. Rapetti 1, and E. Zampieri 2 1 Laboratoire J.-A. Dieudonné, CNRS & Université de Nice et Sophia-Antipolis,
More information