Next Article in Journal
Parametric Study on the Sensitivity of Different Properties of Waves Propagating through an Incompressible Medium
Next Article in Special Issue
The Gravity Force Generated by a Non-Rotating Level Ellipsoid of Revolution with Low Eccentricity as a Series of Spherical Harmonics
Previous Article in Journal
A Posteriori Error Estimators for the Quasi-Newtonian Stokes Problem with a General Boundary Condition
Previous Article in Special Issue
Nonlinear Stability of the Monotone Traveling Wave for the Isothermal Fluid Equations with Viscous and Capillary Terms
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

The Regulator Problem to the Convection–Diffusion Equation

by
Andrés A. Ramírez
and
Francisco Jurado
*
División de Estudios de Posgrado e Investigación, Tecnológico Nacional de México/I.T. La Laguna, Revolución Blvd. and Instituto Tecnológico de La Laguna Av., Torreón 27000, Mexico
*
Author to whom correspondence should be addressed.
Mathematics 2023, 11(8), 1944; https://doi.org/10.3390/math11081944
Submission received: 24 March 2023 / Revised: 12 April 2023 / Accepted: 17 April 2023 / Published: 20 April 2023

Abstract

:
In this paper, from linear operator, semigroup and Sturm–Liouville problem theories, an abstract system model for the convection–diffusion (C–D) equation is proposed. The state operator for this abstract system model is here defined as given in the form of the Sturm–Liouville differential operator (SLDO) plus an integral term of the same SLDO. Our aim is to achieve the trajectory tracking task in the presence of external disturbances to the C–D equation invoking the regulator problem theory, where the state from a finite-dimensional exosystem is the state to the feedback law. In this context, the regulator (Francis) equations, established from the abstract system model for the C–D equation, here are solved; i.e., the state feedback regulator problem (SFRP) for the C–D system has a solution. Our proposal is validated via numerical simulation results.

1. Introduction

Systems whose dynamics evolve in an infinite-dimensional Hilbert space are denominated as infinite-dimensional systems modeled by partial differential equations (PDEs), which are also termed as distributed parameter systems (DPSs), since it reflects the spatial distribution of a physical quantity. The main goal when designing the control system for these class of systems is satisfying stability in the presence of external disturbances.
Dynamical systems given either in the input–output equation form, modeled through ordinary differential equations (ODEs), or in the state-space form can be transformed to the transfer function (transfer matrix) form; this latter description is always rational with real coefficients. Transfer functions from DPSs are non-rational functions which can be analytic in the complex plane and having no poles, such as in the case of the transport equation, namely, a first-order PDE, or having only zeros in their denominator, such as for the diffusion equation with Neumann boundary conditions or for the wave equation [1]. From classical control theory, classical controllers are designed from the knowledge of the transfer function, i.e., from an output/input description of the system. From DPSs, if a closed-form expression of their transfer function is provided, then the direct design of the controller may be possible. This approach is referred to as direct controller design. The primary drawback from this approach is the requirement of an explicit representation of the transfer function. In addition, the controller design will be infinite-dimensional, so this must be approximated by a finite-dimensional system. For some practical applications, when a transfer function for a DPS is not available, then the indirect controller design approach is the most common alternative to be employed. It consists of obtaining a finite-dimensional approximation of the system from which the controller can be designed [2].
The design of a feedback law such that it guarantees the tracking of a reference signal in the presence of an external disturbance, the latter generated through an exosystem, is the main objective when invoking the regulator problem theory [3]. Beyond finite-dimensional systems, the regulator problem theory has been playing an interesting role in the control of infinite-dimensional systems. In this work, we deal with the state feedback regulator problem (SFRP) where the state of the feedback law is from a finite-dimensional exosystem.
From linear finite-dimensional systems, the regulator problem theory has been extended to infinite-dimensional systems also known as DPSs [4,5,6,7,8]. In [6,7], control systems governed by a discrete spectral operator were introduced, where the so-called state operator meets with the property of spectrum decomposition [9,10] from which a controllability condition was determined implying the stabilizability of the control system through a finite-dimensional controller. In [8], the regulator problem was extended to DPSs for bounded input and output operators, with reference and disturbance signals generated through a finite-dimensional exosystem, providing criteria for the solvability of the regulator equations. The linear regulator problem when considering bounded input and output operators but also bounded disturbance operators entering along the entire interval is shown in [11]. In this last work, the linear regulation problem was solved to the heat equation, damped wave equation, harmonic tracking for a coupled wave equation, control of a damped Rayleigh beam, vibration of a 2D plate, thermal control of a 2D fluid flow, thermal regulation in a 3D room, and control of a linearized Stokes flow in 2D. Reviews about the generalization of the regulator problem to infinite-dimensional systems can be found in [12,13]. The output regulation problem for DPSs has been studied extensively for different classes of PDEs systems; a summary is given in [14]. In [14,15], following the methodology from [11] which is based on the derivation of the transfer function from the system model representation in the state-space form to the solvability of the regulator equations, the SFRP was solved for the R–D equation. An abstract model for the R–D equation was derived where the state operator has the form given by the Sturm–Liouville differential operator (SLDO) plus a parametric term. Simulation results validate their proposal showing the achievement of the regulation tasks to a set-point as well as to harmonic tracking under both set-point disturbance and harmonic disturbance rejection.
The phenomena described by a C–D equation exhibit diffusion and convection properties that are common in many scenarios. Diffusion is the mix of a substance through the medium while convection is the movement of the substance by means of the medium, e.g., when considering smoke rising from a chimney, the smoke particles are convected upward with the air and diffuse within the air currents. It is possible for the convection of the substance to contribute more of a movement in the substance than the diffusion itself. In [16], an illustration about the solutions and behavior for diffusion problems when including the convection term is given. A review of different diffusion models, namely, the Maxwell–Stefan model, the generalized Fick’s law, the classical Fick’s law, and the irreversible thermodynamic model, is given in [17]. The importance for an accurate measuring and prediction of the diffusion coefficients as well as the importance of considering the dragging effect is emphasized. So, an accurate method to approximate the gas–oil mass transfer mechanism based on irreversible thermodynamics was proposed. In this last work, molecular diffusion is only discussed since the system was assumed to be an isothermal one. His proposal is validated through numerical examples and experimental test cases when considering convection for some cases. A chemotaxis–diffusion–convection coupling system which describes a form of buoyant convection in which the fluid develops convection cells and plume patterns is studied in [18]. The pattern formation and hydrodynamical stability of the system was investigated through the development of an upwind finite element method based on a two-dimensional convective chemotaxis–fluid model. Numerical results show the influence of the deterministic initial condition on the overall behavior regarding the number of plumes and that the overall system was stabilized by the chemotaxis. To the best of our knowledge, there is no work about the solvability of the SFRP to a C–D system. In fact, there is not much in the literature about works related to the control of a C–D system.
In this work, our proposal is related to solving the SFRP to the C–D equation. Our main contribution is the definition of the state operator in terms of the SLDO plus an integral term, giving rise to an abstract model for the C–D equation from which the regulator equations have solutions.
The organization of the manuscript is as follows. A summary about the properties of the modeling of DPSs through transfer functions, a brief description of the SFRP for finite-dimensional and infinite-dimensional systems and its application to the R–D equation are given in Section 1. In Section 2, we formulated the problem statement; the design of the regulator is carried out in Section 3; in Section 4, we included simulation results; and the conclusion is given at the end.

2. Problem Statement

2.1. Sturm–Liouville System

Typical problems of mathematical physics lead to Sturm–Liouville eigenvalue and boundary value problems (BVPs). The method of separation of variables in initial BVPs for PDEs lead to Sturm–Liouville eigenvalue and BVPs for ODEs [19]. Most of the problems involve the wave equation or heat (diffusion) equation.
Let us consider the following differential equation
( p x y ( x ) ) q x y ( x ) + λ r x y ( x ) = 0 , a x b ,
subject to symmetric (separated) boundary conditions
α 1 y ( a ) + α 2 y ( a ) = 0 , β 1 y ( b ) + β 2 y ( b ) = 0 ,
with α 1 , α 2 , β 1 and β 2 constant values, p ( x ) , q ( x ) and r ( x ) denoting some given functions, and λ representing a separation constant, which is positive in typical applications whose value represents real eigenvalues [20]. The function y ( x ) will be required to vanish at both ends of the interval. Anyway, separated boundary conditions can be specified in such a way that y ( x ) could vanish at one endpoint and its derivative could vanish at the other. The boundary conditions can be interpreted as defining a Hilbert space Z . The boundary conditions are satisfied in many problems in mathematical physics and are determined by the physical application under study. The BVP given by (1)–(2) is the so-called Sturm–Liouville boundary value problem (SLBVP) or Sturm–Liouville system (SLS) [20,21]. In this case, the BVP is said to be regular.
Let us consider the SLDO [22] given by
L = d d x p x d d x + q x ,
if p ( x ) 0 and both p ( x ) and q ( x ) are continuous on [ a , b ] , then it is said that the SLDO (3) is regular.
Considering the linear homogeneous differential operator
L [ y ( x ) ] = ( p x y ( x ) ) + q x y ( x )
the Sturm–Liouville Equation (1) may be rewritten as
L [ y ( x ) ] = λ r x y ( x ) .
Modeling involving linear second-order ODEs can always be put into the so-called self-adjoint form, which for higher-order equations not always is possible. The way to convert linear second-order ODEs into the self-adjoint form is summarized in the next theorem.
Theorem 1
([23]). Assume that f x > 0 , g x and h x are analytic real-valued functions in the finite (or infinite) interval a < x < b ; then, the existing functions p x > 0 , q x and r x > 0 are similarly analytic and real valued in the same interval such that
1 r x ( p x y ) + q x y = f x y + g x y + h x y
identically in y
Proof. 
See [23]. □
The expression from the right-hand side in (6) is referred to as the self-adjoint form, which is also known as the Sturmian form. The Hilbert space definition of self-adjoint not only depends from the shape of (3) but also from the boundary conditions and scalar product for an unweighted integral from a to b that makes (4) a Hermitian operator [24].
The solutions (eigenfunctions) of the SLS have many properties in common, such as the orthogonality property useful in eigenfunctions expansions in terms of Fourier series, Chebyshev polynomials, Laguerre polynomials, Hermite polynomials, spherical Bessel functions, and many others [14,21,25,26]. The most important property of the eigenfunctions of an SLS is that they form a complete set.
The SLS is an infinite dimensional generalization of the finite dimensional matrix eigenvalue problem
M ψ = λ ψ
with M a n × n matrix and ψ a n-dimensional column vector. As in the matrix case, the SLS will have solutions only for certain values of the eigenvalue λ . The solutions ψ corresponding to this λ are the eigenfunctions. For the finite dimensional case with a n × n matrix M, there can be at most n linearly independent eigenfunctions. In general, for the SLBVP, there will be an infinite set of eigenvalues λ with corresponding eigenfunctions ψ [26].
To the second order case,
L = i , j = 1 n x j a i j x x i a x
with a x 0 for all x Ω , Ω R n is a bounded domain with piecewise smooth boundary, a i j x = a j i x , a i j C ( Ω ¯ ) , where Ω ¯ denotes the closure of Ω , and self-adjointness is determined by the boundary conditions from the differential equation. If there exist constants 0 < c 1 < c 2 < , for all ξ R n and x Ω ¯ , by uniform ellipticity, it implies that
c 1 ξ 2 i , j = 1 n a i j x ξ i ξ j c 2 ξ 2 ,
where · represents the Euclidean norm in R n .

2.2. Abstract Control Model

Consider the abstract evolution (differential) equation
z t ( x , t ) = A z ( x , t ) ,
z ( x , 0 ) = z 0 x ,
where A is the state operator on Z , z 0 means the state of the system defined at time zero, z ( x , t ) means the state of the system at time t, and z t denotes the derivative d z ( t ) / d t .
Definition 1.
If positive constants M and α exist and
z ( · , t ) = e A t z 0 M e α t z 0 t 0 ,
with z 0 Z , is satisfied, so the system (10) and (11) is exponentially stable, which means that A is (exponentially) stable, i.e., A generates an exponentially stable C 0 semigroup in Z [27].
In other words, the uncontrolled state z ( · , t ) = e A t z 0 converges exponentially fast to zero as t [28].
Now, let us consider the abstract control model
z t = A z ( t ) + B i n u ( t ) + B d d ( t ) ,
z ( 0 ) = z 0 , z 0 Z ,
y ( t ) = C z ( t ) ,
where A refers to an unbounded densely defined operator, with domain D ( A ) in Z , u ( t ) U means the input, and y ( t ) Y means the output. U , Y Z may be either finite- or infinite-dimensional. B i n denotes the input operator, B d L ( U , Z ) denotes the disturbance operator, C L ( Z , Y ) denotes the output operator, and d ( t ) refers to the disturbance.
The operator C is a set of bounded output operators C i given by
y i ( t ) = C i z = 1 | Ω i | Ω i z ( x , t ) d x
for some Ω j of the domain Ω with Lebesgue measure
| Ω i | = Ω i d x > 0
of the set Ω j Ω .
More generally,
y i ( t ) = C i z = z , Ψ i = Ω z ( x , t ) Ψ i x d x
with
Ψ i x = 1 | Ω i | 1 Ω i x L 2 ( Ω )
and indicator function
1 Ω i x = 1 , x Ω i 0 , elsewhere .
So,
y = C z = [ C 1 ( z ) C 2 ( z ) C n c ( z ) ] T ,
with n c denoting the number of components in the mixture [17].
The input, output and disturbance operators are bounded operators acting in the interior of the domain. The input to the system is spatially uniform over a small interval about a fixed point x i n = x 0 ( 0 , 1 ) , where B i n u = b x u with
b x = 1 2 ν 0 1 [ x 0 ν 0 , x 0 + ν 0 ] x
and
1 [ a , b ] x = 1 , x [ a , b ] 0 , elsewhere .
The input operators B i n and B d are given as
B i n u ( t ) = j = 1 n i n B i n j u j ( t ) , B d d ( t ) = j = 1 n d B d j d j ( t ) ,
where u j ( t ) and d j ( t ) are scalar control inputs and disturbances, respectively. B i n j and B d j x are characteristic functions of a bounded subset of Ω , namely
B i n j x = 1 | Ω j | 1 Ω j x .
To guarantee that B i n j Z , here, it is assumed that | Ω j | > 0 .
For linear infinite-dimensional systems in the form (13) and (14), it is required for A to be an infinitesimal generator of a C 0 semigroup. Conditions about the characterization of infinitesimal generators are established in the following theorem.
Theorem 2.
Hille–Yosida Theorem. A necessary and sufficient condition for a closed, densely defined, linear operator A on a Hilbert space Z to be the infinitesimal generator of a C 0 semigroup is that there exist real numbers M , w such that for all real α > w , α ρ ( A ) , the resolvent set of A , and
R ( α , A ) r M ( α w ) r r 1 ,
where R ( α , A ) = ( α I A ) 1 is the resolvent operator. In this case
T ( t ) M e w t
with T ( t ) a linear operator.
Proof. 
See [28]. □
Definition 2.
A continuously differentiable function z : [ 0 , ) Z is called a classical solution of (10) and (11) if for all t 0 , we have z ( t ) D ( A ) and (10) and (11) is satisfied.
Lemma 1.
Let A be the infinitesimal generator of a strongly continuous semigroup T ( t ) . If z 0 D ( A ) , then the unique classical solution of (10) and (11) is given by
z ( t ) = T ( t ) z 0 .
Proof. 
See [28]. □
It is worth noting that even when z 0 does not belong to D ( A ) , the function (21) is well defined, so it is said that z : [ 0 , ) Z is the mild solution of (10) and (11). From the above, it is clear that the operator T ( t ) plays the role of e A t in finite-dimensional systems. The strongly continuous semigroup theory on Z is a generalization of the concept of e A t for unbounded operators A on abstract spaces.
Converse to Lemma 1, the next theorem, providing that A has a non-empty resolvent, establishes the property for which the existence of unique classical solutions implies the existence of a strongly continuous semigroup.
Theorem 3.
Let A be a linear operator from D ( A ) Z to Z with λ I A boundedly invertible for some λ C , i.e., ( λ I A ) 1 L ( Z ) . If for all z 0 D ( A ) , the abstract differential Equations (10) and (11) possesses a unique classical solution, then A generates a C 0 semigroup.
Proof. 
See [28]. □

2.3. SFRP

Consider a finite-dimensional neutrally stable exosystem, which generates both reference output y r ( t ) and disturbance d ( t ) , given by
w t = S w ( t ) ,
y r ( t ) = Q w ( t ) ,
d ( t ) = P w ( t ) ,
w ( 0 ) = w 0 ,
where W denotes the state space of the exosystem, S L ( W ) , Q L ( W , Y ) and P L ( W , Z ) .
Let us define the error signal
e ( t ) = y ( t ) y r ( t ) ,
or, equivalently
e ( t ) = C z ( t ) Q w ( t ) .
The main task for the regulator consists of forcing the output of the system to track a reference signal in the presence of a disturbance d ( t ) , i.e., e ( t ) 0 as t . Thus, the problem is stated as follows.
Problem 1.
The SFRP consists of finding a control law
u ( t ) = Γ w ( t )
such that for the system
z t = A z ( t ) + ( B d P + B i n Γ ) w ( t ) ,
w t = S w ( t ) ,
with Γ L ( W , U ) , given by interconnection of (13)–(15) with (22)–(25), the norm of the error
e ( t ) = C z ( t ) Q w ( t )
satisfies
e ( t ) 0 as t
for any z 0 Z and w 0 W .
In view of e ( t ) being a finite-dimensional vector, all l p norms in (32) are equivalent. Since it has been assumed exponential stability for the system (10) and (11), a state feedback control law is not required. In what follows, we state the solvability to the SFRP.
Theorem 4.
If there exist mappings Π L ( W , Z ) and Γ L ( W , U ) , with r a n k ( Π ) D ( A ) , satisfying the regulator equations
Π S = A Π + B i n Γ + B d P ,
C Π = Q ,
the feedback control law that solves the SFRP is given by
u ( t ) = Γ w ( t ) .
Proof. 
The proof can be carried out along the same lines as in [11]. □
From (12), and from the fact that the exosystem is neutrally stable, if (33) holds, then e ( t ) 0 as t for all z 0 Z and w 0 W , if and only if [ C Π Q ] = 0 .

3. Regulator Design

Let us consider the C–D system
z t ( x , t ) = D z x x ( x , t ) B z x ( x , t ) + B i n u + B d d ,
z ( 0 , t ) = 0 , ( Dirichlet BC )
z x ( 1 , t ) = 0 , ( Neumann BC )
z ( x , 0 ) = ϕ x ,
y ( t ) = C z ( t ) ,
with D z x x ( x , t ) the diffusion term, with diffusion coefficient D > 0 , and B z x ( x , t ) the convection term, with convection coefficient B > 0 , z x x denotes the second partial derivative and z x denotes the first partial derivative both with respect to space.
In our work, the system (36)–(40) is defined in the abstract form (13)–(15) in the Hilbert state space Z = L 2 ( 0 , 1 ) . The maximal elliptic operator is given by L = d 2 / d x 2 in (8) belonging to D ( L ) = H 2 ( 0 , 1 ) , indicating the Sobolev space of functions in Z with a square integrable second derivative.
So, from (8), here, the state operator is defined in the form of the SLDO as given by
A = A A
with
A = d 2 / d x 2 ,
this latter expression (operator) an infinitesimal generator of a C 0 semigroup for abstract differential equations related with parabolic PDEs as the heat (diffusion) equation [10,11].
Let us assume that the state operator (41) is a self-adjoint (Hermitian) operator in Z , i.e.,
A [ ζ ] = ζ ζ
with
D ( A ) = { ζ H 2 ( 0 , 1 ) ζ ( 0 ) = 0 , ζ ( 1 ) = 0 } Z ,
where, because for z 0 D ( A ) the mild solution is a classical solution, the symmetric boundary conditions (2) are part of the domain of A .
The spectrum of A denoted by
σ ( A ) = { λ k } k = 0 ,
where λ k = k 2 π 2 is purely discrete with a set of orthonormal eigenvectors
ζ k x = 2 cos ( k π x ) , k = 1 , 2 , ,
and A is an infinitesimal generator in terms of the eigenfunction expansion
e A t ζ = j = 0 e λ j t ζ , ζ j ζ j ,
which gives rise to an orthonormal basis for L 2 ( 0 , 1 ) .
The system (36)–(40) is a single-input/single-output system with scalar input B i n , output C, and disturbance B d . The output is the average transport reaction over a small interval about a point x o u t = x 1 ( 0 , 1 ) , i.e.,
C ϕ = 0 1 c x ϕ x d x
with
c x = 1 2 ν 1 1 [ x 1 ν 1 , x 1 + ν 1 ] x .
Since C ζ = ζ , c , C is a bounded linear observation functional on Z . In our work, d ( t ) = A d sin ( β t ) R entering across the entire interval is considered, so B d = 1 .

3.1. Trajectory Tracking in Presence of a Constant Disturbance

In our proposal, the regulator (33) and (34) takes the form
Π S w = A Π w + B i n Γ w + B d P w ,
C Π w = Q w = w 1 ,
with
Π = [ Π 1 Π 2 Π 3 ] , w ( t ) = w 1 w 2 w 3 = A r sin ( α t ) A r cos ( α t ) M d ,
S = 0 α 0 α 0 0 0 0 0 , Q = [ 1 0 0 ] , P = [ 0 0 1 ] .
The block diagonal matrix S allows us to decouple the regulator equations. To achieve this, when trying with the trajectory tracking, the regulator equations can be written as
Π α S α w α = A Π α w α + B i n Γ α w α , C Π α w α = w 1
with
Π α = [ Π 1 , Π 2 ] , S α = 0 α α 0 , w α = w 1 w 2 = A r sin ( α t ) A r cos ( α t ) .
Defining P α = [ 0 0 ] and Q α = [ 1 0 ] , then
y r ( t ) = Q α w α = A r sin ( α t ) .
In this case W = R 2 , so Π α = [ Π 1 Π 2 ] and Γ α = [ Γ 1 Γ 2 ] R 2 with Π j Z . The regulator equations applied to the vector w α = [ w 1 w 2 ] T W are then given as
Π S w = A Π w + B i n Γ w ,
C Π w = Q w .
From these last equations, expanding the regulator equation results in
α Π 1 w 2 α Π 2 w 1 = A Π 1 w 1 + A Π 2 w 2 + B i n Γ 1 w 1 + B i n Γ 2 w 2 .
Since (49) must be satisfied for all w, let us first consider the event for which w 1 = 1 and w 2 = 0 and then for w 1 = 0 and w 2 = 1 yielding
α Π 2 A Π 1 = B i n Γ 1 ,
α Π 1 A Π 2 = B i n Γ 2 .
Recalling that the exosystem (22)–(25) is neutrally stable, multiplying (51) by i, an imaginary number, and adding the result to (50) results in
( i α I A ) Π 1 + i ( i α I A ) Π 2 = B i n i Γ 2 + B i n Γ 1 .
Since i α ρ ( A ) , premultiplying both sides of (52) by ( i α I A ) 1 yields
Π 1 + i Π 2 = ( i α I A ) 1 B i n ( Γ 1 + i Γ 2 ) .
Premultiplying by C both sides of (53) and from the identity C Π w = Q w , with C Π 1 = 1 and C Π 2 = 0 , then Q = [ 1 0 ] . Consequently,
1 = C ( i α I A ) 1 B i n ( i Γ 2 + Γ 1 ) = G ( i α ) ( Γ 1 + i Γ 2 ) .
From the identity G ( i α ) = Re ( G ( i α ) ) + i Im ( G ( i α ) ) , (54) is rewritten as
1 = ( Re ( G ( i α ) ) + i Im ( G ( i α ) ) ) ( Γ 1 + i Γ 2 ) .
From the above, matching real and imaginary parts,
1 = ( Re ( G ( i α ) ) Γ 1 Im ( G ( i α ) ) Γ 2 , 0 = ( Im ( G ( i α ) ) Γ 1 + Re ( G ( i α ) ) Γ 2 .
Thus,
Γ 1 = Re ( G ( i α ) ) | G ( i α ) | 2 , Γ 2 = Im ( G ( i α ) ) | G ( i α ) | 2 .
Accordingly, from the definition
G ( i α ) 1 = 1 G ( i α ) = G ( i α ) ¯ | G ( i α ) | 2 = Re ( G ( i α ) ) i Im ( G ( i α ) ) | G ( i α ) | 2 ,
the control gains are given by
Γ α = [ Γ 1 Γ 2 ] = [ Re ( G ( i α ) 1 ) Im ( G ( i α ) 1 ) ] .
Here, the system has been assumed to be real, i.e., G ( s ¯ ) = G ( s ) ¯ for all s ρ ( A ) . It is worthwhile to mention that G ( i α ) must differ from zero as well as be invertible for solvability.
Trying now with the rejection of a constant disturbance, the regulator equations are given as
0 = A Π 3 w 3 + B d w 3 + B i n Γ 3 w 3 , C Π 3 w 3 = Q 3 w 3 = 0 ,
where Q 3 = 0 and P 3 = 1 . Thus, the regulator equations become
0 = A Π 3 + B d + B i n Γ 3 , C Π 3 = 0 .
Solving for Γ , the last system of equations then
Π 3 = A 1 ( B d + B i n Γ 3 ) , 0 = C ( A ) 1 B d + C ( A ) 1 B i n Γ 3 .
So,
Γ 3 = C ( A ) 1 B d C ( A ) 1 B i n = G B d ( 0 ) G ( 0 ) .
Again, assuming G ( 0 ) 0 and using the following definition
G B d ( s ) = C ( s I A ) 1 B d ,
consequently, combining (55) and (56) results in
Γ = [ Γ α , Γ 3 ] = [ Γ 1 , Γ 2 , Γ 3 ] = [ Re ( G ( i α ) 1 ) , Im ( G ( i α ) 1 ) , G ( 0 ) 1 ( C A 1 B d ) ] .

3.2. Trajectory Tracking in Presence of Harmonic Disturbance

Along the same lines as that for the trajectory tracking control problem with rejection of a constant disturbance but now trying with the rejection of a harmonic disturbance, from (45) and (46) with
Π = [ Π 1 Π 2 Π 3 Π 4 ] , w ( t ) = w 1 w 2 w 3 w 4 = A r sin ( α t ) A r cos ( α t ) A d sin ( β t ) A d cos ( β t ) ,
S = 0 α 0 0 α 0 0 0 0 0 0 β 0 0 β 0 , Q = [ 1 0 0 0 ] , P = [ 0 0 1 0 ] ,
and, from the block diagonal matrix S, decoupling the regulator equations when considering the case of trajectory tracking with rejection of a harmonic disturbance
Π α S α w α = A Π α w α + B i n Γ α w α , C Π α w α = w 1 ,
with
Π α = [ Π 1 , Π 2 ] , S α = 0 α α 0 , w α = w 1 w 2 , Q α = [ 1 , 0 ] , P α = [ 0 , 0 ] ,
thus, the solution is given by (55). In addition, referring us to the blocks in S to try with the rejection of harmonic disturbance, where d = A d sin ( β t ) , then
Π β S β w β = A Π β w β + B d P β w β + B i n Γ β w β , C Π β w β = 0 ,
where
Π β = [ Π 3 , Π 4 ] , S β = 0 β β 0 , w β = w 3 w 4 , Q β = [ 0 , 0 ] , P β = [ 1 , 0 ] .
Again, to this case W = R 2 so, looking for Π = [ Π 3 Π 4 ] , where Π j Z , and Γ = [ Γ 3 Γ 4 ] R 2 , the regulator equations applied to the vector w = [ w 3 w 4 ] T W result in the system
Π S w = A Π w + B d P w + B i n Γ w ,
C Π w = 0 .
Expanding (57) yields
β Π 3 w 4 β Π 4 w 3 = A Π 3 w 3 + A Π 4 w 4 + B d w 3 + B i n Γ 3 w 3 + B i n Γ 4 w 4 .
Since (59) must be satisfied for w, first consider the event for w 3 = 1 and w 4 = 0 to then consider the episode for w 3 = 0 and w 4 = 1 , thus
β Π 4 A Π 3 = B d + B i n Γ 3 ,
β Π 3 A Π 4 = B i n Γ 4 .
Multiplying (61) by the imaginary number i and then adding the result to (60) yields
( i β I A ) Π 3 + i ( i β I A ) Π 4 = B d + B i n i Γ 4 + B i n Γ 3 .
Noting that i β ρ ( A ) , premultiplying (62) by ( i β I A ) 1 , it becomes
Π 3 + i Π 4 = ( i β I A ) 1 B d + ( i β I A ) 1 B i n ( i Γ 4 + Γ 3 ) .
Thus, premultiplying by C both sides of (63) and from the fact that (58) implies
C Π 3 = 0 , C Π 4 = 0 ,
so, results
0 = C ( i β I A ) 1 B d + C ( i β I A ) 1 B i n ( i Γ 4 + Γ 3 ) , = G B d ( i β ) + G ( i β ) ( i Γ 4 + Γ 3 ) ,
where the definition
G B d ( s ) = C ( s I A ) 1 B d
was used. At last, solving for Γ yields
Γ 3 = Re ( G ( i β ) 1 ) Re ( G B d ( i β ) ) + Im ( G ( i β ) 1 ) Im ( G B d ( i β ) ) ,
Γ 4 = Re ( G ( i β ) 1 ) Im ( G B d ( i β ) ) Im ( G ( i β ) 1 ) Re ( G B d ( i β ) ) ,
where
Γ = [ Γ 3 Γ 4 ] .
Hence, combining (55) and (64) yields
Γ = [ Γ α Γ β ] = [ Γ 1 Γ 2 Γ 3 Γ 4 ]
with
Γ α = [ Γ 1 , Γ 2 ] = [ Re ( G ( i α ) 1 ) , Im ( G ( i α ) 1 ) ] , Γ β = [ Γ 3 , Γ 4 ] = [ Re ( G ( i β ) 1 ) Re ( G B d ( i β ) ) + Im ( G ( i β ) 1 ) Im ( G B d ( i β ) ) , Re ( G ( i β ) 1 ) Im ( G B d ( i β ) ) Im ( G ( i β ) 1 ) Re ( G B d ( i β ) ) ] .

4. Simulation Results

In order to validate our proposal via numerical simulation, for the case of trajectory tracking under the influence of a constant disturbance, we have set M d = 10 , A r = 3 , α = 8 , x 0 = 0.75 , x 1 = 0.25 , and ν 0 = ν 1 = 0.5 . Figure 1a shows the tracking of the reference signal y r ( t ) by the output y ( t ) from the initial condition φ x = 4 cos ( π x ) . Figure 1b shows the error signal e ( t ) between the controlled output y ( t ) and the reference signal y r ( t ) from which it can be seen that e ( t ) 0 as t . The solution surface is shown in Figure 2.
To the case of trajectory tracking under the influence of harmonic disturbance, we set M d = 1 , A r = 3 , A d = 3 , α = 12 , β = 12 , x 0 = 0.75 , x 1 = 0.25 and ν 0 = ν 1 = 0.5 . Figure 3a shows the tracking of the reference signal y r ( t ) by the controlled output y ( t ) for the initial condition φ ( x ) = 4 cos ( π x ) . Figure 3b exhibits that e ( t ) 0 as t . The corresponding solution surface is shown in Figure 4.
So, the regulator performs well under the presence of external disturbances for both cases, i.e., in the presence of either a constant disturbance or harmonic disturbance.

5. Conclusions

In our work, the SFRP approach is focused on the trajectory tracking control with the rejection of external disturbances to the C–D equation. The C–D system is modeled through a state operator given in the form of the SLDO plus an integral term involved in an abstract control system model from which the regulator equations are derived and solved. From the simulation results, it is concluded that our proposal performs well since when considering both constant and harmonic disturbances, the regulator is capable of tracking the reference trajectory, showing the rejection of external disturbances. As future work, we are focused on extend our proposal to multiple-input/multiple-output systems.

Author Contributions

Conceptualization, A.A.R. and F.J.; methodology, A.A.R. and F.J.; software, A.A.R.; validation, A.A.R.; formal analysis, A.A.R. and F.J.; investigation, A.A.R. and F.J.; data curation, A.A.R.; writing—original draft preparation, A.A.R. and F.J.; writing—review and editing, F.J.; visualization, A.A.R.; resources, F.J.; supervision, F.J.; funding acquisition, F.J.; project administration, F.J. All authors have read and agreed to the published version of the manuscript.

Funding

This research work was funded by Tecnológico Nacional de México (TecNM) projects and, partially, under grant number 39873 from EDD 2022 program. This work was supported by CONACYT, México, through grant 862135.

Data Availability Statement

Data available on request due to restrictions eg privacy or ethical. The data presented in this study are available on request from the corresponding author. The data are not publicly available due to legal reasons.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:
DPSsDistributed Parameters Systems
ODEsOrdinary Differential Equations
PDEsPartial Differential Equations
C–DConvection–Diffusion
R–DReaction–Diffusion
SFRPState Feedback Regulation Problem
SLBVPSturm–Liouville Boundary Value Problem
SLDOSturm–Liouville Differential Operator
SLSSturm–Liouville System

References

  1. Morris, K.A. Controller Design for Distributed Parameter Systems; Springer: Berlin/Heidelberg, Germany, 2020; pp. 59–67. [Google Scholar]
  2. Morris, K. Control of Systems Governed by Partial Differential Equations. In The Control Handbook, 2nd ed.; Levine, W.S., Ed.; CRC Press: Boca Raton, FL, USA, 2011; pp. 67-1–67-37. [Google Scholar]
  3. Isidori, A. Nonlinear Control Systems; Springer: London, UK, 1995. [Google Scholar]
  4. Pohjolainen, S.A. On the asymptotic regulation problem for distributed parameter systems. In Control of Distributed Parameter Systems; Pergamon: Oxford, UK, 1983; pp. 197–201. [Google Scholar]
  5. Pohjolainen, S.A. Robust multivariable PI–controller for infinite dimensional systems. IEEE Trans. Autom. Control 1982, 27, 17–30. [Google Scholar] [CrossRef]
  6. Schumacher, J.M. Dynamic feedback in finite– and infinite–dimensional linear systems. In Mathematisch Centrum Tracts; Centrum Voor Wiskunde en Informatica: Amsterdam, The Netherlands, 1981; Volume 143. [Google Scholar]
  7. Schumacher, J.M. Finite–dimensional regulators for a class of infinite–dimensional systems. Syst. Control. Lett. 1983, 3, 7–12. [Google Scholar] [CrossRef]
  8. Byrnes, C.I.; Laukó, I.G.; Gilliam, D.S.; Shubov, V.I. Output regulation for linear distributed parameter systems. IEEE Trans. Autom. Control 2000, 45, 2236–2252. [Google Scholar]
  9. Müller, P.H. “T. Kato, Perturbation theory for linear operators.(Die Grundlehren der mathematischen Wissenschaften in Einzeldarstellungen mit besonderer Berücksichtigung der Anwendungsgebiete, Band 132) XX+ 592 S. m. 3 Fig. Berlin/Heidelberg/New York Springer-Verlag. Preis geb. DM 79, 20”. Z. Angew. Math. Mech. 1967, 47, 554. [Google Scholar]
  10. Curtain, R.F.; Zwart, H. An Introduction to Infinite-Dimensional Linear Systems Theory; Springer: New York, NY, USA, 1995. [Google Scholar]
  11. Aulisa, E.; Gilliam, D. A Practical Guide to Geometric Regulation for Distributed Parameter Systems; CRC Press: Boca Raton, FL, USA, 2016. [Google Scholar]
  12. Natarajan, V.; Gilliam, D.S.; Weiss, G. The state feedback regulator problem for regular linear systems. IEEE Trans. Autom. Control 2014, 59, 2708–2723. [Google Scholar] [CrossRef]
  13. Xu, X.; Dubljevic, S. Output and error feedback regulator designs for linear infinite–dimensional systems. Automatica 2017, 83, 170–178. [Google Scholar] [CrossRef]
  14. Jurado, F.; Ramírez, A.A. State feedback regulation problem to the reaction-diffusion equation. Mathematics 2020, 8, 1983. [Google Scholar] [CrossRef]
  15. Ramírez, A.A.; Jurado, F. Tracking regulator design with disturbance rejection to the reaction–diffusion equation. In Proceedings of the 2019 16th International Conference on Electrical Engineering, Computing Science and Automatic Control (CCE), Mexico City, Mexico, 11–13 September 2019; pp. 11–13. [Google Scholar]
  16. Farlow, S.J. Partial Differential Equations for Scientists and Engineers; Dover Publications: New York, NY, USA, 1993. [Google Scholar]
  17. Hoteit, H. Modeling diffusion and gas–oil mass transfer in fractured reservoirs. J. Pet. Sci. Eng. 2013, 105, 1–17. [Google Scholar] [CrossRef]
  18. Deleuze, Y.; Chiang, C.-Y.; Thiriet, M.; Sheu, T.W.H. Numerical study of plume patterns in a chemotaxis–diffusion–convection coupling system. Comput. Fluids 2016, 126, 58–70. [Google Scholar] [CrossRef]
  19. Guenther, R.B.; Lee, J.W. Sturm–Liouville Problems Theory and Numerical Implementation; CRC Press: Boca Raton, FL, USA, 2019. [Google Scholar]
  20. Spiegel, M.R. Fourier Analysis with Applications to Boundary Value Problems; McGraw–Hill: New York, NY, USA, 1974. [Google Scholar]
  21. Kreyszig, E. Advanced Engineering Mathematics; John Wiley & Sons: Hoboken, NJ, USA, 2006. [Google Scholar]
  22. Mathews, J. Walker, R.L. Mathematical Methods of Physics; The Benjamin/Cummings Publishing Co.: Menlo Park, CA, USA, 1970. [Google Scholar]
  23. Holland, S.S., Jr. Applied Analysis by the Hilbert Space Method, An Introduction with Applications to the Wave, Heat, and Schrödinger Equations; Dover Publications: Mineola, NY, USA, 1990. [Google Scholar]
  24. Arfken, G.B.; Weber, H.J.; Harris, F.E. Mathematical Methods for Physicists A Comprehensive Guide; Academic Press: Waltham, MA, USA, 2013. [Google Scholar]
  25. Kaplan, K. Advanced Calculus; Pearson Education: Boston, MA, USA, 2003. [Google Scholar]
  26. Wyld, H.W. Mathematical Methods for Physics; CRC Press: Boca Raton, FL, USA, 2021. [Google Scholar]
  27. Pazy, A. Semigroups of Linear Operators and Applications to Partial Differential Equations; Springer: New York, NY, USA, 1983. [Google Scholar]
  28. Curtain, R.; Zwart, H. Introduction to Infinite-Dimensional Systems Theory; Springer: New York, NY, USA, 2020. [Google Scholar]
Figure 1. Regulator performance with Γ = [ 4.9348 , 8 , 10 ] : (a) Comparison of the output y ( t ) and the reference y r ( t ) , (b) Error between y ( t ) and y r ( t ) .
Figure 1. Regulator performance with Γ = [ 4.9348 , 8 , 10 ] : (a) Comparison of the output y ( t ) and the reference y r ( t ) , (b) Error between y ( t ) and y r ( t ) .
Mathematics 11 01944 g001
Figure 2. Spatial distribution of the solution surface corresponding to the rejection of constant disturbance.
Figure 2. Spatial distribution of the solution surface corresponding to the rejection of constant disturbance.
Mathematics 11 01944 g002
Figure 3. Regulator performance with Γ = [ 4.9348 , 12 , 1 , 0 ] : (a) Comparison of the output y ( t ) and the reference signal y r ( t ) , (b) Error between y ( t ) and y r ( t ) .
Figure 3. Regulator performance with Γ = [ 4.9348 , 12 , 1 , 0 ] : (a) Comparison of the output y ( t ) and the reference signal y r ( t ) , (b) Error between y ( t ) and y r ( t ) .
Mathematics 11 01944 g003
Figure 4. Spatial distribution of the solution surface related with the rejection of harmonic disturbance.
Figure 4. Spatial distribution of the solution surface related with the rejection of harmonic disturbance.
Mathematics 11 01944 g004
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Ramírez, A.A.; Jurado, F. The Regulator Problem to the Convection–Diffusion Equation. Mathematics 2023, 11, 1944. https://doi.org/10.3390/math11081944

AMA Style

Ramírez AA, Jurado F. The Regulator Problem to the Convection–Diffusion Equation. Mathematics. 2023; 11(8):1944. https://doi.org/10.3390/math11081944

Chicago/Turabian Style

Ramírez, Andrés A., and Francisco Jurado. 2023. "The Regulator Problem to the Convection–Diffusion Equation" Mathematics 11, no. 8: 1944. https://doi.org/10.3390/math11081944

APA Style

Ramírez, A. A., & Jurado, F. (2023). The Regulator Problem to the Convection–Diffusion Equation. Mathematics, 11(8), 1944. https://doi.org/10.3390/math11081944

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

Article Metrics

Back to TopTop