DOE PAGES title logo U.S. Department of Energy
Office of Scientific and Technical Information

Title: A hierarchical butterfly LU preconditioner for two-dimensional electromagnetic scattering problems involving open surfaces

Abstract

This paper introduces a hierarchical interpolative decomposition butterfly-LU factorization (H-IDBF-LU) preconditioner for solving two-dimensional electric-field integral equations (EFIEs) in electromagnetic scattering problems of perfect electrically conducting objects with open surfaces. H-IDBF-LU leverages the interpolative decomposition butterfly factorization (IDBF) to compress dense blocks of the discretized EFIE operator to expedite its application; this compressed operator also serves as an approximate LU factorization of the EFIE operator leading to an efficient preconditioner in iterative solvers. Both the memory requirement and computational cost of the H-IDBF-LU solver scale as $$O$$($$N$$ log2 $$N$$) in one iteration; the total number of iterations required for a reasonably good accuracy scales as $$O$$(1) to O(log2$$N$$) in all of our numerical tests. The efficacy and accuracy of the proposed preconditioned iterative solver are demonstrated via its application to a broad range of scatterers involving up to 100 million unknowns.

Authors:
 [1]; ORCiD logo [2]
  1. Lawrence Berkeley National Lab. (LBNL), Berkeley, CA (United States)
  2. National Univ. of Singapore (Singapore). Dept. of Mathematics
Publication Date:
Research Org.:
Lawrence Berkeley National Laboratory (LBNL), Berkeley, CA (United States)
Sponsoring Org.:
USDOE Office of Science (SC), Advanced Scientific Computing Research (ASCR)
OSTI Identifier:
1564075
Alternate Identifier(s):
OSTI ID: 1576807
Grant/Contract Number:  
AC02-05CH11231
Resource Type:
Accepted Manuscript
Journal Name:
Journal of Computational Physics
Additional Journal Information:
Journal Name: Journal of Computational Physics; Journal ID: ISSN 0021-9991
Publisher:
Elsevier
Country of Publication:
United States
Language:
English

Citation Formats

Liu, Yang, and Yang, Haizhao. A hierarchical butterfly LU preconditioner for two-dimensional electromagnetic scattering problems involving open surfaces. United States: N. p., 2019. Web. doi:10.1016/j.jcp.2019.109014.
Liu, Yang, & Yang, Haizhao. A hierarchical butterfly LU preconditioner for two-dimensional electromagnetic scattering problems involving open surfaces. United States. https://doi.org/10.1016/j.jcp.2019.109014
Liu, Yang, and Yang, Haizhao. Thu . "A hierarchical butterfly LU preconditioner for two-dimensional electromagnetic scattering problems involving open surfaces". United States. https://doi.org/10.1016/j.jcp.2019.109014. https://www.osti.gov/servlets/purl/1564075.
@article{osti_1564075,
title = {A hierarchical butterfly LU preconditioner for two-dimensional electromagnetic scattering problems involving open surfaces},
author = {Liu, Yang and Yang, Haizhao},
abstractNote = {This paper introduces a hierarchical interpolative decomposition butterfly-LU factorization (H-IDBF-LU) preconditioner for solving two-dimensional electric-field integral equations (EFIEs) in electromagnetic scattering problems of perfect electrically conducting objects with open surfaces. H-IDBF-LU leverages the interpolative decomposition butterfly factorization (IDBF) to compress dense blocks of the discretized EFIE operator to expedite its application; this compressed operator also serves as an approximate LU factorization of the EFIE operator leading to an efficient preconditioner in iterative solvers. Both the memory requirement and computational cost of the H-IDBF-LU solver scale as $O$($N$ log2 $N$) in one iteration; the total number of iterations required for a reasonably good accuracy scales as $O$(1) to O(log2$N$) in all of our numerical tests. The efficacy and accuracy of the proposed preconditioned iterative solver are demonstrated via its application to a broad range of scatterers involving up to 100 million unknowns.},
doi = {10.1016/j.jcp.2019.109014},
journal = {Journal of Computational Physics},
number = ,
volume = ,
place = {United States},
year = {Thu Oct 10 00:00:00 EDT 2019},
month = {Thu Oct 10 00:00:00 EDT 2019}
}

Journal Article:

Citation Metrics:
Cited by: 5 works
Citation information provided by
Web of Science

Save / Share:

Works referenced in this record:

Rapid solution of integral equations of scattering theory in two dimensions
journal, February 1990


Diagonal Forms of Translation Operators for the Helmholtz Equation in Three Dimensions
journal, December 1993


A wideband fast multipole method for the Helmholtz equation in three dimensions
journal, July 2006

  • Cheng, Hongwei; Crutchfield, William Y.; Gimbutas, Zydrunas
  • Journal of Computational Physics, Vol. 216, Issue 1
  • DOI: 10.1016/j.jcp.2005.12.001

The fast multipole method for the wave equation: a pedestrian prescription
journal, June 1993

  • Coifman, R.; Rokhlin, V.; Wandzura, S.
  • IEEE Antennas and Propagation Magazine, Vol. 35, Issue 3
  • DOI: 10.1109/74.250128

The Fast Multipole Method I: Error Analysis and Asymptotic Complexity
journal, January 2000


Multipole Translation Theory for the Three-Dimensional Laplace and Helmholtz Equations
journal, July 1995

  • Epton, Michael A.; Dembart, Benjamin
  • SIAM Journal on Scientific Computing, Vol. 16, Issue 4
  • DOI: 10.1137/0916051

Multilevel fast multipole algorithm for electromagnetic scattering by large complex objects
journal, January 1997

  • Song, J.
  • IEEE Transactions on Antennas and Propagation, Vol. 45, Issue 10
  • DOI: 10.1109/8.633855

Fast Directional Computation of High Frequency Boundary Integrals via Local FFTs
journal, January 2015

  • Ying, Lexing
  • Multiscale Modeling & Simulation, Vol. 13, Issue 1
  • DOI: 10.1137/140985123

A fast directional algorithm for high frequency acoustic scattering in two dimensions
journal, January 2009


Fast Directional Multilevel Algorithms for Oscillatory Kernels
journal, January 2007

  • Engquist, Björn; Ying, Lexing
  • SIAM Journal on Scientific Computing, Vol. 29, Issue 4
  • DOI: 10.1137/07068583X

Fast directional multilevel summation for oscillatory kernels based on Chebyshev interpolation
journal, February 2012

  • Messner, Matthias; Schanz, Martin; Darve, Eric
  • Journal of Computational Physics, Vol. 231, Issue 4
  • DOI: 10.1016/j.jcp.2011.09.027

Efficient Computation of Oscillatory Integrals via Adaptive Multiscale Local Fourier Bases
journal, July 2000

  • Averbuch, A.; Braverman, E.; Coifman, R.
  • Applied and Computational Harmonic Analysis, Vol. 9, Issue 1
  • DOI: 10.1006/acha.2000.0305

Fast Numerical Computations of Oscillatory Integrals Related to Acoustic Scattering, I
journal, December 1993

  • Bradie, B.; Coifman, R.; Grossmann, A.
  • Applied and Computational Harmonic Analysis, Vol. 1, Issue 1
  • DOI: 10.1006/acha.1993.1007

Sparse Approximation for Solving Integral Equations with Oscillatory Kernels
journal, January 1992

  • Canning, Francis X.
  • SIAM Journal on Scientific and Statistical Computing, Vol. 13, Issue 1
  • DOI: 10.1137/0913004

Wave atoms and sparsity of oscillatory patterns
journal, November 2007


A multilevel matrix decomposition algorithm for analyzing scattering from large structures
journal, January 1996

  • Michielssen, E.; Boag, A.
  • IEEE Transactions on Antennas and Propagation, Vol. 44, Issue 8
  • DOI: 10.1109/8.511816

An algorithm for the rapid evaluation of special function transforms
journal, March 2010

  • O'Neil, Michael; Woolfe, Franco; Rokhlin, Vladimir
  • Applied and Computational Harmonic Analysis, Vol. 28, Issue 2
  • DOI: 10.1016/j.acha.2009.08.005

Butterfly Factorization
journal, January 2015

  • Li, Yingzhou; Yang, Haizhao; Martin, Eileen R.
  • Multiscale Modeling & Simulation, Vol. 13, Issue 2
  • DOI: 10.1137/15M1007173

Analytic preconditioners for the electric field integral equation
journal, October 2004

  • Antoine, X.; Bendali, A.; Darbas, M.
  • International Journal for Numerical Methods in Engineering, Vol. 61, Issue 8
  • DOI: 10.1002/nme.1106

A Preconditioner for the Electric Field Integral Equation Based on Calderon Formulas
journal, January 2002

  • Christiansen, Snorre H.; Nédélec, Jean-Claude
  • SIAM Journal on Numerical Analysis, Vol. 40, Issue 3
  • DOI: 10.1137/S0036142901388731

Second-kind integral solvers for TE and TM problems of diffraction by open arcs: SECOND-KIND OPEN ARC INTEGRAL SOLVERS
journal, December 2012


Directional Preconditioner for 2D High Frequency Obstacle Scattering
journal, January 2015

  • Ying, Lexing
  • Multiscale Modeling & Simulation, Vol. 13, Issue 3
  • DOI: 10.1137/140985135

Fast direct solvers for integral equations in complex three-dimensional domains
journal, May 2009


Multiscale Compressed Block Decomposition for Fast Direct Solution of Method of Moments Linear System
journal, February 2011

  • Heldring, Alex; Rius, Juan M.; Tamayo, José M.
  • IEEE Transactions on Antennas and Propagation, Vol. 59, Issue 2
  • DOI: 10.1109/TAP.2010.2096385

Direct Solve of Electrically Large Integral Equations for Problem Sizes to 1 M Unknowns
journal, August 2008


An ${\cal H}^{2}$-Matrix-Based Integral-Equation Solver of Reduced Complexity and Controlled Accuracy for Solving Electrodynamic Problems
journal, October 2009


Hierarchical LU Decomposition-based Preconditioners for BEM
journal, November 2004


A Recursive Skeletonization Factorization Based on Strong Admissibility
journal, January 2017

  • Minden, Victor; Ho, Kenneth L.; Damle, Anil
  • Multiscale Modeling & Simulation, Vol. 15, Issue 2
  • DOI: 10.1137/16M1095949

An O ( N ) direct solver for integral equations on the plane
journal, March 2015

  • Corona, Eduardo; Martinsson, Per-Gunnar; Zorin, Denis
  • Applied and Computational Harmonic Analysis, Vol. 38, Issue 2
  • DOI: 10.1016/j.acha.2014.04.002

A fast direct solver for scattering problems involving elongated structures
journal, January 2007


Scattering from elongated objects: direct solution in O(N log2 N) operations
journal, January 1996

  • Michielssen, E.; Boag, A.; Chew, W. C.
  • IEE Proceedings - Microwaves, Antennas and Propagation, Vol. 143, Issue 4
  • DOI: 10.1049/ip-map:19960400

Fast Direct Solver for Essentially Convex Scatterers Using Multilevel Non-Uniform Grids
journal, August 2014

  • Brick, Yaniv; Lomakin, Vitaliy; Boag, Amir
  • IEEE Transactions on Antennas and Propagation, Vol. 62, Issue 8
  • DOI: 10.1109/TAP.2014.2327651

An HSS Matrix-Inspired Butterfly-Based Direct Solver for Analyzing Scattering From Two-Dimensional Objects
journal, January 2017

  • Liu, Yang; Guo, Han; Michielssen, Eric
  • IEEE Antennas and Wireless Propagation Letters, Vol. 16
  • DOI: 10.1109/LAWP.2016.2626786

A Butterfly-Based Direct Integral-Equation Solver Using Hierarchical LU Factorization for Analyzing Scattering From Electrically Large Conducting Objects
journal, September 2017

  • Guo, Han; Liu, Yang; Hu, Jun
  • IEEE Transactions on Antennas and Propagation, Vol. 65, Issue 9
  • DOI: 10.1109/TAP.2017.2727511

A butterfly-based direct solver using hierarchical LU factorization for Poggio-Miller-Chang-Harrington-Wu-Tsai equations
journal, April 2018

  • Guo, Han; Liu, Yang; Hu, Jun
  • Microwave and Optical Technology Letters, Vol. 60, Issue 6
  • DOI: 10.1002/mop.31166

Construction and Arithmetics of H -Matrices
journal, August 2003


Fast construction of hierarchical matrix representation from matrix–vector multiplication
journal, May 2011


A Fast Randomized Algorithm for Computing a Hierarchically Semiseparable Representation of a Matrix
journal, October 2011

  • Martinsson, P. G.
  • SIAM Journal on Matrix Analysis and Applications, Vol. 32, Issue 4
  • DOI: 10.1137/100786617

Volterra approximation for low grazing angle shadowing on smooth ocean-like surfaces
journal, November 1995

  • Holliday, D.; DeRaad, L. L.; St-Cyr, G. J.
  • IEEE Transactions on Antennas and Propagation, Vol. 43, Issue 11
  • DOI: 10.1109/8.475091

Forward-backward: a new method for computing low-grazing angle scattering
journal, May 1996

  • Holliday, D.; DeRaad, L. L.; St-Cyr, G. J.
  • IEEE Transactions on Antennas and Propagation, Vol. 44, Issue 5
  • DOI: 10.1109/8.496263

Forward-backward iterative physical optics algorithm for computing the RCS of open-ended cavities
journal, February 2005

  • Burkholder, R. J.; Lundin, T.
  • IEEE Transactions on Antennas and Propagation, Vol. 53, Issue 2
  • DOI: 10.1109/TAP.2004.841317

Sweeping preconditioner for the Helmholtz equation: Hierarchical matrix representation
journal, February 2011

  • Engquist, Björn; Ying, Lexing
  • Communications on Pure and Applied Mathematics, Vol. 64, Issue 5
  • DOI: 10.1002/cpa.20358

On the existence of nonoscillatory phase functions for second order ordinary differential equations in the high-frequency regime
journal, June 2015


An algorithm for the rapid numerical evaluation of Bessel functions of real orders and arguments
journal, May 2018


Numerical Differentiation by High Order Interpolation
journal, November 1987

  • Hoffman, Peter; Reddy, K. C.
  • SIAM Journal on Scientific and Statistical Computing, Vol. 8, Issue 6
  • DOI: 10.1137/0908079

Works referencing / citing this record:

Sparse Approximate Multifrontal Factorization with Butterfly Compression for High-Frequency Wave Equations
journal, January 2021

  • Liu, Yang; Ghysels, Pieter; Claus, Lisa
  • SIAM Journal on Scientific Computing, Vol. 43, Issue 5
  • DOI: 10.1137/20m1349667