{"id":"e74b042b-ce0f-4737-9641-be3092425c6a","arxiv_id":"2606.31152","paper_version":1,"verdict":"UNVERDICTED","confidence":"LOW","novelty_score":4.0,"correctness_risk":"unknown","formal_verification":"none","parameter_count":0,"one_line_summary":"A second-order GSAV consistent-splitting scheme for the 2D perturbed Boussinesq system is proposed with unconditional weak stability and optimal second-order error estimates, though error constants grow exponentially with inverse viscosity and diffusivity.","lead":"The paper proposes and analyzes a second-order consistent-splitting generalized scalar auxiliary variable scheme for the two-dimensional perturbed Boussinesq system. A smart generalist might read it to see how explicit treatment of nonlinear terms can yield unconditionally stable time-stepping for stratified fluid simulations.","discovery_kind":"unclear","skeptic_critique":{"model":"grok-4.3","headline":"Error constants show quadruply-nested exponential dependence on 1/ν and 1/κ, so second-order estimates are non-uniform despite unconditional stability.","rationale":"Reader flagged the explicit treatment of nonlinear terms as the weakest assumption, but the abstract already asserts that the GSAV construction yields unconditional stability. The more load-bearing issue for the headline claim is the non-robustness of the error constants, which the abstract itself surfaces. This moves the verdict from UNVERDICTED to CONDITIONAL pending verification that the exponential nesting cannot be reduced.","tokens_in":1680,"tokens_out":336,"duration_ms":26806,"concrete_test":"Locate the error-analysis section (likely after the stability theorem) and count the successive applications of Gronwall or discrete Gronwall inequalities that produce the exponential factors; recompute the constant dependence symbolically if the proof is written with explicit parameter tracking, and check whether any step can be replaced by a parameter-uniform estimate.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The abstract states that a careful tracing of the error analysis produces constants with quadruply-nested exponential growth in the reciprocals of viscosity and thermal diffusivity. For the central claim of 'optimal second-order error estimates' to be meaningful for the perturbed Boussinesq system (where small ν or κ are physically relevant), the estimates must remain useful as those parameters approach zero; the reported dependence violates uniformity. The unconditional weak stability theorem may hold, but the error result is the load-bearing claim whose practical content is undercut by this constant blow-up. No other internal inconsistency is visible from the given material.","agreement_with_reader":"disagree"},"referee_report":{"model":"grok-4.3","summary":"The manuscript proposes and analyzes a second-order consistent-splitting generalized scalar auxiliary variable (GSAV) scheme for the two-dimensional perturbed Boussinesq system. Nonlinear convection and advection terms, together with linear buoyancy and stratification couplings, are treated explicitly, reducing each time step to decoupled linear systems. The paper proves an unconditional weak stability result and derives optimal second-order error estimates for velocity, pressure, and temperature; a careful tracing of the constants shows quadruply-nested exponential dependence on the reciprocals of viscosity and thermal diffusivity. Numerical experiments confirm second-order convergence and reproduce expected long-time stratified-flow dynamics.","tokens_in":1831,"tokens_out":415,"duration_ms":24213,"significance":"If the stability and error results hold, the scheme offers an efficient, unconditionally stable discretization with decoupled solves for a physically relevant stratified-flow model. The explicit disclosure of the non-uniform error-constant dependence is a strength in transparency. The contribution extends prior consistent-splitting GSAV work on Navier-Stokes equations to the perturbed Boussinesq setting.","major_comments":[{"comment":"Error analysis section (the theorem establishing the second-order estimates): the claimed optimality is formally correct, yet the quadruply-nested exponential dependence of the constant on 1/ν and 1/κ (explicitly traced and stated in the abstract) renders the estimates non-uniform. This dependence is load-bearing for the practical content of the error result in the perturbed Boussinesq regime, where small viscosity and diffusivity are relevant; the manuscript should add a dedicated paragraph discussing the implications for robustness relative to fully implicit schemes.","section":"Error analysis section (the theorem establishing the second-order estimates)"}],"minor_comments":[{"comment":"The abstract states the scheme is for the 'two-dimensional' system but does not repeat this in the introduction; adding a brief sentence on the 2D setting would improve clarity.","section":"Introduction"}],"recommendation":"major_revision","confidential_remarks":null},"author_rebuttal":{"model":"grok-4.3","summary":"We thank the referee for the detailed and constructive report. We address the single major comment below and will incorporate the suggested addition in the revised manuscript.","responses":[{"response":"We agree that the non-uniformity of the error constants, already disclosed in the abstract, merits explicit discussion of its practical implications. In the revised manuscript we will insert a dedicated paragraph (likely in Section 4 or the concluding remarks) that contrasts the robustness of the present explicit-treatment GSAV scheme with that of fully implicit discretizations. The paragraph will note that the exponential dependence on 1/ν and 1/κ arises from the explicit handling of the buoyancy and nonlinear terms, which enables the decoupled linear solves and unconditional weak stability, while fully implicit schemes typically produce more uniform constants at the expense of solving coupled nonlinear systems at each step.","revision_made":"yes","referee_comment":"Error analysis section (the theorem establishing the second-order estimates): the claimed optimality is formally correct, yet the quadruply-nested exponential dependence of the constant on 1/ν and 1/κ (explicitly traced and stated in the abstract) renders the estimates non-uniform. This dependence is load-bearing for the practical content of the error result in the perturbed Boussinesq regime, where small viscosity and diffusivity are relevant; the manuscript should add a dedicated paragraph discussing the implications for robustness relative to fully implicit schemes."}],"tokens_in":1342,"tokens_out":308,"duration_ms":17883,"standing_objections":[]},"desk_editor":{"model":"grok-4.3","letter":"This paper takes the consistent-splitting GSAV-BDF2 scheme that Huang and Shen used for Navier-Stokes and applies it to the perturbed Boussinesq system. The nonlinear convection and advection terms are treated explicitly along with the buoyancy couplings, so each step breaks into a handful of decoupled linear solves. They prove unconditional weak stability and derive optimal second-order error estimates for velocity, pressure, and temperature.\n\nThe numerics confirm second-order convergence and show the expected internal-wave behavior plus relaxation to hydrostatic balance. The abstract is direct about the limitation: the error constant carries a quadruply-nested exponential dependence on the reciprocals of viscosity and thermal diffusivity, so the estimates are not uniform and become useless as either parameter approaches zero.\n\nThat dependence is the real soft spot. The stability result stands, but the error analysis loses practical content for the small-viscosity or small-diffusivity regimes that matter in stratified flows. No other internal problems are visible from the given material.\n\nThe work is aimed at numerical analysts who already follow SAV or consistent-splitting methods for incompressible fluids. A reader looking for another concrete application with proofs will get value from the details. It is a solid, incremental piece that deserves a serious referee.","headline":"Routine extension of GSAV-BDF2 to perturbed Boussinesq gives unconditional weak stability and second-order errors, but the error constants have quadruply-nested exponential blow-up in 1/ν and 1/κ.","tokens_in":2338,"tokens_out":340,"would_cite":false,"duration_ms":15954,"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 consistent-splitting GSAV scheme for the perturbed Boussinesq system is unconditionally weakly stable and yields optimal second-order error estimates.","keywords":["perturbed Boussinesq system","generalized scalar auxiliary variable","consistent-splitting scheme","unconditional stability","error estimates","time discretization","stratified flow","numerical scheme"],"falsifier":"A computation with fixed large time step and successively smaller viscosity (or thermal diffusivity) values that checks whether the discrete solutions remain bounded and converge at the claimed second-order rate or instead grow unbounded.","tokens_in":2577,"feed_emoji":"","tokens_out":720,"duration_ms":27331,"temperature":0.7,"pith_summary":"The paper develops a second-order time-stepping method for the two-dimensional perturbed Boussinesq system, which arises by subtracting a stable linearly stratified equilibrium from the standard Boussinesq equations. The discretization combines the generalized scalar auxiliary variable approach with consistent splitting, treating nonlinear convection and advection explicitly so that each step reduces to decoupled linear systems. An unconditional weak stability theorem is established, followed by proofs of optimal second-order convergence rates for velocity, pressure, and temperature. A reader would care because the method supports long-time simulations of stratified flows without artificial time-step restrictions, although the error bounds grow exponentially with inverse viscosity and inverse thermal diffusivity.","feed_headline":"GSAV scheme unconditionally stable for perturbed Boussinesq","feed_subtitle":"Consistent splitting reduces each step to decoupled linear solves while delivering second-order accuracy for velocity, pressure and temperat","key_machinery":"The consistent-splitting generalized BDF2 framework combined with the generalized scalar auxiliary variable (GSAV) approach, which reformulates the equations around an auxiliary scalar to permit explicit treatment of nonlinear terms while retaining stability.","core_discovery":"We propose and analyze a second-order consistent-splitting scheme based on the generalized scalar auxiliary variable approach for the two-dimensional perturbed Boussinesq system. The system is obtained by subtracting a stable, linearly stratified hydrostatic equilibrium from the standard Boussinesq system. The time discretization extends the consistent-splitting generalized BDF2 framework, treating the nonlinear convection and advection together with the linear buoyancy and stratification couplings explicitly, so that each time step reduces to a small number of decoupled linear systems. We prove an unconditional weak stability theorem for the GSAV scheme and derive optimal second-order error","pith_inferences":["The explicit treatment of convection may become impractical for very small viscosity, suggesting that selective implicit treatment could be needed for robustness at high Reynolds numbers.","The quadruply-nested exponential dependence in the error constant indicates that practical accuracy degrades rapidly as either viscosity or thermal diffusivity approaches zero.","The decoupling property may extend naturally to related systems such as the Navier-Stokes equations with temperature-dependent buoyancy."],"forward_implications":["Each time step reduces to a small number of decoupled linear systems.","Unconditional weak stability holds independently of the time-step size.","Optimal second-order error estimates are obtained for velocity, pressure, and temperature.","The scheme reproduces internal-wave dynamics and exponential relaxation to hydrostatic balance in long-time stratified-flow simulations."],"fun_headline_variants":["GSAV consistent splitting for perturbed Boussinesq","Second-order GSAV splitting for Boussinesq system","Unconditional weak stability for GSAV Boussinesq","Decoupled linear solves with consistent GSAV"],"cache_read_input_tokens":2112,"weakest_assumption_plain":"The nonlinear convection and advection can be treated explicitly together with the linear buoyancy and stratification couplings without destroying unconditional stability or second-order accuracy.","fun_headline_variants_meta":{"raw":{"variants":["GSAV consistent splitting for perturbed Boussinesq","Second-order GSAV splitting for Boussinesq system","Unconditional weak stability for GSAV Boussinesq","Decoupled linear solves with consistent GSAV"]},"model":"grok-4.3","cost_usd":0.005731,"raw_usage":{"total_tokens":2738,"prompt_tokens":676,"num_sources_used":0,"completion_tokens":61,"cost_in_usd_ticks":57312000,"prompt_tokens_details":{"text_tokens":676,"audio_tokens":0,"image_tokens":0,"cached_tokens":256},"completion_tokens_details":{"audio_tokens":0,"reasoning_tokens":2001,"accepted_prediction_tokens":0,"rejected_prediction_tokens":0}},"tokens_in":676,"tokens_out":61,"duration_ms":17512,"temperature":1.0,"reasoning_tokens":2001,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-07-01T05:03:48.894181+00:00","model_set":{"reader":"grok-4.3"},"falsifier":"A computation with fixed large time step and successively smaller viscosity (or thermal diffusivity) values that checks whether the discrete solutions remain bounded and converge at the claimed second-order rate or instead grow unbounded.","supporting_citations":[],"review_version":1}