{"id":"5ed4666d-3a04-4508-9e21-ef7e94c6cc91","arxiv_id":"2608.01925","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"Numerical simulations of the log-strain viscoelastic model show that the viscosity ratio ξ determines how strongly increasing the Weissenberg number flattens velocity profiles, and that flow-type (shear vs. extension) dependence appears even around a cylinder.","lead":"The authors simulate a recently proposed viscoelastic fluid model in channel, cylinder, and contraction flows, showing that the viscosity ratio, not just the Weissenberg number, controls how non-Newtonian the flow profiles become. The work gives computational rheologists a tested finite element method and clarifies which features of the new model a simpler generalized-Newtonian description cannot capture.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Projection Eq. (18) enforces only symmetry and unit determinant, not positive definiteness, so the central claim that computed Bel remains physically admissible is unsupported.","rationale":"The reader's weakest-assumption analysis is well targeted. The continuous log-strain model indeed preserves SPD, but the discrete method's stated projection does not enforce it. This is directly load-bearing because the constitutive equation and the logarithmic stress require a positive-definite Bel; if the numerical solution leaves the SPD cone, the subsequent logarithms and stresses are not physically meaningful, and the reported flow phenomena could be numerical artifacts. The manuscript provides no diagnostic evidence (minimum eigenvalues, positivity checks, or reproducible code) to close this gap. The concern does not require alleging any dishonesty; it is a precise technical point about the projection in Eq. (18). The paper has genuine strengths: it studies an interesting model, uses a reasonable VMS stabilization, and reports plausible qualitative trends. If the authors add a positivity check or use a true SPD projection and confirm the results are unchanged, the central claims would be substantially supported. Therefore the reader's CONDITIONAL verdict is appropriate, and no stronger rejection is justified without evidence that positivity is actually violated in the reported runs.","tokens_in":21271,"tokens_out":3688,"duration_ms":50807,"concrete_test":"Instrument the solver to record the minimum eigenvalue of Bel over all elements and all time steps for the highest-Wi cylinder and contraction cases (e.g., ξ=25, C=150 or C=10 in Figs. 4 and 7). If min λ(Bel) ≤ 0 at any time or location, the projection does not enforce positive definiteness. As a second check, replace projection (18) with a true SPD projection: symmetrize, eigendecompose, clamp negative eigenvalues to a small positive threshold, renormalize det=1, and rerun a representative case; if velocity profiles or drag coefficients change beyond discretization error, the reported results depend on the unenforced SPD constraint.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The abstract and Section 3.1 assert that the computational method enforces symmetry, unit determinant, and positive definiteness of the elastic strain tensor Bel. The actual projection, Eq. (18), symmetrizes the tensor and rescales it by sqrt(det(sym(Bel))). This guarantees symmetry and det=1 when the determinant is positive, but it does not guarantee positive eigenvalues: a symmetric 2x2 matrix with eigenvalues (-1,-1) has det=1 but is negative definite. If det(sym(Bel)) <= 0, the projection is not even real-valued. The discrete update (17) uses log(Bel^n), so if Bel^n ever loses positive definiteness, the logarithm is undefined or nonphysical and the subsequent update is not a valid evolution of the log-strain model. Since the paper's central physical conclusions — plug-like profiles, drag reduction, and the GNF comparison — all depend on these simulations, the lack of any reported check on the minimum eigenvalue of Bel is a load-bearing gap. Increment norms reaching machine precision only show steady state; they do not detect violation of the SPD constraint. Mesh convergence is likewise asserted without quantitative tables, and no code or data are provided, so the positivity property cannot be independently verified from the manuscript.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript reports a numerical study of the log-strain viscoelastic model of Alrashdi and Giusteri, focusing on the separate roles of the Weissenberg number Wi and the viscosity ratio ξ = κτ_r/η. A stabilized mixed finite element method (equal-order P1-P1 with VMS stabilization and a DG0 strain field) is combined with a backward-Euler/semi-Lagrangian treatment of the upper-convected derivative. The method is applied to three benchmark flows: straight channel, flow past a cylinder, and 4:1 planar contraction. The authors find that ξ controls the degree of shear-thinning and hence the appearance of plug-like profiles, that drag reduction occurs in cylinder flow, that no lip vortices appear in the contraction, and that the model cannot be reduced to a Generalized Newtonian Fluid: comparing the log-strain model with a GNF model fitted to the same simple-shear effective viscosity, they report significant profile differences in the cylinder flow that they attribute to flow-type dependence.","tokens_in":21526,"tokens_out":6607,"duration_ms":85510,"significance":"If the numerical results are reliable, the paper makes a useful contribution by demonstrating that Wi and ξ are independent control parameters for the log-strain model and that the model's behavior in non-viscometric flows is not captured by a shear-rate-dependent viscosity alone. The GNF comparison is a good way to isolate flow-type effects, and reaching Wi ~ 35 with a relatively new tensorial model is a practical achievement. However, the paper's central numerical-verification claim is not supported: the projection used in the discrete update does not enforce positive definiteness, and mesh/time-convergence evidence is only qualitative. The physical conclusions are therefore conditional on additional verification.","major_comments":[{"comment":"The abstract and Section 3.1 state that the computational method enforces symmetry, unit determinant, and positive definiteness of B_el. The projection in Eq. (18) only symmetrizes the tensor and rescales it to unit determinant. A symmetric 2x2 matrix with eigenvalues (-1,-1) has determinant 1 but is negative definite, and Eq. (18) leaves it unchanged. If det(sym(B)) ≤ 0 the projection is undefined. Since Eq. (17) requires log B_el^n, loss of positive definiteness makes the evolution ill-defined and the subsequent results non-physical. The manuscript reports no check on the minimum eigenvalue of B_el. This is load-bearing because all reported flow fields and the GNF comparison depend on the computed B_el being admissible. The authors should either prove that the discrete update plus projection preserves positive definiteness under the conditions used, or add numerical monitoring of the s","section":"Section 3.1, Eq. (18)"},{"comment":"Mesh convergence is asserted qualitatively: Section 4.1.1 states a 64x64 mesh is 'sufficient', and Sections 4.2.1 and 4.3.1 state that convergence was 'verified through successive refinements' without quantitative tables or error norms. Similarly, steady state is defined by increment norms decaying to 'machine precision', but no tolerance, history, or time-step dependence is reported. For a numerical study that reaches Wi ~ 35 and includes a re-entrant corner singularity, quantitative mesh-convergence and time-step studies are necessary to support the reliability of the reported steady-state fields. Please add tables or figures with grid sizes, L2 error norms or Richardson extrapolation, and the convergence history of the increment norms.","section":"Section 4 (general numerical verification)"},{"comment":"The central claim that flow-type dependence is significant in the cylinder flow rests on comparing the log-strain model with a GNF model whose viscosity is an approximate Cross fit (β = 0.788, n = 1.458). The authors acknowledge a mismatch in the intermediate shear-rate region and attribute the 2.63% channel profile discrepancy to it. In the cylinder flow, they argue that the profile differences are not a mere scaling and hence cannot be due to fitting error. This is suggestive but not conclusive: a local shear-rate error in the GNF model could produce shape changes in a non-unidirectional flow. To make the identification robust, the comparison should be repeated using the exact tabulated effective viscosity from the log-strain model in simple shear, or the sensitivity of the cylinder profiles to plausible variations of β and n should be quantified.","section":"Section 4.4, Figures 11-12"}],"minor_comments":[{"comment":"The axis/panel labels appear garbled ('b=25' and '8=...' instead of ξ and Wi). Please check the figure rendering.","section":"Figure 7"},{"comment":"The legend in panel (b) repeats the entry 'Wi=0.0, C=10.0' twice; one of the entries should be a different Newtonian reference (e.g., C=15.0).","section":"Figure 9"},{"comment":"The hat notation is used in Eq. (14) before it is defined after Eq. (20). Define the notation earlier.","section":"Section 3.1"},{"comment":"The Cross fit parameters β and n are given without confidence intervals or a goodness-of-fit measure. Given the importance of the fit to the GNF comparison, report at least one measure such as R² or the mean relative error.","section":"Section 4.4.1"},{"comment":"The code and data are only available 'on request'. For a numerical study of this type, a public repository with the FEniCSx scripts would substantially strengthen reproducibility.","section":"Data availability"}],"recommendation":"major_revision","confidential_remarks":"The paper is likely from a capable group and the qualitative results are consistent with the known behavior of the model. The main issue is not the physical interpretation but the verification of the numerical method: the projection in Eq. (18) does not enforce positive definiteness, and this is claimed in the abstract. I would be willing to accept after the authors add eigenvalue monitoring and quantitative convergence studies. The GNF comparison should also be made robust to the fitting error."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Worth reading if you work on viscoelastic FE or on the log-strain model. The novel content is real: first finite-element simulations of Alrashdi–Giusteri's log-strain model in standard benchmarks, first systematic ξ–Wi exploration, and a clean GNF comparison. The method is competently built—VMS stabilization for u/p, semi-Lagrangian Lie-derivative treatment of Bel, DG0 for the tensor. Results are internally consistent with the model's known shear-thinning behavior: plug-like profiles at high ξ, drag reduction with strain rate, no lip vortices in the 4:1 contraction, and a non-monotonic corner vortex. The GNF baseline is fitted to the model's own effective viscosity, but it is used only as a comparison, so no circularity damages the main claims. The demonstration that flow-type dependence is visible around a cylinder is the paper's best point.\n\nThe stress-test concern is on target. Eq. (18) does not enforce positive definiteness. Symmetrizing and rescaling to det = 1 leaves open the possibility of negative eigenvalues, and the projection is undefined if the symmetrized determinant is non-positive. The update (17) needs log B^n, so any loss of SPD makes the next step invalid. The paper says these properties are 'enforced', but the stated algebra doesn't do that. This is not cosmetic: it is the difference between simulating the log-strain model and simulating a nearby projected process. The fix is cheap—report the minimum eigenvalue of Bel over the domain and time, or check positivity after each projection. Mesh convergence is also asserted rather than shown: a 64×64 channel grid is called sufficient, and the cylinder/contraction meshes are said to be verified, but no tables or quantitative convergence data appear. No code is shared, and the data-availability statement is 'on request', so the positivity property cannot be checked independently.\n\nThe qualitative conclusions may survive: the trends match what the model should do in shear-dominated flows, and the GNF comparison uses steady-state fields, so a modest positivity violation would not obviously change the pictures. But because the positivity claim is load-bearing and unverified, the paper needs revision before a practitioner should trust the solver at Wi ≈ 35. I'd send it to peer review—the numerical method and benchmarks are worth refereeing—but the request must include a demonstrable positivity check and mesh-convergence tables, or access to code.","headline":"First FE benchmarks for the log-strain model with a credible ξ–Wi story, but the claimed enforcement of positive definiteness isn't in the equations.","tokens_in":22053,"tokens_out":4306,"would_cite":true,"duration_ms":52086,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["76A10","65M60","76M10"],"pacs":["47.50.-d","83.60.Bc"],"model":"deepseek-v4-flash","headline":"The paper shows that in a log-strain viscoelastic model, the polymer-to-solvent viscosity ratio independently controls how strongly non-Newtonian flow profiles develop with increasing Weissenberg number, and that the model cannot be reduced","keywords":["viscoelastic fluids","log-strain tensorial model","Weissenberg number","viscosity ratio","shear-thinning flows","generalized Newtonian fluid","flow-type dependence","finite element simulation"],"falsifier":"Check the computed Bel field in any of the reported runs for points where the symmetrized tensor has a non-positive determinant or a negative eigenvalue. The existence of even one such point would mean the asserted preservation of positive definiteness does not hold in the discretization, and the steady states reported at high Wi may reflect the projection's repair rather than the log-strain model's dynamics.","tokens_in":21118,"feed_emoji":"🌊","tokens_out":11713,"duration_ms":102919,"temperature":0.7,"pith_summary":"This paper studies a recently proposed viscoelastic model in which the elastic stress is a logarithmic function of the elastic strain tensor, the log-strain model. Using a stabilized finite-element method, it simulates pressure-driven flow in a straight channel, around a cylinder, and through a 4:1 contraction. It finds that the polymer-to-solvent viscosity ratio, not the Weissenberg number alone, controls how far the steady velocity profiles depart from the Newtonian parabolic shape: at small ratios the profiles stay nearly parabolic for all Wi, while at large ratios they become plug-like and shear-thinning. A second result is that a Generalized Newtonian Fluid fitted to the model's simple-shear viscosity reproduces the channel flow but not the cylinder flow, indicating flow-type dependence that a rate-dependent viscosity cannot describe. The paper concludes that the viscosity ratio must be treated as an independent parameter alongside the Weissenberg number in any study of this model.","feed_headline":"Viscosity ratio sets how strongly Weissenberg reshapes flow","feed_subtitle":"In a log-strain model, the polymer/solvent viscosity ratio, not Wi alone, dictates how far non-Newtonian profiles go.","key_machinery":"The central object is the elastic strain tensor Bel, a symmetric positive definite tensor with unit determinant that plays the role of a conformation tensor; its evolution is the upper-convected derivative of Bel balanced by the relaxation term -(1/tau_r) Bel*log(Bel), and the elastic stress is Tel = kappa*log(Bel). The independent dimensionless parameter xi = kappa*tau_r/eta, the ratio of polymeric to solvent viscosity contributions, is the knob that controls how much shear thinning the model displays. The computational machinery is a staggered semi-implicit time-stepping scheme using a generalized Lie derivative with semi-Lagrangian advection for the strain update, a Variational Multiscale","core_discovery":"The central claim, on the paper's own terms, is that the Weissenberg number and the viscosity ratio xi = kappa*tau_r/eta play independent roles in the log-strain model, with xi determining the extent to which raising Wi produces observably non-Newtonian behavior. In straight-channel flow the profiles interpolate between the solvent-only and zero-shear Newtonian limits, developing a plug-like central region at intermediate Wi that shrinks as Wi grows; this plug-like regime requires sufficiently large xi. In flow past a cylinder the drag coefficient transitions from the zero-shear to the solvent baseline as the effective strain rate rises. In the 4:1 contraction the corner vortex size varies n","pith_inferences":["If this xi-dependence persists in other geometries, then reporting the Weissenberg number alone is insufficient: experiments and simulations should quote the viscosity ratio as well, or the same Wi can correspond to qualitatively different flow fields.","The mismatch between the log-strain and Generalized Newtonian solutions in mixed shear/extension flows could serve as a quantitative probe of flow-type-dependent rheology, potentially distinguishing this model from purely shear-thinning constitutive equations in microfluidic experiments.","A direct way to test the numerical enforcement claim is to record the eigenvalues of Bel during the simulations; if positivity is ever lost locally, the high-Wi results would be artifacts of the projection step rather than predictions of the model.","The predicted absence of secondary lip vortices in the 4:1 contraction could be used experimentally to distinguish strongly shear-thinning log-strain behavior from constant-viscosity viscoelastic models that typically produce lip vortices."],"forward_implications":["In straight-channel pressure-driven flow, increasing xi at fixed Wi turns the parabolic velocity profile into a plug-like shear-thinning profile with a central plug region that shrinks as Wi grows; at small xi the profile stays nearly parabolic for all Wi.","In flow past a cylinder, the drag coefficient drops from the zero-shear Newtonian baseline toward the solvent-only baseline as the effective strain rate increases, with the magnitude of the reduction set by xi.","In the 4:1 planar contraction at large xi, the corner vortex size changes non-monotonically with Wi and no secondary lip vortex appears, in line with strongly shear-thinning behavior.","A Generalized Newtonian Fluid fitted to the simple-shear viscosity of the log-strain model reproduces the channel-flow profile to within a few percent but fails around the cylinder, where extensional deformation makes the log-strain model more dissipative; the model is therefore not equivalent to a shear-rate-dependent viscosity.","The stabilized formulation sustains steady-state computations up to Wi around 35 in these benchmark geometries, beyond the typical limit of many unstabilized viscoelastic formulations."],"fun_headline_variants":["Viscosity ratio, not Wi alone, dictates flow reshaping","Polymer-to-solvent ratio sets Weissenberg's flow impact","In log-strain model, viscosity ratio outweighs Wi alone","Weissenberg reshapes flow only when polymer ratio is right","Non-Newtonian profiles hinge on viscosity ratio, not just Wi"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"The simulations assume that the computed elastic strain tensor Bel stays symmetric positive definite with unit determinant; the code's projection enforces only symmetry and unit determinant, not positive eigenvalues, so the physical validity of the reported solutions depends on an unverified positivity property.","fun_headline_variants_meta":{"raw":{"variants":["Viscosity ratio, not Wi alone, dictates flow reshaping","Polymer-to-solvent ratio sets Weissenberg's flow impact","In log-strain model, viscosity ratio outweighs Wi alone","Weissenberg reshapes flow only when polymer ratio is right","Non-Newtonian profiles hinge on viscosity ratio, not just Wi"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000241,"raw_usage":{"total_tokens":1369,"prompt_tokens":765,"completion_tokens":604,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":509,"completion_tokens_details":{"reasoning_tokens":515}},"tokens_in":509,"tokens_out":604,"duration_ms":6041,"temperature":1.0,"reasoning_tokens":515,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-04T18:05:44.911397+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Check the computed Bel field in any of the reported runs for points where the symmetrized tensor has a non-positive determinant or a negative eigenvalue. The existence of even one such point would mean the asserted preservation of positive definiteness does not hold in the discretization, and the steady states reported at high Wi may reflect the projection's repair rather than the log-strain model's dynamics.","supporting_citations":[],"review_version":1}