{"id":"344bb842-1d83-49e8-ba89-cd2e52521a6d","arxiv_id":"2606.17596","paper_version":1,"verdict":"UNVERDICTED","confidence":"LOW","novelty_score":5.0,"correctness_risk":"unknown","formal_verification":"none","parameter_count":0,"one_line_summary":"A posteriori evaluation of twelve AMD-based subgrid models in T-LES identifies a two-layered AMD/Gradient mixed model as providing the best agreement with DNS for solar receiver conditions.","lead":"The study compares twelve subgrid-scale models for thermal large eddy simulation of turbulent channel flows mimicking solar receivers against direct numerical simulation data. A two-layered mixed model using AMD and Gradient approaches shows the best performance.","discovery_kind":"new_application","skeptic_critique":{"model":"grok-4.3","headline":"DNS accuracy as ground truth under solar-receiver conditions is the load-bearing assumption","rationale":"The reader's weakest_assumption directly identifies the single point on which the entire a-posteriori demonstration rests. No other internal inconsistency (e.g., in model formulation or error metric) is visible from the given material, and the paper's own checks on mesh/numerics do not address the reference data quality.","tokens_in":1581,"tokens_out":279,"duration_ms":21595,"concrete_test":"Recompute the global error rate after replacing the reference DNS with a 2× finer mesh run (or published higher-fidelity data at identical parameters); if the ranking of the four selected models changes or the mixed-model advantage drops below the reported margin, the original claim is sensitive to DNS quality.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim (mixed AMD/Gradient model yields best agreement) is established entirely via a posteriori comparison to DNS. For this ranking to be reliable, the DNS must itself be a converged reference solution at the target Re, Pr, and temperature ratios. The abstract provides no information on DNS grid resolution, temporal convergence, or validation against known channel-flow statistics; any under-resolution in the thermal boundary layer would systematically favor or penalize models that alter near-wall dissipation or scalar transport.","agreement_with_reader":"agree"},"referee_report":{"model":"grok-4.3","summary":"The manuscript presents an a posteriori evaluation of 12 subgrid-scale models for Thermal Large Eddy Simulation (T-LES) of anisothermal turbulent channel flows representative of solar receiver conditions. All models are variants of the Anisotropic Minimum Dissipation (AMD) model. The authors compute a global error rate to rank the models against Direct Numerical Simulation (DNS) data and select four for detailed analysis on mesh resolution, numerical schemes, and formulations. They conclude that a two-layered mixed model combining AMD/AMD-scalar with the Gradient model provides the best agreement with DNS.","tokens_in":1678,"tokens_out":418,"duration_ms":18929,"significance":"If the DNS benchmark is a converged reference solution, the work supplies concrete guidance on SGS model selection for T-LES under high temperature-ratio conditions relevant to concentrated solar power. The systematic ranking of AMD-based variants and the identification of a specific mixed formulation constitute a falsifiable contribution to the literature on thermal turbulence modeling.","major_comments":[{"comment":"Abstract: The central claim that the two-layered mixed AMD/Gradient model yields the best agreement is established entirely via a posteriori comparison to DNS. The abstract supplies no information on DNS grid resolution, temporal convergence, or validation against known channel-flow statistics at the target Re, Pr, and temperature ratios. This assumption is load-bearing for the model ranking and must be documented with quantitative evidence.","section":"Abstract"},{"comment":"Abstract: The definition and normalization of the 'global error rate' used to evaluate and select among the 12 models are not stated. Without this, the quantitative basis for declaring one model superior cannot be assessed or reproduced.","section":"Abstract"}],"minor_comments":[{"comment":"The abstract refers to 'solar receivers operating conditions' without listing the specific non-dimensional parameters (Re, Pr, temperature ratio) employed in the simulations.","section":null}],"recommendation":"major_revision","confidential_remarks":null},"author_rebuttal":{"model":"grok-4.3","summary":"We thank the referee for the constructive comments. We address each major comment below and will revise the abstract to improve clarity and self-containment.","responses":[{"response":"We agree that the abstract should include key quantitative details on the DNS to support the model ranking. In the revised manuscript we will add concise information on the DNS grid resolution (number of points and wall-normal spacing), temporal averaging period and convergence criteria, and confirmation that the DNS reproduces established channel-flow statistics at the target Re, Pr and temperature ratio.","revision_made":"yes","referee_comment":"[Abstract] Abstract: The central claim that the two-layered mixed AMD/Gradient model yields the best agreement is established entirely via a posteriori comparison to DNS. The abstract supplies no information on DNS grid resolution, temporal convergence, or validation against known channel-flow statistics at the target Re, Pr, and temperature ratios. This assumption is load-bearing for the model ranking and must be documented with quantitative evidence."},{"response":"The global error rate is defined in the body of the manuscript (Section 3) as a normalized aggregate L2 error across mean and rms profiles of velocity and temperature. To make the abstract self-contained we will insert a brief clause stating its definition and normalization while respecting length constraints.","revision_made":"yes","referee_comment":"[Abstract] Abstract: The definition and normalization of the 'global error rate' used to evaluate and select among the 12 models are not stated. Without this, the quantitative basis for declaring one model superior cannot be assessed or reproduced."}],"tokens_in":1258,"tokens_out":348,"duration_ms":17859,"standing_objections":[]},"desk_editor":{"model":"grok-4.3","letter":"The core result is a head-to-head a posteriori test of twelve AMD variants on low-Mach channel flow at temperature ratios and Reynolds numbers typical of solar receivers. They compute a global error rate, pick the top four, and then check mesh sensitivity and scheme effects. The winner is a two-layer mix of AMD/AMD-scalar with the Gradient model.\n\nWhat the work actually supplies is a concrete ordering for this specific regime. That ordering is new in the sense that prior AMD papers did not target these exact anisothermal solar conditions or run the same twelve-model tournament.\n\nThe soft spot is the DNS itself. The abstract gives no grid size, no check against known channel statistics, and no statement on thermal-boundary-layer resolution. If the DNS under-resolves the near-wall scalar transport, models that change dissipation or scalar flux will be ranked on an imperfect target. The global error metric is also undefined in the abstract, so it is impossible to judge whether the reported advantage is robust or metric-dependent.\n\nThis is a narrow engineering study. Readers who run T-LES for concentrated solar power or similar high-temperature channel flows will find the model ranking useful as a starting point. Readers outside that niche will not.\n\nSend it to review. The comparison is systematic enough to be worth referee time, provided the authors supply the missing DNS convergence data and error definitions in revision.","headline":"The paper ranks twelve AMD-based T-LES models against DNS for anisothermal solar-receiver channel flow and flags one two-layer mixed variant as best, but the ranking rests entirely on the unshown quality of that DNS reference.","tokens_in":2147,"tokens_out":370,"would_cite":false,"duration_ms":11299,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"grok-4.3","headline":"A two-layered mixed subgrid model combining AMD with the Gradient model matches DNS most closely in thermal large-eddy simulations of solar-receiver channel flows.","keywords":["thermal large-eddy simulation","anisothermal channel flow","subgrid-scale modeling","AMD model","solar receiver","a posteriori evaluation","mixed model","low-Mach Navier-Stokes"],"falsifier":"A new direct numerical simulation at the same Reynolds and Prandtl numbers but with a different numerical method or domain size that produces mean profiles or fluctuation statistics differing from the reference DNS used here.","tokens_in":2517,"feed_emoji":"","tokens_out":681,"duration_ms":12213,"temperature":0.7,"pith_summary":"The paper evaluates twelve subgrid-scale models, all derived from the Anisotropic Minimum Dissipation approach, inside low-Mach-number thermal large-eddy simulations of anisothermal turbulent channel flow at conditions typical of solar receivers. It computes a global error rate against existing direct numerical simulation data, then examines four models in greater detail under changes in mesh resolution and numerical scheme. The central finding is that a two-layered mixed formulation that adds the Gradient model on top of AMD for both momentum and scalar transport produces the smallest deviations from the reference DNS across the tested quantities.","feed_headline":"Mixed AMD-Gradient model matches DNS best in solar-receiver T-LES","feed_subtitle":"A two-layer closure combining AMD with the Gradient model reduces error across velocity and temperature statistics under receiver conditions","key_machinery":"The two-layered mixed AMD/AMD-scalar plus Gradient model, which supplies the subgrid stresses and scalar fluxes in the filtered low-Mach Navier-Stokes equations.","core_discovery":"Under solar-receiver operating conditions, a two-layered mixed model that combines the AMD/AMD-scalar formulation with the Gradient model yields the best agreement with DNS data among the twelve AMD-based closures examined.","pith_inferences":["The same mixed closure could be tested in geometries that include the actual absorber tubes or windows of a solar receiver rather than a periodic channel.","If the model remains superior when the flow includes radiation or variable properties, it would reduce the need for finer grids in engineering calculations of receiver efficiency.","Running the identical model set on an independent experimental data set at comparable conditions would confirm whether the ranking holds outside the DNS reference.","Extending the two-layer mixing strategy to other base models besides AMD might produce further gains at higher Reynolds numbers."],"forward_implications":["The mixed model can be adopted directly for production T-LES of solar-receiver geometries at the tested resolution.","Coarser meshes remain usable provided the mixed model is retained, because the global error rate stays lower than for the other closures.","Both convective and diffusive numerical schemes interact with the model choice, so the mixed formulation should be paired with the same schemes used in the evaluation.","The advantage of the mixed model appears in both velocity and temperature statistics, indicating consistent improvement for the coupled thermal problem."],"fun_headline_variants":["Two-layer AMD-Gradient model best matches DNS in solar T-LES","AMD-Gradient hybrid leads twelve AMD models vs DNS in receiver flows","Mixed AMD scalar-Gradient tops T-LES accuracy under solar conditions","Gradient model with AMD-AMD scalar excels in a posteriori solar T-LES"],"cache_read_input_tokens":2112,"weakest_assumption_plain":"The DNS data set provides an accurate and complete ground truth for every quantity compared in the a-posteriori tests.","fun_headline_variants_meta":{"raw":{"variants":["Two-layer AMD-Gradient model best matches DNS in solar T-LES","AMD-Gradient hybrid leads twelve AMD models vs DNS in receiver flows","Mixed AMD scalar-Gradient tops T-LES accuracy under solar conditions","Gradient model with AMD-AMD scalar excels in a posteriori solar T-LES"]},"model":"grok-4.3","cost_usd":0.006136,"raw_usage":{"total_tokens":2827,"prompt_tokens":530,"num_sources_used":0,"completion_tokens":76,"cost_in_usd_ticks":61362000,"prompt_tokens_details":{"text_tokens":530,"audio_tokens":0,"image_tokens":0,"cached_tokens":256},"completion_tokens_details":{"audio_tokens":0,"reasoning_tokens":2221,"accepted_prediction_tokens":0,"rejected_prediction_tokens":0}},"tokens_in":530,"tokens_out":76,"duration_ms":13433,"temperature":1.0,"reasoning_tokens":2221,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-06-26T21:57:44.667104+00:00","model_set":{"reader":"grok-4.3"},"falsifier":"A new direct numerical simulation at the same Reynolds and Prandtl numbers but with a different numerical method or domain size that produces mean profiles or fluctuation statistics differing from the reference DNS used here.","supporting_citations":[],"review_version":1}