Time decay rates for the coupled modified Navier-Stokes and Maxwell equations on a half space

: This paper is concerned with time decay rates of the strong solutions of an incompressible the coupled modified Navier-Stokes and Maxwell equations in a half space R 3 + . With the use of the spectral decomposition of the Stokes operator and L p − L q estimates developed by Borchers and Miyakawa [2], we study the L 2 -decay rate of strong solutions.


Introduction
In this paper, we study the non-Newtonian fluids associated with Maxwell equations: Here u : Q T × (0, T ) → R 3 is the flow velocity vector, b : Q T × (0, T ) → R 3 is the magnetic vector, π : Q T × (0, T ) → R is the total pressure and Du is the symmetric part of the velocity gradient, i.e. Du = D i j u := 1 2 To motivate the conditions on the stress tensor S , we recall the following examples of constitutive laws S (Du) = (µ 0 + µ 1 |Du| p−2 )Du, (1.2) where µ 0 ≥ 0 and µ 1 > 0 are constants (see e.g. [1,16]). We consider the initial-boundary value problem of (1.1), which requires initial conditions u(x, 0) = u 0 (x) and b(x, 0) = b 0 (x) together with the boundary conditions defined as follows: u = 0 and b · n = 0, (∇ × b) × n = 0, (1.4) where n is the outward unit normal vector along boundary ∂R 3 + . The fluids such as coal-water, faint, soaps, etc., is not followed linear relationship between rate of strain and shear stress. In languages of mathematics, commonly, the standard Navier-Stokes equations refers to the equations of motion of an incompressible fluid type of S = µ 0 Du, where µ 0 > 0 is constant. On the other hand, one class of non-Newtonian fluids is defined by S = µ(|Du|)Du (µ(·) a positive nonlinear function). For example, we note that For the existence of solutions to the model with the stress tensor (1.5), we make a brief comment. After the pioneering study done by Ladyzhenskaya [10], the global-in-time existence of strong solutions is proved for q > 11 5 in [11]. On the other hand, they also established the small data global-in-time existence of strong solution for 5 3 < q < ∞ in three dimensional space (see also [3] for weak solutions) The Eqs (1.1), which are the generalized incompressible magnetohydrodynamics equations is regarded as one of the simplest model describing the dynamics of electrically conducting liquid with involved rheological structure in a magnetic field. For the model (1.1), Gunzburger et al. in [6] considered (1.1)-(1.4) for the case of bounded or periodic domains, and they established unique solvability of the initial-boundary value problem. More specifically, assuming that u 0 ∈ H 2 (Ω) and b 0 ∈ H 1 (Ω) with the following boundary conditions (1.4) for a bounded domain, it was shown in [6] that if 5 2 < q ≤ 6, a generalized solution exists, and moreover, it satisfies u ∈ L ∞ (0, T ; L 2 (Ω) ∩ H 1 (Ω)), ∇u ∈ L ∞ (0, T ; L q (Ω)), b ∈ L ∞ (0, T ; L 2 (Ω) ∩ H 1 (Ω)), ∇b ∈ L ∞ (0, T ; L 2 (Ω)), b ∈ L 2 (0, T ; H 2 (Ω)).
(1.6) Furthermore, they shown the uniqueness of solutions. Here strong solutions means that solutions satisfy (1.1) pointwise a.e. and the energy equality holds. Recently, the authors in [8] establish global unique solvability to (1.1)-(1.4) for u 0 ∈ (W 1,2 ∩ W 1,p ) and b 0 ∈ W 1,2 , 5 2 ≤ p in the same class above (see [12] for weak solutions). For a half space, the proof in [6] is also held. And thus, we will not comment further on the existence and uniqueness of strong solutions to (1.1) and (1.2).
For the asymptotic behavior of strong solutions to (1.1) and (1.2), the author in [9] recently examined the L 2 -algebraic decay in the whole space R 3 with respect to the monopolar shear thickening fluids using Fourier splitting method in [13]. Precisely, he shown that We also refer to [5,7] for Navier-Stokes equations of non-Newtonian type.
On the other hand, in the case of a half space R n + , Dong and Chen [4] obtain that the weak solution u(t) of the Naiver-Stokes equations of non-Newtonian type, that is b ≡ 0 in (1.1) enjoys the optimal algebraic decay estimates by using the L p − L q estimates based on the explicit solution formula to the Stokes equation given by Ukai [15] and the spectral decomposition method and fractional powers of the Stokes operator derived by Borchers and Miyakawa [2]. For a domain, we briefly some comments. For a whole space, we obtain the same result like as Theorem 1.1 by the fourier splitting method (see [8]). However, in R 3 + case with boundary conditions, we needs to be handled the solution form due to the boundary effects. To effectively deal with boundary effect, we use the well-known spectral method, we would like to obtain temporal rate of the strong solutions. So far, there are few known results on the time decay problem to (1.1) in R 3 + . In this direction, our main results are as follows: Suppose that u is a strong solution of Naiver-Stokes equations of non-Newtonian type . Remark 1.4. Comparing to [4], since the result of Corollary 1.3 is about the time decay rate for strong solution to Naiver-Stokes equations of non-Newtonian type, the restriction to the range of p can be slightly relaxed.
Let us rewrite the abstract formulation of (1.1), where the Stokes operators A s and Laplacian operator A to be specially considered the boundary conditions (1.4) are defined as follows: (see [14, page 262]). And also the bilinear operator B 1 and B 2 is defined as follows: Here P is the orthogonal projection mapping L q (R 3 + ) onto L q σ (R 3 + ) (see, e.g. Ukai [15]). Lemma 1.5. ([2, Theorem 3.6]). Let either 1 < r ≤ q < ∞ or 1 ≤ r < q ≤ ∞ and v ∈ L 2 σ ∩ L r . Then we have the L r -L q estimate ∥∇ k e −At v∥ L q ≤ ct − k where e −At denotes the analytic semigroup generated by the Stokes operator A. Here P is the orthogonal projection mapping L q (R 3 + ) onto L q σ (R 3 + ) (see, for example, Ukai [12]).

Proof of Theorems
Proof. Following [2, Theorem 3.6], we note that for 1 < r < ∞, α > 0 and 0 < 1 The proof is almost same to that in [4, Lemma 2.3]. For the convenience of readers, we give a proof. Indeed, a proof is based on (2.1) and Gagliardo-Nirenberg inequality as follows: Similarly, we have which is complete of the proof. □ Proof of Theorem 1. Following the argument in [4], we will prove Theorem. From the energy inequality, we know To obtain a lower bound of the second term in (2.3), we get for ρ > 0, and And thus, we have To estimate the right-hand side of (2.4), we consider the integral form of (1.9), Applying the operator E(ρ) to the both sides of this integral equation and integrating by parts, we obtain This together with Lemma 2.1 implies ∥E(ρ)u(t)−E(ρ)e −tA s u 0 ∥ ≤ cρ 5 4 t 0 e −ρ(t−s) ∥u∥ 2 +∥∇u∥ p−1 p−1 +cρ where we use Korn's inequality and the following inequality in the last inequality: For 11 and for p ≥ 3 Hence, we have In the same way, we get Here, we use the following estimate Now let ρ = 3(1 + t) −1 and multiply both sides above by (1 + t) 3 to obtain And thus, we get Since ∥e −tA s u 0 ∥ 2 L 2 → 0 and ∥e −tB b 0 ∥ 2 L 2 → 0 as t → ∞, we conclude that u(t) → 0 and b(t) → 0 as t → ∞ and thus it complete the proof of first part in Theorem 1.1.
For second part, it follows from (2.9) and Lemma 2.2 that Let ρ = 3t −1 and then multiply the both sides of this equation by t 3 to get for t ≥ 1, and complete the proof of Theorem 3.1.
Proof of Theorem 2.
Let ρ = 3t −1 and multiply both sides above by t 3 to obtain The proof is complete.