Nitsche's method for the stationary Boussinesq system under mixed and nonlinear boundary conditions

Aparna Bansal; Dwijendra N. Pandey; Gianmarco Sperone; Nicol\'as A. Barnafi

arxiv: 2604.06560 · v1 · submitted 2026-04-08 · 🧮 math.NA · cs.NA

Nitsche's method for the stationary Boussinesq system under mixed and nonlinear boundary conditions

Aparna Bansal , Nicol\'as A. Barnafi , Gianmarco Sperone , Dwijendra N. Pandey This is my paper

Pith reviewed 2026-05-10 18:30 UTC · model grok-4.3

classification 🧮 math.NA cs.NA

keywords Nitsche's methodBoussinesq systemNavier's slipnonlinear boundary conditionsfinite element methoda posteriori error estimatorsconvergence rateswell-posedness

0 comments

The pith

Nitsche's method yields a well-posed finite element scheme for the stationary Boussinesq system with Navier slip and nonlinear boundary conditions.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper develops a Nitsche-type finite element discretization for the stationary Boussinesq system that incorporates Navier slip and a nonlinear boundary condition. It proves that the resulting discrete nonlinear problem is well-posed when the data satisfies a standard smallness assumption, by applying fixed-point theorems. The analysis also establishes optimal convergence rates for the approximation errors in velocity, pressure, and temperature, together with the reliability and efficiency of residual-based a posteriori error estimators. A reader would care because the method remains robust on domains with arbitrarily complex boundaries, removing the need for body-fitted meshes in buoyancy-driven flow simulations.

Core claim

The central claim is that Nitsche's method applied to the stationary Boussinesq system with Navier's slip and nonlinear boundary conditions produces a robust finite element scheme. Under a standard smallness assumption on the data, fixed-point theorems establish well-posedness of the discrete problem. Optimal convergence rates are proved for the approximation error, and residual-based a posteriori error estimators are shown to be both efficient and reliable. The theoretical results are confirmed by numerical tests.

What carries the argument

Nitsche's method, a symmetric penalty-based finite element technique that weakly enforces the mixed and nonlinear boundary conditions without requiring the mesh to conform to the boundary.

If this is right

The scheme remains stable and convergent on domains with arbitrarily complicated boundaries.
Optimal error rates hold simultaneously for velocity, pressure, and temperature.
Residual-based estimators can be used to drive adaptive refinement while preserving reliability.
The numerical tests confirm that the predicted rates and estimator properties appear in practice.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

The same Nitsche treatment could be tried for time-dependent or non-stationary Boussinesq problems with similar boundary conditions.
The approach suggests a route to handle other coupled nonlinear flow models on non-smooth or multiply-connected domains.
Relaxing the small-data restriction would require a different existence proof, perhaps based on topological degree or monotonicity arguments.

Load-bearing premise

The data must be small enough for a fixed-point argument to guarantee existence and uniqueness of the discrete solution.

What would settle it

A computation in which the data violates the smallness condition yet the nonlinear solver still produces a unique convergent solution, or in which the error estimators lose reliability, would contradict the main claims.

Figures

Figures reproduced from arXiv: 2604.06560 by Aparna Bansal, Dwijendra N. Pandey, Gianmarco Sperone, Nicol\'as A. Barnafi.

**Figure 2.** Figure 2: Test 3: Plots of numerical solutions of the velocity magnitude [PITH_FULL_IMAGE:figures/full_fig_p022_2.png] view at source ↗

**Figure 3.** Figure 3: Test 3: Plots of numerical solutions of the velocity magnitude [PITH_FULL_IMAGE:figures/full_fig_p023_3.png] view at source ↗

**Figure 4.** Figure 4: Test 3: Global indicators for (a) L-shape (b) T-shape. [PITH_FULL_IMAGE:figures/full_fig_p023_4.png] view at source ↗

**Figure 5.** Figure 5: Test 4: Velocity Magnitude, Velocity, pressure isovalues, and temperature isovalues, respectively. [PITH_FULL_IMAGE:figures/full_fig_p024_5.png] view at source ↗

read the original abstract

In this paper we analyze Nitsche's method for the stationary Boussinesq system with Navier's slip and a nonlinear boundary condition. Our analysis of the formulation establishes the robustness of a finite elements scheme in arbitrarily complex boundaries. The well-posedness of the discrete problem is established using fixed-point theorems under a standard smallness assumption on the data. We also provide optimal convergence rates for the approximation error. Furthermore, the efficiency and reliability of residual-based a posteriori error estimators are established. We validate our theory through several numerical tests.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Desk Editor's Note private letter to a colleague

This paper gives a Nitsche scheme for the stationary Boussinesq system with Navier slip and nonlinear boundary conditions, plus well-posedness, optimal rates, and a posteriori estimators under a small-data assumption.

read the letter

The paper extends Nitsche's method to the stationary Boussinesq system with Navier slip plus nonlinear boundary conditions on complex domains. The specific analysis of the discrete problem, including residual-based a posteriori estimators, is new relative to the cited literature on Nitsche for fluids or Boussinesq alone. It does a clean job setting up the formulation, applying fixed-point arguments for existence and uniqueness, and deriving the expected convergence rates. The numerical tests on several examples are a plus because they show the scheme behaves as predicted on irregular boundaries. The robustness claim for arbitrarily complex domains is practically useful in CFD settings. The smallness assumption on the data is the main soft spot. The stress-test note is on target: if the threshold constant in the fixed-point argument depends on the mesh size h or the Nitsche penalty through inverse inequalities or Lipschitz estimates on the convective and buoyancy terms, then existence is not guaranteed uniformly under refinement. The paper uses standard arguments, so the constants are probably independent, but that needs explicit confirmation in the proofs rather than being left implicit. The a posteriori results and efficiency claims rest on the same foundation, so any gap there would propagate. This is aimed at numerical analysts working on finite-element methods for buoyancy-driven flows with slip-type boundaries. Readers who need a ready-to-implement scheme for such problems will get value from the analysis and the tests. It deserves peer review because the core extension is grounded and the numerics are present, even though the uniformity of the smallness condition should be checked carefully by referees.

Referee Report

1 major / 2 minor

Summary. The manuscript analyzes Nitsche's method for the stationary Boussinesq system with Navier slip and nonlinear boundary conditions. It establishes well-posedness of the discrete problem via fixed-point theorems under a standard small-data assumption, derives optimal a priori convergence rates, proves reliability and efficiency of residual-based a posteriori estimators, and validates the results with numerical experiments.

Significance. If the claims hold, the work provides a robust unfitted discretization for buoyancy-driven flows on complex domains together with a posteriori control, which is useful for adaptive simulations. The combination of well-posedness, optimal rates, and estimator analysis under nonlinear boundary conditions is a solid contribution to numerical analysis of nonlinear fluid problems.

major comments (1)

[Well-posedness of the discrete problem] Well-posedness analysis (fixed-point argument for the discrete nonlinear problem): the smallness assumption on the data must be shown to be independent of the mesh size h and the Nitsche penalty parameter. If the threshold constant deteriorates with h, existence and uniqueness of the discrete solution are not guaranteed uniformly under refinement, which is load-bearing for the subsequent optimal convergence rates and a posteriori results that presuppose a discrete solution exists.

minor comments (2)

The abstract and introduction refer to 'several numerical tests' but the validation section would benefit from explicit tables or figures reporting observed convergence rates, effectivity indices for the estimators, and the specific values of the small-data parameter used.
Notation for the nonlinear boundary condition and the Nitsche penalty terms should be introduced with an early equation reference to improve readability for readers unfamiliar with the formulation.

Simulated Author's Rebuttal

1 responses · 0 unresolved

We thank the referee for the positive assessment of our work and for the detailed comment on the well-posedness analysis. We address this point below.

read point-by-point responses

Referee: [Well-posedness of the discrete problem] Well-posedness analysis (fixed-point argument for the discrete nonlinear problem): the smallness assumption on the data must be shown to be independent of the mesh size h and the Nitsche penalty parameter. If the threshold constant deteriorates with h, existence and uniqueness of the discrete solution are not guaranteed uniformly under refinement, which is load-bearing for the subsequent optimal convergence rates and a posteriori results that presuppose a discrete solution exists.

Authors: We agree that this independence is crucial. Upon reviewing our proof of well-posedness (Section 3), the smallness condition is indeed independent of h. The fixed-point map is shown to be a contraction in a ball whose radius depends on the data, with the contraction constant controlled by terms that are bounded uniformly in h thanks to the mesh-independent stability of the Nitsche formulation (the penalty parameter γ is chosen larger than a constant independent of h, and all inverse inequalities and trace theorems yield h-independent constants in this context). We will revise the manuscript to explicitly state this fact and include a short paragraph clarifying the uniformity with respect to h and γ. revision: yes

Circularity Check

0 steps flagged

No circularity: well-posedness, convergence, and a-posteriori results rest on standard fixed-point theorems and residual analysis independent of the paper's own constructs.

full rationale

The derivation proceeds by applying Nitsche's method to the stationary Boussinesq system, then invoking standard fixed-point theorems (Schauder or contraction mapping) under a small-data hypothesis to obtain discrete well-posedness, followed by standard Céa-type arguments for optimal convergence rates and classical residual-based a-posteriori theory for reliability/efficiency. None of these steps reduces by construction to a fitted parameter, self-definition, or self-citation chain; the smallness assumption is an external hypothesis whose uniformity is a separate analytic question, not a tautology. The numerical validation is confirmatory rather than load-bearing for the claims. The analysis is therefore self-contained against external benchmarks.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

Based solely on the abstract, the paper relies on standard functional-analytic assumptions for Navier-Stokes-type systems and a small-data hypothesis; no free parameters or invented entities are explicitly introduced.

axioms (2)

domain assumption Smallness assumption on the data
Invoked to guarantee existence via fixed-point theorem for the nonlinear discrete problem.
standard math Standard Sobolev space setting and trace theorems for the velocity-temperature coupling
Required for the variational formulation and error analysis of the Boussinesq system.

pith-pipeline@v0.9.0 · 5402 in / 1409 out tokens · 43623 ms · 2026-05-10T18:30:26.866217+00:00 · methodology

discussion (0)

Reference graph

Works this paper leans on

3 extracted references · 3 canonical work pages

[1]

Bansal, N

[BBP24] A. Bansal, N. A. Barnafi, and D. N. Pandey. Nitsche method for Navier-Stokes equations with slip boundary conditions: convergence analysis and VMS-LES stabilization. ESAIM. Mathematical Modelling and Numerical Analysis , 58(5):2079–2115,

work page 2079
[2]

[GOGKT25] P. A. Gazca-Orozco, F. Gmeineder, E. M. Kokavcová, and T. Tscherpel. A Nitsche method for in- compressible fluids with general dynamic boundary conditions. ArXiv preprint, arXiv:2502.09550,

work page arXiv
[3]

[Shi84] Z. C. Shi. On the convergence rate of the boundary penalty method. International Journal for Numerical Methods in Engineering , 20(11):2027–2032,

work page 2027

[1] [1]

Bansal, N

[BBP24] A. Bansal, N. A. Barnafi, and D. N. Pandey. Nitsche method for Navier-Stokes equations with slip boundary conditions: convergence analysis and VMS-LES stabilization. ESAIM. Mathematical Modelling and Numerical Analysis , 58(5):2079–2115,

work page 2079

[2] [2]

[GOGKT25] P. A. Gazca-Orozco, F. Gmeineder, E. M. Kokavcová, and T. Tscherpel. A Nitsche method for in- compressible fluids with general dynamic boundary conditions. ArXiv preprint, arXiv:2502.09550,

work page arXiv

[3] [3]

[Shi84] Z. C. Shi. On the convergence rate of the boundary penalty method. International Journal for Numerical Methods in Engineering , 20(11):2027–2032,

work page 2027