Entropy-Compatible Reconstruction for High-Weissenberg Viscoelastic Flow

Sai Peng

REVIEW 2 major objections 1 minor 1 cited by

Entropy-Compatible Reconstruction for High-Weissenberg Viscoelastic Flow

T0 review · 2 major / 1 minor · reviewed 2026-07-01 · grok-4.3

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

Pith's one-line read A corrected logarithmic reconstruction selected by a least-damping entropy constraint ensures compatibility with the discrete free-energy balance in high-Weissenberg viscoelastic flows.

desk verdict The paper introduces an entropy-compatible log reconstruction that targets three specific defect mechanisms beyond positivity, with local proofs, but the global stability when inserted into the coupled solver is unanalyzed. read the letter →

arxiv 2606.04005 v1 pith:IGQZSYC7 submitted 2026-05-26 cond-mat.soft

Sai Peng This is my paper · ORCID

classification cond-mat.soft

keywords viscoelasticflowshighWeissenbergnumberlog-conformationentropycompatibilityfree-energybalancetensorreconstructionpositivedefinitetensorsnumericalstability

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

Log-conformation and square-root reconstructions keep conformation tensors positive definite yet can still violate the discrete free-energy balance through Jensen-type entropy bias, exponential amplification of perturbations at high stretch, and sign-indefinite work defects. The paper formulates an entropy-compatible reconstruction principle that selects a corrected logarithmic reconstruction via a least-damping entropy constraint. The correction is local, positive, and obtained by bisection while remaining compatible with coupled velocity-pressure-conformation stepping. Proofs establish the existence of a maximal admissible parameter, convexity of the entropy profile along the logarithmic path, a compatible free-energy estimate, a defect-budget estimate for noncompatible cases, asymptotic inactivity on high-order defects, and a conditional resolution advantage at high stretch. Reproducible diagnostics confirm reduced entropy and work defects relative to standard reconstructions.

What carries the argument

The entropy-compatible reconstruction principle together with the least-damping entropy constraint that selects the corrected logarithmic reconstruction.

What would settle it

A high-Weissenberg simulation in which the corrected reconstruction produces larger entropy defects or polymeric-work defects than the standard logarithmic reconstruction would falsify the claimed compatibility advantage.

Watch

Extended reading notes

Core claim

The paper establishes an entropy-compatible reconstruction principle. A corrected logarithmic reconstruction is selected by a least-damping entropy constraint. This produces a positive, spectrally controlled tensor that satisfies the discrete free-energy balance. The construction admits proofs of the maximal admissible parameter's existence, convexity of the entropy profile along the logarithmic path, a compatible free-energy estimate, a defect-budget estimate for noncompatible reconstructions, asymptotic inactivity on high-order admissible defects, and a conditional high-stretch resolution advantage in log-relative and entropy metrics.

Load-bearing premise

The local bisection correction integrates into the coupled velocity-pressure-conformation time stepping without introducing new instabilities or violating other unstated physical constraints of the underlying discretization.

Editorial extensions

If this is right

The corrected reconstruction preserves a compatible free-energy estimate.
Noncompatible reconstructions are bounded by a defect-budget estimate.
The correction is asymptotically inactive on high-order admissible defects.
A conditional high-stretch resolution advantage holds in log-relative and entropy metrics.
The method remains compatible with existing velocity-pressure-conformation time stepping.

Reading between the lines

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

The approach may permit stable runs at higher Weissenberg numbers with less reliance on artificial stabilization.
The entropy-constraint idea could be tested on other tensor evolution equations that require both positivity and energy consistency.
Benchmark comparisons on standard viscoelastic test cases would quantify whether the bisection overhead is offset by reduced long-time error accumulation.

Share X Bluesky LinkedIn Reddit HN

Signed reviews

No signed human review yet.

Editorial analysis

A structured set of objections, weighed in public.

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

Desk Editor's Note

read the letter

The main takeaway is that this paper defines an entropy-compatible reconstruction principle for the conformation tensor, using a corrected logarithmic map picked by a least-damping entropy constraint and local bisection. It goes beyond standard log-conformation and square-root approaches by addressing Jensen-type entropy bias, exponential amplification in stretched states, and sign-indefinite work defects.

They do a clear job laying out those three mechanisms and then proving several local properties of the fix: existence of the maximal admissible parameter, convexity of the entropy profile along the log path, a compatible free-energy estimate, a defect-budget estimate, asymptotic inactivity on high-order defects, and a conditional resolution advantage at high stretch. The plan to compare against logarithmic, square-root, and linear reconstructions with reproducible diagnostics on entropy defects, work defects, and stress-force errors is also a practical step.

The soft spot is the missing piece on integration. All the proofs stay inside the reconstruction operator. There is no analysis of how the bisection correction behaves inside the full velocity-pressure-conformation time step, especially with the discrete divergence constraint, pressure Poisson solve, or advection. The stress-test concern holds: local control does not automatically prevent sign-indefinite defects from reappearing at the system level, and the paper offers no stability argument for that coupling. Without the full text the actual derivations cannot be checked either.

This is written for people doing computational rheology and viscoelastic CFD, especially those hitting high-Weissenberg breakdowns in polymer processing simulations. A reader already working on conformation tensor reconstructions will find the defect taxonomy and entropy constraint useful to consider.

It should go to peer review. The formal local claims and the focus on a recognized barrier are enough to justify referee time, even if the global behavior and full proofs need checking.

Referee Report

2 major / 1 minor

Summary. The manuscript formulates an entropy-compatible reconstruction principle for the conformation tensor to address nonphysical behavior in high-Weissenberg viscoelastic flow simulations. It identifies Jensen-type entropy bias, exponential amplification of perturbations, and sign-indefinite work defects as mechanisms that can arise even with positive tensors, then proposes a corrected logarithmic reconstruction selected by a least-damping entropy constraint and realized via local bisection. The authors prove existence of the maximal admissible parameter, convexity of the entropy profile along the logarithmic path, a compatible free-energy estimate, a defect-budget estimate for noncompatible reconstructions, asymptotic inactivity on high-order admissible defects, and a conditional high-stretch resolution advantage; reproducible diagnostics compare logarithmic, square-root, and linear reconstructions and verify predicted defects and accumulation.

Significance. If the local reconstruction properties and free-energy compatibility hold and extend to the coupled system, the work would supply a parameter-free, mathematically controlled method for preserving discrete free-energy balance in viscoelastic simulations. The explicit proofs of convexity, existence, and defect budgets together with the reproducible diagnostics constitute clear strengths that could improve reliability at high Weissenberg numbers.

major comments (2)

[Abstract] Abstract: the assertion that the correction is 'compatible with coupled velocity--pressure--conformation time stepping' is not supported by any analysis of interactions with the discrete divergence constraint, pressure Poisson solve, or advection of the conformation tensor; all stated proofs remain local to the reconstruction operator itself.
[Abstract] Abstract: the defect-budget estimate controls local reconstruction defects, yet no argument or numerical test is supplied showing that sign-indefinite polymeric-work defects cannot reappear globally once the corrected tensor is inserted into the full time-stepping scheme, which is load-bearing for the claim of avoiding new instabilities.

minor comments (1)

The abstract is information-dense; separating the list of proved statements from the description of the method would improve readability.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the careful reading and constructive comments. We agree that the abstract overstates the scope of our results and will revise it to reflect that all proofs and estimates remain local to the reconstruction operator.

read point-by-point responses

Referee: [Abstract] Abstract: the assertion that the correction is 'compatible with coupled velocity--pressure--conformation time stepping' is not supported by any analysis of interactions with the discrete divergence constraint, pressure Poisson solve, or advection of the conformation tensor; all stated proofs remain local to the reconstruction operator itself.

Authors: We agree that the stated proofs and estimates are strictly local to the reconstruction and provide no analysis of interactions with the discrete divergence constraint, pressure Poisson solve, or advection of the conformation tensor. The phrasing in the abstract was meant to indicate that the operator can be inserted into existing coupled schemes without altering their structure, but this does not constitute demonstrated compatibility for the full system. We will revise the abstract to remove or qualify the claim. revision: yes
Referee: [Abstract] Abstract: the defect-budget estimate controls local reconstruction defects, yet no argument or numerical test is supplied showing that sign-indefinite polymeric-work defects cannot reappear globally once the corrected tensor is inserted into the full time-stepping scheme, which is load-bearing for the claim of avoiding new instabilities.

Authors: The referee is correct that the defect-budget estimate is local and that the manuscript supplies neither a global argument nor numerical tests in the full time-stepping scheme to show that sign-indefinite work defects cannot reappear. The reproducible diagnostics verify local behavior but do not address global reappearance or new instabilities in the coupled system. We will revise the abstract to avoid implying global stability guarantees. revision: yes

Circularity Check

0 steps flagged · score 0.0 of 10

No circularity; derivation introduces and analyzes a new explicit principle without reduction to inputs

full rationale

The paper formulates a new entropy-compatible reconstruction principle defined directly by a least-damping entropy constraint and local bisection correction. All listed proofs (existence of maximal admissible parameter, convexity along the log path, compatible free-energy estimate, defect-budget estimate) are derived from this definition. No fitted parameters are renamed as predictions, no self-citations are load-bearing for the central claim, and no ansatz or uniqueness result is smuggled in. The chain is self-contained; the result does not reduce to its inputs by construction.

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

Review based on abstract only; no free parameters, new entities, or ad-hoc axioms are explicitly introduced beyond standard domain assumptions of positive-definite symmetric tensors and entropy convexity in viscoelastic models.

assumptions (1)

domain assumption Conformation tensor remains symmetric positive definite under the reconstruction.
Invoked throughout the discussion of log-conformation, square-root, and the proposed corrected reconstruction.

how reviews work

0 comments

Cite this review

Pith. "Pith review of Entropy-Compatible Reconstruction for High-Weissenberg Viscoelastic Flow." pith.science (2026). https://pith.science/paper/IGQZSYC7

@misc{pith2026260604005,
  author       = {Pith},
  title        = {Pith review of: Entropy-Compatible Reconstruction for High-Weissenberg Viscoelastic Flow},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/IGQZSYC7}},
  note         = {Machine review of arXiv:2606.04005}
}

read the original abstract

Log-conformation and square-root reconstructions preserve positive definiteness in high-Weissenberg viscoelastic simulations, but positivity alone does not guarantee compatibility with the discrete free-energy balance. We identify three reconstruction-level mechanisms by which strictly positive tensors can still generate nonphysical behavior: Jensen-type entropy bias, exponential amplification of logarithmic perturbations in highly stretched states, and sign-indefinite polymeric-work defects caused by using incompatible tensors in stress work and entropy variables. We formulate an entropy-compatible reconstruction principle and a corrected logarithmic reconstruction selected by a least-damping entropy constraint. The correction is local, positive, computable by bisection, spectrally controlled, and compatible with coupled velocity--pressure--conformation time stepping. We prove existence of the maximal admissible parameter, convexity of the entropy profile along the logarithmic path, a compatible free-energy estimate, a defect-budget estimate for noncompatible reconstructions, asymptotic inactivity on high-order admissible defects, and a conditional high-stretch resolution advantage in log-relative and entropy metrics. Reproducible diagnostics compare logarithmic, square-root, and linear reconstructions and verify the predicted entropy defects, work defects, stress-force errors, and high-Weissenberg accumulation.

Figures

Figures reproduced from arXiv: 2606.04005 by the authors.

**Figure 1.** Coupled manufactured diagnostic. Positivitypreserving reconstructions generate entropy and stress-force defects that increase with stretch. The corrected log reconstruction damps the incompatible defect before it enters the momentum equation. 0 0.2 0.4 0.6 0.8 1 0 2 4 6 8 ·10−2 t mean entropy standard log physical predictor corrected [PITH_FULL_IMAGE:figures/full_fig_p014_1.png] view at source ↗

**Figure 2.** Entropy correction diagnostic. Positivity is preserved in all logarithmic updates, but the uncorrected reconstruction accumulates additional entropy from the highfrequency defect [PITH_FULL_IMAGE:figures/full_fig_p014_2.png] view at source ↗

**Figure 3.** Dynamic high-Weissenberg reconstruction loop. Positive standard log reconstruction accumulates entropy and stress-force defects, while the corrected log reconstruction damps the incompatible defect at each step [PITH_FULL_IMAGE:figures/full_fig_p015_3.png] view at source ↗

Figures from the paper (3 more)

**Figure 4.** Figure 4: Scalar logarithmic Jensen bias. Both the mean physical stretch and the elastic entropy increase quadratically with the log perturbation amplitude. Appendix J. Extended coupled diagnostics The coupled manufactured diagnostic measures the stress-force defect that would e…

**Figure 5.** Figure 5: Coupled manufactured reconstruction diagnostic. Positive raw log and square-root reconstructions both produce stretch-amplified force defects; the corrected log reconstruction removes the entropy-incompatible component. 0 0.2 0.4 0 2 4 6 ·10−3 t entropy excess raw lo…

**Figure 6.** Figure 6: Dynamic high-Weissenberg reconstruction loop at Wi = 80 and background stretch Λ = 500. The corrected reconstruction suppresses accumulation of the incompatible entropy and force defects. [2] R. B. Bird, R. C. Armstrong, and O. Hassager. Dynamics of Polymeric Liquids, …

Discussion (0). Continue with ORCID to comment.

Forward citations

Cited by 1 Pith paper

Reviewed papers in the Pith corpus that reference this work. Sorted by Pith novelty score. Full citation record

Loss of positive definiteness is a symptom, not the cause, of high-Weissenberg-number breakdown
physics.comp-ph 2026-07 conditional novelty 7.0 of 10

For viscoelastic simulations with solvent viscosity, loss of positive definiteness is neither necessary nor sufficient for breakdown; the discrete stress-coupling route controls survival.

Reference graph

Works this paper leans on

7 extracted references · 7 canonical work pages · cited by 1 Pith paper

[1]

Balci, B

N. Balci, B. Thomases, M. Renardy, and C. R. Doering. Symmetric factorization of the conformation tensor in viscoelastic fluid models.J. Non-Newtonian Fluid Mech., 166:546–553, 2011. 22 SAI PENG 0 500 1,000 10−5 10−2 101 stretch force defect raw log square-root corrected log 0 500 1,000 10−13 10−8 10−3 stretch entropy defect Figure 5.Coupled manufactured ...

work page 2011
[2]

R. B. Bird, R. C. Armstrong, and O. Hassager.Dynamics of Polymeric Liquids, Volume 1. Wiley, 2nd edition, 1987

work page 1987
[3]

Boyaval, T

S. Boyaval, T. Lelievre, and M. Mangoubi. Free-energy-dissipative schemes for the Oldroyd–B model.ESAIM Math. Model. Numer. Anal., 43:523–561, 2009

work page 2009
[4]

Fattal and R

R. Fattal and R. Kupferman. Constitutive laws for the matrix-logarithm of the con- formation tensor.J. Non-Newtonian Fluid Mech., 123:281–285, 2004

work page 2004
[5]

Fattal and R

R. Fattal and R. Kupferman. Time-dependent simulation of viscoelastic flows at high Weissenberg number using the log-conformation representation.J. Non-Newtonian Fluid Mech., 126:23–37, 2005

work page 2005
[6]

M. A. Hulsen, R. Fattal, and R. Kupferman. Flow of viscoelastic fluids past a cylinder at high Weissenberg number: stabilized simulations using matrix logarithms.J. Non- Newtonian Fluid Mech., 127:27–39, 2005

work page 2005
[7]

Keunings

R. Keunings. On the high Weissenberg number problem.J. Non-Newtonian Fluid Mech., 20:209–226, 1986. ENTROPY-COMPATIBLE RECONSTRUCTION 23 School of Mathematics and Computational Science, Xiangtan University, Xiangtan, Hunan, China Email address:pscfd@xtu.edu.cn

work page 1986

Pith tools

Reviewed July 1, 2026 · model on record in the stance chip above.

[1] [1]

Balci, B

N. Balci, B. Thomases, M. Renardy, and C. R. Doering. Symmetric factorization of the conformation tensor in viscoelastic fluid models.J. Non-Newtonian Fluid Mech., 166:546–553, 2011. 22 SAI PENG 0 500 1,000 10−5 10−2 101 stretch force defect raw log square-root corrected log 0 500 1,000 10−13 10−8 10−3 stretch entropy defect Figure 5.Coupled manufactured ...

work page 2011

[2] [2]

R. B. Bird, R. C. Armstrong, and O. Hassager.Dynamics of Polymeric Liquids, Volume 1. Wiley, 2nd edition, 1987

work page 1987

[3] [3]

Boyaval, T

S. Boyaval, T. Lelievre, and M. Mangoubi. Free-energy-dissipative schemes for the Oldroyd–B model.ESAIM Math. Model. Numer. Anal., 43:523–561, 2009

work page 2009

[4] [4]

Fattal and R

R. Fattal and R. Kupferman. Constitutive laws for the matrix-logarithm of the con- formation tensor.J. Non-Newtonian Fluid Mech., 123:281–285, 2004

work page 2004

[5] [5]

Fattal and R

R. Fattal and R. Kupferman. Time-dependent simulation of viscoelastic flows at high Weissenberg number using the log-conformation representation.J. Non-Newtonian Fluid Mech., 126:23–37, 2005

work page 2005

[6] [6]

M. A. Hulsen, R. Fattal, and R. Kupferman. Flow of viscoelastic fluids past a cylinder at high Weissenberg number: stabilized simulations using matrix logarithms.J. Non- Newtonian Fluid Mech., 127:27–39, 2005

work page 2005

[7] [7]

Keunings

R. Keunings. On the high Weissenberg number problem.J. Non-Newtonian Fluid Mech., 20:209–226, 1986. ENTROPY-COMPATIBLE RECONSTRUCTION 23 School of Mathematics and Computational Science, Xiangtan University, Xiangtan, Hunan, China Email address:pscfd@xtu.edu.cn

work page 1986