Next Article in Journal
Numerical Solution of Third-Order Rosenau–Hyman and Fornberg–Whitham Equations via B-Spline Interpolation Approach
Previous Article in Journal
A Fuzzy Logic for Semi-Overlap Functions and Their Residua
Previous Article in Special Issue
Fractional Sequential Coupled Systems of Hilfer and Caputo Integro-Differential Equations with Non-Separated Boundary Conditions
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Existence Result for a Class of Time-Fractional Nonstationary Incompressible Navier–Stokes–Voigt Equations

1
College of General Education, Guangxi Vocational University of Agriculture, Nanning 530009, China
2
School of Mathematics and Physics, Center for Applied Mathematics of Guangxi, Guangxi Minzu University, Nanning 530006, China
*
Author to whom correspondence should be addressed.
Axioms 2024, 13(8), 499; https://doi.org/10.3390/axioms13080499
Submission received: 22 May 2024 / Revised: 23 July 2024 / Accepted: 23 July 2024 / Published: 25 July 2024

Abstract

:
We are devoted in this work to dealing with a class of time-fractional nonstationary incompressible Navier–Stokes–Voigt equation involving the Caputo fractional derivative. By exploiting the properties of the operators in the equation, we use the Rothe method to show the existence of weak solutions to the equation by verifying all the conditions of the surjectivity theorem for nonlinear weakly continuous operators.

1. Introduction

The Navier–Stokes (NS) equations are fundamental in fluid dynamics, describing the conservation of momentum in viscous, incompressible fluids. However, for fluids with viscoelastic properties, such as polymer solutions or blood, the NS equations alone may not be sufficient to accurately capture the fluid’s behavior. The Navier–Stokes–Voigt (NSV) equations are a modified form of the NS equations that incorporate the effects of fluid viscoelasticity through the introduction of the Voigt model (see, e.g., [1,2,3,4,5]). The NSV equations extend the NS equations by incorporating the Voigt model, which is a simple representation of viscoelasticity. This model combines the elastic and viscous components of the fluid’s stress tensor, allowing for a more accurate description of fluid flows with viscoelastic characteristics. The NSV equations have several applications in fluid dynamics, particularly in the study of complex fluids and their interactions with solid boundaries. They can be used to simulate flows in biological systems, such as blood flow through blood vessels, or to analyze the behavior of polymer solutions in industrial processes.
The fractional Navier–Stokes–Voigt (FNSV) equations are the generalization of the classical NSV equations that incorporate the concept of fractional calculus, and they are also the natural generalization of the fractional Navier–Stokes (FNS) equations. For more details about the theory of FNS equations, we refer to [6,7,8,9,10,11,12].
There are several aspects of the significance of the FNSV equations. First, by incorporating fractional derivatives, the FNSV equations introduce a degree of memory and nonlocality into the fluid dynamics, enabling it to describe these phenomena more precisely. This enhanced description allows researchers to gain deeper insights into the underlying physics of fluid motion. Second, the study of the FNSV equations not only advance fluid mechanics but also contribute to the development of fractional calculus and nonlinear partial differential equations. By delving into the theoretical underpinnings of these equations, researchers can uncover new mathematical properties and relationships that enrich our understanding of fluid dynamics and its applications. Third, the significance of the FNSV equations transcends individual disciplines, fostering interdisciplinary collaboration among mathematicians, physicists, engineers, and other experts. This collaboration enables the integration of diverse perspectives and expertise, leading to more comprehensive and innovative solutions to complex fluid dynamics problems. Fourth, the equations have broad application potentials, including modeling blood flow, modeling drug delivery, investigating the rheological properties of viscoelastic materials, simulating groundwater flow, and analyzing the impact of complex fluids on aircraft design and performance.
Now, we introduce our problem. Let Ω be a bounded and connected domain in R 3 with Lipschitz continuous boundary Γ = Ω , and Q = Ω × ( 0 , T ) . The following time-fractional Navier–Stokes–Voigt equation discussed in this work has the following form:
D t q 0 C u ν 0 Δ u α 2 Δ ( D t q 0 C u ) + ( u · ) u + p = f , in Q , · u = 0 , in Q , u ( x , t ) = 0 , on Ω × ( 0 , T ) , u ( x , 0 ) = u 0 ( x ) , in Ω .
Here, D t q 0 C u is the Caputo time-fractional derivative of velocity field u, ( u · ) v s . = j = 1 3 u j v i x j i = 1 3 is the nonlinear convective term, p = p ( x , t ) is the unknown pressure occupying Ω , ν 0 > 0 stands for the constant viscosity of the fluid, α is a length-scale parameter characterizing the elasticity of the fluid, and f = f ( x , t ) is a given external volume force of density occupying Ω .
In this work, we first give its variational formulation, which can be transferred to an operator equation, in Section 2 and then present an existence result on solutions to the equation in Section 3. The main novelties of this paper are shown in the following: First, instead of the FNS equations and integral NSV equations widely investigated in the literature (see [1,2,13]), we deal with FNSV equations, which are more general than the previous equations. Second, we use Rothe’s method to prove the existence result based on the surjective theorem of nonlinear weakly continuous operators (see [14,15,16]), which are more general and easier to verify than the ones in previous contributions. Third, the existence result for our problem improves many existence results in previous works, for instance [13], and the main problem can be generalized to a more complex form such as a mixed FNSV equation. Note that the weak solution of the integral NSV equation is unique. However, in this work, we do not investigate the uniqueness of the weak solution of the FNSV equation, which can be our future work.

2. Variational Formulation

Let X be a Banach space with the norm · X , the dual space X , and the duality pairing · , · X . For T > 0 , p > 1 , let L p ( 0 , T ; X ) and C ( [ 0 , T ] ; X ) be the spaces of all p-th integral and continuous functions from [ 0 , T ] into X with usual norms, respectively.
Definition 1
([17,18,19,20]). Let g : [ 0 , ) X be a function. The fractional integral and the Caputo fractional derivative of order q ( 0 , 1 ) for g are defined by, respectively,
I t q 0 g ( t ) = 1 Γ ( q ) 0 t ( t s ) q 1 g ( s ) d s , t > 0 ,
D t q 0 C g ( t ) = 1 Γ ( 1 q ) 0 t ( t s ) q g ( s ) d s = I t 1 q 0 g ( t ) , t > 0 ,
where Γ is the gamma function.
Next, we show the variational formulation of problem (1).
Define the following two inner products of ( L 2 ( Ω ) ) 3 and ( H 0 1 ( Ω ) ) 3 , respectively,
( u , v ) 1 : = Ω j = 1 3 u j v j d s , u = ( u 1 , u 2 , u 3 ) , v = ( v 1 , v 2 , v 3 ) ( L 2 ( Ω ) ) 3
( u , v ) 2 : = Ω j = 1 3 u j · v j d s , u = ( u 1 , u 2 , u 3 ) , v = ( v 1 , v 2 , v 3 ) ( H 0 1 ( Ω ) ) 3 ,
and the corresponding norms v 1 : = ( v , v ) 1 , v 2 : = ( v , v ) 2 . Consider the following:
V 0 = { v ( C 0 ( Ω ) ) 3 : · v = 0 } .
Denote by H the closure of V 0 in ( L 2 ( Ω ) ) 3 and by V the closure of V 0 in ( H 0 1 ( Ω ) ) 3 . Clearly, H and V are Hilbert spaces with scalar products ( · , · ) 1 and ( · , · ) 2 . It is easy to see that V H V and the embeddings are dense, continuous, and compact. For q ( 0 , 1 ) , let p > 1 q and
V = L p ( 0 , T ; V ) , V = L q ( 0 , T ; V ) a n d W = { v V : D t q 0 C v V } .
The dual pairing between V and V is given by
v , v V : = 0 T v ( t ) , v ( t ) V d t , v V , v V .
It is essential to define the trilinear form b : V × V × V R by
b ( u , v , w ) : = i , j 3 Ω u i v j x i w j d x .
Lemma 1
([21,22]). The following properties hold:
(i) 
b ( y , z , w ) = b ( y , w , z ) , y , z , w V ;
(ii) 
b ( y , z , z ) = 0 , y , z V ;
(iii) 
| b ( y , z , w ) | C y 1 1 4 y 2 3 4 z 2 w 1 1 4 w 2 3 4 , C y 2 z 2 w 2 , y , z , w V
and in particular,
| b ( y , z , y ) | C y 1 1 2 y 2 3 2 z 2 , y , z V .
We now give the definition of weak solutions to problem (1).
Definition 2
([1,2]). A function u W is called a weak solution to problem (1) if
( D t q 0 C u ( t ) , v ) 1 + ν 0 ( u ( t ) , v ) 2 + α 2 ( D t q 0 C u ( t ) , v ) 2 + b ( u ( t ) , u ( t ) , v ) = f ( t ) , v V , v V , a . e . t ( 0 , T ) , u ( 0 ) = u 0 .
We next give an equivalent formulation for Definition 2, which is in fact an operator equation. To this aim, we introduce two linear bounded operators, D , E : V V , and a nonlinear operator, F : V V , defined by, respectively,
D u , v V : = ( u , v ) 1 , u , v V ,
E u , v V : = ( u , v ) 2 , u , v V ,
F u , v V : = b ( u , u , v ) , u , v V .
Now, we have an equivalent formulation with Definition 2.
Definition 3.
A function u W is called a weak solution to problem (1) on ( 0 , T ) if
D ( D t q 0 C u ( t ) ) + ν 0 E u ( t ) + α 2 E ( D t q 0 C u ( t ) ) + F u ( t ) = f ( t ) , i n V , u ( 0 ) = u 0 , i n V .
Lemma 2
([2,14,22]). The operators D , E are linear and bounded, and there exist constants c i > 0 ( i = 1 , 2 , , 5 ) such that
D v V c 1 v V , E v V c 2 v V , F v V c 3 v V 2 , v V ,
D v , v V c 4 v V 2 , E v , v V c 5 v V 2 , F v , v V = 0 , v V .

3. Existence of Solutions

We first present the main result of this work.
Theorem 1.
If f C ( [ 0 , T ] ; V ) and u 0 V , then problem (1) has a weak solution.
The proof of Theorem 1 will be divided into several steps by applying the Rothe method (see [23,24]). We shall construct a sequence of lemmas. First, we introduce an auxiliary problem associated with problem (1). Then we establish the existence of a solution to the Rothe problem corresponding to the auxiliary problem by exploiting the surjectivity result for nonlinear weakly continuous operators. Next, the Rothe sequences of piecewise approximations are estimated. Finally, we obtain the existence of solutions to problem (1), which are the limit of the Rothe sequences.
Let w ( t ) = D t q 0 C u ( t ) for a.e. t ( 0 , T ) . Then
u ( t ) = u 0 + I t q 0 w ( t ) for   all t [ 0 , T ] .
It is not hard to see that u W defined by (2) is a weak solution to problem (1) if and only if w V solves the following problem:
Problem 1.
Find w V such that
( D + α 2 E ) w ( t ) + ( ν 0 E + F ) ( u 0 + I t q 0 w ( t ) ) = f ( t ) f o r   a . e . t ( 0 , T ) .
Therefore, in the sequel, it is sufficient to establish the result on the existence of solutions to Problem 1.
Now, we introduce the Rothe method for Problem 1 [23,24,25,26]. Let N N be fixed, τ = T N denote the time step, and t k = k τ for k = 1 , 2 , , N . For any g V , define its approximation and the fractional integral I t q 0 g ( t ) by, respectively,
g τ k = 1 τ t k 1 t k g ( s ) d s for k = 1 , , N ,
I ˜ t k q 0 g ( t k ) = 1 Γ ( q ) j = 1 k t j 1 t j ( t k s ) q 1 g ( t j ) d s = τ q Γ ( q + 1 ) j = 1 k g ( t j ) ( k j + 1 ) q ( k j ) q for k = 1 , , N .
The following problem is called the Rothe problem for Problem 1.
Problem 2.
Find a sequence { w τ k } k = 1 N V such that
( D + α 2 E ) w τ k + ( ν 0 E + F ) v τ k = f τ k
for k = 1 , , N , where
v τ k = u 0 + τ q Γ ( q + 1 ) j = 1 k w τ j ( k j + 1 ) q ( k j ) q .
Lemma 3.
There exists τ 0 > 0 such that for all τ ( 0 , τ 0 ) , Problem 2 has a solution { w τ k } k = 1 N V .
Proof. 
Clearly, Problem 2 is equivalent to the following problem:
find w τ k V such that
Φ ( w τ k ) = f τ k for k = 1 , , N ,
where Φ : V V is given, without loss of generality, by
Φ ( v ) = ( D + α 2 E ) v + ( ν 0 E + F ) u 0 + τ q Γ ( q + 1 ) j = 1 k 1 w τ j ( k j + 1 ) q ( k j ) q + τ q Γ ( q + 1 ) v .
We will prove that Equation (4) is surjective for small values of τ by proving that Φ is weakly continuous and coercive (see [27], Theorem 1.2).
Since D , E , F are linear and bounded, they are weakly continuous. From ([27], Proposition 2.6), we deduce that Φ is weakly continuous. It remains to be proved that Φ is coercive. Denote
z τ = u 0 + τ q Γ ( q + 1 ) j = 1 k 1 w τ j ( k j + 1 ) q ( k j ) q ,
K τ = u 0 V + τ q Γ ( q + 1 ) j = 1 k 1 w τ j V ( k j + 1 ) q ( k j ) q ,
C τ = τ q Γ ( q + 1 ) , v τ = z τ + C τ v , v V .
Then we have
( ν 0 E + F ) v τ , v V = ( ν 0 E + F ) v τ , 1 C τ ( v τ z τ ) V = 1 C τ ( ν 0 E + F ) v τ , v τ V 1 C τ ( ν 0 E + F ) v τ , z τ V c 5 ν 0 C τ v τ V 2 1 C τ ( ν 0 E + F ) v τ , z τ V c 5 ν 0 C τ v τ V 2 1 C τ z τ V ( ν 0 E + F ) v τ V c 5 ν 0 C τ v τ V 2 1 C τ z τ V ( c 2 ν 0 v τ V + c 3 v τ V 2 ) = c 5 ν 0 C τ v τ V 2 1 C τ z τ V ( c 2 ν 0 z τ + C τ v V + c 3 z τ + C τ v V 2 ) c 5 ν 0 C τ v τ V 2 1 C τ z τ V ( c 2 ν 0 z τ V + c 2 ν 0 C τ v V + c 3 z τ V 2 + 2 c 3 C τ z τ V v V + c 3 C τ 2 v V 2 ) c 5 ν 0 C τ v τ V 2 c 3 C τ z τ V v V 2 ( c 2 ν 0 + 2 c 3 z τ V ) v V 1 C τ z τ V ( c 2 ν 0 z τ V + c 3 z τ V 2 ) .
This gives
Φ v , v V ( c 4 + c 5 α 2 ) v V 2 + c 5 ν 0 C τ v τ V 2 c 3 C τ z τ V v V 2 ( c 2 ν 0 + 2 c 3 z τ V ) v V 1 C τ z τ V ( c 2 ν 0 z τ V + c 3 z τ V 2 ) ( c 4 + c 5 α 2 ) v V 2 c 3 τ q Γ ( q + 1 ) K τ v V 2 ( c 2 ν 0 + 2 c 3 K τ ) v V Γ ( q + 1 ) τ q K τ ( c 2 ν 0 K τ + c 3 K τ 2 ) ( c 4 + c 5 α 2 c 3 τ q K τ Γ ( q + 1 ) ) v V 2 ( c 2 ν 0 + 2 c 3 K τ ) v V Γ ( q + 1 ) τ q K τ ( c 2 ν 0 K τ + c 3 K τ 2 ) .
It is easily seen that we can find sufficiently small τ 0 > 0 such that
c 4 + c 5 α 2 c 3 τ q K τ Γ ( q + 1 ) > 0 , τ ( 0 , τ 0 ) .
Thus, Φ is coercive for all τ ( 0 , τ 0 ) . The proof is complete.  □
Lemma 4.
There exists τ 0 > 0 such that for all τ ( 0 , τ 0 ) ,
max k = 1 , , N w τ k V C ,
max k = 1 , , N v τ k V C
for a constant C > 0 independent of τ.
Proof. 
We first prove estimate (6). For all 1 k N , we have
( D + α 2 E ) w τ k + ( ν 0 E + F ) v τ k = f τ k .
Multiplying the above equation by w τ k , we can obtain
f τ k , w τ k V = ( D + α 2 E ) w τ k , w τ k V + ( ν 0 E + F ) v τ k , w τ k V ( c 4 + c 5 α 2 c 3 τ q K τ Γ ( q + 1 ) ) w τ k V 2 + α 2 Γ ( q + 1 ) τ q v τ k V 2 ( c 2 ν 0 + 2 c 3 K τ ) w τ k V Γ ( q + 1 ) τ q K τ ( c 2 ν 0 K τ + c 3 K τ 2 ) .
Hence, we obtain
( c 4 + c 5 α 2 c 3 τ q K τ Γ ( q + 1 ) ) w τ k V 2 ( α f + c 2 ν 0 + 2 c 3 K τ ) w τ k V + Γ ( q + 1 ) τ q K τ ( c 2 ν 0 K τ + c 3 K τ 2 ) α 2 v τ k V 2 .
We can choose sufficiently small τ 0 > 0 such that
c 4 + c 5 α 2 c 3 τ q K τ Γ ( q + 1 ) > c 4 + c 5 α 2 2 , τ ( 0 , τ 0 ) .
With the above τ , the following two cases are discussed:
(i)
c 2 ν 0 u 0 V 2 + c 3 u 0 V 3 α 2 v τ k V 2 ;
(ii)
c 2 ν 0 u 0 V 2 + c 3 u 0 V 3 < α 2 v τ k V 2 .
Let (i) hold. We have
( c 4 + c 5 α 2 ) w τ k V 2 ( D + α 2 E ) w τ k , w τ k = f τ k , w τ k ( ν 0 E + F ) v τ k , w τ k f τ k V w τ k V + ( ν 0 E + F ) v τ k V w τ k V α f w τ k V + ( c 2 ν 0 v τ k V + c 3 v τ k V 2 ) w τ k V ,
and hence
( c 4 + c 5 α 2 ) w τ k V α f + c 2 ν 0 v τ k V + c 3 v τ k V 2 ,
which implies that { w τ k } is bounded since { v τ k } is bounded.
Let (ii) hold. We have
w τ k V 2 2 c 4 + c 5 α 2 ( α f + c 2 ν 0 + 2 c 3 K τ ) w τ k V + 2 Γ ( q + 1 ) ( c 4 + c 5 α 2 ) τ q K τ ( c 2 ν 0 K τ + c 3 K τ 2 ) α 2 v τ k V 2 = 2 c 4 + c 5 α 2 ( α f + c 2 ν 0 + 2 c 3 K τ ) w τ k V + 2 Γ ( q + 1 ) ( c 4 + c 5 α 2 ) τ q ( u 0 V + τ q Γ ( q + 1 ) j = 1 k 1 w τ j V ( k j + 1 ) q ( k j ) q × ( c 2 ν 0 u 0 V + τ q Γ ( q + 1 ) j = 1 k 1 w τ j V ( k j + 1 ) q ( k j ) q + c 3 u 0 V + τ q Γ ( q + 1 ) j = 1 k 1 w τ j V ( k j + 1 ) q ( k j ) q 2 ) α 2 v τ k V 2 ) = 2 c 4 + c 5 α 2 ( α f + c 2 ν 0 + 2 c 3 K τ ) w τ k V + 2 Γ ( q + 1 ) ( c 4 + c 5 α 2 ) τ q ( c 2 ν 0 u 0 V 2 + c 3 u 0 V 3 α 2 v τ k V 2 + τ q Γ ( q + 1 ) ( m 1 + τ q Γ ( q + 1 ) ( c 2 ν 0 + 3 c 3 u 0 V ) j = 1 k 1 w τ j ( k j + 1 ) q ( k j ) q V + c 3 τ 2 q Γ 2 ( q + 1 ) j = 1 k 1 w τ j V ( k j + 1 ) q ( k j ) q 2 ) × j = 1 k 1 w τ j V ( k j + 1 ) q ( k j ) q ) 2 c 4 + c 5 α 2 ( α f + c 2 ν 0 + 2 c 3 K τ ) w τ k V + 2 c 4 + c 5 α 2 ( m 1 + τ q Γ ( q + 1 ) ( c 2 ν 0 + 3 c 3 u 0 V ) × j = 1 k 1 w τ j V ( k j + 1 ) q ( k j ) q + c 3 τ 2 q Γ 2 ( q + 1 ) j = 1 k 1 w τ j V ( k j + 1 ) q ( k j ) q 2 ) × j = 1 k 1 w τ j V ( k j + 1 ) q ( k j ) q ,
where m 1 = b ¯ 0 + 2 b ¯ 1 u 0 V + 3 b ¯ 2 u 0 V 2 . It follows that
w τ k V 2 c 4 + c 5 α 2 ( α f + c 2 ν 0 + 2 c 3 K τ ) + ( 2 c 4 + c 5 α 2 ( m 1 + τ q Γ ( q + 1 ) ( c 2 ν 0 + 3 c 3 u 0 V ) × j = 1 k 1 w τ j V ( k j + 1 ) q ( k j ) q + c 3 τ 2 q Γ 2 ( q + 1 ) × j = 1 k 1 w τ j V ( k j + 1 ) q ( k j ) q 2 ) × j = 1 k 1 w τ j V ( k j + 1 ) q ( k j ) q ) 1 2 2 c 4 + c 5 α 2 ( α f + c 2 ν 0 + 2 c 3 K τ ) + 1 c 4 + c 5 α 2 j = 1 k 1 w τ j V ( k j + 1 ) q ( k j ) q + 1 2 ( m 1 + τ q Γ ( q + 1 ) ( c 2 ν 0 + 3 c 3 u 0 V ) j = 1 k 1 w τ j V ( k j + 1 ) q ( k j ) q + c 3 τ 2 q Γ 2 ( q + 1 ) j = 1 k 1 w τ j V ( k j + 1 ) q ( k j ) q 2 ) m 2 + m 3 j = 1 k 1 w τ j V ( k j + 1 ) q ( k j ) q + m 4 j = 1 k 1 w τ j V ( k j + 1 ) q ( k j ) q 2 ,
where
m 2 = 2 c 4 + c 5 α 2 ( α f + c 2 ν 0 + 2 c 3 u 0 V ) + 1 2 m 1 ,
m 3 = 4 c 3 T q ( c 4 + c 5 α 2 ) Γ ( q + 1 ) + T q 2 Γ ( q + 1 ) ( c 2 ν 0 + 3 c 3 u 0 V ) , m 4 = c 3 T 2 q Γ 2 ( q + 1 ) .
Moreover, for w τ 1 , we have
α f w τ 1 V f τ 1 , w τ 1 V = ( D + α 2 E ) w τ 1 , w τ 1 V + ( ν 0 E + F ) v τ 1 , w τ 1 V c 4 + c 5 α 2 w τ 1 V 2 + Γ ( q + 1 ) τ q A v τ 1 , v τ 1 V Γ ( q + 1 ) τ q ( ν 0 E + F ) v τ 1 , u 0 V ( c 4 + c 5 α 2 ) w τ 1 V 2 + α 2 Γ ( q + 1 ) τ q v τ 1 V 2 Γ ( q + 1 ) τ q u 0 V c 2 ν 0 u 0 + τ q Γ ( q + 1 ) w τ 1 V + c 3 u 0 + τ q Γ ( q + 1 ) w τ 1 V 2 ( c 4 + c 5 α 2 c 3 τ q u 0 V Γ ( q + 1 ) ) w τ 1 V 2 ( c 2 ν 0 u 0 V + 2 c 3 u 0 V 2 ) w τ 1 V Γ ( q + 1 ) τ q u 0 V ( c 2 ν 0 u 0 V + c 3 u 0 V 2 ) + α 2 Γ ( q + 1 ) τ q v τ 1 V 2 ,
and hence
( c 4 + c 5 α 2 c 3 τ q u 0 V Γ ( q + 1 ) ) w τ 1 V 2 ( α f + c 2 ν 0 u 0 V + 2 c 3 u 0 V 2 ) w τ 1 V + Γ ( q + 1 ) τ q ( c 2 ν 0 u 0 V 2 + c 3 u 0 V 3 α 2 v τ 1 V 2 ) ( α f + c 2 ν 0 u 0 V + 2 c 3 u 0 V 2 ) w τ 1 V .
Obviously, c 4 + c 5 α 2 c 3 τ q u 0 V Γ ( q + 1 ) > c 4 + c 5 α 2 2 for all τ ( 0 , τ 0 ) . Then
w τ 1 V 2 c 4 + c 5 α 2 ( α f + c 2 ν 0 u 0 V + 2 c 3 u 0 V 2 ) ,
which gives the boundedness of w τ 1 . Therefore, we can obtain (6) by induction.
Furthermore, estimate (7) follows from (6). The proof of Lemma 4 is complete.  □
In the final step of the Rothe method, define w ¯ τ , v ¯ τ V , f ¯ τ V by, respectively,
w ¯ τ ( t ) = w τ k , v ¯ τ ( t ) = v τ k , f ¯ τ ( t ) = f τ k ,
for a.e. t ( t k 1 , t k ] with k = 1 , , N .
Lemma 5.
There exists w V solving Problem 1, satisfying, for a subsequence,
w ¯ τ w w e a k l y   i n V .
Proof. 
Consider the Nemytski operators D , E , : V V , F : W V corresponding to D , E , F , respectively, i.e.,
( D v ) ( t ) = D v ( t ) , ( E v ) ( t ) = E v ( t ) , ( F v ¯ ) ( t ) = F v ¯ ( t ) for v s . V , v ¯ W a . e . t ( 0 , T ) .
Then Equation (8) is transferred to be
( ν 0 D + E ) w ¯ τ + ( ν 0 E + F ) v ¯ τ = f ¯ τ .
It follows from (6) that { w ¯ τ } is bounded in V . Then we obtain, for a subsequence,
w ¯ τ w weakly   in V .
We are left with the task of determining that w solves Problem 1. Since the operator I t q 0 : V C ( [ 0 , T ] ; V ) is linear and continuous, it is also weakly continuous, and hence
I t q 0 w ¯ τ I t q 0 w weakly   in C ( [ 0 , T ] ; V ) , as τ 0 .
Then
I t q 0 w ¯ τ ( t ) I t q 0 w ( t ) weakly   in V , as τ 0 , t [ 0 , T ] .
Moreover, there exists a constant C > 0 such that w ¯ τ ( s ) C for a.e. s ( 0 , T ) . For t ( t k 1 , t k ] , we obtain that
v ¯ τ ( t ) u 0 I t q 0 w ¯ τ ( t ) V = 1 Γ ( q ) 0 t k ( t k s ) q 1 w ¯ τ ( s ) d s 0 t ( t s ) q 1 w ¯ τ ( s ) d s V 1 Γ ( q ) 0 t ( t k s ) q 1 ( t s ) q 1 w ¯ τ ( s ) d s V + 1 Γ ( q ) t t k ( t k s ) q 1 w ¯ τ ( s ) d s V C Γ ( q ) 0 t | ( t k s ) q 1 ( t s ) q 1 | d s + C Γ ( q ) t t k ( t k s ) q 1 d s
C Γ ( q + 1 ) | t k q ( t k t ) q t q | + C Γ ( q ) ( t k t ) q .
We can obtain that v ¯ τ ( t ) u 0 I t q 0 w ¯ τ ( t ) 0 , as τ 0 for all t [ 0 , T ] . Then
v ¯ τ ( t ) u 0 + I t q 0 w ( t ) weakly   in V , as τ 0
for all t [ 0 , T ] .
Furthermore, it is clear that D , E are weakly continuous. Hence
( ν 0 D + E ) w τ ( ν 0 D + E ) w weakly   in V , as τ 0 .
It follows from [2,27] that F is weakly continuous. Then for any v V , we infer that
lim τ 0 ( ν 0 E + F ) v ¯ τ , v V = lim τ 0 0 T ( ν 0 E + F ) v ¯ τ ( t ) , v ( t ) V d t = 0 T lim τ 0 ( ν 0 E + F ) v ¯ τ ( t ) , v ( t ) V d t = 0 T ( ν 0 E + F ) ( u 0 + I t q 0 w ( t ) ) , v ( t ) V d t = ( ν 0 E + F ) ( u 0 + I t q 0 w ) , v V .
Therefore,
( ν 0 E + F ) v ¯ τ ( ν 0 E + F ) ( u 0 + I t q 0 w ) weakly   in V , as τ 0 .
Together with f ¯ τ f in V , as τ 0 (see [28], Lemma 3.3), from (10), (16), and (15), taking the limit in (9), we finally obtain
( ν 0 D + E ) w + ( ν 0 E + F ) ( u 0 + I t q 0 w ) = f ,
which proves that w V is a solution to Problem 1. The proof of Lemma 5 is complete.  □
In the end, combining the above lemmas, we obtain Theorem 1.
By apply the method above, we can also solve the following generalized time-fractional Navier–Stokes–Voigt equation:
ρ 1 D t q 0 C u ρ 2 Δ u ρ 3 Δ ( D t q 0 C u ) + ( u · ) u + p = f , in Q , · u = 0 , in Q , u ( x , t ) = 0 , on Ω × ( 0 , T ) , u ( x , 0 ) = u 0 ( x ) , in Ω .
Here, ρ i L ( Ω ; R + ) ( i = 1 , 2 , 3 ) .

Author Contributions

Writing—original draft preparation, K.X.; writing—review and editing, B.Z. All authors have read and agreed to the published version of the manuscript.

Funding

The first author is supported by the Basic Ability Improvement Project for Middle-Aged and Young Teachers of Universities in Guangxi (No. 2022KY1210). The second author is supported by the Special Fund for Science and Technological Bases and Talents of Guangxi (No. GUIKE AD21220103), the National Natural Science Foundation of China (No. 12201137), and the Start-up Project of Scientific Research on Introducing talents at school level in Guangxi Minzu University (No. 2019KJQD04).

Data Availability Statement

No new data were created or analyzed in this study. Data sharing is not applicable to this article.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Anh, C.T.; Nguyet, T.M. Optimal control of the instationary three dimensional Navier-Stokes-Voigt equations. Numer. Funct. Anal. Optim. 2016, 37, 415–439. [Google Scholar] [CrossRef]
  2. Anh, C.T.; Trang, P.T. Pull-back attractors for three-dimensional Navier-Stokes-Voigt equations in some unbounded domains. Proc. R. Soc. Edinb. Sect. A 2013, 143, 223–251. [Google Scholar] [CrossRef]
  3. García-Luengo, J.; Marín-Rubio, P.; Real, J. Pullback attractors for three-dimensional non-autonomous Navier-Stokes-Voigt equations. Nonlinearity 2012, 25, 905–930. [Google Scholar] [CrossRef]
  4. Kalantarov, V.K.; Titi, E.S. Gevrey regularity for the attractor of the 3D Navier-Stokes-Voight equations. J. Nonlinear Sci. 2009, 19, 133–152. [Google Scholar] [CrossRef]
  5. Yue, G.; Zhong, C.K. Attractors for autonomous and nonautonomous 3D Navier-Stokes-Voigt equations. Discret. Contin. Dyn. Syst. Ser. B 2011, 16, 985–1002. [Google Scholar]
  6. Ahn, J.; Kim, J.; Lee, J. Coriolis effect on temporal decay rates of global solutions to the fractional Navier-Stokes equations. Math. Ann. 2022, 383, 259–289. [Google Scholar] [CrossRef]
  7. Carvalho-Neto, P.M.D.; Planas, G. Mild solutions to the time fractional Navier-Stokes equations in R n . J. Differ. Equ. 2015, 259, 2948–2980. [Google Scholar] [CrossRef]
  8. Li, L.; Liu, J.-G. Some compactness criteria for weak solutions of time fractional PDE’s. SIAM J. Math. Anal. 2018, 50, 3963–3995. [Google Scholar] [CrossRef]
  9. Li, X.C.; Yang, X.Y.; Zhang, Y.H. Error estimates of mixed finite element methods for time-fractional Navier-Stokes equations. J. Sci. Comput. 2017, 70, 500–515. [Google Scholar] [CrossRef]
  10. Peng, L.; Debbouche, A.; Zhou, Y. Existence and approximations of solutions for time-fractional Navier-Stokes equations. Math. Meth. Appl. Sci. 2018, 41, 8973–8984. [Google Scholar] [CrossRef]
  11. Zhou, Y.; Peng, L. On the time-fractional Navier-Stokes equations. Comput. Math. Appl. 2017, 73, 874–891. [Google Scholar] [CrossRef]
  12. Zhou, Y.; Peng, L. Weak solutions of the time-fractional Navier-Stokes equations and optimal control. Comput. Math. Appl. 2017, 73, 1016–1027. [Google Scholar] [CrossRef]
  13. Zeng, B. Feedback control for non-stationary 3D Navier-Stokes-Voigt equations. Math. Mech. Solids 2020, 25, 2210–2221. [Google Scholar] [CrossRef]
  14. Zeng, B. Feedback control for nonlinear evolutionary equations with applications. Nonlinear Anal. RWA 2022, 66, 103535. [Google Scholar] [CrossRef]
  15. Zeng, B. Existence for a class of time-fractional evolutionary equations with applications involving weakly continuous operator. Fract. Calc. Appl. Anal. 2023, 26, 172–192. [Google Scholar] [CrossRef]
  16. Zeng, B.; Migórski, S. Evolutionary subgradient inclusions with nonlinear weakly continuous operators and applications. Comput. Math. Appl. 2018, 75, 89–104. [Google Scholar] [CrossRef]
  17. Hilfer, R. Application of Fractional Calculus in Physics; World Scientific: Singapore, 2000. [Google Scholar]
  18. Kilbas, A.A.; Srivastava, H.M.; Trujillo, J.J. Theory and Applications of Fractional Differential Equations; North-Holland Mathematics Studies 204; Elsevier Science B.V.: Amsterdam, The Netherland, 2006. [Google Scholar]
  19. Lakshmikantham, V.; Vatsala, A.S. Basic theory of fractional differential equations. Nonlinear Anal. 2008, 69, 2677–2682. [Google Scholar] [CrossRef]
  20. Podlubny, I. Fractional Differential Equations; Academic Press: San Diego, CA, USA, 1999. [Google Scholar]
  21. Contantin, P.; Foias, C. Navier-Stokes Equations; Chicago Lectures in Mathematics; University of Chicago Press: Chicago, IL, USA, 1988. [Google Scholar]
  22. Temam, R. Navier-Stokes Equations: Theory and Numerical Analysis; American Mathematical Society: Providence, RI, USA, 2001. [Google Scholar]
  23. Kačur, J. Method of Rothe in Evolution Equations; Teubner-Texte zur Mathematik 80; B.G. Teubner: Leipzig, Germany, 1985. [Google Scholar]
  24. Roubiček, T. Nonlinear Partial Differential Equations with Applications; Birkhäuser: Basel, Switzerland; Boston, MA, USA; Berlin, Germany, 2005. [Google Scholar]
  25. Han, J.F.; Migórski, S.; Zeng, H.D. Weak solvability of a fractional viscoelastic frictionless contact problem. Appl. Math. Comput. 2017, 303, 1–18. [Google Scholar] [CrossRef]
  26. Zeng, S.D.; Migórski, S. A class of time-fractional hemivariational inequalities with application to frictional contact problem. Commun. Nonlinear Sci. Numer. Simul. 2018, 56, 34–48. [Google Scholar] [CrossRef]
  27. Francǔ, J. Weakly continuous operators, Applications to differential equations. Appl. Math. 1994, 39, 45–56. [Google Scholar] [CrossRef]
  28. Carstensen, C.; Gwinner, J. A theory of discretization for nonlinear evolution inequalities applied to parabolic Signorini problems. Ann. Mat. Pura Appl. 1999, 177, 363–394. [Google Scholar] [CrossRef]
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Xu, K.; Zeng, B. Existence Result for a Class of Time-Fractional Nonstationary Incompressible Navier–Stokes–Voigt Equations. Axioms 2024, 13, 499. https://doi.org/10.3390/axioms13080499

AMA Style

Xu K, Zeng B. Existence Result for a Class of Time-Fractional Nonstationary Incompressible Navier–Stokes–Voigt Equations. Axioms. 2024; 13(8):499. https://doi.org/10.3390/axioms13080499

Chicago/Turabian Style

Xu, Keji, and Biao Zeng. 2024. "Existence Result for a Class of Time-Fractional Nonstationary Incompressible Navier–Stokes–Voigt Equations" Axioms 13, no. 8: 499. https://doi.org/10.3390/axioms13080499

APA Style

Xu, K., & Zeng, B. (2024). Existence Result for a Class of Time-Fractional Nonstationary Incompressible Navier–Stokes–Voigt Equations. Axioms, 13(8), 499. https://doi.org/10.3390/axioms13080499

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop