Analysis of Navier–Stokes Models for Flows in Bidisperse Porous Media

2 years ago

(original), ISI/JCR, Nonlinear analysis, paper, PDEs (theory)

Abstract

Having in view a model proposed by Nield and Kuznetsov (Transp Porous Media 59:325–339, 2005; Transp Porous Media 96:495–499, 2013), we consider a more general system of coupled Navier–Stokes type equations in the incompressible case subject to the homogeneous Dirichlet condition in a bounded domain. We provide a deep theoretical analysis for large classes of equations and coupled systems of Navier–Stokes type with various non-homogeneous terms of reaction type. Existence results are obtained by using a variational approach making use of several fixed point principles and matrix theory.

Authors

Mirela Kohr
Department of Mathematics Babes-Bolyai University, Cluj-Napoca, Romania

Radu Precup
Faculty of Mathematics and Computer Science and Institute of Advanced Studies in Science and Technology, Babeş–Bolyai University
Tiberiu Popoviciu Institute of Numerical Analysis, Romanian Academy

Keywords

Navier–Stokes equations; Bidisperse porous media; Fixed point technique.

Paper coordinates

M. Kohr, R. Precup, Analysis of Navier-Stokes models for flows in bidisperse porous media, J. Math. Fluid Mechanics, 25 (2023) art. no. 38, https://doi.org/10.1007/s00021-023-00784-w

PDF

https://link.springer.com/article/10.1007/s00021-023-00784-w

About this paper

Journal

Journal of Mathematical Fluid Mechanics

Publisher Name

Springer

DOI

https://doi.org/10.1007/s00021-023-00784-w

Print ISSN

1422-6928

Online ISSN

1422-6952

google scholar link

Paper (preprint) in HTML form

Analysis of Navier-Stokes models for flows in bidisperse porous media

Mirela Kohr, Radu Precup M. Kohr, Faculty of Mathematics and Computer Science, Babeş–Bolyai University, 400084 Cluj-Napoca, Romania mkohr@math.ubbcluj.ro R. Precup, Faculty of Mathematics and Computer Science and Institute of Advanced Studies in Science and Technology, Babeş–Bolyai University, 400084 Cluj-Napoca, Romania & Tiberiu Popoviciu Institute of Numerical Analysis, Romanian Academy, P.O. Box 68-1, 400110 Cluj-Napoca, Romania r.precup@math.ubbcluj.ro

Abstract.

Having in view a model proposed by Nield and Kuznetsov (2005, 2013), we consider a more general system of coupled Navier-Stokes type equations in the incompressible case subject to the homogeneous Dirichlet condition in a bounded domain. We provide a deep theoretical analysis for large classes of equations and coupled systems of Navier-Stokes type with various non-homogeneous terms of reaction type. Existence results are obtained by using a variational approach making use of several fixed point principles and matrix theory.

Key words and phrases:

Navier-Stokes equations; bidisperse porous media; fixed point technique

1991 Mathematics Subject Classification:

35Q30, 76D05

In memory of Professor Gabriela Kohr, with deep respect

1. Introduction

The Stokes and Navier-Stokes systems play a main role in various areas of fluid mechanics, engineering, biology, chemistry, and there is a huge list of references concerning the mathematical analysis of related boundary value problems and of their applications. Among of them, we mention the books [4], [5], [9], [13], [24], [25], [29], [41], [42], [45], [46].

Let $N\geq 2$ and $\Omega$ be an open set in $\mathbb{R}^{N}$ . Let $\eta,\kappa\geq 0$ and $\mu>0$ be given constants. Let $\mathbf{u}$ and $p$ be unknown vector and scalar fields. Let us assume that $\mathbf{f}$ is a given vector field defined on $\Omega$ . Then the equations

(1.1)

\begin{array}[]{lll}-\mu\Delta\mathbf{u}+\eta\mathbf{u}+\kappa\left(\mathbf{u}% \cdot\nabla\right)\mathbf{u}+\nabla p=\mathbf{f},\ \mathrm{div}\ \mathbf{u}=0% \ \mbox{ in }\ \Omega&&\end{array}

determine a Navier-Stokes type system in the incompressible framework. If $\eta=0$ and $\kappa>0$ we then obtain the well-known Navier-Stokes system in the incompressible case, while for $\eta=\kappa=0$ , (1.1) becomes the Stokes system, which is an Agmon-Douglis-Nirenberg elliptic and linear system (see, e.g., [13], [16], [29], [45] for further details).

Extensions to a more general case of anisotropic Stokes and Navier-Stokes systems with $L^{\infty}$ -variable coefficients, and the analysis of various boundary value problems involving them can be consulted in [18], [19], [20], [21] and the references therein.

Fabes, Kenig and Verchota [12] used a layer potential approach in the analysis of the Dirichlet problem for the Stokes system on Lipschitz domains in the Euclidean setting (see also [10] for applications of the layer potential approach for strongly elliptic differential operators). Dindos̆ and Mitrea [11] proved the well-posedness in Sobolev and Besov spaces for the Dirichlet problem for the Stokes and Navier-Stokes systems with smooth coefficients in Lipschitz domains on compact Riemannian manifolds. Mitrea and Wright [27] obtained well-posedness results in Sobolev and Besov spaces for Dirichlet problems for the Stokes system with constant coefficients in Lipschitz domains in ${\mathbb{R}}^{n}$ for Dirichlet problems for the Stokes, Oseen and Navier-Stokes systems with constant coefficients in a non-solenoidal framework (see also the references therein).

Korobkov, Pileckas and Russo [22] studied the flux problem in the theory of steady Navier-Stokes equations with constant coefficients and non-homogeneous boundary conditions. Amrouche and Rodríguez-Bellido [1] proved the existence of a very weak solution for the non-homogeneous Dirichlet problem for the compressible Navier-Stokes system in a bounded domain of the class $C^{1,1}$ in ${\mathbb{R}}^{3}$ .

Bulí $\check{c}$ ek, Málek and $\check{Z}$ abenský [6] studied a boundary value problem with homogeneous Dirichlet condition associated with a system of nonlinear partial differential equations that generalize the classical fluid flow models of Stokes, Darcy, Forchheimer and Brinkman, by assuming that the viscosity and the drag coefficient depend on the shear rate and the pressure. The authors proved the existence of weak solutions to the problem under a minimal number of conditions, and analyzed relevant examples of viscosities and drag coefficients modeling real physical situations.

The authors in [18] analyzed in $L^{2}$ -based Sobolev spaces, the non-homogeneous boundary value problems of Dirichlet-transmission type for the anisotropic Stokes and Navier-Stokes systems in a compressible framework in a bounded Lipschitz domain with a transversal Lipschitz interface in ${\mathbb{R}}^{n}$ , $n\geq 2$ ( $n=2,3$ for the nonlinear problems). They proved the existence of a weak solution to the Dirichlet problem and the Dirichlet-transmission problem for the nonlinear anisotropic Navier-Stokes system by using the Leray-Schauder fixed point theorem and various results and estimates from the linear case, as well as the Leray-Hopf inequality and some other norm inequalities. Explicit conditions for uniqueness of solutions to the nonlinear problems have been also provided. Mixed problems and mixed-transmission problems for the anisotropic Stokes and Navier-Stokes systems in bounded Lipschitz domains with transversal Lipschitz interfaces have been considered in [19] and analyzed from the variational point of view. The authors in [17] used a layer potential approach and the Leray-Schauder fixed point theorem and proved existence results for a nonlinear Neumann-transmission problem for the Stokes and Brinkman systems in $L^{p}$ , Sobolev, and Besov spaces. Mazzucato and Nistor [26] obtained well-posedness and regularity results in Sobolev spaces for the linear elasticity equations in the anisotropic case with mixed boundary conditions on polyhedral domains.

1.1. Bidisperse porous media and some related models

A bidisperse porous medium (BDPM) may be described as a standard porous medium in which the solid phase is replaced by another porous medium. Thus, a BDPM can be viewed as a medium composed of clusters of large particles that are agglomerations of small particles (cf. [8], see also [32], [33]). The voids between the clusters are macro-pores and the voids within the clusters, which are much smaller in size, are micro-pores. We can then define the $f$ -phase (the macro-pores) and the $p$ -phase (the remainder of the structure). Bidisperse adsorbent or bidisperse capillary wicks in a heat pipe are practical applications related to bidisperse porous media. There are also various biological structures, such as bone regeneration scaffolds, that can be described in terms of bidisperse porous media (see also [33]).

Extending the Brinkman model for a monodisperse porous medium, Nield and Kuznetsov [30] considered a model to describe the steady-state momentum transfer in a BDPM by the following pair of coupled equations for $v_{f}^{\ast}$ and $v_{p}^{\ast}$ ,

\left\{\begin{array}[]{l}\mathbf{G}=\frac{\mu}{K_{f}}\mathbf{v}_{f}^{\ast}+% \zeta(\mathbf{v}_{f}^{\ast}-\mathbf{v}_{p}^{\ast})-\widetilde{\mu}_{f}\Delta^{% \ast}\mathbf{v}_{f}^{\ast}\\ \mathbf{G}=\frac{\mu}{K_{p}}\mathbf{v}_{p}^{\ast}+\zeta(\mathbf{v}_{p}^{\ast}-% \mathbf{v}_{f}^{\ast})-\widetilde{\mu}_{p}\Delta^{\ast}\mathbf{v}_{p}^{\ast}\,% ,\end{array}\right.

where the asterisks denote dimensional variables, $\mathbf{G}$ is the negative of the applied pressure gradient, $\mu$ is the fluid viscosity, $K_{p}$ and $K_{f}$ are the permeabilities of the two phases, $\zeta$ is the coefficient for momentum transfer between the two phases, and $\widetilde{\mu}_{f}$ and $\widetilde{\mu}_{p}$ are the effective viscosities of the two phases (cf. also [32]). In the model proposed by Nield and Kuznetsov [30], the quadratic or Forchheimer terms $|\mathbf{v}_{p}^{\ast}|\mathbf{v}_{p}^{\ast}$ and $\mathbf{v}_{f}^{\ast}|\mathbf{v}_{f}^{\ast}$ have been neglected. This model and various extensions have been considered in many studies related to forced, natural and mixed convection. Among them, we mention [7], [23], [28], [31], [33], [40], [39], [43] (see also the references therein). The thermal convection in an anisotropic bidisperse porous medium has been investigated in [44].

Nield and Kuznetsov [33] extended their linear model proposed in [30] by adding some semilinear terms, called the Forchheimer drag terms, as follows

\left\{\begin{array}[]{l}\mathbf{G}=\frac{\mu}{K_{f}}\mathbf{v}_{f}^{\ast}+% \zeta(\mathbf{v}_{f}^{\ast}-\mathbf{v}_{p}^{\ast})-\widetilde{\mu}_{f}\Delta^{% \ast}\mathbf{v}_{f}^{\ast}+\frac{c_{f}\rho}{K_{f}^{1/2}}|\mathbf{v}_{f}^{\ast}% |\mathbf{v}_{f}^{\ast}\\ \mathbf{G}=\frac{\mu}{K_{p}}\mathbf{v}_{p}^{\ast}+\zeta(\mathbf{v}_{p}^{\ast}-% \mathbf{v}_{f}^{\ast})-\widetilde{\mu}_{p}\Delta^{\ast}\mathbf{v}_{p}^{\ast}+% \frac{c_{p}\rho}{K_{p}^{1/2}}|\mathbf{v}_{p}^{\ast}|\mathbf{v}_{p}^{\ast}\,,% \end{array}\right.

where $\rho$ is the density of the fluid, and $c_{f}$ and $c_{p}$ are the Forchheimer coefficients.

The Nield-Kuznetsov models described above are based on the same pressure gradient $-{\bf G}$ in both phases. Other models consider possible different pressures in the macro and micro phases. For instance, Straughan [44] having analyzed a model of thermal convection in an anisotropic bidispersive porous medium with permeability tensors in the macro and micro phases, considers different velocities ${\bf U}^{f}$ and ${\bf U}^{p}$ and different pressures $p^{f}$ and $p^{p}$ in the macro and micro-pores (see also [7] and the references therein for similar models of bidisperse porous media with different velocities and different pressures in macro and micro phases).

Having in view the model of Nield and Kuznetsov [33], where the steady-state momentum transfer is described by the previous semilinear system, and also the model of Straughan [44], we consider a more general nonlinear coupled type Navier-Stokes system arising in the analysis of fluid flows in bidisperse porous media. Thus, our paper is build around the following system

\left\{\begin{array}[]{l}-\mu_{1}\Delta\mathbf{u}_{1}+\eta_{1}\mathbf{u}_{1}+% \kappa_{1}\left(\mathbf{u}_{1}\cdot\nabla\right)\mathbf{u}_{1}+\nabla p_{1}=% \mathbf{h}_{1}-\ \alpha_{1}\left|\mathbf{u}_{1}\right|^{p-1}\mathbf{u}_{1}\ -% \gamma_{1}\left(\mathbf{u}_{1}-\mathbf{u}_{2}\right)\ \ \ \text{in }\Omega\\ -\mu_{2}\Delta\mathbf{u}_{2}+\eta_{2}\mathbf{u}_{2}+\kappa_{2}\left(\mathbf{u}% _{2}\cdot\nabla\right)\mathbf{u}_{2}+\nabla p_{2}=\mathbf{h}_{2}-\ \alpha_{2}% \left|\mathbf{u}_{2}\right|^{p-1}\mathbf{u}_{2}\ -\gamma_{2}\left(\mathbf{u}_{% 2}-\mathbf{u}_{1}\right)\ \ \ \text{in }\Omega\\ \mathrm{div}\text{\ }\mathbf{u}_{i}=0\ \ \ \text{in }\Omega\\ \mathbf{u}_{i}=0\ \ \ \text{on }\partial\Omega\ \ \ \left(i=1,2\right),\end{% array}\right.

where $\Omega\subset\mathbb{R}^{N}$ is a bounded domain ( $N\leq 3$ ), $p\geq 1$ , $\eta_{i},\,\kappa_{i},\,\alpha_{i}\geq 0$ and $\mu_{i},\,\gamma_{i}>0$ , $i=1,2$ , are given constants whose meaning depends on the physical properties of fluid flow and porous medium, while $\mathbf{h}_{i}$ , $i=1,2$ , are given data in some Sobolev spaces.

In order to analyze this system, we provide a deep analysis of a homogeneous Dirichlet problem of more general coupled Navier-Stokes systems with various non-homogeneous terms of reaction type, and obtain existence results by using a variational approach combined with fixed point theorems, a technique already used for other classes of equations (see [34, Ch. 6], [35], [37, Chs. 9-12], [38]).

The paper is structured as follows. First, we mention some well-known but useful results regarding the stationary Navies-Stokes equations in the incompressible case. The next section is devoted to the analysis of the homogeneous Dirichlet problem for the Navier-Stokes equations with reaction terms. We obtain existence results based on the Schauder fixed point theorem and the Leray-Schauder fixed point theorem. Uniqueness result can be also obtained by using the Banach contraction principle, by imposing additional conditions to the reaction terms. The third section is devoted to the analysis of a coupled system of Navier-Stokes equations. The last section is devoted to a coupled system that could describe a fluid flow in a bidisperse porous medium. We obtain related existence and uniqueness results that follow as consequences of the results obtained in the previous sections.

1.2. Stationary Navier-Stokes type equations

Let $N\leq 3$ and $\Omega\subset\mathbb{R}^{N}$ be an bounded domain.

Next we recall some well-known results about the system

(1.2)

\left\{\begin{array}[]{l}-\mu\Delta\mathbf{u}+\eta\mathbf{u}+\kappa\left(% \mathbf{u}\cdot\nabla\right)\mathbf{u}+\nabla p=\mathbf{f}\ \ \ \text{in }% \Omega\\ \mathrm{div\ }\mathbf{u}=0\ \ \ \text{in }\Omega\\ \mathbf{u}=0\ \ \ \text{on }\partial\Omega\end{array}\right.

where $\mu>0,\,\kappa,\,\eta\geq 0$ are given constants and $\mathbf{f}\in H^{-1}\left(\Omega\right)^{N}$ is a given distribution.

The variational form of system (1.2) with the unknown pair $\left(\mathbf{u},p\right)\in V\times L^{2}\left(\Omega\right)$ is

(1.3)

\displaystyle\mu\left(\mathbf{u},\mathbf{v}\right)_{H_{0}^{1}}+\eta\left(% \mathbf{u},\mathbf{v}\right)_{L^{2}}+\kappa b\left(\mathbf{u},\mathbf{u},% \mathbf{v}\right)-\left(p,\mathrm{div}\ \mathbf{v}\right)_{L^{2}}

\displaystyle=

\displaystyle\left(\mathbf{f},\mathbf{v}\right),\ \ \ \mathbf{v}\in H_{0}^{1}% \left(\Omega\right)^{N},

where

$\displaystyle\left(\mathbf{u},\mathbf{v}\right)_{L^{2}}$	$\displaystyle=$	$\displaystyle\sum\limits_{i=1}^{N}\int_{\Omega}{u}_{i}v_{i}\,,\ \ \ \ \left(% \mathbf{u},\mathbf{v}\right)_{H_{0}^{1}}=\sum\limits_{i=1}^{N}\int_{\Omega}% \nabla{u}_{i}\cdot\nabla v_{i}$
$\displaystyle b\left(\mathbf{u},\mathbf{v},\mathbf{w}\right)$	$\displaystyle=$	$\displaystyle\sum\limits_{i,j=1}^{N}\int_{\Omega}u_{j}\frac{\partial v_{i}}{% \partial x_{j}}w_{i}\,,\ \ \ \ \left(\mathbf{f},\mathbf{v}\right)=\sum\limits_% {i=1}^{N}\left(f_{i},v_{i}\right),$
$\displaystyle V$	$\displaystyle=$	$\displaystyle\left\{\mathbf{v}=\left(v_{1},\ldots,v_{N}\right)\in H_{0}^{1}% \left(\Omega\right)^{N}:\ \mathrm{div\ }\mathbf{v}=0\right\}.$

For $\mathbf{u}\in V,$ equation (1.3) gives

(1.4)

\mu\left(\mathbf{u},\mathbf{v}\right)_{H_{0}^{1}}+\eta\left(\mathbf{u},\mathbf% {v}\right)_{L^{2}}+\kappa b\left(\mathbf{u},\mathbf{u},\mathbf{v}\right)=\left% (\mathbf{f},\mathbf{v}\right),\ \ \ \mathbf{v}\in V.

Once a solution $\mathbf{u}\in V$ to (1.4) is found, the pressure $p\in L^{2}\left(\Omega\right)$ is guaranteed by De Rham’s Theorem (cf., e.g., [45, Proposition 1.1, Chapter 1], [15, Theorem 2.3, Chapter 1], see also [2], [3, Theorem 2.1]).

On $H_{0}^{1}\left(\Omega\right)^{N}$ consider the inner product and norm

\left(\mathbf{u},\mathbf{v}\right)_{0}:=\int_{\Omega}\left(\mu\nabla\mathbf{u}% \cdot\nabla\mathbf{v}+\eta\mathbf{u}\cdot\mathbf{v}\right),\ \ \ \left|\mathbf% {u}\right|^{2}_{0}:=\int_{\Omega}\left(\mu\left|\nabla\mathbf{u}\right|^{2}+% \eta\left|\mathbf{u}\right|^{2}\right),

which when applied to the subspace $V$ will be denoted by $\left(\cdot,\cdot\right)_{V}$ and $\left|\cdot\right|_{V}$ , respectively. Then the embedding constants and the corresponding inequalities for the inclusions $\ V\subset L^{2}\left(\Omega\right)^{N}\subset V^{\prime}$ are

(1.5)

\left|\mathbf{u}\right|_{L^{2}}\leq\frac{1}{\sqrt{\mu\lambda_{1}+\eta}}\left|% \mathbf{u}\right|_{V}\ \ \ \left(\mathbf{u}\in V\right),\ \ \ \left|\mathbf{h}% \right|_{V^{\prime}}\leq\frac{1}{\sqrt{\mu\lambda_{1}+\eta}}\left|\mathbf{h}% \right|_{L^{2}}\ \ \left(\mathbf{h}\in L^{2}\left(\Omega\right)^{N}\right),

where $\lambda_{1}$ is the first eigenvalue of $-\Delta$ with respect to the homogeneous Dirichlet problem. Indeed, knowing that

\left|\mathbf{u}\right|_{L^{2}}\leq\frac{1}{\sqrt{\lambda_{1}}}\left|\mathbf{u% }\right|_{H_{0}^{1}}\ \ \ \left(\mathbf{u}\in H_{0}^{1}\left(\Omega\right)^{N}% \right),\ \ \ \left|\mathbf{h}\right|_{H^{-1}}\leq\frac{1}{\sqrt{\lambda_{1}}}% \left|\mathbf{h}\right|_{L^{2}}\ \ \left(\mathbf{h}\in L^{2}\left(\Omega\right% )^{N}\right),

we have

(1.6)

\left|\mathbf{u}\right|_{V}^{2}=\mu\left|\mathbf{u}\right|_{H_{0}^{1}}^{2}+% \eta\left|\mathbf{u}\right|_{L^{2}}^{2}\geq\left(\mu\lambda_{1}+\eta\right)% \left|\mathbf{u}\right|_{L^{2}}^{2},

whence the first inequality in (1.5). Based on (1.6), the second inequality is obtained as follows:

\left|\mathbf{h}\right|_{V^{\prime}}=\sup_{\mathbf{v}\in V}\frac{\left|\left(% \mathbf{h},\mathbf{v}\right)\right|}{\left|\mathbf{v}\right|_{V}}\leq\left|% \mathbf{h}\right|_{L^{2}}\sup_{\mathbf{v}\in V}\frac{\left|\mathbf{v}\right|_{% L^{2}}}{\left|\mathbf{v}\right|_{V}}\leq\frac{1}{\sqrt{\mu\lambda_{1}+\eta}}% \left|\mathbf{h}\right|_{L^{2}}.

Also $\ V\subset H_{0}^{1}\left(\Omega\right)^{N}\subset H^{-1}\left(\Omega\right)^{% N}\subset V^{\prime}$ and

(1.7)

\left|\mathbf{h}\right|_{V^{\prime}}=\sup_{\mathbf{v}\in V}\frac{\left|\left(% \mathbf{h},\mathbf{v}\right)\right|}{\left|\mathbf{v}\right|_{V}}\leq\sup_{% \mathbf{v}\in H_{0}^{1}\left(\Omega\right)^{N}}\frac{\left|\left(\mathbf{h},% \mathbf{v}\right)\right|}{\left|\mathbf{v}\right|_{V}}=\left|\mathbf{h}\right|% _{H^{-1}}\ \ \left(\mathbf{h}\in H^{-1}\left(\Omega\right)^{N}\right).

Recall that the trilinear functional $\ b:V\times V\times V\rightarrow\mathbb{R}\$ satisfies

	$\displaystyle b\left(\mathbf{u},\mathbf{v},\mathbf{w}\right)+b\left(\mathbf{u}% ,\mathbf{w},\mathbf{v}\right)$	$\displaystyle=$	$\displaystyle 0\ \ \ \ \left(\mathbf{u},\mathbf{v},\mathbf{w}\in V\right),$
	$\displaystyle\left\|b\left(\mathbf{u},\mathbf{v},\mathbf{w}\right)\right\|$	$\displaystyle\leq$	$\displaystyle M\left\|\mathbf{u}\right\|_{V}\left\|\mathbf{v}\right\|_{V}\left\|% \mathbf{w}\right\|_{V},$

where $M>0$ is a constant depending on $\mu$ and $\eta.$ Also, using the Galerkin method, one can prove that for every $\mathbf{f}\in V^{\prime}$ , equation $(\ref{2})$ has at least one solution $\mathbf{u}\in V$ (see, e.g., [13], [45]).

Uniqueness: For every $\mathbf{f}\in V^{\prime}{},$ with $\left|\mathbf{f}\right|_{V^{\prime}}<1/\left(\kappa M\right),$ equation (1.4) has at most one solution $\mathbf{u}\in V.$ Indeed, if $\mathbf{u}_{1},\mathbf{u}_{2}\in V$ are solutions and we let $\mathbf{u}=\mathbf{u}_{1}-\mathbf{u}_{2},$ then using (1.4) one has

$\displaystyle 0$	$\displaystyle=$	$\displaystyle\left(\mathbf{u},\mathbf{v}\right)_{V}+\kappa b\left(\mathbf{u}_{% 1},\mathbf{u}_{1},\mathbf{v}\right)-\kappa b\left(\mathbf{u}_{2},\mathbf{u}_{2% },\mathbf{v}\right)$
	$\displaystyle=$	$\displaystyle\left(\mathbf{u},\mathbf{v}\right)_{V}+\kappa\left(b\left(\mathbf% {u}_{1},\mathbf{u}_{1},\mathbf{v}\right)-b\left(\mathbf{u}_{1},\mathbf{u}_{2},% \mathbf{v}\right)+b\left(\mathbf{u}_{1},\mathbf{u}_{2},\mathbf{v}\right)-b% \left(\mathbf{u}_{2},\mathbf{u}_{2},\mathbf{v}\right)\right)$
	$\displaystyle=$	$\displaystyle\left(\mathbf{u},\mathbf{v}\right)_{V}+\kappa\left(b\left(\mathbf% {u}_{1},\mathbf{u},\mathbf{v}\right)+b\left(\mathbf{u},\mathbf{u}_{2},\mathbf{% v}\right)\right),$

which for $\mathbf{v}=\mathbf{u},$ since $b\left(\mathbf{u}_{1},\mathbf{u},\mathbf{u}\right)=0,$ gives

\left|\mathbf{u}\right|_{V}^{2}=-\kappa b\left(\mathbf{u},\mathbf{u}_{2},% \mathbf{u}\right)\leq\kappa M\left|\mathbf{u}\right|_{V}^{2}\left|\mathbf{u}_{% 2}\right|_{V}.

On the other hand for $\mathbf{u}_{2},$ taking in (1.4) $\mathbf{v}=\mathbf{u}=\mathbf{u}_{2},$ one has

\left|\mathbf{u}_{2}\right|_{V}^{2}=\left(\mathbf{f},\mathbf{u}_{2}\right)\leq% \left|\mathbf{f}\right|_{V^{\prime}}\left|\mathbf{u}_{2}\right|_{V}.

Then

\left|\mathbf{u}\right|_{V}^{2}\left(1-\kappa M\left|\mathbf{f}\right|_{V^{% \prime}}\right)\leq 0,

which for $\left|\mathbf{f}\right|_{V^{\prime}}<1/\left(\kappa M\right)$ yields $\left|\mathbf{u}\right|_{V}=0,$ that is $\mathbf{u}_{1}=\mathbf{u}_{2}.$ Denote the unique solution by $\mathbf{u}_{\mathbf{f}}.$

Thus we may define the solution operator

S:D_{0}\rightarrow V,\ \ \ \ \mathbf{f}\longmapsto\mathbf{u}_{\mathbf{f}}.

Here $D_{0}:=\left\{\mathbf{f}\in V^{\prime}:\ \kappa M\left|\mathbf{f}\right|_{V^{% \prime}}<1\right\}.$ In addition, since the trilinear functional $b$ satisfies

b\left(\mathbf{u},\mathbf{v},\mathbf{w}\right)+b\left(\mathbf{u},\mathbf{w},% \mathbf{v}\right)=0,

one has $b\left(\mathbf{u},\mathbf{u},\mathbf{u}\right)=0,$ whence taking in (1.4) $\mathbf{v}=\mathbf{u}=S\left(\mathbf{f}\right)$ we see that

(1.8)

\left|S\left(\mathbf{f}\right)\right|_{V}\leq\left|\mathbf{f}\right|_{V^{% \prime}}.

Also using the linearity of $b$ in each of its variables gives

			$\displaystyle b\left(\mathbf{u},\mathbf{u},\mathbf{w}\right)-b\left(\mathbf{v}% ,\mathbf{v},\mathbf{w}\right)$
		$\displaystyle=$	$\displaystyle b\left(\mathbf{u},\mathbf{u},\mathbf{w}\right)-b\left(\mathbf{v}% ,\mathbf{u},\mathbf{w}\right)+b\left(\mathbf{v},\mathbf{u},\mathbf{w}\right)-b% \left(\mathbf{v},\mathbf{v},\mathbf{w}\right)$
		$\displaystyle=$	$\displaystyle b\left(\mathbf{u}-\mathbf{v},\mathbf{u},\mathbf{w}\right)+b\left% (\mathbf{v},\mathbf{u}-\mathbf{v},\mathbf{w}\right).$

Then

$\displaystyle\left\|S\left(\mathbf{f}\right)-S\left(\mathbf{g}\right)\right\|_{V% }^{2}$	$\displaystyle=$	$\displaystyle\left(\mathbf{f}-\mathbf{g},S\left(\mathbf{f}\right)-S\left(% \mathbf{g}\right)\right)$
		$\displaystyle+\kappa b\left(S\left(\mathbf{g}\right),S\left(\mathbf{g}\right),% S\left(\mathbf{f}\right)-S\left(\mathbf{g}\right)\right)-\kappa b\left(S\left(% \mathbf{f}\right),S\left(\mathbf{f}\right),S\left(\mathbf{f}\right)-S\left(% \mathbf{g}\right)\right)$
	$\displaystyle=$	$\displaystyle\left(\mathbf{f}-\mathbf{g},S\left(\mathbf{f}\right)-S\left(% \mathbf{g}\right)\right)+\kappa b\left(S\left(\mathbf{g}\right)-S\left(\mathbf% {f}\right),S\left(\mathbf{g}\right),S\left(\mathbf{f}\right)-S\left(\mathbf{g}% \right)\right)$
		$\displaystyle+\kappa b\left(S\left(\mathbf{f}\right),S\left(\mathbf{g}\right)-% S\left(\mathbf{f}\right),S\left(\mathbf{f}\right)-S\left(\mathbf{g}\right)% \right).$

Hence

	$\displaystyle\left\|S\left(\mathbf{f}\right)-S\left(\mathbf{g}\right)\right\|_{V% }^{2}$	$\displaystyle\leq$	$\displaystyle\left\|\mathbf{f}-\mathbf{g}\right\|_{V^{\prime}}\left\|S\left(% \mathbf{f}\right)-S\left(\mathbf{g}\right)\right\|_{V}+\kappa M\left\|S\left(% \mathbf{f}\right)-S\left(\mathbf{g}\right)\right\|_{V}^{2}\left(\left\|S\left(% \mathbf{g}\right)\right\|_{V}+\left\|S\left(\mathbf{f}\right)\right\|_{V}\right)$
		$\displaystyle\leq$	$\displaystyle\left\|\mathbf{f}-\mathbf{g}\right\|_{V^{\prime}}\left\|S\left(% \mathbf{f}\right)-S\left(g\right)\right\|_{V}+\kappa M\left(\left\|{\bf f}\right% \|_{V^{\prime}}+\left\|\mathbf{g}\right\|_{V^{\prime}}\right)\left\|S\left(\mathbf% {f}\right)-S\left(\mathbf{g}\right)\right\|_{V}^{2}.$

As a result, if $\ \left|\mathbf{f}\right|_{V^{\prime}},\left|\mathbf{g}\right|_{V^{\prime}}% \leq\rho,$ then

(1.9)

\left(1-2\kappa M\rho\right)\left|S\left(\mathbf{f}\right)-S\left(\mathbf{g}% \right)\right|_{V}\leq\left|\mathbf{f}-\mathbf{g}\right|_{V^{\prime}}.

Therefore, if $\ 1-2\kappa M\rho>0,$ i.e., $\ \rho<\frac{1}{2\kappa M},$ then the operator $S$ is Lipschitz continuous on the ball of $V^{\prime}$ centered at the origin and of radius $\rho$ . Note that in the case $\kappa=0$ , inequality (1.9) shows that the solution operator $S$ is Lipschitz continuous on the entire space $V^{\prime}$ . Let us introduce a notation for the Lipschitz constant, namely

L\left(\rho\right):=\left(1-2\kappa M\rho\right)^{-1}.

Note that, in the case of the Stokes system, one has $\kappa=\eta=0$ and thus $S$ is well-defined and Lipschitz continuous on the whole space $V^{\prime}$ .

1.3. Notions of vector analysis in fixed point theory

In this paper we assume that the reader is familiar with Banach contraction principle and Schauder and Leray-Schauder fixed point theorems. However, we consider it useful to present some less known elements of vector analysis in fixed point theory. For more details we refer the reader to [36].

A square matrix $\mathcal{A}$ with nonnegative entries is said to be convergent to zero if its power $\mathcal{A}^{k}$ tends to the zero matrix as $k$ tends to infinity. This property is equivalent to the fact that the spectral radius of $\mathcal{A}$ is less than one, and also to the property that $\mathcal{J}-\mathcal{A}$ is invertible and its inverse also has nonnegative entries (here $\mathcal{J}$ is the unit matrix).

A matrix $\mathcal{A}=\left[a_{ij}\right]_{i,j=1,2}$ of size two is convergent to zero if and only if

(1.10)

\text{tr}\ \mathcal{A}<\min\left\{2,\ 1+\text{det\ }\mathcal{A}\right\}.

The property of being convergent to zero of a matrix $\mathcal{A}$ is useful to pass from a matrix inequality of the form $\left(\mathcal{J}-\mathcal{A}\right)\mathbf{u}\leq\mathbf{v},$ where $\mathbf{u},\mathbf{v}$ are column vectors and the inequality is understood on components, to the inequality $\mathbf{u}\leq\left(\mathcal{J}-\mathcal{A}\right)^{-1}\mathbf{v},$ without change of inequality. The notion is even more important since it extends to matrices the situation on real numbers $0\leq a<1,$ asked on the Lipschitz constant in Banach contraction theorem. More exactly we have the following result, a special case of the more general Perov’s fixed point theorem:

Let $k\geq 1$ be a given integer, $\left(X,\left|.\right|_{X}\right)$ be a Banach space, and $X^{k}:=\underbrace{X\times\cdots\times X}_{{\rm{k-times}}}$ . If $D\subset X^{k}$ is closed, and $T:D\rightarrow D,$ $T=\left(T_{1},\cdot\cdot\cdot,T_{k}\right)$ is a mapping satisfying the matrix condition

\left[\begin{array}[]{c}\left|T_{1}\left(\mathbf{u}\right)-T_{1}\left(\mathbf{% v}\right)\right|_{X}\\ \vdots\\ \left|T_{k}\left(\mathbf{u}\right)-T_{k}\left(\mathbf{v}\right)\right|_{X}\end% {array}\right]\leq\mathcal{A}\left[\begin{array}[]{c}\left|u_{1}-v_{1}\right|_% {X}\\ \vdots\\ \left|u_{k}-v_{k}\right|_{X}\end{array}\right]

for all $\mathbf{u}=\left(u_{1},\cdot\cdot\cdot,u_{k}\right),\ \mathbf{v}=\left(v_{1},% \cdot\cdot\cdot,v_{k}\right)\in D$ and some matrix $\mathcal{A}$ , which is convergent to zero, then $T$ has a unique fixed point $\mathbf{u}\in D,$ i.e., $T_{i}\left(\mathbf{u}\right)=u_{i}$ for $i=1,\ldots,k.$

2. Navier-Stokes type equations with reaction terms

Consider now the problem

(2.1)

\left\{\begin{array}[]{l}-\mu\Delta\mathbf{u}+\eta\mathbf{u}+\kappa\left(% \mathbf{u}\cdot\nabla\right)\mathbf{u}+\nabla p=\mathbf{h}+F\left(\mathbf{u}% \right)\ \ \ \text{in }\Omega\\ \mathrm{div\ }\mathbf{u}=0\ \ \ \text{in }\Omega\\ \mathbf{u}=0\ \ \ \text{on }\partial\Omega\end{array}\right.

where $\mathbf{h}\in H^{-1}\left(\Omega\right)^{N}$ and $F:H_{0}^{1}\left(\Omega\right)^{N}\rightarrow H^{-1}\left(\Omega\right)^{N}.$ The problem can be reduced to the fixed point equation in $V,$

(2.2)

\mathbf{u}=S\left(\mathbf{h}+F\left(\mathbf{u}\right)\right)

having in mind that the solution operator $S$ is defined and Lipschitz continuous on the open ball of $H^{-1}\left(\Omega\right)^{N}$ centered at the origin and of radius $1/\left(2\kappa M\right).$

Theorem 2.1.

Assume that $F:H_{0}^{1}\left(\Omega\right)^{N}\rightarrow H^{-1}\left(\Omega\right)^{N},\$ $F\left(0\right)=0\$ and

(2.3)

\text{\ }\left|F\left(\mathbf{u}\right)-F\left(\mathbf{v}\right)\right|_{H^{-1% }}\leq a\left|\mathbf{u}-\mathbf{v}\right|_{0}\ \ \ \left(\left|\mathbf{u}% \right|_{0},\left|\mathbf{v}\right|_{0}<\frac{1}{2\kappa M}\right)

for some constant $a<1.$ Then for each $\mathbf{h}\in H^{-1}\left(\Omega\right)^{N}$ with

(2.4)

\left|\mathbf{h}\right|_{H^{-1}}<\frac{\left(1-a\right)^{2}}{2\kappa M}\,,

equation (2.2) has a unique solution $\mathbf{u}\in V$ such that

(2.5)

\left|\mathbf{u}\right|_{V}<\frac{1-a}{2\kappa M}\,.

Proof.

First, we consider the case $\kappa>0$ . Let $\varepsilon$ be arbitrarily closed to $1/a$ with $1<\varepsilon<1/a,$ such that

\left|\mathbf{h}\right|_{H^{-1}}\leq\left(1-a\right)\frac{\varepsilon-1}{2% \kappa M\varepsilon}.

Denote $R=R(\varepsilon):=\displaystyle\frac{\varepsilon-1}{2\kappa M\varepsilon}.$ Thus $\left|\mathbf{h}\right|_{H^{-1}}\leq\left(1-a\right)R.$ Clearly $R<\frac{1-a}{2\kappa M}<\frac{1}{2\kappa M},$ hence inequality (2.3) holds for all $\mathbf{u},\mathbf{v}$ in the closed ball $B_{R}$ of $\left(V,\left|.\right|_{V}\right)$ centered at the origin and of radius $R .$ In addition, for $\mathbf{u}\in B_{R},$ one has

	$\displaystyle\left\|\mathbf{h}+F\left(\mathbf{u}\right)\right\|_{V^{\prime}}$	$\displaystyle\leq$	$\displaystyle\left\|\mathbf{h}+F\left(\mathbf{u}\right)\right\|_{H^{-1}}\leq% \left\|\mathbf{h}\right\|_{H^{-1}}+a\left\|\mathbf{u}\right\|_{V}$
		$\displaystyle\leq$	$\displaystyle\left\|\mathbf{h}\right\|_{H^{-1}}+aR\leq\left(1-a\right)R+aR=R<% \frac{1}{2\kappa M}.$

Consequently the operator

T\left(\mathbf{u}\right):=S\left(\mathbf{h}+F\left(\mathbf{u}\right)\right)

is well defined in $B_{R}$ and for $\left|\mathbf{u}\right|_{V},\left|\mathbf{v}\right|_{V}\leq R,\$ using (1.9) one has

	$\displaystyle\left\|T\left(\mathbf{u}\right)-T\left(\mathbf{v}\right)\right\|_{V}$	$\displaystyle\leq$	$\displaystyle L\left(R\right)\left\|F\left(\mathbf{u}\right)-F\left(\mathbf{v}% \right)\right\|_{V^{\prime}}\leq L\left(R\right)\left\|F\left(\mathbf{u}\right)-% F\left(\mathbf{v}\right)\right\|_{H^{-1}}$
		$\displaystyle\leq$	$\displaystyle aL\left(R\right)\left\|\mathbf{u}-\mathbf{u}\right\|_{V}.$

Since $\ aL\left(R\right)<1$ (which can be checked easily), we have that $T$ is a contraction on $B_{R}.$ In addition from (1.8) and (2) we have

(2.7)

\left|T\left(\mathbf{u}\right)\right|_{V}\leq\left|\mathbf{h}+F\left(\mathbf{u% }\right)\right|_{H^{-1}}\leq R,

which proves that $T\left(B_{R}\right)\subset B_{R}.$ Thus the Banach contraction principle applies and gives the existence and uniqueness of solution ${\bf u}$ in $B_{R}$ , where $R=R(\varepsilon)=\displaystyle\frac{1-\varepsilon^{-1}}{2\kappa M}$ . Since $R(\varepsilon)\to\frac{1-a}{2\kappa M}$ as $\varepsilon\to\frac{1}{a}$ , we obtain the existence and uniqueness of solution ${\bf u}$ satisfying (2.5). This closes the proof in the case $\kappa>0$ .

In the case $\kappa=0$ , the operator $T$ is a contraction on the entire space $V$ and hence the solution exists and is unique in $V$ . ∎

For the next result, instead of the Lipschitz continuity, we only assume a linear growth of $F .$

Theorem 2.2.

Assume that $F:H_{0}^{1}\left(\Omega\right)^{N}\rightarrow H^{-1}\left(\Omega\right)^{N}$ is completely continuous, satisfies the growth condition

(2.8)

\left|F\left(\mathbf{u}\right)\right|_{H^{-1}}\leq a\left|\mathbf{u}\right|_{V% }\ \ \ \left(\left|\mathbf{u}\right|_{0}<\frac{1}{2\kappa M}\right)

with some constant $a>0$ and $\mathbf{h}$ is as in (2.4). Then equation (2.2) has at least one solution $\mathbf{u}\in V$ satisfying (2.5).

Proof.

As above, $T$ maps the ball $B_{R}$ of $\left(V,\left|\cdot\right|_{V}\right)$ into itself. In addition, the operator $T$ is completely continuous. The result follows from Schauder’s fixed point theorem. ∎

A better result can be derived from the Leray-Schauder fixed point theorem, without the linear growth condition on $F$ .

Theorem 2.3.

Assume that the operator $F:H_{0}^{1}\left(\Omega\right)^{N}\rightarrow H^{-1}\left(\Omega\right)^{N}$ is completely continuous and there exists $c<1$ such that

(2.9)

\left(F\left(\mathbf{u}\right),\mathbf{u}\right)\leq c\left|\mathbf{u}\right|_% {0}^{2}\ \ \ \left(\mathbf{u}\in H_{0}^{1}\left(\Omega\right)^{N}\right).

Then for each $R>0$ satisfying

\left(1-c\right)R+\sigma<\frac{1}{2\kappa M},

where $\sigma=\sup_{\mathbf{u}\in B_{R}}\left|F\left(\mathbf{u}\right)\right|_{H^{-1}}$ and $B_{R}=\left\{\mathbf{u}\in V:\ \left|\mathbf{u}\right|_{V}\leq R\right\},$ equation (2.2) has at least one solution $\mathbf{u}\in B_{R}$ for every $\mathbf{h}\in H^{-1}\left(\Omega\right)^{N}$ with

(2.10)

\left|\mathbf{h}\right|_{H^{-1}}\leq\left(1-c\right)R.

Moreover, any solution $\mathbf{u}$ satisfies

(2.11)

\left|\mathbf{u}\right|_{V}\leq\frac{1}{1-c}\left|\mathbf{h}\right|_{H^{-1}}.

Proof.

For $\mathbf{u}\in B_{R},$ one has

\left|\mathbf{h}+F\left(\mathbf{u}\right)\right|_{H^{-1}}\leq\left|\mathbf{h}% \right|_{H^{-1}}+\sigma\leq\left(1-c\right)R+\sigma<\frac{1}{2\kappa M}\,.

Hence $T$ is well-defined and continuous on the closed ball $B_{R}$ of $V .$ Moreover, the Lipschitz continuity of the solution operator $S$ implies that $S$ is a bounded operator (it maps bounded sets into bounded sets) which together with the complete continuity of $F$ implies that the operator $T$ is also completely continuous on $B_{R}.$ (Recall that $T\left(\mathbf{u}\right):=S\left(\mathbf{h}+F\left(\mathbf{u}\right)\right)$ .)

Next we show that the Leray-Schauder condition holds, namely

\left|\mathbf{u}\right|_{V}=R\ \Rightarrow\ T\left(\mathbf{u}\right)\neq% \lambda\mathbf{u}\ \text{ for all }\lambda>1.

Assume the contrary, i.e. $T\left(\mathbf{u}\right)=\lambda\mathbf{u}$ for some $\mathbf{u}$ with $\left|\mathbf{u}\right|_{V}=R$ and $\lambda>1.$ Then as in (1.4) we obtain the variational equation

\lambda\mu\left(\mathbf{u},\mathbf{v}\right)_{H_{0}^{1}}+\lambda\eta\left(% \mathbf{u},\mathbf{v}\right)_{L^{2}}+\lambda^{2}\kappa b\left(\mathbf{u},% \mathbf{u},\mathbf{v}\right)=\left(\mathbf{h}+F\left(\mathbf{u}\right),\mathbf% {v}\right)\,,\ \ \ {\bf v}\in V\,.

Choosing ${\bf v}=\mathbf{u}$ in the above equation and using inequalities (1.7) and (2.9) one obtains

	$\displaystyle\left\|\mathbf{u}\right\|_{V}^{2}<\lambda\left\|\mathbf{u}\right\|_{V% }^{2}$	$\displaystyle=\left(\mathbf{h}+F\left(\mathbf{u}\right),\mathbf{u}\right)$
		$\displaystyle\leq\left\|\mathbf{h}\right\|_{V^{\prime}}\left\|\mathbf{u}\right\|_{% V}+c\left\|\mathbf{u}\right\|_{V}^{2}\leq\left\|\mathbf{h}\right\|_{H^{-1}}\left\|% \mathbf{u}\right\|_{V}+c\left\|\mathbf{u}\right\|_{V}^{2}\,.$

Since $\left|\mathbf{u}\right|_{V}=R,$ it follows that

R<\left|\mathbf{h}\right|_{H^{-1}}+cR\leq\left(1-c\right)R+cR=R,

a contradiction. Thus, the Leray-Schauder fixed point theorem guarantees the existence of a solution.

To prove (2.11), assume that $\mathbf{u}\in V$ is a solution of equation (2.2). Then, as above, we find

\left|\mathbf{u}\right|_{V}\leq\left|\mathbf{h}\right|_{H^{-1}}+c\left|\mathbf% {u}\right|_{V},

whence the conclusion. ∎

Remark 2.1.

(a):

If $F:L^{2}\left(\Omega\right)^{N}\rightarrow L^{2}\left(\Omega\right)^{N}$ and there exists a constant $a_{0}>0$ such that

\ \left|F\left(\mathbf{u}\right)-F\left(\mathbf{v}\right)\right|_{L^{2}}\leq a% _{0}\left|\mathbf{u}-\mathbf{v}\right|_{L^{2}}\ \ \left(\mathbf{u},\mathbf{v}% \in L^{2}\left(\Omega\right)^{N}\right),

then using twice Poincaré’s inequality gives

	$\displaystyle\left\|F\left(\mathbf{u}\right)-F\left(\mathbf{v}\right)\right\|_{H% ^{-1}}$	$\displaystyle\leq$	$\displaystyle\frac{1}{\sqrt{\lambda_{1}}}\left\|F\left(\mathbf{u}\right)-F\left% (\mathbf{V}\right)\right\|_{L^{2}}\leq\frac{a_{0}}{\sqrt{\lambda_{1}}}\left\|% \mathbf{u}-\mathbf{v}\right\|_{L^{2}}$
		$\displaystyle\leq$	$\displaystyle\frac{a_{0}}{\sqrt{\lambda_{1}\left(\mu\lambda_{1}+\eta\right)}}% \left\|\mathbf{u}-\mathbf{v}\right\|_{V}$

for all $\mathbf{u},\mathbf{v}\in V$ and so condition (2.3) holds with $a=a_{0}/\sqrt{\lambda_{1}\left(\mu\lambda_{1}+\eta\right)}.$

(b):

Analogously, if $F:L^{2}\left(\Omega\right)^{N}\rightarrow L^{2}\left(\Omega\right)^{N}$ and

\left|F\left(\mathbf{u}\right)\right|_{L^{2}}\leq a_{0}\left|\mathbf{u}\right|% _{L^{2}}\ \ \ \ \left(\mathbf{u}\in L^{2}\left(\Omega\right)^{N}\right),

then condition (2.8) holds with $a=a_{0}/\sqrt{\lambda_{1}\left(\mu\lambda_{1}+\eta\right)}.$

(c):

If $F:L^{p}\left(\Omega\right)^{N}\rightarrow L^{q}\left(\Omega\right)^{N}$ is continuous and either $p<2^{\ast},$ $q\geq\left(2^{\ast}\right)^{\prime},$ or $p\leq 2^{\ast},\ q>\left(2^{\ast}\right)^{\prime},$ then $F$ is completely continuous from $H_{0}^{1}\left(\Omega\right)^{N}$ to $H^{-1}\left(\Omega\right)^{N}.$ This follows from the continuous (compact) embeddings $H_{0}^{1}\left(\Omega\right)\subset L^{p}\left(\Omega\right)$ for $p\leq 2^{\ast}\ \left(p<2^{\ast}\right),$ and $\ L^{q}\left(\Omega\right)\subset H^{-1}\left(\Omega\right)$ for $q\geq\left(2^{\ast}\right)^{\prime}$ $\left(q>\left(2^{\ast}\right)^{\prime}\right).$ In particular, if $F:L^{2}\left(\Omega\right)^{N}\rightarrow L^{2}\left(\Omega\right)^{N}$ is continuous, then it is completely continuous from $H_{0}^{1}\left(\Omega\right)^{N}$ to $H^{-1}\left(\Omega\right)^{N}.$

Recall that for $N>2$ , $2^{*}=\frac{2N}{N-2}$ , while for $N=2$ , $2^{*}=+\infty$ . In addition, $\frac{1}{2^{*}}+\frac{1}{\left(2^{*}\right)^{\prime}}=1$ .

3. A coupled system of two Navier-Stokes type equations

Consider now the system

(3.1)

\left\{\begin{array}[]{l}-\mu_{i}\Delta\mathbf{u}_{i}+\eta_{i}\mathbf{u}_{i}+% \kappa_{i}\left(\mathbf{u}_{i}\cdot\nabla\right)\mathbf{u}_{i}+\nabla p_{i}=\ % \mathbf{h}_{i}+F_{i}\left(\mathbf{u}\right)\ \ \text{in }\Omega\\ \mathrm{div\ }\mathbf{u}_{i}=0\ \ \ \text{in }\Omega\\ \mathbf{u}_{i}=0\ \ \ \text{on }\partial\Omega\ \ \ \left(i=1,2\right)\end{% array}\right.

where $\mathbf{u}$ stands for the pair $\left(\mathbf{u}_{1},\mathbf{u}_{2}\right),\ \mathbf{h}_{i}\in H^{-1}\left(% \Omega\right)^{N}$ and $F_{i}:H_{0}^{1}\left(\Omega\right)^{2N}\rightarrow H^{-1}\left(\Omega\right)^{% N}\ \left(i=1,2\right).$ The problem can be reduced to the fixed point equation in $V^{2}$

(3.2)

\mathbf{u}_{i}=S_{i}\left(\mathbf{h}_{i}+F_{i}\left(\mathbf{u}\right)\right),% \ \ i=1,2,

where ${\bf u}=\left({\bf u}_{1},{\bf u}_{2}\right)\in V^{2}$ , and $S_{i}$ stands for the solution operator corresponding to ${\bf u}_{i}\in V$ . (Here $S_{i}$ is the solution operator corresponding to system (3.1) with fixed $i$ , compare also with (2.1).)

Denote by $\left(\cdot,\cdot\right)_{i}$ and $\left|\cdot\right|_{i}$ the inner product $\left(\cdot,\cdot\right)_{0}$ and norm $\left|\cdot\right|_{0}$ corresponding to $\mu_{i}$ and $\eta_{i}.$ Also denote by $M_{i}$ the constant $M$ with respect to the norm $\left|\cdot\right|_{i}$ of $V .$

Assume that the following conditions hold:

(H1):

for $i=1,2,$ one has $\ F_{i}:H_{0}^{1}\left(\Omega\right)^{2N}\rightarrow H^{-1}\left(\Omega\right)% ^{N},\$ $F_{i}\left(0\right)=0\$ and

\left|F_{i}\left(\mathbf{u}\right)-F_{i}\left(\mathbf{v}\right)\right|_{H^{-1}% }\leq a_{i1}\left|\mathbf{u}_{1}-\mathbf{v}_{1}\right|_{1}+a_{i2}\left|\mathbf% {u}_{2}-\mathbf{v}_{2}\right|_{2}\ \ \ \left(\left|\mathbf{u}_{j}\right|_{j},% \left|\mathbf{v}_{j}\right|_{j}<\theta_{j},\ j=1,2\right),

where $a_{i1},a_{i2}$ are nonnegative constants such that $a_{i1}+a_{i2}<1$ , and $\theta_{j}:=1/\left(2\kappa_{j}M_{j}\right)$ .

(H2):

there exist $\varepsilon_{i}>1$ , $i=1,2,$ such that

(3.3)

\varepsilon_{1}a_{11}+\varepsilon_{2}a_{22}<\min\left\{2,\ 1+\varepsilon_{1}% \varepsilon_{2}\left(a_{11}a_{22}-a_{12}a_{21}\right)\right\}.

(H3):

$\mathbf{h}_{i}\in H^{-1}\left(\Omega\right)^{N}$ and

\left|\mathbf{h}_{i}\right|_{H^{-1}}\leq\left(1-a_{i1}-a_{i2}\right)R\ \ \ % \left(i=1,2\right),

where

(3.4)

R:=\min\left\{R_{1},\ R_{2}\right\},\ \ \ R_{i}:=\frac{\varepsilon_{i}-1}{2% \kappa_{i}\varepsilon_{i}M_{i}}.

The next result extends Theorem 2.1 to system (3.1).

Theorem 3.1.

Under assumptions (H1)-(H3), system (3.2) has a unique solution $\mathbf{u}=\left(\mathbf{u}_{1},\mathbf{u}_{2}\right)$ $\in V^{2}$ such that

\left|\mathbf{u}_{i}\right|_{i}\leq R\ \ \ \left(i=1,2\right).

Proof.

First note that the operator $T_{i}\left(\mathbf{u}\right):=S_{i}\left(\mathbf{h}_{i}+F_{i}\left(\mathbf{u}% \right)\right)$ is well-defined on

(3.5)

D_{R}:=\left\{\mathbf{u}=\left(\mathbf{u}_{1},\mathbf{u}_{2}\right)\in V^{2}:% \ \left|\mathbf{u}_{i}\right|_{i}\leq R,\ \ i=1,2\right\}.

Indeed, using (H1) one has

$\displaystyle\left\|\mathbf{h}_{i}+F_{i}\left(\mathbf{u}\right)\right\|_{H^{-1}}$	$\displaystyle\leq$	$\displaystyle\left\|\mathbf{h}_{i}\right\|_{H^{-1}}+\left\|F_{i}\left(\mathbf{u}% \right)\right\|_{H^{-1}}$
	$\displaystyle\leq$	$\displaystyle\left\|\mathbf{h}_{i}\right\|_{H^{-1}}+a_{i1}\left\|\mathbf{u}_{1}% \right\|_{1}+a_{i2}\left\|\mathbf{u}_{2}\right\|_{2}$
	$\displaystyle\leq$	$\displaystyle\left(1-a_{i1}-a_{i2}\right)R+a_{i1}R+a_{i2}R$
	$\displaystyle=$	$\displaystyle R\leq R_{i}=\frac{\varepsilon_{i}-1}{2\kappa_{i}\varepsilon_{i}M% _{i}}<\frac{1}{2\kappa_{i}M_{i}}.$

Also if $\ T=\left(T_{1},T_{2}\right),$ then $\ T\left(D_{R}\right)\subset D_{R}.$ Indeed, if $\mathbf{u}\in D_{R},$ then according to (1.8) and the last estimate,

\left|T_{i}\left(\mathbf{u}\right)\right|_{i}=\left|S_{i}\left(\mathbf{h}_{i}+% F_{i}\left(\mathbf{u}\right)\right)\right|_{i}\leq\left|\mathbf{h}_{i}+F_{i}% \left(\mathbf{u}\right)\right|_{H^{-1}}\leq R,

that is $T\left(\mathbf{u}\right)\in D_{R}.$ Next, again from (H1), we have

\left|T_{i}\left(\mathbf{u}\right)-T_{i}\left(\mathbf{v}\right)\right|_{i}\leq L% _{i}\left(R_{i}\right)\left(a_{i1}\left|\mathbf{u}_{1}-\mathbf{v}_{1}\right|_{% 1}+a_{i2}\left|\mathbf{u}_{2}-\mathbf{v}_{2}\right|_{2}\right)\,,\ \ \ i=1,2,

where

L_{i}(R_{i})=\frac{1}{1-2\kappa_{i}M_{i}R_{i}}\,,\ \ \ i=1,2\,.

These inequalities can be put under the vector form

\left[\begin{array}[]{c}\left|T_{1}\left(\mathbf{u}\right)-T_{1}\left(\mathbf{% v}\right)\right|_{1}\\ \left|T_{2}\left(\mathbf{u}\right)-T_{2}\left(\mathbf{v}\right)\right|_{2}\end% {array}\right]\leq\mathcal{M}\left[\begin{array}[]{c}\left|\mathbf{u}_{1}-% \mathbf{v}_{1}\right|_{1}\\ \left|\mathbf{u}_{2}-\mathbf{v}_{2}\right|_{2}\end{array}\right]

in terms of the matrix

\mathcal{M}=\left[L_{i}\left(R_{i}\right)a_{ij}\right]_{i,j=1,2}.

Here, a direct computation based on the expression of $\ R_{i}$ given by (3.4) shows that

L_{i}\left(R_{i}\right)=\varepsilon_{i},

and hence $\mathcal{M}=\left[\varepsilon_{i}a_{ij}\right]_{i,j=1,2}.$ Moreover, inequality (3.3), in view of (1.10), implies that $\mathcal{M}$ is a matrix that converges to zero, equivalently, whose spectral radius is less than one. Therefore Perov’s fixed point theorem applies to $T$ and guarantees that $T$ has in $D_{R}$ a unique fixed point.

Again we recall that if $\kappa_{i}=0$ for some $i\in\{1,2\}$ , one has that $R_{i}=+\infty$ . Moreover, if $\kappa_{1}=\kappa_{2}=0$ , then $D_{R}=V^{2}$ . ∎

Remark 3.1.

$(a)$ A sufficient condition for (H2) to hold is that the matrix $\mathcal{A}:=\left[a_{ij}\right]_{i,j=1,2}$ is convergent to zero, i.e.

(3.6)

a_{11}+a_{22}<\min\left\{2,\ 1+a_{11}a_{22}-a_{12}a_{21}\right\}.

Indeed, the strict inequality in (3.6) remains true if we insert, as condition (3.3) requires, the coefficients $\varepsilon_{1},\,\varepsilon_{2}>1$ sufficiently closed to $1$ .

$(b)$ A sufficient condition for (3.6) to hold is that $a_{i1}+a_{i2}<1$ for $i=1,2.$ Indeed, summing up the two inequalities gives $\sum\limits_{i,j=1}^{2}a_{ij}<2,$ whence one has even more $a_{11}+a_{22}<2.$ In addition, the inequality $a_{11}+a_{22}<1+a_{11}a_{22}-a_{12}a_{21}$ also holds in view of its equivalent form $a_{12}a_{21}<\left(1-a_{11}\right)\left(1-a_{22}\right)$ and of the assumptions $a_{12}<1-a_{11}$ and $a_{21}<1-a_{22}$ .

Remark 3.2.

In the case of the Stokes type equations, when $\kappa_{1}=\kappa_{2}=0$ , under condition (H1) with $a_{i1}+a_{i2}<1\$ for $i=1,2,$ problem (3.1) has a unique solution for every $\left(\mathbf{h}_{1},\mathbf{h}_{2}\right)\in H^{-1}\left(\Omega\right)^{2N}.$

If instead of the Lipschitz condition (H1) and of condition (H3) we consider the more relaxed hypotheses:

(H1^∗):

for $i=1,2,$ $F_{i}:H_{0}^{1}\left(\Omega\right)^{2N}\rightarrow H^{-1}\left(\Omega\right)^{N}$ is completely continuous and

\left|F_{i}\left(\mathbf{u}\right)\right|_{H^{-1}}\leq a_{i1}\left|\mathbf{u}_% {1}\right|_{1}+a_{i2}\left|\mathbf{u}_{2}\right|_{2}\ \ \ \left(\left|u_{j}% \right|_{j}<\theta_{j},\ j=1,2\right);

(H3^∗):

$1>a_{i1}+a_{i2},$ $\ \mathbf{h}_{i}\in H^{-1}\left(\Omega\right)^{N}$ and

(3.7)

\left|\mathbf{h}_{i}\right|_{H^{-1}}<\left(1-a_{i1}-a_{i2}\right)\theta\ \ \ % \left(i=1,2\right),

where

\theta:=\min\left\{\theta_{1},\ \theta_{2}\right\},

then we can prove the existence of at least one solution.

Theorem 3.2.

Under assumptions (H1^∗) and (H3^∗), system (3.2) has at least one solution $\ u=\left(\mathbf{u}_{1},\mathbf{u}_{2}\right)$ $\in V^{2}\$ such that

\left|\mathbf{u}_{i}\right|_{i}<\theta\ \ \ \left(i=1,2\right).

Proof.

In view of the strict inequalities (3.7), we may chose a number $R$ with $0<R<\theta$ such that

\left|\mathbf{h}_{i}\right|_{H^{-1}}\leq\left(1-a_{i1}-a_{i2}\right)R\ \ \ % \left(i=1,2\right).

Then, as in the proof of Theorem 3.1, the operators $T_{i}$ are well-defined and for $T=\left(T_{1},T_{2}\right)$ and $D_{R}$ given by (3.5) with $R$ as above, one has $T\left(D_{R}\right)\subset D_{R}.$ Here in addition $T$ is completely continuous, hence Schauder’s fixed point theorem applies and gives the conclusion. ∎

The next result is the version for systems of Theorem 2.3.

Theorem 3.3.

Assume that the operators $F_{i}:H_{0}^{1}\left(\Omega\right)^{2N}\rightarrow H^{-1}\left(\Omega\right)^{N}$ $\left(i=1,2\right)$ are completely continuous and there exist $c_{i},\ d_{i}\geq 0$ with $c_{i}+d_{i}<1$ such that that for each $i,$

(3.8)

\left(F_{i}\left(\mathbf{u}\right),\mathbf{u}_{i}\right)\leq c_{i}\left|% \mathbf{u}_{i}\right|_{i}^{2}+d_{i}\left|\mathbf{u}_{1}\right|_{1}\left|% \mathbf{u}_{2}\right|_{2}\ \ \ \left(\mathbf{u}\in H_{0}^{1}\left(\Omega\right% )^{2N}\right).

Then for each $R>0$ satisfying

(3.9)

\sigma_{i}\left(R\right)<\theta_{i}-\left(1-c_{i}-d_{i}\right)R,

where $\ \sigma_{i}\left(R\right):=\sup_{\mathbf{u}\in D_{R}}\left|F_{i}\left(\mathbf% {u}\right)\right|_{H^{-1}},$ system (3.2) has at least one solution $\mathbf{u}\in D_{R}$ for every $\left(\mathbf{h}_{1},\mathbf{h}_{2}\right)\in H^{-1}\left(\Omega\right)^{2N}$ with

(3.10)

\ \left|\mathbf{h}_{i}\right|_{H^{-1}}\leq\left(1-c_{i}-d_{i}\right)R.

Moreover, for any solution $\mathbf{u}=\left(\mathbf{u}_{1},\mathbf{u}_{2}\right)$ one has the matrix estimate

(3.11)

\left[\begin{array}[]{c}\left|\mathbf{u}_{1}\right|_{1}\\ \left|\mathbf{u}_{2}\right|_{2}\end{array}\right]\leq\left(\mathcal{I}-% \mathcal{A}\right)^{-1}\left[\begin{array}[]{c}\left|\mathbf{h}_{1}\right|_{H^% {-1}}\\ \left|\mathbf{h}_{2}\right|_{H^{-1}}\end{array}\right],

where

\mathcal{A}=\left[\begin{array}[]{ll}c_{1}&d_{1}\\ d_{2}&c_{2}\end{array}\right]

and $\mathcal{I}$ is the unit matrix.

Proof.

We follow the same ideas as in the proof of Theorem 2.3. Here the Leray-Schauder theorem applies on the set $\ D_{R}.$

For $\mathbf{u}\in D_{R},$ using (3.9) one has

	$\displaystyle\left\|\mathbf{h}_{i}+F_{i}\left(\mathbf{u}\right)\right\|_{H^{-1}}$	$\displaystyle\leq$	$\displaystyle\left\|\mathbf{h}_{i}\right\|_{H^{-1}}+\left\|F_{i}\left(\mathbf{u}% \right)\right\|_{H^{-1}}\leq\left\|\mathbf{h}_{i}\right\|_{H^{-1}}+\sigma_{i}(R)$
		$\displaystyle\leq$	$\displaystyle\left(1-c_{i}-d_{i}\right)R+\sigma_{i}(R)<\frac{1}{2\kappa_{i}M_{% i}}=\theta_{i}\,.$

Hence $T$ is well-defined and continuous on the closed subset $D_{R}$ of $V^{2}.$ Moreover, the Lipschitz continuity of $S_{i}$ implies that $S_{i}$ is a bounded operator (maps bounded sets into bounded sets) which together with the complete continuity of $F_{i}$ implies that $T=\left(T_{1},T_{2}\right),$ $T_{i}\left(\mathbf{u}\right)=S_{i}\left(\mathbf{h}_{i}+F_{i}\left(\mathbf{u}% \right)\right),$ is completely continuous on $D_{R}.$ Next we show that the Leray-Schauder condition holds.

Assume the contrary, i.e. $T\left(\mathbf{u}\right)=\lambda\mathbf{u}$ for some $\mathbf{u}\in D_{R}$ with $\left|\mathbf{u}_{1}\right|_{1}=R$ or $\left|\mathbf{u}_{2}\right|_{2}=R$ and some $\lambda>1.$ Assume that $\left|\mathbf{u}_{1}\right|_{1}=R.$ Then multiplying by $\mathbf{u}_{1}$ the equation

-\lambda\mu_{1}\Delta\mathbf{u}_{1}+\lambda\eta_{1}\mathbf{u}_{1}+\lambda^{2}% \kappa_{1}\left(\mathbf{u}_{1}\cdot\nabla\right)\mathbf{u}_{1}+\nabla p_{1}=% \mathbf{h}_{1}+F_{1}\left(\mathbf{u}\right),

and using (1.7) gives

(3.12)

\left|\mathbf{u}_{1}\right|_{1}^{2}<\lambda\left|\mathbf{u}_{1}\right|_{1}^{2}% =\left(\mathbf{h}_{1}+F_{1}\left(\mathbf{u}\right),\mathbf{u}_{1}\right)\leq% \left|\mathbf{h}_{1}\right|_{H^{-1}}\left|\mathbf{u}_{1}\right|_{1}+\left(F_{1% }\left(\mathbf{u}\right),\mathbf{u}_{1}\right).

Furthermore,using (3.8) we obtain

\left(F_{1}\left(\mathbf{u}\right),\mathbf{u}_{1}\right)\leq c_{1}\left|% \mathbf{u}_{1}\right|_{1}^{2}+d_{1}\left|\mathbf{u}_{1}\right|_{1}\left|% \mathbf{u}_{2}\right|_{2}\leq\left(c_{1}+d_{1}\right)R^{2}.

Hence

R^{2}<\left|\mathbf{h}_{1}\right|_{H^{-1}}R+\left(c_{1}+d_{1}\right)R^{2}.

Similarly, in case that $\left|\mathbf{u}_{2}\right|_{2}=R$ we derive

R^{2}<\left|\mathbf{h}_{2}\right|_{H^{-1}}R+\left(c_{2}+d_{2}\right)R^{2}.

Hence in view of (3.10), we arrive to the contradiction $R<R.$

To prove (3.11) we start with

\left|\mathbf{u}_{i}\right|_{i}^{2}=\left(\mathbf{h}_{i}+F_{i}\left(\mathbf{u}% \right),\mathbf{u}_{i}\right)\leq\left|\mathbf{h}_{i}\right|_{H^{-1}}\left|% \mathbf{u}_{i}\right|_{i}+c_{i}\left|\mathbf{u}_{i}\right|_{i}^{2}+d_{i}\left|% \mathbf{u}_{1}\right|_{1}\left|\mathbf{u}_{2}\right|_{2}\

which give

	$\displaystyle\left\|\mathbf{u}_{1}\right\|_{1}\leq\left\|\mathbf{h}_{1}\right\|_{H% ^{-1}}+c_{1}\left\|\mathbf{u}_{1}\right\|_{1}+d_{1}\left\|\mathbf{u}_{2}\right\|_{% 2}\,,$
	$\displaystyle\left\|\mathbf{u}_{2}\right\|_{2}\leq\left\|\mathbf{h}_{2}\right\|_{H% ^{-1}}+c_{2}\left\|\mathbf{u}_{2}\right\|_{2}+d_{2}\left\|\mathbf{u}_{1}\right\|_{% 1}\,.$

These can be put together under the matrix form

\left(\mathcal{I}-\mathcal{A}\right)\left[\begin{array}[]{c}\left|\mathbf{u}_{% 1}\right|_{1}\\ \left|\mathbf{u}_{2}\right|_{2}\end{array}\right]\leq\left[\begin{array}[]{c}% \left|\mathbf{h}_{1}\right|_{H^{-1}}\\ \left|\mathbf{h}_{2}\right|_{H^{-1}}\end{array}\right].

From the assumption $c_{i}+d_{i}<1\ \left(i=1,2\right)$ we easily can check that matrix $\mathcal{A}$ is convergent to zero, which guarantees that $\mathcal{I}-\mathcal{A}$ has an inverse whose entries are nonnegative. Thus the multiplication by $\left(\mathcal{I}-\mathcal{A}\right)^{-1}$ does not change inequality and yields (3.11). ∎

Remark 3.3.

Simple computation of $\left(\mathcal{I}-\mathcal{A}\right)^{-1}$ shows that estimate (3.11) means explicitly:

	$\displaystyle\left\|\mathbf{u}_{1}\right\|_{1}$	$\displaystyle\leq$	$\displaystyle\frac{\left(1-c_{2}\right)\left\|\mathbf{h}_{1}\right\|_{H^{-1}}+d_% {1}\left\|\mathbf{h}_{2}\right\|_{H^{-1}}}{\left(1-c_{1}\right)\left(1-c_{2}% \right)-d_{1}d_{2}},$
	$\displaystyle\left\|\mathbf{u}_{2}\right\|_{2}$	$\displaystyle\leq$	$\displaystyle\frac{d_{2}\left\|\mathbf{h}_{1}\right\|_{H^{-1}}+\left(1-c_{1}% \right)\left\|\mathbf{h}_{2}\right\|_{H^{-1}}}{\left(1-c_{1}\right)\left(1-c_{2}% \right)-d_{1}d_{2}}.$

Remark 3.4.

In the case of the Stokes type equations, when $\kappa_{1}=\kappa_{2}=0,$ under condition (3.8) with the coefficients $c_{i},\ d_{i}$ such that the matrix $\mathcal{A}$ is convergent to zero, the problem has a unique solution for every $\left(\mathbf{h}_{1},\mathbf{h}_{2}\right)\in H^{-1}\left(\Omega\right)^{2N}.$ Indeed, in this case the operator $T$ is defined and completely continuous on the whole space $V^{2}$ and (3.11) gives the a priori bounds of solutions.

4. Navier-Stokes type model for fluid flow in bidisperse porous media

Now we come back to the specific model of Navier-Stokes type for bidisperse porous media. Thus we consider the following system

(4.1)

\left\{\begin{array}[]{l}-\mu_{1}\Delta\mathbf{u}_{1}+\eta_{1}\mathbf{u}_{1}+% \kappa_{1}\left(\mathbf{u}_{1}\cdot\nabla\right)\mathbf{u}_{1}+\nabla p_{1}=% \mathbf{h}_{1}-\ \alpha_{1}\left|\mathbf{u}_{1}\right|^{p-1}\mathbf{u}_{1}\ -% \gamma_{1}\left(\mathbf{u}_{1}-\mathbf{u}_{2}\right)\ \ \ \text{in }\Omega\\ -\mu_{2}\Delta\mathbf{u}_{2}+\eta_{2}\mathbf{u}_{2}+\kappa_{2}\left(\mathbf{u}% _{2}\cdot\nabla\right)\mathbf{u}_{2}+\nabla p_{2}=\mathbf{h}_{2}-\ \alpha_{2}% \left|\mathbf{u}_{2}\right|^{p-1}\mathbf{u}_{2}\ -\gamma_{2}\left(\mathbf{u}_{% 2}-\mathbf{u}_{1}\right)\ \ \ \text{in }\Omega\\ \mathrm{div}\text{\ }\mathbf{u}_{i}=0\ \ \ \text{in }\Omega\\ \mathbf{u}_{i}=0\ \ \ \text{on }\partial\Omega\ \ \ \left(i=1,2\right),\end{% array}\right.

where $\ 1\leq p<2^{\ast}-1$ $(N\leq 3)$ . Here for

F_{1}\left(\mathbf{u}\right)=-\ \alpha_{1}\left|\mathbf{u}_{1}\right|^{p-1}% \mathbf{u}_{1}\ -\gamma_{1}\left(\mathbf{u}_{1}-\mathbf{u}_{2}\right),\ \ \ F_% {2}\left(\mathbf{u}\right)=-\alpha_{2}\left|\mathbf{u}_{2}\right|^{p-1}\mathbf% {u}_{2}-\gamma_{2}\left(\mathbf{u}_{2}-\mathbf{u}_{1}\right),

we have $F_{i}:H_{0}^{1}\left(\Omega\right)^{2N}\rightarrow L^{2^{\ast}/p}\left(\Omega% \right)^{N}.$ Indeed, from $H_{0}^{1}\left(\Omega\right)\subset L^{2^{\ast}}\left(\Omega\right)$ we have $\left|\mathbf{u}\right|^{p-1}\mathbf{u}\in L^{2^{\ast}/p}\left(\Omega\right)^{% N}.$

Wishing to prove the existence of solutions we shall guarantee that all the assumptions of Theorem 3.3 can be satisfied.

First, since $p<2^{\ast}-1,$ $2^{\ast}/p>\left(2^{\ast}\right)^{\prime}$ , and so one has $L^{2^{\ast}/p}\left(\Omega\right)\subset H^{-1}\left(\Omega\right)$ compactly. It follows that $F_{i}$ is completely continuous from $H_{0}^{1}\left(\Omega\right)^{2N}$ to $H^{-1}\left(\Omega\right)^{N}.$

Secondly, if for some $R>0$ we take $\ \left|\mathbf{u}_{1}\right|_{1},\ \left|\mathbf{u}_{2}\right|_{2}\leq R,$ then we have

$\displaystyle\left\|F_{1}\left(\mathbf{u}\right)\right\|_{H^{-1}}$	$\displaystyle\leq$	$\displaystyle\alpha_{1}\left\|\left\|\mathbf{u}_{1}\right\|^{p-1}\mathbf{u}_{1}% \right\|_{H^{-1}}+\gamma_{1}\left(\left\|\mathbf{u}_{1}\right\|_{H^{-1}}+\left\|% \mathbf{u}_{2}\right\|_{H^{-1}}\right)$
	$\displaystyle\leq$	$\displaystyle\alpha_{1}c_{1}\left\|\mathbf{u}_{1}\right\|_{L^{2^{\ast}}}^{p}+% \frac{\gamma_{1}}{\sqrt{\lambda_{1}}}\left(\left\|\mathbf{u}_{1}\right\|_{L^{2}}% +\left\|\mathbf{u}_{2}\right\|_{L^{2}}\right)$
	$\displaystyle\leq$	$\displaystyle\alpha_{1}\widetilde{c_{1}}R^{p}+2\gamma_{1}\widetilde{c_{2}}R,$

where $c_{1}$ is the embedding constant of the inclusion $L^{2^{\ast}/p}\left(\Omega\right)\subset H^{-1}\left(\Omega\right)$ and $\widetilde{c_{1}},\ \widetilde{c_{2}}$ take into account the embedding constants for $V\subset L^{2^{\ast}}\left(\Omega\right)^{N}$ and $V\subset L^{2}\left(\Omega\right)^{N},$ respectively. Similarly

\left|F_{2}\left(\mathbf{u}\right)\right|_{H^{-1}}\leq\alpha_{2}\widetilde{c_{% 1}}R^{p}+2\gamma_{2}\widetilde{c_{2}}R.

Therefore

\sigma_{i}\left(R\right)=\sup_{\mathbf{u}\in D}\left|F_{i}\left(\mathbf{u}% \right)\right|_{H^{-1}}\leq\alpha_{i}\widetilde{c_{1}}R^{p}+2\gamma_{i}% \widetilde{c_{2}}R\ \ \ \left(i=1,2\right).

Third, we have

	$\displaystyle\left(F_{1}\left(\mathbf{u}\right),\mathbf{u}_{1}\right)$	$\displaystyle=$	$\displaystyle\left(-\ \alpha_{1}\left\|\mathbf{u}_{1}\right\|^{p-1}\mathbf{u}_{1% }\ -\gamma_{1}\left(\mathbf{u}_{1}-\mathbf{u}_{2}\right),\ \mathbf{u}_{1}% \right)\leq\gamma_{1}\left(\mathbf{u}_{2},\mathbf{u}_{1}\right)$
		$\displaystyle\leq$	$\displaystyle\gamma_{1}\left\|\mathbf{u}_{1}\right\|_{L^{2}}\left\|\mathbf{u}_{2}% \right\|_{L^{2}}\leq C\gamma_{1}\left\|\mathbf{u}_{1}\right\|_{1}\left\|\mathbf{u}% _{2}\right\|_{2}\leq C\gamma_{1}R^{2},$

where

(4.2)

C=\frac{1}{\sqrt{\left(\lambda_{1}\mu_{1}+\eta_{1}\right)\left(\lambda_{1}\mu_% {2}+\eta_{2}\right)}},

and analogously

\left(F_{2}\left(\mathbf{u}\right),\mathbf{u}_{2}\right)\leq C\gamma_{2}R^{2}.

Thus, in this example, one has $\ c_{i}=0,$ $d_{i}=C\gamma_{i}$ , $i=1,2,$ and condition (3.9) holds if

(4.3)

\left(1-C\gamma_{i}\right)R+\alpha_{i}\widetilde{c_{1}}R^{p}+2\gamma_{i}% \widetilde{c_{2}}R<\theta_{i}.

Therefore Theorem 3.3 yields the following existence result.

Theorem 4.1.

Let the coefficients $\mu_{i},\,\alpha_{i},\,\gamma_{i}>0$ and $\kappa_{i},\,\eta_{i}\geq 0$ be given with $\gamma_{i}<1/C$ , $i=1,2$ , where $C$ is the constant given by (4.2). Then for each $R>0$ satisfying (4.3), system (4.1) has at least one solution $\mathbf{u}\in D_{R}$ for every $\left(\mathbf{h}_{1},\mathbf{h}_{2}\right)\in H^{-1}\left(\Omega\right)^{2N}$ with

(4.4)

\ \left|\mathbf{h}_{i}\right|_{H^{-1}}\leq\left(1-C\gamma_{i}\right)R.

Moreover, for any solution $\mathbf{u}=\left(\mathbf{u}_{1},\mathbf{u}_{2}\right)\in V^{2},$ one has

\left|\mathbf{u}_{1}\right|_{1}\leq\frac{\left|\mathbf{h}_{1}\right|_{H^{-1}}+% C\gamma_{1}\left|\mathbf{h}_{2}\right|_{H^{-1}}}{1-C^{2}\gamma_{1}\gamma_{2}},% \ \ \ \ \ \left|\mathbf{u}_{2}\right|_{2}\leq\frac{C\gamma_{2}\left|\mathbf{h}% _{1}\right|_{H^{-1}}+\left|\mathbf{h}_{2}\right|_{H^{-1}}}{1-C^{2}\gamma_{1}% \gamma_{2}}.

Notice that condition (4.3) is fulfilled if $R$ is sufficiently small and that it can be chosen like this, according to (4.4), provided that $\left|\mathbf{h}_{i}\right|_{H^{-1}}$ are small enough. Thus we have

Corollary 4.2.

Let the coefficients $\mu_{i},\,\alpha_{i},\,\gamma_{i}>0$ and $\kappa_{i},\,\eta_{i}\geq 0$ be given with $\gamma_{i}<1/C$ . Then problem (4.1) has solutions for each $\left(\mathbf{h}_{1},\mathbf{h}_{2}\right)\in H^{-1}\left(\Omega\right)^{2N}$ with small $\left|\mathbf{h}_{i}\right|_{H^{-1}}.$

If now $\mathbf{h}_{i}$ are given arbitrarily, then (4.4) is fulfilled by some large enough $R .$ Under this number $R,$ condition (4.3) holds provided that $\theta_{i}$ are large (equivalently, $\kappa_{i}$ are small, that is system (4.1) is close to a Stokes type system). Thus we have

Corollary 4.3.

For every coefficients $\mu_{i},\,\alpha_{i},\ \gamma_{i}>0$ and $\kappa_{i},\,\eta_{i}\geq 0$ with $\gamma_{i}<1/C$ , and every $\left(\mathbf{h}_{1},\mathbf{h}_{2}\right)\in H^{-1}\left(\Omega\right)^{2N},$ problem (4.1) has solutions provided that $\ \kappa_{i}$ are small enough.

We note that in the case of the Stokes equations, when $\kappa_{1}=\kappa_{2}=0,$ the above theorem guarantees the existence of a solution for every $\left(\mathbf{h}_{1},\mathbf{h}_{2}\right)\in H^{-1}\left(\Omega\right)^{2N},$ under the only condition that $\gamma_{i}<1/C$ for $i=1,2.$ Recall that $1/C=\sqrt{\left(\lambda_{1}\mu_{1}+\eta_{1}\right)\left(\lambda_{1}\mu_{2}+% \eta_{2}\right)}.$ However, a better result is the following:

Theorem 4.4.

Let $\kappa_{1}=\kappa_{2}=0.$ Then problem (4.1) has a solution for every $\left(\mathbf{h}_{1},\mathbf{h}_{2}\right)\in H^{-1}\left(\Omega\right)^{2N}$ provided that

(4.5)

\gamma_{1}\gamma_{2}<\left(\lambda_{1}\mu_{1}+\eta_{1}\right)\left(\lambda_{1}% \mu_{2}+\eta_{2}\right).

Proof.

In this case, the operator $T$ being well-defined and completely continuous on the whole space $V^{2},$ it remains to guarantee the a priori boundedness of the solutions of the equations $T\left(\mathbf{u}\right)=\lambda\mathbf{u}$ for $\lambda\geq 1.$ For any such a solution, one has

\left|\mathbf{u}_{i}\right|_{i}^{2}\leq\left|\mathbf{h}_{i}\right|_{H^{-1}}% \left|\mathbf{u}_{i}\right|_{i}+C\gamma_{i}\left|\mathbf{u}_{1}\right|_{1}% \left|\mathbf{u}_{2}\right|_{2},\ \ \ i=1,2.

These can be put under a matrix form

\left[\begin{array}[]{c}\left|\mathbf{u}_{1}\right|_{1}\\ \left|\mathbf{u}_{2}\right|_{2}\end{array}\right]\leq\left[\begin{array}[]{c}% \left|\mathbf{h}_{1}\right|_{H^{-1}}\\ \left|\mathbf{h}_{2}\right|_{H^{-1}}\end{array}\right]+\left[\begin{array}[]{% ll}0&C\gamma_{1}\\ C\gamma_{2}&0\end{array}\right]\left[\begin{array}[]{c}\left|\mathbf{u}_{1}% \right|_{1}\\ \left|\mathbf{u}_{2}\right|_{2}\end{array}\right].

The involved square matrix $\mathcal{A}$ is convergent to zero if $\gamma_{1}\gamma_{2}<1/C^{2},$ which is our assumption (4.5). Thus the above matrix inequality is equivalent to

\left[\begin{array}[]{c}\left|\mathbf{u}_{1}\right|_{1}\\ \left|\mathbf{u}_{2}\right|_{2}\end{array}\right]\leq\left(\mathcal{I}-% \mathcal{A}\right)^{-1}\left[\begin{array}[]{c}\left|\mathbf{h}_{1}\right|_{H^% {-1}}\\ \left|\mathbf{h}_{2}\right|_{H^{-1}}\end{array}\right]

which gives the desired a priori bounds. The result follows from the Leray-Schauder fixed point theorem. ∎

Remark 4.1.

If in (4.1), $\alpha_{1}=\alpha_{2}=0$ , under some suitable conditions on $\mathbf{h}_{1}$ and $\mathbf{h}_{2}$ , one can obtain an existence and uniqueness result by using Banach contraction principle and a similar argument as in [19].

Finally we note that our analysis can be developed in order to treat some models of fluid flows in tridisperse porous media. Moreover, the analysis can be adapted to treat some models with variable coefficients in the anisotropic case. For a numerical approach related to such models we refer the reader to the paper [14] and the references therein.

Acknowledgements

M. Kohr acknowledges the support of the grant PN-III-P4-PCE-2021-0993 (cod PCE 69/2022), UEFSCDI, Romania.

Data Availability Statements

Data sharing not applicable to this article as no datasets were generated or analysed during the current study.

Conflict of interest

On behalf of all authors, the corresponding author states that there is no conflict of interest

References

[1] C. Amrouche, M.A. Rodríguez-Bellido, The Oseen and Navier-Stokes equations in a non-solenoidal framework. Math. Meth. Appl. Sci. 39 (2016), 5066–5090.
[2] C. Amrouche, P.G. Ciarlet, C. Mardare, On a Lemma of Jacques-Louis Lions and its relation to other fundamental results. J. Math. Pures Appl. 104 (2015), 207–226.
[3] C. Amrouche, V. Girault, Decomposition of vector spaces and application to the Stokes problem in arbitrary dimension. Czechoslovak Math. J. 44 (1994), 109–140.
[4] F. Boyer, P. Fabrie, Mathematical Tools for the Study of the Incompressible Navier-Stokes Equations and Related Models. Springer, New York, 2013.
[5] C.E. Brennen, Fundamentals of Multiphase Flows. Cambridge University Press, 2005.
[6] M. Bulí $\check{c}$ ek, J. Málek, J. $\check{Z}$ abenský, On generalized Stokes’ and Brinkman’s equations with a pressure- and shear-dependent viscosity and drag coefficient. Nonlinear Analysis: Real World Applications, 26 (2015), 109–132.
[7] F. Capone, M. Gentile, G. Massa, The onset of thermal convection in anisotropic and rotating bidisperse porous media. Z. Angew. Math. Phys. 72:169 (2021).
[8] Z.Q. Chen, P. Cheng, C.T. Hsu, A theoretical and experimental study on stagnant thermal conductivity of bi-dispersed porous media. Int. Commun. Heat Mass Transfer. 27 (2000), 601–610.
[9] P. Constantin, C. Foias, Navier-Stokes Equations. The University of Chicago Press, Chicago, London, 1988.
[10] M. Costabel, Boundary integral operators on Lipschitz domains: Elementary results. SIAM J. Math. Anal. 19 (1988), 613–626.
[11] M. Dindos̆, M. Mitrea, The stationary Navier-Stokes system in nonsmooth manifolds: The Poisson problem in Lipschitz and $C^{1}$ domains. Arch. Rational Mech. Anal. 174 (2004), 1–47.
[12] E. Fabes, C. Kenig, G. Verchota, The Dirichlet problem for the Stokes system on Lipschitz domains. Duke Math. J. 57 (1988), 769–793.
[13] G.P. Galdi, An Introduction to the Mathematical Theory of the Navier–Stokes Equations. Steady-State Problems. Second Edition, Springer, New York, 2011.
[14] M. Ghalambaz, H. Hendizadeh, H. Zargartalebi, I. Pop, Free convection in a square cavity filled with a tridisperse porous medium. Transp. Porous Media 116 (2017), 379–392.
[15] V. Girault, P.A. Raviart, Finite Element Methods for Navier–Stokes Equations, Springer, Berlin, 1986.
[16] G.C. Hsiao, W.L. Wendland, Boundary Integral Equations. Springer-Verlag, Heidelberg, First Edition 2008, Second Edition 2021.
[17] M. Kohr, M. Lanza de Cristoforis, W.L. Wendland, Nonlinear Neumann-transmission problems for Stokes and Brinkman equations on Euclidean Lipschitz domains. Potential Anal. 38 (2013), 1123–1171.
[18] M. Kohr, S.E. Mikhailov, W.L. Wendland, Non-homogeneous Dirichlet-transmission problems for the anisotropic Stokes and Navier-Stokes systems in Lipschitz domains with transversal interfaces. Calc. Var. Partial Differ. Equ. 61:198 (2022), 47 pp.
[19] M. Kohr, S.E. Mikhailov, W.L. Wendland, On some mixed-transmission problems for the anisotropic Stokes and Navier-Stokes systems in Lipschitz domains with transversal interfaces. J. Math. Anal. Appl. 516 (2022), 126464.
[20] M. Kohr, W.L. Wendland, Variational approach for the Stokes and Navier-Stokes systems with nonsmooth coefficients in Lipschitz domains on compact Riemannian manifolds. Calc. Var. Partial Differ. Equ. 57:165 (2018), 1–41.
[21] M. Kohr, W. L. Wendland, Boundary value problems for the Brinkman system with $L^{\infty}$ coefficients in Lipschitz domains on compact Riemannian manifolds. A variational approach. J. Math. Pures Appl. 131 (2019), 17–63.
[22] M.V. Korobkov, K. Pileckas, R. Russo, On the flux problem in the theory of steady Navier-Stokes equations with non-homogeneous boundary conditions. Arch. Rational Mech. Anal. 207 (2013), 185–213.
[23] A.V. Kuznetsov, D.A. Nield, Thermally developing forced convection in a bidisperse porous medium. J. Porous Media 9 (2006), 393–402.
[24] P.G. Lemarié-Rieusset, The Navier-Stokes Problem in the 21st Century. CRC Press, Boca Raton, London, New York, 2016.
[25] G. Łukaszewicz, P. Kalita, Navier-Stokes Equations. An Introduction with Applications. Advances in Mechanics and Mathematics, 34. Springer, Cham, 2016.
[26] A.L. Mazzucato, V. Nistor, Well-posedness and regularity for the elasticity equation with mixed boundary conditions on polyhedral domains and domains with cracks. Arch. Rational Mech. Anal. 195 (2010), 25–73.
[27] M. Mitrea, M. Wright, Boundary value problems for the Stokes system in arbitrary Lipschitz domains. Astérisque. 344 (2012), viii+241 pp.
[28] D.A. Nield, A note on the modelling of bidisperse porous media. Transp. Porous. Media 111 (2016), 517–520.
[29] D.A. Nield, A. Bejan, Convection in Porous Media. Third Edition, Springer, New York, 2013.
[30] D.A. Nield, A.V. Kuznetsov, A two-velocity two-temperature model for a bi-dispersed porous medium: forced convection in a channel. Transp. Porous Media 59 (2005), 325–339.
[31] D.A. Nield, A.V. Kuznetsov, Heat transfer in bidisperse porous media. In: D.B. Ingham, I. Pop (Eds.), Transport Phenomena in Porous Media III, Elsevier, Oxford, 2005, pp. 34–59.
[32] D.A. Nield, A.V. Kuznetsov, The onset of convection in a bidisperse porous medium. Int. J. Heat Mass Transfer 49 (2006), 3068–3074.
[33] D.A. Nield, A.V. Kuznetsov, A note on modeling high speed flow in a bidisperse porous medium. Transp. Porous Med. 96 (2013), 495–499.
[34] D. O’Regan, R. Precup, Theorems of Leray-Schauder Type and Applications. Gordon and Breach, Amsterdam, 2001.
[35] R. Precup, Existence results for nonlinear boundary value problems under nonresonance conditions. In: Qualitative Problems for Differential Equations and Control Theory, C. Corduneanu (Ed.), World Scientific, Singapore, 1995, pp. 263–273.
[36] R. Precup, Methods in Nonlinear Integral Equations. Kluwer, Dordrecht, 2002.
[37] R. Precup, Linear and Semilinear Partial Differential Equations. De Gruyter, Berlin, 2013.
[38] R. Precup, P. Rubbioni, Stationary solutions of Fokker-Planck equations with nonlinear reaction terms in bounded domains. Potential Anal. 57 (2022), 181–199.
[39] F.O. Pătrulescu, T. Groşsan, I. Pop, Natural convection from a vertical plate embedded in a non-Darcy bidisperse porous medium, J. Heat Transfer 142 (2020): 012504.
[40] C. Revnic, T. Groşan, I. Pop, D.B. Ingham, Free convection in a square cavity filled with a bidisperse porous medium, Int. J. Thermal Sci. 48 (2009), 1876–1883.
[41] G. Seregin, Lecture Notes on Regularity Theory for the Navier-Stokes Equations. World Scientific, London, 2015.
[42] H. Sohr, The Navier-Stokes Equations: An Elementary Functional Analytic Approach. Birkhäuser Verlag, Basel 2001.
[43] B. Straughan, Bidispersive Porous Media. In: Convection with Local Thermal Non-Equilibrium and Microfluidic Effects. Advances in Mechanics and Mathematics, vol. 32, 2015. Springer, Cham.
[44] B. Straughan, Anisotropic bidispersive convection. Proc. R. Soc. A 475:20190206, 2019.
[45] R. Temam, Navier-Stokes Equations. Theory and Numerical Analysis. AMS Chelsea Edition, American Mathematical Society, 2001.
[46] W. Varnhorn, The Stokes Equations. Akademie Verlag, Berlin, 1994.

[1] Amrouche, C., Rodríguez-Bellido, M.A.: The Oseen and Navier–Stokes equations in a non-solenoidal framework. Math. Methods Appl. Sci. 39, 5066–5090 (2016) ADS MathSciNet MATH Google Scholar
[2] Amrouche, C., Ciarlet, P.G., Mardare, C.: On a Lemma of Jacques–Louis Lions and its relation to other fundamental results. J. Math. Pures Appl. 104, 207–226 (2015) MathSciNet MATH Google Scholar
[3] Amrouche, C., Girault, V.: Decomposition of vector spaces and application to the Stokes problem in arbitrary dimension. Czechoslov. Math. J. 44, 109–140 (1994) MathSciNet MATH Google Scholar
[4] Boyer, F., Fabrie, P.: Mathematical Tools for the Study of the Incompressible Navier–Stokes Equations and Related Models. Springer, New York (2013) MATH Google Scholar
[5] Brennen, C.E.: Fundamentals of Multiphase Flows. Cambridge University Press, Cambridge (2005) MATH Google Scholar
[6] Bulíček, M., Málek, J., Žabenský, J.: On generalized Stokes’ and Brinkman’s equations with a pressure- and shear-dependent viscosity and drag coefficient. Nonlinear Anal. Real World Appl. 26, 109–132 (2015) MathSciNet MATH Google Scholar
[7] Capone, F., Gentile, M., Massa, G.: The onset of thermal convection in anisotropic and rotating bidisperse porous media. Z. Angew. Math. Phys. 72, 169 (2021) MathSciNet MATH Google Scholar
[8] Chen, Z.Q., Cheng, P., Hsu, C.T.: A theoretical and experimental study on stagnant thermal conductivity of bi-dispersed porous media. Int. Commun. Heat Mass Transf. 27, 601–610 (2000) Google Scholar
[9] Constantin, P., Foias, C.: Navier–Stokes Equations. The University of Chicago Press, Chicago (1988) MATH Google Scholar
[10] Costabel, M.: Boundary integral operators on Lipschitz domains: elementary results. SIAM J. Math. Anal. 19, 613–626 (1988) MathSciNet MATH Google Scholar
[11] Dindos, M., Mitrea, M.: The stationary Navier–Stokes system in nonsmooth manifolds: the Poisson problem in Lipschitz and C¹ domains. Arch. Ration. Mech. Anal. 174, 1–47 (2004)
[12] Fabes, E., Kenig, C., Verchota, G.: The Dirichlet problem for the Stokes system on Lipschitz domains. Duke Math. J. 57, 769–793 (1988) MathSciNet MATH Google Scholar
[13] Galdi, G.P.: An Introduction to the Mathematical Theory of the Navier–Stokes Equations. Steady-State Problems, 2nd edn. Springer, New York (2011) MATH Google Scholar
[14] Ghalambaz, M., Hendizadeh, H., Zargartalebi, H., Pop, I.: Free convection in a square cavity filled with a tridisperse porous medium. Transp. Porous Media 116, 379–392 (2017) MathSciNet Google Scholar
[15] Girault, V., Raviart, P.A.: Finite Element Methods for Navier–Stokes Equations. Springer, Berlin (1986) MATH Google Scholar
[16] Hsiao, G.C., Wendland, W.L.: Boundary Integral Equations. Springer, Heidelberg, First Edition 2008, Second Edition 2021[17] Kohr, M., Lanza de Cristoforis, M., Wendland, W.L.: Nonlinear Neumann-transmission problems for Stokes and Brinkman equations on Euclidean Lipschitz domains. Potential Anal. 38, 1123–1171 (2013) MathSciNet MATH Google Scholar
[18] Kohr, M., Mikhailov, S.E., Wendland, W.L.: Non-homogeneous Dirichlet-transmission problems for the anisotropic Stokes and Navier–Stokes systems in Lipschitz domains with transversal interfaces. Calc. Var. Partial Differ. Equ. 61, 198 (2022) MathSciNet MATH Google Scholar
[19] Kohr, M., Mikhailov, S.E., Wendland, W.L.: On some mixed-transmission problems for the anisotropic Stokes and Navier–Stokes systems in Lipschitz domains with transversal interfaces. J. Math. Anal. Appl. 516, 126464 (2022) MathSciNet MATH Google Scholar
[20] Kohr, M., Wendland, W.L.: Variational approach for the Stokes and Navier–Stokes systems with nonsmooth coefficients in Lipschitz domains on compact Riemannian manifolds. Calc. Var. Partial Differ. Equ. 57(165), 1–41 (2018) MathSciNet MATH Google Scholar
[21] Kohr, M., Wendland, W.L.: Boundary value problems for the Brinkman system with $L^{\infty}$ coefficients in Lipschitz domains on compact Riemannian manifolds.
A variational approach. J. Math. Pures Appl. 131, 17–63 (2019) MathSciNet
[22] Korobkov, M.V., Pileckas, K., Russo, R.: On the flux problem in the theory of steady Navier–Stokes equations with non-homogeneous boundary conditions. Arch. Ration. Mech. Anal. 207, 185–213 (2013) MathSciNet MATH Google Scholar
[23] Kuznetsov, A.V., Nield, D.A.: Thermally developing forced convection in a bidisperse porous medium. J. Porous Media 9, 393–402 (2006) Google Scholar
[24] Lemarié-Rieusset, P.G.: The Navier–Stokes Problem in the 21st Century. CRC Press, Boca Raton (2016) MATH Google Scholar
[25] Łukaszewicz, G., Kalita, P.: Navier–Stokes Equations. An Introduction with Applications. Advances in Mechanics and Mathematics, vol. 34. Springer, Cham (2016) MATH Google Scholar
[26] Mazzucato, A.L., Nistor, V.: Well-posedness and regularity for the elasticity equation with mixed boundary conditions on polyhedral domains and domains with cracks. Arch. Ration. Mech. Anal. 195, 25–73 (2010) MathSciNet MATH Google Scholar
[27] Mitrea, M., Wright, M.: Boundary value problems for the Stokes system in arbitrary Lipschitz domains. Astérisque. 344, viii+241 pp (2012)
[28] Nield, D.A.: A note on the modelling of bidisperse porous media. Transp. Porous Media 111, 517–520 (2016) MathSciNet Google Scholar
[29] Nield, D.A., Bejan, A.: Convection in Porous Media, 3rd edn. Springer, New York (2013) MATH Google Scholar
[30] Nield, D.A., Kuznetsov, A.V.: A two-velocity two-temperature model for a bi-dispersed porous medium: forced convection in a channel. Transp. Porous Media 59, 325–339 (2005) Google Scholar
[31] Nield, D.A., Kuznetsov, A.V.: Heat transfer in bidisperse porous media. In: Ingham, D.B., Pop, I. (eds.) Transport Phenomena in Porous Media III, pp. 34–59. Elsevier, Oxford (2005) Google Scholar
[32] Nield, D.A., Kuznetsov, A.V.: The onset of convection in a bidisperse porous medium. Int. J. Heat Mass Transf. 49, 3068–3074 (2006) MATH Google Scholar
[33] Nield, D.A., Kuznetsov, A.V.: A note on modeling high speed flow in a bidisperse porous medium. Transp. Porous Media 96, 495–499 (2013) MathSciNet Google Scholar
[34] O’Regan, D., Precup, R.: Theorems of Leray–Schauder Type and Applications. Gordon and Breach, Amsterdam (2001) MATH Google Scholar
[35] Precup, R.: Existence results for nonlinear boundary value problems under nonresonance conditions. In: Corduneanu, C. (ed.) Qualitative Problems for Differential Equations and Control Theory, pp. 263–273. World Scientific, Singapore (1995) MATH Google Scholar
[36] Precup, R.: Methods in Nonlinear Integral Equations. Kluwer, Dordrecht (2002) MATH Google Scholar
[37] Precup, R.: Linear and Semilinear Partial Differential Equations. De Gruyter, Berlin (2013) MATH Google Scholar
[38] Precup, R., Rubbioni, P.: Stationary solutions of Fokker–Planck equations with nonlinear reaction terms in bounded domains. Potential Anal. 57, 181–199 (2022) MathSciNet MATH Google Scholar
[39] Pătrulescu, F.O., Groşsan, T., Pop, I.: Natural convection from a vertical plate embedded in a non-Darcy bidisperse porous medium. J. Heat Transf. 142, 012504 (2020) Google Scholar
[40] Revnic, C., Groşan, T., Pop, I., Ingham, D.B.: Free convection in a square cavity filled with a bidisperse porous medium. Int. J. Therm. Sci. 48, 1876–1883 (2009) MATH Google Scholar
[41] Seregin, G.: Lecture Notes on Regularity Theory for the Navier–Stokes Equations. World Scientific, London (2015) mMATH Google Scholar
[42] Sohr, H.: The Navier–Stokes Equations: An Elementary Functional Analytic Approach. Birkhäuser, Basel (2001) MATH Google Scholar
[43] Straughan, B.: Bidispersive porous media. In: Straughan, B. (ed.) Convection with Local Thermal Non-Equilibrium and Microfluidic Effects. Advances in Mechanics and Mathematics, vol. 32. Springer, Cham (2015) MATH Google Scholar

[44] Straughan, B.: Anisotropic bidispersive convection. Proc. R. Soc. A 475, 20190206 (2019) ADS MathSciNet MATH Google Scholar

[45] Temam, R.: Navier–Stokes Equations. Theory and Numerical Analysis, AMS Chelsea edn. American Mathematical Society (2001)

[46] Varnhorn, W.: The Stokes Equations. Akademie, Berlin (1994) MATH Google Scholar

2023

A method of lower and upper solutions for control problems and application to a model of bone marrow transplantation

October 24, 2023

Componentwise Dinkelbach algorithm for nonlinear fractional optimization problems

September 20, 2023

Linking Methods for Componentwise Variational Systems

October 23, 2023

	$\displaystyle\left\|S\left(\mathbf{f}\right)-S\left(\mathbf{g}\right)\right\|_{V% }^{2}$	$\displaystyle\leq$	$\displaystyle\left\|\mathbf{f}-\mathbf{g}\right\|_{V^{\prime}}\left\|S\left(% \mathbf{f}\right)-S\left(\mathbf{g}\right)\right\|_{V}+\kappa M\left\|S\left(% \mathbf{f}\right)-S\left(\mathbf{g}\right)\right\|_{V}^{2}\left(\left\|S\left(% \mathbf{g}\right)\right\|_{V}+\left\|S\left(\mathbf{f}\right)\right\|_{V}\right)$
		$\displaystyle\leq$	$\displaystyle\left\|\mathbf{f}-\mathbf{g}\right\|_{V^{\prime}}\left\|S\left(% \mathbf{f}\right)-S\left(g\right)\right\|_{V}+\kappa M\left(\left\|{\bf f}\right% \|_{V^{\prime}}+\left\|\mathbf{g}\right\|_{V^{\prime}}\right)\left\|S\left(\mathbf% {f}\right)-S\left(\mathbf{g}\right)\right\|_{V}^{2}.$

	$\displaystyle\left\|\mathbf{h}+F\left(\mathbf{u}\right)\right\|_{V^{\prime}}$	$\displaystyle\leq$	$\displaystyle\left\|\mathbf{h}+F\left(\mathbf{u}\right)\right\|_{H^{-1}}\leq% \left\|\mathbf{h}\right\|_{H^{-1}}+a\left\|\mathbf{u}\right\|_{V}$
		$\displaystyle\leq$	$\displaystyle\left\|\mathbf{h}\right\|_{H^{-1}}+aR\leq\left(1-a\right)R+aR=R<% \frac{1}{2\kappa M}.$

	$\displaystyle\left\|T\left(\mathbf{u}\right)-T\left(\mathbf{v}\right)\right\|_{V}$	$\displaystyle\leq$	$\displaystyle L\left(R\right)\left\|F\left(\mathbf{u}\right)-F\left(\mathbf{v}% \right)\right\|_{V^{\prime}}\leq L\left(R\right)\left\|F\left(\mathbf{u}\right)-% F\left(\mathbf{v}\right)\right\|_{H^{-1}}$
		$\displaystyle\leq$	$\displaystyle aL\left(R\right)\left\|\mathbf{u}-\mathbf{u}\right\|_{V}.$

	$\displaystyle\left\|F\left(\mathbf{u}\right)-F\left(\mathbf{v}\right)\right\|_{H% ^{-1}}$	$\displaystyle\leq$	$\displaystyle\frac{1}{\sqrt{\lambda_{1}}}\left\|F\left(\mathbf{u}\right)-F\left% (\mathbf{V}\right)\right\|_{L^{2}}\leq\frac{a_{0}}{\sqrt{\lambda_{1}}}\left\|% \mathbf{u}-\mathbf{v}\right\|_{L^{2}}$
		$\displaystyle\leq$	$\displaystyle\frac{a_{0}}{\sqrt{\lambda_{1}\left(\mu\lambda_{1}+\eta\right)}}% \left\|\mathbf{u}-\mathbf{v}\right\|_{V}$

	$\displaystyle\left\|\mathbf{u}_{1}\right\|_{1}$	$\displaystyle\leq$	$\displaystyle\frac{\left(1-c_{2}\right)\left\|\mathbf{h}_{1}\right\|_{H^{-1}}+d_% {1}\left\|\mathbf{h}_{2}\right\|_{H^{-1}}}{\left(1-c_{1}\right)\left(1-c_{2}% \right)-d_{1}d_{2}},$
	$\displaystyle\left\|\mathbf{u}_{2}\right\|_{2}$	$\displaystyle\leq$	$\displaystyle\frac{d_{2}\left\|\mathbf{h}_{1}\right\|_{H^{-1}}+\left(1-c_{1}% \right)\left\|\mathbf{h}_{2}\right\|_{H^{-1}}}{\left(1-c_{1}\right)\left(1-c_{2}% \right)-d_{1}d_{2}}.$

Analysis of Navier–Stokes Models for Flows in Bidisperse Porous Media

Abstract

Authors

Keywords

Paper coordinates

PDF

About this paper

Journal

Publisher Name

DOI

Print ISSN

Online ISSN

Paper (preprint) in HTML form

Analysis of Navier-Stokes models for flows in bidisperse porous media

Abstract.

Key words and phrases:

1991 Mathematics Subject Classification:

1. Introduction

1.1. Bidisperse porous media and some related models

1.2. Stationary Navier-Stokes type equations

1.3. Notions of vector analysis in fixed point theory

2. Navier-Stokes type equations with reaction terms

Theorem 2.1.

Proof.

Theorem 2.2.

Proof.

Theorem 2.3.

Proof.

Remark 2.1.

3. A coupled system of two Navier-Stokes type equations

Theorem 3.1.

Proof.

Remark 3.1.

Remark 3.2.

Theorem 3.2.

Proof.

Theorem 3.3.

Proof.

Remark 3.3.

Remark 3.4.

4. Navier-Stokes type model for fluid flow in bidisperse porous media

Theorem 4.1.

Corollary 4.2.

Corollary 4.3.

Theorem 4.4.

Proof.

Remark 4.1.

Acknowledgements

Data Availability Statements

Conflict of interest

References

References

Related Posts