On three-dimensional flows of thermo-viscoelastic fluids of Giesekus type
Pith reviewed 2026-05-07 06:51 UTC · model grok-4.3
The pith
Global weak solutions exist for the three-dimensional system of thermo-viscoelastic Giesekus fluids under natural energy and entropy bounds.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
The central claim is that global weak solutions exist in the full three-dimensional setting for the system describing heat-conducting viscoelastic rate-type fluids of Giesekus type. In contrast to existing literature, no smallness, regularity, or structural restrictions on the initial data are imposed beyond natural energy and entropy bounds, and no additional regularising mechanisms such as artificial stress diffusion are required. The analysis relies on a weak-strong framework combining energy and entropy balances with compactness tools to handle the full nonlinear coupling between the fluid velocity, the temperature, and the elastic stress tensor.
What carries the argument
A weak-strong framework that combines energy and entropy balances with compactness tools to control the nonlinear coupling between velocity, temperature, and the Giesekus elastic stress tensor.
If this is right
- The thermodynamically consistent model admits global weak solutions for arbitrary initial data meeting the energy and entropy bounds.
- Heat conduction and temperature-dependent laws can be included without losing the existence result or requiring extra regularization.
- The Giesekus elastic stress remains controllable in the full three-dimensional setting through the weak-strong approach.
- The framework avoids smallness or regularity assumptions that limited earlier analyses of similar systems.
Where Pith is reading between the lines
- The same combination of energy-entropy balances and compactness may extend to other rate-type viscoelastic models with thermal coupling.
- The result supports further study of long-time behavior, stability, and possible asymptotic regimes for these flows.
- Numerical methods for simulating three-dimensional viscoelastic flows with heat transfer can now be grounded in this existence theorem without artificial diffusion.
Load-bearing premise
The energy and entropy balances together with compactness arguments suffice to control the full nonlinear coupling between velocity, temperature, and the Giesekus stress tensor in three dimensions.
What would settle it
Constructing initial data with finite energy and entropy for which no global weak solution exists or the stress tensor loses the required integrability would falsify the existence result.
read the original abstract
Viscoelastic rate-type fluid models constitute a fundamental framework for the mathematical description of complex materials exhibiting coupled elastic and viscous effects, with a wide range of applications in engineering, biomaterials, and medicine. In realistic regimes, thermal effects are essential and lead to strongly coupled systems in which heat conduction and temperature-dependent constitutive laws play a decisive role. In this paper, we develop a thermodynamically consistent model for heat-conducting viscoelastic rate-type fluids. We establish the existence of a global weak solution in the full three-dimensional setting. In contrast to the existing literature, no smallness, regularity, or structural restrictions on the initial data are imposed beyond natural energy and entropy bounds, and no additional regularising mechanisms such as artificial stress diffusion are required. The analysis is based on a weak-strong framework combining energy and entropy balances with compactness tools, allowing us to treat the full nonlinear coupling between the fluid velocity, the temperature, and the elastic stress.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper develops a thermodynamically consistent system of PDEs for incompressible, heat-conducting viscoelastic fluids of Giesekus type in three space dimensions. The model couples the Navier-Stokes equations for the velocity field with a temperature equation and a nonlinear constitutive relation for the symmetric elastic stress tensor that includes an upper-convected derivative and a quadratic Giesekus term, with temperature-dependent viscosity and relaxation time. The central result is the existence of a global weak solution for arbitrary initial data satisfying only natural energy and entropy bounds, obtained via a weak-strong framework that combines integrated energy/entropy inequalities with compactness arguments; no smallness assumptions, artificial stress diffusion, or extra structural hypotheses on the data are imposed.
Significance. If the technical estimates close, the result would constitute a substantial advance over existing literature on thermo-viscoelastic fluids, which typically requires small data, additional regularisation, or restrictive growth conditions on the constitutive functions. The combination of entropy methods with compactness to control the full nonlinear coupling (including temperature-dependent coefficients multiplying the stress transport and quadratic terms) is a technically demanding approach whose success would strengthen the applicability of weak-solution theory to realistic complex-fluid models.
major comments (3)
- [§4.3] §4.3 (passage to the limit in the stress equation): the argument for strong convergence of the Giesekus quadratic term (τ:τ) and the upper-convected derivative relies on uniform integrability of μ(T) and λ(T) times derivatives of τ. The entropy inequality (3.7) yields only ∫|∇√T|² dx dt < ∞ together with control on log T, which does not preclude T oscillating between values where μ(T) or λ(T) become unbounded or degenerate. No explicit growth or Lipschitz assumptions on μ(T), λ(T) are stated in the constitutive hypotheses of §2.2; without them the claimed compactness cannot be justified from the given a-priori bounds.
- [Theorem 3.1] Theorem 3.1 (definition of weak solution) and the integrability class for τ: the space in which τ is sought (L^{4/3} or similar) appears to presuppose that the temperature-dependent coefficients remain bounded, yet the only temperature control comes from the entropy balance. If μ(T) or λ(T) are allowed to grow arbitrarily (consistent with the claim of “no structural restrictions”), the product μ(T)τ may fail to belong to the dual space needed to pass to the limit in the weak formulation of the stress equation.
- [§5.2] §5.2 (compactness via Aubin-Lions or similar): the time-compactness for the stress is obtained from the energy/entropy bounds, but the temperature-dependent transport coefficients destroy the uniform bound on the time derivative of τ unless additional integrability of ∇T or 1/T is available. The entropy estimate does not supply this; a concrete counter-example or additional structural hypothesis on the constitutive functions is needed to close the argument.
minor comments (2)
- [§2] Notation for the Giesekus parameter α and the temperature dependence of the mobility function should be introduced once in §2 and used consistently; several places in §4 switch between α and α(T) without explicit redefinition.
- [Theorem 1.1] The statement of the main existence theorem (Theorem 1.1) lists the initial-data assumptions but omits the precise function spaces for the initial stress τ₀; this should be aligned with the integrability class used in the weak formulation.
Simulated Author's Rebuttal
We thank the referee for the careful reading and constructive major comments, which correctly identify gaps in the stated assumptions on the constitutive functions μ(T) and λ(T). These points are essential for closing the compactness arguments. We will revise the manuscript by adding explicit hypotheses in §2.2 (positive continuous functions with at most polynomial growth compatible with the entropy bounds) and by expanding the justifications in §§4.3 and 5.2. The claim of “no structural restrictions” applies only to the initial data; the constitutive laws may and will now be equipped with the minimal growth conditions needed for the estimates. Below we respond point by point.
read point-by-point responses
-
Referee: [§4.3] §4.3 (passage to the limit in the stress equation): the argument for strong convergence of the Giesekus quadratic term (τ:τ) and the upper-convected derivative relies on uniform integrability of μ(T) and λ(T) times derivatives of τ. The entropy inequality (3.7) yields only ∫|∇√T|² dx dt < ∞ together with control on log T, which does not preclude T oscillating between values where μ(T) or λ(T) become unbounded or degenerate. No explicit growth or Lipschitz assumptions on μ(T), λ(T) are stated in the constitutive hypotheses of §2.2; without them the claimed compactness cannot be justified from the given a-priori bounds.
Authors: We agree that the current statement of §2.2 lacks explicit growth conditions on μ(T) and λ(T), which is a genuine omission. The entropy bound controls log T in L¹ and ∇√T in L², which by Sobolev embedding yields T ∈ L³(Ω×(0,T)) and prevents T from reaching zero too rapidly, but does not automatically give uniform integrability of μ(T)·(τ:τ) or μ(T)·(upper-convected terms) if μ grows arbitrarily. In the revised version we will add the hypothesis that μ and λ are positive, continuous, and satisfy m ≤ μ(T), λ(T) ≤ M(1+T^α) with α small enough that the products remain in L¹ (e.g., α<1). Under this assumption the uniform integrability follows from the already-established L² bound on τ and the L³ integrability of T, allowing the strong convergence of the quadratic and convective terms via Vitali’s theorem. This revision does not alter the main existence result or the absence of restrictions on initial data. revision: yes
-
Referee: [Theorem 3.1] Theorem 3.1 (definition of weak solution) and the integrability class for τ: the space in which τ is sought (L^{4/3} or similar) appears to presuppose that the temperature-dependent coefficients remain bounded, yet the only temperature control comes from the entropy balance. If μ(T) or λ(T) are allowed to grow arbitrarily (consistent with the claim of “no structural restrictions”), the product μ(T)τ may fail to belong to the dual space needed to pass to the limit in the weak formulation of the stress equation.
Authors: The integrability class for τ is derived from the basic energy inequality and is independent of μ and λ. However, the weak formulation of the stress equation requires that μ(T)τ and λ(T)·(upper-convected τ) belong to the dual of the test-function space. We acknowledge that without growth control this may fail. In the revision we will (i) state explicitly in §2.2 the polynomial-growth assumption described above, (ii) verify that under this assumption μ(T)τ ∈ L^{4/3} (or the precise dual space used in Theorem 3.1), and (iii) update the statement of Theorem 3.1 to list the constitutive hypotheses on μ and λ. The phrase “no structural restrictions” in the abstract and introduction refers exclusively to the initial data; it was never intended to apply to the constitutive functions themselves. revision: yes
-
Referee: [§5.2] §5.2 (compactness via Aubin-Lions or similar): the time-compactness for the stress is obtained from the energy/entropy bounds, but the temperature-dependent transport coefficients destroy the uniform bound on the time derivative of τ unless additional integrability of ∇T or 1/T is available. The entropy estimate does not supply this; a concrete counter-example or additional structural hypothesis on the constitutive functions is needed to close the argument.
Authors: The referee is right that the entropy inequality alone does not directly furnish a uniform bound on ∂t τ in a space strong enough for Aubin-Lions when the coefficients depend on T. However, the combination of ∇√T ∈ L² and T ∈ L³ (from Sobolev) already yields ∇T ∈ L^{3/2} and, together with the L² bound on the stress derivatives, produces an integrable right-hand side for the stress equation. To make the argument fully rigorous we will add a short paragraph in §5.2 that (a) recalls the precise integrability of the coefficients under the new growth assumption, (b) shows that the resulting ∂t τ lies in L^{4/3}(0,T; W^{-1,4/3}) or the space required by the chosen compactness lemma, and (c) cites the analogous estimates used in related works on temperature-dependent viscoelastic models. No counter-example is needed once the growth condition is imposed; the revision will therefore close the gap without altering the overall strategy. revision: yes
Circularity Check
No circularity: standard PDE existence proof via energy/entropy estimates and compactness
full rationale
The paper is a pure mathematical existence result for a system of nonlinear PDEs. Its derivation chain consists of a priori energy and entropy bounds, followed by compactness arguments to pass to the limit in a weak-strong framework. No parameters are fitted to data, no quantities are defined in terms of themselves, and no load-bearing step reduces by construction to an input or to a self-citation whose content is the target result. The abstract and described method rely on standard functional-analytic tools (energy/entropy balances plus compactness) that are independent of the specific conclusion being proved. The skeptic's concern addresses whether the stated bounds suffice for the compactness step, but that is a question of proof correctness, not circularity. No quoted equation or self-citation in the provided material exhibits the forbidden patterns.
Axiom & Free-Parameter Ledger
axioms (1)
- standard math A priori energy and entropy estimates plus compactness arguments suffice to obtain global weak solutions for the coupled system
Reference graph
Works this paper leans on
-
[1]
J. M. Ball and F. Murat. Remarks on Chacon’s biting lemma.Proc. Amer. Math. Soc., 107(3):655–663, 1989
work page 1989
-
[2]
J. W. Barrett and S. Boyaval. Existence and approximation of a (regularized) Oldroyd-B model. Math. Models Methods Appl. Sci., 21(9):1783–1837, 2011. ISSN 0218-2025,1793-6314. doi: 10.1142/ S0218202511005581. URLhttps://doi.org/10.1142/S0218202511005581
-
[3]
M. Bathory, M. Bulíček, and J. Málek. Large data existence theory for three-dimensional unsteady flows of rate-type viscoelastic fluids with stress diffusion.Adv. Nonlinear Anal., 10(1):501–521, 2021
work page 2021
-
[4]
Michal Bathory, Miroslav Bulíček, and Josef Málek. Coupling the Navier-Stokes-Fourier equations with the Johnson-Segalman stress-diffusive viscoelastic model: global-in-time and large-data anal- ysis.Math. Models Methods Appl. Sci., 34(3):417–476, 2024. ISSN 0218-2025,1793-6314. doi: 10.1142/S0218202524500064. URLhttps://doi.org/10.1142/S0218202524500064. 48
-
[5]
L. Boccardo and F. Murat. Almost everywhere convergence of the gradients of solutions to elliptic and parabolic equations.Nonlinear Anal., 19(6):581–597, 1992. ISSN 0362-546X,1873-5215. doi: 10.1016/ 0362-546X(92)90023-8. URLhttps://doi.org/10.1016/0362-546X(92)90023-8
-
[6]
L. Boccardo, A. Dall’Aglio, T. Gallouët, and L. Orsina. Nonlinear parabolic equations with measure data.J. Funct. Anal., 147(1):237–258, 1997. doi: 10.1006/jfan.1996.3040. URLhttps://doi.org/10. 1006/jfan.1996.3040
-
[7]
Miroslav Bulíˇ cek and Jakub Wóznicki. On various models describing behaviour of thermoviscoelastic rate-type fluids.Advances in Nonlinear Analysis, accepted, 2026
work page 2026
-
[8]
M. Bulíček, E. Feireisl, and J. Málek. A Navier-Stokes-Fourier system for incompressible fluids with temperature dependent material coefficients.Nonlinear Anal. Real World Appl., 10(2):992–1015, 2009. ISSN 1468-1218,1878-5719. doi: 10.1016/j.nonrwa.2007.11.018. URLhttps://doi.org/10.1016/j. nonrwa.2007.11.018
-
[9]
M. Bulíček, J. Málek, and K. R. Rajagopal. Mathematical analysis of unsteady flows of fluids with pressure, shear-rate, and temperature dependent material moduli that slip at solid boundaries.SIAM J. Math. Anal., 41(2):665–707, 2009. ISSN 0036-1410,1095-7154. doi: 10.1137/07069540X. URLhttps: //doi.org/10.1137/07069540X
-
[10]
M. Bulíček, E. Feireisl, and J. Málek. On a class of compressible viscoelastic rate-type fluids with stress-diffusion.Nonlinearity, 32(12):4665–4681, 2019
work page 2019
-
[11]
M. Bulíček, J. Málek, V. Průša, and E. Süli. On incompressible heat-conducting viscoelastic rate-type fluids with stress-diffusion and purely spherical elastic response.SIAM J. Math. Anal., 53(4):3985–4030, 2021
work page 2021
-
[12]
M. Bulíček, A. Jüngel, M. Pokorný, and N. Zamponi. Existence analysis of a stationary compressible fluid model for heat-conducting and chemically reacting mixtures.J. Math. Phys., 63(5):Paper No. 051501, 48,
-
[13]
ISSN 0022-2488,1089-7658. doi: 10.1063/5.0041053. URLhttps://doi.org/10.1063/5.0041053
-
[14]
M. Bulíček, T. Los, Y. Lu, and J. Málek. On planar flows of viscoelastic fluids of Giesekus type. Nonlinearity, 35(12):6557–6604, 2022
work page 2022
-
[15]
M. Bulíček, T. Los, Y. Lu, and J. Málek. On three=dimensional flows of viscoelastic fluids of Giesekus type.Preprint submitted to Nonlinearity, 2024
work page 2024
-
[16]
Miroslav Bulíček, Josef Málek, and Endre Süli. Weak sequential stability of solutions to a nonisothermal kinetic model for incompressible dilute polymeric fluids.Mathematical Models and Methods in Applied Sciences, accepted(0):1–83, 2026. doi: 10.1142/S0218202526500260. URLhttps://doi.org/10.1142/ S0218202526500260
-
[17]
J. M. Burgers. Mechanical considerations - model systems - phenomenological theories of relaxations and viscosity. InFirst report on viscosity and plasticity, pages 5–67. Nordemann Publishing Company, New York, 1939
work page 1939
-
[18]
The div-curl lemma for sequences whose divergence and curl are compact inW−1,1.C
Sergio Conti, Georg Dolzmann, and Stefan Müller. The div-curl lemma for sequences whose divergence and curl are compact inW−1,1.C. R. Math. Acad. Sci. Paris, 349(3-4):175–178, 2011. ISSN 1631-073X. doi: 10.1016/j.crma.2010.11.013. URLhttps://doi.org/10.1016/j.crma.2010.11.013
-
[19]
T. Debiec and E. Süli. On a class of generalised solutions to the kinetic Hookean dumbbell model for incompressible dilute polymeric fluids: existence and macroscopic closure.Arch. Ration. Mech. Anal., 249(4):Paper No. 43, 54, 2025. ISSN 0003-9527,1432-0673. doi: 10.1007/s00205-025-02115-x. URL https://doi.org/10.1007/s00205-025-02115-x
-
[20]
R. J. DiPerna and P.-L. Lions. Ordinary differential equations, transport theory and Sobolev spaces. Inventiones Mathematicae, 98(3):511–547, 1989
work page 1989
-
[21]
M. Dressler, B. J. Edwards, and H. Ch. Öttinger. Macroscopic thermodynamics of flowing polymeric liquids.Rheologica Acta, 38(2):117–136, 1999. doi: 10.1007/s003970050162. URLhttps://doi.org/ 10.1007/s003970050162. 49
-
[22]
E. Feireisl.Dynamics of viscous compressible fluids, volume 26 ofOxford Lecture Series in Mathematics and its Applications. Oxford University Press, Oxford, 2004
work page 2004
-
[23]
E. Feireisl and J. Málek. On the Navier-Stokes equations with temperature-dependent transport coef- ficients.Differ. Equ. Nonlinear Mech., pages Art. ID 90616, 14, 2006. ISSN 1687-4099,1687-4102. doi: 10.1155/denm/2006/90616. URLhttps://doi.org/10.1155/denm/2006/90616
-
[24]
NečasCenterSeries.Birkhäuser/Springer, Cham, 2022
E.FeireislandA.Novotný.Mathematics of open fluid systems. NečasCenterSeries.Birkhäuser/Springer, Cham, 2022. ISBN 978-3-030-94792-7; 978-3-030-94793-4. doi: 10.1007/978-3-030-94793-4. URLhttps: //doi.org/10.1007/978-3-030-94793-4
-
[25]
E. Feireisl and A. Novotný.Singular Limits in Thermodynamics of Viscous Fluids. Advances in Math- ematical Fluid Mechanics. Birkhäuser/Springer, Cham, 2017
work page 2017
- [26]
-
[27]
J. Hron, V. Miloš, V. Průša, O. Souček, and K. Tůma. On thermodynamics of incompressible viscoelastic rate type fluids with temperature dependent material coefficients.International Journal of Non-Linear Mechanics, 95:193–208, 2017. ISSN 0020-7462. doi: 10.1016/j.ijnonlinmec.2017.06.011. URLhttps: //doi.org/10.1016/j.ijnonlinmec.2017.06.011
-
[28]
D. Hu and T. Lelièvre. New entropy estimates for Oldroyd-B and related models.Commun. Math. Sci., 5(4):909–916, 2007. ISSN 1539-6746,1945-0796. doi: 10.4310/cms.2007.v5.n4.a9. URLhttps: //doi.org/10.4310/cms.2007.v5.n4.a9
-
[29]
M.W. Johnson and D. Segalman. A model for viscoelastic fluid behavior which allows non-affine defor- mation.Journal of Non-Newtonian Fluid Mechanics, 2(3):255–270, 1977. ISSN 0377-0257. doi: https: //doi.org/10.1016/0377-0257(77)80003-7. URLhttps://www.sciencedirect.com/science/article/ pii/0377025777800037
- [30]
-
[31]
P.-L. Lions and N. Masmoudi. Global solutions for some Oldroyd models of Non-Newtonian flows. Chinese Annals of Mathematics. Series B, 21(2):131–146, 2000
work page 2000
-
[32]
J. Málek, V. Průša, T. Skřivan, and E. Süli. Thermodynamics of viscoelastic rate-type fluids with stress diffusion.Physics of Fluids, 30(2):023101, 02 2018. doi: 10.1063/1.5018172. URLhttps://doi.org/ 10.1063/1.5018172
- [33]
-
[34]
Global existence of weak solutions to macroscopic models of polymeric flows
Nader Masmoudi. Global existence of weak solutions to macroscopic models of polymeric flows. Journal de Mathématiques Pures et Appliquées, 96(5):502–520, 2011. ISSN 0021-7824. doi: https: //doi.org/10.1016/j.matpur.2011.04.008. URLhttps://www.sciencedirect.com/science/article/ pii/S0021782411000547
-
[35]
J. Málek and V. Průša. Derivation of equations for continuum mechanics and thermodynamics of fluids. InHandbook of Mathematical Analysis in Mechanics of Viscous Fluids, pages 3–72. Springer, Cham, 2018
work page 2018
- [36]
- [37]
-
[38]
J.G. Oldroyd. On the formulation of rheological equations of state.Proceedings of the Royal Society of London. Series A. Mathematical and Physical Sciences, 200:523–541, 1950. 50
work page 1950
-
[39]
K. R. Rajagopal and A. R. Srinivasa. A thermodynamic frame work for rate type fluid models.Journal of Non-Newtonian Fluid Mechanics, 88(3):207–227, 2000
work page 2000
-
[40]
K. R. Rajagopal and A. R. Srinivasa. On thermomechanical restrictions of continua.Proceedings of the Royal Society A: Mathematical, Physical and Engineering Sciences, 460(2042):631–651, 2004. 51
work page 2042
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.