Existence and uniqueness of classical solution to an initial-boundary value problem for the unsteady general planar Broadwell model with four velocities

Amah S\'ena D'Almeida; Koudzo Togb\'evi Selom Sobah

REVIEW 4 major objections 5 minor 19 references

Existence and uniqueness of classical solution to an initial-boundary value problem for the unsteady general planar Broadwell model with four velocities

T0 review · 4 major / 5 minor · reviewed 2026-08-07 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗

Pith's one-line read The unsteady general planar Broadwell model with four velocities admits a unique non-negative classical solution in a rectangle with Dirichlet data, provided a geometric-data smallness condition holds.

desk verdict New result for the theta≠0 planar Broadwell IBVP, but the N2 characteristic partition is wrong and self-contradictory, so the proof as written does not establish the theorem. read the letter →

arxiv 2506.05658 v1 pith:4R3F4SZU submitted 2025-06-06 math.AP

Koudzo Togbévi Selom Sobah , Amah Séna d'Almeida This is my paper

classification math.AP MSC 76A0276M28

keywords discretevelocityBoltzmannmodelsBroadwellmodelinitial-boundaryvalueproblemDirichletboundaryconditionsexistenceanduniquenessSchauderfixedpointtheoremmethodofcharacteristicsclassicalsolution

verification ladder T0 review T1 audit T2 compute T3 formal

The pith

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

The reading

The paper aims to show that the unsteady general planar Broadwell model with four velocities has a unique non-negative classical solution when posed on a rectangular spatial domain with Dirichlet boundary conditions. The main theorem gives two positive parameters, $p$ depending on the dimensions of the space-time domain and $q$ depending on the data, such that if $pq \leq 1/4$, the system admits a unique solution that is continuous, differentiable except possibly on finitely many planes, and has bounded first-order derivatives with explicit estimates. This matters because multidimensional initial-boundary value problems for discrete velocity models have been largely open, and the result extends the earlier $\theta=0$ case to every angle $0 \leq \theta < \pi/2$.

What carries the argument

The central object is the fixed-point operator $T$, defined by the explicit formulas (A.1)-(A.15), which solve the linear system obtained by freezing the nonlinear collision term $Q(M)=2cS(M_2M_3-M_1M_4)$ at a candidate function $M$. For each of the four velocities, the integration is carried backward along the characteristic until it hits either the initial plane or one of the two relevant sides of the rectangle, with the space-time domain partitioned into three regions by characteristic planes. The proof shows $T$ is continuous, lands in the space $E^4$ of functions with bounded derivatives off finitely many planes, satisfies $V(T(M)) \le p (V(M))^2 + q$, and is compact on the balls $B_R$; when $pq \le 1/4$, a radius $R$ exists with $T(B_R) \subset B_R$, so Schauder's theorem supplies the fixed point.

What would settle it

Choose a point in the region where the printed $N_2$ formula uses the $N_1$-type partition yet the backward characteristic of velocity $(-c\sin\theta, c\cos\theta)$ exits through the right boundary $x=b_1$; differentiating the printed $T_2(M)$ along that characteristic and substituting into (2.4) will show whether the residual is zero. A nonzero residual localized near the plane $x+ct\sin\theta=b_1$ would demonstrate that the printed operator does not solve the linearized system, so a fixed point need not satisfy the original PDE.

Watch

Extended reading notes

Core claim

The central claim, stated as Theorem 2.1, is that under the smallness condition $pq \leq 1/4$ the system $\Sigma$ has exactly one non-negative continuous solution $N=(N_1,N_2,N_3,N_4)$ whose first-order partial derivatives exist, are bounded, and are continuous away from finitely many characteristic planes. The existence follows by constructing an operator $T$ from explicit integration along backward characteristics of the linearized system and applying Schauder's fixed point theorem to the ball $B_R$; the uniqueness follows from a uniform Lipschitz estimate showing that the map is a contraction on the set of fixed points with the stated bound. The proof also shows that any such solution is automatically non-negative.

Load-bearing premise

The proof rests on the assumption that the explicit integral formulas defining $T$ really solve the linearized transport equations, which requires the characteristic partition for each velocity to be correct; the printed partition for $N_2$ appears copied from $N_1$'s geometry, and the proof also applies Schauder's fixed point theorem to a relatively compact ball rather than a closed one.

Editorial extensions

If this is right

The system $\Sigma$ is well-posed within the class of non-negative initial and Dirichlet data with continuous and bounded first derivatives, under $pq \le 1/4$.
Every solution carries the explicit a priori bound $V(N) \le (1+\sqrt{1-4pq})/(2p)$, so the densities and their derivatives cannot exceed a quantity fixed by the domain and the data.
Non-negativity is a consequence of the dynamics, not an extra assumption: the fixed point of $T$ is automatically non-negative because the modified linear systems used in the proof have non-negative solutions.
The solution is classical up to finitely many planes, with the same regularity as the data; this gives a precise sense in which the multi-dimensional IBVP is solvable.
The uniqueness proof shows the Lipschitz constant $p'(1+\sqrt{1-4pq})/p$ is strictly below $1$, so no second solution with the stated bounds can exist.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

If the characteristic partition printed for $N_2$ is genuinely copied from $N_1$'s geometry—the backward characteristics for velocity $(-c\sin\theta, c\cos\theta)$ start at $(x+ct\sin\theta, y-ct\cos\theta)$, so the boundary region should be $x+ct\sin\theta \le b_1$, $y-ct\cos\theta \ge a_2$—then the printed operator $T_2$ does not integrate the correct transport equation, and a fixed point of the
A direct repair would recompute the partition for each velocity independently; the same fixed-point estimates would then likely yield the same existence and uniqueness conclusion with $pq \le 1/4$.
The use of Schauder's theorem on the relatively compact ball $B_R$ can be made airtight by applying the theorem to the closed convex hull of $T(B_R)$; this is a standard adjustment that would preserve the conclusion if the estimates hold.
The characteristic-integral formulas give a constructive iteration $M \mapsto T(M)$ that could serve as the basis of a numerical scheme for the unsteady Broadwell system, with the explicit bounds providing a stopping criterion.

Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Desk Editor's Note

read the letter

Short version: the paper really does treat theta ≠ 0 for the general four-velocity planar Broadwell model, which is new compared with the authors' earlier theta=0 paper and the rest of the cited literature. The architecture is standard — characteristic integration, a fixed point operator, Schauder plus contraction — and the bounds involving p and q are explicit and self-contained. I want to give credit for that.

The problem is that the operator T2 is built from the wrong characteristic regions. For N2, the velocity is (−c sinθ, c cosθ), so a backward characteristic from (t, x, y) lands at (x + c t sinθ, y − c t cosθ) at time zero, and the inflow sides are x=b1 and y=a2. Appendix A.5 and Section 3.15 copy the N1 regions instead: they use {x − c t cosθ ≥ a1, y − c t sinθ ≥ a2} and the two N1 boundary sets. The integrands in A.7 and A.8 use N2's boundary data, so the formula is internally inconsistent. Appendix B.11 then gives the correct planes for N2, contradicting A.5. A fixed point of the displayed T is therefore not known to satisfy equation (4.2). This is a load-bearing defect, not a typo, because the Schauder argument only produces a fixed point of this misdefined operator.

Two more issues, in decreasing order of severity. Schauder is invoked on B_R, which is convex and relatively compact but not shown to be closed in C(P; R^4); E is not closed under the uniform norm, so this needs an extra step. And uniqueness is proved only for solutions with V(N) below the upper root of p r^2 − r + q; the inequality allows solutions larger than that, so the claimed global uniqueness does not follow.

I cannot recommend accepting the paper as is. But I would not desk reject it either: the result is new, the strategy is coherent, and a referee can pinpoint the exact repair. Send it to review, and let the authors fix the N2 partition and the closedness gap.

Referee Report

4 major / 5 minor

Summary. The paper studies the initial-boundary value problem for the general planar four-velocity Broadwell model on a rectangle, with non-negative bounded initial and Dirichlet data having first-order derivatives continuous and bounded. The main theorem (Theorem 2.1) asserts the existence of positive parameters p and q such that if pq ≤ 1/4, the system admits a unique non-negative continuous solution, differentiable except possibly on finitely many planes, with bounded derivatives and explicit bounds. The proof constructs an operator T by solving the linearized system along characteristics (Appendix A), proves continuity and mapping estimates for T, and invokes Schauder's fixed point theorem to obtain a fixed point in an invariant ball, followed by a contraction-type uniqueness argument.

Significance. If correct, the result would extend the authors' earlier work for θ = 0 to the general four-velocity planar Broadwell model, contributing to the sparse literature on multi-dimensional discrete-velocity initial-boundary value problems. The paper is self-contained and provides explicit bounds, which is a positive feature. However, the proof as written contains a fundamental error in the characteristic decomposition for the second velocity component, which invalidates the identification of fixed points with solutions. Additional gaps in the Schauder argument and in the uniqueness claim prevent the main theorem from being established in its present form.

major comments (4)

[Appendix A, Eq. (A.5); Section 3, Eq. (3.15)] The characteristic partition for T2(M) is identical to that for T1(M), but N2 satisfies ∂tN2 − c sinθ ∂xN2 + c cosθ ∂yN2 = −Q, whose backward characteristics are (x + c s sinθ, y − c s cosθ). The correct partition, which coincides with the one used in Appendix B, Eq. (B.11), is based on the planes x + ct sinθ = b1 and y − ct cosθ = a2, not on x − ct cosθ = a1 and y − ct sinθ = a2. As written, T2(M) is not the solution of the linearized equation (4.2) with boundary data (2.10)–(2.11), so a fixed point of the displayed operator need not satisfy problem Σ. Since the Schauder fixed point obtained in Section 4.2 is only a fixed point of this misdefined operator, the existence claim of Theorem 2.1 is not supported.
[Section 4.2, Lemma 2; Proposition 4.2] Schauder's theorem is applied to the set BR = {N ∈ E^4 : V(N) ≤ R}, which is not closed in C(P; R^4). The subspace E is defined by differentiability and boundedness of derivatives on the complement of finitely many planes; these conditions are not stable under uniform convergence. Proposition 4.2 proves only that BR is relatively compact, and Lemma 2 asserts compactness of T on BR based on T(BR) lying in a relatively compact set. Neither statement supplies the closedness required by the standard Schauder theorem (or by the version quoted in Theorem 4.1). Without taking the closure of BR or an equivalent argument, the existence of a fixed point does not follow.
[Section 4.2, final paragraph] The uniqueness argument proves that two solutions M and N satisfying ∥M∥, ∥N∥ ≤ (1 + √(1−4pq))/(2p) coincide. Theorem 2.1, however, asserts uniqueness of the solution without any such size restriction. The paper does not show that every solution lies in this invariant ball, so the uniqueness statement as written goes beyond what is proved. The theorem either needs to be weakened to uniqueness among solutions in the ball or supplemented with a global a priori bound.
[Eqs. (4.44)–(4.48)] The inequality β > 2α is stated without proof, and it is not evident from the definitions of α and β that it holds for all parameter values. For example, if the spatial domain is large or c is large, the inequality may fail. This inequality is needed to conclude p′/p ≤ 1/2 and hence to complete the contraction argument. Without a proof, or an explicit condition under which it holds, the uniqueness step is incomplete.

minor comments (5)

[Abstract] The abstract contains a typo: 'fourh velocities' should read 'four velocities.'
[Section 2.2] The notation for the indicator function I used in formulas (3.11)–(3.26) is not defined in the text; the reader is left to infer its meaning from the displayed conditions.
[Section 4.1] The statement 'It is immediate that the solutions of Σ are the fixed points of the operator T' is not immediate in light of the characteristic partition issue identified above; the equivalence requires the operator to be the correct solution operator of the linearized system.
[Throughout] The text alternates between 'positive solution' and 'non-negative solution'; the theorem states non-negative, so the terminology should be made consistent.
[References] Reference [16] is cited as having a DOI and a HAL identifier; if it is not yet published, the citation should be updated or marked as forthcoming.

Circularity Check

0 steps flagged · score 0.0 of 10

No significant circularity: p and q are explicit data/domain constants, T is built from characteristic integrals, and fixed points are verified against the system rather than being fitted to the conclusion.

full rationale

The derivation is self-contained rather than circular. The existence proof constructs an operator T by explicit characteristic integration formulas (A.1)-(A.15) for the linearized system (4.1)-(4.4), and then applies Schauder's fixed-point theorem to a convex, relatively compact set B_R. The parameters p and q are not fitted to the conclusion: they are defined by closed-form expressions (4.39)-(4.40) involving only the domain dimensions, the time horizon T, the physical constants c and S, and norms of the initial and boundary data. The inequality pq <= 1/4 is then used to find an invariant ball, which is a standard fixed-point argument, not a redefinition of the solution. The uniqueness proof is an independent contraction estimate (4.41)-(4.43) using the explicit Lipschitz constant p'. The paper cites its own prior work [16] for the theta=0 case and for 'classical arguments', but the positivity proof in Section 3 and the continuity proof in Lemma 1 are reproduced in the present paper, so the self-citations are not load-bearing reductions of the central claim. The skeptical observation about Appendix A.5 is a genuine correctness concern: the characteristic regions displayed for N2 repeat the N1 partition, whereas the backward characteristics for (4.2) point in the direction (x + c t sin(theta), y - c t cos(theta)), and Appendix B.11 uses the correct planes x + c t sin(theta) = b1 and y - c t cos(theta) = a2. If A.5 is wrong, then T(M) as written may not be the true solution operator of (4.2) with data (2.10)-(2.11), and a fixed point of the displayed T need not satisfy the PDE. However, this is a mathematical error in the verification that fixed points solve the system, not an instance of circular reasoning: no prediction is defined in terms of the conclusion, no fitted parameter is renamed as a prediction, and no load-bearing argument reduces to a self-citation chain. Similarly, applying Schauder's theorem to the relatively compact but not closed set B_R is a technical gap, not circularity. Therefore the paper's central derivation is not circular, and the appropriate circularity score is 0.

Assumptions & free parameters 0 free parameters · 7 assumptions · 0 invented entities

The theorem rests on standard fixed point and compactness tools plus a smallness condition; no fitted parameters or new postulated entities are introduced.

assumptions (7)

standard math Schauder fixed point theorem
Used in Section 4.1, Theorem 4.1, to obtain a fixed point of the operator T.
standard math Arzelà-Ascoli theorem
Used in Proposition 4.2 to conclude that B_R is relatively compact in C(P;R^4).
standard math Existence and uniqueness for linear transport equations with continuous coefficients via characteristic integration
Used throughout Section 3 and Appendix A to define the solution T(M) of the linearized system.
domain assumption Compatibility conditions (2.16)-(2.27) ensure the piecewise characteristic formulas define a continuous function
Invoked in Section 2.2 when defining the initial-boundary value problem and in the explicit formulas for T.
domain assumption Non-negative initial and boundary data with bounded continuous first derivatives
Assumed in Section 2.2 and used in Theorem 3.1 for positivity and in the estimates of Section 4.
domain assumption The velocity parameters satisfy theta in (0, pi/2), c > 0, S > 0, and the domain is a rectangle over a finite time interval
Stated in Section 2.1 and 2.2; theta=0 is covered by the authors' prior paper [16].
domain assumption Smallness condition pq <= 1/4
This is the sufficient condition in Theorem 2.1 guaranteeing that the quadratic bound V(T(M)) <= p V(M)^2 + q admits an invariant ball.

how reviews work

0 comments

Cite this review

Pith. "Pith review of Existence and uniqueness of classical solution to an initial-boundary value problem for the unsteady general planar Broadwell model with four velocities." pith.science (2026). https://pith.science/paper/4R3F4SZU

@misc{pith2026250605658,
  author       = {Pith},
  title        = {Pith review of: Existence and uniqueness of classical solution to an initial-boundary value problem for the unsteady general planar Broadwell model with four velocities},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/4R3F4SZU}},
  note         = {Machine review of arXiv:2506.05658}
}

read the original abstract

We consider the unsteady problem for the general planar Broadwell model with four velocities in a rectangular spatial domain over a finite time interval. We impose a class of non-negative initial and Dirichlet boundary data that are bounded and continuous, along with their first-order partial derivatives. We then prove the existence and uniqueness of a non-negative continuous solution, bounded together with its first-order partial derivatives, to the initial-boundary value problem.

Discussion (0). Sign in to comment.

Reference graph

Works this paper leans on

19 extracted references · 17 canonical work pages

[1]

Broadwell, Shock structure in a simple discrete velocity gas

J.E. Broadwell, Shock structure in a simple discrete velocity gas. Phys. Fluids. 7 (1964), 1243–1247

work page 1964
[2]

Broadwell, Study of rarefied shear flow by the discrete velocity method

J.E. Broadwell, Study of rarefied shear flow by the discrete velocity method. Journal of Fluid Mechanics. 19 (1964), 401–414

work page 1964
[3]

Gatignol, Théorie cinétique d’un gaz à répartition discrète de vitesses

R. Gatignol, Théorie cinétique d’un gaz à répartition discrète de vitesses. Zeit-Fur Flugwissenschaften, 18,(1970), 93–97

work page 1970
[4]

Platkowski, R

T. Platkowski, R. Illner, Discrete velocity models of the Boltzmann equation: a survey on the mathematical aspects of the theory. SIAM. Review. 30 (1988), 213–255

work page 1988
[5]

J. M. Bony, Existence globale à données de Cauchy petites pour les modèles discrets de l’équation de Boltzmann. Communications in Par- tial Diffferential Equations 16 (1991), 533–545

work page 1991
[6]

Cornille, Exact solutions of the broadwell model in 1+1 dimensions J

H. Cornille, Exact solutions of the broadwell model in 1+1 dimensions J. Phys. Fluids. A20 (1987), 1973–1988

work page 1987
[7]

Illner, M

R. Illner, M. Shinbrot, C. Cercignani, A boundary value problem for discrete velocity models. Duke Mathematical Journal 55 (1987), 889– 900

work page 1987
[8]

Cercignani, R

C. Cercignani, R. Illner, M. Shinbrot, A boundary value problem for the two dimensional broadwell model. Commun. Math. Phys. 114 (1988), 687–698

work page 1988

Show all 19 references

[9]

d’Almeida, Exact solutions for discrete velocity models

A. d’Almeida, Exact solutions for discrete velocity models. Mech. Res. Com. 34 (2007), 405–409. 37 K.T. S. SOBAH and A. S. d’ALMEIDA

2007
[10]

d’Almeida, transition of unsteady flows of evaporation to steady state

A. d’Almeida, transition of unsteady flows of evaporation to steady state. C. R. Mecanique 336 (2008), 612–615

2008
[11]

D. R. Smart, Fixed point theorems. Cambridge University Press, 1974

1974
[12]

Nicoué & A.S.d’Almeida, Existence of solutions of a boundary value problem for the four velocity Broadwell model, DOI:10.37418/amsj.10.11.4,November 2021

A.E. Nicoué & A.S.d’Almeida, Existence of solutions of a boundary value problem for the four velocity Broadwell model, DOI:10.37418/amsj.10.11.4,November 2021

2021 doi
[13]

A.S.d’Almeida, Exact Steady Solutions for a Fifteen Velocity Model of GasAugust2020DOI:10.1007/978-3-030-57336-2_9Inbook: Nonlinear Analysis, Geometry and Applications (pp.231-261)

doi
[14]

Defoou & K.K.L

P.L. Defoou & K.K.L. Sossou & A.S. d’Almeida, On the steady solutions of the general four velocity Broadwell, DOI:10.37418/amsj.12.5.2, May 2023

2023 doi
[15]

Temam, Sur la résolution exacte et approchée d’un problème hyper- bolique non linéaire de T

R. Temam, Sur la résolution exacte et approchée d’un problème hyper- bolique non linéaire de T. Carleman. Arch. Rat. Mech. Anal. 35 (1969), 351–362

1969
[16]

Sobah & A.S

K.T.S. Sobah & A.S. D’almeida, On the existence and uniqueness of classical solution for an initial-boundary value problem for a discrete Boltzmann system in two space dimensions. Advances in Mathematics: Scientific Journal, 2025, 14, pp.73 - 102.〈10.37418/amsj.14.1.5〉.〈hal- 05050184〉

2025 doi
[17]

and Luol, L.S., 1980

Cabannes, H., Gatignol, R. and Luol, L.S., 1980. The discrete Boltz- mann equation. Lecture Notes at University of California, Berkley, pp.1- 65

1980
[18]

& Kawashima S., Le problème aux valeurs initiales en théorie cinétique discrète, C

Cabannes H. & Kawashima S., Le problème aux valeurs initiales en théorie cinétique discrète, C. r. hebd. Séanc. Acad. Sci. Paris 307, 507- 511 (1988)

1988
[19]

Publications scientifiques de l’Institut Mittag-Leffler, Uppsala (1957), p

Carleman, T., Problèmes mathématiques dans la théorie cinétique des gaz. Publications scientifiques de l’Institut Mittag-Leffler, Uppsala (1957), p. 104–106. 38

1957

Pith tools

Reviewed August 7, 2026 · model on record in the stance chip above.

[1] [1]

Broadwell, Shock structure in a simple discrete velocity gas

J.E. Broadwell, Shock structure in a simple discrete velocity gas. Phys. Fluids. 7 (1964), 1243–1247

work page 1964

[2] [2]

Broadwell, Study of rarefied shear flow by the discrete velocity method

J.E. Broadwell, Study of rarefied shear flow by the discrete velocity method. Journal of Fluid Mechanics. 19 (1964), 401–414

work page 1964

[3] [3]

Gatignol, Théorie cinétique d’un gaz à répartition discrète de vitesses

R. Gatignol, Théorie cinétique d’un gaz à répartition discrète de vitesses. Zeit-Fur Flugwissenschaften, 18,(1970), 93–97

work page 1970

[4] [4]

Platkowski, R

T. Platkowski, R. Illner, Discrete velocity models of the Boltzmann equation: a survey on the mathematical aspects of the theory. SIAM. Review. 30 (1988), 213–255

work page 1988

[5] [5]

J. M. Bony, Existence globale à données de Cauchy petites pour les modèles discrets de l’équation de Boltzmann. Communications in Par- tial Diffferential Equations 16 (1991), 533–545

work page 1991

[6] [6]

Cornille, Exact solutions of the broadwell model in 1+1 dimensions J

H. Cornille, Exact solutions of the broadwell model in 1+1 dimensions J. Phys. Fluids. A20 (1987), 1973–1988

work page 1987

[7] [7]

Illner, M

R. Illner, M. Shinbrot, C. Cercignani, A boundary value problem for discrete velocity models. Duke Mathematical Journal 55 (1987), 889– 900

work page 1987

[8] [8]

Cercignani, R

C. Cercignani, R. Illner, M. Shinbrot, A boundary value problem for the two dimensional broadwell model. Commun. Math. Phys. 114 (1988), 687–698

work page 1988