1. Introduction
Mathematical justifications for modeling digital simulators of automation equipment and mechanisms for regulating the stability of engineering systems are becoming more important due to their increasing complexity. The theoretical foundations of this field were laid out by R.A. DeVore, G.G. Lorentz, M. von Golitschek, Y. Makovoz, G. Mastroianni, G.V. Milovanovi, J. von Neumann, and A.H. Taub [
1,
2,
3,
4], to name a few. At the turn of the 21st century, the need to ensure the goal of 3-E (energy, environmental, and economic efficiency) for innovative transformations grew. The mathematical modeling methods developed by S.V. Alyukov et al. [
5,
6,
7,
8,
9,
10,
11] make it possible to predict the influence of individual factors on the operation of an enterprise as a whole. V. Balakrishnan et al., A. Alabugin, and A. Nikitin showed that new control algorithms with increased accuracy significantly reduce the impact of forecasting errors [
12,
13,
14]. In the economic ecosystems of a number of developing countries, researchers have identified the insufficient decarbonization of industrial territories. This reduces the opportunities for sustainable development, increasing the imbalance in the goals of ensuring 3-E efficiency and innovation. These issues are not being solved because the industrial, urban, and agricultural territories (including households) of a number of developing countries do not meet the requirements on sustainable development: municipalities often use low-tech landfills for industrial and household waste, while available agricultural lands suitable for high-tech production of high-quality organic products constantly decrease.
Structural changes are needed to foster constant interaction between enterprises and institutions of education and science as part of a diversified EPS. This form of economy is effective with the wide-spread and constant use of high-tech “green” energy resources, which is characteristic of technologically developed countries. Without a subsystem of centralized environmental monitoring, it is impossible to ensure the transition of enterprises away from traditional energy resources. Increasing the intensity of integration and the combination of diversified resources using the goal-foresight-controlling mechanism (FCM) in the digital platform of the Center for Coordinating the Interaction of Diversified Enterprises (hereinafter, “the Center”) can provide the high-tech development of enterprises [
15]. This led to the concept of developing a set of mathematical methods to approximate step functions which reflect the modeling of the effects of factorial parameters in relation to technical and thermodynamic controllers of the 3-E efficiency of enterprises as part of one or more complexes of an EPS. Thus, in this study, we propose the use of a theoretical model of the methodology of integration-balancing management (MIBM) for the sustainable development of industrial enterprises in an emerging EPS [
16,
17].
The proposed methodology has improved models and methods for assessing the efficiency and speed of enterprises transitioning away from low-tech levels of the industrial economy, leading to the predominance of “creative destruction” scenarios. Elimination of obsolete technologies, or their minor modernization by factors of volatility, uncertainty, complexity, and ambiguity (VUCA), are commonplace [
18]. Regulation based on deviations and the results of operational control reduces the circle of supporters of long-term scenario planning [
19]. In conditions of environmental uncertainty and crises, foresight-management methods are needed to predict sustainable high-tech development. We must also take into account risky, complex alternatives aimed at radical technological changes [
20]. Thus, methods for visualizing development trends and normative forecasts are becoming more widespread.
The need for continuous correction of the existing indicators of monitoring weak signals has been substantiated. This proposal has been implemented in collaborative networks by screening robust portfolio models to identify promising innovative ideas [
21,
22,
23,
24,
25,
26,
27,
28]. The composition of external VUCA factors was taken as a foundation. VUCA factors should be accounted for in the design of a successful foresight. The proposals of the aforementioned authors improve control based on cyber-physical systems and computer calculations, using programming tools and artificial intelligence. This will contribute to machine-to-machine substantiation of the variety of options for organizational, economic, and physical impacts and the interaction of subsystems.
These works substantiate the relevance of applying the MIBM organizational approach to ensuring the integration of resources and combining methods of mathematical and thermodynamic modeling of high-tech system development. This allows for an extended set of target characteristics of effective strategic development (the aforementioned VUCA factors) to be implemented [
29,
30,
31]. Thus, the possibilities of digital controlling and foresight methods are expanded, since non-linear processes are taken into account. We applied the integration of big data and data science tools to balance the interests of the strategic development of enterprises and elements of the EPS in agglomerations. The algorithmic scheme will make it possible to carry out research, design, and development in the innovation cycle and implement strategic plans and projects while balancing the goals of the 3-E efficiency of the EPS and the knowledge economy [
32,
33].
The possibility of combining new and old functions when modeling evolutionary modernization and abrupt high-tech transformations has determined the expanded interaction of enterprises within the ecosystem. This substantiated the goal of this study to analyze mathematical models and assess the possibilities of their thermodynamic and organizational interpretation in terms of 3-E efficiency indicators and the stability criterion. The evolutionary and abrupt transition towards the use of technologies with the possibility of integrating traditional and alternative energy resources must be mathematically modeled. To this end, we developed new methods of organizational design of a digital simulator of a single FCM, combining the capabilities of extended environmental monitoring, controlling, and foresight, and the goals of the strategic development of enterprises in the EPS. In modern conditions, these subsystems function separately in a number of enterprises, without high-tech interaction with the knowledge economy.
The goal of expanding the space for control and regulation of decarbonization determined the need to implement extended monitoring of environmental factors. To achieve this goal, strategic management tools for decarbonization goals must be combined with foresight management of long-term development trends. In addition, decarbonization technologies should be developed in the digital platform of EPS regulatory structures across the applicable region. We substantiate methods for displaying singular processes of improving the balance of enterprise goals while coordinating the impact on the efficiency of standard and additional management functions. Impacts are evaluated using mathematical and thermodynamic models that reflect the processes of regulation of the heat and power equipment of industrial enterprises.
An analysis revealed the possibility of adapting and upgrading general and special mathematical models to improve environmental monitoring and regulation of technologies at enterprises in the EPS. Models for the formation of the FCM in the structure of the Center, to increase the efficiency and sustainability of decarbonization within a region, were also determined. For ecological stability, this regional regulator should offer forecasts to enterprises on 3-E goals on the efficiency of processes and results of decarbonization in the medium and long term. We propose using the previously developed theoretical FCM model and the methods of organizational design, which were improved based on this same model. Factor indicators of the effectiveness of managing the integration of resources and combining technologies of enterprises with the potential of the knowledge economy are developed when combining the capabilities of these subsystems in a single mechanism. The indicators allow the governing structures of a region to correct the imbalance of interests detected in the target recommendations for changes to 3-E efficiency indicators and to improve sustainability.
The methods for achieving the set goal are presented as an algorithmic scheme of an FCM in the Center’s digital platform with the cyclic use of big data information tools, data science mathematical and software tools, and strategic management. The scheme is distinguished by a special evidence base for the analysis and application of mathematical possibilities for modeling evolutionary and abrupt processes of 3-E efficiency regulation. The effects of regulators are described by thermodynamic formulas, which make it possible to evaluate the efficiency and speed of the transition towards a zero-carbon footprint.
These innovations are necessary for a number of reasons. The residents of major cities are exposed to increased atmospheric emissions of carbon dioxide and other more harmful components. To reduce emissions, it is necessary to modernize the systems for automation, control, and forecasting of greenhouse gases and waste for disposal. Industrial production will only increase in the coming years, which will lead to an increase in pollution.
We used big data information tools to integrate distributed data and take into account environmental quality parameters over a long period. The automatic control of environmental quality monitoring and diagnostics of the functional and technical state of the elements of the micro- and macro-network of the EPS are implemented using sensors, the internet of things, and neural networks of artificial intelligence.
We developed a concept of expanding the space for control and regulation of decarbonization processes in an EPS. This task is solved by combining strategic management tools for decarbonization goals with foresight management of long-term trends. We substantiate methods for displaying singular processes of improving the balance of enterprise goals while coordinating the impact on the efficiency of standard and additional management functions. We describe methods for mathematical modeling of greenhouse gas emissions and for approximating functions in the study of the thermal power industry and economies. We established that new research methods and techniques increase the efficiency of designing an EPS and can be used in practice to improve the environment. The processes of the digital twin should be modeled with consideration of assessments on ensuring the stability of the balance area with the goals of developing the EPS and the influence of diversified VUCA factors.
2. Materials and Methods
Due to the shortcomings in the use of existing signal discretization methods, we propose a combination of numerical and analytical methods of mathematical modeling based on new mathematical relationships that account for the non-stationary nature of the processes being studied. In other words, we revealed the necessity of accounting for conjugation conditions and the variety of boundary conditions. The spasmodic nature of a number of processes and the implicit oscillations within the model due to the inertia of data transmission necessitated the use of stepwise display functions.
2.1. New Methods for Approximating Step Processes: Functions by the Nyquist Theorem
The Gibbs effect, seen when approximating functions by trigonometric expressions, also forces us to be critical of the proof of some important theorems, which reduce the accuracy of calculations. In particular, Kotelnikov’s theorem, known in the English literature as the Nyquist theorem, is widely used in the theory of signal transmission. When proving the Nyquist theorem, the so-called integral sine is used to approximate functions, which is defined by the expression
Based on the integral sine, Kotelnikov constructs a function to prove the theorem
where ω is an argument and
T and ω
1 are some parameters. At the same time, he claims that with an increase in
T, this function tends to the limits shown in
Figure 1a, and is equal to zero for ω > ω
1 and π for ω < ω
1.
S.V. Alyukov proved that the graph of a limit function will appear as seen in
Figure 1b. In other words, for any values of the parameter (even arbitrarily large, but finite values) there will always be
for which the values of the function constructed by Kotelnikov will differ from π, and there will always be ω > ω
1, for which its values will be non-zero. Moreover, it is important to note that the indicated difference with increasing
T does not tend to zero, but to approximately 0.281 [
6,
7,
8,
9,
10,
11,
12,
13,
14,
15,
16,
17,
18,
19,
20]. Mathematical tools for approximating functions make it possible to evaluate and more accurately control jump-like processes.
Expanding function (1) in a Fourier series has all the above disadvantages, namely the Gibbs effect and a significant approximation error. To eliminate these disadvantages, we propose approximating the original step function by a sequence of recursive periodic functions
Since, due to the limited functions
and the relation
due to the monotonicity of functions
on the interval [0, π/2], we obtain that
Based on the Helly theorem, in a sequence
we can extract a subsequence converging at every point [0, π/2] to some function
f, and
Let us show that the original function can act as such a function, f.
Statement 1. The sequence of functions converges to the original function , and the convergence is pointwise, but not uniform.
Proof. At points
and
Therefore, at these points
since
We can take, for example,
.
Since
then the condition
is satisfied for any
.
Then, the sequence
is positive, increasing, and bounded, and therefore has a finite limit, which we denote by
We obtain
finding that
or
.
Since the sequence is positive and increasing, then
. Then, in the considered interval
Taking into account the previous conclusion on the convergence of the sequence at
and
, we conclude that
This convergence is only pointwise and not uniform, since the function is not continuous on the segment . □
Statement 2. In the space of measurable functions , and in the Hilbert space , the sequence of approximating functions converges in the norm to the original function .
Proof. We introduce a sequence of minorants with respect to a sequence
of functions
The set of discontinuity points of the function
is equal to zero. Then, taking into account the non-negative sign and boundedness of the functions
and
, and of the segment under consideration, in space
we obtain
Similarly, we can prove that the sequence converges in the norm to a function in the space .
Thus, the sequence of approximating functions in spaces and is fundamental. In space , the sequence is not fundamental. □
2.2. The Use of New Methods for Approximating Stepwise Processes in the Energy Sector
Analytical methods must be used to provide a greater depth of process analysis, accuracy, and convergence of solutions. In the considered non-stationary processes, we propose the application of step functions. Their traditional areas of application are technical and mathematical disciplines, including the theory of automatic control, the theory of information and signal transmission, equations of mathematical physics, differential equations, and many others. Currently, only the initial steps are being applied to systems with stepped characteristics and functions in thermal power engineering. These systems are referred to as essentially non-linear structures, the solutions of which are difficult to obtain.
Organic Rankine cycles (ORCs), which increase the 3-E efficiency of thermal power plants and industrial plants, can be described with a combination of the above proven mathematical possibilities and thermodynamic dependencies. For example, a group of researchers led by Y. Zhu analyzed the operation of a gas piston unit using thermodynamic and economic methods to introduce the organic Rankine cycle. They identified the most obvious heat losses, which can be used as a heat source for the operation of a heat pump. Other plants produce compressed air. The proposed work combines ORC and the compressed-air generation cycle in a single EPS using an improved thermo-transformer, namely a refrigeration and heating unit operating on the heating water of a thermal power plant.
In addition to the described methods of decarbonization, the same effect is possible in the power industry. For example, another class of parallel mechanism and an electric drive based on switched reluctance motors with electromagnetic speed reduction and an extreme angular position control system, can be combined in a single EPS. Each of these elements solves one of the tasks assigned to it: a parallel mechanism increases rigidity and resistance to wind loads, an electric drive enables operation at low speeds. Extreme control is required; such installations must be installed to achieve the highest possible energy efficiency. For example, when managing collectors, the position is regulated when the negative impact of the environment is minimal. The aforementioned methods of control require the movement of individual parts of the mechanized system to be described. In this case, Alyukov’s method is the most applicable.
2.3. Methods for Substantiating the Composition of the EPS in Terms of Heat and Power Indicators and Efficiency Criteria
To select the thermal power of the source and the power supply system, we use the functional of energy costs from the thermal power of the cooling system
The heat balance (24)–(29) for the air conditioning system of an industrial facility, due to the use of thermal energy in an EPS, can be written as
where
are the temperatures at coordinate x on the thickness of an external wall, on an internal surface, and on an external surface of an external wall of the building at timepoint
τ, respectively; and
a is the thermal diffusion.
where
is the coefficient of heat conductivity of the material of an external wall;
is the thermolysis coefficient for internal and external surfaces of a wall; and
and
are the temperatures of internal air at the initial timepoint and timepoint
τ, respectively.
is the thermolysis coefficient for the external surfaces of a wall; and
is the temperature of external air at timepoint
τ;
where
is the management heating system power at timepoint
τ;
and
are the specific heat and mass of the internal air in the building, respectively;
and
are the areas of the exterior walls of the building and windows, respectively; and
is the window heat transfer coefficient.
Functional
I consists of two parts—temperature difference and electricity costs:
where
is the set period;
is the function estimating the current cost of heating the building; and
is a weight coefficient.
To formulate optimality conditions, we introduce the function (31)
satisfies the following system of Equations (32) and (33)
with initial conditions (34), (35):
The boundary conditions for the function
are given as:
2.4. Methods for Calculating Greenhouse Gas Emissions of an Industrial Facility Due to the Use of Thermal Energy in an EPS
The dependence for the emissivity of carbon dioxide is represented by three functionals
where
are functions derived from experimental data.
where
,
are functions derived from experimental data.
where power
where power
Similarly to (38)–(44), for water vapor
:
where
are functions derived from experimental data.
where
are functions derived from experimental data.
where power
where power
where powers
where powers
When calculating the total radiation of a mixture of gases, we must account for the mutual overlap of the CO2 and H2O absorption bands. The calculation of the emissivity of CO2 is given as an example, using the parameters lp = 0.04, θ = 1373 K (t = 1100 °C), εCO2 = 0.149, and ECO2 = εCO2·5.67θ4 = 3 × 104 W/m2.
2.5. Methods for Forming the Foresight-Controlling Mechanism for the Decarbonization of Enterprises in the EPS
The proposed MIBM model substantiates a set of methods for modeling complex socio-technical systems in the EPS with a transition to a zero-carbon footprint. These methods make it possible to diagnose the conditions for the inclusion of low- and medium-tech enterprises in the EPS to increase interaction with educational and scientific institutions in a post-industrial agglomeration. This determines the parameters for the synthesis of a single FCM to increase the 3-E efficiency of enterprise development in the system being formed. We understand the mechanism as a complex analog–digital toolkit for the systematic organization of the coordinated actions of enterprises and knowledge economy based on the results of expert communications (at least as part of an EPS). Experts should determine the initial estimates and adjust subsequent ones to follow the ideal vision of the long-term goals of enterprises in the EPS. To this end, organizational design of evolutionary and spasmodic processes of the coordinated development of enterprises was carried out, keeping in mind the goals of decarbonization.
Enterprise managers and the regulatory regional body of the Center for forecasting the target indicators of the vision in the assessments of 3-E efficiency make decisions on adjusting target indicators. Designing is carried out according to the stages of sustainable and planned development and takes into account possible bifurcations of trends of the indicator under the conditions of maximum uncertainty using an analogy with the “just in time” logistics system, the forecasts of ideal indicators, the existing target indicators of the resulting property. Factor indicators of management quality are adjusted based on the results of periodic expert assessments in iterative processes of applying the FCM to improve the 3-E efficiency. Therefore, new monitoring subsystems and mechanisms should be included in the digital platform of the Center. The Center can correct imbalances in the interests of enterprises in an industrialized region. We propose considering the organizational models of cyclical regulation of multivariate processes using the big data and data science tools outlined below.
The stability properties of increasing 3-E efficiency justify our proposed principles for regulating the imbalance of enterprise interests in the FCM mechanism for their coordinated development in an EPS. An expanded list of external (EVUCA) factors is accounted for: variability in environmental quality indicators; uncertainty in indicators; complexity of influencing factors; simultaneous changes to three types of efficiency; ambiguity of effects; and situations where energy efficiency increases while environmental and economic efficiency decreases.
The internal (IVUCA) view of factors as organizational impacts of the Center, characterizing the environment of enterprises and elements of micro- and macro-networks in the EPS, include: a shared vision among the staff on the future state of the EPS, to counteract instability in assessing the availability of 3-E efficiency targets, keeping decarbonization in mind; personnel’s understanding of the goals of joint development in the EPS; operational adjustment of strategic development trends in response to a clearly established discrepancy between actual and target indicators of 3-E efficiency as an imbalance of enterprise interests with the sustainable development and the speed of decarbonization; the speed and flexibility of management decisions to counteract the ambiguity of results based on the continuous study of the possibilities of innovative development.
Methods for combining and integrating the resources of microgrid objects, macrogrids, and industrial production technologies lead to more complex processes and require the use of more complex data science tools. This in turn requires the use of more sophisticated data processing tools.
Figure 2 shows the scheme of a typical EPS in an industrial region in terms of the interaction among institutions using diversified energy sources.
The model of the foresight-controlling mechanism for the facilities of the EPS is similar. Therefore, we understand the mechanism as a complex analog–digital toolkit for the systematic organization of the processes of coordinated actions of enterprises and the knowledge economy based on the results of expert communications (at least as part of the EPS). Experts determine the initial and corrected subsequent estimates of the ideal long-term goals for the energy system of the region at different stages of an enterprise’s innovation cycle. The mechanism makes it possible to implement the evolutionary and spasmodic processes of regulating 3-E efficiency of enterprises in the EPS based on the results of the joint actions of experts. Enterprise managers and the regulatory regional body of the Center, for forecasting the target indicators of the vision in the assessments of 3-E efficiency, make decisions on adjusting target indicators. The design is carried out according to the stages of sustainable and planned development and takes into account possible bifurcations of trends of the indicator under the conditions of maximum uncertainty.
Figure 2 shows the main elements of an EPS that uses gas piston units, waste heat boilers, solar photo panels, refrigeration equipment, and heat pumps. The greenhouse gas emissions of the equipment and thermal energy were calculated according to (23)–(58).
Figure 3 shows a diagram of interactions between EPS enterprises based on the FCM in the regulatory structure of the regional coordination center.
Figure 3 shows the interaction of the elements of the EPS. Energy flows were calculated using mathematical methods according to Formulas (1)–(22).
To achieve investment results to implement high-tech development projects, it is necessary to use existing software products that ensure their variability. Regression analysis and the considered mathematical models expand the possibilities of multivariate and multicriteria approaches on digital platforms, without primary filtering of the source data.
3. Practical Use of New Methods of Thermodynamic Calculations in the Design of Power Supply Sources
Elements of an IoT subsystem are needed to organize the interaction of facilities in the electric power system. These include autonomous sensors with self-diagnostic functions connected to a system of automatic regulators of technical parameters based on the results of monitoring environmental quality indicators. We investigate the modeling of processes and systems of a zero-carbon EPS in terms of energy and environmental sustainability of heat and cold supply processes. The proposed models use controlled cooling models and organizational models which account for the air temperature in the cooled device and the temperature of the cases. They also utilize advanced environmental monitoring and control capabilities to continuously assess energy savings and the carbon footprint. Using past experience and calculations, the energy consumption and cost savings associated with decarbonization can be quantified. In addition, the proposed models take into account factors such as power requirements and technological advances. These big data and data science tools can identify and evaluate possible trends and options for improving 3-E efficiency in energy use, costs, and carbon emissions. Applying this level of monitoring in innovative decarbonization cycles allows for the best strategies for reducing energy and carbon emissions to be identified and implemented.
The proposed EPS models have the potential to significantly reduce Jaro’s energy consumption and carbon emissions. According to the experience of technologically advanced countries, energy savings associated with the decarbonization of the energy sector can range from 2% to 5% in the short term to 50–75% in the long term. For example, replacing electric lighting with LED lighting can lead to energy savings of up to 50–60%. Similarly, implementing energy management systems (EMS) can lead to 10–20% energy savings. The quantitative analysis of the proposed models also showed that the introduction of energy efficiency measures alongside renewable energy sources can lead to emission reductions of 70–100%. In addition, modeling the processes and systems of a zero-carbon EPS can lead to cost savings of up to 15–25%. It can also lead to a local reduction in air pollution of up to 40–60%. Finally, the implementation of the proposed models can lead to an increase in jobs in the domestic energy sector, with a potential employment growth of up to 20–30%.
Figure 4 shows a diagram of the use of coal–water fuel as an energy source in an EPS. The designations of the objects of the EPS using coal–water technologies are introduced. The technologies reduce the environmental impact of greenhouse gases. The proposed scheme can reduce the carbon footprint of EPS enterprises and increase the 3-E efficiency of energy production and consumption processes by integrating industrial and domestic energy installations in one ecosystem. Organizing a single mechanism of the FCM increases the speed of the effective and sustainable development of enterprises when it is included in the digital platform of the center for regional coordination.
For the technical implementation of decarbonization, it is necessary to develop a scheme of automatic control regulators according to the process stability parameter. An ecosystem can be created using the capabilities of sensors such as autonomous wireless sensors with self-diagnostic capabilities in the industrial IoT. This corresponds to the goals of the green transition based on monitoring of environmental quality parameters. To increase the speed and reliability of regulatory processes, it is necessary to add tools for intellectual and analytical management based on neural networks of artificial intelligence. This also provides a variety of options for management decisions to improve the quality of forecasting and planning the development of EPS facilities. This is ensured by the regulation of technical parameters based on the results of automatic remote evaluation of the dynamics of environmental quality indicators in decarbonization processes. These include an increased area of the rhizosphere, an increased volume of emissions of greenhouse gases in a territory, the burning of pumped out volumes of biogas from landfills, and other indicators.
The results of regulation were evaluated by the criterion of minimizing the deviations of the actual parameters of environmental quality and efficiency from the standards of their ideal or target level of forecasts/plans. This made it possible to carry out a systematic analysis of the insufficient efficiency and sustainability of processes in a post-industrial ecosystem.
The speed of transportation is selected based on technical and economic considerations, the minimum energy consumption for the movement of solid material, and the minimum wear of pipelines. In particular, A.E. Smoldyrev recommends the following transport rates for coal–water suspensions [
34] (
Table 1).
The proposed integrated forecasting and control system will increase the speed and efficiency of meeting the post-industrial demand for reclamation of existing landfills and conversion into high-tech carbon landfills. In this regard, it is vital to create an effective regulatory framework to guarantee the security and sustainability of energy supply processes, including the environmental, economic, and social (3-E) elements of energy consumption.
Modernization options, formulas, algorithms and calculation results are of great relevance to the sustainability of the EPS with a zero-carbon footprint. They can help to model and optimize complex processes and systems and identify the most efficient energy sources and technologies and the most cost-effective operation and management strategies. Ultimately, this can lead to lower energy fuel costs, reduced environmental pollution, and increased sustainability of the complex of energy technologies.
Figure 5 shows a diagram of the use of heat, including from water–coal fuel and a gas piston plant, using traditional energy resources.
Figure 5 shows the most promising scheme for using a gas piston installation based on the results of calculations using Formulas (23)–(58).
Integrating all available energy sources into the regional macrogrid ecosystem allows for greater energy efficiency. Turning sources on increases the accuracy of control systems, allowing for more efficient management of heating, cooling, and the power supply. Using a control scheme and balance integration methods increases 3-E efficiency under the conditions of advanced control of parameters and environmental factors.
It should be noted that a computer simulation was also carried out in a special program, in which the boundary and initial conditions were taken as those during the combustion of coal dust in the furnace of the E-160 boiler unit,
Figure 6. Analysis of the computer simulation data showed a discrepancy of 3–4% with the calculated data and 5–6% with the experimental data.
Since we approximated the generalized functions by analytic functions, we can differentiate these approximating functions and find their derivatives of any order. The main functions of the parameters for k-omega computer simulation of the quality and speed of decarbonization processes can be approximated by a principally new method. The resulting approximations more accurately describe the original function. The suggestions for improving the mathematical apparatus of the standard k-omega model will help researchers to develop the theory of approximation by the recursive functions with an effect-corrected coefficient estimate of 0.281, determined by S.V. Alyukov. This coefficient solves one of the important problems in achieving more accurate modeling, by reducing calculation time and improving suitability of the constructed model for similar installations and devices.
We determined that, with the help of new methods of approximation and analysis of deviations of the air flow rate in the channel, it is possible to obtain more accurate results of experimental and numerical studies (
Figure 7).
Figure 7 is a 3D distribution of deviations in the calculation of velocity of slurry movement through innovative pipelines. Deviations in the calculation of velocity may be explained by the fact that measurement results are affected by a large number of sources of small random errors. The whole set of measurements has a symmetric bell-shaped Gauss function as a limit distribution. The distribution center, which coincides with its maximum, is the true measured value. We use S.V. Alyukov’s method of approximation, which is considered for the first time in relation to numerical methods and experimental data for automatic control applied to foresight controlling of 3-E efficiency and speed of reclamation of decarbonization (see
Figure 7).
When using these new approximation methods, only small relative errors occur when transmitting experimental data. The obtained modeling results in ANSYS show a good fit to the theoretical and experimental data. The differences are insignificant, since errors in experimental measurements were leveled by using modern tools certified by the European Union. The error in the mathematical modeling was reduced by using empirical data obtained by processing signals from analogous equipment and subsequent sampling, according to the technology developed by us in previous studies. This can be seen, for example, in the decarbonization velocity (
Table 2).
Table 2 shows some indicators of improving decarbonization processes according to the increase in the 3-E efficiency criterion for facilities of the EPS.
The indicators were achieved by combining automatic control according to the scheme in
Figure 7, the technology of a standard trigeneration plant (
Figure 5), and by strengthening the relationship between an industrial enterprise and a consumer (
Figure 4). In addition, resources are implemented through the integration of the relevant energy facilities in the EPS and are shown in
Figure 2. The organizational conditions for combining and integrating diversified resources are shown in
Figure 3. They are implemented in the algorithmic scheme of interaction between enterprises as part of the foresight-controlling mechanism.
The maximum differences are obtained between the experimental data and computer modeling data when our developments are not used. The maximum coincidence is seen between the average values according to the modeling results, taking into account the sampled signals and the experiment (
Figure 8).
4. Discussion
Sen and Lalu studied holacratic flexible (agile) methods for creating self-managing structures [
35,
36,
37]. In comparison with Adizes, Sen and Lalu quite unambiguously determine the types of such structures—they are spiral dynamics models [
37]. However, the quality of control in terms of heat engineering in such structures greatly narrows the applicability of the methods in the formation of an EPS. Real forecasts for the development of intellectual resources in heat engineering are provided by Frumin [
38]. The growth in high-tech engineering services in the innovation ecosystem is in line with the mentioned 3-E concept. The transition of small enterprises of the energy complex to educational and scientific systems 4.0 is growing, along with universities distinguished by high academic achievements and the use of new educational methods using big data. Currently, the use of training digital twins in virtual or real cognitive methods that adhere to data science is recognized for solving practical problems of the post-industrial economy and energy. In the post-industrial digital economy, it is necessary to use open educational platforms. When studying the quality of the results, Azizi et al. took into account the necessary organizational tools for the development of complex logistics systems [
39]. Hoegl and Schulze studied the initial stage of the genesis of the innovation cycle [
40]. Liu et al. revealed the importance of prognostic methods [
41]. The modeling of adaptive systems was supplemented by Lafond [
42] and McCarthy et al. [
43]. Kirsch proposed a solution of complex control problems, including those applicable in the energy sector [
44]. Puschke et al. presented scientific approaches to dynamic optimization of processes [
45]. Promising methods for analytical modeling are presented in [
46,
47]. An overview of these papers is included in this section so that readers can become acquainted with the further development of this area of research.
Technical regulators and process automation tools should have the properties of foresight controlling of 3-E efficiency goals. To do this, they must have the following elements and properties:
Analysis of the motion of water–coal suspension flows in
Figure 8 can be described by the following mathematical model. According to the perturbation method, velocity can be represented as the sum of the stationary and perturbed components:
Let us introduce the dimensionless stream function:
Then, the perturbations of the velocity and the interface are expressed as:
where
c is a parameter and
is the wavelength.
Substituting the relations into the Navier–Stokes equation, we obtain
In this case, the boundary conditions are expressed as
We use the Orr–Sommerfeld equation and the corresponding boundary conditions
where
is the dimensionless stream function,
is maximum speed,
where
is the dimensionless parameter of the normal direction; ρ″ is the density of the gas, (kg/m
3); g is the gravitational acceleration, (m/s
2); ν is the coefficient of kinematic viscosity of the coolant in the liquid state, (m
2/s); and β is the coefficient of volumetric expansion, (K
−1), selected from the reference book;
where
is the coefficient of surface tension,
where
is the surface viscosity parameter; and
and
are coefficients,
where
is the parameter of the inhomogeneity of the surface,
where
is the dimensionless shear stress parameter and
is the dimensionless parameter of dynamic viscosity;
where
and
are the dynamic viscosities of the liquid and gas.
The dependencies obtained by Benjamin in Craik’s modification can be written as:
where
is the coefficient of friction and
is the gas flow velocity on the film surface.
We represent the function as an expansion in a small parameter
and the parameter
appears as
The boundary of film flow stability is determined by the relation
where
5. Conclusions
We substantiated the concept of expanding the control and regulation of decarbonization processes in an EPS. This determined the purpose of the development of methods, including extended monitoring of environmental factors. The set problem was solved by combining strategic controlling tools for decarbonization goals with foresight management of long-term trends to develop the use of decarbonization technologies in a formed digital platform of EPS regulatory structures across a region.
We substantiate methods for displaying singular processes of improving the balance of enterprise goals while coordinating the impact on the efficiency of standard and additional management functions.
We noted the undoubted relevance of using mathematical methods and models for studying the regulation of factorial parameters of the processes of decarbonization of a region and reducing environmental pollution. We described methods for mathematical modeling of greenhouse gas emissions and for approximating functions in the study of processes in the thermal power industry and the economy. The established new research methods and techniques increase the efficiency of designing an EPS and can be used to improve the environment.
We developed a methodology in accordance with the algorithmic scheme of the foresight-controlling mechanism for the efficiency and speed of decarbonization. Digital twins should be modeled with consideration of assessments on ensuring a balance between the goals of developing the EPS and the influence of diversified VUCA factors.
The simulation results can be interpreted in accordance with the distribution of target performance indicators by levels of the Harrington scale (
Table 2). We proposed and used mathematical constructions and models to assess their growth when using the controlling and foresight-management subsystems.