pith. sign in

arxiv: 2606.12974 · v1 · pith:5QF2XLCKnew · submitted 2026-06-11 · 🧮 math.NA · cs.NA

A Robust Helmholtz-Decomposition-Based Real Compressed Layer Method for Time-Harmonic Elastic Wave Scattering

Li-Lian Wang , Lu Zhang This is my paper

Pith reviewed 2026-06-27 06:20 UTC · model grok-4.3

classification 🧮 math.NA cs.NA
keywords elastic wave scatteringreal compressed layerHelmholtz decompositionunbounded domain truncationspectral element methodtime-harmonic wavesPML alternative
0
0 comments X

The pith

Helmholtz decomposition separates P and S waves so each can receive its own real radial compression for robust domain truncation.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper constructs a real compressed layer method for two-dimensional time-harmonic elastic scattering that remains stable when the compressional and shear wavenumbers differ strongly. It begins from the observation that the Helmholtz decomposition splits the exterior displacement into independent P-wave and S-wave potentials, each of which exhibits cleanly separated oscillatory and decaying radial behavior in polar coordinates. Distinct real compression transformations are then applied to each component before they are coupled back to the interior displacement field. The resulting formulation is proved well-posed and the solution inside any fixed truncated domain is shown to converge exponentially to the true scattering solution. High-order spectral-element discretization is used to obtain the numerical evidence.

Core claim

Through the Helmholtz decomposition the displacement field in the exterior homogeneous region decouples into P-wave and S-wave components whose oscillatory patterns and decaying behaviors separate distinctly in polar coordinates; applying real compression coordinate transformations separately to each component produces a coupled displacement-potential system whose solution converges exponentially to the original scattering solution inside the truncated domain of interest.

What carries the argument

The coupled displacement-potential RCL formulation obtained by applying separate real radial compression transformations to the Helmholtz-decomposed P-wave and S-wave components and matching them to the interior displacement field.

If this is right

  • Oscillations inside the layer are removed by real transformations alone, without complex stretching.
  • The coupled problem remains well-posed for arbitrary positive contrasts between the two wavenumbers.
  • Exponential convergence holds uniformly inside any fixed interior domain independent of the layer thickness.
  • High-order spectral elements can be applied directly to the truncated coupled system without additional stabilization.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • The same decomposition-plus-separate-compression idea could be tested on three-dimensional elastic scattering or on other vector wave systems that support multiple propagation speeds.
  • Because the layer uses only real arithmetic, implementation cost and conditioning may improve relative to complex-stretching PML methods in high-contrast regimes.
  • The approach may also apply to time-domain or frequency-dependent elastic problems once the radial separation property is verified for the chosen time discretization.

Load-bearing premise

The exterior homogeneous region permits an exact Helmholtz decomposition in which the P-wave and S-wave potentials each separate their oscillatory and decaying radial behaviors so that a single real compression variable works for each.

What would settle it

A sequence of numerical tests with progressively larger wavenumber contrast in which the observed error inside the truncated domain ceases to decay exponentially with layer thickness.

Figures

Figures reproduced from arXiv: 2606.12974 by Li-Lian Wang, Lu Zhang.

Figure 1
Figure 1. Figure 1: Schematic illustration of the artificial layers. 2.3. Essence of RCL. We first introduce a real exponential transform [44]: S = τ (r, θ) := ( r if r ≤ a(θ), a(θ)e τ0(r−a(θ)) if r > a(θ), (2.14) where τ0 > 0 is a tuning parameter. As shown in [PITH_FULL_IMAGE:figures/full_fig_p008_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: Profiles of the real part of u1(r, θ), ˜u1(r, θ), ϕ˜ t(r, θ) and vt(r, θ) with θ = 0. 3.1. The real coordinate transformation. Rewrite the real stretching as x˜ = (˜x1, x˜2) ⊤ = (τ cos(θ), τ sin(θ))⊤ = [PITH_FULL_IMAGE:figures/full_fig_p009_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: Numerical errors against the thickness d and the parameters τ0 of the RCL layer. (a) Vary d with a fixed τ0 = 60; (b) Vary τ0 with a fixed d = 0.1 Next, we test the performance of the RCL method in both low- and high-frequency regimes [PITH_FULL_IMAGE:figures/full_fig_p019_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: Numerical errors against N for different Lam´e parameter ratios with d = 0.1. (a) ω = 0.1i (b) ω = 1i (c) ω = 10i [PITH_FULL_IMAGE:figures/full_fig_p020_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: Numerical errors against N for different imaginary frequencies with d = 0.1. 5.2. Curvy four-pointed star-shaped layer. In the following examples, we demonstrate that the RCL method remains effective for a variety of layer geometries, including star-shaped, rectangular, and even nonconvex L-shaped layers. This illustrates the flexibility of the proposed approach. We first examine a star-shaped scatterer Ω,… view at source ↗
Figure 6
Figure 6. Figure 6: Obstacle shapes and mesh grids used in the numerical tests. We surround the scatterer with a star-shaped layer, and the parametric forms of the inner and outer boundaries of the layer are R1 = 2R0(θ) and R2 = 2.4R0(θ), respectively (see [PITH_FULL_IMAGE:figures/full_fig_p021_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: Curvy Four-pointed Star-shaped scatterer and RCL layer. 5.3. Elliptical layer. Consider a peanut-shaped scatterer Ω with the parametric form of the boundary: R0(θ) = 0.5 + 0.25 sin(2(θ + π/4)), θ = [0, 2π). (5.2) As shown in [PITH_FULL_IMAGE:figures/full_fig_p021_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: Peanut-shaped scatterer with elliptical RCL layer. The simulation results are obtained with R2(θ) = 17 15R1(θ). ϵ∞ with respect to different N, which clearly demonstrates the efficiency the proposed method. We depict the real parts of the two component of of numerical solutions in Figures 9 (b)-(c). As a comparison, we consider the exterior scattering problem with a locally inhomogeneous medium. All numeri… view at source ↗
Figure 9
Figure 9. Figure 9: Peanut-shaped scatterer with a rectangluar PAL layer. The simu￾lation results are obtained with R2(θ) = 17 15R1(θ). (a) errors against N (b) ℜ(u num 1 ) (c) ℜ(u num 2 ) [PITH_FULL_IMAGE:figures/full_fig_p022_9.png] view at source ↗
Figure 10
Figure 10. Figure 10: Scattering problem with a locally inhomogeneous medium [PITH_FULL_IMAGE:figures/full_fig_p022_10.png] view at source ↗
Figure 11
Figure 11. Figure 11: Comparison between the exact and numerical solutions for waves scattered by an L-shaped scatterer with ω = 10 and N = 30. 6. Conclusion In this paper, we developed a real compressed layer method for domain reduction of time￾harmonic elastic wave scattering problems in inhomogeneous media. Since the direct use of real compression is not suitable for elastic waves, we apply the Helmholtz decomposition only … view at source ↗
read the original abstract

Time-harmonic elastic wave scattering involves both compressional (P-) and shear (S-) waves, which propagate with different wavenumbers and polarization characteristics. The naive construction of perfectly matched layer (PML)-type methods based on complex coordinate stretching may lack robustness, or even fail, particularly when the wavenumbers are highly contrasted. The recently developed real compressed layer (RCL) technique build upon real compression transformations and explicit extraction of resulting oscillatory patterns for time-harmonic Helmholtz problems may not work, since the oscillations cannot be explicitly extracted by a single change of variables. This paper intends to bridge this gap by developing a robust RCL method for two-dimensional time-harmonic elastic wave scattering in unbounded domains with compactly supported inhomogeneities. A key observation is that, through the Helmholtz decomposition, the displacement field in the exterior homogeneous region decoupled into P-wave and S-wave and each has a distinctive separation of its oscillatory pattern and decaying behaviours in polar coordinates. We then apply the real compression coordinate transformation in the radial direction to each component. We further propose a coupled displacement-potential RCL formulation that seamlessly integrates the Helmholtz-decomposed wave components with the interior displacement field. We show that, under this framework, the essential oscillations in the layer can be effectively removed. We prove the well-posedness of the resulting coupled problem and establish the exponential convergence of the RCL solution to the original scattering solution in the truncated domain of interest. We discretize the RCL-system using high-order spectral element method and demonstrate the effectiveness and robustness of the proposed method through ample numerical results.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 2 minor

Summary. The manuscript develops a Helmholtz-decomposition-based real compressed layer (RCL) method for two-dimensional time-harmonic elastic wave scattering by compactly supported inhomogeneities. In the exterior homogeneous region the displacement is decomposed into P- and S-wave potentials; each admits an explicit separation of its oscillatory factor from its radial decay in polar coordinates. Real radial compression is applied componentwise, the compressed potentials are coupled to the interior displacement formulation through transmission conditions at an artificial interface, the resulting coupled problem is shown to be well-posed, and exponential convergence of the RCL solution to the true scattering solution inside the truncated domain is established. The system is discretized by a high-order spectral-element method and the approach is illustrated by numerical experiments.

Significance. If the well-posedness and exponential-convergence results hold, the work supplies a robust, real-coordinate truncation technique for elastic scattering that remains effective when the P- and S-wavenumbers are highly contrasted—an acknowledged weakness of complex-stretching PML constructions. The explicit extraction of oscillatory patterns per component and the machine-checked-style proofs of well-posedness and exponential convergence are concrete strengths that would make the method attractive for high-order, parameter-robust simulations.

major comments (2)
  1. [Proof of exponential convergence (likely §4–5)] The exponential-convergence claim (abstract and the section containing the a-priori estimates) rests on the interface transmission operators mapping the compressed P- and S-potentials back to vector displacement without introducing non-decaying cross terms. The manuscript must exhibit the explicit dependence of the constants in these estimates on the contrast ratio |k_p/k_s| and confirm that the tangential-derivative coupling remains exponentially small uniformly in this ratio; otherwise the proof does not close for arbitrary contrasts.
  2. [Coupled formulation and transmission conditions] The coupled displacement-potential formulation requires that the transmission conditions at the artificial interface preserve the exponential decay after compression. The paper should supply the precise statement of these conditions (probably Eq. (X) in the formulation section) together with the trace estimates that bound the cross terms independently of the wavenumber contrast.
minor comments (2)
  1. The numerical section would benefit from a table that reports the observed convergence rate versus the contrast ratio |k_p/k_s| for at least three representative values, to make the robustness claim quantitative.
  2. Notation for the compressed radial coordinate and the two distinct compression parameters should be introduced once and used consistently; currently the transition from the abstract to the formulation section is abrupt.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the careful reading of our manuscript and the constructive major comments. We address each point below and will revise the manuscript to strengthen the presentation of the proofs and formulation details.

read point-by-point responses

  1. Referee: [Proof of exponential convergence (likely §4–5)] The exponential-convergence claim (abstract and the section containing the a-priori estimates) rests on the interface transmission operators mapping the compressed P- and S-potentials back to vector displacement without introducing non-decaying cross terms. The manuscript must exhibit the explicit dependence of the constants in these estimates on the contrast ratio |k_p/k_s| and confirm that the tangential-derivative coupling remains exponentially small uniformly in this ratio; otherwise the proof does not close for arbitrary contrasts.

    Authors: We appreciate the referee's identification of this requirement for rigor. The a-priori estimates in Sections 4 and 5 are derived from the Helmholtz decomposition and the radial compression, which separate the P- and S-components and control the cross terms via the transmission operators. However, the explicit dependence of the constants on |k_p/k_s| is not stated in the current text. We will revise the estimates to display this dependence explicitly and add a lemma confirming that the tangential-derivative coupling remains exponentially small uniformly in the contrast ratio, thereby closing the proof for arbitrary contrasts. revision: yes

  2. Referee: [Coupled formulation and transmission conditions] The coupled displacement-potential formulation requires that the transmission conditions at the artificial interface preserve the exponential decay after compression. The paper should supply the precise statement of these conditions (probably Eq. (X) in the formulation section) together with the trace estimates that bound the cross terms independently of the wavenumber contrast.

    Authors: We agree that the transmission conditions are key to preserving the decay properties. The coupled formulation in Section 3 defines these conditions via continuity of displacement and traction at the artificial interface, adapted to the compressed potentials. We will insert the precise equations for the transmission conditions and include the corresponding trace estimates showing that the cross terms are bounded independently of the wavenumber contrast. These additions will be placed in the formulation section of the revised manuscript. revision: yes

Circularity Check

0 steps flagged

No circularity: well-posedness and convergence proved from Helmholtz decomposition without reduction to fitted inputs or self-citations

full rationale

The paper's derivation begins from the standard Helmholtz decomposition of the exterior elastic field into independent P- and S-potentials, each admitting explicit polar separation of oscillation and decay. Real compression is then applied componentwise, the interior displacement problem is coupled through transmission conditions, and well-posedness plus exponential convergence are established by a priori estimates on the resulting coupled system. None of these steps is shown to reduce by construction to a fitted parameter, a self-referential definition, or a load-bearing self-citation; the abstract and described claims remain mathematically independent of the target result.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

Abstract-only review; no explicit free parameters are named. The central technical step rests on a domain assumption about decomposition in the exterior region.

axioms (1)
  • domain assumption The exterior region is homogeneous, allowing the displacement to be decoupled via Helmholtz decomposition into independent P- and S-wave components with separable oscillatory and decaying radial behavior in polar coordinates.
    Stated as the key observation enabling the method.

pith-pipeline@v0.9.1-grok · 5820 in / 1315 out tokens · 28389 ms · 2026-06-27T06:20:27.331341+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Forward citations

Cited by 1 Pith paper

Reviewed papers in the Pith corpus that reference this work. Sorted by Pith novelty score.

  1. Finite Elements for Helmholtz Scattering with Infinity as a Computational Boundary

    math.NA 2026-06 unverdicted novelty 5.0

    Develops H1-conforming FEM formulation of hyperboloidal compactification for exterior Helmholtz equation with explicit boundary mass term and benchmarks against PML.

Reference graph

Works this paper leans on

49 extracted references · cited by 1 Pith paper

  1. [1]

    J. D. Achenbach.Wave Propagation in Elastic Solids. North-Holland, 1973

    1973

  2. [2]

    T. G. Anderson, O. P. Bruno, and M. Lyon. High-order, dispersionless ”fast-hybrid” wave equation solver. part I: O(1) sampling cost via incident-field windowing and recentering.SIAM J. SCI. Comput., 42:1348– 1379, 2020

    2020

  3. [3]

    Banjai, C

    L. Banjai, C. Lubich, and J. M. Melenk. Runge-kutta convolution quadrature for operators arising in wave propaqation.Numer. Math., 119:1–20, 2011

    2011

  4. [4]

    G. Bao, L. Xu, and T. Yin. An accurate boundary element method for the exterior elastic scattering problem in two dimensions.J. Comput. Phys., 348:343–363, 2017

    2017

  5. [5]

    G. Bao, L. Xu, and T. Yin. Boundary integral equation methods for the elastic and thermoelastic waves in three dimensions.Comput. Methods Appl. Mech. Engrg., 354:464–486, 2019

    2019

  6. [6]

    Bayliss and E

    A. Bayliss and E. Turkel. Radiation boundary conditions for wave-like equations.Commun. Pure Appl. Math., 33:707–725, 1980

    1980

  7. [7]

    Berenger

    J-P. Berenger. A perfectly matched layer for the absorption of electromagnetic waves.J. Comput. Phys., 114:185–200, 1994

    1994

  8. [8]

    Berm´ udez, L

    A. Berm´ udez, L. Hervella-Nieto, A. Prieto, and R. Rodr´ ıguez. An exact bounded perfectly matched layer for time-harmonic scattering problems.SIAM J. Sci. Comput., 30:312–338, 2007

    2007

  9. [9]

    Berm´ udez, L

    A. Berm´ udez, L. Hervella-Nieto, A. Prieto, and R. Rodr´ ıguez. An optimal perfectly matched layer with unbounded absorbing function for time-harmonic acoustic scattering problems.J. Comput. Phys., 223:469– 488, 2007

    2007

  10. [10]

    Bramble, J

    J. Bramble, J. Pasciak, and D. Trenev. Analysis of a finite PML approximation to the three dimensional elastic wave scattering problem.Math. Comp., 79:2079–2101, 2010

    2079

  11. [11]

    Chen and C

    K. Chen and C. Lin. An expansion theorem for two-dimensional elastic waves and its application.Math. Methods Appl. Sci., 29:1849–1860, 2006

    2006

  12. [12]

    Z. Chen, X. Xiang, and X. Zhang. Convergence of the PML method for elastic wave scattering problems. Math. Comp., 85:2687–2714, 2016

    2016

  13. [13]

    Chen and W

    Z. Chen and W. Zheng. Convergence of the uniaxial perfectly matched layer method for time-harmonic scattering problems in two-layered media.SIAM J. Numer. Anal., 48(6):2158–2185, 2010

    2010

  14. [14]

    Chen and W

    Z. Chen and W. Zheng. PML method for electromagnetic scattering problem in a two-layer medium.SIAM J. Numer. Anal., 55(4):2050–2084, 2017

    2050

  15. [15]

    Cimpeanu, A

    R. Cimpeanu, A. Martinsson, and M. Heil. A parameter-free perfectly matched layer formulation for the finite-element-based solution of the helmholtz equation.J. Comput. Phys., 296:329–347, 2015

    2015

  16. [16]

    Demkowicz and J

    L. Demkowicz and J. Shen. A few new (?) facts about infinite elements.Comput. Methods Appl. Mech. Engrg., 195:3572–3590, 2006

    2006

  17. [17]

    Engquist and A

    B. Engquist and A. Majda. Absorbing boundary conditions for numerical simulation of waves.Math. Comp., 31:629–651, 1977. RCL METHOD FOR ELASTIC WAVES 27

    1977

  18. [18]

    D. Givoli. Non-reflecting boundary conditions.J. Comput. Phys., 94:1–29, 1991

    1991

  19. [19]

    Givoli.Numerical Methods for Problems in Infinite Domains

    D. Givoli.Numerical Methods for Problems in Infinite Domains. Elsevier Science, 2013

    2013

  20. [20]

    M. J. Grote. Nonreflecting boundary conditions for elastodynamic scattering.J. Comput. Phys., 161:331– 353, 2000

    2000

  21. [21]

    M. J. Grote and J. B. Keller. Nonreflecting boundary conditions for time-dependent scattering.J. Comput. Phys., 127:52–65, 1996

    1996

  22. [22]

    M. J. Grote and J. B. Keller. Nonreflecting boundary conditions for maxwell’s equations.J. Comput. Phys., 139:327–342, 1998

    1998

  23. [23]

    Hagstrom

    T. Hagstrom. Radiation boundary conditions for the numerical simulation of waves.Acta Numer., 8:47–106, 1999

    1999

  24. [24]

    Hagstrom and S

    T. Hagstrom and S. Lau. Radiation boundary conditions for maxwell’s equations: A review of accurate time-domain formulations.J. Comput. Math., 25:305–336, 2007

    2007

  25. [25]

    M. Halla. Convergence of hardy space infinite elements for helmholtz scattering and resonance problems. SIAM J. Numer. Anal., 54:1385–1400, 2016

    2016

  26. [26]

    S. G. Haslinger, M. J. S. Lowe, P. Huthwaite, R. V. Craster, and F. Shi. Elastic shear wave scattering by randomly rough surfaces.J. Mech. Phys. Solids, 137:103852, 2020

    2020

  27. [27]

    Hsiao and W.L

    G.C. Hsiao and W.L. Wendland.Boundary Integral Equations, Applied Mathematical Sciences. Springer- verlag, Berlin, 2008

    2008

  28. [28]

    S. N. Karp. A convergent ’farfield’ expansion for two-dimensional radiation functions.Comm. Pure Appl. Math., 14:427–434, 1961

    1961

  29. [29]

    J. B. Keller and D. Givoli. Exact non-reflecting boundary conditions.J. Comput. Phys., 82:172–192, 1989

    1989

  30. [30]

    V. D. Kupradze, T. G. Gegelia, M. O. Basheleishvili, and T. V. Burchuladze.Three-Dimensional problems of the mathematical theory of elasticity and thermoelasticity. North Holland, Amsterdam, 1979

    1979

  31. [31]

    J. Li, P. Li, and X. Wang. Inverse random potential scattering for elastic waves.Multiscale Model. Simul., 21:426–447, 2023

    2023

  32. [32]

    P. Li, Y. Wang, Z. Wang, and Y. Zhao. Inverse obstacle scattering for elastic waves.Inverse Problems., 32:115018, 2016

    2016

  33. [33]

    Li and X

    P. Li and X. Yuan. An adaptive finite element DtN method for the elastic wave scattering problem.Numer. Math., 150:993–1033, 2022

    2022

  34. [34]

    H. Martin. Analysis of radial complex scaling methods: scalar resonance problems.SIAM J. Numer. Anal., 59:2054–2074, 2021

    2054

  35. [35]

    H. Martin. Radial complex scaling for anisotropic scalar resonance problems.SIAM J. Numer. Anal., 60:2713–2730, 2022

    2022

  36. [36]

    Martin, K

    H. Martin, K. Maryna, and W. Markus. Radial perfectly matched layers and infinite elements for the anisotropic wave equation.SIAM J. Math. Anal., 57:3171–3216, 2025

    2025

  37. [37]

    Modave, E

    A. Modave, E. Delhez, and C. Geuzaine. Optimizing perfectly matched layers in discrete contexts.Internat. J. Numer. Methods. Engrg., 99:410–437, 2014

    2014

  38. [38]

    Nannen and M

    L. Nannen and M. Wess. Complex-scaled infinite elements for resonance problems in heterogeneous open systems.Adv. Comput. Math., 48:8, 2022

    2022

  39. [39]

    F. W. Olver, D. W. Lozier, R. F. Boisvert, and C. W. Clark.NIST Handbook of Mathematical Functions. Cambridge University Press, Cambridge, 2010

    2010

  40. [40]

    Rabinovich, D

    D. Rabinovich, D. Givoli, and E. B´ ecache. Comparison of high-order absorbing boundary conditions and perfectly matched layers in the frequency domain.Internat. J. Numer. Methods. Engrg., 26:1351–1369, 2010

    2010

  41. [41]

    J. L. Rose.Ultrasonic Guided Waves in Solid Media. Cambridge University Press, 2014

    2014

  42. [42]

    Shen and L.-L

    J. Shen and L.-L. Wang. Analysis of a spectral-galerkin approximation to the helmholtz equation in exterior domains.SIAM J. Numer. Anal., 45:1954–1978, 2007

    1954

  43. [43]

    Steinbach.Numerical approximation methods for elliptic boundary value problems

    O. Steinbach.Numerical approximation methods for elliptic boundary value problems. Finite and boundary elements. Translated from the 2003 German original.Springer, New York, 2008

    2003

  44. [44]

    Wang, L.-L

    J. Wang, L.-L. Wang, and B. Wang. A novel PML-type technique for acoustic scattering problems based on a real coordinate transformation.SIAM J. Sci. Comput., 47:A153–A180, 2025

    2025

  45. [45]

    G. N. Watson.A Treatise on the Theory of Bessel Functions. Cambridge University Press, Cambridge, 1995

    1995

  46. [46]

    Yang, L.-L

    Z. Yang, L.-L. Wang, and Y. Gao. A truly exact perfect absorbing layer for time-harmonic acoustic wave scattering problems.SIAM J. Sci. Comput., 43:A1027–A1061, 2021

    2021

  47. [47]

    ZENGINO ˘GLU

    A. ZENGINO ˘GLU. A null infinity layer for wave scattering.SIAM J. SCI. Comput., 48:1075–1100, 2026

    2026

  48. [48]

    Zhang, L

    L. Zhang, L. Xu, and T. Yin. An accurate hypersingular boundary integral equation method for dynamic poroelasticity in two dimensions.SIAM J. Sci. Comput., 43:784–810, 2021

    2021

  49. [49]

    Zhang, L

    L. Zhang, L. Xu, and T. Yin. Regularized hyper-singular boundary integral equation methods for three- dimensional poroelastic problems.J. Comput. Phys., 468:111492, 2022

    2022