The theory of contrasting structures in singularly perturbed boundary problems for nonlinear parabolic partial differential equations is applied to the research of formation of steady state distributions of power within the nonlinear “power-society” model. The interpretations of the solutions to the equation are presented in terms of applied model. The possibility theorem for the problem of getting the solution having some preassigned properties by means of parametric control is proved.
1. Introduction
Since the work [1], the theory of contrasting structures has become one of the most booming areas of research of the singularly perturbed differential equations [2–4].
The contrasting structures having the form of nonstationary fronts for parabolic partial differential equations were studied in [5]. The theory was applied to propagation of magnetic fronts in spiral galaxies [6–8]. Here we consider the nonstationary fronts in the Mikhailov “power-society” model [9–12] and the possibility to control them.
In the most general case the “power-society” model has the form of a Neumann boundary value problem for nonlinear parabolic integrodifferential equation. In the absence of some political mechanisms the model is reduced to singularly perturbed parabolic Neumann boundary value problem:
(1)∂p∂t=ɛ2∂2p∂x2+F(p,x),∂p∂x|x=0=∂p∂x|x=1=0,p(x,0)=p0(x).
Here ɛ≪1 is a small positive parameter. Within the “power-society” model this parameter is small if the hierarchy is long or if society is strong.
The steady-state problem has the form
(2)ɛ2∂2p∂x2+F(p,x)=0,∂p∂x|x=0=∂p∂x|x=1=0.
Let the following conditions hold [2–4].
The function F(p,x) has continuous partial derivatives for 0≤x≤1 and p∈(-∞,+∞).
The degenerate equation F(p,x)=0 has three roots p=φ1(x), p=φ2(x), and p=φ3(x) such that φ1(x)<φ2(x)<φ3(x), 0≤x≤1.
The following inequalities
(3)∂F∂p|p=φ1(x)<0,∂F∂p|p=φ2(x)>0,∂F∂p|p=φ3(x)<0
take place.
The equation
(4)Φ(x)=∫φ1(x)φ3(x)F(p,x)dp=0
has isolated root x=x0 on the interval 0<x<1.
Under these conditions,
if Φ′(x0)<0, then the solution p(x,ɛ) of the problem (2) exists such that
(5)limɛ→0p(x,ɛ)={φ3(x),0<x<x0,φ1(x),x0<x<1,
and it is an asymptotically stable stationary solution of problem (1);
if Φ′(x0)>0, solution p(x,ɛ) of problem (2) exists such that
(6)limɛ→0p(x,ɛ)={φ1(x),0<x<x0,φ3(x),x0<x<1,
and it is an asymptotically stable stationary solution of problem (1).
The solutions that satisfy (5) or (6) are called the step-like contrasting structures or stationary fronts. There are also other stable stationary solutions of the problem (1). In particular, under Conditions 1–3 the existence of two more solutions, one of which is close to φ1(x):
(7)limɛ→0p(x,ɛ)=φ1(x),0<x<x0
and the other one is close to φ3(x):
(8)limɛ→0p(x,ɛ)=φ3(x),0<x<x0,
is guaranteed.
There is an important problem of correspondence between a set of initial functions and a set of steady stationary solutions: given initial function p0(x), what steady-state solution will we have at t→+∞? And there is the inverse problem: if one of the steady states is more desirable than others, which conditions on p0(x) guarantee approach to this desirable steady state?
At last, when studying mathematical models of particular processes there is the following question which arises: if the existing p0(x) does not correspond to the desirably steady state, is it possible to change the right-hand part of (1) so that the solution would evolve to the desirable steady state?
This work is aimed at considering these problems for the “power-society” model which describes the dynamics of the power distribution in a hierarchy.
We base our study on the theory of contrasting structures [2–5], especially on the Butuzov-Nedelko theorem [13]. Some other issues related to nonstationary fronts were studied in [14–16].
2. Nonstationary Fronts and Interpretation in the Nonlinear Singularly Perturbed “Power-Society” Model
This section deals with mathematical modeling of the processes of power dynamics in the hierarchical structures. The model was firstly introduced by Mikhailov, 1994, and the books by Samarskii and Mikhailov 1997 and Mikhailov 2005 should also be mentioned.
Here the hierarchy is a ranked set of instances. Each instance has a particular set of powers. The amount of powers changes with time, and we call such variability the power dynamics. We suppose that there exists a numerical variable which specifies the amount of powers of a particular instance. The power dynamics appear through (a) the self-streamlining of the hierarchy and (b) the influence of the society.
Let us denote the rank of the instance in the hierarchy by x so that x=0 at the top of the hierarchy and x=1 at the bottom. Denote by p(x,t) the amount of powers of instance at time t.
The equation of the “power-society” model [1, 2] has the form (1), and F(p,x) is called the reaction of a civil society. The paper [1] has shown that if F(p,x)=-k1(p-p0(x)) (where k1=const>0 and the function p0(x) is the attractive power profile), then the solution p=p0(x) of the stationary degenerated equation F(p,x)=0 is stable. This means that the solution p(x,t,ɛ) of (1) tends to p0(x) when t→+∞, 0<x<1. So for sufficiently large values of t the power profile is close to p0(x).
It was very important in the paper [1] that only one attractive profile is supposed to exist. Here we consider the case of two stable power profiles φ1(x) and φ3(x), and each of them is attractive. We call φ1(x) the participatory profile and φ3(x) the iron-hand profile. Both of them are stable due to inequalities (3).
Henceforth we consider the function F(p,x) having the cubic nonlinearity. So we consider the equation
(9)∂p∂t=ɛ2∂2p∂x2-k1(x)(p-φ1(x))(p-φ2(x))(p-φ3(x)),
with boundary value conditions
(10)∂p∂x|x=0=∂p∂x|x=1=0.
The following conditions are supposed to hold.
Condition 1.
The functions k1(x), φi(x), i=1,2,3, have the continuous derivatives for 0≤x≤1.
Condition 2.
Inequalities k1(x)>0, φ1(x)<φ2(x)<φ3(x) hold true for 0≤x≤1.
We also notice that though φ1(x)>0 because of the politological meaning of the function φi(x), this condition must not be required from the mathematical point of view.
Making the substitution (see [17])
(11)p=qφ3(x)-φ1(x)2+φ3(x)+φ1(x)2,
we obtain the equation for the function u(x,t,ɛ):
(12)∂q∂t=ɛ2∂2q∂x2-γ(x)(q2-1)(q-φ(x))+ɛ2∂q∂xΦ1(x)+ɛ2Φ2(q,x).
Here
(13)γ(x)=k1(x)(φ3(x)-φ1(x)2)2>0,φ(x)=2φ2(x)-φ3(x)-φ1(x)φ3(x)-φ1(x)
(notice that -1<φ(x)<1),
(14)Φ1(x)=2(φ3′(x)-φ1′(x))φ3(x)-φ1(x),Φ2(q,x)=qφ3″(x)-φ1″(x)φ3(x)-φ1(x)+φ3″(x)+φ1″(x)φ3(x)-φ1(x).
Function u(x,t,ɛ) satisfies boundary conditions
(15)[φ3(0)-φ1(0)]∂q∂x|x=0+[φ3′(0)-φ1′(0)]q|x=0+φ3′(0)+φ1′(0)=0,[φ3(1)-φ1(1)]∂q∂x|x=1+[φ3′(1)-φ1′(1)]q|x=1+φ3′(1)+φ1′(1)=0.
Consider stationary (∂/∂t=0) equation related to (12):
(16)ɛ2q′′=γ(x)(q2-1)(q-φ(x))-ɛ2∂q∂xΦ1(x)-ɛ2Φ2(q,x).
Using the boundary functions method [18] we construct the asymptotic contrast solution of the problem (16) and (15). The first-order asymptotic expansion has the form
(17)q(x,t,ɛ)=q¯0(x)+ɛq¯1(x)+Π0q(τ)+ɛΠ1q(τ)+ɛQq(τ0,τ1),
where q¯0(x) and q¯1(x) are the regular terms of asymptotic expansion, Π0q(τ) and Π1q(τ) are zero- and first-order transition layer functions, τ=(x-x*)/ɛ is a stretched variable, x*=x*(ɛ) is a transition point in a small vicinity of which the transition layer is localized, and Qq(τ0,τ1) is function describing the boundary layers near the points x=0, x=1 and τ0=x/ɛ, τ1=(1-x)/ɛ. The transition point has the following asymptotic form:
(18)x*=x0+ɛx1+ɛ2x2+⋯.
Using the boundary functions method procedure [3, 4] we obtain that the principal term x0 of the expansion (18) can be found from the equation
(19)φ(x0)=0.
The full principal order function q~(τ)=q¯0(x0)+Π0q(τ) can be found from equation
(20)q~′′=γ(x0)q~(q~2-1),(21)q~(-∞)=1,q~(+∞)=-1.
From (20) and (21) we have
(22)q~=1-exp[2γ(x0)τ]1+exp[2γ(x0)τ].
So the principal term of the stationary power profile has the form
(23)pst(x,ɛ)=1-exp[2γ(x0)(x-x0)/ɛ]1+exp[2γ(x0)(x-x0)/ɛ]×φ3(x)-φ1(x)2+φ3(x)+φ1(x)2.
The power profile pst(x,ɛ) is close to the iron-hand profile φ3(x) when 0≤x<x0 and to the participatory profile φ1(x) when x0<x≤1. In the vicinity of the transition point x0 we have ∂p/∂x≅ɛ-1. We call such power profiles the contrast power profiles.
Equation (19) can be written in the form
(24)φ2(x0)=φ3(x0)+φ1(x0)2.
We call the function h1(x)=φ2(x)-φ1(x) the participatory domain's width and function h3(x)=φ3(x)-φ2(x) the iron-hand domain's width. Then (20) can be interpreted in the following way: at the transition point x0 of the stationary contrast power profile (SCPP) the participatory domain's width is equal to the iron-hand domain's width h1(x0)=h3(x0).
The stability of contrast structures of (9) was investigated by Bozhevol’nov and Nefëdov [5] and Vasil’eva et al. [6]. In terms of the “power-society” model the stability result can be interpreted as follows.
SCPP, which are close to the iron-hand profile at the top ranks of the hierarchy (p≈φ3(x) when 0≤x<x0) and to the participatory profile at the bottom ranks (p≈φ1(x) when x0<x≤1), are stable if the iron-hand domain's width is greater than the participatory domain's width at the top ranks of the hierarchy (h3(x)>h1(x) when 0≤x<x0) and less at the bottom ranks (h3(x)<h1(x) when x0<x≤1). If h3(x)<h1(x) when 0≤x<x0 and h3(x)>h1(x) when x0<x≤1 then the SCPP is unstable.
Remark 1.
Similar statement holds for the so-called confederative SCPP which are close to the participatory profile at the top ranks of the hierarchy and to the iron-hand profile at the bottom ranks (p≈φ1(x) when 0≤x<x0 and p≈φ3(x) when x0<x≤1). They are stable if h3(x)<h1(x) when 0≤x<x0 and h3(x)>h1(x) when x0<x≤1 and unstable if h3(x)>h1(x) when 0≤x<x0 and h3(x)<h1(x) when x0<x≤1.
Consider again nonstationary equation (9). Suppose that at time t=t0 contrasting structure has appeared with the transition layer at the vicinity of the point x=ξ. Then for t>t0 the solution is a nonstationary contrast structure: p≈φ3(x) when x<R(t,ɛ) and p≈φ1(x) when x>R(t,ɛ), where the transition point R(t,ɛ) depends on time. We call such power profile the nonstationary contrast power profile (NCPP).
Let us construct the asymptotic NCPP.
Like in Section 3, make the substitution (6) and consider (7). It was shown in [7] that the principal term of the nonstationary contrast structure looks similar to one of the stationary contrast structure (18):
(25)q^=1-exp[2γ(R)(x-R)/ɛ]1+exp[2γ(R)(x-R)/ɛ],
where the function R=R(t,ɛ) can be found from the equation
(26)dRdt=-ɛ2γ(R)φ(R).
So the principal term of NCPP has the form
(27)p(x,t,ɛ)=1-exp[2γ(R)(x-R)/ɛ]1+exp[2γ(R)(x-R)/ɛ]×φ3(x)-φ1(x)2+φ3(x)+φ1(x)2.
Power profile p(x,t,ɛ) is close to the iron-hand profile φ3(x) when 0≤x<R and to the participatory profile φ1(x) when R<x≤1. The value of dR/dt represents the speed of the transition layer. In terms of the “power-society” model the expression for dR/dt has the form
(28)dRdt=-ɛk1(R)2[2φ2(R)-φ3(R)-φ1(R)]
or
(29)dRdt=ɛk1(R)2[h3(R)-h1(R)].
Consider now some important cases of using formula (29); see also [17].
2.1. Attraction to the “Iron-Hand” Profile (1)
Let the “iron-hand” domain’s width be larger than participatory domain's width: h3(x)>h1(x) for any x∈[0,1]. This means that the iron-hand profile looks more attractive from the society's point of view. Then SCPP do not exist because (24) has no roots. After appearing at time t0 the contrast structure begins to move according to formula (29). Evidently dR/dt>0, and after small time of order ɛ-1 transition point R(t,ɛ) comes to the right end of the segment [0,1]. So the power profile appears close to the iron-hand profile for any x∈[0,1]. Notice that if at time t=0 function p(x,0,ɛ) is entirely in the participatory domain then for any t the power profile is close to the participatory profile even if h3(x)>h1(x). For appearing the power profile close to the iron-hand profile function p(x,0,ɛ) must be located in the iron-hand domain on at least one point in the interval (0,1). This statement is based on the theorem proved by Bozhevol’nov and Nefedov [5].
2.2. Attraction to the “Iron-Hand” Profile (2)
Let point x0∈(0,1) exist such that h1(x0)=h3(x0), h1(x)>h3(x) when x<x0 and h1(x)<h3(x) when x>x0. Then unstable SCPP exist having transition layer in the vicinity of the point x0. Let function p(x,0,ɛ) be in the iron-hand domain for x<ξ and in the participatory domain for x>ξ where x0<ξ<1. Then the power profile p(x,t,ɛ) is attracted to the iron-hand profile for x<ξ and to the participatory profile for x>ξ. After appearing at time t0 the contrast structure begins to move according to formula (29) and initial condition R(0,ɛ)=ξ. As ξ>x0 then dR/dt>0. So after small time of order ɛ-1 transition point R(t,ɛ) comes to the right end of the segment [0,1]. So the power profile appears close to the iron-hand profile for any x∈[0,1].
2.3. Attraction to the Participatory Profile (1)
Let h3(x)<h1(x) for any x∈[0,1]. This means that the participatory profile looks more attractive from the society's point of view. Then SCPP do not exist because (24) has no roots. After appearing at time t0 the contrast structure begins to move according to formula (29). Evidently dR/dt<0 and after small time of order ɛ-1 transition point R(t,ɛ) comes to the left end of the segment [0,1]. So the power profile appears close to the participatory profile for any x∈[0,1]. Notice that if at time t=0 function p(x,0,ɛ) is entirely in the iron-hand domain then for any t the power profile is close to the iron-hand profile even if h3(x)>h1(x). For appearing the power profile close to the participatory profile function p(x,0,ɛ) must be smooth and located in the participatory domain on at least one point in the interval (0,1). This statement is based on the theorem proved by Bozhevol’nov and Nefedov [5].
2.4. Attraction to the Participatory Profile (2)
Let point x0∈(0,1) exist such that h1(x0)=h3(x0), h1(x)>h3(x) when x<x0 and h1(x)<h3(x) when x>x0. Then unstable SCPP exist having transition layer in the vicinity of the point x0. Let function p(x,0,ɛ) be in the iron-hand domain for x<ξ and in the participatory domain for x>ξ where x0<ξ<1. Then according to (3) power profile p(x,t,ɛ) is attracted to the iron-hand profile for x<ξ and to the participatory profile for x>ξ. After appearing at time t0 the contrast structure begins to move according to formula (26) and initial condition R(0,ɛ)=ξ. As ξ<x0 then dR/dt<0. So after small time of order ɛ-1 transition point R(t,ɛ) comes to the left end of the segment [0,1]. So the power profile appears close to the participatory profile.
2.5. Attraction to SCPP
Let point x0∈(0,1) exist such that h1(x0)=h3(x0), h1(x)<h3(x) when x<x0 and h1(x)>h3(x) if x>x0. Then the stable SCPP exist having transition layer in the vicinity of the point x0. If there is ξ∈(0,1) such that the initial function p0(x) satisfies p0(x)>φ2(x), 0<x<ξ, and p0(x)<φ2(x), ξ<x<1, then the contrasting structure appears after a short time t0 such that p(x)≈φ3(x) for 0<x<ξ and p(x)≈φ1(x) for ξ<x<1. Then the contrasting structure begins to move according to formula (29) and initial condition R(t0,ɛ)=ξ. So if ξ<x0 then dR/dt>0 and if ξ>x0 then dR/dt<0. So when t→∞, the transition point R(t,ɛ) tends to the stationary transition point x0. Thus the SCPP having transition layer in the vicinity of point x0 appears.
3. Parametric Optimization
The total amount of power of the hierarchy is P~(t,ɛ)=∫01p(x,t,ɛ). It was shown in [19] that there exists the optimal value P0 of the total power which provides a maximum of steady-state consumption per capita (in frame of the “power-society-economics” model [19]). So we should introduce the control parameter into the “power-society” model to make it controllable. So the problem would be to find the value of the control parameter under which P~(t,ɛ)→P0, when t→∞, ɛ→0.
Generally speaking, the model could be formulated such that the control is considered to be a function of time or x. In any case, the control describes the exogenous impact on the political system, such as a political pressure through media and political institutions. We restrict ourselves to the parametric control.
Definition 2.
The value P0 is called the asymptotically achievable amount of the total power if there exists an admissible value of control parameter u such that the steady-state total power Pu(ɛ) satisfies Pu(ɛ)→P0 when t→∞, ɛ→0.
So consider the “power-society” model with nonlinear reaction of civil society:
(30)f(p,x,t,u)=-k1(x)(p-φ1(x))×(p-(φ2(x)+γu))(p-φ3(x)).
Here k1(x)>0, the functions k1(x), φ1(x), φ2(x), and φ3(x) have continuous derivatives, and γ is a constant. Thus the lowest and the biggest roots φ1(x), φ3(x) of the degenerate equation F(p,x,u)=0 do not depend on the control, but there is an impact from the control to the “middle” root φ2(x). So the model has the form
(31)∂p∂t=ɛ2∂2p∂x2-k1(x)(p-φ1(x))×(p-(φ2(x)+γu))(p-φ3(x)),(32)∂p∂x|x=0=∂p∂x|x=1=0,p|t=0=p0(x).
The initial function p0(x) is supposed to be smooth and satisfying φ1(x)<p0(x)<φ3(x). In other words, the initial distribution of power is between the iron-hand and participatory profiles.
The steady-state equation for (31) has the form
(33)ɛ2∂2p∂x2=k1(x)(p-φ1(x))(p-(φ2(x)+γu))(p-φ3(x)),∂p∂x|x=0=∂p∂x|x=1=0.
Let us stress here that the control influences the relation between the width of the iron-hand domain and the width of the participatory domain.
Let the following conditions be fulfilled.
Condition 3.
We have -1≤u≤1.
Condition 4.
We have φ1(x)<φ2(x)-γ<φ2(x)+γ<φ3(x) for 0≤x≤1.
Condition 5.
We have H′(x)<0 for 0≤x≤1, and here H(x)=φ1(x)+φ3(x)-2φ2(x).
Conditions 3 and 4 introduce the normalization of the control such that for any admissible control u∈[-1;1] the root φ2(x)+γu is between the φ1(x) and φ3(x). The steady-state solution has no more than one transition point due to Condition 5.
If the control parameter u is increased, the root of the equation φ1(x)+φ3(x)-2(φ2(x)+γu)=0 will move to the left. This means greater support to the participation ideas. If it exists. Analogically, the less the value of u is, the more to the right the root of this equation is.
Consider the following problem. Let the desirable (optimal) value of total power be P0. Is there a value of control parameter u, under which the steady-state solution (33) is such that Pu(ɛ)=∫01p(x,ɛ)dx→P0 when ɛ→0?
From the practical point of view, such a formulation of the problem can be justified in the following way. We know from the “power-society-economics” model that the optimal value of the total power is some P0, so we should try to tune the political system to provide this optimal value of power for the steady-state regime.
Several cases should be distinguished.
Let us start our consideration from the situation in which both equations
(34)φ1(x)+φ3(x)-2(φ2(x)+γ)=0,(35)φ1(x)+φ3(x)-2(φ2(x)-γ)=0
have roots in the interval (0;1). Let us denote these roots by a and b, respectively. Here we have a<b in view of Condition 5.
The points x=a and x=b are the main asymptotic terms for the boundaries of the range within which the transition point of the stationary front is located.
Therefore, the value
(36)Pas(u)=limɛ→0P(u,ɛ)=limɛ→0∫01p(x,u,ɛ)dx
satisfies the inequality
(37)∫0aφ3(x)dx+∫a1φ1(x)dx≤Pas(u)≤∫0bφ3(x)dx+∫b1φ1(x)dx.
Thus, in this case, P0 is asymptotically achievable, if inequality
(38)∫0aφ3(x)dx+∫a1φ1(x)dx≤P0≤∫0bφ3(x)dx+∫b1φ1(x)dx
holds true.
The steady-state problem (32) and (33) has also solutions without transition layers: the iron-hand profile and the participatory one. So the values of total power
(39)P0=∫01φ3(x)dx,P0=∫01φ1(x)dx,
are also asymptotically achievable.
So, if both (34) and (35) have roots in the interval (0;1) then the set of asymptotically achievable values comprises the closed interval (38) and two isolated values (39): one of them is to the left of this closed interval, and the other one is to the right of it.
Now let us consider the situation in which (34) has a root x=a∈(0;1) and (35) has no roots on (0;1).
In other words, φ2(x)-γ<(φ1(x)+φ3(x))/2 for any x. So the values of P0 in the closed interval
(40)∫0aφ3(x)dx+∫a1φ1(x)dx≤P0≤∫01φ3(x)dx
are asymptotically achievable.
Let now (35) have the root x=b∈(0;1), and let (34) have no roots in the interval. That is, at u=1, for any x, the middle root φ2(x)+γ is larger than half-sum of φ1(x) and φ3(x). Then the asymptotically achievable values are given by the inequality
(41)∫01φ1(x)dx≤P0≤∫0bφ3(x)dx+∫b1φ1(x)dx.
At last, consider the case in which neither of the equations has a root in the interval (0;1). Then the following three subcases are possible:
(42)φ1(x)+φ3(x)-2(φ2(x)+γ)<0,φ1(x)+φ3(x)-2(φ2(x)-γ)<0foranyx,(43)φ1(x)+φ3(x)-2(φ2(x)+γ)>0,φ1(x)+φ3(x)-2(φ2(x)-γ)>0foranyx,(44)φ1(x)+φ3(x)-2(φ2(x)+γ)>0,φ1(x)+φ3(x)-2(φ2(x)-γ)<0foranyx.
It can be easily shown (see [2], e.g.) that in subcases (42) and (43) a steady-state front does not exist for any control parameter. So only the values of P0 given by (39) are asymptotically achievable.
In the subcase (44) for any given x0∈(0;1), such u exists that the problem (32) and (33) has the stationary front with the transition point in the ɛ-vicinity of x0. Therefore, any P0 from closed interval
(45)∫01φ1(x)dx≤P0≤∫01φ3(x)dx
is asymptotically achievable.
The above speculations can be summarized as follows.
Theorem 3.
Consider the problem (32) and (33) with parametric control. Let the following conditions be satisfied:
k1(x), φ1(x), φ2(x), φ3(x)∈C1[0,1],
-1≤u≤1,
φ1(x)<φ2(x)-γ<φ2(x)+γ<φ3(x),
H′(x)<0, where H(x)=φ1(x)+φ3(x)-2φ2(x).
Then the set of asymptotically achievable values is not empty.
After the provided analysis of a steady-state problem (32) and (33), we go back to the initial parabolic partial problem (31) and (32).
Let some value P0 be asymptotically achievable in the corresponding stationary problem. It means that there is a value of parametric control u, at which the problem (31), (32) has the steady-state solution for which the total power Pu(ɛ) of the hierarchy asymptotically tends to P0 when ɛ→0. However for the P~u(t,ɛ)=∫01pu(x,t,ɛ)dx we have P~u(t,ɛ)→Pu(ɛ) just for some class of initial functions p0(x).
Thus, there is a problem to determine the class of initial power distributions for which, under the found value parametric control, the solution of the parabolic partial problem converges to the proper steady-state solution at t→∞.
The answer is given by the following theorem.
Theorem 4.
(1) Let all the conditions of Theorem 3 be satisfied;
(2) let P0 be an asymptotically achievable value of total power (denote by u=u0 the corresponding value of parametric control);
(3) let the point x0∈(0;1) exist such that
(46)P0=∫0x0φ3(x)dx+∫x01φ1(x)dx;
(4) let the points x1∈(0;x0), x2∈(x0;1) exist such that
(47)p0(x1)>φ2(x1)+γu0,p0(x2)<φ2(x2)+γu0.
Then the solution pu0(x,t,ɛ) of the parabolic partial problem (31) and (32) is such that the total power
(48)Qu0(t,ɛ)=∫01pu0(x,t,ɛ)dx
converges to
(49)P0:limɛ→∞limt→∞P(t,u0,ɛ)=P0.
Proof.
It is easy to see that under these conditions the Butuzov-Nedelko theorem is fulfilled [13]. Therefore, at the chosen value of control u=u0, the solution pu0(x,t,ɛ) of the problem (31) and (32) has a passage to the limit
(50)limt→∞pu0(x,t,ɛ)=pst,u0(x,ɛ),
where pst(x,ɛ) is the steady state solution for which
(51)limt→∞pst,u0(x,ɛ)={φ3(x),x<x0,φ1(x),x>x0.
By integrating (50) from x=0 to x=1 we get
(52)limt→∞∫01pu0(x,t,ɛ)dx=∫01pst,u0(x,ɛ)dx.
Passing to a limit ɛ→0, we get
(53)limɛ→0limt→∞∫01pu0(x,t,ɛ)dx=limɛ→0∫01pst,u0(x,ɛ)dx=∫0x0φ3(x)dx+∫x01φ1(x)dx=P0.
Thus limɛ→0limt→∞P(t,u0,ɛ)=P0. Theorem 4 is proved.
4. Conclusion
It is shown that the theory of contrasting structures in singularly perturbed boundary value problems allows for investigating the properties of nonstationary fronts in the singularly perturbed “power-society” model. Depending on the initial condition, these fronts evolve to one of the asymptotically stable steady-state distributions of power within a government hierarchy.
There are some reasons to introduce a concept of desirable steady-state total amount of power of the hierarchy. The possibility theorem is proved for the problem of getting this amount by means of parametric control. The results can be used in investigating governing hierarchical systems.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
Acknowledgment
This paper was partially supported by Russian Foundation for Basic Research (project 12-06-00205).
ButuzovV. F.Vasil’evaA. B.Asymptotic behavior of a solution of contrasting structure type1987426831841MR934815ButuzovV. F.Vasil’evaA. B.NefedovN. N.Asymptotic theory of contrasting structures. A survey199758710681091Vasil’evaA. B.ButuzovV. F.NefedovN. N.Singularly perturbed problems with boundary and internal layers20102681258273Vasil’evaA. B.ButuzovV. F.NefëdovN. N.Contrast structures in singularly perturbed problems199843799851MR1800038Bozhevol’novYu. V.NefëdovN. N.Front motion in a parabolic reaction-diffusion problem201050226427310.1134/S0965542510020089MR2681153Vasil’evaA.NikitinA.PetrovA.Stability of contrasting solutions of nonlinear hydromagnetic dynamo equations and magnetic fields reversals in galaxies1994781–426127910.1080/03091929408226582MR1335773MossD.PetrovA.SokoloffD.The motion of magnetic fronts in spiral galaxies2000921-212914910.1080/03091920008203714MR1792617Vasil’evaA. B.PetrovA. P.PlotnikovÀ. À.Alternating contrast structures199838914711480MikhailovA. P.Mathematical modeling of the authority in the hierarchical structures19946article 6SamarskiiA. A.MikhailovA. P.20023London, UKTaylor & FrancisNumerical InsightsMR1920084MikhailovA. P.Mathematical modeling of power distribution in state ierarchical structures interacting with civil society2Proceedings of 14th IMACS World Congress1994Atlanta, Ga, USA828MikhailovA. P.The research of the “Power-Society” systemFizmatlit, 2006, (in Russian)ButuzovV. F.NedelkoI. V.On the global influence domain of stable solutions with internal layers20011925135210.1070/SM2001v192n05ABEH000563MR1859701NefëdovN. N.NikitinA. G.The method of differential inequalities for step-type contrast structures in singularly perturbed integro-differential equations in the spatially two-dimensional case200642573974910.1134/S0012266106050132MR2292167NefedovN. N.NikitinA. G.PetrovaM. A.RekkeL.Moving fronts in integro-parabolic reaction-advection-diffusion equations20114791318133210.1134/S0012266111090096MR2918823ZBL1269.45008NefëdovN. N.NikitinA. G.Boundary and internal layers in the reaction-diffusion problem with a nonlocal inhibitor20115161011101910.1134/S0965542511060157MR2859176ZBL1249.35159DmitrievM.PetrovA.ZhukovaG.The nonlinear “Authority-Society” model200321-23136Vasil’evaA. B.ButuzovV. F.The method of boundary layer functions1985211011071112MR814556DmitrievM. G.PavlovA. A.PetrovA. P.Optimal power volume in social-economical hierarchy by consumption per capita criterion20074411