Source author record

Min Hyung Cho

Min Hyung Cho appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

math.NA Numerical Analysis

Catalog footprint

What is connected

4works

2topics

4close collaborators

Actions

Connect this record

Open graph Browse works

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

preprint2022arXiv

Quadrature by Two Expansions for Evaluating Helmholtz Layer Potentials

In this paper, a Quadrature by Two Expansions (QB2X) numerical integration technique is developed for the single and double layer potentials of the Helmholtz equation in two dimensions. The QB2X method uses both local complex Taylor expansions and plane wave type expansions to achieve a resulting representation which is numerically accurate for all target points (interior, exterior, or exactly on the boundary) inside a leaf box in the fast multipole method (FMM) hierarchical tree structure. Compared to the original Quadrature by Expansion (QBX) method, the QB2X method explicitly includes the nonlinearity from the boundary geometry in the plane wave expansions, thereby providing for higher order representations of both the boundary geometry and density functions in the integrand, with its convergence following standard FMM error analysis. Numerical results are presented to demonstrate the performance of the QB2X method for Helmholtz layer potentials and its comparison with the original QBX method for both flat and curved boundaries with various densities. The QB2X method overcomes the challenges of the original QBX method, and is better suited for efficient solutions of the Helmholtz equation with complex geometries.

preprint2019arXiv

Adapting free-space fast multipole method for layered media Green's function: algorithm and analysis

In this paper, we present a numerical algorithm for the accurate and efficient computation of the convolution of the frequency domain layered media Green's function with a given density function. Instead of compressing the convolution matrix directly as in the classical fast multipole method, fast direct solvers, and fast H-matrix algorithms, the new algorithm considers a translated form of the original matrix so that many existing building blocks from the highly optimized free-space fast multipole method can be easily adapted to the Sommerfeld integral representations of the layered media Green's function. An asymptotic analysis is performed on the Sommerfeld integrals for large orders to provide an estimate of the decay rate in the new "multipole" and "local" expansions. In order to avoid the highly oscillatory integrand in the original Sommerfeld integral representations when the source and target are close to each other, or when they are both close to the interface in the scattered field, mathematically equivalent alternative direction integral representations are introduced. The convergence of the multipole and local expansions and formulas and quadrature rules for the original and alternative direction integral representations are numerically validated.

preprint2015arXiv

Accurate and Efficient Nystrom Volume Integral Equation Method for the Maxwell equations for Multiple 3-D Scatterers

In this paper, we develop an accurate and efficient Nyström volume integral equation (VIE) method for the Maxwell equations for large number of 3-D scatterers. The Cauchy Principal Values that arise from the VIE are computed accurately using a finite size exclusion volume together with explicit correction integrals consisting of removable singularities. Also, the hyper-singular integrals are computed using interpolated quadrature formulae with tensor-product quadrature nodes for several objects, such as cubes and spheres, that are frequently encountered in the design of meta-materials . The resulting Nyström VIE method is shown to have high accuracy with a minimum number of collocation points and demonstrate $p$-convergence for computing the electromagnetic scattering of these objects. Numerical calculations of multiple scatterers of cubic and spherical shapes validate the efficiency and accuracy of the proposed method.

preprint2014arXiv

Robust fast direct integral equation solver for quasi-periodic scattering problems with a large number of layers

We present a new boundary integral formulation for time-harmonic wave diffraction from two-dimensional structures with many layers of arbitrary periodic shape, such as multilayer dielectric gratings in TM polarization. Our scheme is robust at all scattering parameters, unlike the conventional quasi-periodic Green's function method which fails whenever any of the layers approaches a Wood anomaly. We achieve this by a decomposition into near- and far-field contributions. The former uses the free-space Green's function in a second-kind integral equation on one period of the material interfaces and their immediate left and right neighbors; the latter uses proxy point sources and small least-squares solves (Schur complements) to represent the remaining contribution from distant copies. By using high-order discretization on interfaces (including those with corners), the number of unknowns per layer is kept small. We achieve overall linear complexity in the number of layers, by direct solution of the resulting block tridiagonal system. For device characterization we present an efficient method to sweep over multiple incident angles, and show a $25\times$ speedup over solving each angle independently. We solve the scattering from a 1000-layer structure with $3\times 10^5$ unknowns to 9-digit accuracy in 2.5 minutes on a desktop workstation.