Research article

Non-autonomous 3D Brinkman-Forchheimer equation with singularly oscillating external force and its uniform attractor

  • Received: 10 October 2019 Accepted: 09 January 2020 Published: 21 January 2020
  • MSC : 35B40, 35B41, 35Q35

  • For ρ[0,1) and ε>0, the non-autonomous 3D Brinkman-Forchheimer equation with singularly oscillating external force utνΔu+au+b|u|u+c|u|βu+p=f0(x,t)+ερf1(x,tε),divu=0 are considered, together with the averaged equation utνΔu+au+b|u|u+c|u|βu+p=f0(x,t),divu=0 formally corresponding to the limiting case ε=0. First, within the restriction ρ<1 and under suitable translation-compactness assumptions on the external forces, the uniform (w.r.t.ε) boundedness of the related uniform attractors Aε is established when 1<β4/3. This fact is not at all intuitive, since in principle the blow up of the oscillation amplitude might overcome the averaging effect due to the term tε in f1. Next, the convergence of the attractor Aε of the first equation to the attractor A0 of the second one as ε0+ is established.

    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

    Related Papers:

    [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, the non-autonomous 3D Brinkman-Forchheimer equation with singularly oscillating external force utνΔu+au+b|u|u+c|u|βu+p=f0(x,t)+ερf1(x,tε),divu=0 are considered, together with the averaged equation utνΔu+au+b|u|u+c|u|βu+p=f0(x,t),divu=0 formally corresponding to the limiting case ε=0. First, within the restriction ρ<1 and under suitable translation-compactness assumptions on the external forces, the uniform (w.r.t.ε) boundedness of the related uniform attractors Aε is established when 1<β4/3. This fact is not at all intuitive, since in principle the blow up of the oscillation amplitude might overcome the averaging effect due to the term tε in f1. Next, the convergence of the attractor Aε of the first equation to the attractor A0 of the second one as ε0+ is established.


    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(N3) 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
    ufluid velocity vector(m/s)ppressure
    R3three-dimensional whole spaceΩopen, bounded domain of R3
    u1,u2,u3velocity componentsttime
    Ωboundary of a domainνBrinkman coefficient
    aDarcy coefficientb,cForchheimer coefficient
    βpower of nonlinear termfexternal force

    Given a space X, we usually denote the norm in X by X, and we indicate by

    distX(B1,B2)=supb1B1infb2B2b1b2X,

    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(C0(Ω))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(Ω))3V, 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,uHσ:=∥Aσ2u(L2(Ω))3.

    In particular,

    H0=H,H1=V,H2=D(A),

    and we have the generalized Poincarˊe inequality

    uHσ+1λ121uHσ,uHσ+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)=Ωuvdx,u,vH, u∥=(u,u)12,

    ((,)) and V denote the product and norm of V, i.e.,

    ((u,v))=3i=1Ωuividx,u,vV,uV=((u,u))12.

    In this paper, Lp(Ω)=(Lp(Ω))3, and we use p to denote the norm in Lp(Ω) (for p2). 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

    f02L2b:=suptRt+1tf0(s)2ds=M20, (2.2)
    f12L2b:=suptRt+1tf1(s)2ds=M21, (2.3)

    for some constants M0,M10.

    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εL2bQε:={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,bR,

    (|a|p2a|b|p2b)(ab)22p3p/2|ab|p.

    Lemma 2.2. ([27]) For every τR, every nonnegative locally summable function φ on Rτ and every β>0, we have

    tτφ(s)eβ(ts)ds11eβsupθτθ+1θφ(s)ds, (2.5)

    for all tτ.

    Lemma 2.3. ([23]) Let a real function z(t),t0 be uniformly continuous and satisfy the inequality

    dzdt+γz(t)f(t),t0,

    where γ>0,f(t)0 for all t0 and fL1loc(R+). Suppose also that

    t+1tf(s)dsM,t0.

    Then

    z(t)z(0)eγt+M(1+γ1),t0.

    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

    uC(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ε={uH|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)|hR}]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 ˆf0H(f0) and ˆf1H(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])

    ˆf0L2b(R;H)≤∥f0L2b(R;H),ˆf0H(f0), (3.5)
    ˆf1L2b(R;H)≤∥f1L2b(R;H),ˆf1H(f1). (3.6)

    It follows that

    ˆfεL2b(R;H)≤∥f0L2b(R;H)+Cερf1L2b(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ε|HC+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τnuτ in H, (3.10)

    then

    Uˆfεn(t,τ)uτnUˆ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τnUˆ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

    12ddtwn2+νwn2+awn2+(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|βˆunb|ˆu|ˆuc|ˆu|βˆu,ˆunˆu)0, (3.14)

    then (3.13) gives

    ddtwn2+awn21aˆfεnˆfε2. (3.15)

    Applying Gronwall Lemma to (3.15) we have

    wn(t)2≤∥wn(τ)2ea(tτ)+1atτˆfεnˆfε2ea(ts)ds≤∥wn(τ)2+1atτˆfεnˆfε2ds,tτ. (3.16)

    Note that

    ˆfεnˆfε in L2loc(R;H) and uτnuτ 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),tR} of (3.8), which are uniformly bounded in H. The set

    Kˆfε(s)={ˆu(s)|ˆuKˆ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 KL2loc(R;V), then the above problem has a unique solution

    VC(Rτ;H2)L2loc(Rτ;H3).

    Moreover, the inequalities

    V(t)2Ctτea(ts)K(s)2ds, (4.2)
    AV(t)2Ctτe2a(ts)K(s)2Vds, (4.3)
    t+1tV(s)2dsC(V(t)2+t+1tK(s)2ds) (4.4)
    t+1tA32V(s)2dsC(AV(t)2+t+1tK(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

    12ddtV2+νV2+aV2=(K(t),V)a2V2+12aK(t)2, (4.6)
    12ddtAV2+νA32V2+aAV2=(K(t),A2V)12νK(t)2V+ν2A32V2. (4.7)

    It follows from (4.6) and (4.7) that

    ddtV2+aV21aK(t)2,
    ddtAV2+2aAV21νK(t)2V.

    Applying Gronwall lemma it yields

    V(t)21atτea(ts)K(s)2ds,
    AV(t)21νtτe2a(ts)K(s)2Vds.

    From (4.6) and (4.7) we also can get

    ddtV2+2νV21aK(t)2, (4.8)
    ddtAV2+νA32V21νK(t)2V. (4.9)

    Integrating (4.8) and (4.9) on [t,t+1] respectively, we obtain

    2νt+1tV(s)2ds≤∥V(t)2+1at+1tK(s)2ds, (4.10)
    νt+1tA32V(s)2ds≤∥AV(t)2+1νt+1tK(s)2Vds.

    The proof is finished.

    Setting F(t,τ)=tτf1(s)ds,tτ, we assume that

    suptτ,τR(F(t,τ)2+t+1tF(s,τ)2Vds)l2. (4.11)

    Lemma 4.2. Assume that f1L2loc(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+1tv(s)2dsCl2ε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 t0,

    tV(t)=v(t)=t0tv(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

    supt0Fε(t)∥≤lε

    and

    t+1tFε(s)2Vds=ε3t+1εtεF(s,0)2Vds2ε2supt0{t+1tF(s,0)2Vds}2l2ε2.

    By (2.5) we have

    t0ea(ts)Fε(s)2dsCl2ε2,t0e2a(ts)Fε(s)2VdsCl2ε2.

    So applying Lemma 5.1, we obtain

    V(t)2+AV(t)2+t+1tV(s)2ds+t+1tA32V(s)2dsCl2ε2.

    Hence, on account of (4.14) we have

    v(t)∥=∥tV(t)∥≤∥Fε(t)+νAV(t)+aV(t)∥≤Clε

    and

    v(s)2=∥(tV(s))23Fε(s)2V+3ν2A32V(s)2+3a2V(s)2,

    from which we derive the integral estimate

    t+1tv(s)2dsCl2ε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+1tv(s)2dscl2ε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

    12ddtw2+νw2+aw2+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

    12ddtw2+νw2+aw2+b25332w33+c24β3β+22wβ+2β+2b(|v|v,w)c(|v|βv,w)+(f0,w). (4.19)

    Noticing

    (f0,w)a2w2+12af02, (4.20)
    b|(|v|v,w)|=b|Ω|v|vwdx|b(Ω|w|6dx)1/6(Ω(|v|v)6/5dx)5/6=bw6v2125ν2d20w26+Cv4125, (4.21)

    and

    c|(|v|βv,w)|=c|Ω|v|βvwdx|c(Ω|w|6dx)16(Ω(|v|βv)6/5dx)5/6=cw6vβ+165(β+1)ν2d20w26+Cv2(β+1)65(β+1), (4.22)

    according to Sobolev inequality

    vpd0v,1p6, (4.23)

    combining (4.20)-(4.23) with (4.19) we have

    ddtw2+aw2Cv4125+Cv2(β+1)65(β+1)+1af02. (4.24)

    Case Ⅰ. 1<β<43.

    Now we use Gagliardo-Nirenberg inequality to obtain

    v125Cv1/4v3/4, (4.25)
    v65(β+1)Cv3β22(β+1)v4β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

    ddtw2+aw2Cv∥∥v3+Cv3β2v4β+1af02≤∥v2+Cv6+Cv2(4β)43β+1af02≤∥v2+Cl6ε6(1ρ)+C(l2ε2(1ρ))4β43β+1af02. (4.27)

    Let g(s)=∥v(s)2+Cl6ε6(1ρ)+C(l2ε2(1ρ))4β43β+1af02. Noticing (4.16), we have

    t+1tg(s)ds=t+1t[v(s)2+Cl6ε6(1ρ)+C(l2ε2(1ρ))4β43β+1af02]dsC(l2ε2(1ρ)+l6ε6(1ρ)+(l2ε2(1ρ))4β43β+M20),tτ. (4.28)

    Applying Lemma 2.3 to (4.27) we have

    w(t)2≤∥uτ2ea(tτ)+C(1+1a)(l2ε2(1ρ)+l6ε6(1ρ)+(l2ε2(1ρ))4β43β+M20)≤∥uτ2ea(tτ)+C(l2+l6+l2(4β)43β+M20),tτ. (4.29)

    Case Ⅱ. β=43.

    From (4.24)-(4.26) we have

    ddtw2+aw2Cv∥∥v3+Cv2v83+1af02≤∥v2+Cv6+Cv2v83+1af02[1+Cv83]v2+Cv6+1af02(1+Cl83ε83(1ρ))v2+Cl6ε6(1ρ)+1af02. (4.30)

    Similar to the derivation of (4.29), we get

    w(t)2≤∥uτ2ea(tτ)+C(1+1a)[(1+l83ε83(1ρ))l2ε2(1ρ)+l6ε6(1ρ)+M20],≤∥uτ2ea(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τ2ea(tτ)+C(l2+l6+l2(4β)43β+l143+M20),tτ. (4.32)

    Thus, for every εε0, the process {Ufε(t,τ)} has the absorbing set

    B0:={uH|u2C(l2+l6+l2(4β)43β+l143+M20)}.

    On the other hand, if ε0<ε1, the process {Ufε(t,τ)} possesses also the absorbing set (cf (3.2))

    Bε0={uH|u∥≤C0Qε0}.

    In conclusion, for every ε[0,1], the bounded set

    B=B0Bε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)2R21, (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)2R20, (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)2Cl2ε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

    12ddtq2+νq2+aq2+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|dxbΩ|uε|(|v|+|q|)|v|dxbΩ|uε||v|2dx+bΩ|uε||q||v|dx, (5.8)

    and

    bΩ|uε||q||v|dxb(Ω|q|6dx)16(Ω|uε|2dx)12(Ω|v|3dx)13=bq6uε∥∥v3ν4d20q26+Cuε2v23ν4q2+Cuε2v2, (5.9)
    bΩ|uε||v|2dxb(Ω|uε|2dx)12(Ω|v|6dx)16(Ω|v|3dx)13=buε∥∥v6v312d20v26+Cuε2v2312v2+Cuε2v2, (5.10)

    it follows from (5.8)-(5.10) that

    b(|uε|˜w,v)ν4q2+12v2+Cuε2v2. (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|dxbΩ(|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ν4q2+Cu02v2, (5.13)

    and

    bΩ|v|2|u0|dx12v2+Cu02v2. (5.14)

    Hence, from (5.12)-(5.14) we get

    b((|uε||u0|)u0,v)ν4q2+12v2+Cu02v2. (5.15)

    Combining (5.11), (5.15) with (5.7), we have

    b(|uε|uε|u0|u0,v)ν2q2+v2+C(uε2+u02)v2. (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|dxcΩ|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|dxc(Ω|q|6dx)16(Ω|uε|2dx)β2(Ω|v|653βdx)53β6=cq6uεβv653βν4d20q26+Cuε2βv2653βν4q2+Cuε2βv2 (5.19)

    and

    cΩ|uε|β|v|2dxc(Ω|v|6dx)16(Ω|uε|2dx)β2(Ω|v|653βdx)53β6=cv6uεβv653β12v2+Cuε2βv2, (5.20)

    where the last inequalities in (5.19) and (5.20) are valid only if 653β6, i.e. β43, so combining (5.19), (5.20) with (5.18), we have

    c(|uε|β˜w,v)ν4q2+12v2+Cuε2βv2. (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|dxCΩ||uε|β1+|u0|β1||˜w||u0||v|dxCΩ|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

    |xpyp|Cp(xp1+yp1)|xy|

    for any x,y0, 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|dxCΩ|uε|β1|q||u0||v|dx+CΩ|uε|β1|v|2|u0|dx, (5.23)

    and

    CΩ|uε|β1|q||u0||v|dxC(Ω|q|6dx)16(Ω|uε|2dx)β12(Ω|u0|2dx)12(Ω|v|653β)53β6=Cq6uεβ1u0∥∥v653βν8q2+Cuε2(β1)u02v2653βν8q2+Cuε2(β1)u02v2, (5.24)

    similarly,

    CΩ|uε|β1|v|2|u0|dx14v2+Cuε2(β1)u02v2. (5.25)

    It follows from (5.23)-(5.25) that

    CΩ|uε|β1|˜w||u0||v|dxν8q2+14v2+Cuε2(β1)u02v2. (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ν8q2+14v2+Cu02(β1)u02v2. (5.27)

    Combining (5.26), (5.27) with (5.22), we have

    c((|uε|β|u0|β)u0,v)ν4q2+12v2+C(uε2(β1)+u02(β1))u02v2. (5.28)

    So it follows from (5.17), (5.21) and (5.28) that

    c(|uε|βuε|u0|βu0,v)ν2q2+v2+Cuε2βv2   +C(uε2(β1)+u02(β1))u02v2. (5.29)

    Now, considering (5.6), (5.16) and (5.29), and according to (5.1) and (5.2), it yields

    ddtq2+2aq24v2+C(R21+R20)v2+CR2β1v2   +C(R2(β1)1+R2(β1)0)R20v2C(1+R20+R21+R2β1+(R2(β1)0+R2(β1)1)R20)v2. (5.30)

    Recalling that q(τ)∥=0 and (4.16), Lemma 2.3 entails

    q(t)2C(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 ˆf1H(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)2Cl2ε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ε),ˆf0H(f0),ˆf1H(f1)

    such that ˆuε(0)=uε. For every L0 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 ˆf0H(f0). Then, for every δ>0, there is T=T(δ)0, independent of L, such that

    distH(Uˆf0(TL,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.
  • This article has been cited by:

    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
  • Reader Comments
  • © 2020 the Author(s), licensee AIMS Press. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0)
通讯作者: 陈斌, bchen63@163.com
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

Metrics

Article views(3906) PDF downloads(373) Cited by(6)

Other Articles By Authors

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog