1. Introduction
Many results in the theory of asymptotic approximations have been obtained from 1930 onwards. Indeed, there were a lot of results on integral manifolds, equations with retarded argument, quasi- or almost-periodic equations etc. Earlier work on this theory has been presented in the famous book [
1].
Averaging is a valuable method to understand the long-term evolution of dynamical systems characterized by slow dynamics and fast periodic or quasi-periodic dynamics. In [
2], a transparent proof of the validity of averaging in the periodic case is presented. Different proofs for both the periodic and the general case are provided by [
3,
4]. In the last paper, moreover, the relation between averaging and the multiple time-scales method is established.
The averaging method for constructing approximate solutions in the theory of ODEs is presented in [
5,
6]. In [
7], the asymptotic analysis of nonlinear dynamical systems is developed.
The work [
8] is devoted to using an asymptotic method for studying the Cauchy problem for a 1D Euler–Poisson system, which represents a physically relevant hydrodynamic model but also a challenging case for a bipolar semiconductor device by considering two different pressure functions. In [
9], the averaging results for ordinary differential equations perturbed by a small parameter are proved. Here, authors assume only that the right-hand sides of the equations are bounded by some locally Lebesgue integrable functions with the property that their indefinite integrals satisfy a Lipschitz-type condition.
In [
10], the authors prove that averaging can be applied to the extremal flow of optimal control problems with two fast variables, that is considerably more complex because of resonances.
The averaging method is one of the most effective tools for constructing approximate solutions, including optimal control problems for ODEs [
11] and PDEs [
12], where the autors consider the optimal control problem in coefficients in the so-called class of H-admissible solutions.
The Krasnoselski–Krein theorem and its various modifications [
13,
14,
15] play an essential role in all such considerations, since it guarantees the limit transition in perturbed problem with fast-oscillating coefficients of the form
as
.
The typical averaging problem may be defined as follows: one considers an unperturbed problem in which the slow variables remain fixed. Upon perturbation, a slow drift appears in these variables which one would like to approximate independently of the fast variables.
In the present paper we use this approach to nonlinear parabolic system with fast-oscillating (w.r.t. time variable) coefficients on an infinite time interval. We prove that the optimal control of the problem with averaging coefficients can be considered to be ”approximately” optimal for the initial perturbed system.
2. Statement of the Problem
Let
be a bounded domain. In cylinders
we consider an initial boundary-value problem for a parabolic system [
16,
17]
Here
is a small parameter,
A is a real
matrix,
f is a given vector-valued mapping,
g is a given matrix-valued mapping,
is an unknown state function,
is an unknown control function, which are determined by requirements
where
are positive constants.
Under the natural assumptions on
we prove, that the optimal control problem (
1)–(
3) has a solution
, i.e., for every
and for any solution
of (
1) with control
u we have
In what follows we consider the problem of finding an approximate solution of (
1)–(
3) by transition to averaged coefficients. For this purpose we assume that uniformly w.r.t.
there exists
We consider the following optimal control problem
It should be noted that the transition to the averaging parameters can essentially simplify the problem. In particular, if
does not depend on
y then in some cases exact solution of (
1)–(
3) can be found [
18,
19]. Another approaches for finding exact solutions of optimal control problems and approximate procedures can be found in [
20,
21].
Assume that
is a solution of (
5)–(
7). The main goal of the paper is to prove the limit equality
As a consequence of (
8) we will prove that the control
is approximately optimal for the problem (
1)–(
3) in the following sense:
where
is a solution of (
1) with control
.
3. Assumptions, Notations and Basic Results
We assume the following conditions hold.
Assumption 1. , where and I is a unit matrix;
Assumption 2. is continuous and bounded: Assumption 3. is continuous and bounded: Assumption 4. is closed and convex,
Assumption 5. is a Carathéodory function, , such that , Here, denotes the Euclidean norm of .
For
and
we understand solution of (
1) in weak (or generalized) sense on every finite time interval, i.e.,
y is a solution of (
1) if
such that
,
,
the following equality holds:
Here and after we denote by and the classical norm and scalar product in , by and the classical norm and scalar product in , by the norm in , and by the dual space to V.
Due to the Assumptions 1–3, every solution of (
1) satisfies
It means that
every solution of (
1) is an absolutely continuous function from
to
H, and equality (
9) is equivalent to the following one [
16]:
for almost all (a.a.)
.
It is known [
16,
17] that, under Assumptions 1–3, for every
,
there exists at least one solution of (
1), and for a.a.
Remark 1. Uniqueness of solution of (1) is not guaranteed. This can be done under some additional assumptions, e.g., [16] , , In the sequel, we denote by
(or
) a set of all pairs
, where
y is a solution of (
1) (or (
5)) with control
u.
The following Lemma gives us a result about the solvability of the optimal control problem (
1)–(
3) and it also provides some useful inequalities.
Lemma 1. Under the Assumptions 1–5 for every the problem (1)–(3) has a solution , that is, Proof of Lemma 1. Fix
and suppress index
throughout the proof. The idea of the proof is to derive a priori estimates for the minimizing sequence. Obtained estimates allow us to pass to the limit in problem (
1)–(
3).
From (
11), Poincare inequality
,
, and Young inequality we derive that for some
,
(not depending on
) for every
for a.a.
Therefore using Gronwall inequality we get for all
From the inequality (
13) and the first inequality from the Assumption 5 we have that for some
Now let
be a minimizing sequence, that is,
Note that due to the Assumption 5
From (
15) for sufficiently large
nOn the other hand, for
Inequalities (
16) and (
17) imply that
is bounded in
, so for subsequence
Due to convexity of
U we have inclusion
. From (
11) over
and using (
13) we we obtain from (
5) that
is bounded in
is bounded in
. Using Compactness Lemma [
22] we conclude that up to subsequence
From (
19) and Lebesgue’s Dominated Convergence Theorem we can pass to the limit in the equality (
9) applied to
, and obtain that
. Due to pointwise convergence
Fatou’s lemma and weak convergence (
18) we have
Therefore
is a solution of (
1)–(
3). □
4. Main Results
We assume that
,
Assumption (
20) implies that the averaged function
from (
4) is a continuous function and the Assumption 2 holds. It means that under conditions (
4), (
20) the optimal control problem (
5)–(
7) has a solution
.
The main result of the paper is the following
Theorem 1. Suppose that the Assumptions 1–5 and (4), (20) hold and, moreover, the problem (5) has a unique solution for every . Let be a solution of (1)–(3). Thenand up to subsequencewhere is a solution of (5)–(7). Proof. Let
,
be a solution of (
1)–(
3) for
. Due to the optimality of
we have
where
is a solution of (
1) with
and
. Then from (
14)
Repeating arguments used in the prof of Lemma 1 conclude that on subsequence for some
,
:
Let us prove that
, i.e.,
is a solution of the averaged problem (
5) with control
. For this purpose it is sufficient to make a limit transition in the equality
for arbitrary
and
.
Limit transition in the left part of (
24) is a direct consequence of (
23). From the Dominated Convergence Theorem we see that
Then (
23) implies convergence in the last term of (
24).
Let us prove that
,
where
. Due to the Dominated Convergence Theorem
,
Due to Egorov’s theorem
such that
and
Here
is Lebesgue’s measure on
. On the other hand there exists a sequence of step functions
is a covering of
such that
Moreover
such that
and
Due to (
20) for a given
,
Therefore, choosing
such that
we get from (
28) that
On the other hand, for every step function
we have due to (
26):
So
,
,
Furthermore,
,
,
Combining (
30)–(
33), we obtain
,
Inequalities (
29), (
34) imply (
25). So we can pass to the limit in (
24) and obtain that
. Now let us prove that
is an optimal process in (
5)–(
7).
On the other hand, for every
and any
–solution of (
1) with control
u and
we get
Using the same arguments as in proof of the Lemma 1 for
we derive that
where
y is a unique solution of (
5) with control
u.
Indeed due to the Assumption 5 and (
13) we have
As
in
and a.e. in
Q, we deduce from Lebesgue’s Dominated Convergence theorem:
On the other hand, from (
12) and (
36)
where
does not depend on
T and
n. The last inequality together with with (
37) leads to (
35).
From (
35) we conclude the following inequality:
This means that
is a solution of (
5)–(
7).
Now we substitute
in previous arguments. Then
due to uniqueness. So from (
39), we obtain
These inequalities mean that up to subsequence
Since
, then convergence in (
40) holds for the whole sequence. Therefore (
21) is proved.
Moreover, up to subsequence
tends to
in
. So, repeating arguments (
37) and (
38) for
, and using boundness of
, we have
Then from (
40) and weak convergence we deduce (
22) □
Corollary 1. An optimal control of the averaged problem (5)–(7) can serve as an ”approximate” optimal control in the initial problem (1), that is:where is a solution of (1) with control . Indeed, for
,
, we can repeat arguments of the proof of the Theorem, and due to the uniqueness of the solution of (
5) for
we have up to subsequence
Then (
35) holds and taking into account strong convergence (
22), we obtain (
41).
5. Conclusions and Future Research
We sought to obtain a theoretical result that demonstrates the effectiveness of the averaging method of finding an approximate solution of the optimal control problem for a non-linear parabolic system with fast-oscillating coefficients with respect to a time variable. We proved that the optimal control of the problem with averaging coefficients can be considered as an ”approximately” optimal for the initial perturbed system. To demonstrate effectiveness of the method we plan to continue research focusing on the practical applications and simulation results using in particular genetic algorithms.
Author Contributions
Conceptualization, O.A.K., O.V.K., A.R. and V.S.; methodology, O.A.K., O.V.K., A.R. and V.S.; formal analysis, O.A.K., O.V.K., A.R. and V.S.; investigation, O.A.K., O.V.K., A.R. and V.S.; writing—original draft preparation, O.A.K., O.V.K., A.R. and V.S.; writing—review and editing, O.A.K., O.V.K., A.R. and V.S. All authors have read and agreed to the published version of the manuscript.
Funding
This research received no external funding.
Institutional Review Board Statement
Not applicable.
Informed Consent Statement
Not applicable.
Data Availability Statement
Not applicable.
Conflicts of Interest
The authors declare no conflict of interest.
References
- Bogoliubov, N.N.; Mitropolsky, Y.A. Asymptotic Methods in the Theory of Non-Linear Oscillations; Gordon and Breach: New York, NY, USA, 1961. [Google Scholar]
- Roseau, M. Vibrations nonlinéaires et théorie de la stabilité; Springer: Berlin, Germany, 1966. [Google Scholar]
- Besjes, J.G. On the asymptotic methods for non-linear differential equations. J. Mécanique 1969, 8, 357–373. [Google Scholar]
- Perko, L.M. Higher order averaging and related methods for perturbed periodic and quasi-periodic systems. SIAM J. Appl. Math. 1969, 17, 698–724. [Google Scholar] [CrossRef]
- Lochak, P.; Meunier, C. Multiphase Averaging for Classical Systems; Springer: New York, NY, USA, 1988. [Google Scholar]
- Samoilenko, A.M.; Stanzhitskyi, A.N. On averaging differential equations on an infinite interval. Differ. Uravn. 2006, 42, 476–482. [Google Scholar] [CrossRef]
- Sanders, J.A.; Verhulst, F. Averaging Methods in Nonlinear Dynamical Systems; Springer: New York, NY, USA, 1985. [Google Scholar]
- Donatella Donatelli, D.; Mei, M.; Rubino, B.; Sampalmieri, R. Asymptotic behavior of solutions to Euler–Poisson equations for bipolar hydrodynamic model of semiconductors. J. Differ. Equ. 2013, 255, 3150–3184. [Google Scholar] [CrossRef]
- Lakrib, M.; Kherraz, T.; Bourada, A. Averaging for ordinary differential equations perturbed by a small parameter. Math. Bohem. 2016, 141, 143–151. [Google Scholar] [CrossRef] [Green Version]
- Dell’Elce, L.; Caillau, J.-B.; Pomet, J.-B. Considerations on Two-Phase Averaging of Time-Optimal Control Systems. Available online: https://hal.inria.fr/hal-01793704v3 (accessed on 31 March 2022).
- Nosenko, T.V.; Stanzhytskyi, O.M. Averaging method in some problems of optimal control. Nonlin. Osc. 2008, 11, 539–547. [Google Scholar] [CrossRef]
- Ciro D’Apice, C.; De Maio, U.; Kogut, O.P. Optimal Control Problems in Coefficients for Degenerate Equations of Monotone Type: Shape Stability and Attainability Problems. SIAM J. Control Optim. 2012, 50, 1174–1199. [Google Scholar] [CrossRef]
- Kichmarenko, O.; Stanzhytskyi, O. Sufficient conditions for the existence of optimal control for some classes of functional-differential equations. Nonlin. Dyn. Syst. Theory 2018, 18, 196–211. [Google Scholar]
- Plotnikova, N.V. The Krasnoselskii-Krein theorem for differential inclusions. Differ. Uravn. 2005, 41, 997–1000. [Google Scholar]
- Gamma, R.; Guerman, A.; Smirnov, G. On the asymptotic stability of discontinuous systems via the averaging method. Nonlin. Ann. 2011, 74, 1513–1522. [Google Scholar] [CrossRef] [Green Version]
- Chepyzhov, V.V.; Visnik, M.I. Attractors for Equations of Mathematical Physics; AMS: Providence, RI, USA, 2002. [Google Scholar]
- Kapustyan, O.V.; Kasyanov, P.O.; Valero, J. Structure of the global attractor for weak solutions of a reaction-diffusion equation. Appl. Math. Inf. Sci. 2015, 9, 2257–2264. [Google Scholar]
- Kapustian, O.A.; Sobchuk, V.V. Approximate homogenized synthesis for disturbed optimal control problem with superposition type cost functional. Stat. Opt. Inf. Comp. 2018, 6, 233–239. [Google Scholar]
- Kapustian, E.A.; Nakonechny, A.G. The minimax problems of pointwise observation for a parabolic boundary-value problem. J. Autom. Inf. Sci. 2002, 34, 52–63. [Google Scholar]
- Pichkur, V.V.; Sobchuk, V.V. Mathematical Model and Control Design of a Functionally Stable Technological Process. Diff. Eq. App. 2021, 29, 32–41. [Google Scholar] [CrossRef]
- Garashchenko, F.G.; Pichkur, V.V. Structural optimization of dynamic systems by use of generalized Bellman’s principle. J. Autom. Inf. Sci. 2000, 32, 1–6. [Google Scholar]
- Sell, G.R.; You, Y. Dynamics of Evolutionary Equations; Springer: New York, NY, USA, 2002. [Google Scholar]
| Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations. |
© 2022 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).