Quantitative critical homogenization of a 3D non-homogeneous thermoviscous fluid in perforated domains

In this paper, we study a three-dimensional non-homogeneous incompressible fluid in a bounded periodically perforated domain at the critical Stokes-capacity scale. The density is transported by the flow, while the viscosity depends on a quasi-static temperature governed by a uniformly elliptic transmission problem with a homogeneous Neumann boundary condition and zero spatial mean. In the critical regime, the homogenized momentum equation contains a Brinkman resistance term. Beyond the existing qualitative critical Brinkman limits, we prove a quantitative stability estimate for the coupled density--velocity--temperature system. For each fixed perforated geometry, global finite-energy weak solutions exist. Given a regular effective solution, every microscopic weak solution has squared relative error bounded by the initial mismatch plus an $O({\text{eps}}^2)$ remainder. For well-prepared data, this estimate yields order-$O(\text{eps})$ convergence of the density and uncorrected velocity in $L^{\infty }(0,T;L^2)$, of the temperature in $L^2(0,T;H^1)$, and of the corrected velocity in $L^2(0,T;H^1)$. The proof combines a solenoidal restriction operator, an $O(\text{eps})$ cell-capacity residual in the dual energy norm, an $O(\text{eps})$ $L^{6/5}$ corrector-gradient estimate, and fixed-domain temperature stability. The critical boundary layers retain order-one viscous energy, which is represented in the limit by the Brinkman dissipation. Finally, local regular solvability of the effective system is established for smooth prescribed forces satisfying a finite hierarchy of initial boundary compatibility conditions.