{"id":"bc3fedaa-f85f-4535-8714-57526cd61b60","arxiv_id":"2506.21065","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"New entropy-stable in- and outflow boundary conditions for the compressible Navier-Stokes equations come with nonlinear a priori bounds and a finite-volume implementation.","lead":"Researchers propose inflow and outflow boundary conditions for the compressible Navier-Stokes equations and prove they keep entropy, mass, and total energy bounded. The work gives numerical analysts a stable way to let fluid enter and leave computational domains, for example in channel or open-air simulations.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Boundary entropy estimate is imported from [23] without proof; for the Navier-Stokes total-flux boundary conditions this transfer is asserted, not demonstrated, and Theorem 2.2 stands or falls on it.","rationale":"The reader identified the transfer of the boundary entropy estimate from the Euler paper [23] as the weakest assumption, and that is exactly the load-bearing point. The paper's Theorem 2.2 and the semi-discrete stability result both depend on the unstated lemma that the boundary integrand w^T f^b - psi.n is bounded below using only (12), (14), and positivity. The authors do not reproduce the argument and do not address the Navier-Stokes-specific feature that the boundary condition prescribes the total flux, so the solution pressure at a subsonic outflow is shifted by the normal viscous stress (Sec. 2.5). While I found no internal contradiction or obvious counterexample, the proof is not self-contained at the critical step, and the claim should be treated as conditional until the imported estimate is either supplied or independently verified. Since the reader already reached CONDITIONAL and my concern matches theirs, the verdict remains UNCHANGED.","tokens_in":15370,"tokens_out":14670,"duration_ms":179313,"concrete_test":"Independently derive the lower bound for the entropy boundary integrand for the four data fluxes (5)-(8), in normal/tangential coordinates, for a generic admissible Navier-Stokes state (rho, V, W, p). Verify, in the style of Eq. (30) of [23], that for each regime the integrand is bounded below by a constant plus terms controlled by the mass/energy outflow integrals (12), (14) and the inflow data, using only rho,p,T > 0. For the subsonic outflow, substitute the general boundary relation p = p_b + n^T tau n (the general form of Sec. 2.5) and confirm that the viscous normal stress cancels from the integrand or is otherwise harmless; if a term proportional to n^T tau n survives without control, Theorem 2.2 is false for viscous flows. Also check the corresponding discrete boundary term (29) on the Cartesian dual mesh.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim is Theorem 2.2, whose proof in Sec. 2.6 reduces the entropy boundary term in (15) to the assertion that the integrand w^T f^b - (psi_x,psi_y).n is bounded below using (12) and (14). The paper explicitly states that this integrand is 'exactly the same as (30) in [23]' and then says 'we omit the remaining details' (Sec. 2.6). This is a genuine reliance on an external result. In the Navier-Stokes problem, the boundary condition (4) prescribes the total flux n_1(f^I - f^V) + n_2(g^I - g^V), not the inviscid flux as in the Euler paper, and Sec. 2.5 shows this changes the relation between solution pressure and data pressure (for an x-normal boundary, p - f^V_2 = p_b). Thus the state w at the boundary is not the Euler state, and the boundedness argument in [23] does not transfer merely because the algebraic form of f^b is the same; one must verify that the estimate uses only the bounds (12), (14), and positivity, and that no uncontrolled boundary viscous stress term survives in w^T f^b - psi.n. The same imported estimate is used for the semi-discrete scheme in (29). If this external lemma fails, the entropy bound (15), the L^infinity(0,T;L^1) conclusion, and the discrete stability claim all fail. No internal inconsistency is apparent, but the proof is not self-contained at the decisive step.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper proposes in- and outflow boundary conditions for the compressible Navier-Stokes equations, generalizing the authors' earlier Euler boundary conditions. The new conditions prescribe the total (inviscid minus viscous) normal flux at the boundary in terms of data-dependent flux vectors, with distinct forms for supersonic/subsonic inflow and outflow. The main theoretical result (Theorem 2.2) asserts that admissible solutions satisfy a priori L∞(0,T;L1(Ω)) bounds on density, momentum, pressure, ρ log ρ, and ρ log T, derived from mass, energy, and entropy balances. A semi-discrete node-centered finite-volume scheme is presented that reproduces these estimates for the numerical solution, and numerical experiments with a convected vortex and a blast wave demonstrate robustness. The paper also openly states that linear well-posedness is not established for the subsonic inflow case, and that the decisive boundary entropy estimate is imported from the companion Euler paper [23] without a full derivation.","tokens_in":15755,"tokens_out":12497,"duration_ms":128230,"significance":"If the main theorem is correct, the paper provides a practically useful route to nonlinear a priori estimates for compressible Navier-Stokes initial-boundary value problems and to entropy-stable numerical schemes with open boundaries. The boundary conditions are new, conceptually simple, and they reduce to previously derived Euler conditions in the inviscid limit. The finite-volume construction is thoughtful and the numerical tests indicate that the method is stable in several nontrivial flow configurations. However, the central entropy estimate is not self-contained: it relies on an unstated lemma from [23], and the paper does not verify that the hypotheses of that lemma hold for the Navier-Stokes boundary state. Because this estimate is load-bearing for Theorem 2.2 and for the discrete stability claim, the paper is not yet ready for publication in its present form.","major_comments":[{"comment":"The entropy bound is not proven self-containedly. The proof asserts that the boundary integrand w^T f^b - (ψ_x, ψ_y)·n is \"exactly the same\" as (30) in [23] and omits the remaining details. However, the algebraic identity of f^b alone is not sufficient to transfer the boundedness argument. In the present problem the boundary condition (4) prescribes the total flux, so for an x-normal subsonic outflow the boundary state satisfies p - f^V_2 = p_b (Section 2.5), not p = p_b as in the Euler case. The entropy variables w and the potentials ψ are evaluated at the Navier-Stokes state, and the paper does not verify that the estimate in [23] uses only positivity, (12), (14), and the algebraic form of f^b. Without this verification, Theorem 2.2 is not established. Please supply the full argument or state and prove the required lemma explicitly.","section":"Section 2.6, Eq. (15)"},{"comment":"The semi-discrete entropy stability claim inherits the same gap. The boundary terms in (29) are asserted to be bounded below \"as shown in [23]\", but the discrete boundary flux (19) also imposes the total flux, and the discrete entropy variables are evaluated at the numerical solution, which satisfies a different pressure-data relation than in the Euler case. The same missing verification therefore affects the semi-discrete a priori bounds, and the proof should be completed or the needed result stated and proved.","section":"Section 3.3.3, Eq. (29)"}],"minor_comments":[{"comment":"The notation switches from lowercase ψ_x, ψ_y (defined just above) to capital Ψ_x, Ψ_y in (15); please use a consistent notation throughout.","section":"Section 2.6, Eq. (15)"},{"comment":"The summation in the entropy balance is written as \"NX i=N\", which appears to be a typo; it should presumably be the sum over all grid points i ∈ N.","section":"Section 3.3.3, Eq. (28)"},{"comment":"The numerical experiments measure reflections against a freestream solution but do not include a manufactured-solution grid-convergence study. The current tests show stability, but an accuracy verification would substantially strengthen the numerical section.","section":"Section 4.1, Table 2"},{"comment":"The paper transparently notes that linear well-posedness is not established for subsonic inflow. Given that this is one of the main boundary types, it would be helpful to briefly discuss the practical implications and whether the alternative supersonic-inflow data flux is actually recommended in that case.","section":"Section 2.7"}],"recommendation":"major_revision","confidential_remarks":"The decisive estimate in Theorem 2.2 is imported from the first author's companion paper [23]. Self-citation is not inappropriate here, but the editorial decision should hinge on whether the authors can supply the transferred lemma in full detail and verify its hypotheses for the Navier-Stokes total-flux boundary state. The numerical section is also thin on accuracy verification, and the admitted lack of linear well-posedness for subsonic inflow should be taken into account in judging the maturity of the proposal."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The genuinely new thing here is that the authors take the inflow/outflow data flux vectors from the Euler paper [23] and impose them on the total (inviscid minus viscous) flux of the Navier-Stokes equations. That changes the boundary state: at a subsonic outflow, for example, the condition becomes p - f^V_2 = p_b, not the inviscid pressure condition. The paper then proves mass, energy, and entropy bounds for the continuous problem and for a semi-discrete finite-volume scheme, and shows that the scheme is entropy stable with the new boundary treatment. The numerical experiments, mostly a vortex crossing the boundary and a blast wave, demonstrate stability and report reflections below 0.3% of the initial perturbation. All of this is a substantive and honest extension, not a repackaging of known results.\n\nThe soft spot is the proof of the entropy bound. In Section 2.6, the critical boundary integrand w^T f^b - psi.n is declared to be \"exactly the same as (30) in [23]\" and the remaining details are omitted. That is true as an algebraic expression in w, but the w that solves the Navier-Stokes IBVP is not the Euler w: the total-flux boundary condition changes how w relates to the data. It is entirely plausible that the boundedness argument in [23] only uses the bounds (12) and (14) plus positivity, in which case the transfer is valid. But the paper does not demonstrate that, and the theorem depends on it. This is a genuine gap in the proof, though not an obvious error. The authors also disclose that they have not proved linear well-posedness for subsonic inflow, which is a real limitation but one they flag openly. The numerical sections are stability-focused; there is no grid-convergence or manufactured-solution test for the boundary treatment itself, so accuracy of the boundary closure is not directly assessed.\n\nWho should read this: anyone working on entropy-stable schemes, open-boundary conditions, or nonlinear stability of compressible Navier-Stokes solvers. The paper belongs in the peer-reviewed literature, but a serious referee should send it back with a request to either supply the omitted boundary-entropy estimate or state explicitly which result in [23] covers the Navier-Stokes case. My recommendation: send to peer review, and ask for that proof before acceptance.","headline":"A credible extension of the Euler in/outflow entropy-stable BCs to Navier-Stokes, but the decisive boundary entropy estimate is imported from [23] rather than proved here.","tokens_in":16200,"tokens_out":2963,"would_cite":true,"duration_ms":32873,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["65M08","65M12","35Q30","76N10"],"pacs":[],"model":"deepseek-v4-flash","headline":"Entropy-stable in- and outflow boundary conditions are extended from inviscid to viscous compressible flow, with provable nonlinear bounds on mass, energy, and entropy.","keywords":["entropy stability","boundary conditions","compressible Navier-Stokes","finite volume","a priori estimates","inflow outflow","summation-by-parts","nonlinear stability"],"falsifier":"Run the proposed finite-volume boundary treatment on a subsonic outflow with an imposed strong boundary layer so that the viscous normal stress in the momentum flux is comparable to the pressure, and check whether the discrete entropy balance (28) stays bounded from below in time; an unbounded boundary term would show that the estimate imported from the Euler case does not transfer.","tokens_in":15185,"feed_emoji":"🌀","tokens_out":4942,"duration_ms":50682,"temperature":0.7,"pith_summary":"This paper proposes inflow and outflow boundary conditions for the compressible Navier-Stokes equations that are enforced by a boundary data flux vector rather than by setting solution values. It proves that admissible solutions satisfying these conditions have bounded mass, total energy, and entropy over finite time intervals, with density, momentum, pressure, and the logarithmic entropy terms in L∞(0,T;L1(Ω)). The same flux is then built into a node-centred finite-volume scheme using entropy-stable interior fluxes, and the discrete version of the three a priori bounds is shown to hold. The significance is that previously, entropy-stable open-boundary conditions of this type existed for the Euler equations and for far-field Navier-Stokes settings on infinite domains, while finite-domain inlet/outlet conditions for viscous flows lacked a nonlinear stability proof. The result would make boundary closures for viscous compressible simulations a matter of provable nonlinear bounds rather than empirical robustness.","feed_headline":"Entropy-stable in/outflow boundaries proven for Navier-Stokes","feed_subtitle":"Weakly imposed boundary fluxes give nonlinear a priori estimates for mass, energy, and entropy in finite-volume simulations.","key_machinery":"The load-bearing object is the boundary data flux vector [n1 fᵇ + n2 gᵇ], defined separately for supersonic inflow, subsonic inflow, subsonic outflow, and supersonic outflow, and inserted directly into the continuous boundary condition and into the finite-volume boundary fluxes. The entropy argument uses the entropy variables w and entropy flux potentials ψˣ, ψᵜ to rewrite the boundary term in the entropy balance as wᵀ fᵇ − ψ·n; since fᵇ is the same vector as in the Euler case [23], the right-hand side inherits the bound proven there, while the viscous term contributes non-negative entropy dissipation. In the scheme, the shuffle condition (18) on interior fluxes and the result that the discrete viscous term DIF F ≤ 0 leave only the same boundary expression to control.","core_discovery":"The central claim is that the four boundary conditions (4)-(8), written as a boundary data flux [n1 fᵇ + n2 gᵇ] that depends on the local flow regime, give a priori estimates for the Navier-Stokes initial-boundary value problem. Theorem 2.2 states that every admissible solution has {ρ, ρ|v|², p, ρ log ρ, ρ log T} in L∞(0,T;L¹(Ω)). The boundary conditions are designed so that at a subsonic outflow the momentum condition is p − fᵛ₂ = p_b, a pressure condition modified by the viscous normal stress, which reduces to the Euler pressure condition as viscosity tends to zero. The semi-discrete finite-volume scheme, with entropy-stable interior inviscid fluxes and the entropy-dissipative viscous treatment of [31], inherits the same bounds by reproducing the same boundary terms in its discrete entropy balance.","pith_inferences":["The omitted transfer of the Euler boundary entropy estimate to the viscous boundary pressure term is the point to scrutinise: if the estimate does not survive the replacement of p by p − fᵛ₂ at subsonic outflow, Theorem 2.2 would need modification.","Reflection levels near 0.3% for a weak vortex suggest the conditions are nearly transparent at low Mach numbers; whether this degrades with increasing Mach number or vortex strength is a natural testable extension.","The discrete bounds require positivity of ρ and p and a conservation-form scheme; relaxing those, as the paper notes for smooth solutions, would trade away the nonlinear estimates for linear ones, which is a design choice users would have to weigh."],"forward_implications":["The proposed boundary conditions can close entropy-stable finite-volume and other entropy-stable discretisations of the Navier-Stokes equations at inlets, outlets, and far-field boundaries of finite domains.","Simulations with strong non-smooth features, such as the circular blast wave in the paper, run stably under these boundary conditions without tuning.","The subsonic outflow condition is exactly the pressure condition proven linearly well-posed in [29], and the supersonic cases match that theory, so the linear well-posedness gap is confined to subsonic inflow.","The boundary treatment composes with entropy-stable no-slip wall conditions, letting a domain carry walls and open boundaries simultaneously while preserving the discrete bounds."],"supporting_citations":[{"why":"Supplies the Euler-equation boundary conditions and the key boundary entropy estimate that the present proof transfers to the Navier-Stokes setting.","marker":"[23]"},{"why":"Establishes linear well-posedness of the pressure outflow condition to which the subsonic outflow boundary condition reduces.","marker":"[29]"},{"why":"Provides the entropy-dissipative viscous flux approximation used in the finite-volume scheme, giving the DIF F ≤ 0 bound.","marker":"[31]"},{"why":"Supplies the entropy-stable inviscid fluxes used in the numerical experiments and the argument turning bounds on entropy, energy, and density into bounds on ρ log ρ and ρ log T.","marker":"[21]"},{"why":"Defines the entropy-stable no-slip wall condition shown to be compatible with the proposed open-boundary treatment.","marker":"[7]"},{"why":"Provides the earlier far-field Navier-Stokes boundary conditions that the new finite-domain conditions supersede.","marker":"[24]"}],"fun_headline_variants":["Entropy-stable boundaries proved for compressible NS","A priori estimates for Navier-Stokes inflow/outflow","Robust finite-volume scheme with stable boundary fluxes","Boundary conditions that guarantee entropy bounds","NS boundary fluxes with proven mass-energy-entropy control"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The proof borrows a boundary entropy bound from the Euler equations, where the boundary flux vector is identical, and leaves the corresponding viscous-case details out.","fun_headline_variants_meta":{"raw":{"variants":["Entropy-stable boundaries proved for compressible NS","A priori estimates for Navier-Stokes inflow/outflow","Robust finite-volume scheme with stable boundary fluxes","Boundary conditions that guarantee entropy bounds","NS boundary fluxes with proven mass-energy-entropy control"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000237,"raw_usage":{"total_tokens":1436,"prompt_tokens":805,"completion_tokens":631,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":421,"completion_tokens_details":{"reasoning_tokens":557}},"tokens_in":421,"tokens_out":631,"duration_ms":7490,"temperature":1.0,"reasoning_tokens":557,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-06T22:34:40.860423+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Run the proposed finite-volume boundary treatment on a subsonic outflow with an imposed strong boundary layer so that the viscous normal stress in the momentum flux is comparable to the pressure, and check whether the discrete entropy balance (28) stays bounded from below in time; an unbounded boundary term would show that the estimate imported from the Euler case does not transfer.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the Euler-equation boundary conditions and the key boundary entropy estimate that the present proof transfers to the Navier-Stokes setting."},{"cited_title":"Sv¨ ard, M","cited_arxiv_id":null,"evidence_quote":"Establishes linear well-posedness of the pressure outflow condition to which the subsonic outflow boundary condition reduces."},{"cited_title":"Sv¨ ard, M","cited_arxiv_id":null,"evidence_quote":"Provides the entropy-dissipative viscous flux approximation used in the finite-volume scheme, giving the DIF F ≤ 0 bound."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the entropy-stable inviscid fluxes used in the numerical experiments and the argument turning bounds on entropy, energy, and density into bounds on ρ log ρ and ρ log T."},{"cited_title":"Gjesteland and M","cited_arxiv_id":null,"evidence_quote":"Defines the entropy-stable no-slip wall condition shown to be compatible with the proposed open-boundary treatment."},{"cited_title":"Sv¨ ard and A","cited_arxiv_id":null,"evidence_quote":"Provides the earlier far-field Navier-Stokes boundary conditions that the new finite-domain conditions supersede."}],"review_version":1}