Regularity of Weak Solutions to the Inhomogeneous Stationary Navier–Stokes Equations

: One of the most intriguing issues in the mathematical theory of the stationary Navier– Stokes equations is the regularity of weak solutions. This problem has been deeply investigated for homogeneous ﬂuids. In this paper, the regularity of the solutions in the case of not constant viscosity is analyzed. Precisely, it is proved that for a bounded domain Ω ⊂ R 2 , a weak solution u ∈ W 1, q ( Ω ) is locally Hölder continuous if q = 2, and Hölder continuous around x , if q ∈ ( 1,2 ) and | µ ( x ) − µ 0 | is suitably small, with µ 0 positive constant; an analogous result holds true for a bounded domain Ω ⊂ R n ( n > 2 ) and weak solutions in W 1, n /2 ( Ω ) .


Introduction and Statement of the Results
The stationary motions of an inhomogeneous, incompressible viscous fluid in a solid (bounded domain with connected boundary) Ω of R n (n ≥ 2) are governed by the Navier-Stokes equations 2div (µ∇u) − div (u ⊗ u) − ∇p = 0, div u = 0, where u : Ω → R n , p : Ω → R, are the velocity and pressure field, respectively, ∇u = 1 2 (∇u + ∇u ) the symmetric part of ∇u, and µ : Ω → R + is the viscosity coefficient, satisfying the natural assumption almost everywhere in Ω. Unless otherwise specified, we will essentially use a standard notation as in [1]. In particular, B R (x) = {y ∈ R n : |y − x| < R}; if F is a function space the subscript σ in F σ stands for the (weakly) divergence-free condition; the integral mean of the field u over A will be denoted by u A = 1 |A| A u. Throughout the paper, the symbol c will denote a positive constant whose numerical value is unessential for our purposes and may change from line to line and in the same line, too.
Observe that we are considering the case of viscosity exclusively depending on the position x (see [2,3], for unsteady flows). This is the case, for example, of immiscible fluids or fluid mixtures where, in the first instance, we can neglect the natural dependence of viscosity on pressure and temperature (see, e.g., [4] and the references therein). Limiting ourselves to isothermal flows, the dependence of viscosity on pressure was already highlighted by G.G. Stokes in the seminal work [5] of 1845. In the last decades, many efforts have been done to take into account the effects of a pressure-dependent viscosity. In particular, we quote [6,7], where a global in time existence of solutions for the evolution problem was proved (a previous result [8] regarded a local in time existence and uniqueness result).
In [6,7], besides on pressure, a further dependence of the viscosity on the shear rate was required. More recently, a pure dependence of viscosity on pressure was instead supposed in [9], where an existence and uniqueness result for the stationary Stokes system, obtained from (1) by neglecting the inertial term, was presented by requiring general assumptions on µ that are satisfied by the Barus formula and other empiric laws; in the case of a thin straight pipe with variable cross-section the effective behavior of the flow was found via a rigorous asymptotic analysis with respect to the pipe's thickness [10].
In order to deal with the problem of regularity of weak solutions to system (1) with variable viscosity µ, in this paper, we consider the simplified situation of viscosity depending on the point µ = µ(x). Besides including meaningful cases as immiscible fluids in the model, it is worth pointing out that the problem here analyzed is also interesting from a purely mathematical point of view. Indeed, as far as we are aware, up to now the regularization problem has been considered only for variational solutions to (1) (see, e.g., [11], where it is considered a more general non linear system of the type of the Navier-Stokes one; however, it is not difficult to see that our technique can be used to deal with the equations here considered).
A weak solution (variational if q = 2) to (1) is a field u ∈ W 1,q σ (Ω), q ≥ 2n/(n + 2), which satisfies the relation To a solution to (3) is associated a scalar field p ∈ L q (Ω) such that [1] 2 It is natural to ask whether and under what condition a weak solution to (1) is more regular, for instance Hölder continuous. For homogeneous fluids, this issue has been the object of several papers, starting from the classical one of D. Serre (1983) [12], who proved regularity of a weak solution u ∈ W 1,q (Ω), q > 3/2, in the physical case n = 3. The regularity of variational solutions for n = 2, 3 was proved by O.A. Ladyzhenskaia [13] and for n = 4 by C. Gerhardt [14] (see also [11]). Regularity of variational solutions for n > 4 is an open problem [1]. More recently, Serre's results were generalized to distributional solutions in L n (Ω) by G. P. Galdi [1] and also in L n w (Ω) (weak L n (Ω) space) by H. Kim and H. Kozono [15], under the assumption of smallness of the norm u L n w (Ω) . The main purpose of this paper is to extend the results of [1] to system (1). To be precise, starting from the technique introduced in [16], we shall prove From Theorem 2 it easily follows Corollary 1. Let u ∈ W 1,n/2 (Ω), n > 2, be a weak solution to Equation (1). If µ ∈ C(Ω), then u ∈ C o,α loc (Ω) for every α ∈ (0, 1).

Preliminary Results
We collect in this section the main tools we need to prove Theorems 1 and 2 in the next sections. Consider the linear system [17] (see also [18,19] It is well-known that for constant µ a weak solution to (5) is real analytic [1].
, then there is a positive constant γ depending only on n, µ i and µ e such that for all 0 < ρ ≤ R R .
Proof. Let g(|y − x|) be a regular cut-off function, vanishing for |y − x| > 2R, equal to 1 for |y − x| < R and such that |∇g| ≤ cR −1 . The system where c is independent of R. In writing (8), we used the divergence free condition (1) 2 .
Since the field h + g 2 (u − κ) is permissible in (6), one has Therefore, taking into account (2), Since where we also took into account that, applying the divergence theorem and the condition g = 0 outside B 2R (x), By the arithmetic-geometric mean inequality and (8), for every α 1 , Therefore, choosing α 1 , α 2 1, (11) implies Using the properties of the function g and, in particular, that g = 1 in B R (x) and |∇g| ≤ cR −1 , it follows Hence, choosing k = u T R (x) , and using Poincaré's inequality ) (see, e.g., [20]), it follows Adding c B R (x) |∇u| 2 to both sides of (14) yields Now, to (15) we can apply a classical result in the theory of regularity of weak solutions to elliptic systems (see, e.g., [21], Lemma 8.23), to conclude that (7) holds, with γ = − log τ/ log 2.

Lemma 2. (Campanato's inequality)
-If µ is constant and u ∈ W 1,q (B R (x)) is a weak solution to (5), then for every p ∈ [1, +∞) for all 0 < ρ ≤ R << R and for some positive constant c depending only on n and p.

Lemma 3.
If F ∈ L q (B R ), q ∈ (1, +∞), then the system has a unique solution u ∈ W 1,q 0 (B R (x)) and for some constant c independent of R.
for every 0 ≤ ρ ≤ R =R − |x − y| and for every y ∈ BR(x), then u is Hölder continuous around x with exponent α.

Lemma 5.
Let Φ be a non negative and non decreasing function. Then there exists 0 such that if for all ρ ≤ R ≤R with c 0 , c 1 , α, β positive constants, β < α and < 0 , then for all ρ ≤ R ≤R, where c is a constant depending on α, β, c 0 .
• Let u ∈ W 1,q (Ω), q ∈ (1, 2), and let |µ − µ 0 | be suitably small in a neighborhood of x. Now u = u 1 + u 2 , with u 1 and u 2 solutions to and for every constant κ (this follows from the divergence free condition (1) 2 ). By Lemma 2 and by Lemma 3 By Hölder's inequality, Therefore, choosing κ = u B R (x) and using Sobolev's inequality, , one has Moreover, since with suitably small. Putting together (28) and (31) and using the inequality |a + b| q ≤ 2 q−1 (|a| q + |b| q ), for all a, b ∈ R, one gets Hence, Hölder continuity of u around x follows from Lemmas 4 and 5.

Proof of Theorem 2
Let u 1 and u 2 be the solutions to (25) and (26), respectively. By Lemma 2 and by Lemma 3 By Hölder's inequality, Therefore, by Sobolev's inequality, Moreover, since with suitably small. Putting together (33) and (36) and repeating the steps yielding (32), one has |∇u| n/2 , for small positive . Hence the desired result follows from Lemmas 4 and 5.
It should be of some interest to detect whether for n > 2 the above regularity results can be extended to function spaces larger than W 1,n/2 (Ω), like the grand Sobolev spaces introduced by T. Iwaniec and C. Sbordone [25]. In this connection see [26,27].
The interest in considering models with variable viscosity goes back to the researches by G.G. Stokes [5] and all the derived studies where, for instance, a pressure-dependent viscosity is considered (see [6][7][8][9][10] and the references therein for a qualitative analysis of the solutions). The idea is to use the results obtained in this paper, regarding the simplified case of viscosity depending on the point-that is, however, interesting itself, as pointed out in the introduction-to tackle the problem of regularity of solutions in the case of more general inhomogeneities, as the one due to the variation of viscosity with pressure. Obviously, a pressure-dependent viscosity brings an additional nonlinearity to system (1) and makes the pressure no more, merely a Lagrange multiplier, so that the mathematical analysis of the system is more complicated and this work can be considered a starting point in this perspective.

Conflicts of Interest:
The author declares no conflict of interest.