{"id":"eaa0be31-c25a-4cf4-a18a-749bf729e76e","arxiv_id":"2506.02469","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"A forward quasi-1D stagnation line model using measured free-jet temperature and velocity reproduces measured cold-wall heat fluxes in the VKI Plasmatron for a range of conditions, supporting LTE-based characterization.","lead":"This paper combines optical emission spectroscopy, heat flux and Pitot probes, and simulations to characterize the plasma jet in the VKI Plasmatron wind tunnel. It proposes a forward procedure using measured temperature and velocity to compute stagnation line heat flux, agreeing with experiments within a range of conditions and literature catalytic recombination values.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Agreement hinges on the unquantified CF-ICP-derived inlet location δ*u(us); a sensitivity study is needed before the 'measured data only' claim is secure.","rationale":"The paper's central claim is that a forward procedure, using only measured free-jet temperature and dynamic pressure, reproduces intrusive cold-wall heat flux measurements. The weakest point in that argument is the auxiliary input δ*u(us), which is taken from CF-ICP, a model the paper itself shows to be velocity-biased (Sec. 6.2). The authors acknowledge this (Sec. 7.1) but do not quantify how strongly the predicted heat-flux envelope depends on δ*u. Appendix A is a valuable consistency check, but it only shows STAGLINE can reproduce CF-ICP, not that the experimental agreement survives a biased correlation. A perturbative sensitivity study is straightforward and would settle the question. The low-enthalpy failure at 100 mbar (hs < 20 MJ/kg) is another limitation, but the central claim is already qualified as valid 'through a range of test conditions,' so it does not by itself overturn the claim. I concur with the reader's CONDITIONAL verdict: the procedure is plausible and well-supported by internal consistency, but the unquantified δ*u sensitivity keeps the strongest claim from being fully established. Thus the verdict should remain UNCHANGED (CONDITIONAL).","tokens_in":27884,"tokens_out":5881,"duration_ms":55909,"concrete_test":"Rerun the STAGLINE forward procedure for the conditions in Tables B.2 and B.3 (at least the endpoints: 50 mbar/150 kW, 50 mbar/300 kW, 100 mbar/150 kW, 100 mbar/290 kW) with δ*u perturbed by −50%, −20%, +20%, +50% (and, if feasible, with an alternative δ*u obtained from a velocity-bias-corrected CF-ICP correlation), keeping Ts, us and pc fixed. If the computed cold-wall heat-flux envelope for 0.01 ≤ γ_ref ≤ 0.1 fails to contain the measured ˙qcw for any perturbed case, the agreement is conditional on the CF-ICP correlation; if the envelope remains essentially unchanged, the concern is resolved.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The forward procedure in Sec. 7 is presented as using only measured free-jet temperature and dynamic pressure, but it also requires the inlet-point correlation δ*u(us) from CF-ICP (Fig. 12c). Section 6.2 shows CF-ICP underpredicts axial velocity for a given enthalpy, and Sec. 7.1 explicitly admits δ*u may be mispredicted. Appendix A verifies only that STAGLINE reproduces CF-ICP (temperature within ±1%, velocity gradient within ±10%), not that the procedure is insensitive to the bias in δ*u. Since STAGLINE is initialized with Ts and us at δ*u, a biased δ*u shifts the upstream boundary condition along the jet and changes the inviscid deceleration length and, potentially, the BL-edge velocity gradient and cold-wall heat flux. The paper never quantifies this sensitivity. The reported agreement at 100 mbar is limited to hs > 20 MJ/kg, and the attributed cause (jet decay) is the same regime where δ*u is largest and most uncertain, so the two effects are confounded. Without a δ*u sensitivity study, the central claim that the forward method is consistent with intrusive heat flux measurements across test conditions is not fully established.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper presents a combined experimental and numerical characterization of the subsonic air plasma jet in the VKI Plasmatron at 50 and 100 mbar and 150-300 kW. Spatially resolved OES yields radial temperature profiles from atomic lines, Boltzmann plots, and spectral fits, together with electron densities from H-beta Stark broadening; the diagnostics are mutually consistent below about 7000 K, supporting LTE and allowing inference of the free-jet enthalpy. The OES data are compared with CF-ICP simulations, whose power efficiency is tuned to match the measured temperature, and with the traditional inverse heat-flux rebuilding procedure, exposing discrepancies in velocity and enthalpy. A forward procedure is then proposed: measured free-jet temperature and dynamic pressure are imposed as inlet conditions for the quasi-1D STAGLINE solver, with the inlet distance delta*u taken from a CF-ICP correlation. The computed cold-wall heat-flux envelopes agree with intrusive measurements for enthalpies above about 20 MJ/kg at 100 mbar and over the full tested range at 50 mbar, for recombination coefficients between 0.01 and 0.1.","tokens_in":27931,"tokens_out":8284,"duration_ms":73627,"significance":"If the forward procedure is robust, it provides a valuable alternative to inverse enthalpy rebuilding for ICP ground testing, reducing reliance on assumed catalytic efficiencies and yielding boundary-layer edge velocity gradients useful for LHTS scaling. The OES dataset is a genuine strength: the calibration chain, Abel inversion, and internal consistency checks (line-to-line agreement, oxygen Boltzmann plots, H-beta electron densities) are carefully documented. The main risk is that the heat-flux agreement is demonstrated with one particular CF-ICP-derived inlet-location correlation, and the paper does not quantify the sensitivity of the result to that correlation; the 'measured data only' claim is therefore not yet fully secured. A sensitivity study is within the scope of the manuscript and should be required before publication.","major_comments":[{"comment":"The claim that the forward procedure uses only measured free-jet temperature and dynamic pressure is not supported as written, because the inlet location delta*u(us) comes from a CF-ICP correlation. Section 6.2 reports that CF-ICP underpredicts the axial velocity for a given free-jet enthalpy, and Section 7.1 explicitly concedes that delta*u may be mispredicted; a biased delta*u shifts the STAGLINE inlet along the jet, altering the inviscid deceleration length and potentially the BL-edge velocity gradient and cold-wall heat flux. Appendix A verifies only that STAGLINE reproduces CF-ICP at one condition (50 mbar, 100 kW), which is not a sensitivity test of the delta*u choice. Please add a sensitivity study varying delta*u by an amount consistent with the CF-ICP velocity mismatch and the dynamic-pressure uncertainty, and report the resulting change in the Fig. 13 heat-flux envelope.","section":"Section 7.1 and Fig. 12(c)"},{"comment":"The claimed agreement at 100 mbar is restricted to hs > 20 MJ/kg, and the paper attributes the failure below this threshold to jet-decay effects not represented in the quasi-1D formulation. This is the same regime in which Fig. 12(c) shows delta*u to be largest (low us), so the limitation and the delta*u uncertainty are confounded; without a quantitative bound on how much the heat-flux envelope shifts when delta*u is varied, the conclusion that the procedure agrees with intrusive heat-flux measurements across the reported range is not fully established.","section":"Section 7.2 and Fig. 13"},{"comment":"The +/-5 Pa uncertainty in pdyn is a large relative error at the low end of the data (e.g., pdyn = 15 Pa at 100 mbar, 148 kW in Table B.3), yet the paper does not propagate this uncertainty through Eq. (3) into us and hence into delta*u and the heat-flux envelope. The sensitivity study should include this propagation, as it directly affects the inlet-location correlation and the low-enthalpy regime where agreement is already limited.","section":"Section 3.2 and Eq. (3)"},{"comment":"The delta*u(us) correlation is generated with CF-ICP simulations whose numerical power efficiency is tuned to match the OES free-jet temperature; because the same OES temperature defines the inferred enthalpy hs in the forward procedure, the heat-flux comparison in Fig. 13 is not a fully independent test of the measured-only method. This indirect circularity should be addressed, for example by recomputing delta*u with a fixed power efficiency or by measuring the location experimentally, e.g., from Pitot pressure profiles at two axial stations.","section":"Section 6.1 and Section 7.1"}],"minor_comments":[{"comment":"The phrase 'region region' is duplicated in the sentence describing the boundary layer; remove the repetition.","section":"Section 4.1"},{"comment":"The word 'affliation' should be 'affiliation'.","section":"Footnote 1"},{"comment":"The words 'Reseach Center' should be 'Research Center'.","section":"Section 3.4"},{"comment":"The word 'multidimentional' should be 'multidimensional'.","section":"Appendix A"},{"comment":"The figure shows only fitted trend lines without uncertainty bands; adding representative error bars for the experimental heat flux and enthalpy would make the claimed compatibility easier to assess.","section":"Fig. 13"},{"comment":"The sentence 'the numerical power efficiency is not considered here' is confusing because Fig. 9 explicitly uses eta_sim to match the OES temperature; clarify that Fig. 10 compares data and simulation at matched free-jet enthalpy rather than matched input power.","section":"Section 6.2"},{"comment":"The header 'Acc.˚A' appears garbled; the accuracy column should be formatted properly.","section":"Table C.4"}],"recommendation":"major_revision","confidential_remarks":null},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Fairly strong experimental paper. The new OES setup with absolute calibration and spatial resolution gives the community a solid dataset: radial temperature profiles, electron density from H-beta, and enthalpy/heat-flux/dynamic-pressure maps at 50 and 100 mbar. The internal consistency checks are convincing: different atomic lines agree, Boltzmann plots are linear, and electron density matches LTE predictions from the O777 line. The forward STAGLINE procedure is genuinely new. Instead of inverse rebuilding with an assumed catalytic efficiency, it starts from measured Ts and us and computes the stagnation line with quasi-1D Navier-Stokes. The heat-flux comparison is not a fit to the data, and both HS50 and HS30 geometries falling mostly inside the gamma=0.01-0.1 envelope is a meaningful validation.\n\nThe soft spot is exactly the one the stress-test flags. The inlet point delta*u(us) comes from a CF-ICP correlation, and Sec. 6.2 shows CF-ICP underpredicts axial velocity. Sec. 7.1 admits delta*u may be mispredicted. Appendix A only verifies that STAGLINE reproduces CF-ICP when fed CF-ICP inputs, not that the heat flux is insensitive to a biased delta*u. Since a shifted inlet changes the inviscid deceleration length and therefore the boundary-layer-edge velocity gradient, the unquantified sensitivity is a real gap. It is not fatal, but it does make the \"measured data only\" phrasing slightly strong. The low-enthalpy disagreement at 100 mbar (below 20 MJ/kg) is also confounded with this, since that is where delta*u is largest and jet decay effects matter.\n\nOther limitations are stated honestly: LTE is supported experimentally only below 7000 K, and the high-power spectral fits deviate. The paper provides no code or raw data, so reproducibility is moderate. The citation pattern is appropriate; the work builds on the authors' own prior OES paper and the relevant VKI literature, which is not a problem when the results genuinely extend that line.\n\nVerdict: this deserves a serious referee. A sensitivity study on delta*u and ideally a lower-pressure test point would address the main worry. The experimental core is careful, the uncertainties are documented, and the forward procedure is a useful step for plasma wind tunnel characterization. I would bring it to a reading group and would cite it if I worked on aerothermodynamic ground testing.","headline":"Careful OES-based characterization with a genuinely forward stagnation-line procedure, but the CF-ICP-derived inlet-location correlation needs a sensitivity study before the 'measured data only' claim is fully secure.","tokens_in":28651,"tokens_out":2093,"would_cite":true,"duration_ms":22284,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"Two measured free-jet quantities, temperature and dynamic pressure, plus a quasi-one-dimensional stagnation-line solver, reproduce intrusive heat-flux measurements across a range of plasma wind tunnel conditions.","keywords":["plasma wind tunnel","inductively coupled plasma","optical emission spectroscopy","free-jet enthalpy","stagnation line flow","cold-wall heat flux","local thermodynamic equilibrium","catalytic recombination"],"falsifier":"A direct test would be to measure the axial velocity profile along the stagnation line upstream of a probe without relying on the Pitot conversion, for example by laser Doppler or particle imaging velocimetry in the plasma, and compare it with the profile the forward computation produces from the $\\delta_u^*$ correlation; if the true $\\delta_u^*$ at a given free-jet velocity differs substantially from the correlation, the computed heat flux should move outside the experimental envelope. A cheaper version is a sensitivity sweep: re-run the procedure with $\\delta_u^*$ varied by the known velocity uncertainty and check whether the predicted heat-flux range still contains the measured values, especially at 100 mbar and enthalpies below 20 MJ/kg where the paper already reports degraded agreement.","tokens_in":2042,"feed_emoji":"🔥","tokens_out":2786,"duration_ms":89613,"temperature":0.7,"pith_summary":"This paper proposes that the subsonic plasma jet in a 1.2 MW inductively coupled plasma wind tunnel can be characterized without the usual inverse heat-transfer rebuilding or tunable catalytic-efficiency parameters. The authors use spatially resolved ultraviolet-to-near-infrared optical emission spectroscopy to show that the free jet is in local thermodynamic equilibrium below about 7000 K at the probe location, so the flow enthalpy can be read directly from the measured temperature. Combining that temperature with a Pitot-derived dynamic pressure, they impose two boundary conditions on a quasi-one-dimensional Navier-Stokes computation of the stagnation line and find that the predicted cold-wall heat flux brackets the intrusive measurements for both probe radii. The agreement holds for free-jet enthalpies between roughly 20 and 40 MJ/kg and for copper recombination coefficients in the range commonly reported in the literature, offering a framework that is independent of the catalyst assumptions that limit the traditional method.","feed_headline":"Two measurements rebuild the plasma jet and match heat flux","feed_subtitle":"Forward model from measured temperature and Pitot pressure predicts heat flux without tuning catalyst parameters.","key_machinery":"The load-bearing object is a forward stagnation-line computation built on the quasi-one-dimensional form of the Navier-Stokes equations in conservation form, which needs fewer inputs than the boundary-layer approach: at an inlet point chosen far enough upstream that the radial velocity gradient and its axial derivative are negligible, the only boundary conditions are the free-jet temperature, chamber pressure, and free-jet velocity. The temperature comes from absolute emission spectra after Abel inversion, using the oxygen 777 nm triplet as the reference line, and the velocity comes from the measured dynamic pressure with a low-Reynolds-number correction. The location of the inlet point, $\\delta_u^*$, is taken from a correlation computed with the two-dimensional magnetohydrodynamics solver, and the procedure is verified against that solver, showing the thermal boundary layer is reproduced within about 1 percent and the velocity gradient within about 10 percent. The computation also yields a boundary-layer-edge radial velocity gradient that is nearly insensitive to the wall recombination coefficient, which the paper identifies as the quantity needed for local heat-transfer simulation.","core_discovery":"On its own terms, the central discovery is that the stagnation-line flow and the cold-wall heat flux on a hemispherical copper probe can be computed forward from only two experimentally measured free-jet quantities, the centerline temperature inferred from the oxygen 777 nm emission and the dynamic pressure from a Pitot probe, together with the chamber pressure and probe geometry. The paper argues that the traditional inverse boundary-layer rebuilding method overpredicts the gas enthalpy at 100 mbar for the standard catalytic-efficiency assumptions, while the forward quasi-one-dimensional computation reproduces the measured heat-flux-versus-enthalpy trends at 50 and 100 mbar, for 30 and 50 mm probes, with the data mostly falling inside the 0.01 to 0.1 recombination-coefficient envelope. A supporting observation is that the equilibrium heat flux nearly equals the non-equilibrium fully catalytic heat flux at these conditions, so the earlier disagreement with the two-dimensional magnetohydrodynamics solver is attributed to that solver under-predicting the axial velocity rather than to wall catalysis.","pith_inferences":["If the $\\delta_u^*$ correlation from the two-dimensional solver carries a bias, the imposed inlet location shifts along the jet; a sensitivity study of the predicted heat flux to $\\delta_u^*$ would show how much of the reported agreement depends on that correlation.","The same two-input forward framework could be applied to other test gases and probe geometries, provided the equilibrium assumption remains valid at the chosen free-jet point; this is a direct extension the paper does not demonstrate.","Combining the forward stagnation-line computation with a velocity measurement technique that does not rely on the Pitot correction, for example laser Doppler or particle imaging velocimetry, would provide an independent check of the dynamic-pressure-to-velocity conversion and could extend the method to lower enthalpies where jet decay degraded the agreement.","The paper's restriction to temperatures below 7000 K is tied to the LTE spectral fits, so extending the characterization to higher-power conditions will require modeling the observed non-equilibrium molecular emission rather than a single-temperature fit."],"forward_implications":["Free-jet enthalpy maps can be built from optical emission spectroscopy temperature alone, giving experimental $h_s$ versus $\\dot{q}_{cw}$ and $h_s$ versus $p_{dyn}$ data for benchmarking models.","The strong dependence of the traditional inverse rebuilding procedure on the reference copper recombination coefficient is bypassed; the forward method brackets measurements with the literature range $0.01<\\gamma_{ref}<0.1$.","The boundary-layer-edge radial velocity gradient, a key quantity for local heat-transfer simulation, can be extracted without knowing the wall catalytic efficiency.","The observed under-prediction of dynamic pressure by the two-dimensional magnetohydrodynamics solver implies that improving the velocity prediction in that solver is the main route to extending this characterization away from the stagnation line.","Within the tested envelope the approach works for two probe radii, so it can be used to define consistent test conditions for material response studies."],"supporting_citations":[{"why":"Describes the traditional inverse boundary-layer enthalpy rebuilding procedure that the paper contrasts with the forward method.","marker":"[9]"},{"why":"Supplies the literature range of copper recombination coefficients used to bracket the measured heat fluxes.","marker":"[11]"},{"why":"Provides the reference catalytic efficiencies traditionally assumed at 50 and 100 mbar, against which the rebuilt enthalpies are compared.","marker":"[25]"},{"why":"Describes the two-dimensional magnetohydrodynamics solver whose temperature fields and $\\delta_u^*$ correlation the forward procedure relies on.","marker":"[62]"},{"why":"Supplies the stagnation-line Navier-Stokes solver used for the forward heat-flux computations.","marker":"[68]"},{"why":"Gives the conservation-form stagnation-line equations implemented in the solver.","marker":"[69]"},{"why":"Establishes the conditions $\\beta\\approx0$ and $\\beta'\\approx0$ and the local heat-transfer simulation context used to reduce the boundary conditions to temperature, pressure, and velocity.","marker":"[6]"},{"why":"Supplies thermodynamic, transport, and chemical-equilibrium properties for the gas mixture.","marker":"[55]"},{"why":"Provides Stark broadening tables used to infer electron density from the hydrogen-beta line, supporting the LTE diagnosis.","marker":"[58]"},{"why":"Provides the atomic line data used to extract the O777 reference temperature.","marker":"[54]"}],"fun_headline_variants":["Forward model from T and Pitot predicts heat flux","Two measurements rebuild plasma jet, match heat flux","Plasma jet forward modeling: no catalyst tuning","Heat flux from temperature and dynamic pressure only","Improved ICP jet model: forward method beats inverse"],"cache_read_input_tokens":30592,"weakest_assumption_plain":"The forward procedure assumes a distance from the probe, $\\delta_u^*$, where the measured free-jet conditions are imposed, and that distance is taken from a numerical solver that the paper itself finds under-predicts the jet velocity; if this distance is biased, the whole computation shifts and the reported agreement with heat-flux measurements could be lost.","fun_headline_variants_meta":{"raw":{"variants":["Forward model from T and Pitot predicts heat flux","Two measurements rebuild plasma jet, match heat flux","Plasma jet forward modeling: no catalyst tuning","Heat flux from temperature and dynamic pressure only","Improved ICP jet model: forward method beats inverse"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.00023,"raw_usage":{"total_tokens":1504,"prompt_tokens":988,"completion_tokens":516,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":604,"completion_tokens_details":{"reasoning_tokens":445}},"tokens_in":604,"tokens_out":516,"duration_ms":5324,"temperature":1.0,"reasoning_tokens":445,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-07T11:22:49.771473+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A direct test would be to measure the axial velocity profile along the stagnation line upstream of a probe without relying on the Pitot conversion, for example by laser Doppler or particle imaging velocimetry in the plasma, and compare it with the profile the forward computation produces from the $\\delta_u^*$ correlation; if the true $\\delta_u^*$ at a given free-jet velocity differs substantially from the correlation, the computed heat flux should move outside the experimental envelope. A cheaper version is a sensitivity sweep: re-run the procedure with $\\delta_u^*$ varied by the known velocity uncertainty and check whether the predicted heat-flux range still contains the measured values, especially at 100 mbar and enthalpies below 20 MJ/kg where the paper already reports degraded agreement.","supporting_citations":[{"cited_title":"Degrez, P","cited_arxiv_id":null,"evidence_quote":"Describes the traditional inverse boundary-layer enthalpy rebuilding procedure that the paper contrasts with the forward method."},{"cited_title":"Viladegut, O","cited_arxiv_id":null,"evidence_quote":"Supplies the literature range of copper recombination coefficients used to bracket the measured heat fluxes."},{"cited_title":"Panerai, Aerothermochemistry Characterization of Thermal Protection Systems, Ph.D","cited_arxiv_id":null,"evidence_quote":"Provides the reference catalytic efficiencies traditionally assumed at 50 and 100 mbar, against which the rebuilt enthalpies are compared."},{"cited_title":"Degrez, D","cited_arxiv_id":null,"evidence_quote":"Describes the two-dimensional magnetohydrodynamics solver whose temperature fields and $\\delta_u^*$ correlation the forward procedure relies on."},{"cited_title":"Munafo, Multi-Scale Models and Computational Methods for Aerother- modynamics, Ph.D","cited_arxiv_id":null,"evidence_quote":"Supplies the stagnation-line Navier-Stokes solver used for the forward heat-flux computations."},{"cited_title":"Klomfass, S","cited_arxiv_id":null,"evidence_quote":"Gives the conservation-form stagnation-line equations implemented in the solver."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides Stark broadening tables used to infer electron density from the hydrogen-beta line, supporting the LTE diagnosis."}],"review_version":1}