Consider the elastic scattering of a time-harmonic wave by multiple well-separated rigid particles with smooth boundaries in two dimensions. Instead of using the complex Green's tensor of the elastic wave equation, we utilize the Helmholtz decomposition to convert the boundary value problem of the elastic wave equation into a coupled boundary value problem of the Helmholtz equation. Based on single, double, and combined layer potentials with the simpler Green's function of the Helmholtz equation, we present three different boundary integral equations for the coupled boundary value problem. The well-posedness of the new integral equations is established. Computationally, a scattering matrix based method is proposed to evaluate the elastic wave for arbitrarily shaped particles. The method uses the local expansion for the incident wave and the multipole expansion for the scattered wave. The linear system of algebraic equations is solved by GMRES with fast multipole method (FMM) acceleration. Numerical results show that the method is fast and highly accurate for solving elastic scattering problems with multiple particles.
more »
« less
Integral equation methods for the Morse-Ingard equations
We present two (a decoupled and a coupled) integral-equation-based methods for the Morse-Ingard equations subject to Neumann boundary conditions on the exterior domain. Both methods are based on second-kind integral equation (SKIE) formulations. The coupled method is well-conditioned and can achieve high accuracy. The decoupled method has lower computational cost and more flexibility in dealing with the boundary layer; however, it is prone to the ill-conditioning of the decoupling transform and cannot achieve as high accuracy as the coupled method. We show numerical examples using a Nyström method based on quadrature-by-expansion (QBX) with fast-multipole acceleration. We demonstrate the accuracy and efficiency of the solvers in both two and three dimensions with complex geometry.
more »
« less
- PAR ID:
- 10473715
- Publisher / Repository:
- Elsevier
- Date Published:
- Journal Name:
- Journal of Computational Physics
- Volume:
- 492
- Issue:
- C
- ISSN:
- 0021-9991
- Page Range / eLocation ID:
- 112416
- Subject(s) / Keyword(s):
- The Morse-Ingard Equations, Fast Multipole Method, Integral Equation Method, Quadrature-by-Expansion
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
more » « less
Summary In this paper, we propose and analyze two stabilized mixed finite element methods for the dual‐porosity‐Stokes model, which couples the free flow region and microfracture‐matrix system through four interface conditions on an interface. The first stabilized mixed finite element method is a coupled method in the traditional format. Based on the idea of partitioned time stepping, the four interface conditions, and the mass exchange terms in the dual‐porosity model, the second stabilized mixed finite element method is decoupled in two levels and allows a noniterative splitting of the coupled problem into three subproblems. Due to their superior conservation properties and convenience of the computation of flux, mixed finite element methods have been widely developed for different types of subsurface flow problems in porous media. For the mixed finite element methods developed in this article, no Lagrange multiplier is used, but an interface stabilization term with a penalty parameter is added in the temporal discretization. This stabilization term ensures the numerical stability of both the coupled and decoupled schemes. The stability and the convergence analysis are carried out for both the coupled and decoupled schemes. Three numerical experiments are provided to demonstrate the accuracy, efficiency, and applicability of the proposed methods.
-
Frequency-domain wave scattering problems that arise in acoustics and electromagnetism can be often described by the Helmholtz equation. This work presents a boundary integral equation method for the Helmholtz equation in 3-D multilayered media with many doubly periodic smooth layers. Compared with conventional quasi-periodic Green’s function method, the new method is robust at all scattering parameters. A periodizing scheme is used to decompose the solution into near-and far-field contributions. The near-field contribution uses the free-space Green’s function in an integral equation on the interface in the unit cell and its immediate eight neighbors; the far-field contribution uses proxy point sources that enclose the unit cell. A specialized high-order quadrature is developed to discretize the underlying surface integral operators to keep the number of unknowns per layer small. We achieve overall linear computational complexity in the number of layers by reducing the linear system into block tridiagonal form and then solving the system directly via block LU decomposition. The new solver is capable of handling a 100-interface structure with 961.3k unknowns to 10−5 accuracy in less than 2 hours on a desktop workstation.more » « less
-
A path integral Monte Carlo method (PIMC) based on a Feynman-Kac formula for the Laplace equation with mixed boundary conditions is proposed to solve the forward problem of the electrical impedance tomography (EIT). The forward problem is an important part of iterative algorithms of the inverse EIT problem, and the proposed PIMC provides a local solution to find the potentials and currents on individual electrodes. Improved techniques are proposed to compute with better accuracy both the local time of reflecting Brownian motions (RBMs) and the Feynman-Kac formula for mixed boundary problems of the Laplace equation. Accurate voltage-to-current maps on the electrodes of a model 3-D EIT problem with eight electrodes are obtained by solving a mixed boundary problem with the proposed PIMC method.more » « less
-
Abstract. The states of coupling between clouds andsurface or boundary layer have been investigated much more extensively formarine stratocumulus clouds than for continental low clouds, partly due tomore complex thermodynamic structures over land. A manifestation is a lackof robust remote sensing methods to identify coupled and decoupled cloudsover land. Following the idea for determining cloud coupling over the ocean,we have generalized the concept of coupling and decoupling to low cloudsover land, based on potential temperature profiles. Furthermore, by usingample measurements from lidar and a suite of surface meteorologicalinstruments at the U.S. Department of Energy's Atmospheric RadiationMeasurement Program's Southern Great Plains site from 1998 to 2019, we havedeveloped a method to simultaneously retrieve the planetary boundary layer(PBL) height (PBLH) and coupled states under cloudy conditions during thedaytime. The new lidar-based method relies on the PBLH, the liftedcondensation level, and the cloud base to diagnose the cloud coupling. Thecoupled states derived from this method are highly consistent with thosederived from radiosondes. Retrieving the PBLH under cloudy conditions, whichhas been a persistent problem in lidar remote sensing, is resolved in thisstudy. Our method can lead to high-quality retrievals of the PBLH undercloudy conditions and the determination of cloud coupling states. With thenew method, we find that coupled clouds are sensitive to changes in the PBLwith a strong diurnal cycle, whereas decoupled clouds and the PBL are weaklyrelated. Since coupled and decoupled clouds have distinct features, our newmethod offers an advanced tool to separately investigate them in climatesystems.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1016/j.jcp.2023.112416
