pith. sign in

arxiv: 2505.09495 · v1 · submitted 2025-05-14 · 🧮 math.AP

A direct imaging method for inverse scattering problem of biharmonic wave with phased and phaseless data

Tielei Zhu , Zhihao Ge This is my paper

Pith reviewed 2026-05-22 15:15 UTC · model grok-4.3

classification 🧮 math.AP MSC 35R3065N21
keywords inverse biharmonic scatteringdirect imagingreverse time migrationphaseless dataobstacle reconstructionboundary conditionsscattered field
0
0 comments X

The pith

Reverse time migration constructs imaging functions to locate obstacles from biharmonic wave data including phaseless measurements.

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

The paper proposes a direct imaging method using reverse time migration to reconstruct the shape and location of extended obstacles in biharmonic wave scattering problems. It develops imaging functions that rely on only one type of measurement data, such as the scattered field, its normal derivative, bending moment, transverse force, far-field pattern, or the phaseless total field. These functions apply to obstacles with any of four boundary conditions. Resolution analysis indicates that the functions exhibit contrast near the obstacle boundary, and numerical tests show accurate reconstruction of complex shapes with robustness to noise.

Core claim

A direct imaging method based on reverse time migration is proposed for reconstructing the extended obstacle with one of four types of boundary conditions on the obstacle. The newly developed imaging functions are constructed by utilizing merely one of various measurement data, including the scattered field, its normal derivative, the bending moment, the transverse force, its far-field and the phaseless total field. Resolution analysis demonstrates that these imaging functions have a contrast when sampling points are near or far from the boundary of the obstacle.

What carries the argument

Reverse time migration imaging functions built from a single type of biharmonic measurement data to generate a contrast map that identifies obstacle boundaries.

If this is right

  • The imaging functions achieve sufficient contrast for accurate boundary location both near and far from the obstacle.
  • The method works with phaseless total field data as well as phased measurements.
  • Numerical experiments demonstrate the algorithm's efficiency for complex scatter geometries and robustness to noise.
  • The approach applies to four types of boundary conditions on the obstacle.

Where Pith is reading between the lines

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

  • This direct method could be adapted for other types of wave scattering problems beyond biharmonic equations.
  • It may enable faster reconstructions in applications like non-destructive testing where only intensity data is available.
  • Combining data from multiple sources could further improve accuracy in challenging environments.

Load-bearing premise

The resolution analysis and numerical experiments assume that the constructed imaging functions produce sufficient contrast to accurately locate the boundary when sampling points are near or far from the obstacle, relying on the specific properties of biharmonic scattering and the chosen boundary conditions without post-hoc adjustments.

What would settle it

A numerical simulation or physical experiment with a known simple obstacle shape showing that the imaging function fails to produce a clear peak at the true boundary for phaseless data under moderate noise would falsify the reconstruction claim.

Figures

Figures reproduced from arXiv: 2505.09495 by Tielei Zhu, Zhihao Ge.

Figure 1
Figure 1. Figure 1: Example 1: The reconstruction results with no noisy data [PITH_FULL_IMAGE:figures/full_fig_p017_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: Example 2: The j−th row presents the reconstruction results of the indicator function Ij (j = 1, · · · , 11). fig2:kite 19 [PITH_FULL_IMAGE:figures/full_fig_p019_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: Example 3: The j−th row presents the reconstruction results of the indicator function Ij (j = 1, · · · , 11). fig3:obstalces Acknowledgments The research of TZ is supported by the Natural Science Foundation of Henan Province (No. 242300421677). The resaearch of ZG is supported by the National Natural Science Foundation of China (No. 12371393) and Natural Science Foundation of Henan Province (No. 2423004210… view at source ↗
read the original abstract

This paper investigates the inverse biharmonic scattering problems of identifying the shape and location of the obstacle with phased and phaseless measurement data. A direct imaging method based on reverse time migration is proposed for reconstructing the extended obstacle with one of four types of boundary conditions on the obstacle. The newly developed imaging functions are constructed by utilizing merely one of various measurement data, including the scattered field, its normal derivative, the bending moment, the transverse force, its far-field and the phaseless total field. Our resolution analysis demonstrates that these imaging functions have a contrast when sampling points are near or far from the boundary of the obstacle. Numerical experiments are further presented to show the algorithm's efficiency to accurately reconstruct complex scatter geometries and its robustness to noise.

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 paper proposes a direct imaging method based on reverse time migration for inverse biharmonic scattering problems to reconstruct the shape and location of an extended obstacle subject to one of four boundary conditions. Imaging functions are constructed from single types of measurement data (scattered field, normal derivative, bending moment, transverse force, far-field pattern, or phaseless total field). Resolution analysis is claimed to establish contrast in the imaging functions when sampling points are near or far from the obstacle boundary, with numerical experiments demonstrating reconstruction of complex geometries and robustness to noise.

Significance. If the resolution analysis holds uniformly, the work provides an efficient, non-iterative reconstruction technique for biharmonic inverse scattering, extending RTM methods to higher-order PDEs and multiple data modalities including phaseless measurements. This has potential applications in structural acoustics and plate vibration imaging. The parameter-free nature of the imaging functions and the coverage of four boundary conditions represent strengths, though the practical impact depends on the rigor of the contrast proofs and the breadth of the numerical validation.

major comments (2)
  1. [Resolution analysis] Resolution analysis section: the demonstration of contrast for the imaging functions relies on asymptotic properties of the biharmonic Green's function and boundary integral representations, but it is unclear whether the argument is uniform across all four boundary conditions (particularly the free boundary condition, where higher-order boundary terms may modify the leading-order behavior) and for the phaseless total field case; if the contrast fails to hold uniformly, the central reconstruction claim does not follow for the full range of data types advertised.
  2. [Numerical experiments] Numerical experiments section: the reported reconstructions for complex scatterers and noisy data lack quantitative metrics (e.g., Hausdorff distance or L2 boundary error) and direct comparisons to existing methods; without these, the efficiency and robustness claims cannot be fully assessed against the resolution analysis predictions.
minor comments (2)
  1. [Introduction] The introduction should explicitly list the four boundary conditions considered (clamped, simply-supported, free, etc.) rather than referring to them generically.
  2. Notation for the imaging functions should be standardized across the phased and phaseless cases to improve readability.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the careful reading and constructive comments on our manuscript. We address the major comments point by point below and will revise the manuscript to improve clarity and strengthen the presentation of results.

read point-by-point responses

  1. Referee: [Resolution analysis] Resolution analysis section: the demonstration of contrast for the imaging functions relies on asymptotic properties of the biharmonic Green's function and boundary integral representations, but it is unclear whether the argument is uniform across all four boundary conditions (particularly the free boundary condition, where higher-order boundary terms may modify the leading-order behavior) and for the phaseless total field case; if the contrast fails to hold uniformly, the central reconstruction claim does not follow for the full range of data types advertised.

    Authors: The resolution analysis relies on the leading singularity of the biharmonic fundamental solution (Green's function in free space), whose asymptotic behavior as the sampling point approaches the obstacle boundary is independent of the boundary conditions imposed on the obstacle itself. The boundary integral representations express the scattered data in terms of this fundamental solution plus a regular remainder; the contrast in each imaging function is generated by the singular term when the sampling point crosses the boundary. For the free boundary condition the higher-order boundary operators affect only the regular part of the representation and do not alter the leading singular contribution. The phaseless imaging function is constructed from the modulus of the total field; its contrast follows from the same asymptotic analysis once the phase information is recovered via the migration operator. To remove any ambiguity we will add a short remark (or subsection) that explicitly verifies uniformity for the free boundary condition and the phaseless case. revision: partial

  2. Referee: [Numerical experiments] Numerical experiments section: the reported reconstructions for complex scatterers and noisy data lack quantitative metrics (e.g., Hausdorff distance or L2 boundary error) and direct comparisons to existing methods; without these, the efficiency and robustness claims cannot be fully assessed against the resolution analysis predictions.

    Authors: We agree that quantitative error measures would allow a more precise assessment of the numerical results. In the revised version we will include tables reporting the Hausdorff distance and relative L2 boundary error for each reconstructed geometry under several noise levels. Direct comparisons with other reconstruction algorithms are more delicate because our method is non-iterative while most competing approaches are optimization-based; we will nevertheless add a brief discussion of computational cost relative to representative iterative methods and cite relevant literature on RTM-type techniques for other wave equations. These additions will be placed in a new subsection of the numerical experiments. revision: yes

Circularity Check

0 steps flagged

No significant circularity; imaging functions and resolution analysis are independently constructed from scattering data

full rationale

The paper develops RTM-based imaging functions directly from one of several measurement types (scattered field, normal derivative, bending moment, transverse force, far-field pattern, or phaseless total field) and states that resolution analysis shows contrast at the obstacle boundary for the four boundary conditions. No quoted step reduces a central claim to a fitted parameter, self-definition, or unverified self-citation chain. The derivation remains self-contained against the biharmonic Green's function and boundary conditions without the result being forced by construction or renaming.

Axiom & Free-Parameter Ledger

0 free parameters · 0 axioms · 0 invented entities

Based on the abstract alone, the paper does not introduce or fit new free parameters, axioms beyond standard scattering theory, or invented entities; it relies on existing biharmonic wave models and reverse time migration concepts.

pith-pipeline@v0.9.0 · 5653 in / 1176 out tokens · 28050 ms · 2026-05-22T15:15:04.412104+00:00 · methodology

discussion (0)

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

Reference graph

Works this paper leans on

22 extracted references · 22 canonical work pages

  1. [1]

    Bourgeois L and Hazard C 2020 On well-posedness of scattering problems in a Kirchhoff–Love infinite plateSIAM J. Appl. Math.801546-66 bourgeois_linear_2020

    work page 2020

  2. [2]

    Imaging14363-84 colton_linear_2011

    Bourgeois L and Recoquillay A 2020 The linear sampling method for Kirchhoff-Love infinite plates Inverse Probl. Imaging14363-84 colton_linear_2011

    work page 2020

  3. [3]

    Chang Y and Guo Y 2023 An optimization method for the inverse scattering problem of the bihar- monic waveCommun. Anal. Comput.1168-82 chen_reverse_2013

    work page 2023

  4. [4]

    Chen J, Chen Z and Huang G 2013 Reverse time migration for extended obstacles: Acoustic waves Inverse Problems29085005 chen_reverse_2013b

    work page 2013

  5. [5]

    Chen J, Chen Z and Huang G 2013 Reverse time migration for extended obstacles: Electromagnetic wavesInverse Problems29085006 chen_reverse_2015

    work page 2013

  6. [6]

    Chen Z and Huang G 2015 Reverse time migration for extended obstacles: Elastic wavesSci. Sin. Math.451103–14 22 chen_phaseless_2017

    work page 2015

  7. [7]

    Chen Z and Huang G 2017 Phaseless imaging by reverse time migration: Acoustic wavesNumer. Math. Theory Methods Appl.101-21 chen_boundary_2010

    work page 2017

  8. [8]

    Dong H and Li P 2024 A novel boundary integral formulation for the biharmonic wave scattering problemJ. Sci. Comput.9842 dong_uniqueness_2024

    work page 2024

  9. [9]

    Dong H and Li P 2024 Uniqueness of an inverse cavity scattering problem for the time-harmonic biharmonic wave equationInverse Problems40065011 gao_inverse_2023

    work page 2024

  10. [10]

    Gao Y, Liu H and Liu Y 2023 On an inverse problem for the plate equation with passive measurement SIAM J. Appl. Math.831196-214 guo_direct_2025

    work page 2023

  11. [11]

    Harris I, Lee H and Li P 2025 Direct sampling for recovering a clamped cavity from the biharmonic far-field dataInverse Problems4135013 Harris_direct_2022

    work page 2025

  12. [12]

    Imaging161137-62 Hohage_newton_1998

    Harris I, Nguyen D-L and Nguyen T-P 2022 Direct sampling methods for isotropic and anisotropic scatterers with point source measurementsInverse Probl. Imaging161137-62 Hohage_newton_1998

    work page 2022

  13. [13]

    Hohage T and Schormann C 1998 A Newton-type method for a transmission problem in inverse scatteringInverse Problems141207-27 hsiao_bounday_2021

    work page 1998

  14. [14]

    Ito K, Jin B and Zou J 2012 A direct sampling method to an inverse medium scattering problem Inverse Problems28025003 ji_acoustic_2019

    work page 2012

  15. [15]

    Ji X, Liu X and Zhang B 2019 Inverse acoustic scattering with phaseless far field data: Uniqueness, phase retrieval, and direct sampling methodsSIAM J. Imag. Sci.121163-89 kirsch_factorization_2008

    work page 2019

  16. [16]

    Lan T H and Nguyen D-L 2023 Stable reconstruction of anisotropic objects from near-field electro- magnetic data (arXiv:2309.12606) liu_novel_2017

    work page arXiv 2023

  17. [17]

    Liu X 2017 A novel sampling method for multiple multiscale targets from scattering amplitudes at a fixed frequencyInverse Problems33085011 Nguyen_numerical_2023

    work page 2017

  18. [18]

    Vietnam.48551-66 olver_asymptotic_1974

    Nguyen D-L and Truong T 2023 On the numerical solution to an inverse medium scattering problem Acta Math. Vietnam.48551-66 olver_asymptotic_1974

    work page 2023

  19. [19]

    Wu C and Yang J 2024 The obstacle scattering for the biharmonic equation (arXiv:2406.06126) potthast_study_2010

    work page arXiv 2024

  20. [20]

    Potthast R 2010 A study on orthogonality samplingInverse Problems26074015 yang_bayesian_2020

    work page 2010

  21. [21]

    Yang Z, Gui X, Ming J and Hu G 2020 Bayesian approach to inverse time-harmonic acoustic scat- tering with phaseless far-field dataInverse Problems3665012 zhang_approximate_2020

    work page 2020

  22. [22]

    Zhang B and Zhang H 2020 An approximate factorization method for inverse acoustic scattering with phaseless total-field dataSIAM J. Appl. Math.802271-98 23

    work page 2020