Processing math: 100%
Research article Topical Sections

Dynamic analysis and comparison of the performance of linear and nonlinear controllers applied to a nonlinear non-interactive and interactive process


  • This article presents an in-depth dynamic analysis and comparative evaluation of three distinct control strategies—proportional-integral (PI) compensator, linear quadratic regulator (LQR), and sliding mode control (SMC)—applied to a nonlinear process in two configurations: non-interactive system (NIS) and interactive system (IS). The primary objective was to optimize the regulation of fluid levels in a dual-tank system subject to external disturbances and varying operational conditions. The process dynamics were initially modeled using nonlinear differential equations, which were subsequently linearized to facilitate the design of the PI and LQR controllers. The PI compensator design was rooted in state-space representation and was tuned using the Ziegler-Nichols method to achieve the desired transient and steady-state performance. The LQR design employed optimal control theory, minimizing a quadratic cost function to derive the state feedback gain matrix, ensuring system stability by shifting the eigenvalues of the closed-loop system matrix into the left half of the complex plane. In contrast, the SMC leveraged the full nonlinear dynamics of the process, establishing a sliding surface to drive the system states toward a desired trajectory with robustness against model uncertainties and external disturbances. The SMC's performance was evaluated by analyzing the existence and stability of the sliding mode using the derived switching laws for the actuation signal. The comparative study was conducted through simulations in MATLAB/Simulink environments, where each controller's performance was assessed based on transient response, robustness to disturbances, and computational complexity. The results indicate that while the PI compensator and LQR provide satisfactory performance under linearized assumptions, the SMC demonstrates superior robustness and precision in managing the nonlinearities inherent in the IS configuration. This comprehensive analysis underscores the critical trade-offs between simplicity, computational overhead, and control efficacy when selecting appropriate control strategies for nonlinear, multi-variable processes.

    Citation: José M. Campos-Salazar, Pablo Lecaros, Rodrigo Sandoval-García. Dynamic analysis and comparison of the performance of linear and nonlinear controllers applied to a nonlinear non-interactive and interactive process[J]. AIMS Electronics and Electrical Engineering, 2024, 8(4): 441-465. doi: 10.3934/electreng.2024021

    Related Papers:

    [1] Wei Shi, Xinguang Yang, Xingjie Yan . Determination of the 3D Navier-Stokes equations with damping. Electronic Research Archive, 2022, 30(10): 3872-3886. doi: 10.3934/era.2022197
    [2] Jie Qi, Weike Wang . Global solutions to the Cauchy problem of BNSP equations in some classes of large data. Electronic Research Archive, 2024, 32(9): 5496-5541. doi: 10.3934/era.2024255
    [3] Guoliang Ju, Can Chen, Rongliang Chen, Jingzhi Li, Kaitai Li, Shaohui Zhang . Numerical simulation for 3D flow in flow channel of aeroengine turbine fan based on dimension splitting method. Electronic Research Archive, 2020, 28(2): 837-851. doi: 10.3934/era.2020043
    [4] Jie Zhang, Gaoli Huang, Fan Wu . Energy equality in the isentropic compressible Navier-Stokes-Maxwell equations. Electronic Research Archive, 2023, 31(10): 6412-6424. doi: 10.3934/era.2023324
    [5] Keqin Su, Rong Yang . Pullback dynamics and robustness for the 3D Navier-Stokes-Voigt equations with memory. Electronic Research Archive, 2023, 31(2): 928-946. doi: 10.3934/era.2023046
    [6] Linlin Tan, Bianru Cheng . Global well-posedness of 2D incompressible Navier–Stokes–Darcy flow in a type of generalized time-dependent porosity media. Electronic Research Archive, 2024, 32(10): 5649-5681. doi: 10.3934/era.2024262
    [7] Xiao Su, Hongwei Zhang . On the global existence and blow-up for the double dispersion equation with exponential term. Electronic Research Archive, 2023, 31(1): 467-491. doi: 10.3934/era.2023023
    [8] Guochun Wu, Han Wang, Yinghui Zhang . Optimal time-decay rates of the compressible Navier–Stokes–Poisson system in R3. Electronic Research Archive, 2021, 29(6): 3889-3908. doi: 10.3934/era.2021067
    [9] Xu Zhao, Wenshu Zhou . Vanishing diffusion limit and boundary layers for a nonlinear hyperbolic system with damping and diffusion. Electronic Research Archive, 2023, 31(10): 6505-6524. doi: 10.3934/era.2023329
    [10] Jianxia He, Qingyan Li . On the global well-posedness and exponential stability of 3D heat conducting incompressible Navier-Stokes equations with temperature-dependent coefficients and vacuum. Electronic Research Archive, 2024, 32(9): 5451-5477. doi: 10.3934/era.2024253
  • This article presents an in-depth dynamic analysis and comparative evaluation of three distinct control strategies—proportional-integral (PI) compensator, linear quadratic regulator (LQR), and sliding mode control (SMC)—applied to a nonlinear process in two configurations: non-interactive system (NIS) and interactive system (IS). The primary objective was to optimize the regulation of fluid levels in a dual-tank system subject to external disturbances and varying operational conditions. The process dynamics were initially modeled using nonlinear differential equations, which were subsequently linearized to facilitate the design of the PI and LQR controllers. The PI compensator design was rooted in state-space representation and was tuned using the Ziegler-Nichols method to achieve the desired transient and steady-state performance. The LQR design employed optimal control theory, minimizing a quadratic cost function to derive the state feedback gain matrix, ensuring system stability by shifting the eigenvalues of the closed-loop system matrix into the left half of the complex plane. In contrast, the SMC leveraged the full nonlinear dynamics of the process, establishing a sliding surface to drive the system states toward a desired trajectory with robustness against model uncertainties and external disturbances. The SMC's performance was evaluated by analyzing the existence and stability of the sliding mode using the derived switching laws for the actuation signal. The comparative study was conducted through simulations in MATLAB/Simulink environments, where each controller's performance was assessed based on transient response, robustness to disturbances, and computational complexity. The results indicate that while the PI compensator and LQR provide satisfactory performance under linearized assumptions, the SMC demonstrates superior robustness and precision in managing the nonlinearities inherent in the IS configuration. This comprehensive analysis underscores the critical trade-offs between simplicity, computational overhead, and control efficacy when selecting appropriate control strategies for nonlinear, multi-variable processes.



    It is well-known that the Navier-Stokes equations are the typical evolution equations and widely used in the field of science and engineering. The attractors of Navier-Stokes equations are studied by many scholars in the fields of dynamical systems for a long time (see [1,2,3,4,5,6,7,8,9,10,11,12,13,14,15] and reference therein). Especially in recent years, there are many research achievements on g-Navier-Stokes equation. In [16,17,18], Roh deduced the 2D g-Navier-Stokes equations from 3D Navier-Stokes equations on thin region. It can be viewed as a perturbation of the usual Navier-Stokes equations. Bae et al. studied the well-posedness of weak solution for the 2D g-Navier-Stokes equations. Kwak et al. researched the global attractor and its fractal dimension of 2D g-Navier-Stokes equations in [19]. In [20,21,22,23,24], Jiang et al. studied global and the pullback attractor for g-Navier-Stokes equation. Moreover, the long-time behavior for 2D non-autonomous g-Navier-Stokes equations and the stability of solutions to stochastic 2D g-Navier-Stokes equations were studied by Anh in [25,26], The stationary solutions and its pullback attractor are researched in [27]. On the basis of the above research, we have studied the long time properties for g-Navier-Stokes equation with weakly dampness and time delay in [28] recently.

    In this manuscript, the uniform attractor of the g-Navier-Stokes equations with nonlinear dampness is researched. Its usual form is as follows:

    utνΔu+(u)u+c|u|β1u+p=f(x,t)in[τ,+)×Ω(gu)=0in[τ,+)×Ωu(x,t)=0in[τ,+)×Ωu(τ,x)=uτ(x)xΩ (1.1)

    In (1.1), we can see that u(t,x)R2 and p(t,x)R denote the velocity and pressure respectively. ν>0 is the viscosity coefficient, c|u|β1u denotes nonlinear dampness. c>0 and β1 are positive constant. f=f(x,t) is the external force term, 0<m0g=g(x1,x2)M0 and g=g(x1,x2) is a suitable smooth function, Let c=0 and g=1, the Eq (1.1) will become the usual 2D Navier-Stokes equations.

    This manuscript is organized as follows. In Section 2, we recall some basic results of 2D g-Navier-Stokes equations, then we give the concept about process families and uniform attractor. In Section 3, the global well-posedness of weak solutions for 2D g-Navier-Stokes equations with nonlinear dampness is studied. In Section 4, by the energy equation method, the existence of the uniform attractor of 2D g-Navier-Stokes equation with nonlinear dampness is proved on the unbounded domain. In Section 5, the dimension estimation of the uniform attractor in the quasi-periodic case is obtained.

    We assume Ω is a smooth unbounded domian of R2, Let L2(g)=(L2(Ω))2 and we denote (u,v)=Ωuνgdx and ||=(,)1/2,u,vL2(g). Let H10(g)=(H10(Ω))2, Set

    ((u,v))=Ω2j=1ujvjgdx,

    and ||||=((,))1/2,u=(u1,u2),v=(v1,v2)H10(g). We denote D(Ω) be the space of C functions with compact support contained in Ω. So we have the following spaces

    H={v(D(Ω))2:gv=0inΩ};
    Hg=closureofHinL2(g);
    Vg=closureofHinH10(g).

    where Hg and Vg endowed with the inner product and norm of L2(g) and H10(g) respectively.

    We assume that there exists λ1>0, such that

    |u|21λ1||u||2,uVg. (2.1)

    This Poincarˊe-type inequality imposes some restrictions on the geometry of the domain Ω.

    The g-Laplacian operator is defined as follows:

    Δgu=1g(g)u=Δu1ggu.

    The first equation of (1.1) can be rewritten as follows:

    utνΔgu+νggu+c|u|β1u+(u,)u+p=f. (2.2)

    In [16], g-orthogonal projection and g-Stokes operator are defined respectively by Pg:L2(g)Hg and Agu=Pg(1g((gu))). Applying the projection Pg on the Eq (2.2), we have the following weak formulation of (1.1).

    ddt(u,v)+ν((u,v))+c(|u|β1u,v)+bg(u,u,v)+ν(Ru,v)=f,vvVg,t>0, (2.3)
    u(0)=u0, (2.4)

    where bg:Vg×Vg×VgR and

    bg(u,v,w)=2i,j=1uivjxwjgdx,

    we have

    Ru=Pg[1g(g)u],uVg.

    Then the formula (2.3) and (2.4) are equivalent to the following functional equations

    dudt+νAgu+c|u|β1u+Bu+νRu=f (2.5)
    u(0)=u0 (2.6)

    We denote

    Agu,v=((u,v)),u,vVg. (2.7)

    From [16,17,19], we have

    ||B(u)||Vgc|u|||u||,||Ru||Vg|g|m0λ1/21||u||,uVg.

    where B(u)=B(u,u)=Pg(u)u is defined by

    B(u,v),w=bg(u,v,w),u,v,wVg.

    A family of two parametric maps {Uf(t,τ)}={Uf(t,τ)|tτ,τR} is defined in Hg as follows:

    Uf(t,τ):EE,tτ,τR.

    The following concepts and conclusions are given from [7]. fL(R+;Vg), the translation operator is defined in L(R+;Vg) as follows.

    T(h)f(s)=f(s+h),h0,sR.

    Obviously

    ||T(h)f||L(R+;Vg)||f||L(R+;Vg),h0,fL(R+;Vg).

    We set Σ={T(h)f(x,s)=f(x,s+h),hR}, where T() is the positive invariant semigroups which act on Σ and satify T(h)ΣΣ,h0 and

    UT(h)f(t,τ)=Uf(t+h,τ+h),h0,tτ0.

    Let ρF>0 be constant, Σ{fL(R+;Vg):||f||L(R+;Vg)ρF}. For {Uf(t,τ)} with fΣ, we call the parameter f as the symbols of the process family {Uf(t,τ)}, and Σ as the symbol space.

    Definition 2.1 [7] A family of two-parametric maps {U(t,τ)} is called a process in Hg, if

    (1)Uf(t,s)Uf(s,τ)=Uf(t,τ),tsτ,τR,

    (2)Uf(τ,τ)=Id,τR.

    Let E be the Banach space, B(E) is denoted the set of all bounded sets on E, then

    Definition 2.2 [7] A set B0E is said to be uniformly absorbing for the family of processes {Uf(t,τ)},fΣ}, if for any τR and each BB(E), there exists t0=t0(τ,B)τ, such that for all tt0,

    fΣUf(t,τ)BB0.

    Definition 2.3 [7] A set PE is said uniformly atttracting set of {Uf(t,τ)},fΣ}, if for any τR, there is

    limt+(supfΣdistE(Uf(t,τ)B,P))=0.

    Definition 2.4 [7] A closed set AΣE is said to be the uniform attractor of the family of processes {Uf(t,τ)},fΣ}, if

    (1) AΣE is uniformly attractive;

    (2) AΣE is included in any uniformly attracting set A of {Uf(t,τ)},fΣ}, that is AΣA.

    In the section we will prove the well-posedness of the solution for 2D g-Navier-Stokes equations with nonlinear dampness by the Faedo-Galerkin method.

    Definition 3.1 Let u0Hg,fL2Loc(R;Vg), For any τR, uL(τ,T;Vg)L2(τ,T;Vg)Lβ+1(τ,T;Lβ+1(Ω)),T>τ is called a weak solution of problem (1.1) if it fulfils

    ddtu(t)+νAgu(t)+B(u(t))+c|u|β1u+νR(u(t))=f(x,t)onD(τ,+;Vg),u(τ)=u0.

    Theorem 3.1 Let β1, fL2Loc(R;Vg), Then for every uτVg, the equations (1.1) have a unique weak solution u(t)=u(t;τ,uτ)L(τ,T;Vg)L2(τ,T;Vg)Lβ+1(τ,T;Lβ+1(Ω)), and u(t) is continuously depending on the initial value in Vg.

    Proof. Let {wj}j1 be the eigenfunctions of Δ on Ω with homogeneous Dirichlet boundary conditions, Its corresponding eigenvalues are 0<λ1λ2, Obviously, {wj}j1Vg forms a Hilbert basis in Hg, given uτVg and fL2Loc(R;Vg).

    For any positive integer n1, we structure the Galerkin approximate solutions as un(t)=un(t;T,uτ), It has the following form

    un(t,T;uτ)=nj=1γn,j(t)wj.

    where γn,j(t) is determined from the initial values of the following system of nonlinear ordinary differential equations.

    (un(t),wj)+ν((un(t),wj))+c(|un(t)|β1un(t),wj)+b(un(t),un(t),wj)+b(gg,un(t),wj)
    =f(x,t),wj,t>τ,j=1,2,n (3.1)
    ((un(t),wj))=((uτ,wj)).

    where is dual product of Vg and Vg.

    According to the results of the initial value problems of ordinary differential equations, we have that there exists a unique local solution of (3.1). In the following, we prove that the time interval of the solution can be extended to [τ,).

    12ddt|un(t)|22+ν||un(t)||2+c|un(t)|β+1β+1+b((gg)un(t),un(t))=f(x,t),un(t) (3.2)

    Using Cauchy's inequality and Young's inequality, we have

    f(x,t),un(t)||f(x,t)||||un(t)||
    ν2||un||2+12ν||f(x,t)||2 (3.3)

    where |||| is norm of Vg. We take (3.3) into (3.2) to obtain

    12ddt|un(t)|22+ν||un(t)||2+c|un(t)|β+1β+1+b((gg)un(t),un(t))
    ν2||un||2+12ν||f(x,t)||2
    ddt|un(t)|22+2ν||un(t)||2+2c|un(t)|β+1β+1+2b((gg)un(t),un(t))
    ν||un||2+1ν||f(x,t)||2
    ddt|un(t)|22+ν||un(t)||2+2c|un(t)|β+1β+1+2b((gg)un(t),un(t))1ν||f(x,t)||2 (3.4)

    That is

    ddt|un(t)|22+ν||un(t)||2+2c|un(t)|β+1β+11ν||f(x,t)||2+2ν|g|m0λ1/21||un(t)||2
    ddt|un(t)|22+ν(12|g|m0λ1/21)||un(t)||2+2c|un(t)|β+1β+11ν||f(x,t)||2 (3.5)

    By integrating (3.5) from τ to t, we have

    |un(t)|2+ν(12|g|m0λ1/21)tτ||un(s)||2ds+2ctτ|un(s)|β+1β+1ds
    |un(τ)|2+1νtτ||f(x,s)||2ds.

    For any T>0 and β1, we obtain

    supτtT(|un(t)|2)+ν(12|g|m0λ1/21)tτ||un(s)||2ds+2ctτ|un(s)|β+1β+1ds
    |un(τ)||2+1νtτ||f(x,s)||2dsC.
    Sowecanobtainthat{un(t)}isboundedinL(τ,T;Vg), (3.6)
    {un(t)}isboundedinL2(τ,T;Vg), (3.7)
    and{un(t)}isboundedinLβ+1(τ,T;Lβ+1(Ω)). (3.8)

    So un(t)L(τ,T;Vg). Therefore B(un(t))L(τ,T;Vg), |un(t)|β1un(t)Lβ+1(τ,T;Lβ+1(Ω)). As a result,

    ddtun(t),v=f(x,t)c|un(t)|β1un(t)νAun(t)B(un(t))νR(un(t)),v,vVg.

    so {un(t)} is bounded in L2(τ,T;Vg).

    Then we deduce that there is a subsequence in {un(t)}, which is still denoted by {un(t)}. We obtain un(t)L2(τ,T;Vg) and un(t)L2(τ,T;Vg) such that

    (i)un(t)u(t) is weakly convergent in L(τ,T;Vg);

    (ii)un(t)u(t) is weakly convergent in L2(τ,T;Vg);

    (iii)|un(t)|β1un(t)ξ is weakly convergent in Lβ+1(τ,T;Lβ+1(Ω));

    (iv)un(t)u(t) is weakly convergent in L2(τ,T;Vg);

    (v)un(t)u(t) is strongly convergent in L2(τ,T;Hg);

    (vi)un(t)u(t),ae(x,t)Ω×[τ,T].

    From Lemma 1.3 of [29], we can see ξ=|u|β1u. Since nN+span{w1,w2,,wn} is denseness in Vg, Taking the limit n on both sides of Eq (3.1), we can obtain that u is a weak solution of (1.1).

    In the following, the solution is proved to be unique and continuously dependent on initial values. Let u1,u2 be two weak solutions of (1.1) corresponding to the initial values u1τ,u2τVg, We take u=u1u2, From (2.3) we have

    12ddt(|u|2)+ν||u||2+c(|u1|β1u1|u2|β1u2,u)+ν(Ru,u)=B(u2)B(u1),u. (3.9)

    Using Hölder inequality and Sobolev embedding theorem, we obtain

    (|u1|β1u1|u2|β1u2,u)=Ω(|u1|β1u1|u2|β1u2)(u1u2)dx
    Ω(|u1|β+1|u1|β|u2||u2|βu1+|u2|β+1)dx
    =Ω(|u1|β|u2|β)(|u1||u2|)dx0. (3.10)

    we have

    |B(u2)B(u1),u=|B(u2,u2u1)B(u1u2,u1),u|
    C1||u2||||u2u1||||u||+C1||u1u2||||u1||||u||
    =C1||u||2(||u1+||u2||)
    C1||u||2 (3.11)

    where C1>0 is any constant.

    ν(Ru,u)ν||g||m0λ1/21||u|||u|
    ν||g||2m0λ1/21(||u||2+|u|2)
    =α(||u||2+|u|2). (3.12)

    where α=ν||g||2m0λ1/21. so

    12ddt|u|2+ν||u||2C1||u||2+α(||u||2+|u|2).
    ddt|u|2+2(νC1α)λ1|u|2α|u|2.

    Thus

    ddt|u|2[α2λ1(νC1α)]|u|2.

    Let C be a constant and C=α2λ1(νC1α)>0, then

    ddt|u|2C|u|2.

    Therefore

    |u|22eC(tτ)|uτ|22.

    So we prove the continuous dependence on the initial value. When u1τ=u2τ, that is uτ=0, then the uniqueness of the solution holds.

    In the following we have that the family of processes {Uf(t,τ)},fΣ is uniformly bounded (w.r.t.fΣ) and it has uniform absorbing sets.

    Firstly, the existence of uniformly absorbing sets is proved. Taking the inner product of (2.5) with u, we have

    ddt|u|2+2ν||u||2+2c|u|β+1=2f,u2ν((gg)u,u),

    Then

    ddt|u|2+2ν||u||2+2c|u|β+1||f||2L(R+;Vg)ν+ν||u||2+2ν|g|m0λ1/21||u||2,

    For β1, we obtain

    ddt|u|2+νλ1γ|u|2ddt|u|2+νγ||u||2||f||2L(R+;Vg)ν,

    where γ=12|g|m0λ1/21>0 for sufficiently small |g|. Using the Gronwall inequality, we have

    |u(t)|2|u0|2eνλ1γt+||f||2L(R+;Vg)ν2λ1γ,t>0.

    and from

    ddt|u|2+2ν||u||2+2c|u|β+1||f||2L(R+;Vg)ν+ν||u||2+2ν|g|m0λ1/21||u||2,

    we have

    ddt|u|2+ν||u||2+2c|u|β+1||f||2L(R+;Vg)ν+2ν|g|m0λ1/21||u||2.

    So

    ddt|u|2+ν(12|g|m0λ1/21)||u||2||f||2L(R+;Vg)ν. (4.1)

    Integrating (4.1) in s from 0 to t, we have

    1tt0||u(s)||2ds|u0|2tνγ+||f||2L(R+;Vg)ν2γ,t>0.

    then we know that the family of processes corresponding to u is uniformly bounded, and

    B0={uHg:|u|ρ0=1ν2λ1γ||f||L(R+;Vg)}

    is uniformly absorbing set in Hg. Then the following lemma holds.

    Lemma 4.1 Let Σ be symbolic space, The process family corresponding to Eq (1.1) is uniformly bounded in L(R+;Hg)L2(τ,T;Vg) and there is a uniform absorbing set in Hg.

    Lemma 4.2 Let τ0, uτn be the sequence in Hg that weakly converges to uτHg, fnΣ is the sequence in L(R+;Vg) that weakly converges to f, then

    (1) For t>τ, Ufn(t,τ)uτn is weakly converges to Uf(t,τ)uτ in Hg;

    (2) For T>τ, Ufn(,τ)uτn is weakly converges to Uf(,τ)uτ in L2(τ,T;Vg).

    The proof is similar to Lemma 3.2 of [7], so it is omitted.

    As we know, when uτn is bounded in Hg, fnΣ, tn+. If {Ufn(tn,τ)uτn} is precompact in Hg, then the family of processes {Uf(t,τ)},fΣ is asymptotically compact. So we construct an energy functional [,]:Vg×VgR as follows:

    [u,v]=ν((u,v))+ν2((gg,)u,v)+ν2((gg,)v,u)νλ14(u,v)+c(|u|β1u,v),u,vVg.

    Obviously [,] is bilinear and symmetric, and

    [u]2=[u,u]=ν||u||2+ν((gg)u,u)νλ14|u|2+c|u|β+1
    ν||u||2ν(|g|m0λ1/21+14)||u||2
    ν2||u||2. (4.2)

    Let |g| be sufficiently small in (4.2), such that |g|m0λ1/21<14. Hence

    ν2||u||2[u]232ν||u||2,uVg.

    Since

    ddt|u|2+vλ12|u|2+2[u]2=2(f,u),

    Given u=u(t)=Uf(t,τ)uτ,uτHg,tτ0, Then we have

    |Uf(t,τ)uτ|2=|uτ|2eνλ1(tτ)/2+2tτeνλ1(ts)/2((f,Uf(s,τ)uτ)[Uf(s,τ)uτ]2)ds.

    That is uτHg,tτ0, we obtain

    |Uf(t,τ)uτ|2=|uτ|2eνλ1(tτ)/2
    +2tτ0eνλ1(tτs)/2((T(τ)f(s),UT(τ)f(s,0)uτ)[UT(τ)f(s,0)uτ]2)ds.

    Lemma 4.3 Let {Uf(t,τ)}fΣ is the family of processes of Eq (1.1), then {Uf(t,τ)}fΣ is uniformly asymptotically compact.

    Proof. Let BHg is bounded, uτnB,fnΣ and tnR+ is satisfied tn+(n+). From Lemma 4.1, we have a constant M(B,τ)>τ and

    Uf(t,τ)BB0,tM(B,τ),fΣ.

    There exists sufficiently large tnM(B,τ), such that Ufn(tn,τ)BB0. then {Ufn(tn,τ)uτn} is weakly precompact in Hg. For wB0Hg, we can deduce that Ufns(tns,τ)uτns is weakly convergent to w in Hg. Similarly T>0 and tnsT+M(B,τ), we obtain Ufns(tnsT,τ)uτnsB0. The same to wTB0, we can take ns, T>0, so we have utns=Ufns(tnsT,τ)uτns is weakly convergent to wT in Hg. According to the definition of process and translation operator, we have

    Ufns(tns,τ)=UT(tnsT)fns(T,0)Ufns(tnsT,τ).

    Let gT,ns=T(tnsT)fns, we denote limnsHw as weak limit in Hg, then

    w=limnsHwUfns(tns,τ)uτns=limnsHwUgT,ns(T,0)utns=UgT(T,0)wT,

    thus

    |w|lim infns|Ufns(tns,τ)uτns|=lim infns|UgT,ns(T,0)utns|.

    Now we will prove

    lim supns|Ufns(tns,τ)uτns||w|.

    T>0, we have wk=UgkT(T,0)wT. When tnsT+M(B,τ), we obtain

    UgkT,ns(T,0)utns
    =2t0eνλ1(Ts)/2((gkT,ns(s),UgkT,ns(s,0)utsn)[UgkT,ns(s,0)utns]2)ds+|utns|2eνλ1T/2

    Obviously

    lim supns(eνλ1T/2|utns|2)ρ20eνλ1T/2.

    From Lemma 4.2, we obtain

    T0eνλ1T/2[UgkT(s,0)wT]2dslim infnsT0eνλ1T/2[UgkT,ns(s,0)utns]2ds.

    So

    lim supns2T0eνλ1T/2[UgkT,ns(s,0)utns]2ds
    =2lim infnsT0eνλ1T/2[UgkT,ns(s,0)utns]2ds
    2T0eνλ1T/2[UgkT(s,0)wT]2ds.

    For

    limnsT0eνλ1(Ts)/2(gkT,ns(s),UgkT,ns(s,0)utns)ds=T0eνλ1(Ts)/2(gkT(s),UgkT(s,0)wT)ds

    thus

    lim supns|UgkT,ns(T,0)utns|22T0eνλ1(Ts)/2((gkT(s),UgsT(s,0)wT)[UgsT(s,0)wT]2)ds
    +ρ20eνλ1T/2.

    From wk=UgsT(s,0)wT, we have

    |wk|2=|UgsT(s,0)wT|2
    =eνλ1T/2|wT|2+2T0eνλ1(Ts)/2((gkt(s),UgkT(s,0)wT)[UgkT(s,0)wT]2)ds.

    T>0, we have

    lim supns|UgkT,ns(T,0)utns|2|wk|2+(ρ20|wT|2)eνλ1T/2|wk|2+ρ20eνλ1T/2.

    From w=UgT(T,0)wT, by the Lemma 3.3 of [7] and Lemma 4.2, we can obtain wkw in Hg. So there exists any sufficiently small ε>0, such that |wk|2|w|2+ε. Since

    supns+|Ufns(tns,τ)uτns|2=supns+|UgT,ns(T,0)utns|2|w|2+ε+ρ20eνλ1T/2.

    When ε0, T, we have

    lim supns|Ufns(tns,τ)uτns|2|w|2.

    Let BHg be any bounded set, we have

    ωτ,Σ(B)=tτ¯fΣstUf(s,τ)B.

    and vωτ,Σ(B) iff there exists a sequence vnB,fnΣ,tn[τ,+). When n, we have tn+ and Ufn(tn,τ)vnv in Hg. When {Uf(t,τ)},fΣ} is uniformly asymptotically compact, t+, we have

    supfΣdistHg(Uf(t,τ)B,ωτ,Σ(B))0.

    We will obtain the minimization of the uniform attractor in the following.

    Lemma 4.4 Let {Uf(t,τ)},fΣ is any the family of processes, B0 is uniformly absorbing set, AΣ=ω0,Σ(B0). then AΣ is contained in any uniform absorbing set of {Uf(t,τ)},fΣ.

    Proof. τ>0,BHg, Suppose there is another bounded closed set PHg which satisfies

    limtsupfΣdistHg(Uf(t,τ)B,P)=0,

    where AΣ is not contained in the P. We deduce there is at least one vAΣ and vP. Since vAΣ=ω0,F(B0), From the definition of the uniform ω limit set, there is a sequence vnB,fnΣ,tn[τ,+), as n, we have tn+, then Ufn(tn,0)vnv is obtained in Hg. Given ˜vn=Ufn(tn,0)vn, when n+, There must be Ufn(tn,τ)˜vnv. Let ˜vnB, then we obtain vP, It is contradiction, so AΣP.

    Theorem 4.1 Let {Uf(t,τ),fΣ} is a family of processes of Eq (1.1), Then the process family has a unique compact uniform attractor AΣHg. where AΣ=ω0,Σ(B0), B0 is any uniform absorbing set corresponding to a family of processes.

    When f(x,t)=f(x,w1(t),w2(t),,wk(t)) is a quasi-periodic function, That is, there exists a set of rational independent real numbers α1,,αk which satisfies f(x,α1t,,αit+2π,,αkt)=f(x,α1t,,αit,,αkt)(1ik). Here w1(t+α1)=w1(t),w2(t+α2)=w2(t),,wk(t+αk)=wk(t) and α1,α2,,αk are rational independent.

    Let αt=(α1t,,αkt), α=(α1,,αk), w(t)=(w1(t),,wk(t))=[αt+w0]=(αt+w0)mod(2π)k, w0=(w01,,w0k)Tk=[0,2π]k,F(x,w(t))=f(x,t). we can obtain the following conclusion.

    Theorem 5.1 Let A is the uniform attractor of (1.1), then its Hausdorff and Fractal dimensions are estimated as follows:

    dH(A)4νm1λ1(C2q|f|22ν3λ1m1m2+G2kνm1λ1)+k+1.
    dF(A)16νm1λ1(C2q|f|22ν3λ1m1m2+G2kνm1λ1)+2k+2.

    where

    G=(ki=1|Fwi|2BC(Tk,Hg))12,
    m1=1+|g|m0λ1/21,m2=1|g|m0λ1/21.

    Proof. We transform the Eq (1.1) into the following forms of autonomous systems by semigroup S(t)(u0,w0)=(Uw0(t,0)u0,T1(t)w0),

    utνΔgu+B(u,u)+ν(gg)u+c|u|β1u+p=F(x,w(t)) (5.1)
    ddtw(t)=α (5.2)
    u|t=0=u0,u(t)|t=0=w0,u0Hg,w0Tk. (5.3)

    Let y(t)=(u(x,t),w(t))T, M(y(t))=(νΔguB(u,u)ν(gg)uc|u|β1up+F(x,w(t)),α)T.

    Then we can write the Eqs (5.1) and (5.2) as follows,

    y(t)t=M(y(t)) (5.4)
    y(t)|t=0=y0=(u0,w0) (5.5)

    y0A, where y(t)=(u(t),w(t))T is the solution of Eqs (5.1) and (5.2) and y0 as initial value. The linearized equation of (5.1) in y(t) is

    z(t)t=M(y(t))z (5.6)
    z(t)|t=0=z0 (5.7)

    In the equation of (5.6), z(t)=(v(t),w(t))T,μ(t)=(μ1(t),,μk(t)),z0=(v0,μ0)THg×Tk,

    M(y(t))z=(νΔgvB(v,u)B(u,v)ν(gg)vc|u|β1v˜p+Fwu,0)T.

    while

    (M(y(t))z,z)=ν||v||2b(v,u,v)ν|g|m0λ1/21||v||2c|u|β1|v|2ΩFwuvdx
    ν||v||2+Ω|u||v|2dx+ν|g|m0λ1/21||v||2+c|u|β1λ1||v||2+Ω|Fwu||v|dx

    Let

    G=(ki=1|Fwi|2BC(Tk,H))1/2,

    then

    (M(y(t))z,z)ν(1|g|m0λ1/21c|u|β1λ1)||v||2+Ω|u||v|2dx+bG2|v|2+G2b|u|2.

    b is any positive constant, Let (M(y(t))z,z)=(M1v,v)+(M2μ,μ),

    M1v=ν(1|g|m0λ1/21c|u|β1λ1)Δgv+(|u|+bG2)v,M2μ=G2bIkμ.

    Ik is identity operator in Rk. Thus operator

    ~M1=(M100M2)

    is block operator, (v1,0),,(vnk,0) are respectively solution of (2.7) with (ξ1,0),,(ξnk,0) as the initial value, where ξ1,,ξnk is linearly independent basis in Hg, ϕ1,,ϕnk is unit orthogonal basis of span{v1,,vnk}, ˜μnk+1,,˜μn is unit orthogonal basis of Rk, then (ϕ1,0),,(ϕnk,0), (0,˜μnk+1),,(0,˜μn) is unit orthogonal basis of Hg×Rk. Let θi=(ϕi,0)(i=1,,nk),vi=(0,~μi),(nk+1in).

    qn=limTinfsupy0A(1TT0n1(M(y(s))θi,θi)ds).
    n1(M(y(s))θi,θi)ν(1|g|m0λ1/21c|u|β1λ1)nki=1||φi||2+||u||(Ω(nki=1|ϕi|2)2dx)1/2+bG2(nk)+G2bk.

    We make m1=1|g|m0λ1/21c|u|β1λ1, then

    n1(M(y(s))θi,θi)νm1nki=1||φi||2+||u||(Ω(nki=1|ϕi|2)2dx)1/2+bG2(nk)+G2bk
    νm12nki=1||φi||2+C2q2νm1||u||2+bG2(nk)+G2bk.

    We take b=νm1λ12G, then

    n1(M(y(s))θi,θi)νm1λ14(nk)+C2q2νm1||u||2+G2kνm1λ1.

    Since

    12ddt|u|2+ν||u||2=(f,u)ν((gg)u,u)ν(|u|β1u,u),
    ddt|u|2+2ν||u||2|f|2νλ1+νλ1|u|2+2ν|g|m0λ1/21||u||2+ν|u|β1λ1||u||2.

    That is

    ddt|u|2+ν(12|g|m0λ1/21|u|β1λ1)||u||2|f|2νλ1.

    We take m2=12|g|m0λ1/21|u|β1λ1, then

    ddt|u|2+νm2||u||2|f|2νλ1.

    So ||u||2|f|2ν2λ1m2. therefore

    qnνm1λ14(nk)+C2q|f|22ν3λ1m1m2+G2kνm1λ1.

    Let

    n0=[4νm1λ1(C2q|f|22ν3λ1m1m2+G2kνm1λ1)]+k+1,

    where [] denotes trunc, then we have qn0<0. so

    dH(A)4νm1λ1(C2q|f|22ν3λ1m1m2+G2kνm1λ1)+k+1.

    We take another

    n1=[8νm1λ1(C2q|f|22ν3λ1m1m2+G2kνm1λ1)]+k+1,

    then qn1<0, and max1jn11(qj)+qn1<1, thus

    dF(A)16νm1λ1(C2q|f|22ν3λ1m1m2+G2kνm1λ1)+2k+2.

    In this paper, using a priori estimates of the solutions and the energy equation method, we show how to control the nonlinear dampness and obtain the uniform attractor of the g-Navier-Stokes equation on unbounded domain. Meanwhile, the dimension of the uniform attractor is estimated in the quasi-periodic case. The methods in this paper can bring some inspiration for the research of 3D Navier-Stokes equations in the future.

    From a theoretical point of view, it is important to analysis the connection between Navier-Stokes equations and g-Navier-Stokes equations. So it is of great significance to study the dynamics for the g-Navier-Stokes equations. To obtain more research results for the study of g-Navier-Stokes equations in the next research, we may continue the research in this line, extending the case of Lebesgue space L2 to the case of L2,λ, for suitable 0<λ<2. On the other hand, we may consider that the pullback asymptotic behavior of solutions for 2D g-Navier-Stokes equations with nonlinear dampness on the unbounded domain.

    The author would like to thank the referees for the helpful suggestions. This work is supported by the National Natural Science Foundation of China (NSFC) (Grant No. 11971378).

    The authors declare this work does not have any conflicts of interest.



    [1] Saju S, Revathi R, Suganya KP (2014) Modeling and Control of Liquid Level Non-Linear Interacting and Non-Interacting System. International Journal of Advanced Research in Electrical, Electronics and Instrumentation Engineering 3: 8003-8013.
    [2] Vijayachitra S, Pooja S (2017) Modeling and Simulation of Level Control Phenomena in a Non-Linear System. International Journal of Intellectual Advancements and Research in Engineering Computations 5: 465-469.
    [3] Naşcu I, De Keyser R, Naşcu I, Buzdugan T (2010) Modeling and Simulation of a Level Control System. Proceedings of the 2010 IEEE International Conference on Automation, Quality and Testing, Robotics (AQTR) 1: 1–6. https://doi.org/10.1109/AQTR.2010.5520894 doi: 10.1109/AQTR.2010.5520894
    [4] Mukherjee D, Kundu PK, Ghosh A (2016) PID Controller Design for an Interacting Tank Level Process with Time Delay Using MATLAB FOMCON Toolbox. 2016 2nd International Conference on Control, Instrumentation, Energy & Communication (CIEC), 1-5. https://doi.org/10.1109/CIEC.2016.7513803
    [5] Brahman K, Baruah L, Deori H, Brahma S (2015) Study of Interacting and Non-Interacting with Disturbance and PID Controller Design. International Journal of Advanced Computing and Electronics Technology (IJACET), 2.
    [6] Meghna PR; Saranya V, Pandian BJ (2017) Design of Linear-Quadratic-Regulator for a CSTR Process. IOP Conf Ser Mater Sci Eng 263: 052013. https://doi.org/10.1088/1757-899X/263/5/052013 doi: 10.1088/1757-899X/263/5/052013
    [7] Lengyel K, Dulf ÉH, Kovács L (2020) Linear Quadratic Control on a Cascaded Multitank System. Proceedings of the 2020 IEEE 24th International Conference on Intelligent Engineering Systems (INES), 185–190. https://doi.org/10.1109/INES49302.2020.9147194
    [8] Anbumani K, Ranihemamalini R (2019) Linear Quadratic Regulator for Three Interacting Cylindrical Tank Control. Int J Resent Technol Eng 8: 19.
    [9] Aktaş M, Altun Y, Erol O (2017) LQR Control of Liquid Level and Temperature Control for Coupled-Tank System. International Conference on Hydraulics and Pneumatics, 79.
    [10] Anbumani K, Hemamalini RR (2020) Optimal State Feedback Controller for Three Tank Cylindrical Interacting System Using Grey Wolf Algorithm. Microprocessors and Microsystems 79: 103269. https://doi.org/10.1016/j.micpro.2020.103269 doi: 10.1016/j.micpro.2020.103269
    [11] Parvat B, Ratnaparkhi SD (2015) A Second Order Sliding Mode Controller Applications in Industrial Process. International Journal of Engineering Trends and Technology 19: 217–222. https://doi.org/10.14445/22315381/IJETT-V19P238 doi: 10.14445/22315381/IJETT-V19P238
    [12] Toms T, Hepsiba D (2014) Comparison of PID Controller with a Sliding Mode Controller for a Coupled Tank System. International Journal of Engineering Research & Technology 3: 151-154.
    [13] Arunshankar J (2018) Control of Nonlinear Two-Tank Hybrid System Using Sliding Mode Controller with Fractional-Order PI-D Sliding Surface. Computers & Electrical Engineering 71: 953–965. https://doi.org/10.1016/j.compeleceng.2017.10.005 doi: 10.1016/j.compeleceng.2017.10.005
    [14] Parvat BJ, Jadhav VK, Lokhande NN (2012) Design and Implementation of Sliding Mode Controller for Level Control. IOSR Journal of Electronics and Communication Engineering (IOSR-JECE) 2: 51-54.
    [15] Efe MÖ (2007) MIMO Variable Structure Controller Design for a Bioreactor Benchmark Process. ISA Transactions 46: 459–469. https://doi.org/10.1016/j.isatra.2007.03.005 doi: 10.1016/j.isatra.2007.03.005
    [16] Biswas PP, Srivastava R, Ray S, Samanta AN (2009) Sliding Mode Control of Quadruple Tank Process. Mechatronics 19: 548–561. https://doi.org/10.1016/j.mechatronics.2009.01.001 doi: 10.1016/j.mechatronics.2009.01.001
    [17] Smith CA, Corripio AB (2005) Principles and Practice of Automatic Process Control, 2nd edition, John Wiley & Sons.
    [18] Almutairi NB, Zribi M (2006) Sliding Mode Control of Coupled Tanks. Mechatronics 16: 427–441. https://doi.org/10.1016/j.mechatronics.2006.03.001 doi: 10.1016/j.mechatronics.2006.03.001
    [19] Ogata K (2009) Modern Control Engineering, 5th edition, Pearson: Boston.
    [20] Brogan WL (1991) Modern Control Theory, Prentice Hall.
    [21] Houpis CH, Sheldon SN, D'Azzo JJ (2003) Linear Control System Analysis and Design: Fifth Edition, Revised and Expanded, CRC Press. https://doi.org/10.1201/9780203911426
    [22] Dorato P, Abdallah CT, Cerone V (2000) Linear Quadratic Control: An Introduction, Krieger Pub Co: Melbourne, FL.
    [23] Khalil H (2014) Nonlinear Control, 1st edition, Pearson: Boston.
  • Reader Comments
  • © 2024 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(1083) PDF downloads(131) Cited by(0)

Figures and Tables

Figures(9)  /  Tables(4)

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog