Related papers: A Fast Butterfly-compressed Hadamard-Babich Integr…
Starting from the Kirchhoff-Huygens representation and Duhamel's principle of time-domain wave equations, we propose novel butterfly-compressed Hadamard integrators for self-adjoint wave equations in both time and frequency domain in an…
We present a fast and approximate multifrontal solver for large-scale sparse linear systems arising from finite-difference, finite-volume or finite-element discretization of high-frequency wave equations. The proposed solver leverages the…
In this paper, we will introduce a new heterogeneous fast multipole method (H-FMM) for 2-D Helmholtz equation in layered media. To illustrate the main algorithm ideas, we focus on the case of two and three layers in this work. The key…
A butterfly-accelerated volume integral equation (VIE) solver is proposed for fast and accurate electromagnetic (EM) analysis of scattering from heterogeneous objects. The proposed solver leverages the hierarchical off-diagonal butterfly…
This paper proposes a new Helmholtz decomposition based windowed Green function (HD-WGF) method for solving the time-harmonic elastic scattering problems on a half-space with Dirichlet boundary conditions in both 2D and 3D. The Helmholtz…
We develop an effective computational tool for simulating the scattering of 1D waves by a composite layer architected in an otherwise homogeneous medium. The layer is designed as the union of segments cut from various mother periodic media,…
A butterfly-based fast direct integral equation solver for analyzing high-frequency scattering from two-dimensional objects is presented. The solver leverages a randomized butterfly scheme to compress blocks corresponding to near- and…
A boundary integral equation method for the 3-D Helmholtz equation in multilayered media with many quasi-periodic layers is presented. Compared with conventional quasi-periodic Green's function method, the new method is robust at all…
We propose a novel Hadamard integrator for the self-adjoint time-dependent wave equation in an inhomogeneous medium. First, we create a new asymptotic series based on the Gelfand-Shilov function, dubbed Hadamard's ansatz, to approximate the…
Model-based treatment planning and exposimetry for high-intensity focused ultrasound (HIFU) requires the numerical simulation of nonlinear ultrasound propagation through heterogeneous and absorbing media. This is a computationally demanding…
Hybrid quantum-classical algorithms have been proposed to circumvent noise limitations in quantum computers. Such algorithms delegate only a calculation of the expectation value to the quantum computer. Among them, the Variational Quantum…
Reductions of the self-consistent mean field theory model of amphiphilic molecules in solvent can lead to a singular family of functionalized Cahn-Hilliard energies. We modify these energies, mollifying the singularities to stabilize the…
This paper is concerned with the asymptotic description of high-frequency waves in locally periodic media. A key issue is that the Bloch-dispersion curves vary with the local microstructure, giving rise to hidden singularities associated…
In this paper, a fast multipole method (FMM) is proposed to compute long-range interactions of wave sources embedded in 3-D layered media. The layered media Green's function for the Helmholtz equation, which satisfies the transmission…
Efficient numerical solution of the acoustic Helmholtz equation in heterogeneous media remains challenging, particularly for large-scale problems with spatially-varying density - a limitation that restricts applications in biomedical…
Boundary integral equations lead to dense system matrices when discretized, yet they are data-sparse. Using the $\mathcal{H}$-matrix format, this sparsity is exploited to achieve $\mathcal{O}(N\log N)$ complexity for storage and…
Developing efficient path integral (PI) methods for atomistic simulations of vibrational spectra in heterogeneous condensed phases and interfaces has long been a challenging task. Here, we present the h-CMD method, short for hybrid centroid…
Hybrid spectral-spatial representations are introduced to rapidly calculate periodic scalar and dyadic Green's functions of the Helmholtz equation for 2D and 3D configurations with a 1D (linear) periodicity. The presented schemes work…
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…
In this paper, we present a fast boundary integral method accelerated by the fast multipole method (FMM) for acoustic wave scattering governed by the scalar Helmholtz equation in multi-layered two-dimensional media. Multiple scatterers are…