{"id":"e2b55ee9-a90f-4708-89c1-3ee1dcbaf4fc","arxiv_id":"2411.14286","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"Direct numerical simulations of K-type boundary-layer transition for two supercritical-fluid wall-temperature regimes show delayed, less violent breakdown when the wall exceeds the pseudo-boiling temperature.","lead":"This paper uses detailed computer simulations to study how a flat-plate boundary layer becomes turbulent when the flowing fluid is in a supercritical state. It finds that in the transcritical regime, the breakdown to turbulence is gentler and happens later than in the ideal-gas or subcritical liquid-like regimes.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Causal claim of streak secondary instability is inferred, not demonstrated; delay lacks quantitative metric.","rationale":"The reader's verdict is CONDITIONAL and focuses on the idealized equation of state. I agree the real-fluid transferability is uncertain, but I see a more immediate internal gap: the paper's novel causal mechanism (LSI) is inferred from snapshots and mode amplitudes without a linear secondary-stability validation, and the 'significantly delayed' claim is not accompanied by a defined transition metric. These are exactly the places where the abstract's strongest statements go beyond the presented evidence. The suggested checks are feasible and would not require rerunning the full DNS if the flow fields are archived: one local eigenvalue analysis and one Cf-based onset calculation. The EOS concern remains legitimate but is secondary because the vdW model is the standard tool in this line of work and the qualitative mode-II mechanism has been established with it; the internal causal support is the weaker link. Therefore the CONDITIONAL verdict is retained, with the added requirement to substantiate the LSI mechanism and quantify the delay.","tokens_in":5174,"tokens_out":7585,"duration_ms":74743,"concrete_test":"For case Tw110, extract the streamwise velocity profile in the y-z plane at x/δ99,0=238 (Rex/10^5≈4.35), form the local streak base state, and solve a secondary linear stability problem for varicose and sinuous perturbations scanning spanwise wavenumber and frequency. If no unstable mode is found, the claim that LSI produces the observed hairpin vortices fails. Separately, compute a quantitative transition-onset metric (e.g., the first streamwise location where Cf deviates by 5% from the laminar value) for TadIG, Tw095, and Tw110, and report the values; if the delay is within the grid or forcing-amplitude uncertainty, the 'significantly delayed' claim is not established.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central assertion that localized streak secondary instabilities 'lead to' the simultaneous growth of hairpin vortices and near-wall streaks (abstract; Sec. 4) is supported only by instantaneous flow visualization in Figs. 3-4 and by amplitude thresholds, not by a secondary-instability analysis. In Fig. 4 the symmetric varicose and antisymmetric sinuous modes are identified from the symmetry of u′ contours at a single station (x/δ=238); their causal role is inferred from the appearance of hairpins further downstream in the same instantaneous snapshot. A spatial DNS can suggest instability, but claiming LSI as the mechanism requires showing that the measured streak profile is linearly unstable to those modes at the observed spanwise scales. In addition, the conclusion that transition is 'significantly delayed' is not quantified: no transition-onset Reynolds number based on Cf or amplitude is reported for Tw095/Tw110 against TadIG, and the subcritical-vs-transcritical comparison is made with different Q thresholds and x-ranges in Fig. 3. Because the abstract's mechanism and delay statements are the core of the paper, these omissions are the most load-bearing weaknesses.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper reports direct numerical simulations of controlled K-type breakdown in a flat-plate boundary layer with a supercritical fluid at a reduced pressure of 1.10. Two wall-temperature cases are considered: a liquid-like subcritical case (Tw095) and a transcritical case with wall temperature above the pseudo-boiling line (Tw110), together with an ideal-gas reference case (TadIG) from Sayadi et al. Using a blowing/suction disturbance strip with a fundamental frequency and spanwise wavenumber, the authors compare disturbance-amplitude evolution, flow structures via the Q-criterion, and local vortex interactions. They find that in the subcritical case the K-type breakdown is slightly delayed relative to the ideal-gas reference, with a steady longitudinal mode (0,1) dominating the late transitional stage. In the transcritical case, higher harmonic (2,0) dominates the early nonlinear stage, transition is delayed further, and near the legs of aligned Lambda-vortices the authors observe symmetric varicose and antisymmetric sinuous modes, which they interpret as localized streak secondary instabilities that lead to additional hairpin vortices and near-wall streaky structures. They conclude that transcritical K-type breakdown is more gradual and less violent than the H-type breakdown under the same flow conditions.","tokens_in":5265,"tokens_out":3850,"duration_ms":33905,"significance":"If the main claims hold, this is an important contribution to the understanding of laminar-turbulent transition in supercritical fluids, documenting a qualitatively different breakdown route from the previously reported H-type case. The paper's strengths include a controlled DNS setup matched to an ideal-gas baseline, comparison with linear stability theory, and clear visual documentation of the flow structures. The observation that transcritical K-type transition is delayed and gradual, and that a localized streak secondary instability may be involved, is a falsifiable statement that can guide future stability analyses and engineering predictions. However, the absence of quantitative metrics for the delay and of direct evidence for the instability mechanism currently limits the significance to a descriptive level.","major_comments":[{"comment":"The central claim that transition is 'significantly delayed' in the transcritical case is not quantified. The paper reports no transition-onset Reynolds number, no skin-friction-based criterion, and no direct comparison of Cf or amplitude curves between Tw095, Tw110, and TadIG; Figure 2 only shows Tw095 and Tw110, and the TadIG curves are not shown. In addition, Figure 3 compares the cases with different Q thresholds (0.06 vs 0.016) and different x-ranges, making the visual assessment of delay unreliable. The authors should report a quantitative measure such as Rex at min{Cf}, or a defined transition-onset location, for all three cases, and use consistent visualization parameters for a fair comparison.","section":"Section 4, Figures 2 and 3"},{"comment":"The causal statement that localized streak secondary instabilities 'lead to' the simultaneous development of hairpin vortices and near-wall streaky structures is inferred from a single instantaneous flow field and from an amplitude threshold (the (0,1) mode exceeding 10%), but no secondary-instability analysis is performed. The symmetric varicose and antisymmetric sinuous modes are identified from the symmetry of u′ contours at one station, x/δ=238, and the downstream hairpins are observed in the same snapshot. To support the mechanism, the authors should either (i) perform a local linear stability analysis of the extracted streak profile to show that it is linearly unstable at the observed spanwise scales and frequencies, or (ii) explicitly reword the claim as a plausible conjecture rather than an established mechanism, downplaying it in the abstract and conclusions.","section":"Section 4, Figure 4 and following paragraph"},{"comment":"No grid-convergence or resolution-sensitivity study is presented. The transcritical case uses a substantially different grid (Nx=8400, Δx+=6.1) and a much longer domain (Lx/δ0=570) than the subcritical and ideal-gas cases (Nx=3000, Δx+≈9–10, Lx/δ0≈347–352). Since the conclusions depend on resolving the breakdown processes and on comparing the transition location across cases, the reader cannot rule out that part of the observed delay or change in breakdown route is numerical. A concise grid-convergence test, for instance comparing Cf and amplitude evolution for Tw110 at a coarser resolution, should be included.","section":"Sections 3 and 4, Table 1 and Figures 1–3"}],"minor_comments":[{"comment":"The caption states that 'LST results are shown in black circles', but in the reproduced figure no circles are visible; please verify and correct the caption or the figure.","section":"Caption of Figure 2"},{"comment":"The boundary-layer thickness is introduced as δ∗_99(x0) = δ∗_99,0 and later the local Blasius length scale δ∗ is defined; the notation is close but could be unified to avoid confusion between the inlet thickness and the local scale.","section":"Section 2"},{"comment":"The comparison with the H-type transition relies on Ref. [6], which is a symposium paper; please provide more specific quantitative statements from that comparison (e.g., transition-onset locations, Cf values) rather than only the qualitative remark that the transcritical K-type breakdown is 'not violent'.","section":"Section 3, references to Ref. [6]"},{"comment":"The phrase 'significantly delayed' is used without a quantitative context; consider giving the delay in terms of Rex or x/δ99,0 units to make the claim more precise.","section":"Abstract and Conclusions"}],"recommendation":"major_revision","confidential_remarks":"The paper is a concise DNS study with interesting and potentially important results, but the absence of a quantitative delay metric and of direct evidence for the localized secondary instability mechanism will require a major revision before it can satisfy journal standards. The authors should be asked to add a grid-convergence test, to provide quantitative transition-onset data, and to either perform a local stability analysis or temper the causal language. The strong point is the controlled setup and the contrast between K-type and H-type behavior, which deserves careful, quantified support."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nThis is the first DNS of controlled K-type breakdown in a supercritical-fluid boundary layer, and it does what a good exploratory DNS should: it reuses the well-established Sayadi et al. setup, runs one subcritical and one transcritical wall-temperature case against the ideal-gas reference, and reports a clean contrast. The subcritical case behaves like a delayed version of ideal-gas K-type; the transcritical case breaks down through higher harmonics and a more gradual route, without the violent valley structures seen in the group's earlier H-type study. That contrast is the real contribution.\n\nCredit where due: the amplitude curves and Q-criterion visualizations are internally consistent. The LST diagrams in Fig. 1 anchor the chosen disturbance frequency in the unstable range for the transcritical case, and using the same disturbance strip as Sayadi et al. makes the comparison meaningful. The setup is described in enough detail to be reproduced by someone with access to the CUBENS solver.\n\nThe stress-test note has a fair point. The abstract's statement that streak secondary instabilities 'lead to' the hairpins and near-wall streaks is inferred from instantaneous snapshots and the symmetry of u′ contours at a single station, not from a secondary-instability analysis. It is a plausible mechanism, but calling it LSI is a stretch without showing that the measured streak profile is linearly unstable at those scales. The delay claim is also not quantified: no transition-onset Reynolds number based on Cf or amplitude is given for Tw095 versus Tw110 versus TadIG, and Fig. 3 uses different Q thresholds and x-ranges, so part of the 'significantly delayed' impression could be visual. These are fixable in revision—add a quantitative onset metric and soften the causal language or add a stability analysis.\n\nAlso minor: there is no grid-convergence study, though the reported Δx+ values are reasonable, and the idealized van der Waals/Jossi-Stiel-Thodos transport model limits quantitative transfer to real fluids like CO2. Neither undermines the qualitative findings.\n\nOverall, the core observation—transcritical K-type is delayed and gradual, unlike H-type—is supported by the amplitude data and visualizations. The gaps are real but not load-bearing. I would send this out for peer review with a request for quantification of the onset, a more cautious mechanism claim, and ideally a grid-convergence check. If I worked in non-ideal-fluid transition, I would cite it. Worth bringing to a reading group focused on transition or supercritical flows.","headline":"First DNS of K-type transition in a supercritical boundary layer, with a plausible but not fully proven mechanism and an unquantified delay claim; deserves serious review with revisions.","tokens_in":5898,"tokens_out":2058,"would_cite":true,"duration_ms":20129,"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":"Heated supercritical wall delays K-type transition and makes the breakdown gradual.","keywords":["K-type transition","boundary-layer transition","supercritical fluids","pseudo-boiling","direct numerical simulation","localized secondary instability","mode II instability","flat-plate boundary layer"],"falsifier":"Run the same K-type forcing with a more accurate real-fluid equation of state for carbon dioxide, or in an experiment, at a reduced pressure of 1.10 and a wall-to-free-stream temperature ratio of 1.222, and locate where the disturbance amplitudes saturate and where the skin-friction maximum occurs; if transition is not delayed relative to the subcritical case, or if a violent breakdown appears at the Lambda-vortex valleys, the central claim is refuted.","tokens_in":4881,"feed_emoji":"🌊","tokens_out":7312,"duration_ms":56031,"temperature":0.7,"pith_summary":"This paper uses direct numerical simulation to compare controlled K-type boundary-layer transition for a supercritical fluid at a reduced pressure of 1.10 under two wall-temperature conditions: one liquid-like and one above the pseudo-boiling (Widom) temperature. It finds that when the wall is heated into the vapour-like regime, the aligned Lambda-vortices of the K-type route do not break down violently; instead, transition is delayed and proceeds through streak secondary instabilities near the legs of the primary vortices, producing extra hairpin vortices and near-wall streaks. The paper argues that this transcritical K-type path is qualitatively different from the violent H-type breakdown reported for the same flow conditions, and also different from the subcritical case, where steady longitudinal modes dominate and transition is only slightly delayed relative to the ideal-gas reference. If correct, the result matters because supercritical-fluid systems such as heat exchangers and rocket engines may experience transition later and more gradually than ideal-gas correlations predict.","feed_headline":"Supercritical wall heat delays and softens K-type transition","feed_subtitle":"In DNS, wall temperature above the pseudo-boiling line postpones turbulence and turns breakdown into streak-driven vortices.","key_machinery":"The central mechanism is the K-type (fundamental) breakdown, in which a two-dimensional fundamental wave and an oblique three-dimensional wave pair produce streamwise-aligned Lambda-vortices whose nonlinear growth leads to turbulence. In the transcritical flow, the load-bearing new element is the localized secondary instability of low-speed streaks: once the steady longitudinal mode (0,1) exceeds roughly 10% amplitude, co-rotating and counter-rotating vortex legs create low-speed streaks that destabilize in a symmetric varicose mode, making hairpin vortices, and an antisymmetric sinuous mode, making the streaks meander and shed near-wall vortices. The paper also relies on mode-II instability, a supercritical-flow instability tied to the kinematic-viscosity minimum at the Widom line, which is the only unstable mode for the fundamental frequency within the DNS domain.","core_discovery":"At supercritical pressure with a wall temperature above the pseudo-boiling temperature, K-type breakdown of a flat-plate boundary layer follows a gradual route: large-amplitude higher harmonics dominate the early breakdown stage, the fundamental resonance is delayed and weakened, and just before transition the steady longitudinal mode (0,1) triggers a localized secondary instability on low-speed streaks near the legs of the aligned Lambda-vortices. The symmetric varicose component evolves into hairpin vortices while the antisymmetric sinuous component causes the streak pair to meander and generates additional near-wall vortices. As a result, transition to turbulence is not violent and occurs significantly later than in the subcritical heating case, and it lacks the strong secondary vortices at the spanwise valleys that characterise the H-type route.","pith_inferences":["The paper leaves implicit that the delaying effect depends on the wall temperature crossing the pseudo-boiling line; systematically varying the reduced wall temperature between 0.95 and 1.10 would reveal how the (0,1) streak amplitude threshold and the transition-onset shift scale with stratification.","By analogy with the H-type comparison, one might expect that the noise-receptivity scenario determines whether supercritical boundary layers transition gradually or violently, so real flows with broad-spectrum disturbances may show a mix of both routes.","A practical corollary is that transition-prediction correlations calibrated on ideal-gas or subcritical DNS may misestimate heat-transfer enhancement in supercritical heat exchangers, since the transcritical K-type route places the turbulent region farther downstream."],"forward_implications":["In the transcritical regime the fundamental (K-type) resonance is delayed and weakened, and higher harmonics such as (2,0) can dominate the early breakdown stage.","Subharmonic resonance between mode (2,0) and the oblique pair (1,±1) becomes active just before transition onset, so the route loses its purely fundamental character.","Steady longitudinal streaks can reach amplitudes above 10% and trigger localised secondary instability, making streak dynamics a necessary ingredient in transcritical transition prediction.","Because the transcritical K-type breakdown is gradual whereas the H-type under the same conditions is violent, the disturbance scenario determines whether supercritical boundary layers transition abruptly or softly."],"supporting_citations":[{"why":"Supplies the ideal-gas K-type DNS baseline against which the delayed transition is measured.","marker":"[9]"},{"why":"Provides the controlled K-type transition setup with aligned Lambda-vortices that the paper reproduces and extends to supercritical fluids.","marker":"[7]"},{"why":"Is the companion H-type DNS under identical flow conditions whose violent breakdown is the direct comparison for the gradual transcritical route.","marker":"[6]"},{"why":"Identifies the mode-II instability tied to the kinematic-viscosity minimum at the Widom line, the instability that drives the fundamental disturbance here.","marker":"[2]"},{"why":"Documents the high-order GPU-accelerated finite-difference solver used for the DNS.","marker":"[8]"}],"fun_headline_variants":["Pseudo-boiling wall delays and softens K-type transition","Supercritical heat above pseudo-boiling softens K-type transition","K-type transition delayed and softened by supercritical wall heat","Streak-driven vortices soften supercritical K-type breakdown","Wall above pseudo-boiling makes K-type transition streak-driven"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"All supercritical simulations use the reduced van der Waals equation of state and analytical transport correlations, so the claim that transcritical K-type breakdown is delayed and gradual depends on those models reproducing the real property variations of a supercritical fluid near the pseudo-boiling line.","fun_headline_variants_meta":{"raw":{"variants":["Pseudo-boiling wall delays and softens K-type transition","Supercritical heat above pseudo-boiling softens K-type transition","K-type transition delayed and softened by supercritical wall heat","Streak-driven vortices soften supercritical K-type breakdown","Wall above pseudo-boiling makes K-type transition streak-driven"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001242,"raw_usage":{"total_tokens":5073,"prompt_tokens":895,"completion_tokens":4178,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":511,"completion_tokens_details":{"reasoning_tokens":4095}},"tokens_in":511,"tokens_out":4178,"duration_ms":27373,"temperature":1.0,"reasoning_tokens":4095,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T15:20:12.856990+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Run the same K-type forcing with a more accurate real-fluid equation of state for carbon dioxide, or in an experiment, at a reduced pressure of 1.10 and a wall-to-free-stream temperature ratio of 1.222, and locate where the disturbance amplitudes saturate and where the skin-friction maximum occurs; if transition is not delayed relative to the subcritical case, or if a violent breakdown appears at the Lambda-vortex valleys, the central claim is refuted.","supporting_citations":[{"cited_title":"and Moin, P.: Direct numerical simulation of complete H-type and K-type transitions with implications for the dynamics of turbulent boundary layers","cited_arxiv_id":null,"evidence_quote":"Supplies the ideal-gas K-type DNS baseline against which the delayed transition is measured."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the controlled K-type transition setup with aligned Lambda-vortices that the paper reproduces and extends to supercritical fluids."},{"cited_title":"C., Gaspar, R., Bugeat, B., Costa, P., Peeters, J","cited_arxiv_id":null,"evidence_quote":"Is the companion H-type DNS under identical flow conditions whose violent breakdown is the direct comparison for the gradual transcritical route."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Identifies the mode-II instability tied to the kinematic-viscosity minimum at the Widom line, the instability that drives the fundamental disturbance here."},{"cited_title":"C., Hirai, R., Costa, P., Peeters J","cited_arxiv_id":null,"evidence_quote":"Documents the high-order GPU-accelerated finite-difference solver used for the DNS."}],"review_version":1}