Abstract and Applied Analysis

Volume 2013 (2013), Article ID 172654, 8 pages

http://dx.doi.org/10.1155/2013/172654

## Nonstationary Fronts in the Singularly Perturbed Power-Society Model

^{1}Institute of System Analysis of RAS, Higher School of Economics, Moscow 101000, Russia^{2}Russian State Social University, Moscow 129226, Russia^{3}Keldysh Institute of Applied Mathematics RAS, Moscow 125047, Russia

Received 23 July 2013; Accepted 15 September 2013

Academic Editor: Allaberen Ashyralyev

Copyright © 2013 M. G. Dmitriev et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

#### Abstract

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: Here 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

Let the following conditions hold [2–4].(1) The function has continuous partial derivatives for and .(2) The degenerate equation has three roots , , and such that , .(3) The following inequalities take place.(4) The equation has isolated root on the interval .

Under these conditions,(i) if , then the solution of the problem (2) exists such that and it is an asymptotically stable stationary solution of problem (1);(ii) if , solution of problem (2) exists such that 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 : and the other one is close to : 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 , what steady-state solution will we have at ? And there is the inverse problem: if one of the steady states is more desirable than others, which conditions on 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 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 so that at the top of the hierarchy and at the bottom. Denote by the amount of powers of instance at time .

The equation of the “power-society” model [1, 2] has the form (1), and is called the reaction of a civil society. The paper [1] has shown that if (where and the function is the attractive power profile), then the solution of the stationary degenerated equation is stable. This means that the solution of (1) tends to when , . So for sufficiently large values of the power profile is close to .

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 and , and each of them is attractive. We call the participatory profile and the iron-hand profile. Both of them are stable due to inequalities (3).

Henceforth we consider the function having the cubic nonlinearity. So we consider the equation with boundary value conditions The following conditions are supposed to hold.

*Condition 1. *The functions , , , have the continuous derivatives for .

*Condition 2. *Inequalities , hold true for .

We also notice that though because of the politological meaning of the function , this condition must not be required from the mathematical point of view.

Making the substitution (see [17]) we obtain the equation for the function : Here (notice that ), Function satisfies boundary conditions Consider stationary () equation related to (12): 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 where and are the regular terms of asymptotic expansion, and are zero- and first-order transition layer functions, is a stretched variable, is a transition point in a small vicinity of which the transition layer is localized, and is function describing the boundary layers near the points , and , . The transition point has the following asymptotic form: Using the boundary functions method procedure [3, 4] we obtain that the principal term of the expansion (18) can be found from the equation The full principal order function can be found from equation From (20) and (21) we have So the principal term of the stationary power profile has the form The power profile is close to the iron-hand profile when and to the participatory profile when . In the vicinity of the transition point we have . We call such power profiles the contrast power profiles.

Equation (19) can be written in the form
We call the function the participatory domain's width and function the iron-hand domain's width. Then (20) can be interpreted in the following way: *at the transition point ** of the stationary contrast power profile (SCPP) the participatory domain's width is equal to the iron-hand domain's width *.

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 ( when ) and to the participatory profile at the bottom ranks ( when ), are stable if the iron-hand domain's width is greater than the participatory domain's width at the top ranks of the hierarchy ( when ) and less at the bottom ranks ( when ). If when and when 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 ( when and when ). They are stable if when and when and unstable if when and when .

Consider again nonstationary equation (9). Suppose that at time contrasting structure has appeared with the transition layer at the vicinity of the point . Then for the solution is a nonstationary contrast structure: when and when , where the transition point 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): where the function can be found from the equation So the principal term of NCPP has the form Power profile is close to the iron-hand profile when and to the participatory profile when . The value of represents the speed of the transition layer. In terms of the “power-society” model the expression for has the form or 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: for any . 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 the contrast structure begins to move according to formula (29). Evidently , and after small time of order transition point comes to the right end of the segment . So the power profile appears close to the iron-hand profile for any . Notice that if at time function is entirely in the participatory domain then for any the power profile is close to the participatory profile even if . For appearing the power profile close to the iron-hand profile function must be located in the iron-hand domain on at least one point in the interval . 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 exist such that , when and when . Then unstable SCPP exist having transition layer in the vicinity of the point . Let function be in the iron-hand domain for and in the participatory domain for where . Then the power profile is attracted to the iron-hand profile for and to the participatory profile for . After appearing at time the contrast structure begins to move according to formula (29) and initial condition . As then . So after small time of order transition point comes to the right end of the segment . So the power profile appears close to the iron-hand profile for any .

##### 2.3. Attraction to the Participatory Profile (1)

Let for any . 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 the contrast structure begins to move according to formula (29). Evidently and after small time of order transition point comes to the left end of the segment . So the power profile appears close to the participatory profile for any . Notice that if at time function is entirely in the iron-hand domain then for any the power profile is close to the iron-hand profile even if . For appearing the power profile close to the participatory profile function must be smooth and located in the participatory domain on at least one point in the interval . This statement is based on the theorem proved by Bozhevol’nov and Nefedov [5].

##### 2.4. Attraction to the Participatory Profile (2)

Let point exist such that , when and when . Then unstable SCPP exist having transition layer in the vicinity of the point . Let function be in the iron-hand domain for and in the participatory domain for where . Then according to (3) power profile is attracted to the iron-hand profile for and to the participatory profile for . After appearing at time the contrast structure begins to move according to formula (26) and initial condition . As then . So after small time of order transition point comes to the left end of the segment . So the power profile appears close to the participatory profile.

##### 2.5. Attraction to SCPP

Let point exist such that , when and if . Then the stable SCPP exist having transition layer in the vicinity of the point . If there is such that the initial function satisfies , , and , , then the contrasting structure appears after a short time such that for and for . Then the contrasting structure begins to move according to formula (29) and initial condition . So if then and if then . So when , the transition point tends to the stationary transition point . Thus the SCPP having transition layer in the vicinity of point appears.

#### 3. Parametric Optimization

The total amount of power of the hierarchy is . It was shown in [19] that there exists the optimal value 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 , when , .

Generally speaking, the model could be formulated such that the control is considered to be a function of time or . 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 is called the asymptotically achievable amount of the total power if there exists an admissible value of control parameter such that the steady-state total power satisfies when , .

So consider the “power-society” model with nonlinear reaction of civil society: Here , the functions , , , and have continuous derivatives, and is a constant. Thus the lowest and the biggest roots , of the degenerate equation do not depend on the control, but there is an impact from the control to the “middle” root . So the model has the form The initial function is supposed to be smooth and satisfying . 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

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 .

*Condition 4. *We have for .

*Condition 5. *We have for , and here .

Conditions 3 and 4 introduce the normalization of the control such that for any admissible control the root is between the and . The steady-state solution has no more than one transition point due to Condition 5.

If the control parameter is increased, the root of the equation will move to the left. This means greater support to the participation ideas. If it exists. Analogically, the less the value of 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 . Is there a value of control parameter , under which the steady-state solution (33) is such that when ?

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 , 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 have roots in the interval . Let us denote these roots by and , respectively. Here we have in view of Condition 5.

The points and 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 satisfies the inequality

Thus, in this case, is asymptotically achievable, if inequality 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 are also asymptotically achievable.

So, if both (34) and (35) have roots in the interval 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 and (35) has no roots on .

In other words, for any . So the values of in the closed interval are asymptotically achievable.

Let now (35) have the root , and let (34) have no roots in the interval. That is, at , for any , the middle root is larger than half-sum of and . Then the asymptotically achievable values are given by the inequality At last, consider the case in which neither of the equations has a root in the interval . Then the following three subcases are possible:

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 given by (39) are asymptotically achievable.

In the subcase (44) for any given , such exists that the problem (32) and (33) has the stationary front with the transition point in the -vicinity of . Therefore, any from closed interval 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: *(1)*, , , ,*(2)*,*(3)*, *(4)*, where . ** 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 be asymptotically achievable in the corresponding stationary problem. It means that there is a value of parametric control , at which the problem (31), (32) has the steady-state solution for which the total power of the hierarchy asymptotically tends to when . However for the we have just for some class of initial functions .

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 .

The answer is given by the following theorem.

Theorem 4. *(1) Let all the conditions of Theorem 3 be satisfied;**(2) let be an asymptotically achievable value of total power (denote by the corresponding value of parametric control);**(3) let the point exist such that
**(4) let the points , exist such that
**
Then the solution of the parabolic partial problem (31) and (32) is such that the total power
**
converges to
*

*Proof. *It is easy to see that under these conditions the Butuzov-Nedelko theorem is fulfilled [13]. Therefore, at the chosen value of control , the solution of the problem (31) and (32) has a passage to the limit
where is the steady state solution for which
By integrating (50) from to we get
Passing to a limit , we get
Thus . 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).

#### References

- V. F. Butuzov and A. B. Vasil’eva, “Asymptotic behavior of a solution of contrasting structure type,”
*Akademiya Nauk SSSR*, vol. 42, no. 6, pp. 831–841, 1987. View at Google Scholar · View at MathSciNet - V. F. Butuzov, A. B. Vasil’eva, and N. N. Nefedov, “Asymptotic theory of contrasting structures. A survey,”
*Automation and Remote Control*, vol. 58, no. 7, pp. 1068–1091, 1997. View at Google Scholar - A. B. Vasil’eva, V. F. Butuzov, and N. N. Nefedov, “Singularly perturbed problems with boundary and internal layers,”
*Proceedings of the Steklov Institute of Mathematics*, vol. 268, no. 1, pp. 258–273, 2010. View at Google Scholar - A. B. Vasil’eva, V. F. Butuzov, and N. N. Nefëdov, “Contrast structures in singularly perturbed problems,”
*Fundamental'naya i Prikladnaya Matematika*, vol. 4, no. 3, pp. 799–851, 1998. View at Google Scholar · View at MathSciNet - Yu. V. Bozhevol’nov and N. N. Nefëdov, “Front motion in a parabolic reaction-diffusion problem,”
*Computational Mathematics and Mathematical Physics*, vol. 50, no. 2, pp. 264–273, 2010. View at Publisher · View at Google Scholar · View at MathSciNet - A. Vasil’eva, A. Nikitin, and A. Petrov, “Stability of contrasting solutions of nonlinear hydromagnetic dynamo equations and magnetic fields reversals in galaxies,”
*Geophysical and Astrophysical Fluid Dynamics*, vol. 78, no. 1–4, pp. 261–279, 1994. View at Publisher · View at Google Scholar · View at MathSciNet - D. Moss, A. Petrov, and D. Sokoloff, “The motion of magnetic fronts in spiral galaxies,”
*Geophysical and Astrophysical Fluid Dynamics*, vol. 92, no. 1-2, pp. 129–149, 2000. View at Publisher · View at Google Scholar · View at MathSciNet - A. B. Vasil’eva, A. P. Petrov, and À. À. Plotnikov, “Alternating contrast structures,”
*Computational Mathematics and Mathematical Physics*, vol. 38, no. 9, pp. 1471–1480, 1998. View at Google Scholar - A. P. Mikhailov, “Mathematical modeling of the authority in the hierarchical structures,”
*Mathematical Modeling*, vol. 6, article 6, 1994 (Russian). View at Google Scholar - A. A. Samarskii and A. P. Mikhailov,
*Principles of mathematical modeling, Ideas, methods, examples*, vol. 3 of*Numerical Insights*, Taylor & Francis, London, UK, 2002. View at MathSciNet - A. P. Mikhailov, “Mathematical modeling of power distribution in state ierarchical structures interacting with civil society,” in
*Proceedings of 14th IMACS World Congress*, vol. 2, p. 828, Atlanta, Ga, USA, 1994. - A. P. Mikhailov, “The research of the “Power-Society” system,” Fizmatlit, 2006, (in Russian).
- V. F. Butuzov and I. V. Nedelko, “On the global influence domain of stable solutions with internal layers,”
*Differential Equations*, vol. 192, no. 5, pp. 13–52, 2001. View at Publisher · View at Google Scholar · View at MathSciNet - N. N. Nefëdov and A. G. Nikitin, “The method of differential inequalities for step-type contrast structures in singularly perturbed integro-differential equations in the spatially two-dimensional case,”
*Differential Equations*, vol. 42, no. 5, pp. 739–749, 2006. View at Publisher · View at Google Scholar · View at MathSciNet - N. N. Nefedov, A. G. Nikitin, M. A. Petrova, and L. Rekke, “Moving fronts in integro-parabolic reaction-advection-diffusion equations,”
*Differential Equations*, vol. 47, no. 9, pp. 1318–1332, 2011. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - N. N. Nefëdov and A. G. Nikitin, “Boundary and internal layers in the reaction-diffusion problem with a nonlocal inhibitor,”
*Computational Mathematics and Mathematical Physics*, vol. 51, no. 6, pp. 1011–1019, 2011. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - M. Dmitriev, A. Petrov, and G. Zhukova, “The nonlinear “Authority-Society” model,”
*WSEAS Transactions on Mathematics*, vol. 2, no. 1-2, pp. 31–36, 2003. View at Google Scholar - A. B. Vasil’eva and V. F. Butuzov, “The method of boundary layer functions,”
*Differential Equations*, vol. 21, no. 10, pp. 1107–1112, 1985. View at Google Scholar · View at MathSciNet - M. G. Dmitriev, A. A. Pavlov, and A. P. Petrov, “Optimal power volume in social-economical hierarchy by consumption per capita criterion,”
*Informacionnye Technologii i Vychislitelnye Sistemy*, vol. 4, pp. 4–11, 2007 (Russian). View at Google Scholar