{"id":"f41b63c5-3b13-428a-80c8-21525a1cf7ce","arxiv_id":"2504.15177","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":7,"one_line_summary":"An rp-adaptive implicit shock tracking method resolves viscous shocks and boundary layers on coarse grids by coupling optimization-driven mesh compression with local p-refinement.","lead":"This paper presents a computational method that moves mesh points and raises polynomial order to resolve thin shocks and boundary layers in high-speed viscous flow without artificial stabilization. It reports accurate heat flux predictions on a hypersonic cylinder using fewer degrees of freedom than standard adaptive refinement in that test.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Heat-flux accuracy is asserted from self-convergence and an undisclosed reference rather than an independent converged solution.","rationale":"The reader's conditional verdict is appropriate and I do not recommend changing it. However, I identify a slightly different load-bearing concern than the reader's stated weakest assumption: rather than parameter transferability being the primary risk, the central claim of accurate heat-flux prediction is directly under-evidenced because the cylinder results lack any external converged reference. Parameter sensitivity (lambda, Delta, continuation) is a real generality concern that the paper itself documents, but it would affect the breadth of the method's applicability; the missing reference affects whether the flagship accuracy claim is demonstrated even for the single presented case. The paper deserves credit for its honest ablations, side-by-side h-adaptation comparisons, and the external Blasius check for the flat plate; those are genuine evidence. The concrete test proposed would settle the heat-flux accuracy question directly, and the existing material should be amended with the reference solution and the origin of the error bands before the central claim can be taken as fully established.","tokens_in":30401,"tokens_out":6049,"duration_ms":63652,"concrete_test":"Run the Section 4.2.2 cylinder case with an independent, grid-converged solver and compare the pointwise heat flux profile to Figure 25; also disclose the reference used for the Figure 27 sub-1% error band. If the pointwise deviation exceeds 5% or the integrated value falls outside the quoted band, the accurate heat-flux claim is unsupported.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The abstract's central quantitative claim is that the rp-adaptive HOIST method yields \"accurate prediction of heat flux profiles\" for hypersonic cylinder flow. In Section 4.2.2, Figure 25 shows heat flux profiles at successive p-adaptation iterations, which become smoother and stabilize; however, no converged reference solution is overlaid, so the profile is only shown to stop changing, not to be accurate. Figure 27 reports that both HOIST and the h-adaptive comparator reach \"sub-1% error\" in integrated heat flux, but the reference defining that 1% is not specified; if it is itself produced by the same IPDG/DWR-h-adaptive family, apparent agreement may reflect shared discretization bias rather than independent accuracy. This matters because the headline claim is precisely heat-flux accuracy, and internal p-convergence alone is consistency, not validation. The flat-plate case has an external Blasius reference (Figure 18), which supports that configuration, but the cylinder—the paper's flagship case—has no analogous external check. The lack of released code or data compounds the difficulty of settling this independently.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper proposes an rp-adaptive extension of the high-order implicit shock tracking (HOIST) method for viscous, shock-dominated flows. The method minimizes an enriched IPDG residual simultaneously over flow variables and mesh nodal coordinates, thereby compressing elements into viscous shocks and boundary layers; robustness is achieved through residual weighting, step constraints and modifications, viscosity continuation, and Hessian regularization. A p-adaptivity loop locally increases the polynomial degree using residual-based, dual-weighted, or feature-based indicators. The method is tested on viscous Burgers problems, a laminar flat-plate boundary layer, and hypersonic flow over a cylinder, and is compared with anisotropic h-adaptation in terms of accuracy per degree of freedom.","tokens_in":30590,"tokens_out":5515,"duration_ms":54150,"significance":"If the central claims hold, this is a substantial contribution: it extends implicit shock tracking from inviscid discontinuities to viscous transitions, removes the need for nonlinear stabilization in the final DG solve, and couples aggressive r-adaptivity with p-adaptivity in an optimization-based framework. The paper contains several genuinely independent checks: the Burgers solutions are compared with reference finite volume solutions, the flat-plate boundary-layer profile and skin friction are compared with the Blasius solution, and the stability/accuracy benefit of each proposed component is demonstrated through ablation studies. These elements give credibility to the algorithmic framework. However, the headline quantitative claim — accurate heat flux prediction for hypersonic cylinder flow — is supported mainly by self-convergence of the heat flux profile and by a 'sub-1% error' assertion whose reference solution is not disclosed. That gap, together with the apparent sensitivity of the method to hand-set parameters, prevents the paper from being accepted in its present form.","major_comments":[{"comment":"The headline claim of accurate heat flux prediction is not backed by an independent converged reference solution. Figure 25 shows that successive p-adaptation iterates of the pointwise heat flux become smoother and stabilize, which demonstrates self-consistency but not accuracy relative to the true heat flux. Figure 27 claims that both HOIST and the h-adaptive comparator reach 'sub-1% error' in integrated heat flux, but the reference value defining that 1% is not identified anywhere in the text. If the reference is generated by the same IPDG/DWR-h-adaptive method family, agreement could reflect shared discretization bias rather than independent accuracy. I request that the authors overlay a grid-converged external reference (e.g., an independently computed fine-mesh solution) in Figure 25 and report the integrated heat flux error relative to that reference in Figure 27.","section":"§4.2.2, Figures 25 and 27"},{"comment":"The method's central behavior depends on user-set parameters — enrichment degree Δ=2 and boundary-layer residual scaling λ=10 for the cylinder and λ=100 for the flat plate — and the ablation studies in §4.2.4 and §4.2.5 show that Δ=1 or λ=1 substantially degrades boundary-layer resolution and heat-flux prediction. No criterion or automatic procedure is provided for selecting these values for a new geometry or flow regime, and the values were evidently tuned on the cylinder problem used for the headline heat flux claim. This makes the 'stabilization-free' resolution claim contingent on parameters that may not transfer. I ask for either a priori selection rules or a sensitivity study over Δ and λ on an independent configuration.","section":"§3.5.1, §3.5.3, §4.2.4, §4.2.5"},{"comment":"The p-adaptivity indicator comparison shows that a dual-weighted-residual indicator tailored to the integrated heat flux produces a trial space for which the PTC solver cannot converge without stabilization, while the feature-based indicator yields inaccurate integrated heat flux; the enriched-residual indicator is then declared to combine the best aspects of both. This is a plausible conclusion, but the supporting evidence is mostly qualitative degree distributions. A quantitative comparison of the three indicators in terms of heat-flux profile error against the same external reference requested above would strengthen the claim that the enriched-residual indicator is the right choice for the proposed rp-adaptive framework.","section":"§4.2.6"}],"minor_comments":[{"comment":"There is a typo on page 2: 'whiche amounts to a savings factor' should be 'which amounts to a savings factor.'","section":"§1"},{"comment":"References [53] and [54] are duplicates of the same paper (Zahr, Shi, and Persson, J. Comput. Phys. 410:109385, 2020); one should be removed and the citations renumbered. Also, reference [9] lists 'AIAA Paper 2025-153' despite the 2009 conference date; the paper number appears incorrect.","section":"References"},{"comment":"In the degree-of-freedom comparison, flow and mesh degrees of freedom are counted for HOIST while only flow degrees of freedom are counted for the h-adaptive method. This is conservative for the HOIST efficiency claim, but the asymmetry should be stated explicitly in the caption or text so readers do not misread the comparison as apples-to-apples.","section":"§4.2.1, Figure 19"},{"comment":"The initialization for the cylinder problem uses a shock-capturing DG solution with PDE-based artificial viscosity, yet the abstract claims the method operates 'without nonlinear stabilization.' This is not inconsistent — the final solution phase is stabilization-free — but the sentence should be qualified to clarify that an artificial-viscosity solution is used only to initialize the optimization.","section":"§4.2.2"},{"comment":"The notation in (48)–(50) overloads C: it denotes both the scalar function of the residual logarithm and the range parameter in the pair (Cl, Cu). Using a different symbol for the log-residual field would avoid ambiguity.","section":"§3.5.6"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is technically interesting and the Burgers and flat-plate validations are credible. The main obstacle is the missing external reference for the cylinder heat-flux claim, which is the paper's flagship result; this is fixable with additional computational experiments. The parameter-sensitivity issue is also fixable with a robustness or sensitivity study."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Here's my read. The paper is a credible extension of HOIST from inviscid shocks to viscous shock-dominated flows. What's actually new: a primal IPDG form of implicit shock tracking that avoids the auxiliary variables of viscous MDG-ICE, aggressive r-adaptation that compresses elements into viscous shocks and boundary layers, residual weighting so weak boundary layers get equal priority in the objective, viscosity continuation to avoid carbuncles, and p-adaptivity driven by the enriched residual. The Burgers tests are checked against refined finite volume solutions, the flat plate against Blasius, and the authors honestly report that the flat plate case is not competitive with h-adaptation. That earns trust.\n\nThe main soft spot is the flagship claim of accurate heat flux prediction for the hypersonic cylinder. Figure 25 shows the heat flux profile stabilizing and losing inter-element jumps as p increases, but no external converged reference is overlaid. Figure 27's 'sub-1% error' line is not tied to an independent reference, so if both HOIST and the h-adaptive comparator use the same IPDG/DWR family, apparent agreement could share discretization bias. Internal p-convergence is consistency, not validation. This matters because the abstract's headline is precisely heat-flux accuracy. The flat plate Blasius check is a genuine external anchor for the method's boundary layer behavior, so the circularity is partial, not total, but the cylinder case still needs an independent reference.\n\nThe method also leans on several hand-set parameters (lambda, Delta, continuation schedule, Hessian regularization constants), and the sensitivity studies for lambda and Delta run on the same cylinder problem used for the headline claim. The paper is honest that changing them degrades the solution (carbuncle without continuation, poor boundary layers with Delta=1 or lambda=1), which confirms the tuned character. No code or data is released, so independent reproduction is harder than it should be. I don't see a load-bearing flaw in the mathematics or the optimization formulation; the claims are scoped carefully except for the cylinder reference issue.\n\nThis paper deserves a serious referee. The useful outcome would be revision asking the authors to specify the 1% reference, add an external converged heat flux comparison for the cylinder, and release code or data. I'd cite it for the r-versus-p trade-off study and the viscous HOIST formulation.","headline":"A solid, honest extension of HOIST to viscous flows; the cylinder heat-flux claim is credible but lacks an independent reference.","tokens_in":31184,"tokens_out":3534,"would_cite":true,"duration_ms":33047,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["65N30","65N50","76M10","76K05"],"pacs":[],"model":"deepseek-v4-flash","headline":"Viscous shock-dominated flows can be resolved without nonlinear stabilization by letting the mesh nodes move as part of the discrete solve, compressing elements into shocks and boundary layers.","keywords":["implicit shock tracking","rp-adaptivity","discontinuous Galerkin","hypersonic flow","viscous shock","boundary layer","heat flux prediction","mesh optimization"],"falsifier":"Apply the method with the paper's fixed parameters ($\\lambda=10$, $\\Delta=2$, same continuation schedule) to hypersonic flow over a different blunt geometry or Mach number: if the final pseudo-transient solve cannot drive the DG residual to about $10^{-9}$ without stabilization, or the heat flux profile departs by more than 1% from a high-fidelity reference, the central claim fails outside the demonstrated case.","tokens_in":30101,"feed_emoji":"🔥","tokens_out":12168,"duration_ms":101013,"temperature":0.7,"pith_summary":"The paper claims that viscous, shock-dominated flows can be resolved accurately without nonlinear stabilization if the computational mesh is allowed to move while the flow is being solved. The method treats the mesh nodal coordinates and the discrete flow variables as unknowns of a constrained optimization: the objective is an enriched discontinuous Galerkin residual, and the constraint is the IPDG equations, so minimizing drives elements into the viscous shocks and boundary layers. A p-adaptivity loop then raises the polynomial degree where the enriched residual is largest, reducing how much mesh compression is required. On hypersonic flow over a cylinder, the method predicts pointwise heat flux profiles, and for the integrated heat flux it reaches sub-1% error with fewer degrees of freedom than high-order h-adaptation.","feed_headline":"rp-adaptivity captures hypersonic heat flux without stabilization","feed_subtitle":"Moves mesh nodes into shocks and boundary layers, then p-refines where error remains; beats h-adaptation on heat-flux.","key_machinery":"The load-bearing object is a PDE-constrained optimization problem whose objective is the two-norm of the enriched IPDG residual plus a small mesh-distortion penalty, with the IPDG equations imposed as constraints. The enriched residual uses a test space two degrees higher than the trial space (enrichment degree $\\Delta=2$), making the objective sensitive to under-resolved viscous transitions. A constant amplification weight $\\rho=\\lambda$ on wall-adjacent elements forces boundary layers to compete with stronger shocks in the objective; viscosity and Reynolds-number continuation supplies robust initializations that prevent carbuncles; and a sequential quadratic programming solver with step-length limits, pseudo-transient step modifications, and enriched-residual-based Hessian regularization drives elements to extreme compression. The p-adaptive loop uses the enriched residual as an error indicator to raise polynomial degree locally in under-resolved regions.","core_discovery":"The central claim is that aggressive r-adaptation via implicit shock tracking can resolve the internal structure of viscous shocks and boundary layers on coarse meshes, making nonlinear stabilization unnecessary. The discrete solution and nodal coordinates are found together as the minimizer of the norm of an enriched interior-penalty DG residual (a test space two polynomial degrees above the trial space), subject to the DG equations; the minimizer compresses elements into the steep transitions. With viscosity continuation to avoid carbuncles, boundary-layer residual weighting to balance weak and strong features, and p-refinement driven by the same enriched residual, the method reproduces the Blasius boundary-layer profile and skin friction on a flat plate, and yields accurate heat flux profiles for $M_\\infty=5$ flow over a cylinder. In the cylinder test, the rp-adaptive method achieves sub-1% integrated heat-flux error with fewer degrees of freedom than h-adaptation at constant polynomial degrees 1, 2, and 3.","pith_inferences":["If the hand-set parameters ($\\lambda=10$ for the cylinder, $\\lambda=100$ for the flat plate, $\\Delta=2$, and the continuation schedule) could be chosen adaptively, the method would transfer more readily to new geometries; the paper's sensitivity experiments show these choices are load-bearing.","The reference-domain visualization, where compressed features are expanded, suggests the same machinery could be coupled with output-based error estimation to target specific quantities such as heating rates rather than the residual norm.","The approach may extend to three-dimensional and turbulent hypersonic flows, provided the optimizer can compress elements in three dimensions and the continuation path remains carbuncle-free; these extensions are not demonstrated here.","A stricter cost comparison that counts the expense of the continuation and optimization iterations, rather than only degrees of freedom, would place the rp-adaptive method and h-adaptation on more equal footing."],"forward_implications":["The hypersonic cylinder example indicates that heat flux, a notoriously sensitive output, can be predicted on coarse grids without artificial viscosity or limiting after several p-refinement rounds.","Because p-adaptivity shares the resolution burden, the required mesh compression is less extreme than in pure r-adaptation, so element stretching and polynomial degree can be balanced to save degrees of freedom.","Reynolds-number continuation is a necessary ingredient of the stabilization-free claim: without it, carbuncles form and the moving mesh tracks them.","Using the enriched residual as the p-adaptation indicator also keeps the shock-aligned mesh that allows deep pseudo-transient convergence without stabilization, something the dual-weighted-residual and feature-based indicators do not achieve.","On the flat-plate drag comparison, the rp-adaptive method is not more efficient than h-adaptation; its efficiency advantage appears in the cylinder heat-flux problem."],"supporting_citations":[{"why":"It supplies the original optimization-based implicit shock tracking formulation, minimizing an enriched residual over flow variables and mesh coordinates with a boundary-preserving mesh parametrization.","marker":"[54]"},{"why":"It provides the robust full-space SQP solver, merit function, and mesh-quality regularization that the viscous method modifies.","marker":"[26]"},{"why":"It introduces viscosity continuation and the viscous moving-DG context that this method builds on to compress elements into shocks and layers.","marker":"[31]"},{"why":"It defines the interior penalty DG discretization of the compressible Navier-Stokes equations used as the PDE constraint.","marker":"[23]"},{"why":"It supplies the PDE-based artificial viscosity shock capturing DG solution used to initialize the continuation process.","marker":"[4]"},{"why":"It provides the prior p-adaptive HOIST framework whose refinement loop is embedded and extended here.","marker":"[15]"},{"why":"It supports use of the lifted viscous flux for accurate and less oscillatory heat flux output via dual-consistent output evaluation.","marker":"[34, 22]"},{"why":"It supplies the dual-weighted residual error estimators used for the h-adaptive comparisons and studied as p-adaptation indicators.","marker":"[48, 49]"}],"fun_headline_variants":["rp-adaptive shock tracking resolves hypersonic heat flux","No stabilization needed: rp-adaptive method tracks shocks","Implicit shock tracking rp-adaptivity for viscous flows","Coordinate and solution optimized to resolve shock layers","rp-adaptivity with shock tracking outperforms h-adaptation"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The stabilization-free resolution rests on the hand-tuned choices of the boundary-layer amplification factor $\\lambda$, enrichment degree $\\Delta=2$, the viscosity continuation schedule, and the enriched-residual Hessian regularization producing sufficient mesh compression on any new viscous, shock-dominated problem; the paper demonstrates failures when these are changed, including a carbuncle without continuation and poor boundary-layer resolution with $\\Delta=1$ or $\\lambda=1$.","fun_headline_variants_meta":{"raw":{"variants":["rp-adaptive shock tracking resolves hypersonic heat flux","No stabilization needed: rp-adaptive method tracks shocks","Implicit shock tracking rp-adaptivity for viscous flows","Coordinate and solution optimized to resolve shock layers","rp-adaptivity with shock tracking outperforms h-adaptation"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000313,"raw_usage":{"total_tokens":1826,"prompt_tokens":1037,"completion_tokens":789,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":653,"completion_tokens_details":{"reasoning_tokens":708}},"tokens_in":653,"tokens_out":789,"duration_ms":6861,"temperature":1.0,"reasoning_tokens":708,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-16T11:30:41.147784+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Apply the method with the paper's fixed parameters ($\\lambda=10$, $\\Delta=2$, same continuation schedule) to hypersonic flow over a different blunt geometry or Mach number: if the final pseudo-transient solve cannot drive the DG residual to about $10^{-9}$ without stabilization, or the heat flux profile departs by more than 1% from a high-fidelity reference, the central claim fails outside the demonstrated case.","supporting_citations":[{"cited_title":"Zahr, Andrew Shi, and Per-Olof Persson","cited_arxiv_id":null,"evidence_quote":"It supplies the original optimization-based implicit shock tracking formulation, minimizing an enriched residual over flow variables and mesh coordinates with a boundary-preserving mesh parametrization."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"It provides the robust full-space SQP solver, merit function, and mesh-quality regularization that the viscous method modifies."},{"cited_title":"Kercher, Andrew Corrigan, and David A","cited_arxiv_id":null,"evidence_quote":"It introduces viscosity continuation and the viscous moving-DG context that this method builds on to compress elements into shocks and layers."},{"cited_title":"An optimal order interior penalty discontinuous Galerkin discretiza- tion of the compressible Navier–Stokes equations.Journal of Computational Physics, 227(22):9670–9685, 2008","cited_arxiv_id":null,"evidence_quote":"It defines the interior penalty DG discretization of the compressible Navier-Stokes equations used as the PDE constraint."},{"cited_title":"Barter and David L","cited_arxiv_id":null,"evidence_quote":"It supplies the PDE-based artificial viscosity shock capturing DG solution used to initialize the continuation process."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"It provides the prior p-adaptive HOIST framework whose refinement loop is embedded and extended here."}],"review_version":1}