pith. sign in

arxiv: 2410.05176 · v2 · submitted 2024-10-07 · 🧮 math.AP

Homogenized Equations for Isentropic Gas in a Pipe with Periodically-Varying Cross-Section

Laila S. Busaleh , David I. Ketcheson This is my paper

Pith reviewed 2026-05-23 19:58 UTC · model grok-4.3

classification 🧮 math.AP
keywords homogenized equationsmultiple-scale perturbationisentropic gasperiodic cross-sectionsolitary wavesdispersive termsRiemann solverpipe flow
0
0 comments X

The pith

Multiple-scale perturbation yields homogenized equations with dispersion for isentropic gas in periodically varying pipes, leading to solitary waves rather than shocks.

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

The paper applies multiple-scale perturbation theory to derive effective equations for isentropic gas flowing in a narrow pipe whose cross-section varies periodically. These homogenized equations form a constant-coefficient system that includes higher-order dispersive derivative terms. An approximate Riemann solver is developed for the original system to enable numerical comparisons. The solutions of the homogenized system consistently appear as solitary waves instead of the shock waves typical in standard gas dynamics.

Core claim

Using multiple-scale perturbation theory, the authors derive homogenized effective equations for an isentropic gas in a pipe with periodically-varying cross-section. These take the form of a constant-coefficient system of evolution equations including dispersive higher-order derivative terms. Numerical comparisons using an approximate Riemann solver for the variable-cross-section equations show that the resulting solutions take the form of solitary waves rather than shock waves under fairly general conditions.

What carries the argument

Multiple-scale perturbation theory that produces a constant-coefficient homogenized system with dispersive terms from the variable-geometry isentropic gas equations.

Load-bearing premise

The derivation requires clear scale separation between the periodic pipe variation and the flow features of interest, plus the validity of the isentropic assumption for the gas.

What would settle it

A direct numerical simulation of the original equations with periodic cross-section where the period is comparable to the wavelength, showing persistent shocks instead of solitary waves, would falsify the homogenized model's applicability.

Figures

Figures reproduced from arXiv: 2410.05176 by David I. Ketcheson, Laila S. Busaleh.

Figure 1
Figure 1. Figure 1: Flow in a pipe of periodically-varying cross-section. [PITH_FULL_IMAGE:figures/full_fig_p003_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: Comparison of the homogenized system (dashed line) and the variable-coefficient system (solid [PITH_FULL_IMAGE:figures/full_fig_p012_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: Comparison of the homogenized system (dashed line) and the variable coefficients system (solid [PITH_FULL_IMAGE:figures/full_fig_p013_3.png] view at source ↗
read the original abstract

We analyze the behavior of an isentropic gas in a narrow pipe with periodically-varying cross-sectional area. Using multiple-scale perturbation theory, we derive homogenized effective equations, which take the form of a constant-coefficient system of evolution equations, including dispersive higher-order derivative terms. We provide an approximate Riemann solver for the variable-cross-section isentropic gas equations, and compare numerical solutions of the original system with those of the homogenized system. We observe that the resulting solutions take the form of solitary waves, rather than shock waves, under fairly general conditions.

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 applies multiple-scale perturbation theory to derive homogenized constant-coefficient evolution equations, including dispersive higher-order terms, for isentropic gas flow in a narrow pipe with periodically varying cross-section. It also provides an approximate Riemann solver for the original variable-cross-section equations and compares numerical solutions of both systems, observing that the homogenized model produces solitary waves rather than shocks under fairly general conditions.

Significance. If the multiple-scale derivation is valid and the numerical observations hold, the work provides an effective dispersive model that explains the regularization of shocks into solitary waves due to the periodic geometry. This could have implications for modeling wave propagation in variable-geometry conduits. The numerical comparisons add practical value, but the lack of rigorous error analysis reduces the overall significance for mathematical analysis of the approximation.

major comments (2)
  1. [§2] §2 (Multiple-scale analysis): The derivation of the homogenized system with dispersive terms relies on the multiple-scale ansatz, but no a priori error estimate or rigorous justification is provided for the approximation when the scale separation parameter is not asymptotically small.
  2. [§4] §4 (Numerical comparisons): The numerical tests use small values of the period length ε and small amplitudes; this does not address whether the original system can form shocks when scale separation fails or amplitudes are O(1), undermining the claim that solitary waves occur under 'fairly general conditions'.
minor comments (2)
  1. [Abstract] The abstract mentions 'fairly general conditions' without specifying what they are; this should be clarified in the introduction or conclusion.
  2. [Notation] The definition of the cross-section variation function A(x) should be stated explicitly with its periodicity assumption early in the paper.

Simulated Author's Rebuttal

2 responses · 1 unresolved

We thank the referee for the detailed report and constructive feedback. We address the two major comments point by point below. The derivation is formal perturbation theory, and we agree that the numerical regime is limited to small ε; we will revise statements accordingly while maintaining the paper's focus on the effective model and observations within its validity range.

read point-by-point responses

  1. Referee: §2 (Multiple-scale analysis): The derivation of the homogenized system with dispersive terms relies on the multiple-scale ansatz, but no a priori error estimate or rigorous justification is provided for the approximation when the scale separation parameter is not asymptotically small.

    Authors: We agree that the multiple-scale derivation is formal and lacks rigorous a priori error estimates or justification for finite ε. This is standard for such asymptotic analyses in the literature on homogenization of hyperbolic systems. The manuscript emphasizes the formal derivation of the effective dispersive system and its numerical comparison rather than a rigorous approximation theorem. We will add an explicit remark in §2 clarifying the formal nature of the expansion and that it is expected to hold asymptotically as ε → 0. revision: partial

  2. Referee: §4 (Numerical comparisons): The numerical tests use small values of the period length ε and small amplitudes; this does not address whether the original system can form shocks when scale separation fails or amplitudes are O(1), undermining the claim that solitary waves occur under 'fairly general conditions'.

    Authors: The numerical tests are performed in the small-ε regime consistent with the homogenization ansatz, using amplitudes that remain within the perturbative setting. We acknowledge that the phrasing 'under fairly general conditions' in the abstract and introduction is imprecise and could be read as applying outside the small-ε, small-amplitude regime. We will revise the abstract, introduction, and conclusion to state that solitary-wave behavior is observed in the regime of small period length and moderate amplitudes where the homogenized model applies. No additional large-amplitude tests are planned for this revision, as they fall outside the paper's scope. revision: yes

standing simulated objections not resolved
  • Providing a rigorous a priori error estimate or full justification of the multiple-scale approximation for non-asymptotically small ε, which would require a separate analytical effort beyond the formal derivation and numerical validation presented.

Circularity Check

0 steps flagged

Multiple-scale derivation is independent; no fitted inputs or self-referential steps

full rationale

The paper applies standard multiple-scale perturbation theory directly to the original variable-cross-section isentropic gas system to obtain the homogenized constant-coefficient equations with dispersive terms. No coefficients are obtained by fitting to solutions of the target homogenized system, no self-citations supply load-bearing uniqueness theorems or ansatzes, and the numerical comparisons are presented as post-derivation validation rather than inputs. The observation of solitary waves follows from solving the derived system, which is constructed from the original PDEs without reduction to its own outputs. This is the most common honest non-finding for a first-principles asymptotic derivation.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The central claim rests on the validity of multiple-scale perturbation theory for the given scale separation and on the isentropic assumption; no free parameters, invented entities, or additional axioms are stated in the abstract.

axioms (2)
  • domain assumption Multiple-scale perturbation theory applies when the pipe-period length is small compared with the wavelength of interest.
    Invoked by the choice of homogenization method in the abstract.
  • domain assumption The gas obeys the isentropic closure.
    Stated in the title and abstract as the physical model under study.

pith-pipeline@v0.9.0 · 5622 in / 1411 out tokens · 19911 ms · 2026-05-23T19:58:36.459696+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

23 extracted references · 23 canonical work pages

  1. [1]

    Oller, Numerical Simulation of Mechanical Behavior of Composite Materials

    S. Oller, Numerical Simulation of Mechanical Behavior of Composite Materials. Springer Cham, first ed., 2014

    work page 2014

  2. [2]

    A brief introduction to homogenization and miscellaneous applications,

    G. Allaire, “A brief introduction to homogenization and miscellaneous applications,” in ESAIM: Pro- ceedings, vol. 37, pp. 1–49, EDP Sciences, September 2012

    work page 2012

  3. [3]

    Numerical simulation of cylindrical solitary waves in periodic media,

    M. Quezada de Luna and D. I. Ketcheson, “Numerical simulation of cylindrical solitary waves in periodic media,” Journal of Scientific Computing , pp. 672–689, 2013

    work page 2013

  4. [4]

    A dispersive effective medium for wave propagation in periodic composites,

    F. Santosa and W. Symes, “A dispersive effective medium for wave propagation in periodic composites,” SIAM Journal on Applied Mathematics , vol. 51, no. 4, p. 984–1005, 1991

    work page 1991

  5. [5]

    Two-dimensional wave propagation in layered periodic media,

    M. Quezada de Luna and D. I. Ketcheson, “Two-dimensional wave propagation in layered periodic media,” SIAM Journal on Applied Mathematics , vol. 74, no. 6, pp. 1852–1869, 2014

    work page 2014

  6. [6]

    One-dimensional wave dispersion in layered media,

    H. Ockendon, J. Ockendon, C. Farmer, and D. Allwright, “One-dimensional wave dispersion in layered media,” SIAM Journal on Applied Mathematics , vol. 75, no. 5, pp. 2128–2146, 2015

    work page 2015

  7. [7]

    Crime pays; homogenization for long times,

    G. Allaire, A. Lamacz, and J. Rauch, “Crime pays; homogenization for long times,” S´ eminaire Laurent Schwartz—EDP et applications , pp. 1–9, 2022

    work page 2022

  8. [8]

    Solitary waves in layered nonlinear media,

    R. J. LeVeque, Yong, and H. Darryl, “Solitary waves in layered nonlinear media,” SIAM Journal on Applied Mathematics, vol. 63, pp. 1539–1560, 2003

    work page 2003

  9. [9]

    Longitudinal strain solitary wave in a two-layered polymeric bar,

    G. Dreiden, K. Khusnutdinova, A. Samsonov, and I. Semenova, “Longitudinal strain solitary wave in a two-layered polymeric bar,” Strain, vol. 46, no. 6, pp. 589–598, 2010

    work page 2010

  10. [10]

    Coherent structures and carrier shocks in the nonlinear periodic maxwell equations,

    G. Simpson and M. I. Weinstein, “Coherent structures and carrier shocks in the nonlinear periodic maxwell equations,” Multiscale Modeling & Simulation , vol. 9, no. 3, pp. 955–990, 2011

    work page 2011

  11. [11]

    Shock dynamics in layered periodic media,

    D. Ketcheson and R. J. LeVeque, “Shock dynamics in layered periodic media,” Communications in Mathematical Sciences, vol. 10, no. 3, pp. 859–874, 2012

    work page 2012

  12. [12]

    Diffractons: Solitary waves created by diffraction in periodic media,

    D. I. Ketcheson and M. Q. de Luna, “Diffractons: Solitary waves created by diffraction in periodic media,” Multiscale Modeling and Simulation , vol. 13, no. 1, pp. 440–458, 2015

    work page 2015

  13. [13]

    Numerical modeling of nonlinear acoustic waves in a tube connected with helmholtz resonators,

    B. Lombard and J.-F. Mercier, “Numerical modeling of nonlinear acoustic waves in a tube connected with helmholtz resonators,” Journal of Computational Physics , vol. 259, pp. 421–443, 2014

    work page 2014

  14. [14]

    Dispersive and dispersive-like bores in channels with sloping banks,

    R. Chassagne, A. G. Filippini, M. Ricchiuto, and P. Bonneton, “Dispersive and dispersive-like bores in channels with sloping banks,” Journal of Fluid Mechanics , vol. 870, pp. 595–616, 2019

    work page 2019

  15. [15]

    Solitary water waves created by variations in bathymetry,

    M. Quezada de Luna and D. I. Ketcheson, “Solitary water waves created by variations in bathymetry,” Journal of Fluid Mechanics , vol. 917, p. A45, 2021

    work page 2021

  16. [16]

    Strongly nonlinear wave dynamics of continuum phononic materials with periodic rough contacts,

    G. U. Patil and K. H. Matlack, “Strongly nonlinear wave dynamics of continuum phononic materials with periodic rough contacts,” Physical Review E, vol. 105, no. 2, p. 024201, 2022

    work page 2022

  17. [17]

    A multiscale model for weakly nonlinear shallow water waves over periodic bathymetry,

    D. I. Ketcheson, L. L´ oczi, and G. Russo, “A multiscale model for weakly nonlinear shallow water waves over periodic bathymetry,” arXiv preprint arXiv:2311.02603 , 2023

    work page arXiv 2023

  18. [18]

    Solving Boundary-Value Problems for Systems of Hyperbolic Con- servation Laws with Rapidly Varying Coefficients,

    Yong, D. H., and J. Kevorkian, “Solving Boundary-Value Problems for Systems of Hyperbolic Con- servation Laws with Rapidly Varying Coefficients,” Studies in Applied Mathematics , vol. 108, no. 3, pp. 259–303, 2002

    work page 2002

  19. [19]

    The Riemann problem for fluid flows in a nozzle with discontinuous cross-section,

    P. G. Lefloch and M. D. Thanh, “The Riemann problem for fluid flows in a nozzle with discontinuous cross-section,” Communications in Mathematical Sciences , vol. 1, no. 4, pp. 763–797, 2003. 16

    work page 2003

  20. [20]

    Nonlinear waves

    B. Deconinck, “Nonlinear waves.” Department of Applied Mathematics, University of Washington, February 2018

    work page 2018

  21. [21]

    Efficient implementation of essentially non-oscillatory shock-capturing schemes, ii,

    C.-W. Shu and S. Osher, “Efficient implementation of essentially non-oscillatory shock-capturing schemes, ii,” Journal of Computational Physics , vol. 83, no. 1, pp. 32–78, 1989

    work page 1989

  22. [22]

    R. J. LeVeque, Finite Volume Methods for Hyperbolic Problems . Cambridge University Press, 1 ed., 2002

    work page 2002

  23. [23]

    A wave propagation method for con- servation laws and balance laws with spatially varying flux functions,

    D. S. Bale, R. J. LeVeque, S. Mitran, and J. A. Rossmanith, “A wave propagation method for con- servation laws and balance laws with spatially varying flux functions,” SIAM Journal on Scientific Computing, vol. 24, no. 3, pp. 955–978, 2003. 17

    work page 2003