Citation: Xueli Song, Jianhua Wu. Non-autonomous 3D Brinkman-Forchheimer equation with singularly oscillating external force and its uniform attractor[J]. AIMS Mathematics, 2020, 5(2): 1484-1504. doi: 10.3934/math.2020102
[1] | Dagmar Medková . Classical solutions of the Dirichlet problem for the Darcy-Forchheimer-Brinkman system. AIMS Mathematics, 2019, 4(6): 1540-1553. doi: 10.3934/math.2019.6.1540 |
[2] | Ailing Ban . Asymptotic behavior of non-autonomous stochastic Boussinesq lattice system. AIMS Mathematics, 2025, 10(1): 839-857. doi: 10.3934/math.2025040 |
[3] | Xue-li Song, Yuan-yuan Liu, Xiao-tian Xie . The existence of uniform attractors for the 3D micropolar equations with nonlinear damping term. AIMS Mathematics, 2024, 9(4): 9608-9630. doi: 10.3934/math.2024470 |
[4] | Xin Liu . Stability of random attractors for non-autonomous stochastic $ p $-Laplacian lattice equations with random viscosity. AIMS Mathematics, 2025, 10(3): 7396-7413. doi: 10.3934/math.2025339 |
[5] | Rou Lin, Min Zhao, Jinlu Zhang . Random uniform exponential attractors for non-autonomous stochastic Schrödinger lattice systems in weighted space. AIMS Mathematics, 2023, 8(2): 2871-2890. doi: 10.3934/math.2023150 |
[6] | Bengisen Pekmen Geridönmez . Numerical investigation of ferrofluid convection with Kelvin forces and non-Darcy effects. AIMS Mathematics, 2018, 3(1): 195-210. doi: 10.3934/Math.2018.1.195 |
[7] | Xiaobin Yao . Random attractors for non-autonomous stochastic plate equations with multiplicative noise and nonlinear damping. AIMS Mathematics, 2020, 5(3): 2577-2607. doi: 10.3934/math.2020169 |
[8] | Rinko Miyazaki, Dohan Kim, Jong Son Shin . Uniform boundedness of solutions to linear difference equations with periodic forcing functions. AIMS Mathematics, 2023, 8(10): 24116-24131. doi: 10.3934/math.20231229 |
[9] | Xiao Bin Yao, Chan Yue . Asymptotic behavior of plate equations with memory driven by colored noise on unbounded domains. AIMS Mathematics, 2022, 7(10): 18497-18531. doi: 10.3934/math.20221017 |
[10] | Li Yang . Pullback random attractors of stochastic strongly damped wave equations with variable delays on unbounded domains. AIMS Mathematics, 2021, 6(12): 13634-13664. doi: 10.3934/math.2021793 |
For ρ∈[0,1) and ε>0, we consider the following non-autonomous Brinkman-Forchheimer equation with a singularly oscillating external force on Ω⊂R3:
{ut−νΔu+au+b|u|u+c|u|βu+∇p=f0(x,t)+ε−ρf1(x,tε),in Ω×(τ,T),divu=0,in Ω×(τ,T),u|t=τ=uτ,in Ω,u=0,on ∂Ω×(τ,T), | (1.1) |
where Ω is an open, bounded domain of R3 with sufficiently smooth boundary ∂Ω. u=(u1,u2,u3) is the fluid velocity vector, ν is the Brinkman coefficient, a>0 is the Darcy coefficient, b>0,c>0 are the Forchheimer coefficients, p is the pressure, β∈(1,43] is a constant.
Along with (1.1), we consider the averaged Brinkman-Forchheimer equation
{ut−νΔu+au+b|u|u+c|u|βu+∇p=f0(x,t),in Ω×(τ,T),divu=0,in Ω×(τ,T),u|t=τ=uτ,in Ω,u=0,on ∂Ω×(τ,T), | (1.2) |
without by rapid and singular oscillations, which formally corresponds to ε=0.
The Brinkman-Forchheimer equation describes the motion of fluid flow in a saturated porous medium and has been studied by many researchers. We should note that most of the previous studies have been focused on physical viewpoint or numerical simulation viewpoint(see [1,2,3,4,5,6,7,8]), or continuous dependence of solutions on the coefficients ν,b and c (see [9,10,11,12,13,14,15]). The asymptotic behavior of solutions was examined in [16,17,18,19,20,21,22], where [16,17,18,19,20,21] were mainly for the case of the parameter β=2. In [16], using condition (C) method, U˘gurlu showed the existence of global attractor in H10(Ω). In [17], Wang and Lin showed that Brinkman-Forchheimer equation has a global attractor in H2(Ω) by a very clever way. In [18], using the method of regularization of solutions and the compact embedding to deduce the uniformly ω-limit compactness of the associated evolutionary process, You, Zhao and Zhou proved the existence of uniform attractor in H10(Ω). In [19], the existence of D-pullback attractors was deduced by establishing the D-pullback asymptotical compactness of θ-cocycle. In [20], Song and Qiao proved the existence and structure of the uniform attractor in H10(Ω) for the processes associated to the fluid when the external force f0(x,t) is translation compact in L2loc(R,(L2(Ω))3) and investigated the averaging problems of the equations with oscillating external forces. In [21], Zhang, Su and Wen investigated the existence of global attractor and uniform attractor for the 3D autonomous and nonautonomous Brinkman-Forchheimer equations. For β≠2, in [22], using condition (C) method, Ouyang and Yang proved the existence of global attractor in H10(Ω) when 1<β≤43.
Stability of attractors for a dynamical system with some oscillating (or perturbed) external forces is very important in natural phenomenon. Indeed, this issue has been considered by some mathematicians and engineers. Chepyzhov et al. [23] studied the non-autonomous sine-Gordon type equations with rapidly oscillating external force. Efendiev and Zelik [24,25] considered the reaction-diffusion systems with rapidly oscillating coefficients and nonlinear rapidly oscillating in time. Chepyzhov and Vishik [26,27,28] investigated the Navier-Stokes equations with terms that rapidly oscillate with respect to spatial and time variables. Qin et al. [29] investigated the uniform attractors for a 3D non-autonomous Navier-Stokes-Voight equation with singularly oscillating forces. Anh and Toan [30] considered the nonclassical diffusion equation on RN(N≥3) with a singularly oscillating external force. Medjo [31] and [32] used different methods to discuss non-autonomous planetary 3D geostrophic equation with oscillating external force and its global attractor. Medjo [33] investigated a non-autonomous two-phase flow model with oscillating external force and its global attractor. As far as we know, there is almost no paper dealing with 3D Brinkman-Forchheimer equations with rapidly oscillating terms have been published.
Motivated by [22] and [23,24,25,26,27,28,29,30,31,32,33], we consider the properties of (1.1), depending on the small parameter ε, which reflects the rate of fast time oscillation in the term ε−ρf1(tε,x), having the growing amplitude of order ε−ρ. By using the method in [27,29,30,32,33], under suitable assumptions on the external force, we prove the stability of the uniform attractor Aε(0<ε≤1) associated to problem (1.1)-(1.2) as ε→0+ in space H. The uncertainty of parameter β brings a lot of trouble to our proof, because when β is too large, some Sobolev embedding inequlaities can not be used. In the proof process of this paper, we finally determine that 1<β≤43, just like in [22].
The main purpose of this paper is to show:
(1) the uniform (w.r.t.ε) boundedness of the family Aε in H which is defined in Section 3:
supε∈[0,1]∥Aε∥H<+∞; |
(2) the convergence of Aε to A0 as ε→0+ in the standard Hausdorff semidistance in H, i.e., :
limε→0+distH(Aε,A0)=0. |
This paper is organized as follows. In Section 2, we present the notations and preliminaries that are required for this study. In Section 3, we show the existence of uniform attractor Aε and we demonstrate the structure of the uniform attractor. In Section 4, we verify the uniform boundedness of the uniform attractor Aε. In Section 5, we prove the convergence Aε→A0 as ε→0+.
Nomenclature | |||
u | fluid velocity vector(m/s) | p | pressure |
R3 | three-dimensional whole space | Ω | open, bounded domain of R3 |
u1,u2,u3 | velocity components | t | time |
∂Ω | boundary of a domain | ν | Brinkman coefficient |
a | Darcy coefficient | b,c | Forchheimer coefficient |
β | power of nonlinear term | f | external force |
Given a space X, we usually denote the norm in X by ∥⋅∥X, and we indicate by
distX(B1,B2)=supb1∈B1infb2∈B2∥b1−b2∥X, |
the Hausdorff semidistance in X from a set B1 to a set B2. Throughout this paper, we set Rτ=[τ,+∞),τ∈R. C will stand for a generic positive constant, which is different from line to line or even in the same line.
The mathematical setting of our problem is similar to that of the Navier-Stokes equations. Let us introduce the following spaces
V={u∈(C∞0(Ω))3:divu=0},H=cl(L2(Ω))3V,V=cl(H10(Ω))3V, |
where clX denotes the closure in the space X. Operator P is the Helmholtz-Leray orthogonal projection from (L2(Ω))3 onto H. A:=−PΔ is the Stokes operator subject to the nonslip homogeneous Dirichlet boundary with the domain (H2(Ω))3∩V, and A is a self-adjoint positively defined operator on H. We define, for σ∈R, the scale of Hilbert spaces
Hσ:=D(Aσ2) |
with inner products and norms
⟨u,v⟩σ:=⟨Aσ2u,Aσ2v⟩(L2(Ω))3,∥u∥Hσ:=∥Aσ2u∥(L2(Ω))3. |
In particular,
H0=H,H1=V,H2=D(A), |
and we have the generalized Poincarˊe inequality
∥u∥Hσ+1≥λ121∥u∥Hσ,∀u∈Hσ+1, | (2.1) |
where λ1 is the first eigenvalue of the Stokes operator A.
In this paper, we use (⋅,⋅) and ∥⋅∥ denote the product and the norm of H, i.e.,
(u,v)=∫Ωu⋅vdx,∀u,v∈H, ∥u∥=(u,u)12, |
((⋅,⋅)) and ∥⋅∥V denote the product and norm of V, i.e.,
((u,v))=3∑i=1∫Ω∇ui⋅∇vidx,∀u,v∈V,∥u∥V=((u,u))12. |
In this paper, Lp(Ω)=(Lp(Ω))3, and we use ∥⋅∥p to denote the norm in Lp(Ω) (for p≠2). Assumptions on the external forces The function f0(x,t) and f1(x,t) are taken from the space L2b(R;H) of translation bounded functions in L2loc(R;H), with
∥f0∥2L2b:=supt∈R∫t+1t∥f0(s)∥2ds=M20, | (2.2) |
∥f1∥2L2b:=supt∈R∫t+1t∥f1(s)∥2ds=M21, | (2.3) |
for some constants M0,M1≥0.
Putting
fε(x,t):={f0(x,t)+ε−ρf1(x,tε),ε>0,f0(x,t),ε=0. | (2.4) |
It is easy to check that fε∈L2b(R;H), and
∥fε∥L2b≤Qε:={M0+√2M1ε−ρ,ε>0,M0,ε=0. |
Now we recall some inequalities and a Gronwall-type lemma that will be needed in the sequel.
Lemma 2.1. ([34]) Let p∈[2,∞). Then for every a,b∈R,
(|a|p−2a−|b|p−2b)(a−b)≥2−2−p3−p/2|a−b|p. |
Lemma 2.2. ([27]) For every τ∈R, every nonnegative locally summable function φ on Rτ and every β>0, we have
∫tτφ(s)e−β(t−s)ds≤11−e−βsupθ≥τ∫θ+1θφ(s)ds, | (2.5) |
for all t≥τ.
Lemma 2.3. ([23]) Let a real function z(t),t≥0 be uniformly continuous and satisfy the inequality
dzdt+γz(t)≤f(t),∀t≥0, |
where γ>0,f(t)≥0 for all t≥0 and f∈L1loc(R+). Suppose also that
∫t+1tf(s)ds≤M,∀t≥0. |
Then
z(t)≤z(0)e−γt+M(1+γ−1),∀t≥0. |
We rewrite (1.1) and (1.2) in the abstract form
{ut+νAu+au+G(u)=fε(x,t),u|t=τ=uτ, | (3.1) |
where the pressure p has disappeared by force of the application of the Helmholtz-Leray projection P, and G(u)=P(b|u|u+c|u|βu), fε(x,t)=f0(x,t)+ε−ρf1(x,tε) and ε>0 is fixed.
Proposition 3.1. Given fε∈L2b(R;H) and uτ∈H. Then the system (3.1) has a unique weak solution u(t) satisfying
u∈C(Rτ;H)∩L∞(Rτ;H)∩L2loc(Rτ;V). |
Proof. We can prove the global existence and uniqueness result by using the Faedo-Galerkin method (see [35,36]).
If the functions f0(t) and f1(t) are translation bounded, i.e. conditions (2.2) and (2.3) hold, equation (3.1) generates the dynamical process
{Ufε(t,τ),t≥τ,τ∈R} |
acting on H by the formula
Ufε(t,τ)uτ=u(t),t≥τ, |
where u(t) is the solution to (3.1).
Proposition 3.2. For any ε>0, the process {Ufε(t,τ)} associated to (3.1) is uniformly compact in H and it has a uniform (with respect to τ∈R) absorbing set Bε in H,
Bε={u∈H|∥u∥≤C0Qε}, | (3.2) |
where the constant C0 depends on ν and λ1. Moreover, there exists a uniform attractor Aε in H.
Proof. The proof process is similar to the proof in [21], so we omit it here.
We consider the hull H(fε) of fε(x,t) in the space L2loc(R;H):
H(fε)=[{fε(⋅,t+h)|h∈R}]L2loc(R;H). | (3.3) |
Recall that H(fε) is compact in L2loc(R;H) and each element ˆfε∈H(fε) can be written as
ˆfε(x,t)=ˆf0(x,t)+ε−ρˆf1(x,tε), | (3.4) |
with ˆf0∈H(f0) and ˆf1∈H(f1), where H(f0) and H(f1) are the hulls of f0 and f1 in L2loc(R;H) respectively.
We also note that (see [23])
∥ˆf0∥L2b(R;H)≤∥f0∥L2b(R;H),∀ˆf0∈H(f0), | (3.5) |
∥ˆf1∥L2b(R;H)≤∥f1∥L2b(R;H),∀ˆf1∈H(f1). | (3.6) |
It follows that
∥ˆfε∥L2b(R;H)≤∥f0∥L2b(R;H)+Cερ∥f1∥L2b(R;H),∀ˆfε∈H(fε), | (3.7) |
where C is independent of f0,f1,ρ and ε.
To describe the structure of the attractor Aˆfε, we consider the family of equations
dˆudt+νAˆu+aˆu+G(ˆu)=ˆfε(x,t), | (3.8) |
with the external force ˆfε∈H(fε).
For ˆfε∈H(fε), the Eq.(3.8) generates a process {Uˆfε(t,τ)} that satisfies the same properties as {Ufε(t,τ)}. Moreover, the process {Uˆfε(t,τ)} has a uniform attractor Aˆfε that satisfies Aˆfε⊂Afε.
Proposition 3.3. Let f0(x,t),f1(x,t) be translation compact in the space L2loc(R;H). Then for any fixed ε,0<ε≤1, the family of processes {Uˆfε(t,τ)},ˆfε∈H(fε) corresponding to (3.8) has an absorbing set Bε, which is bounded in H and satisfies
|Bε|H≤C+Cε−ρ. | (3.9) |
The family {Uˆfε(t,τ)},ˆfε∈H(fε) is (H×H(fε);H)-continuous. That is, if
ˆfεn→ˆfε in L2loc(R;H),uτn→uτ in H, | (3.10) |
then
Uˆfεn(t,τ)uτn→Uˆfε(t,τ)uτ in H. | (3.11) |
Proof. The first part of the proposition is proved in Proposition 3.2. Now, let us prove the second part. Let wn=ˆun−ˆu=Uˆfεn(t,τ)uτn−Uˆfε(t,τ)uτ. Then wn satisfies
dwndt+νAwn+awn+G(ˆun)−G(ˆu)=ˆfεn−ˆfε. | (3.12) |
Multiplying (3.12) by wn we have
12ddt∥wn∥2+ν∥∇wn∥2+a∥wn∥2+(G(ˆun)−G(ˆu),wn)=(ˆfεn−ˆfε,wn). | (3.13) |
Noticing b|u|u+c|u|βu is monotonic, i.e.,
(G(ˆun)−G(ˆu),wn)=(b|ˆun|ˆun+c|ˆun|βˆun−b|ˆu|ˆu−c|ˆu|βˆu,ˆun−ˆu)≥0, | (3.14) |
then (3.13) gives
ddt∥wn∥2+a∥wn∥2≤1a∥ˆfεn−ˆfε∥2. | (3.15) |
Applying Gronwall Lemma to (3.15) we have
∥wn(t)∥2≤∥wn(τ)∥2e−a(t−τ)+1a∫tτ∥ˆfεn−ˆfε∥2e−a(t−s)ds≤∥wn(τ)∥2+1a∫tτ∥ˆfεn−ˆfε∥2ds,∀t≥τ. | (3.16) |
Note that
ˆfεn→ˆfε in L2loc(R;H) and uτn→uτ in H as n→∞, | (3.17) |
therefore, it follows from (3.16) that
∥wn(t)∥=∥ˆun(t)−ˆu(t)∥→0 as n→∞, |
and (3.11) is proved, i.e., the family of processes {Uˆfε(t,τ)},ˆfε∈H(fε) is (H×H(fε);H)-continuous.
We denote by Kˆfε the kernel of (3.8) with the external force ˆfε∈H(fε). Let us recall that Kˆfε is the family of all complete solutions {ˆu(t),t∈R} of (3.8), which are uniformly bounded in H. The set
Kˆfε(s)={ˆu(s)|ˆu∈Kˆfε}⊂H |
is called the kernel section of Kˆfε at time t=s.
For every ε∈[0,1], the following representation of the uniform attractor Aε of equation (3.1) holds:
Aε=⋃ˆfε∈H(fε)Kˆfε(0). | (3.18) |
Actually, Kˆfε(0) can be replaced by Kˆfε(τ), for an arbitrary τ∈R.
First, we consider the auxiliary linear equation with nonautonomous external force and give some useful estimates and then prove the uniform boundedness of Aε in H.
Considering the linear equation
Vt+νAV+aV=K(t), V|t=τ=0, | (4.1) |
we get the following lemma.
Lemma 4.1. If K∈L2loc(R;V), then the above problem has a unique solution
V∈C(Rτ;H2)∩L2loc(Rτ;H3). |
Moreover, the inequalities
∥V(t)∥2≤C∫tτe−a(t−s)∥K(s)∥2ds, | (4.2) |
∥AV(t)∥2≤C∫tτe−2a(t−s)∥K(s)∥2Vds, | (4.3) |
∫t+1t∥∇V(s)∥2ds≤C(∥V(t)∥2+∫t+1t∥K(s)∥2ds) | (4.4) |
∫t+1t∥A32V(s)∥2ds≤C(∥AV(t)∥2+∫t+1t∥K(s)∥2Vds) | (4.5) |
hold for every t≥τ and some constant C>0, independent of the initial time τ∈R.
Proof. Multiplying the equation (4.1) by V and A2V respectively, we have
12ddt∥V∥2+ν∥∇V∥2+a∥V∥2=(K(t),V)≤a2∥V∥2+12a∥K(t)∥2, | (4.6) |
12ddt∥AV∥2+ν∥A32V∥2+a∥AV∥2=(K(t),A2V)≤12ν∥K(t)∥2V+ν2∥A32V∥2. | (4.7) |
It follows from (4.6) and (4.7) that
ddt∥V∥2+a∥V∥2≤1a∥K(t)∥2, |
ddt∥AV∥2+2a∥AV∥2≤1ν∥K(t)∥2V. |
Applying Gronwall lemma it yields
∥V(t)∥2≤1a∫tτe−a(t−s)∥K(s)∥2ds, |
∥AV(t)∥2≤1ν∫tτe−2a(t−s)∥K(s)∥2Vds. |
From (4.6) and (4.7) we also can get
ddt∥V∥2+2ν∥∇V∥2≤1a∥K(t)∥2, | (4.8) |
ddt∥AV∥2+ν∥A32V∥2≤1ν∥K(t)∥2V. | (4.9) |
Integrating (4.8) and (4.9) on [t,t+1] respectively, we obtain
2ν∫t+1t∥∇V(s)∥2ds≤∥V(t)∥2+1a∫t+1t∥K(s)∥2ds, | (4.10) |
ν∫t+1t∥A32V(s)∥2ds≤∥AV(t)∥2+1ν∫t+1t∥K(s)∥2Vds. |
The proof is finished.
Setting F(t,τ)=∫tτf1(s)ds,t≥τ, we assume that
supt≥τ,τ∈R(∥F(t,τ)∥2+∫t+1t∥F(s,τ)∥2Vds)≤l2. | (4.11) |
Lemma 4.2. Assume that f1∈L2loc(R;H) and satisfies (4.11). Then the solution v(t) to the Cauchy problem
vt+νAv+av=f1(tε), v|t=τ=0 | (4.12) |
with ε∈(0,1], satisfies the inequality
∥v(t)∥2+∫t+1t∥∇v(s)∥2ds≤Cl2ε2,∀t≥τ, | (4.13) |
where C>0 is a constant independent of f1.
Proof. Without loss of generality, we may assume τ=0. Denoting V(t)=∫t0v(s)ds, we have, for any t≥0,
∂tV(t)=v(t)=∫t0∂tv(s)ds, |
as v(0)=0. Integrating (4.12) in time, we see that the function V(t) solves the problem
∂tV+νAV+aV=Fε(t),V|t=0=0, | (4.14) |
with external force
Fε(t)=∫t0f1(sε)ds=ε∫tε0f1(s)ds=εF(tε,0). |
It follows from (4.11) that
supt≥0∥Fε(t)∥≤lε |
and
∫t+1t∥Fε(s)∥2Vds=ε3∫t+1εtε∥F(s,0)∥2Vds≤2ε2supt≥0{∫t+1t∥F(s,0)∥2Vds}≤2l2ε2. |
By (2.5) we have
∫t0e−a(t−s)∥Fε(s)∥2ds≤Cl2ε2,∫t0e−2a(t−s)∥Fε(s)∥2Vds≤Cl2ε2. |
So applying Lemma 5.1, we obtain
∥V(t)∥2+∥AV(t)∥2+∫t+1t∥∇V(s)∥2ds+∫t+1t∥A32V(s)∥2ds≤Cl2ε2. |
Hence, on account of (4.14) we have
∥v(t)∥=∥∂tV(t)∥≤∥Fε(t)∥+ν∥AV(t)∥+a∥V(t)∥≤Clε |
and
∥∇v(s)∥2=∥∇(∂tV(s))∥2≤3∥Fε(s)∥2V+3ν2∥A32V(s)∥2+3a2∥∇V(s)∥2, |
from which we derive the integral estimate
∫t+1t∥∇v(s)∥2ds≤Cl2ε2. |
This finishes the proof.
Theorem 4.1. Let (4.11) holds true. Then the uniform attractors Aε are uniformly (w.r.t. ε) bounded in H, that is,
supε∈[0,1]∥Aε∥<∞. |
Proof. Let u be the solution to (3.1) with initial data uτ∈H. For ε>0, we consider the problem
vt+νAv+av=ε−ρf1(tε), v|t=τ=0. | (4.15) |
Lemma 4.2 provides the estimate
∥v(t)∥2+∫t+1t∥∇v(s)∥2ds≤cl2ε2(1−ρ),∀t≥τ. | (4.16) |
Then, the function w(t)=u(t)−v(t) clearly satisfies the equation
wt−νΔw+aw+b|w+v|(w+v)+c|w+v|β(w+v)+∇p=f0 | (4.17) |
with initial condition w|t=τ=uτ. Taking the inner product of (4.17) with w in H, we obtain
12ddt∥w∥2+ν∥∇w∥2+a∥w∥2+b(|w+v|(w+v)−|v|v,w)+c(|w+v|β(w+v)−|v|βv,w)=−b(|v|v,w)−c(|v|βv,w)+(f0,w). | (4.18) |
By Lemma 2.1, we have
12ddt∥w∥2+ν∥∇w∥2+a∥w∥2+b⋅2−53−32∥w∥33+c⋅2−4−β3−β+22∥w∥β+2β+2≤−b(|v|v,w)−c(|v|βv,w)+(f0,w). | (4.19) |
Noticing
(f0,w)≤a2∥w∥2+12a∥f0∥2, | (4.20) |
b|(|v|v,w)|=b|∫Ω|v|vwdx|≤b(∫Ω|w|6dx)1/6(∫Ω(|v|v)6/5dx)5/6=b∥w∥6∥v∥2125≤ν2d20∥w∥26+C∥v∥4125, | (4.21) |
and
c|(|v|βv,w)|=c|∫Ω|v|βvwdx|≤c(∫Ω|w|6dx)16⋅(∫Ω(|v|βv)6/5dx)5/6=c∥w∥6∥v∥β+165(β+1)≤ν2d20∥w∥26+C∥v∥2(β+1)65(β+1), | (4.22) |
according to Sobolev inequality
∥v∥p≤d0∥∇v∥,1≤p≤6, | (4.23) |
combining (4.20)-(4.23) with (4.19) we have
ddt∥w∥2+a∥w∥2≤C∥v∥4125+C∥v∥2(β+1)65(β+1)+1a∥f0∥2. | (4.24) |
Case Ⅰ. 1<β<43.
Now we use Gagliardo-Nirenberg inequality to obtain
∥v∥125≤C∥∇v∥1/4∥v∥3/4, | (4.25) |
∥v∥65(β+1)≤C∥∇v∥3β−22(β+1)∥v∥4−β2(β+1). | (4.26) |
Considering 1<β<43, so 3β−2<2. Combining (4.25), (4.26) with (4.24), and according to (4.16), we obtain
ddt∥w∥2+a∥w∥2≤C∥∇v∥∥v∥3+C∥∇v∥3β−2∥v∥4−β+1a∥f0∥2≤∥∇v∥2+C∥v∥6+C∥v∥2(4−β)4−3β+1a∥f0∥2≤∥∇v∥2+Cl6ε6(1−ρ)+C(l2ε2(1−ρ))4−β4−3β+1a∥f0∥2. | (4.27) |
Let g(s)=∥∇v(s)∥2+Cl6ε6(1−ρ)+C(l2ε2(1−ρ))4−β4−3β+1a∥f0∥2. Noticing (4.16), we have
∫t+1tg(s)ds=∫t+1t[∥∇v(s)∥2+Cl6ε6(1−ρ)+C(l2ε2(1−ρ))4−β4−3β+1a∥f0∥2]ds≤C(l2ε2(1−ρ)+l6ε6(1−ρ)+(l2ε2(1−ρ))4−β4−3β+M20),∀t≥τ. | (4.28) |
Applying Lemma 2.3 to (4.27) we have
∥w(t)∥2≤∥uτ∥2e−a(t−τ)+C(1+1a)(l2ε2(1−ρ)+l6ε6(1−ρ)+(l2ε2(1−ρ))4−β4−3β+M20)≤∥uτ∥2e−a(t−τ)+C(l2+l6+l2(4−β)4−3β+M20),∀t≥τ. | (4.29) |
Case Ⅱ. β=43.
From (4.24)-(4.26) we have
ddt∥w∥2+a∥w∥2≤C∥∇v∥∥v∥3+C∥∇v∥2∥v∥83+1a∥f0∥2≤∥∇v∥2+C∥v∥6+C∥∇v∥2∥v∥83+1a∥f0∥2≤[1+C∥v∥83]∥∇v∥2+C∥v∥6+1a∥f0∥2≤(1+Cl83ε83(1−ρ))∥∇v∥2+Cl6ε6(1−ρ)+1a∥f0∥2. | (4.30) |
Similar to the derivation of (4.29), we get
∥w(t)∥2≤∥uτ∥2e−a(t−τ)+C(1+1a)[(1+l83ε83(1−ρ))l2ε2(1−ρ)+l6ε6(1−ρ)+M20],≤∥uτ∥2e−a(t−τ)+C(l2+l143+l6+M20),∀t≥τ. | (4.31) |
Recalling that u=w+v, using (4.16), (4.29) and (4.31), we end up with
∥u(t)∥2≤∥uτ∥2e−a(t−τ)+C(l2+l6+l2(4−β)4−3β+l143+M20),∀t≥τ. | (4.32) |
Thus, for every ε≤ε0, the process {Ufε(t,τ)} has the absorbing set
B0:={u∈H|∥u∥2≤C(l2+l6+l2(4−β)4−3β+l143+M20)}. |
On the other hand, if ε0<ε≤1, the process {Ufε(t,τ)} possesses also the absorbing set (cf (3.2))
Bε0={u∈H|∥u∥≤C0Qε0}. |
In conclusion, for every ε∈[0,1], the bounded set
B∗=B0∪Bε0 |
is an absorbing set for {Ufε(t,τ)} which is independent of ε. Since Aε⊂B∗, the proof is completed.
The main result of this section is the following.
Theorem 5.1. Let (4.11) hold. Then, the uniform attractor Aε converges to A0 as ε→0+ in the following sense:
limε→0+distH(Aε,A0)=0. |
The proof of this theorem requires some steps. Now, we shall study the difference of two solutions to (3.1) with ε>0 and ε=0, respectively, sharing the same initial data. We denote
uε(t)=Ufε(t,τ)uτ, |
with uτ belonging to the absorbing set B∗ found in the previous section. Owing to (4.32), we have the uniform bound:
∥uε(t)∥2≤R21, | (5.1) |
for some R1=R1(l,M0) because the size of B∗ depends on l and M0. In particular, for ε=0, since uτ∈B∗, we have the bound
∥u0(t)∥2≤R20, | (5.2) |
for some R0=R0(l,M0).
Lemma 5.1. For every ε∈(0,1], every τ∈R and every uτ∈B∗, the deviation ˜w(t)=uε(t)−u0(t) with uε(0)=u0(0)=uτ, fulfills the estimate
∥˜w(t)∥2≤Cl2ε2(1−ρ),∀t≥τ, | (5.3) |
for some positive constant C independent of ε.
Proof. Since the deviation ˜w(t) solves
˜wt−νΔ˜w+a˜w+b|uε|uε−b|u0|u0+c|uε|βuε−c|u0|βu0=ε−ρf1(x,tε),˜w|t=τ=0, | (5.4) |
the difference q(t)=˜w(t)−v(t), where v(t) is the solution to (4.15), fulfills the Cauchy problem
qt−νΔq+aq+b|uε|uε−b|u0|u0+c|uε|βuε−c|u0|βu0=0,q|t=τ=0. | (5.5) |
At this point, we take the scalar product in H of (5.5) with q, so getting
12ddt∥q∥2+ν∥∇q∥2+a∥q∥2+b(|uε|uε−|u0|u0,˜w)+c(|uε|βuε−|u0|βu0,˜w)=b(|uε|uε−|u0|u0,v)+c(|uε|βuε−|u0|βu0,v). | (5.6) |
Noting the first term on the right-hand side of (5.6) is given by
b(|uε|uε−|u0|u0,v)=b(|uε|˜w,v)+b((|uε|−|u0|)u0,v), | (5.7) |
we now proceed to estimate the first term on the right-hand side of (5.7). Since
b(|uε|˜w,v)≤b∫Ω|uε||˜w||v|dx≤b∫Ω|uε|(|v|+|q|)|v|dx≤b∫Ω|uε||v|2dx+b∫Ω|uε||q||v|dx, | (5.8) |
and
b∫Ω|uε||q||v|dx≤b(∫Ω|q|6dx)16(∫Ω|uε|2dx)12(∫Ω|v|3dx)13=b∥q∥6∥uε∥∥v∥3≤ν4d20∥q∥26+C∥uε∥2∥v∥23≤ν4∥∇q∥2+C∥uε∥2∥∇v∥2, | (5.9) |
b∫Ω|uε||v|2dx≤b(∫Ω|uε|2dx)12(∫Ω|v|6dx)16(∫Ω|v|3dx)13=b∥uε∥∥v∥6∥v∥3≤12d20∥v∥26+C∥uε∥2∥v∥23≤12∥∇v∥2+C∥uε∥2∥∇v∥2, | (5.10) |
it follows from (5.8)-(5.10) that
b(|uε|˜w,v)≤ν4∥∇q∥2+12∥∇v∥2+C∥uε∥2∥∇v∥2. | (5.11) |
Now, let us estimate the second term on the right-hand side of (5.7). Noting
b((|uε|−|u0|)u0,v)≤b∫Ω|uε−u0||u0||v|dx=b∫Ω|˜w||u0||v|dx≤b∫Ω(|q|+|v|)|u0||v|dx=b∫Ω|u0||q||v|dx+b∫Ω|v|2|u0|dx, | (5.12) |
similar arguments as (5.9) and (5.10), we have
b∫Ω|u0||q||v|dx≤ν4∥∇q∥2+C∥u0∥2∥∇v∥2, | (5.13) |
and
b∫Ω|v|2|u0|dx≤12∥∇v∥2+C∥u0∥2∥∇v∥2. | (5.14) |
Hence, from (5.12)-(5.14) we get
b((|uε|−|u0|)u0,v)≤ν4∥∇q∥2+12∥∇v∥2+C∥u0∥2∥∇v∥2. | (5.15) |
Combining (5.11), (5.15) with (5.7), we have
b(|uε|uε−|u0|u0,v)≤ν2∥∇q∥2+∥∇v∥2+C(∥uε∥2+∥u0∥2)∥∇v∥2. | (5.16) |
Noting the second term on the right-hand side of (5.6) is given by
c(|uε|βuε−|u0|βu0,v)=c(|uε|β˜w,v)+c((|uε|β−|u0|β)u0,v), | (5.17) |
we now proceed to estimate the first term on the right-hand side of (5.17). Since
c(|uε|β˜w,v)≤c∫Ω|uε|β(|q|+|v|)|v|dx≤c∫Ω|uε|β|q||v|dx+c∫Ω|uε|β|v|2dx, | (5.18) |
so we should estimate the right-hand side of the last inequality in (5.18) term by term. Because
c∫Ω|uε|β|q||v|dx≤c(∫Ω|q|6dx)16(∫Ω|uε|2dx)β2(∫Ω|v|65−3βdx)5−3β6=c∥q∥6∥uε∥β∥v∥65−3β≤ν4d20∥q∥26+C∥uε∥2β∥v∥265−3β≤ν4∥∇q∥2+C∥uε∥2β∥∇v∥2 | (5.19) |
and
c∫Ω|uε|β|v|2dx≤c(∫Ω|v|6dx)16(∫Ω|uε|2dx)β2(∫Ω|v|65−3βdx)5−3β6=c∥v∥6⋅∥uε∥β⋅∥v∥65−3β≤12∥∇v∥2+C∥uε∥2β∥∇v∥2, | (5.20) |
where the last inequalities in (5.19) and (5.20) are valid only if 65−3β≤6, i.e. β≤43, so combining (5.19), (5.20) with (5.18), we have
c(|uε|β˜w,v)≤ν4∥∇q∥2+12∥∇v∥2+C∥uε∥2β∥∇v∥2. | (5.21) |
Now let us estimate the second term on the right-hand side of (5.17). Since
c((|uε|β−|u0|β)u0,v)≤c∫Ω||uε|β−|u0|β||u0||v|dx≤C∫Ω||uε|β−1+|u0|β−1||˜w||u0||v|dx≤C∫Ω|uε|β−1|˜w||u0||v|dx+C∫Ω|u0|β−1|˜w||u0||v|dx, | (5.22) |
in the second inequality of (5.22) we used the fact that
|xp−yp|≤Cp(xp−1+yp−1)|x−y| |
for any x,y≥0, where C is an absolute constant. So let us estimate the right-hand side of the last inequality in (5.22) term by term. For the first term, we have
C∫Ω|uε|β−1|˜w||u0||v|dx≤C∫Ω|uε|β−1|q||u0||v|dx+C∫Ω|uε|β−1|v|2|u0|dx, | (5.23) |
and
C∫Ω|uε|β−1|q||u0||v|dx≤C(∫Ω|q|6dx)16(∫Ω|uε|2dx)β−12(∫Ω|u0|2dx)12(∫Ω|v|65−3β)5−3β6=C∥q∥6∥uε∥β−1∥u0∥∥v∥65−3β≤ν8∥∇q∥2+C∥uε∥2(β−1)∥u0∥2∥v∥265−3β≤ν8∥∇q∥2+C∥uε∥2(β−1)∥u0∥2∥∇v∥2, | (5.24) |
similarly,
C∫Ω|uε|β−1|v|2|u0|dx≤14∥∇v∥2+C∥uε∥2(β−1)∥u0∥2∥∇v∥2. | (5.25) |
It follows from (5.23)-(5.25) that
C∫Ω|uε|β−1|˜w||u0||v|dx≤ν8∥∇q∥2+14∥∇v∥2+C∥uε∥2(β−1)∥u0∥2∥∇v∥2. | (5.26) |
Similar arguments as (5.23)-(5.26), for the second term on the right-hand side of inequality (5.22), we have
C∫Ω|u0|β−1|˜w||u0||v|dx≤ν8∥∇q∥2+14∥∇v∥2+C∥u0∥2(β−1)∥u0∥2∥∇v∥2. | (5.27) |
Combining (5.26), (5.27) with (5.22), we have
c((|uε|β−|u0|β)u0,v)≤ν4∥∇q∥2+12∥∇v∥2+C(∥uε∥2(β−1)+∥u0∥2(β−1))∥u0∥2∥∇v∥2. | (5.28) |
So it follows from (5.17), (5.21) and (5.28) that
c(|uε|βuε−|u0|βu0,v)≤ν2∥∇q∥2+∥∇v∥2+C∥uε∥2β∥∇v∥2 +C(∥uε∥2(β−1)+∥u0∥2(β−1))∥u0∥2∥∇v∥2. | (5.29) |
Now, considering (5.6), (5.16) and (5.29), and according to (5.1) and (5.2), it yields
ddt∥q∥2+2a∥q∥2≤4∥∇v∥2+C(R21+R20)∥∇v∥2+CR2β1∥∇v∥2 +C(R2(β−1)1+R2(β−1)0)R20∥∇v∥2≤C(1+R20+R21+R2β1+(R2(β−1)0+R2(β−1)1)R20)∥∇v∥2. | (5.30) |
Recalling that ∥q(τ)∥=0 and (4.16), Lemma 2.3 entails
∥q(t)∥2≤C(1+12a)(1+R20+R21+R2β1+(R2(β−1)0+R2(β−1)1)R20)l2ε2(1−ρ)≤Cl2ε2(1−ρ). |
Finally, as ˜w=q+v, using (4.16) to control ∥v∥, we obtain the desired conclusion (5.3).
In order to study the convergence of the uniform attractors, we actually need a generalization of Lemma 5.1, which applies to the whole family of equations
ˆut+νAˆu+aˆu+G(ˆu)=ˆfε,ˆfε∈H(fε), | (5.31) |
with the external force ˆf=ˆfε∈H(fε). To this end, we observe that every function ˆf1∈H(f1) fulfills the inequality (4.11). Defining
ˆF1(t,τ)=∫tτˆf1(s)ds,t≥τ, |
we have
supt≥τ,τ∈R{∥ˆF1(t,τ)∥2+∫t+1t∥ˆF1(s,τ)∥2Vds}≤l2. | (5.32) |
For any ε∈[0,1], let ˆuε(t)=Uˆfε(t,τ)uτ be the solution to (5.31) with external force ˆfε=ˆf0+ε−ρˆf1(⋅/ε)∈H(fε) and uτ∈B∗. For ε>0, we consider the deviation ˆw(t)=ˆuε(t)−ˆu0(t).
Lemma 5.2. The inequality
∥ˆw(t)∥2≤Cl2ε2(1−ρ),∀t≥τ, | (5.33) |
holds, where C is independent of ε.
Proof. As the similar argument to the proof of Lemma 5.1, with ˆuε,ˆf0 and ˆf1 in place of uε,f0 and f1, respectively. Noting that (5.2) still holds for ˆu0, and the family {Uˆfε(t,τ)}(ˆfε∈H(fε)) is (H×H(fε),H)-continuous, and using (5.32) in place of (4.11), finally complete the proof of the lemma.
We can now complete the proof of Theorem 5.1, using the following argument from [27], which we report in some detail for the reader's convenience.
Proof of Theorem 5.1 Let ε>0 and uε∈Aε. Thus, in view of (3.18), there exists a complete bounded trajectory ˆuε(t) of (5.31), with the external force
ˆfε=ˆf0+ε−ρˆf1(⋅/ε)∈H(fε),ˆf0∈H(f0),ˆf1∈H(f1) |
such that ˆuε(0)=uε. For every L≥0 to be specified later, consider the vector
ˆuε(−L)∈Aε⊂B∗. |
From the straightforward equality
uε=Uˆfε(0,−L)ˆuε(−L), |
by applying Lemma 5.2, we have that
∥uε−Uˆf0(0,−L)ˆuε(−L)∥≤Clε1−ρ. | (5.34) |
On the other hand, the set A0 attracts Uˆf0(t,−L)B∗, uniformly as ˆf0∈H(f0). Then, for every δ>0, there is T=T(δ)≥0, independent of L, such that
distH(Uˆf0(T−L,−L)ˆuε(−L),A0)≤δ. | (5.35) |
Setting L=T, and collecting the two above inequalities, we readily get
distH(uε,A0)≤Clε1−ρ+δ. |
Since uε∈Aε and δ>0 are arbitrary, taking the limit ε→0+, the conclusion follows.
In this paper, we investigated a class of three-dimensional Brinkman-Forchheimer equation with oscillating external forces fε(x,t)=f0(x,t)+ε−ρf1(x,tε). Based on some translation-compactness assumptions on the external forces, we obtained the uniform boundedness of the uniform attractor Aε of the system (1.1) in (L2(Ω))3, and the convergence of Aε to the attractor A0 of the system (1.2) as ε→0+. To prove the uniform boundedness and the convergence of the uniform attractors, the Gagliardo-Nirenberg inequality is needed. In the proof process, we concluded that the parameter β∈(1,43].
The authors are thankful to the editors and the anonymous reviewers for their valuable suggestions and comments on the manuscript. The first author was supported in part by the National Natural Science Foundation of China (No. 11601417), Natural Science Basic Research Plan in Shaanxi Province of China (Nos. 2018JM1047, 2019JM-283) and Postdoctoral Fund in Shaanxi Province of China (No. 2016BSHEDZZ112). The second author was supported by the National Natural Science Foundation of China (No. 11771262).
The authors declare no conflict of interest in this paper.
[1] |
R. C. Gilver, S. A. Altobelli, A determination of effective viscosity for the Brinkman-Forchheimer flow model, J. Fluid Mech., 258 (1994), 355-370. doi: 10.1017/S0022112094003368
![]() |
[2] |
D. A. Nield, The limitations of the Brinkman-Forchheimer equation in modeling flow in a saturated porous medium and at an interface, Int. J. Heat Fluid Fl., 12 (1991), 269-272. doi: 10.1016/0142-727X(91)90062-Z
![]() |
[3] |
K. Vafai, S. J. Kim, Fluid mechanics of the interface region between a porous medium and a fluid layer-an exact solution, Int. J. Heat Fluid Fl., 11 (1990), 254-256. doi: 10.1016/0142-727X(90)90045-D
![]() |
[4] |
K. Vafai, S. J. Kim, On the limitations of the Brinkman-Forchheimer-extended Darcy equation, Int. J. Heat Fluid Fl., 16 (1995), 11-15. doi: 10.1016/0142-727X(94)00002-T
![]() |
[5] |
S. Whitaker, The Forchheimer equation: a theoretical development, Transp. Porous Media, 25 (1996), 27-62. doi: 10.1007/BF00141261
![]() |
[6] |
M. M. Bhatti, A. Zeeshan, R. Ellahi, et al. Mathematical modelling of heat and mass transfer effects on MHD peristaltic propulsion of two-phase flow through a Darcy-Brinkman-Forchheimer porous medium, Adv. Powder Technol., 29 (2018), 1189-1197. doi: 10.1016/j.apt.2018.02.010
![]() |
[7] | M. Marin, S. Vlase, R. Ellahi, et al. On the Partition of Energies for the Backward in Time Problem of Thermoelastic Materials with a Dipolar Structure, Symmetry, 11 (2019), 863. |
[8] | H. M. Zhang, C. Yuan, G. S. Yang, et al. A novel constitutive modelling approach measured under simulated freeze-thaw cycles for the rock failure, Engineering with Computers, 2019. |
[9] |
A. O. Celebi, V. Kalantarov, Continuous dependence for the covective Brinkman-Forchheimer equations, Appl. Anal., 84 (2005), 877-888. doi: 10.1080/00036810500148911
![]() |
[10] |
A. O. Celebi, V. Kalantarov, D. Uğurlu, On continuous dependence on coefficients of the BrinkmanForchheimer equations, Appl. Math. Lett., 19 (2006), 801-807. doi: 10.1016/j.aml.2005.11.002
![]() |
[11] |
F. Franchi, B. Straughan, Continuous dependence and decay for the Forchheimer equations, Proc. R. Soc. Lond. A, 459 (2003), 3195-3202. doi: 10.1098/rspa.2003.1169
![]() |
[12] |
L. E. Payne, B. Straughan, Convergence and continuous dependence for the BrinkmanForchheimer equations, Stud. Appl. Math., 102 (1999), 419-439. doi: 10.1111/1467-9590.00116
![]() |
[13] |
Y. Liu, Convergence and continuous dependence for the Brinkman-Forchheimer equations, Math. Comput. Model., 49 (2009), 1401-1415. doi: 10.1016/j.mcm.2008.11.010
![]() |
[14] |
Y. Liu, S. Z. Xiao, Y. W. Lin, Continuous dependence for the Brinkman-Forchheimer fluid interfacing with a Darcy fluid in a bounded domain, Math. Comput. Simulat., 150 (2018), 66-82. doi: 10.1016/j.matcom.2018.02.009
![]() |
[15] | Y. F. Li, C. H. Lin, Continuous dependence for the nonhomogeneous Brinkman-Forchheimer equations in a semi-infinite pipe, Appl. Math. Comput., 244 (2014), 201-208. |
[16] |
D. Uğurlu, On the existence of a global attractor for the Brinkman-Forchheimer equation, Nonlinear Anal., 68 (2008), 1986-1992. doi: 10.1016/j.na.2007.01.025
![]() |
[17] |
B. Wang, S. Lin, Existence of global attractors for the three-dimensional Brinkman-Forchheimer equation, Math. Method. Appl. Sci., 31 (2008), 1479-1495. doi: 10.1002/mma.985
![]() |
[18] |
Y. You, C. Zhao, S. Zhou, The existence of uniform attractors for 3D Brinkman-Forchheimer equations, Discrete Contin. Dyn. Syst., 32 (2012), 3787-3800. doi: 10.3934/dcds.2012.32.3787
![]() |
[19] | X. Song, Pullback D-attractors for a non-autonomous Brinkman-Forchheimer system, J. Math. Res. Application, 33 (2013), 90-100. |
[20] | X. Song, Q. Bao, Uniform attractors for three-dimensional Brinkman-Forchheimer system and some averaging problems, Far East J. Dyn. Syst., 25 (2014), 99-122. |
[21] | L. Zhang, K. Su, S. Wen, Attractors for the 3D autonomous and nonautunomous BrinkmanForchheimer equations, Bound. Value. Probl., 17 (2016), 1-18. |
[22] |
Y. Ouyang, L. Yang, A note on the existence of a global attractor for the Brinkman-Forchheimer equations, Nonlinear Anal., 70 (2009), 2054-2059. doi: 10.1016/j.na.2008.02.121
![]() |
[23] |
V. V. Chepyzhov, M. I. Vishik, W. L. Wendland, On non-autonomous sine-Gorden type equations with a simple global attractor and some averaging, Discrete Contin. Dyn. Syst., 12 (2005), 27-38. doi: 10.3934/dcds.2005.12.27
![]() |
[24] | M. Efendiev, S. V. Zelik, Attractors of the reaction-diffusion systems with rapidly oscillating coefficients and their homogenenization, Ann. I. H. Poincaré-An, 19 (2002), 961-989. |
[25] | M. Efendiev, S. V. Zelik, The regular attractor for the reaction-diddusion systems with rapidly oscillating in time and its averaging, Adv. Differential Equ., 8 (2003), 673-732. |
[26] |
V. V. Chepyzhov, M. I. Vishik, Non-autonomous 2D Navier-Stokes system with singularly oscillating external force and its global attractor, J. Dyn. Differ. Equ., 19 (2007), 655-684. doi: 10.1007/s10884-007-9077-y
![]() |
[27] |
V. V. Chepyzhov, V. Pata, M. I. Vishik, Averaging of 2D Navier-Stokes equations with singularly oscillating forces, Nonlinearity, 22 (2009), 351-370. doi: 10.1088/0951-7715/22/2/006
![]() |
[28] |
V. V. Chepyzhov, M. I. Vishik, Non-autonomous Navier-Stokes system with a simple global attractor and some averaging problems, ESAIM: Control, Optimisation and Calculus of Variations, 8 (2002), 467-487. doi: 10.1051/cocv:2002056
![]() |
[29] |
Y. Qin, X. Yang, X. Liu, Averaging of a 3D Navier-Stokes-Voight equation with singularly oscillating forces, Nonlinear Anal-Real, 13 (2012), 893-904. doi: 10.1016/j.nonrwa.2011.08.025
![]() |
[30] |
C. T. Anh, N. D. Toan, Nonclassical diffusion equations on RN with singularly oscillating external forces, Appl. Math. Lett., 38 (2014), 20-26. doi: 10.1016/j.aml.2014.06.008
![]() |
[31] |
T. T. Medjo, Non-autonomous planetary 3D geostrophic equations with oscillating external force and its global attractor, Nonlinear Anal-Real, 12 (2011), 1437-1452. doi: 10.1016/j.nonrwa.2010.10.004
![]() |
[32] | T. T. Medjo, Averaging of the planetary 3D geostrophic equations with oscillating external forces, Appl. Math. Comput., 218 (2012), 5910-5928. |
[33] |
T. T. Medjo, A non-autonomous two-phase flow model with oscillating external force and its global attractor, Nonlinear Anal., 75 (2012), 226-243. doi: 10.1016/j.na.2011.08.024
![]() |
[34] |
M. Bigert, A priori estimates for the difference of solutions to quasi-linear elliptic equations, Manuscripta Math., 133 (2010), 273-306. doi: 10.1007/s00229-010-0367-z
![]() |
[35] | P. Constantin, C. Foias, Navier-Stokes equations, Chicago and London: Univ. Chicago Press, 1989. |
[36] | R. Temam, Infinite Dimensional Dynamical Systems in Mechanics and Physics, second ed., in: Appl. Math. Sci., vol. 68, New York: Springe-Verlag, 1988. |
1. | Lei Zhang, Ping Gao, Exploring the Independent Learning of Music Students in Higher Teacher Training Colleges in a Pluralistic Network Environment, 2022, 2022, 1687-9813, 1, 10.1155/2022/5024006 | |
2. | Zhengwang Tao, Xin-Guang Yang, Yan Lin, Chunxiao Guo, Determination of Three-Dimensional Brinkman—Forchheimer-Extended Darcy Flow, 2023, 7, 2504-3110, 146, 10.3390/fractalfract7020146 | |
3. | Chengfei Ai, Zhong Tan, Gromov–Hausdorff stability of global attractors for 3D Brinkman–Forchheimer equations, 2022, 45, 0170-4214, 11117, 10.1002/mma.8440 | |
4. | Xueli SONG, Xi DENG, Baoming QIAO, Dimension Estimate of the Global Attractor for a 3D Brinkman- Forchheimer Equation, 2023, 28, 1007-1202, 1, 10.1051/wujns/2023281001 | |
5. | Xueli SONG, Shuang XU, Baoming QIAO, Existence of Global Attractor for a 3D Brinkman-Forchheimer Equfation in Some Poincaré Unbounded Domains, 2023, 28, 1007-1202, 282, 10.1051/wujns/2023284282 | |
6. | Qiaozhen Ma, Lulu Wang, Uniform attractors of non-autonomous wave equations with singularly oscillating external forces and displacement-dependent damping, 2025, 426, 00220396, 388, 10.1016/j.jde.2025.01.061 |