“…In this article, we focus on the spatially global estimates, i.e., on the entire domain. Since the equation for pressure is of degenerate parabolic type, finding L ∞ -bounds for its gradient, in general, is a difficult task and requires much work, see, for e.g., [5,18,19,24]. We will show that the gradient estimates for the particular flows under the current study can be obtained in a relatively simple way thanks to their equations' specific structure.…”