pith. sign in

arxiv: 2605.11788 · v2 · pith:7XQS7H74new · submitted 2026-05-12 · 🧮 math.AP · math-ph· math.MP

The unified transform for Burgers' equation: Application to unsaturated flow in finite interval

Konstantinos Kalimeris , Leonidas Mindrinos , Athanasios Paraskevopoulos This is my paper

Pith reviewed 2026-05-13 05:47 UTC · model grok-4.3

classification 🧮 math.AP math-phmath.MP
keywords Burgers equationunified transform methodunsaturated flowRichards equationHopf-Cole transformationfinite intervalintegral representationdiffusion equation
0
0 comments X

The pith

The Unified Transform Method supplies an explicit integral representation for solutions of the linearized Burgers equation on a finite interval.

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

The paper reduces Richards' equation for unsaturated flow to Burgers' equation under constant diffusivity and quadratic conductivity assumptions, then applies the Hopf-Cole transformation to obtain a diffusion equation with mixed boundary conditions. It solves this initial-boundary value problem using the Unified Transform Method to derive an explicit integral form for the solution. Numerical evaluation of this form reproduces the results of classical Fourier series solutions exactly while exhibiting improved convergence and stability. The approach is demonstrated on two hydrological application examples involving one-dimensional vertical infiltration.

Core claim

The Unified Transform Method yields an explicit integral representation of the solution to the diffusion equation on a finite interval with mixed boundary conditions that arises from linearizing Burgers' equation for modeling unsaturated flow.

What carries the argument

The Unified Transform Method, which derives a global relation from the PDE and boundary conditions to construct an explicit integral representation of the solution.

If this is right

  • The method allows direct computation of water content profiles in finite soil columns without relying on series expansions.
  • Improved numerical stability facilitates more reliable simulations of infiltration processes over longer times or with sharper gradients.
  • Applications in hydrology can leverage this for better handling of boundary conditions in bounded domains.

Where Pith is reading between the lines

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

  • Similar integral representations could apply to other nonlinear PDEs in hydrology after suitable transformations.
  • Testing the method on problems with variable diffusivity would reveal its robustness beyond the constant case assumed here.
  • Integration with existing hydrological software might improve efficiency in predicting soil moisture dynamics.

Load-bearing premise

The assumptions that diffusivity is constant and hydraulic conductivity depends quadratically on water content, which permit the exact reduction of Richards' equation to Burgers' equation and its linearization.

What would settle it

Numerical tests on the finite interval problem where the integral representation diverges or fails to match the Fourier series solution for chosen initial and boundary data.

Figures

Figures reproduced from arXiv: 2605.11788 by Athanasios Paraskevopoulos, Konstantinos Kalimeris, Leonidas Mindrinos.

Figure 1
Figure 1. Figure 1: The boundary ∂D+ (black curve) of the domain D+ and the deformed contour C + (blue curve) given by (20). Let V1(λ, x, t) denote the integrand of w1. Then, the term w1 is expressed as: Z C+ V1(λ, x, t)dλ = Z −ℓ −∞ V1(z1(s), x, t)z ′ 1 (s)ds + Z ℓ −ℓ V1(z2(s), x, t)z ′ 2 (s)ds + Z ∞ ℓ V1(z3(s), x, t)z ′ 3 (s)ds = Z ∞ ℓ  V1(z3(s), x, t)e i π 6 − V1(z1(−s), x, t)e i 5π 6  ds + Z ℓ −ℓ V1(z2(s), x, t)ds . 4 Nu… view at source ↗
Figure 2
Figure 2. Figure 2: The two different contour deformations C + considered in the first example. Note that the series had to be truncated at N = 2000 to obtain an illustration of the solution with a sufficiently small deviation from the actual solution. A visualization of the truncation error for different numbers of Fourier nodes is shown in [PITH_FULL_IMAGE:figures/full_fig_p008_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: The water content θ (6) using the analytical solution (14) (left) and the Fourier series solution, see [16], (right) for the setup of the first example. We tested our solution on two numerical experiments modeling real-world applications of sub￾surface water flow under different ponding conditions. By exploiting tools from complex analysis, we derived efficient numerical schemes that are easy to implement.… view at source ↗
Figure 4
Figure 4. Figure 4: Convergence of the series solution to the analytical one (blue curve) for [PITH_FULL_IMAGE:figures/full_fig_p010_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: Example 1: The three-dimensional visualization of the analytical solution [PITH_FULL_IMAGE:figures/full_fig_p010_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: The water content θ (6) using the analytical solution (14) (left) and the Fourier series solution, see [16], (right) for the setup of the second example. [9] F. P. J. De Barros, M. J. Colbrook, and A. S. Fokas. A hybrid analytical-numerical method for solving advection-dispersion problems on a half-line. International Journal of Heat and Mass Transfer, 139:482–491, 2019. [10] A. Fokas and S. De Lillo. The … view at source ↗
Figure 7
Figure 7. Figure 7: Example 2: The three-dimensional visualization of the analytical solution [PITH_FULL_IMAGE:figures/full_fig_p012_7.png] view at source ↗
read the original abstract

In this paper, we focus on one-dimensional vertical infiltration, assuming constant diffusivity and a quadratic relationship between hydraulic conductivity and water content. Under these assumptions, Richards' equation reduces to Burgers' equation, which we then linearize via the Hopf-Cole transformation. This turns the initial boundary value problem into a diffusion equation on a finite interval with mixed boundary conditions. To solve it, we use the Unified Transform Method (also known as the Fokas method). This approach gives an explicit integral representation of the solution, and when evaluated numerically, the results match classical Fourier series solutions exactly, but with better convergence and stability. Two examples from hydrological applications are examined.

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 applies the Unified Transform Method (Fokas method) to the heat equation obtained after Hopf-Cole linearization of Burgers' equation, which arises from Richards' equation for one-dimensional vertical infiltration under the assumptions of constant diffusivity and quadratic hydraulic conductivity. It derives an explicit integral representation for the solution on a finite interval subject to mixed boundary conditions and presents numerical evaluations that are asserted to match classical Fourier-series solutions exactly while showing improved convergence and stability. Two hydrological application examples are examined.

Significance. If the derivation and numerical claims hold, the work supplies an alternative explicit representation for linearized infiltration problems that may avoid Gibbs phenomena or slow decay associated with eigenfunction expansions, offering potential advantages for stable long-time simulations in hydrology. The approach aligns with known strengths of contour-integral methods for linear IBVPs and could extend to related nonlinear diffusion models once the linearization assumptions are satisfied.

major comments (2)
  1. [§3] §3, Eq. (3.8): the global relation is stated but the explicit elimination of the unknown boundary values (arising from the mixed Dirichlet-Neumann conditions after Hopf-Cole) is not carried out in detail; without this step the claimed explicit integral representation cannot be verified independently.
  2. [§4.2] §4.2, Figure 3: the reported numerical agreement with the Fourier solution is shown only via overlaid plots; no L²-error tables, convergence rates versus number of quadrature points, or contour-deformation parameters are supplied, so the asserted superiority in convergence and stability remains unquantified.
minor comments (2)
  1. [Abstract] The abstract and §1 state that the numerical results 'match exactly'; this should be rephrased as 'agree to machine precision' to reflect floating-point evaluation of the contour integrals.
  2. [§2] Notation for the transformed boundary functions (e.g., the functions appearing in the global relation) is introduced without a dedicated table or list of symbols, reducing 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 each major comment below and indicate the revisions planned to improve clarity and support for the claims.

read point-by-point responses

  1. Referee: [§3] §3, Eq. (3.8): the global relation is stated but the explicit elimination of the unknown boundary values (arising from the mixed Dirichlet-Neumann conditions after Hopf-Cole) is not carried out in detail; without this step the claimed explicit integral representation cannot be verified independently.

    Authors: We appreciate this observation. While the global relation and its role in the unified transform are presented in Section 3, we agree that the elimination of the unknown boundary values for the mixed conditions could be shown with additional explicit steps to facilitate independent verification. In the revised manuscript we will expand this derivation, inserting the intermediate algebraic manipulations that solve for the unknown transforms and yield the final contour-integral formula. revision: yes

  2. Referee: [§4.2] §4.2, Figure 3: the reported numerical agreement with the Fourier solution is shown only via overlaid plots; no L²-error tables, convergence rates versus number of quadrature points, or contour-deformation parameters are supplied, so the asserted superiority in convergence and stability remains unquantified.

    Authors: We agree that quantitative evidence is needed to substantiate the asserted advantages in convergence and stability. In the revised version we will add L²-error tables for representative test cases, report observed convergence rates as a function of quadrature points, and document the contour-deformation parameters used in the numerical evaluations of the integral representation. revision: yes

Circularity Check

0 steps flagged

No significant circularity identified

full rationale

The paper reduces Richards' equation to Burgers' equation under constant diffusivity and quadratic conductivity assumptions, applies the Hopf-Cole transformation to obtain the linear heat equation on a finite interval with mixed boundary conditions, and then invokes the established Unified Transform Method (Fokas method) to derive an explicit contour-integral representation. This representation is shown to agree numerically with the classical Fourier-series solution of the identical IBVP, which is expected once both methods correctly enforce the global relation and eliminate unknown boundary values. No step reduces by definition to its own inputs, no fitted parameter is relabeled as a prediction, and no load-bearing premise rests on a self-citation chain; the derivation is self-contained against the external benchmark of the standard eigenfunction expansion.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The central claim rests on two standard elements: the Hopf-Cole transformation (standard math) and the domain-specific reduction of Richards' equation under the stated constitutive assumptions (domain assumption). No free parameters or invented entities are introduced.

axioms (2)
  • standard math Hopf-Cole transformation converts Burgers' equation into a linear diffusion equation
    Standard change of variables routinely used for this nonlinear PDE.
  • domain assumption Constant diffusivity and quadratic hydraulic conductivity-water content relation reduce Richards' equation to Burgers' equation
    Specific modeling choice for the infiltration problem stated in the abstract.

pith-pipeline@v0.9.0 · 5415 in / 1261 out tokens · 53541 ms · 2026-05-13T05:47:02.010168+00:00 · methodology

discussion (0)

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