{"id":"8424ffdb-e476-4318-9a86-8dbe8c43ff61","arxiv_id":"2512.10623","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":3,"one_line_summary":"A spectral rescaling method inflates low-Reynolds-number velocity slices to high Reynolds numbers, producing turbulent boundary layer inflow whose mean-flow quantities match a precursor simulation right from the inlet.","lead":"This paper introduces a way to generate turbulent inflow for simulations at high Reynolds numbers by rescaling velocity slices from a cheaper, lower-Reynolds-number flow. The method reportedly shortens the streamwise distance required for the flow to become realistic by about an order of magnitude compared with current approaches.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The right-from-inlet cf/H12 agreement is by construction: the target mean profile in Eq. (2.16) is taken from the same reference database, so the headline metric does not test the scaling procedure.","rationale":"The reader's verdict is already CONDITIONAL, and I agree with that overall assessment. The reader's weakest_assumption focuses on the self-similar scaling of outer modes; that is a legitimate concern, especially given the Reθ=790 failure. However, I find a more load-bearing and more easily demonstrated problem: the validation's headline metric (cf and H12 right from the inlet) is determined primarily by the target mean profile, which is supplied verbatim from the reference database. This does not invalidate the method's core idea—the scaled fluctuations and their downstream development are still tested—but it means the paper's strongest quantitative claim ('within ±3.5% and ±0.5% right from the inlet') is not evidence for the scaling procedure. The concrete test of replacing the database mean profile with a universal correlation would cleanly separate the contribution of the mean-profile input from the contribution of the fluctuation-scaling method. Until that test is run, the claim should be presented as conditional on having an accurate target mean profile, which is a substantially weaker statement for a method whose selling point is avoiding expensive precursor simulations. Since the paper already acknowledges the Reθ=790 limitation, the 'regardless of base Re' phrasing is also overstated, but this is secondary to the circularity concern. The verdict remains CONDITIONAL: the method shows promise, but the central quantitative claim needs revalidation with an independent mean profile.","tokens_in":17793,"tokens_out":6819,"duration_ms":78217,"concrete_test":"Repeat the inflow DNS for re2k_sc and re4k_sc exactly as in §3.2, but replace U_i,target(y) in Eq. (2.16) with a universal composite mean profile (e.g., Monkewitz, Chauhan & Nagib 2007 or Coles–Fernholz) at Reθ=8000 instead of the Eitel-Amor et al. (2014) database profile. Keep all scaling operations and parameters identical. If cf and H12 at the inlet (x=0) remain within ±3.5% and ±0.5% of the re8k reference, the claim stands; if they deviate beyond these bounds, the headline agreement is an artifact of supplying the exact reference mean profile, and the order-of-magnitude development-length claim must be revised to the case where an accurate a priori mean profile is available.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim is that cf and H12 are within ±3.5% and ±0.5% of a precursor simulation right from the inlet, supporting an order-of-magnitude reduction in development length. However, §2 makes clear that the method requires as input the mean velocity profile at the target Re, and the validation uses the exact mean profile from the Eitel-Amor et al. (2014) database: 'since the database contained the information at the target Re, we use the mean velocity profiles as well from the database.' Equation (2.16) directly adds this target mean profile to the scaled fluctuations. cf is the wall shear (determined by the imposed near-wall mean profile) and H12 is an integral of the mean profile. Therefore the inlet cf and H12 agreement is essentially guaranteed by the input, not produced by the scaling procedure. The actual method's contribution is the fluctuating field, whose recovery after ~8δ99 is the real evidence; but the 'right from the inlet' phrasing and the order-of-magnitude gain attributed to it are confounded. A trivial method that only imposed the reference mean profile would achieve the same cf/H12 metrics. The abstract also claims 'regardless of the base Re tested,' yet §4 admits base Reθ=790 required hand-tuning of space and energy scaling parameters—so the claim is overstated for at least one tested case. The method may still work, but the validation as presented does not substantiate the headline quantitative claim.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript proposes a spectral-scaling inflow-generation method for spatially developing turbulent boundary layers at high Reynolds number. Starting from two-dimensional cross-stream velocity slices at a lower base Re, the method Fourier-decomposes the fields in span, separates inner and outer wavenumber bands, applies POD in the wall-normal direction, shifts the outer wavenumbers by the ratio sc = δ99,target/δ99,base, rescales the wall-normal coordinate and time, applies an energy scaling factor, and finally adds the target-Re mean velocity profile (Eq. 2.16). The generated fields are used as Dirichlet inflow in Neko DNS for a TBL from Reθ ≈ 8000 to 9000. The paper reports that cf and H12 agree with the precursor reference within ±3.5% and ±0.5% right from the inlet, that Reynolds stresses recover after about 8δ99_0, and that this yields an order-of-magnitude reduction in development length relative to previous high-Re TBL simulations.","tokens_in":18156,"tokens_out":4894,"duration_ms":54336,"significance":"If fully supported, the method would be a substantial practical contribution: it would allow high-Re TBL DNS to start from a small, low-Re precursor database instead of a long, expensive precursor domain, and the presented computational experiments are large and carefully executed. The spectral and POD-based scaling description is clear, the algorithm is explicitly formulated, and the downstream development of Reynolds stresses provides genuine evidence that the scaled outer structures survive and the missing spectral gap fills in. However, the headline inlet cf/H12 agreement is largely a consequence of imposing the target-Re mean profile taken from the same database used as reference, and the 'regardless of base Re' claim is not supported by the evidence: only two base Re values work without adjustment, while the Reθ=790 case required hand-tuning. The core idea is promising but the validation as presented does not substantiate the strongest quantitative claims.","major_comments":[{"comment":"The 'right from the inlet' cf and H12 agreement is not a test of the scaling procedure. Equation (2.16) explicitly adds Ui,target(y), and §2 states that, because the database contained the information at the target Re, the mean velocity profiles are taken from the same database. Since cf is determined primarily by the near-wall mean profile and H12 by an integral of the mean profile, the inlet agreement is imposed by construction. The genuine validation of the method is the downstream recovery of Reynolds stresses and spectra. The abstract should be revised so that the inlet cf/H12 statement is attributed to the prescribed mean profile, and a test using a target mean profile obtained from an external correlation (rather than from the reference database) should be reported or at least discussed.","section":"Abstract; §2 and §2.3, Eq. (2.16)"},{"comment":"The claim 'regardless of the base Re tested' is stronger than the evidence. Only Reθ=2240 and Reθ=4430 scaled to Reθ=8000 are presented as working seamlessly. The Reθ=790 case, discussed in §4, required 'a couple of iterations of trial and error' to find the correct scaling parameters for space and energy, and no DNS inflow test is shown for that case. The base-Re independence claim should therefore be restricted to the range of base Re for which outer scaling is valid, and the 790 case should be presented as a limitation rather than as confirming the general claim.","section":"Abstract and §4"},{"comment":"The method depends on several user-selected parameters whose sensitivity is not assessed: the inner/outer wavenumber cutoffs (30 ≤ λz+ ≤ 500 and the Re-dependent outer band), the POD truncation nmodes = 40, and the energy scaling factor Esc obtained from the Alfredsson et al. correlation (Eq. 2.15) and then applied isotropically to all three velocity components. The paper should either provide a sensitivity analysis or explicitly state these as part of the method's input. Without such an analysis, the claimed robustness of the method cannot be separated from the hand-tuning that was needed at low base Re.","section":"§2.1, §2.2.1, §2.2.5"},{"comment":"The 'order-of-magnitude reduction in development length' is based on comparing the development of Reynolds stresses up to about 8δ99_0 with the development length of 3–4δ99U∞+ reported by Sillero et al. for H12 to approach an empirical fit. These are different diagnostics, and the criterion for 'development length' is not stated precisely. The comparison should use the same quantitative criterion (for example, H12 within a specified tolerance of a correlation, or the same Reynolds-stress measure) and should report the result for the current method with the same non-dimensionalization. This would make the order-of-magnitude claim directly verifiable.","section":"§3.2.2 and Fig. 8"}],"minor_comments":[{"comment":"Typo: 'unitray normalization' should be 'unitary normalization'.","section":"§2.1"},{"comment":"The notation Kouter, K_outer, and K_outer is confusing. Use a single consistent symbol, e.g., K_outer, and define it once before Eq. (2.8).","section":"§2.2.1"},{"comment":"The text says the square of the scaling factor Esc^2 is found between the base and target intensities, but then Esc is applied to the POD modes in Eq. (2.16). Please clarify whether the factor applied to the modes is Esc or sqrt(Esc^2), and make the convention consistent.","section":"§2.2.5"},{"comment":"The first term is written as Ui,target(y,z) for kz=0, n=0, but the mean profile is a function of y only. Writing Ui,target(y) would be clearer.","section":"Eq. (2.16)"},{"comment":"The caption does not mention that the right panel shows median-filtered data; this is described only in the text. Either add it to the caption or make the distinction clearer in the figure itself.","section":"Fig. 5 caption"}],"recommendation":"major_revision","confidential_remarks":"The paper has a genuinely interesting and potentially useful method, and the DNS investment is substantial. The main concern is not the scaling operation itself but the framing of the validation: the inlet cf/H12 agreement is largely built into the input, and the base-Re independence claim is overstated. If the authors revise the abstract and the conclusions to separate the imposed mean-flow part from the fluctuating-field part, perform at least one test with a mean profile obtained from a correlation, and report a sensitivity study for the user-chosen parameters, the paper could become a strong contribution. As it stands, the strongest quantitative claims outrun the evidence."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Quick take: the paper has a genuinely new idea for generating high-Re TBL inflows — rescale the outer-layer spectral content from a lower-Re database in spanwise wavenumber, wall-normal coordinate, and time — and the DNS evidence for the fluctuating field is genuinely encouraging. But the headline claim that cf and H12 match a precursor 'right from the inlet' is largely by construction, because the target mean profile is taken from the same database and imposed in the reconstruction. The method deserves peer review; the abstract overstates what is actually tested.\n\nWhat's new: the specific rescaling recipe — shifting the outer wavenumbers, stretching the POD modes in y, slowing the time coefficients, boosting energy — is not in the cited recycling, precursor, or synthetic-eddy literature. The re8k_onlyIO control case is a good move: it shows that cutting the same spanwise wavenumbers without the rescaling and energy treatment leaves the flow under-energized, so the scaling is doing real work. The authors are honest about the coarse y-sampling, the energy over-prediction at the inlet, and the Reθ=790 case needing hand-tuning. The DNS itself — 3.2e9 grid points in Neko, three inflow variants plus a reference case — is serious computation.\n\nSoft spots. First, the cf/H12 agreement at the inlet is essentially forced. cf is wall shear and H12 is an integral of the mean profile; both come from Ui,target(y), read straight from the Eitel-Amor database (the paper says so in §2). The 'right from the inlet' numbers don't test the scaling procedure — imposing the reference mean profile alone would give the same result. The real test is the Reynolds-stress recovery after roughly 8δ99, which looks good in the figures but is never quantified with an error metric. Second, 'regardless of the base Re tested' is too strong: two favorable base Re values (2240 and 4430) are shown, and the 790 case required trial-and-error parameter tuning to work at all. Third, the order-of-magnitude development-length claim compares Reynolds-stress recovery here with an H12-based development length in Sillero et al., so it is not apples-to-apples.\n\nNone of this sinks the method. The fluctuation evidence is real, and the issues are fixable: test with a universal mean profile as input, report a quantitative mismatch measure for the Reynolds stresses, and show sensitivity to the free parameters (Esc, nmodes, the wavenumber cutoffs).\n\nThis is a paper for people who run high-Re TBL DNS and want to cut precursor costs. Practical advance, not a conceptual breakthrough. Send it to peer review — a good referee can push for the clean mean-profile test and the sensitivity study, and the work will be stronger for it.","headline":"Genuinely new inflow method with a plausible DNS demonstration, but the headline cf/H12 agreement is by construction (target mean profile imposed from the same database) and the 'regardless of base Re' claim is only supported by two favorable cases.","tokens_in":18620,"tokens_out":6357,"would_cite":true,"duration_ms":61523,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["47.27.E-","47.27.nb"],"model":"deepseek-v4-flash","headline":"A single scaling factor applied to low-Reynolds-number velocity slices generates high-Reynolds-number turbulent inflow conditions with nearly immediate accuracy, cutting the required development length by roughly an order of magnitude.","keywords":["turbulent boundary layer","inflow generation","Reynolds number scaling","direct numerical simulation","pre-multiplied spectra","proper orthogonal decomposition","skin friction coefficient","development length"],"falsifier":"Apply the scaling to a base Re_theta=2240 and a target Re_theta=16000 (well beyond the calibrated range) using only the base slices and a mean profile, without any tuning. If the pre-multiplied spectra at the inlet show the outer peak shifted by more than the spectral resolution from a reference DNS at that same Re, or if the developed skin friction deviates by more than 3.5% from the reference, the single-sc scaling assumption is falsified.","tokens_in":17704,"feed_emoji":"🌀","tokens_out":4982,"duration_ms":48394,"temperature":0.7,"pith_summary":"The paper proposes a method to generate turbulent inflow conditions for high-Reynolds-number (Re) boundary layer simulations by scaling low-Re velocity slices instead of simulating a long development region. The method splits each spanwise-wavenumber slice into inner and outer regions, then stretches the outer-region modes in span, wall-normal distance, and time by a factor equal to the ratio of boundary-layer thicknesses between target and base Re, and adjusts their energy. In direct numerical simulations targeting Re_theta=8000 from base Re_theta=2240 and 4430, skin friction and shape factor match a full precursor simulation within 3.5% and 0.5% right from the inlet, and Reynolds stresses converge after about eight inlet boundary-layer thicknesses. This reduces the required development length by roughly an order of magnitude compared to existing methods, lowering the cost of high-Re simulation.","feed_headline":"Scaling leap cuts high-Re turbulence simulation cost tenfold","feed_subtitle":"Low-Re velocity slices, stretched by one thickness ratio, reproduce high-Re boundary layers right from the inlet.","key_machinery":"The machinery is a spectral-space scaling pipeline. Cross-stream velocity slices from a low-Re precursor are Fourier-transformed in the periodic spanwise direction; wavenumbers with lambda_z+ between about 30 and 500 are classified as inner and left unchanged, while the outer wavenumbers (beyond roughly 500 to 1700 in lambda_z+ depending on base Re) are processed. Singular value decomposition in the wall-normal direction yields POD modes and time coefficients. For each outer mode, the spanwise wavenumber is shifted to round(k_z/sc), the wall-normal axis is stretched by sc and re-interpolated, the time axis is compressed by sc, and the mode amplitude is scaled by an energy factor E_sc. Recons","core_discovery":"The central discovery is that large-scale outer-region structures in a turbulent boundary layer can be transported across Reynolds numbers by a single geometric scaling factor, sc = delta99_target/delta99_base, applied to the spanwise wavenumbers, wall-normal coordinate, and time of the outer-region POD modes extracted from low-Re precursor slices. Because the inner region is nearly invariant in inner units, it is left untouched. The scaled slices, when seeded with a target-Re mean profile and an energy correction, behave like genuine turbulence at the target Re: the flow equilibrates in about 8 delta99_0, an order of magnitude faster than methods that must let large scales develop on their","pith_inferences":["If the single-factor scaling is genuinely Re-independent, then the method should work for significantly higher targets than the tested Re_theta=8000, e.g., 16000 or beyond, using the same base data; a successful demonstration would considerably strengthen the cost-saving claim.","The observed equilibration time being governed by the spectral gap (missing intermediate wavenumbers) points to a specific physical bottleneck: it is the filling of the mid-range scales, not the survival of the large scales, that sets the development length. Seeding those mid-range scales directly might shorten the adaptation further.","The failure at Re_theta=790 is a caution: the outer-scaling assumption is only valid above some threshold Re. A practical version of the method should state that threshold as a measurable condition (e.g., existence of a clear outer peak in the pre-multiplied spectrum), rather than an ad-hoc base-Re value.","Since the inner region is left untouched, the method relies on the near-wall structures regenerating themselves quickly; an extreme test would be to use a base Re so low that the inner and outer spectral regions overlap, which is exactly where the method is expected to break down."],"forward_implications":["The method removes the need for a large precursor domain: a small, fixed-size precursor at a low base Re can serve any higher target Re, since the scaling leap avoids simulating the expensive development of outer structures.","Inflow quality is high immediately: skin friction coefficient and shape factor stay within ±3.5% and ±0.5% of the reference precursor right from the inlet for both base Re tested.","Two-point statistics equilibrate in roughly 8 inlet boundary-layer thicknesses, which is an order of magnitude shorter than the development lengths reported for other high-Re TBL simulations.","The method naturally combines with autoregressive time-series extension of the POD coefficients, so the inflow signal can be made arbitrarily long without resampling.","The same spectral scaling idea should transfer to other wall-bounded flows (pipes, channels) and to flows such as jets and mixing layers, as long as a scaling law for the spectra is known."],"fun_headline_variants":["Scale outer turbulence to jump Re cheaply","Low-Re slices stretch into high-Re turbulence","One scaling factor reproduces high-Re boundary layers","10x faster inflow: stretch low-Re eddies to high-Re","Outer-layer scaling speeds high-Re inflow generation"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"The method assumes that the large-scale outer structures of a low-Re boundary layer can be faithfully rescaled to a higher Re by a single factor sc determined by the boundary-layer thickness ratio; if this self-similarity of the outer region is not exact, the claimed inlet accuracy fails.","fun_headline_variants_meta":{"raw":{"variants":["Scale outer turbulence to jump Re cheaply","Low-Re slices stretch into high-Re turbulence","One scaling factor reproduces high-Re boundary layers","10x faster inflow: stretch low-Re eddies to high-Re","Outer-layer scaling speeds high-Re inflow generation"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001006,"raw_usage":{"total_tokens":4152,"prompt_tokens":868,"completion_tokens":3284,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":612,"completion_tokens_details":{"reasoning_tokens":3211}},"tokens_in":612,"tokens_out":3284,"duration_ms":20766,"temperature":1.0,"reasoning_tokens":3211,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-03T17:04:11.629989+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Apply the scaling to a base Re_theta=2240 and a target Re_theta=16000 (well beyond the calibrated range) using only the base slices and a mean profile, without any tuning. If the pre-multiplied spectra at the inlet show the outer peak shifted by more than the spectral resolution from a reference DNS at that same Re, or if the developed skin friction deviates by more than 3.5% from the reference, the single-sc scaling assumption is falsified.","supporting_citations":[],"review_version":1}