An Omni-Temporal Theory for Hydrodynamic Dispersion and Reaction in Porous Media

Kyle C. Smith; Md Abdul Hamid

arxiv: 2505.06063 · v2 · pith:23HYQETHnew · submitted 2025-05-09 · ⚛️ physics.flu-dyn · math-ph· math.MP

An Omni-Temporal Theory for Hydrodynamic Dispersion and Reaction in Porous Media

Md Abdul Hamid , Kyle C. Smith This is my paper

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

classification ⚛️ physics.flu-dyn math-phmath.MP

keywords hydrodynamic dispersionporous mediafrequency domainsolute transportupscalingadvection-diffusionreactive transportbreakthrough curves

0 comments

The pith

A frequency-based theory captures both fast and slow dispersion in porous media through volume averaging of Fourier-transformed pore-scale equations.

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

The paper sets out an omni-temporal framework that avoids the long-time limits of classical dispersion theory. Instead it works in the frequency domain to keep track of both rapid early transients and slower later spreading of solutes. Volume averaging is applied to the Fourier-transformed advection-diffusion-reaction equation inside a periodic unit cell. This step produces four frequency-dependent coefficients that function as transfer functions. These can be inverted via Fourier transform to recover time-domain predictions such as breakthrough curves that match direct simulations even for fast pulses.

Core claim

Volume averaging of the Fourier-transformed pore-scale advection-diffusion equation over a periodic unit cell produces four frequency-dependent upscaled transport coefficients: a dispersion tensor, an advection-suppression transfer function, a spectral Sherwood number, and a reactivity-bias vector. These coefficients serve as transfer functions that relate microscopic driving forces to effective fluxes in the frequency domain, so that inverse Fourier transformation yields predictions of transient transport dynamics.

What carries the argument

Volume averaging in the frequency domain of the pore-scale advection-diffusion-reaction equation over a periodic unit cell, which generates four frequency-dependent transfer functions for upscaled transport.

If this is right

Analytical expressions for the four transfer functions become available for Poiseuille flow in parallel plates and circular tubes.
The resulting transfer functions can be inserted into an FFT framework to generate breakthrough curves for both reactive and non-reactive cases.
For fast solute pulses the theory matches direct numerical simulations while conventional long-time theories overpredict propagation speeds by orders of magnitude.
The same coefficients apply to cross-flow through periodic arrays of square rods, showing the method works for both active and inactive solid phases.

Where Pith is reading between the lines

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

The transfer functions could be precomputed once for a given unit-cell geometry and then reused for many different inlet conditions or reaction rates.
Numerical evaluation of the averaging step for random rather than periodic media would turn the theory into a practical upscaling tool for field-scale models.
Coupling the frequency-domain coefficients with existing reservoir simulators might improve early-time forecasts of contaminant plumes or drug delivery.
Testing the inverted curves against laboratory tracer experiments in packed beds would check whether the periodic-cell assumption holds under real flow conditions.

Load-bearing premise

Volume averaging the Fourier-transformed pore-scale advection-diffusion equation over a periodic unit cell produces accurate effective transfer functions whose inverse Fourier transforms recover correct time-domain behavior.

What would settle it

A direct numerical simulation of breakthrough curves for a fast solute pulse between inactive parallel plates; if the theory's FFT-inverted curves deviate markedly from the simulated early arrival times while asymptotic theory also fails, the averaging step would be falsified.

read the original abstract

A frequency-based omni-temporal dispersion theory is developed to capture the transient interplay between diffusion, advection, and reaction during solute transport through porous media. Unlike classical asymptotic dispersion theories, which commonly rely on long-time approximation, the proposed framework simultaneously captures both fast and slow components of dispersion. The theory is formulated by volume averaging the Fourier-transformed pore-scale advection-diffusion equation, yielding four frequency-dependent upscaled transport coefficients for a periodic unit cell: a dispersion tensor, an advection-suppression transfer function, a spectral Sherwood number, and a reactivity-bias vector. These coefficients act as transfer functions that relate microscopic driving forces to corresponding effective fluxes in the frequency domain, enabling prediction of transient transport dynamics in the time domain through inverse Fourier transformation. The utility of the proposed framework is demonstrated by deriving analytical expressions for the transfer functions in Poiseuille flow between parallel plates and through circular tubes, and subsequently using them within a Fast Fourier Transform framework to obtain breakthrough curves. For fast solute pulses between inactive parallel plates, the proposed theory produces breakthrough curves in close agreement with direct numerical simulations, whereas conventional asymptotic theory overpredicts propagation rates by orders of magnitude. Finally, the framework is applied to reactive and non-reactive porous media consisting of periodic arrays of square rods under cross flow, demonstrating the generality and versatility of the proposed omni-temporal theory.

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.

Desk Editor's Note private letter to a colleague

This paper gives a frequency-domain volume-averaging route to transient dispersion and reaction coefficients that handles fast pulses better than standard long-time theories in the cases they check.

read the letter

The main thing to know is that the authors volume-average the Fourier-transformed pore-scale advection-diffusion-reaction equation over a periodic cell and extract four frequency-dependent transfer functions: a dispersion tensor, an advection-suppression term, a spectral Sherwood number, and a reactivity-bias vector. These are then inverted via FFT to produce time-domain breakthrough curves that match direct simulations for fast solute pulses between inactive parallel plates, where the usual asymptotic theory overpredicts speeds by a lot. They also derive closed analytical forms for Poiseuille flow and circular tubes and show results for reactive square-rod arrays under cross flow. That agreement on the simple inactive case is a concrete check, and starting from the pore-scale equation without fitted parameters keeps the approach grounded. The extension beyond long-time limits is the clearest novelty here. The soft spot is the limited validation. The abstract reports agreement for one fast-pulse scenario but gives no error bars, no sweep across frequencies or reaction strengths, and no detail on how the deviation concentration problem is closed in frequency space. If that closure leaves residual non-local memory at some frequencies, the inverse transform could lose accuracy even if the averaged equation looks closed on paper. The stress-test note about unstated approximations for the deviation fields lines up with what is visible from the abstract and the reported checks. This is aimed at readers in hydrology or chemical engineering who need transient upscaled models for porous-media transport. Someone already working on volume averaging or frequency-domain methods would see the practical payoff in the transfer-function route. It is worth sending to peer review because the core construction is new, the simple-geometry test is encouraging, and the gaps are fixable with more comparisons rather than fatal to the idea.

Referee Report

2 major / 3 minor

Summary. The manuscript develops a frequency-based omni-temporal theory for hydrodynamic dispersion and reaction in porous media. It applies volume averaging to the Fourier-transformed pore-scale advection-diffusion-reaction equation over a periodic unit cell, producing four frequency-dependent upscaled coefficients: a dispersion tensor, an advection-suppression transfer function, a spectral Sherwood number, and a reactivity-bias vector. These act as transfer functions relating microscopic forces to effective fluxes in frequency space, which are inverted via FFT to recover time-domain breakthrough curves capturing both fast and slow transients. Analytical expressions are derived for Poiseuille flow between parallel plates and in circular tubes; the framework is then applied to non-reactive and reactive periodic arrays of square rods under cross flow, with reported agreement to direct numerical simulations for a fast-pulse case where asymptotic theory fails.

Significance. If the central claims hold after addressing the noted issues, the work would represent a meaningful advance by providing a unified approach to transient dispersion and reaction that avoids long-time asymptotic restrictions. The derivation of explicit analytical transfer functions for canonical geometries and their use in FFT-based predictions are clear strengths, as is the demonstration of improved performance over classical methods in at least one fast-transient regime. This could have utility in environmental and chemical engineering applications involving reactive transport in porous media.

major comments (2)

[Theory derivation] The volume-averaging closure for the deviation concentration field in frequency space (theory section following the Fourier transform of the pore-scale ADE) is presented as directly furnishing the four exact frequency-dependent coefficients without residual memory or non-local effects. This closure is load-bearing for the invertibility claim via FFT, yet the manuscript provides no explicit cell-problem solution details, bounds on validity across frequencies, or tests confirming absence of non-invertible artifacts, particularly for reactive cases.
[Numerical results and comparisons] Validation is restricted to a single fast-pulse non-reactive case between parallel plates (results section) with no error bars, quantitative error metrics (e.g., L2 norms or breakthrough time deviations), or comparisons across slow-transient, high-Pe, or reactive regimes. This leaves the claim of close agreement with DNS and superiority over asymptotic theory only partially supported.

minor comments (3)

[Abstract and results] The abstract and results text describe 'close agreement' with DNS; adding specific quantitative measures (e.g., relative L2 error or peak arrival time differences) would improve clarity.
[Figures] Figure captions for breakthrough curves should explicitly label all three curves (proposed theory, asymptotic theory, DNS) and note the Fourier inversion parameters used.
[Theory] Notation for the four coefficients could be introduced with a summary table early in the theory section to aid readability.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for their constructive and detailed review of our manuscript. We have carefully addressed each major comment below and will incorporate revisions to enhance the clarity of the theoretical derivation and the robustness of the numerical validation.

read point-by-point responses

Referee: [Theory derivation] The volume-averaging closure for the deviation concentration field in frequency space (theory section following the Fourier transform of the pore-scale ADE) is presented as directly furnishing the four exact frequency-dependent coefficients without residual memory or non-local effects. This closure is load-bearing for the invertibility claim via FFT, yet the manuscript provides no explicit cell-problem solution details, bounds on validity across frequencies, or tests confirming absence of non-invertible artifacts, particularly for reactive cases.

Authors: We agree that the manuscript would benefit from greater explicitness in the closure procedure. In the revised version, we will expand the relevant theory section to present the full cell problems for the deviation concentration and the resulting analytical expressions for the four frequency-dependent coefficients. We will also add a discussion of validity bounds, noting that the derivation holds for all frequencies under the assumptions of linear reactions and periodic media. To confirm the absence of non-invertible artifacts, we will include a supplementary numerical test of FFT inversion stability for a reactive case. revision: yes
Referee: [Numerical results and comparisons] Validation is restricted to a single fast-pulse non-reactive case between parallel plates (results section) with no error bars, quantitative error metrics (e.g., L2 norms or breakthrough time deviations), or comparisons across slow-transient, high-Pe, or reactive regimes. This leaves the claim of close agreement with DNS and superiority over asymptotic theory only partially supported.

Authors: The present manuscript emphasizes the fast-pulse regime to illustrate the theory's advantage where classical asymptotics break down. We acknowledge that broader validation strengthens the claims. In the revision, we will add quantitative metrics such as L2 norms and breakthrough-time deviations, include error bars from repeated simulations, and provide additional comparisons for slow-transient, high-Pe, and reactive regimes using the square-rod arrays. These additions will more fully substantiate agreement with DNS and superiority over asymptotic theory. revision: yes

Circularity Check

0 steps flagged

Derivation from pore-scale ADE via Fourier transform and volume averaging is self-contained with no reduction to inputs

full rationale

The paper starts from the standard pore-scale advection-diffusion-reaction equation, applies a Fourier transform in time, performs volume averaging over a periodic unit cell, and solves the resulting cell problem to obtain four frequency-dependent coefficients (dispersion tensor, advection-suppression transfer function, spectral Sherwood number, reactivity-bias vector). These are then used as transfer functions whose inverse FFT yields time-domain predictions. This chain is a direct homogenization procedure; the effective coefficients are outputs of the cell problem rather than fitted parameters or self-referential definitions. No load-bearing self-citations, uniqueness theorems, or ansatzes imported from prior author work are indicated in the abstract or described steps. Analytical solutions for Poiseuille flow and square-rod arrays are derived explicitly and validated against DNS, confirming the derivation does not collapse to its inputs by construction.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The central claim rests on standard mathematical operations and one domain assumption; no free parameters or invented entities are introduced in the abstract.

axioms (1)

domain assumption Volume averaging of the Fourier-transformed pore-scale advection-diffusion equation is valid for periodic unit cells.
This step is the explicit starting point for deriving the four frequency-dependent coefficients.

pith-pipeline@v0.9.0 · 5770 in / 1231 out tokens · 72283 ms · 2026-05-22T16:01:23.897498+00:00 · methodology

discussion (0)

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/AlphaCoordinateFixation.lean costAlphaLog_high_calibrated_iff unclear

?

unclear
Relation between the paper passage and the cited Recognition theorem.

ansatz … c̃ = a + b·∇⟨c⟩ … closure problems for a and b

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.