pith. sign in

arxiv: 2510.14088 · v2 · pith:DXM4Y7I2new · submitted 2025-10-15 · 🧮 math.NA · cs.NA

Geometric local parameterization for solving Hele-Shaw problems with surface tension

Zengyan Zhang , Wenrui Hao , John Harlim This is my paper

Pith reviewed 2026-05-21 20:23 UTC · model grok-4.3

classification 🧮 math.NA cs.NA
keywords Hele-Shaw problemfree boundarypoint cloudlocal parameterizationsurface tensionconvergence analysisnumerical discretizationmoving boundary
0
0 comments X

The pith

A geometric local parameterization from point clouds enables high-order discretization and convergence analysis for the Hele-Shaw free boundary problem with surface tension.

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

The paper presents a computational approach to the two-dimensional Hele-Shaw problem that tracks the moving boundary using only a set of points instead of a fixed coordinate system. Local approximations built from these points compute the necessary geometric information such as curvature and allow the boundary integral equation to be discretized with an exact treatment of singular parts. Rigorous analysis establishes that the method is consistent and stable when the boundary remains smooth and the points stay evenly spaced, producing error estimates that scale with the point spacing and the order of the integration rule. Tests show the expected accuracy rates and demonstrate evolution of various initial shapes into circles under the smoothing effect of surface tension.

Core claim

The paper establishes a spatial discretization for the Hele-Shaw problem with surface tension based on local geometric parameterization of the boundary from a uniformly sampled point cloud. This parameterization supplies high-order approximations to curvature and other quantities and supports an analytical formula for the singular integrals appearing in the governing equation. Under assumptions of sufficient smoothness and uniform sampling, consistency and stability are proved, and an explicit error bound is derived in terms of the point cloud resolution, boundary smoothness, and quadrature order.

What carries the argument

The geometric local parameterization, which constructs local coordinate systems from nearby point data on the boundary to approximate differential quantities like curvature without requiring a global mapping of the entire interface.

If this is right

  • The numerical error is bounded by a quantity that decreases as the distance between neighboring boundary points is reduced.
  • High-order spatial accuracy is obtained when higher-order quadrature rules are employed for the integral discretization.
  • Long-time simulations remain reliable provided the point distribution on the boundary stays uniform.
  • Interfaces with anisotropic surface tension evolve toward circular equilibrium shapes as expected from the physics.

Where Pith is reading between the lines

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

  • This local-chart technique could be adapted to other free-boundary problems in fluid mechanics that involve curvature-driven motion.
  • Combining the point-cloud representation with dynamic resampling might extend the method to cases where the boundary develops regions of high curvature.
  • Similar local parameterization ideas may apply to three-dimensional versions of the Hele-Shaw problem or to related moving-boundary models in materials science.

Load-bearing premise

The moving boundary must remain sufficiently smooth and the sampled points must stay uniformly distributed for the local approximations to retain their designed accuracy order.

What would settle it

Numerical experiments in which the boundary develops a singularity or the point cloud becomes clustered, resulting in a loss of the predicted convergence rate.

Figures

Figures reproduced from arXiv: 2510.14088 by John Harlim, Wenrui Hao, Zengyan Zhang.

Figure 1
Figure 1. Figure 1: Root-Mean-Square errors of the curvature as defined in [PITH_FULL_IMAGE:figures/full_fig_p006_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: Mesh refinement tests. (a) Evolution of the circle boundary over time. (b) Error for the circle radius over time [PITH_FULL_IMAGE:figures/full_fig_p016_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: The evolution dynamics for the circular case. In (a), the profiles of [PITH_FULL_IMAGE:figures/full_fig_p016_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: The evolution dynamics for perturbed circles with different [PITH_FULL_IMAGE:figures/full_fig_p018_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: The evolution dynamics for smooth closed curves. In (a) and (c), the profiles of [PITH_FULL_IMAGE:figures/full_fig_p019_5.png] view at source ↗
read the original abstract

In this work, we introduce a novel computational framework for solving the two-dimensional Hele-Shaw free boundary problem with surface tension. The moving boundary is represented by point clouds, eliminating the need for a global parameterization. Our approach leverages Generalized Moving Least Squares (GMLS) to construct local geometric charts, enabling high-order approximations of geometric quantities such as curvature directly from the point cloud data. This local parameterization is systematically employed to discretize the governing boundary integral equation, including an analytical formula of the singular integrals. We provide a rigorous convergence analysis for the proposed spatial discretization, establishing consistency and stability under certain conditions. The resulting error bound is derived in terms of the size of the uniformly sampled point cloud data on the moving boundary, the smoothness of the boundary, and the order of the numerical quadrature rule. Numerical experiments confirm the theoretical findings, demonstrating high-order spatial convergence and the expected temporal convergence rates. The method's effectiveness is further illustrated through simulations of complex initial shapes, including interfaces driven by anisotropic surface tension, which correctly evolve towards circular equilibrium states under the influence of surface tension, highlighting the versatility of the method for complex geometry-dependent interface dynamics.

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

1 major / 2 minor

Summary. The paper introduces a computational framework for the two-dimensional Hele-Shaw free-boundary problem with surface tension. The moving interface is represented by a point cloud rather than a global parameterization; Generalized Moving Least Squares (GMLS) is used to construct local geometric charts that furnish high-order approximations to curvature and other geometric quantities. These local charts are then employed to discretize the boundary-integral formulation of the problem, including an analytic treatment of the singular integrals. A rigorous consistency-and-stability analysis is supplied for the spatial discretization, producing an error bound controlled by the point-cloud spacing h, the smoothness of the interface, and the quadrature order. Numerical experiments are reported to confirm the predicted high-order spatial and temporal rates, and the method is demonstrated on complex initial shapes, including those driven by anisotropic surface tension, that relax toward circular equilibria.

Significance. If the global error estimate remains valid under the moving-boundary evolution, the work supplies a useful alternative to traditional parameterized or level-set approaches for free-boundary problems that must accommodate complex or topology-changing geometries. The combination of GMLS local charts with an analytic singular-integral treatment and an attempt at a rigorous error analysis constitutes a clear technical contribution to numerical methods for interface dynamics.

major comments (1)
  1. [Convergence analysis section] Convergence analysis (the section containing the consistency/stability theorem and the statement of the error bound): the derived bound is stated to hold for a uniformly sampled point cloud whose spacing is controlled by h. The normal velocity of the Hele-Shaw flow can stretch or cluster points, yet the manuscript supplies no explicit estimate showing that any redistribution or interpolation step preserves the O(h^k) consistency error already present in the spatial discretization. Without such a bound, the global error after O(1/h) time steps is not guaranteed to remain of the claimed order.
minor comments (2)
  1. [Numerical experiments / Implementation details] Add a short paragraph or subsection describing the precise redistribution algorithm (including any tolerance or frequency) used to maintain approximate uniformity of the point cloud.
  2. [Convergence analysis section] Clarify the precise Sobolev or Hölder class assumed for the interface in the error analysis and state the quadrature order explicitly in the theorem statement.

Simulated Author's Rebuttal

1 responses · 0 unresolved

We thank the referee for the careful reading of the manuscript and the constructive comment on the convergence analysis. The observation regarding the preservation of the spatial error bound under interface evolution is well taken, and we address it directly below.

read point-by-point responses

  1. Referee: [Convergence analysis section] Convergence analysis (the section containing the consistency/stability theorem and the statement of the error bound): the derived bound is stated to hold for a uniformly sampled point cloud whose spacing is controlled by h. The normal velocity of the Hele-Shaw flow can stretch or cluster points, yet the manuscript supplies no explicit estimate showing that any redistribution or interpolation step preserves the O(h^k) consistency error already present in the spatial discretization. Without such a bound, the global error after O(1/h) time steps is not guaranteed to remain of the claimed order.

    Authors: We agree that the consistency/stability theorem is stated for a fixed, uniformly spaced point cloud. In the numerical experiments the point cloud is periodically redistributed to maintain quasi-uniform spacing, but the manuscript does not supply an explicit a-priori bound showing that this redistribution preserves the O(h^k) consistency error. We will revise the convergence-analysis section to include a short additional remark (or lemma) that quantifies the redistribution error. Under the assumption that the re-sampling operator is accurate to order k+1 and that the new point spacing remains O(h), the perturbation to the GMLS operators remains O(h^k). With this estimate the spatial error bound continues to hold at each time step; combined with the existing stability result, the global error after O(1/h) steps remains of the claimed order. The revised text will make this dependence explicit. revision: yes

Circularity Check

0 steps flagged

No circularity in derivation chain; analysis is conditional but self-contained

full rationale

The paper introduces a GMLS-based local parameterization for discretizing the Hele-Shaw boundary integral equation on point clouds and states a convergence analysis establishing consistency and stability, with the error bound expressed in terms of point-cloud size h, boundary smoothness, and quadrature order. No quoted equations, self-citations, or fitted parameters reduce any claimed result to its own inputs by construction; the bound is derived under explicitly stated standing assumptions of persistent uniform sampling and smoothness rather than being tautological or statistically forced. The derivation chain therefore remains independent of the target quantities and does not exhibit any of the enumerated circular patterns.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The central claim rests on standard numerical-analysis assumptions about boundary smoothness and uniform point sampling; no new physical entities or heavily fitted parameters are introduced in the abstract.

axioms (2)
  • domain assumption The moving boundary remains sufficiently smooth for the GMLS local charts to achieve the stated high-order approximation of curvature and geometric quantities.
    Invoked to derive the error bound in terms of point-cloud size and boundary smoothness.
  • domain assumption The point cloud remains uniformly sampled on the boundary throughout the evolution.
    Required for the stated dependence of the error bound on point-cloud size.

pith-pipeline@v0.9.0 · 5732 in / 1473 out tokens · 67947 ms · 2026-05-21T20:23:37.212098+00:00 · methodology

discussion (0)

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

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

  • IndisputableMonolith/Foundation/AbsoluteFloorClosure.lean reality_from_one_distinction unclear
    ?
    unclear

    Relation between the paper passage and the cited Recognition theorem.

    We provide a rigorous convergence analysis for the proposed spatial discretization, establishing consistency and stability under certain conditions. The resulting error bound is derived in terms of the size of the uniformly sampled point cloud data on the moving boundary, the smoothness of the boundary, and the order of the numerical quadrature rule.

What do these tags mean?
matches
The paper's claim is directly supported by a theorem in the formal canon.
supports
The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends
The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses
The paper appears to rely on the theorem as machinery.
contradicts
The paper's claim conflicts with a theorem or certificate in the canon.
unclear
Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.

Reference graph

Works this paper leans on

37 extracted references · 37 canonical work pages

  1. [1]

    Routledge, 2018

    Vasilios Alexiades.Mathematical modeling of melting and freezing processes. Routledge, 2018

    work page 2018

  2. [2]

    Chen and A

    X. Chen and A. Friedman. A free boundary problem for an elliptic-hyperbolic system: an application to tumor growth.SIAM Journal on Mathematical Analysis, 35(4):974–986, 2003

    work page 2003

  3. [3]

    Diffusion maps.Applied and computational harmonic analysis, 21(1):5–30, 2006

    Ronald R Coifman and Stéphane Lafon. Diffusion maps.Applied and computational harmonic analysis, 21(1):5–30, 2006

    work page 2006

  4. [4]

    Constantin and M Pugh

    P . Constantin and M Pugh. Global solutions for small data to the hele-shaw problem.Nonlinearity, 6(3):393, 1993

    work page 1993

  5. [5]

    Nonlinear simulation of tumor growth.Journal of mathemati- cal biology, 46(3):191–224, 2003

    Vittorio Cristini, John Lowengrub, and Qing Nie. Nonlinear simulation of tumor growth.Journal of mathemati- cal biology, 46(3):191–224, 2003

    work page 2003

  6. [6]

    DiBenedetto and A

    E. DiBenedetto and A. Friedman. The ill-posed hele-shaw model and the stefan problem for supercooled water. Transactions of the American Mathematical Society, 282(1):183–204, 1984

    work page 1984

  7. [7]

    Friedman

    A. Friedman. Free boundary problems in biology.Phil. Trans. R. Soc. A, 373(2050):20140368, 2015

    work page 2050

  8. [8]

    Friedman and W

    A. Friedman and W . Hao. A mathematical model of atherosclerosis with reverse cholesterol transport and associated risk factors.Bulletin of mathematical biology, 77(5):758–781, 2015

    work page 2015

  9. [9]

    Friedman and J

    A. Friedman and J. Spruck.Variational and free boundary problems, volume 53. Springer Science & Business Media, 2012

    work page 2012

  10. [10]

    Prentice Hall, 1963

    Avner Friedman.Generalized functions and partial differential equations. Prentice Hall, 1963

    work page 1963

  11. [11]

    Meshfree methods on manifolds for hydro- dynamic flows on curved surfaces: A Generalized Moving Least-Squares (GMLS) approach.J

    Ben J Gross, Nathaniel Trask, Paul Kuberry, and Paul J Atzberger. Meshfree methods on manifolds for hydro- dynamic flows on curved surfaces: A Generalized Moving Least-Squares (GMLS) approach.J. Comput. Phys., 409:109340, 2020

    work page 2020

  12. [12]

    Elsevier, 2017

    Sushil Chandra Gupta.The classical Stefan problem: basic concepts, modelling and analysis with quasi- analytical solutions and methods, volume 45. Elsevier, 2017

    work page 2017

  13. [13]

    W . Hao, E. Crouser, and A. Friedman. Mathematical model of sarcoidosis.Proceedings of the National Academy of Sciences, 111(45):16065–16070, 2014

    work page 2014

  14. [14]

    Hao and A

    W . Hao and A. Friedman. The LDL-HDL profile determines the risk of atherosclerosis: a mathematical model. PLoS ONE, 9(3):e90497, 2014

    work page 2014

  15. [15]

    W . Hao, J. Hauenstein, B. Hu, Y. Liu, A. Sommese, and Y.-T . Zhang. Bifurcation for a free boundary problem modeling the growth of a tumor with a necrotic core.Nonlinear Analysis: Real World Applications, 13(2):694– 709, 2012. 21 APREPRINT- OCTOBER17, 2025

    work page 2012

  16. [16]

    W . Hao, J. Hauenstein, B. Hu, T . McCoy, and A. Sommese. Computing steady-state solutions for a free boundary problem modeling tumor growth by stokes equation.Journal of Computational and Applied Mathematics, 237(1):326–334, 2013

    work page 2013

  17. [17]

    W . Hao, J. Hauenstein, B. Hu, and A. Sommese. A three-dimensional steady-state tumor system.Applied Mathematics and Computation, 218(6):2661–2669, 2011

    work page 2011

  18. [18]

    W . Hao, B. Hu, and A. Sommese. Numerical algebraic geometry and differential equations. InFuture Vision and Trends on Shapes, Geometry and Algebra, pages 39–53. Springer, 2014

    work page 2014

  19. [19]

    A homotopy method for parameter estimation of nonlinear differential equations with multiple optima.Journal of Scientific Computing, 74(3):1314–1324, 2018

    Wenrui Hao. A homotopy method for parameter estimation of nonlinear differential equations with multiple optima.Journal of Scientific Computing, 74(3):1314–1324, 2018

    work page 2018

  20. [20]

    Radial basis approximation of tensor fields on manifolds: from operator estimation to manifold learning.Journal of Machine Learning Research, 24(345):1–85, 2023

    John Harlim, Shixiao Willing Jiang, and John Wilson Peoples. Radial basis approximation of tensor fields on manifolds: from operator estimation to manifold learning.Journal of Machine Learning Research, 24(345):1–85, 2023

    work page 2023

  21. [21]

    Hele-Shaw

    H. Hele-Shaw. Flow of water.Nature, 58:520, 1898

    work page

  22. [22]

    A hybrid method for moving interface problems with application to the hele–shaw flow.Journal of Computational Physics, 134(2):236–252, 1997

    Thomas Y Hou, Zhilin Li, Stanley Osher, and Hongkai Zhao. A hybrid method for moving interface problems with application to the hele–shaw flow.Journal of Computational Physics, 134(2):236–252, 1997

    work page 1997

  23. [23]

    Ghost point diffusion maps for solving elliptic PDEs on manifolds with classical boundary conditions.Communications on Pure and Applied Mathematics, 76(2):337–405, 2023

    Shixiao Willing Jiang and John Harlim. Ghost point diffusion maps for solving elliptic PDEs on manifolds with classical boundary conditions.Communications on Pure and Applied Mathematics, 76(2):337–405, 2023

    work page 2023

  24. [24]

    Generalized finite difference method on unknown manifolds.Journal of Computational Physics, 502:112812, 2024

    Shixiao Willing Jiang, Rongji Li, Qile Yan, and John Harlim. Generalized finite difference method on unknown manifolds.Journal of Computational Physics, 502:112812, 2024

    work page 2024

  25. [25]

    Springer, 1989

    Rainer Kress.Linear integral equations, volume 82. Springer, 1989

    work page 1989

  26. [26]

    Springer, 1999

    Rainer Kress.Linear integral equations, volume 82. Springer, 1999

    work page 1999

  27. [27]

    Eigenvalues of tridiagonal pseudo-toeplitz matrices

    Devadatta Kulkarni, Darrell Schmidt, and Sze-Kai Tsui. Eigenvalues of tridiagonal pseudo-toeplitz matrices. Linear Algebra and its Applications, 297:63–80, 1999

    work page 1999

  28. [28]

    A rescaling scheme with application to the long-time simulation of viscous fingering in a hele–shaw cell.Journal of Computational Physics, 225(1):554–567, 2007

    Shuwang Li, John S Lowengrub, and Perry H Leo. A rescaling scheme with application to the long-time simulation of viscous fingering in a hele–shaw cell.Journal of Computational Physics, 225(1):554–567, 2007

    work page 2007

  29. [29]

    Solving partial differential equations on point clouds.SIAM J

    Jian Liang and Hongkai Zhao. Solving partial differential equations on point clouds.SIAM J. Sci. Comput., 35(3):A1461–A1486, 2013

    work page 2013

  30. [30]

    Bifurcation analysis of a free boundary model of vascular tumor growth with a necrotic core and chemotaxis.Journal of mathematical biology, 86(1):19, 2023

    Min-Jhe Lu, Wenrui Hao, Bei Hu, and Shuwang Li. Bifurcation analysis of a free boundary model of vascular tumor growth with a necrotic core and chemotaxis.Journal of mathematical biology, 86(1):19, 2023

    work page 2023

  31. [31]

    An improved geometry-aware curvature discretization for level set methods: application to tumor growth.journal of Computational Physics, 215(2):392–401, 2006

    Paul Macklin and John Lowengrub. An improved geometry-aware curvature discretization for level set methods: application to tumor growth.journal of Computational Physics, 215(2):392–401, 2006

    work page 2006

  32. [32]

    Variation properties of some simple boundary integral equations

    W McLean. Variation properties of some simple boundary integral equations. InTheoretical and Numer- ical Aspects of Geometric Variational Problems, volume 26, pages 168–179. Australian National University, Mathematical Sciences Institute, 1991

    work page 1991

  33. [33]

    On generalized moving least squares and diffuse derivatives.IMA Journal of Numerical Analysis, 32(3):983–1000, 2012

    Davoud Mirzaei, Robert Schaback, and Mehdi Dehghan. On generalized moving least squares and diffuse derivatives.IMA Journal of Numerical Analysis, 32(3):983–1000, 2012

    work page 2012

  34. [34]

    A review of physics of moving contact line dynamics models and its applications in interfacial science.Journal of Applied Physics, 132(8), 2022

    Alireza Mohammad Karim. A review of physics of moving contact line dynamics models and its applications in interfacial science.Journal of Applied Physics, 132(8), 2022

    work page 2022

  35. [35]

    Saffman and G

    P . Saffman and G. Taylor. The penetration of a fluid into a porous medium or hele-shaw cell containing a more viscous liquid. InProceedings of the Royal Society of London A: Mathematical, Physical and Engineering Sciences, volume 245, pages 312–329. The Royal Society, 1958

    work page 1958

  36. [36]

    Princeton university press, 1996

    Steven Vogel.Life in moving fluids: the physical biology of flow. Princeton university press, 1996

    work page 1996

  37. [37]

    Spectral methods for solving elliptic pdes on unknown manifolds.Journal of Computational Physics, 486:112132, 2023

    Qile Yan, Shixiao Willing Jiang, and John Harlim. Spectral methods for solving elliptic pdes on unknown manifolds.Journal of Computational Physics, 486:112132, 2023. 22

    work page 2023