{"id":"97bcd2d4-7356-421c-834b-597c1ff5ef3e","arxiv_id":"2412.12779","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":7,"one_line_summary":"A new OpenFOAM solver couples variable-compressibility, thermal, cavitating lubricant flow with a structural solver to simulate elastohydrodynamically lubricated line contacts under different slip conditions.","lead":"This paper describes a new OpenFOAM solver for lubricant flow between surfaces under gigapascal pressures, adding variable liquid compressibility and a temperature equation to an existing cavitation solver. A smart generalist might read it to see how open-source CFD, structural mechanics, and coupling tools can be combined to simulate elastohydrodynamic lubrication in bearings and gears.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The only validation of the new thermal/variable-compressibility features compares the solver against a co-authored near-identical model; the independent comparison is isothermal with different rheology, so the distinctive validation claim lacks an independent anchor.","rationale":"The reader's verdict is CONDITIONAL, and my analysis supports that verdict rather than moving it. The reader's stated weakest assumption concerned the accuracy of the external constitutive parameters (Bair, Björling) in the extreme EHL regime. My concern is related but distinct: even if those parameters are correct, the paper's thermal/variable-compressibility validation is not independent. The only case exercising the new features is compared with Havaej et al., a model using almost identical constitutive choices and sharing two co-authors with the present paper, so the agreement in Fig. 9 could reflect shared implementation choices rather than physical correctness. The independent comparison in Fig. 10 is valuable and should be credited, but it deliberately uses different lubricant models (Dowson/Roelands) and is isothermal, so it does not test the new temperature equation or variable liquid compressibility. Thus there is a real gap between the evidence presented and the broad validation claim in Section 5. I do not see an internal inconsistency or a fatal numerical flaw in the derivation; the paper is a plausible engineering contribution. However, the conditional element should be explicit: acceptance of the full claim should require an independent thermal/variable-compressibility benchmark or experimental comparison. The reader's CONDITIONAL verdict already captures this, so I recommend no change to the verdict label.","tokens_in":20097,"tokens_out":5014,"duration_ms":52982,"concrete_test":"Run the Section 4.2 thermal case (100 kN/m, SRR 1, squalane with the Table 1 parameters) against an independent reference: either a TEHL solver from another research group using the same constitutive models, or published experimental measurements of central film thickness and friction for squalane under comparable loading and sliding conditions. Compare central film thickness h_c, minimum film thickness h_min, pressure spike amplitude, and friction coefficient. If the independent agreement is not within the tolerance implied by Fig. 9 (e.g., more than a few percent in h_c and h_min), then the Section 5 'successfully validated' claim is not supported for the thermal/variable-compressibility regime.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's central conclusion in Section 5 ('successfully validated against results found in literature') rests on two comparisons in Section 4.2. The first comparison, which alone exercises the genuinely new features of the solver — variable liquid compressibility, the temperature equation, shear-thinning and the squalane parameters from Table 1 — is against Havaej et al. [34], a reference explicitly described as using 'almost identical lubricant models and parameters' and co-authored by two of the present authors (P. Havaej and D. Fauconnier). The second comparison (Srirattayawong [45] and Tošić et al. [19]) is independent, but it is an isothermal pure-rolling case run with Dowson compressibility and Roelands piezoviscosity, so it validates the base flow solver and partitioned coupling machinery rather than the thermal and variable-compressibility extensions that distinguish this work. Consequently, the distinctive claim — that the solver accurately models thermal EHL with variable liquid compressibility — is currently supported only by comparing a new code against a near-twin code from the same group. Agreement in Figure 9 is therefore much weaker evidence than agreement against an independent numerical implementation or experimental data. This is the load-bearing gap: the distinctive validation claim is anchored to a non-independent reference, and the only independent comparison does not test the novel physics. The absence of released code and of a grid/time-step convergence study further prevents an independent check of the presented results.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper presents a new OpenFOAM-based two-phase flow solver for elastohrodynamic lubrication (EHL) line contacts. The solver extends cavitatingFoam with variable liquid compressibility via the Tait equation, a temperature equation, Doolittle/Carreau rheology, thermal-conductivity and heat-capacity models, and a homogeneous equilibrium cavitation model. It is coupled in a partitioned way to the Kratos Multiphysics structural solver through the CoCoNuT coupling tool. The manuscript describes the governing equations, the PIMPLE-based implementation, and a line-contact test case with squalane. Validation in Section 4.2 compares pressure and film-thickness profiles against Havaej et al. [34], Srirattayawong [45], and Tošić et al. [19], and Section 4.3 presents results for three slip-to-roll ratios.","tokens_in":20514,"tokens_out":6554,"duration_ms":65455,"significance":"If the solver is adequately validated, it would be a useful open-source tool for CFD-based EHL simulation, combining modular constitutive models, cavitation, thermal effects, and partitioned FSI with entirely open-source components. A clear strength is that all constitutive parameters are taken from external literature sources (Bair [42], Björling et al. [27]); no target result is fitted in this paper, so the core derivation is not circular in a parameter-fitting sense. The governing equations and implementation steps are presented in sufficient detail to be reproduced. However, the significance is currently limited by a validation strategy whose distinctive thermal/variable-compressibility claim rests on comparison with a near-identical, co-authored solver, with no quantitative error metrics or convergence study.","major_comments":[{"comment":"The only case exercising the new solver's distinguishing features—variable liquid compressibility, the temperature equation, shear thinning, and the squalane parameters of Table 1—is compared against Havaej et al. [34], a model that shares two authors with the present paper and is explicitly described as using 'almost identical lubricant models and parameters.' Agreement between two implementations built on the same models by the same group is not independent evidence that the new physics is captured correctly, and no quantitative error metric is reported; the comparison is visual only. Please add a comparison against an independent numerical solver or experimental data for a thermal EHL case, or at minimum report error norms and a sensitivity study of the boundary-condition differences (solid energy equation, Carslaw–Jaeger condition) that are invoked to explain the discrepancies.","section":"§4.2, Fig. 9"},{"comment":"No grid-refinement or time-step convergence study is reported for this solver. The mesh size and time step are justified by reference to Tošić et al. [19], but that work used different discretization techniques and did not include the new thermal and variable-compressibility terms. Since the central claim is quantitative agreement of pressure and film thickness in a thin-film contact with very high pressure gradients, the absence of a convergence check leaves open the possibility that discretization error contributes to the observed agreement or disagreement. Please add a convergence study (e.g., central and minimum film thickness, maximum pressure, and the Petrusevich spike) or provide a quantitative argument that the present resolution is sufficient for the new terms.","section":"§3.3, §4.2"},{"comment":"The validation set does not directly test the temperature solution. Figure 9 plots only pressure and film thickness, despite the thermal nature of that case, and Figure 10 is isothermal. Because the temperature equation is one of the two main novelties of the solver, the temperature field itself, or a directly temperature-sensitive integrated quantity such as friction, should be compared with an independent result. Without this, the thermal model is not independently anchored.","section":"§4.2, Figs. 9–10"},{"comment":"The conclusion states that the 'resulting solution setup is successfully validated against results found in literature.' As it stands, this overstates the evidence: the thermal/variable-compressibility validation is anchored to a non-independent reference, and the independent comparison is an isothermal case with different rheology that does not exercise the novel physics. Please either add the additional independent validation suggested above or substantially narrow the conclusion to describe what is actually demonstrated.","section":"§5"}],"minor_comments":[{"comment":"The first validation case is described as 'SSR 1' and the caption of Table 2 uses 'SSR'; these should be 'SRR' for slip-to-roll ratio.","section":"§4.2, Table 2"},{"comment":"The sentence 'and at the same the illustration' is missing the word 'time'; it should read 'and at the same time the illustration.'","section":"Abstract"},{"comment":"The temperature equation is presented with correction terms that vanish upon convergence. A short explanation of why the correction terms are needed, and how their magnitude is monitored, would improve readability.","section":"§2.2.7, Eq. (30)"},{"comment":"The fixed-temperature wall condition on the roller is justified by a time-scale argument, but the paper does not quantify the resulting error relative to the Carslaw–Jaeger condition used in the reference case. A quantitative estimate or a sensitivity check would strengthen the comparison in Fig. 9.","section":"§3.5"},{"comment":"The footnote noting that other sources use an exponent of −4 instead of −3 in the heat-capacity model is helpful, but the implications of this difference for the thermal results are not discussed. A one-sentence comment on the sensitivity would be useful.","section":"§2.1.3, Eq. (17)"},{"comment":"The paper states that the solver is intended to benefit the OpenFOAM community, but no repository or code-availability statement is given. A link to a public repository or a clear statement of availability would improve reproducibility.","section":"Generally"}],"recommendation":"major_revision","confidential_remarks":"The main issue for the editor is the independence of the validation: the thermal comparison in Fig. 9 is against a model from the same group with nearly identical constitutive choices, and the genuinely independent comparison in Fig. 10 does not test the novel physics. This is not a sign of bad faith—the overlap is explicitly stated—but it is a load-bearing weakness for a paper whose central claim is validation. The authors may also want to consider whether the journal's standards require quantitative error metrics and a convergence study for a new CFD solver. I recommend major revision rather than rejection because the underlying equations and implementation are coherent and the weakness is addressable by adding independent validation and convergence evidence."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Good to see this. The genuinely new contribution is the OpenFOAM solver itself: cavitatingFoam extended with Tait compressibility, Doolittle viscosity, Carreau shear thinning, thermal conductivity/heat-capacity models, and a temperature equation, all selectable at runtime, and coupled to Kratos via CoCoNuT for partitioned FSI. That is real engineering work, and the paper is clearly written. The implementation details—flux correction, HEM pressure correction, relative-flux handling, even the OpenFOAM version-8 bug note—are useful and show the authors know what they're doing. The SRR 1 and 2 results, with shear bands and 8–29 K temperature rises, look plausible and illustrate why the extensions matter.\n\nThe soft spots are about validation. The distinctive claim—thermal, variable-compressibility EHL with slip—is checked only against Havaej et al., a co-authored reference that uses almost identical lubricant models and parameters. That is a consistency check, not an independent validation. The other comparison, against Tošić et al. and Srirattayawong, is independent but isothermal pure rolling with Dowson/Roelands; it tests the base flow solver and coupling machinery, not the new physics. So the Section 5 sentence 'successfully validated against results found in literature' overstates what is actually shown. There are also no quantitative error metrics, no grid/time-step convergence study, and no released code or data, so readers cannot independently check the curves. These issues are fixable and do not undermine the basic soundness of the solver, but they materially weaken the paper's main claim as written.\n\nCitation pattern is fine: constitutive parameters come from external sources (Bair, Björling et al.) and no target result is fitted here. The Havaej et al. reference is relevant and acknowledged as near-identical; the problem is not that they cite it, but that it is their only thermal validation.\n\nBottom line: a solid engineering contribution that deserves referee time. I would send it to review, but I would ask for supplementary material (code, mesh-convergence data) and for a rewritten validation section. If no independent thermal EHL case can be added, the authors should state explicitly that the thermal extension is self-consistent with a similar in-house solver rather than independently validated.","headline":"A genuinely useful OpenFOAM EHL solver, well implemented and clearly written, but the thermal/variable-compressibility validation rests on a near-twin co-authored reference, so 'validated against literature' is too strong as written.","tokens_in":21027,"tokens_out":2932,"would_cite":true,"duration_ms":26655,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"A two-phase CFD solver with variable liquid compressibility and a temperature equation, coupled to a structural solver, reproduces reference pressure and film-thickness profiles in elastohydrodynamic lubrication and predicts…","keywords":["elastohydrodynamic lubrication","two-phase flow","cavitation","variable compressibility","thermal effects","shear thinning","partitioned fluid-structure interaction","squalane"],"falsifier":"Measure central film thickness and friction coefficient in a controlled rolling-sliding EHL contact at the same nominal load and slip ratios 0, 1 and 2, and compare with the predicted values of about 0.33 μm and friction coefficients from 0.001 to 0.033; deviations well beyond the spread between the two literature benchmarks would show that the constitutive closure fails.","tokens_in":19917,"feed_emoji":"⚙️","tokens_out":10908,"duration_ms":85403,"temperature":0.7,"pith_summary":"Elastohydrodynamically lubricated contacts — the narrow, high-pressure oil films that separate rollers, gears and cams — are usually modeled with reduced Reynolds and Boussinesq equations. This paper argues that a full two-phase CFD solver with variable liquid compressibility, cavitation, shear-thinning rheology and a temperature equation can replace those reductions and still reproduce accepted reference solutions. The solver is validated against two independent numerical benchmarks, one thermal sliding and one isothermal pure rolling, with close agreement in pressure and film-thickness profiles. If the validation holds, it gives engineers a modular route to predict film thickness, load and friction in contacts where slip, temperature and compressibility matter.","feed_headline":"New CFD solver reproduces oil films under gigapascal loads","feed_subtitle":"Bearings, gears and cams depend on sub-micron oil films; a validated compressible, thermal solver now predicts them.","key_machinery":"The load-bearing mechanism is the homogeneous equilibrium model (HEM) for cavitation, in which liquid and vapor share velocity, pressure and temperature and the pressure is pinned at the saturation pressure while the vapor fraction adjusts. Around that core, the solver couples the pressure equation with a compressibility relation from the Tait equation, a Doolittle viscosity law, Carreau shear thinning, a limiting shear stress, and an enthalpy-based temperature equation. The segregated PIMPLE loop — an iterative momentum-prediction/pressure-correction sequence — carries the pressure-velocity coupling, and the cavitation step enforces saturation through a linearized pressure-density relation. The whole sequence runs on a moving mesh and is driven by a quasi-Newton coupling iteration with a structural solver, which is what turns the flow solver into an elastohydrodynamic lubrication solver.","core_discovery":"The paper claims that a pressure-based finite-volume solver for two-phase lubricant flow, built by adding a variable liquid compressibility and a temperature equation to an existing homogeneous-equilibrium cavitation solver, can simulate elastohydrodynamic line contacts at gigapascal pressures without relying on the Reynolds and Boussinesq reductions. The liquid phase is closed by the Tait equation for compressibility, the Doolittle equation for piezoviscosity, the Carreau model for shear thinning and a limiting-shear-stress cutoff, with thermal conductivity and heat capacity depending on the thermodynamic state. Coupled in a partitioned manner to a structural finite-element solver, the framework produces pressure and film-thickness profiles that closely match published numerical results for a thermal sliding contact and for an isothermal pure-rolling contact, and it resolves cavitation, a sevenfold variation of liquid compressibility across the contact, and temperature rises up to about 29 K under slip.","pith_inferences":["The partitioned moving-mesh architecture should extend to point contacts, wavy or rough surfaces, and transient load or speed reversal, although those cases are not demonstrated in the paper.","A combined validation at high slip with independently measured traction would close the gap between the thermal-sliding and pure-rolling benchmark cases, which the paper compares separately.","The fixed-temperature boundary on the rolling solid is justified by slow thermal time scales, but repeated or long-duration contacts would likely need a conjugate heat model.","Because constitutive models are interchangeable, the same solver could be used to quantify model-form uncertainty by comparing Tait-Doolittle-Carreau closures against classical Dowson-Roelands closures, which the paper does not do."],"forward_implications":["Film thickness and pressure profiles for line contacts can be obtained without the Reynolds and Boussinesq assumptions, so inertial and large-deformation effects at the inlet are captured.","Cavitation appears where the pressure drops to the saturation value; the vapor pocket location and density field are outputs of the solver rather than imposed.","The liquid compressibility varies by about a factor of 7 across a pure-rolling contact, so constant-compressibility cavitation solvers would misrepresent the high-pressure region.","Under sliding, viscous heating raises the central temperature by roughly 8 K at slip-to-roll ratio 1 and 29 K at ratio 2, dropping viscosity by an order of magnitude and forming a shear band; the friction coefficient rises from 0.00115 to 0.03331.","Constitutive models can be selected at run time, so the same solver can be retargeted to other lubricants or to other compressibility and viscosity closures without rewriting the flow solver."],"supporting_citations":[{"why":"Supplies the Tait, Doolittle and Carreau model constants for squalane used in the liquid closure.","marker":"[42]"},{"why":"Supplies the thermal conductivity and heat capacity model forms and constants.","marker":"[27]"},{"why":"Provides the thermal sliding EHL reference solution used for the first validation.","marker":"[34]"},{"why":"Provides the isothermal pure-rolling EHL reference solution used for the second validation.","marker":"[45]"},{"why":"Independent CFD EHL study whose mesh size and approach guide the discretization and whose solution is aligned with the pure-rolling validation.","marker":"[19]"},{"why":"Motivates the Tait-Doolittle pairing as a free-volume framework for liquid compressibility and piezoviscosity.","marker":"[26]"},{"why":"Shows partitioned fluid-structure interaction modeling of EHL line contacts, the approach the solver adopts.","marker":"[10]"},{"why":"Supplies the coupling tool and quasi-Newton stabilization used to iterate the partitioned EHL solution.","marker":"[40]"}],"fun_headline_variants":["Solver predicts 29 K rise in sliding EHL contacts","New solver models compressible lubricants with cavitation and heat","Partitioned CFD-FEM solver tackles EHL with compressibility and temperature","Two-phase flow solver for gigapascal EHL contacts validated"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The central claim depends on the constitutive equations and constants for squalane — Tait compressibility, Doolittle viscosity, Carreau shear thinning and the thermal property laws — remaining faithful at pressures up to several gigapascals and high shear rates; if those models or fitted constants are off in that regime, the validated agreement and predicted film thickness and friction would not transfer to other conditions.","fun_headline_variants_meta":{"raw":{"variants":["Solver predicts 29 K rise in sliding EHL contacts","New solver models compressible lubricants with cavitation and heat","Partitioned CFD-FEM solver tackles EHL with compressibility and temperature","Two-phase flow solver for gigapascal EHL contacts validated"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.00081,"raw_usage":{"total_tokens":3536,"prompt_tokens":908,"completion_tokens":2628,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":524,"completion_tokens_details":{"reasoning_tokens":2554}},"tokens_in":524,"tokens_out":2628,"duration_ms":17055,"temperature":1.0,"reasoning_tokens":2554,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-11T13:44:22.702443+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure central film thickness and friction coefficient in a controlled rolling-sliding EHL contact at the same nominal load and slip ratios 0, 1 and 2, and compare with the predicted values of about 0.33 μm and friction coefficients from 0.001 to 0.033; deviations well beyond the spread between the two literature benchmarks would show that the constitutive closure fails.","supporting_citations":[{"cited_title":"Reference liquids for quantitative elastohydrodynamics: selection and rheological characterization,","cited_arxiv_id":null,"evidence_quote":"Supplies the Tait, Doolittle and Carreau model constants for squalane used in the liquid closure."},{"cited_title":"Friction reduction in elastohydrodynamic contacts by thin-layer thermal insulation,","cited_arxiv_id":null,"evidence_quote":"Supplies the thermal conductivity and heat capacity model forms and constants."},{"cited_title":"CFD study of surface roughness effects on the thermo-elastohydrodynamic lubrication line contact problem,","cited_arxiv_id":null,"evidence_quote":"Provides the isothermal pure-rolling EHL reference solution used for the second validation."},{"cited_title":"A computational fluid dynamics study on shearing mechanisms in thermal elastohydrodynamic line contacts,","cited_arxiv_id":null,"evidence_quote":"Independent CFD EHL study whose mesh size and approach guide the discretization and whose solution is aligned with the pure-rolling validation."},{"cited_title":"Revisiting the ASME pressure-viscosity report using the tait-doolittle correlations,","cited_arxiv_id":null,"evidence_quote":"Motivates the Tait-Doolittle pairing as a free-volume framework for liquid compressibility and piezoviscosity."}],"review_version":1}