{"id":"976f075a-a4c6-4909-b9e3-797f55476e8c","arxiv_id":"2411.14900","paper_version":1,"verdict":"ACCEPT","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"Temperature-dependent thermo-viscoelastic wave equations admit global generalized solutions for arbitrary large finite-energy initial data when the elasticity tensors are bounded and uniformly elliptic.","lead":"The authors prove that a standard model of heat generation by acoustic waves in solids, with stiffness that depends on temperature, has solutions that exist for all time even from large initial data. They introduce a new generalized solution concept to handle the loss of the classical energy identity.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"No significant objection identified: Theorem 1.1's proof is internally coherent, and the boundedness assumptions, while restrictive for applications, are explicit and not contradicted by the numerical section.","rationale":"The reader's weakest_assumption identifies boundedness and uniform ellipticity of gamma and Gamma as the least secure premise, and I agree that this is the assumption most likely to fail for physically motivated constitutive laws such as the power law in Section 8. However, that is a scope limitation, not an internal flaw: Theorem 1.1 is explicitly conditional on (1.8) and (1.12)-(1.13), and the paper openly marks the numerical power law as outside those hypotheses. My structural check of the proof found the main argument coherent: the regularized system is globally solvable, the a priori estimates are sufficient for compactness, and the lower-semicontinuity lemmas are applied with consistent signs and admissible parameters. I therefore do not see a reason to alter the ACCEPT verdict. The agreement is partial because the reader's concern is relevant to applicability rather than to the validity of the central theorem.","tokens_in":38280,"tokens_out":29205,"duration_ms":282784,"concrete_test":"Verify Lemma 7.4 by recomputing identity (7.30) and then checking the three liminf applications: Lemma 6.6 with beta=gamma, w=v+a/2 u, eta=kappa/a; Lemma 6.4 with beta=Gamma, w=v, eta=a lambda/kappa; and Lemma 6.4 with beta=gamma, w=u, eta=a^2/(kappa(a+2 mu)). If all three hypotheses (7.16)-(7.18) are satisfied and the signs match, the localized energy inequality (2.5) follows; if any liminf inequality has the wrong sign, the proof fails at that point.","verdict_should_be":"UNCHANGED","load_bearing_attack":"I reconsidered the chain leading to Theorem 1.1. Lemma 3.4 provides the epsilon-independent L-infinity(L2) and L2(H1) bounds for (v_epsilon,u_epsilon); Lemmas 5.1-5.5 supply the Lq bounds on Theta_epsilon and Lr bounds on grad Theta_epsilon plus the dual-time-derivative bound needed for the Aubin-Lions extraction in Lemma 7.3. The weak lower-semicontinuity machinery in Lemmas 6.3, 6.4 and 6.6 is consistent: the sign of each liminf inequality matches the corresponding term in rearrangement (7.30), and the parameter choices (7.16)-(7.18) make the three applications admissible. I found no gap in the passage from (7.20) to (7.31), and no circularity: the regularized problems are global for fixed epsilon by Section 4, and all subsequent estimates are epsilon-independent. The numerical power law (8.33) indeed violates (1.8) and (1.12)-(1.13), but the authors explicitly state that this model is outside Theorem 1.1, and even the bounded exponential law cases that overflow are accompanied by the caveat that overflow may be numerical. The boundedness and uniform-ellipticity assumptions are the most fragile with respect to physical coverage, but they are explicit hypotheses, not hidden steps.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper studies the initial-boundary value problem (1.6) for a Kelvin-Voigt thermoviscoelastic system with temperature-dependent elasticity tensors γ(Θ), Γ(Θ). The main result, Theorem 1.1, asserts global existence of generalized solutions for arbitrary large finite-energy initial data in bounded smooth domains, under the assumptions that γ and Γ are bounded, C², symmetric, and uniformly positive definite. The generalized solution concept in Definition 2.1 consists of a weak formulation for the displacement, a one-sided inequality for the heat equation, and a localized energy dissipation inequality. The proof regularizes the system as (2.12), derives ε-independent estimates, proves global existence of approximate solutions, extracts limits via Aubin-Lions and compactness, and passes to the limit using the lower-semicontinuity lemmas of Section 6. Section 8 contains illustrative FDTD simulations of a one-dimensional resonator, including temperature-dependent elasticity laws that lie outside the theorem's hypotheses.","tokens_in":38431,"tokens_out":42734,"duration_ms":385408,"significance":"The result is significant: it provides a large-data global solvability statement for a model in which the classical Lyapunov structure (1.7) fails, and it introduces a generalized solution concept that is shown to coincide with classical solvability for smooth data (Proposition 2.2). The proof is detailed and largely self-contained; the a priori estimates in Sections 3–5 and the lower-semicontinuity machinery in Section 6 are coherent, and the parameter choices (7.16)–(7.18) are admissible. The paper is also honest about the limitations of the numerical section: it explicitly states that the power law (8.33) lies outside the hypotheses of Theorem 1.1. The main reservation is an algebraic sign inconsistency in the proof of the key energy inequality (2.5), which is correctable but must be fixed.","major_comments":[{"comment":"There is an algebraic sign inconsistency in the central derivation of (2.5). In (7.20), the coefficient of ∫∫ |∇uε|² ζe^{-μt}ψ is −κμ/2, which is the coefficient correctly obtained from the time derivative of Fεζe^{-μt}. However, in (7.30) and in the liminf expression displayed before (7.31), the same coefficient appears as +κμ/2. The final inequality (2.5) requires the plus sign, and the proof can be repaired by multiplying (7.20) by −1 before the rearrangement, but this step is not stated and the chain of displayed equations as written is inconsistent. This is a load-bearing point in the proof of Theorem 1.1 and must be corrected.","section":"Lemma 7.4, Eqs. (7.20)–(7.31)"}],"minor_comments":[{"comment":"The displayed operator A has third row −⟨Γ(Θ):∇sv:∇s·⟩−D∆, which does not match the Θ-equation in (2.12), whose source is quadratic in ∇sv and contains no ∇sΘ term; please clarify the notation or correct the row.","section":"Lemma 2.3"},{"comment":"The sentence 'That u = 0 and Θ = 0 on Ω×(0,∞) immediately results...' should read 'on ∂Ω×(0,∞)' rather than 'on Ω×(0,∞)'.","section":"Proposition 2.2"},{"comment":"The power law (8.33) and the exponential law (8.34) with large b violate the boundedness assumption (1.8); although the text acknowledges this, it would be helpful to state explicitly that the numerical experiments are heuristic and are not intended as a numerical test of Theorem 1.1.","section":"Section 8"},{"comment":"A plus sign appears to be missing before the term κ(a+2μ)/4 |∇u|² ζ(t)e^{-μt}ψ on the left-hand side of the displayed liminf inequality.","section":"Eq. (7.31)"}],"recommendation":"major_revision","confidential_remarks":"The sign inconsistency in §7 appears to be a typographical slip rather than a fatal flaw: the proof can be repaired by multiplying (7.20) by −1 before the rearrangement and adjusting the words 'first five summands' accordingly. The manuscript is otherwise within scope for the journal and the result is substantial. I would support acceptance after the correction is made."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Punchline: the core result is real. For bounded, uniformly elliptic temperature-dependent tensors gamma and Gamma, the coupled wave-heat system has a global generalized solution for any finite-energy data, with no smallness condition. I read through the main estimates and the limit passage, and I did not find a load-bearing gap. The generalized solution concept, joining a weak heat inequality with a localized energy inequality, is a legitimate way to compensate for the loss of the energy identity (1.7), and Proposition 2.2 confirms classical consistency.\n\nWhat is new: previous global existence theory for these systems essentially requires constant coefficients or a favorable structural identity. This paper handles gamma = gamma(Theta) and Gamma = Gamma(Theta) in a genuinely new way. The regularized system is locally well-posed, the a priori bounds in Sections 3-5 are epsilon-independent, and the lower-semicontinuity lemmas in Section 6 are assembled correctly; the parameter choices in (7.16)-(7.18) match the hypotheses of Lemmas 6.4 and 6.6. The numerical section is honest: it uses a power law that violates the theorem's assumptions and says so explicitly, and it flags the overflow ambiguity rather than hiding it.\n\nSoft spots: the global boundedness and uniform positive definiteness of gamma and Gamma are restrictive. Real piezoceramics likely have bounded elasticity, but the experimentally motivated power law (8.33) is outside the theorem. Because that is an explicit hypothesis rather than a hidden step, I treat it as a limitation, not a flaw. Uniqueness and higher regularity are not addressed, which is unsurprising for this kind of generalized framework. The proof invokes some standard tools as black boxes, but they are standard and likely checkable.\n\nThis paper deserves a serious referee. The result is substantive, the argument is carefully structured, and the authors are clear about the gap between the mathematical assumptions and the numerical experiments. I would engage with it.","headline":"Global large-data existence for thermoviscoelasticity with temperature-dependent coefficients is real; the proof is coherent and the boundedness assumptions are explicit, so send it to peer review.","tokens_in":39063,"tokens_out":2560,"would_cite":true,"duration_ms":28407,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["35D99","35L05","74F05","74J10"],"pacs":[],"model":"deepseek-v4-flash","headline":"This paper proves that temperature-dependent elasticity in a Kelvin-Voigt thermo-acoustic system does not obstruct global-in-time generalized solutions for arbitrary initial data.","keywords":["thermoviscoelasticity","Kelvin-Voigt model","acoustic wave heating","temperature-dependent elasticity","global generalized solutions","large data","piezoelectric ceramics","quadratic heat source"],"falsifier":"Run the regularized scheme (2.12) in $n=1$ with a bounded, uniformly elliptic, $C^2$ law such as $\\gamma(\\Theta)=\\Gamma(\\Theta)=1+0.5\\sin^2(\\Theta)$, large initial data such as $u_0=0$, $u_{0t}=M\\sin(\\pi x)$, $\\Theta_0=1$, and check whether the uniform estimates (3.8)-(3.12) and (5.1)-(5.6) hold on a fixed time interval; any violation, or finite-time overflow of the temperature in the limit, would contradict Theorem 1.1. The same test with the bounded exponential law (8.34) and parameters chosen so that $C$ stays within a fixed positive range would show whether the observed overflow is numerical or physical.","tokens_in":37964,"feed_emoji":"🌡️","tokens_out":8213,"duration_ms":77199,"temperature":0.7,"pith_summary":"This paper proves that a Kelvin-Voigt model of acoustic waves generating heat in a solid, with elasticity tensors that depend on temperature, admits global-in-time solutions for arbitrarily large initial data. Previous global results for such thermoviscoelastic systems required the elastic coefficients to be constant, because constancy gives an exact energy-dissipation identity. Here that identity is lost, so the paper introduces a generalized notion of solution in which the temperature equation is replaced by two inequalities, and it shows that smooth approximating solutions converge to such a generalized solution. The result matters for piezoelectric ceramics, whose stiffness measurably changes with temperature and whose resonant heating can push them toward their Curie point. If correct, it says that temperature-dependent stiffness alone does not prevent a global solution; the price is a weaker solution concept and boundedness and uniform-ellipticity assumptions on the tensors.","feed_headline":"Global thermo-acoustic solutions with temperature-dependent stiffness","feed_subtitle":"Large initial data still evolve globally when stiffness depends on temperature, in a generalized sense.","key_machinery":"The argument is carried by the localized energy functional $F=\\frac12|u_t|^2+\\frac{\\kappa}{2}|\\nabla u|^2+\\lambda\\Theta$, for suitably chosen constants $\\kappa>0$, $\\lambda>0$, and $\\mu>0$. Around this functional the authors build a generalized solution concept whose third requirement, (2.5), is a one-sided integrated inequality that is an identity along classical trajectories. To pass from parabolic regularizations to the limit, the proof uses weak lower-semicontinuity lemmas (Lemmas 6.3, 6.4, and 6.6) for expressions of the form $\\langle B(z):\\nabla w,\\nabla w\\rangle$, obtained by representing the positive tensor $B$ through its square root and applying a Lebesgue-type product-convergence lemma. The regularized problem (2.12) adds $-\\varepsilon\\Delta^2 v$ to the velocity equation and $\\varepsilon\\Delta u$ to the displacement equation, which supplies enough smoothing for each approximate solution to be global (Lemma 4.2) and for the compactness extraction in Lemma 7.3.","core_discovery":"The central claim is Theorem 1.1: for bounded $C^2$ symmetric tensors $\\gamma$ and $\\Gamma$ that are uniformly positive definite at every temperature, and for initial data $u_0\\in W^{1,2}_0(\\Omega;\\mathbb{R}^n)$, $u_{0t}\\in L^2(\\Omega;\\mathbb{R}^n)$, and nonnegative $\\Theta_0\\in L^1(\\Omega;\\mathbb{R})$, the problem (1.6) has a global generalized solution in the sense of Definition 2.1 on any bounded smooth domain. The solution satisfies $u\\in L^\\infty_{\\mathrm{loc}}([0,\\infty);W^{1,2}_0(\\Omega;\\mathbb{R}^n))$, $\\Theta\\in L^\\infty_{\\mathrm{loc}}([0,\\infty);L^1(\\Omega))$ with the stated additional $L^q$ and $W^{1,r}$ integrability, and $u_t\\in L^\\infty_{\\mathrm{loc}}([0,\\infty);L^2(\\Omega;\\mathbb{R}^n))\\cap L^2_{\\mathrm{loc}}([0,\\infty);W^{1,2}_0(\\Omega;\\mathbb{R}^n))$. The discovery is that global solvability survives the loss of the energy identity (1.7) caused by temperature dependence: the generalized concept requires the wave equation weakly, the heat equation as a one-sided inequality (2.4), and a localized energy-dissipation inequality (2.5) built around the coupled functional $F=\\frac12|u_t|^2+\\frac{\\kappa}{2}|\\nabla u|^2+\\lambda\\Theta$.","pith_inferences":["If the boundedness or uniform-ellipticity assumptions fail, as they do for the unbounded power law (8.33), Theorem 1.1 gives no information; a natural next step is to ask whether generalized solutions persist under one-sided or unbounded stiffness laws, or whether singularities can actually occur.","The one-sided inequality technique may transfer to other wave-heat systems with $L^1$ temperatures and quadratic sources, such as thermoelasticity including thermal expansion, where the favorable energy identity is also lost.","The hot-spot patterns reported for steep temperature dependence suggest a mechanism worth testing experimentally: if real, they would imply that even ideal cooling cannot stabilize high-power piezoelectric transducers once stiffness rises steeply with temperature; if purely numerical, they still motivate better-posed regularizations."],"forward_implications":["For any bounded smooth domain and any finite-energy initial data, temperature-dependent elasticity tensors that stay bounded and uniformly elliptic do not cause finite-time blow-up in the generalized sense.","The generalized solution concept is consistent with classical solvability: for sufficiently regular data and solutions, the two inequalities collapse to the original system in the classical sense (Proposition 2.2).","The heat-production rate $\\langle\\Gamma(\\Theta):\\nabla_s u_t,\\nabla_s u_t\\rangle$ only needs to be controlled in $L^1$; the theory does not rely on temperature-dependent heat capacity or on additional inelastic variables, unlike several earlier constructions.","In the one-dimensional resonator simulations, temperature-dependent stiffness shifts the resonance frequency and changes the temperature growth from superlinear to sublinear; steep temperature dependence can produce hot spots or numerical overflow, which the paper leaves as an open question.","If the bounded exponential law (8.34) is used, forcing the elasticity to remain within a fixed positive range, the same qualitative detuning behavior appears, suggesting that the core mechanism is robust as long as the boundedness assumptions hold."],"supporting_citations":[{"why":"Provides the previous global weak-solution result for a related thermoviscoelastic system with temperature-independent coefficients, whose favorable energy structure the present paper replaces.","marker":"[6]"},{"why":"Supplies the L1-data framework and Neumann-boundary treatment that the new generalized solution concept builds on and modifies.","marker":"[50]"},{"why":"Gives the general quasilinear parabolic existence theory used to obtain local classical solutions of the regularized problem (2.12).","marker":"[1]"},{"why":"Supplies the Aubin-Lions lemma used to extract convergent subsequences of the regularized solutions.","marker":"[57]"},{"why":"Provides analytic semigroup smoothing estimates used in Lemma 4.1 to make the regularized solutions global.","marker":"[17]"},{"why":"Supplies the matrix square-root construction used in Lemma 6.2 to express dissipative quadratic forms and prove weak lower semicontinuity.","marker":"[24]"},{"why":"Provides Korn's inequality, used in Lemma 3.4 to convert coercivity of $\\gamma$ into control of $\\nabla v_\\varepsilon$.","marker":"[35]"},{"why":"Supplies the experimental measurement of temperature-dependent elastic stiffness that motivates the model and Figure 1.","marker":"[18]"}],"fun_headline_variants":["Global solutions with temperature-dependent stiffness for acoustic heat","Large-data global existence for thermo-acoustic waves","Acoustic heat: global solutions without smallness assumptions","Temperature-dependent elasticity still yields global thermo-acoustic solutions","Generalized global solutions for heat from acoustic waves"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The proof needs both elasticity tensors $\\gamma$ and $\\Gamma$ to be bounded, $C^2$, and uniformly positive definite at every temperature, with the symmetries in (1.9)-(1.11); if stiffness can grow without bound or lose positive definiteness, the theorem says nothing, and the paper's own unbounded power-law simulations overflow for large $k$.","fun_headline_variants_meta":{"raw":{"variants":["Global solutions with temperature-dependent stiffness for acoustic heat","Large-data global existence for thermo-acoustic waves","Acoustic heat: global solutions without smallness assumptions","Temperature-dependent elasticity still yields global thermo-acoustic solutions","Generalized global solutions for heat from acoustic waves"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000825,"raw_usage":{"total_tokens":3620,"prompt_tokens":973,"completion_tokens":2647,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":589,"completion_tokens_details":{"reasoning_tokens":2574}},"tokens_in":589,"tokens_out":2647,"duration_ms":18819,"temperature":1.0,"reasoning_tokens":2574,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T14:46:11.803029+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Run the regularized scheme (2.12) in $n=1$ with a bounded, uniformly elliptic, $C^2$ law such as $\\gamma(\\Theta)=\\Gamma(\\Theta)=1+0.5\\sin^2(\\Theta)$, large initial data such as $u_0=0$, $u_{0t}=M\\sin(\\pi x)$, $\\Theta_0=1$, and check whether the uniform estimates (3.8)-(3.12) and (5.1)-(5.6) hold on a fixed time interval; any violation, or finite-time overflow of the temperature in the limit, would contradict Theorem 1.1. The same test with the bounded exponential law (8.34) and parameters chosen so that $C$ stays within a fixed positive range would show whether the observed overflow is numerical or physical.","supporting_citations":[{"cited_title":"Blanchard and O","cited_arxiv_id":null,"evidence_quote":"Provides the previous global weak-solution result for a related thermoviscoelastic system with temperature-independent coefficients, whose favorable energy structure the present paper replaces."},{"cited_title":"Roubíček","cited_arxiv_id":null,"evidence_quote":"Supplies the L1-data framework and Neumann-boundary treatment that the new generalized solution concept builds on and modifies."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Gives the general quasilinear parabolic existence theory used to obtain local classical solutions of the regularized problem (2.12)."},{"cited_title":"Friedman","cited_arxiv_id":null,"evidence_quote":"Provides analytic semigroup smoothing estimates used in Lemma 4.1 to make the regularized solutions global."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the matrix square-root construction used in Lemma 6.2 to express dissipative quadratic forms and prove weak lower semicontinuity."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides Korn's inequality, used in Lemma 3.4 to convert coercivity of $\\gamma$ into control of $\\nabla v_\\varepsilon$."},{"cited_title":"Friesen, L","cited_arxiv_id":null,"evidence_quote":"Supplies the experimental measurement of temperature-dependent elastic stiffness that motivates the model and Figure 1."}],"review_version":1}