Sustainability of Management Decisions in a Digital Logistics Network

: Globalization has given a powerful impetus to the development of international commercial activity and logistics management systems taking full advantage of cross-border networking. The solution lies at the intersection of information technologies, technical means of machine-to-machine (M2M) interaction, mobile high-speed networks, geolocation, cloud services, and a number of international standards. The current trend towards creating digital logistics platforms has set a number of serious challenges for developers. The most important requirement is the condition of sustainability of the obtained solutions with respect to disturbances in the conditions of logistics activities caused not only by market uncertainty but also by a whole set of unfavorable factors accompanying the transportation process. Within the framework of the presented research, the problem of obtaining the conditions for the stability of solutions obtained on the basis of mathematical models is set. At the same time, the processes of transferring not only discrete but also continuous material ﬂows through complex structured networks are taken into account. This study contains the results of the analysis of the stability of solutions of differential systems of various types that simulate the transfer processes in network media. Initial boundary value problems for evolutionary equations and differential-difference systems are relevant in logistics, both for the discrete transportation of a wide range of goods and for the quasi-continuous transportation of, for example, liquid hydrocarbons. The criterion for the work of a logistics operator is the integral functional. For the mathematical description of the transport process of continuous and discrete media, a wide class of integrable functions are used, which adequately describe the transport of media with a complex internal rheological structure.


Introduction
Different factors could facilitate the transformation of economies to sustainable economies [1,2] as well as bioeconomies [3].
This study presents a thematic description of the process of transfer of continuous and discrete media by formalisms of differential systems with continuously varying time and differential-difference systems with discretely varying time. In 2010, the authors completed a qualitative analysis of evolutionary transport equations with a spatial variable changing on a network or a network-like domain (the existence and uniqueness of the solution, the correctness of the stated initial-boundary value problems). The latter led to the emergence of further research by the authors that in 2020 ended with scientific results in the direction of studying the properties of solutions (stability, stabilization, and controllability of system solutions) and the practical use of the obtained results, the logical conclusions of which were the results of this work. It should be noted that the use of the mentioned mathematical results is a distinctive feature of this work and represents an advantage over most studies of an economic orientation, in which the issues of the stability of mathematical models of an economic nature, and therefore the stability of processes, are extremely rarely and superficially studied or are not considered at all. The authors also note that the results obtained are new-the authors are not aware of similar publications in the studied direction.
The development of correct mathematical models reflecting the evolution of material flows in networks is a complex task and requires the involvement of a serious mathematical apparatus [4]. For the found theoretical solutions and algorithms, real-world business puts forward the additional condition of the study of their sustainability. It is necessary due to such factors as the instability of market conditions, the spread of parameters of mathematical models, and the inaccuracy of the initial data.

The Purpose of the Research and Methods
The goal of the study, due to its complexity, includes two tasks. The first task is to find the stability of the solution of an evolutionary parabolic system with continuous time. The system has distributed parameters on the graph used as a prototype of the logistics transport network [5,6]. The mathematical model itself describes the process of transferring material flows through a spatial network. The parabolic system is studied using the class of integrable functions reflecting mass transfer with a complex internal rheology. A generalized solution or weak solution of the system in this case is an integrable function that determines the variational formulation of the initial boundary value problem.
The second problem involves the analysis of an analogue of the initial boundary value problem for a differential system of parabolic type with a discretely varying time variable: a differential-difference system of equations with a spatial variable changing on a graph (network) that defines a mathematical model of the transfer process in discrete time. To construct such an analogue, we use the E. Rote method of semi-sampling over a time variable. The conditions for the existence of a weak solution of a differential-difference system and the analysis of its stability are obtained. The results of the study can be used in problems, such as initial boundary value problems and differential-difference systems of equations, of the optimization of network-like transfer processes.
One of the most relevant problems is the discrepancy between the development of hardware and brainware. The problem is that although hardware for logistics platforms is represented by a full line of modern technical means for tracking the flow of goods and cargo, as well as the exchange of information between all participants in the logistics process, the brainware algorithms have only just begun to develop in the direction of network logistics. For practice, we need algorithms that, first of all, are optimal by economic criteria and acceptable for embedding in the software of servers of digital logistics platforms.
The following parameters are noted separately: the traffic level; the capacity of the routes; the spread of the waiting time and the process of loading and unloading operations; customs services when crossing state borders; the influence of weather conditions on the speed of movement; technical failures of carriers; and other factors known to logistics management. In addition, the costs are affected by fluctuations in the cost of stevedoring operations, mileage rates in multimodal logistics, a difficult-to-predict range of prices for the rental of port and station infrastructure, etc.
Their impact can be assessed using the theory of random functions and the theory of queuing (or queuing theory, QT). Traditionally, for these purposes, the methods of A. M. Lyapunov have been used to assess the sustainability of systems described by differential equations for small perturbations of the initial condition. The development of this approach lies in the area of logistics networks modeling taking into account the sustainability of material flows. This improvement considers the methods that use maps, graph theory, and systems of equations from the apparatus of mathematical physics and applies methods of the theory of optimal decision-making.
Within the framework of the presented research, the problem of obtaining the conditions for the stability of solutions obtained on the basis of mathematical models is set. At the same time, the processes of transferring not only discrete but also continuous material flows through complex structured networks are taken into account. This study contains the results of the analysis of the stability of solutions of differential systems of various types that simulate the transfer processes in network media. The problems under consideration are relevant to logistics, both for the discrete transportation of a wide range of goods and for the quasi-continuous transportation of, for example, liquid hydrocarbons.
A generalized solution is a function that defines the variational formulation of an initial boundary value problem or a differential-difference system. Going beyond the classical continuously differentiable solutions is dictated by the need to describe more precisely the physical nature of the processes of transport of continuous media (e.g., gas, oil, petroleum products) or discrete flows (commercial goods and cargo) through logistics networks.
The analysis of the influence of perturbations of the source data on the behavior of a real logistics system and an abstract mathematical model (the solutions obtained either in the class of initial boundary value problems or the system of differential-difference equations used) is critical primarily from the point of view of economic indicators of commercial activity. This makes it possible to assess the risks for investors as well as analyze the sustainability of the proposed business models.
Note that in the class of ordinary differential equations, the results of the mathematical theory of stability have been obtained and are widely applied. The new tasks, dictated by the need to manage economic processes, require optimization according to efficiency criteria containing at least two variables, one of which is the current time of the process. In such economic applications, due to the increasing complexity of algorithms for modeling evolutionary transfer processes, there is a natural need to use formalisms of equations and systems of partial differential equations that are more appropriate to the specifics of evolutionary processes.
In the present paper, the abovementioned formalisms are in practice taken to be the main ones in the mathematical description of evolutionary processes on networks, and their analysis is the main goal of the study. Namely, we will obtain the sustainability conditions for solutions for two classes of evolutionary differential systems-systems with a continuously varying time variable and systems with a discretely varying time variable.

Main Definitions
We use the classical notation and concepts introduced by the authors to describe differential systems with a spatial variable changing on a graph Γ (in applications on a logistics network). At the same time, the internal structure of the graph can be arbitrary (the graph can contain cycles) while remaining connected and bounded.
Due to the use of functions with a carrier on the graph (and, consequently, differential equations with a variable changing on the graph), at each edge of γ ⊂ Γ a parametrization is introduced by the segment [0, 1].
The case of parametrization by different segments is not excluded. In accordance with the notation introduced in [3,4], ∂Γ, J(Γ) are the sets of boundary ζ and internal ξ nodes of Γ, respectively.
By Γ 0 ⊂ Γ we denote a graph that does not include boundary and end nodes: T < ∞-constant). We use the concept of the Lebesgue integral on the graph Γ and on the domain Γ t : f (·) γ is the narrowing of the function f (·) to a fixed edge γ. We use the classical space C(Γ) of continuous functions on a graph and the Lebesgue space L p (Γ) (p = 1, 2) of pth-degree measurable and summable functions with norm 1/2 and the space L p (Γ T ) with norm: as well as a space L 2,1 (Γ T ) of functions from L 1 (Γ T ) with norm: In the analysis of differential equations and systems of equations with a spatial variable changing on the graph Γ, we use the main Sobolev function spaces: W 1 2 (Γ) with elements from L 2 (Γ), for which the derivative of the first order (the generalized derivative) belongs to L 2 (Γ); the space of functions W 1,0 2 (Γ T ) with elements from L 2 (Γ T ), for which the derivative of the first order in the spatial variable x belongs to L 2 (Γ T ), and the space We further introduce the state space with auxiliary spaces. This is necessary not only for the decision-making process, but also for the development of the conditions for the existence of these decisions. To do this, in W 1 2 (Γ), consider the expression: defining a bilinear form with fixed measurable and bounded on Γ 0 coefficients a(x),b(x) from the space L 2 (Γ):

Theoretical Fundamentals
The authors consider a statement that is proved in the work [5].
is continuous for the set of internal nodes of the graph under consideration for an arbitrary edge γ ⊂ Γ.
We introduce a set of Ω a (Γ) functions u(x) that satisfy the statement of Lemma 1 and the following relations, called balance conditions in applications: for any internal node ξ ∈ J(Γ) of the graph (in applications, the nodes of accumulation and distribution of the transported medium of the logistics network); here, R(ξ) and r(ξ) denote the sets of edges γ oriented "to node ξ" and "from node ξ", respectively. The operation of closing Ω a (Γ) by the norm of the space W 1 2 (Γ) leads to the definition of the space W 1 (a, Γ). If we assume that for u(x) from the set Ω a (Γ) there is a boundary condition u(x )| ∂Γ = 0, then another space , whose traces exist on the sections of the domain Γ T with the plane t = t 0 (t 0 ∈ [0, T]) and belong to the space W 1 0 (a, Γ). The following relations are also valid for them: is an arbitrary node of the logistics network). The operation of closing Ω a (Γ T ) by norm (1) leads to the definition of the space We define another subspace for W 1,0 2 (Γ T ) by the closure in the norm W 1,0 2 (Γ T ) of the set of functions that are differentiable and meet the conditions (3) in ξ ∈ J(Γ) and in t ∈ [0, T].
We introduce the notation W 1,0 (a, Γ T ) for it and also define it similarly W 1 (a, Γ T ); in this case, is performed (in applications using this simulation, this corresponds to the absorption condition).
We have defined the V 1,0 (a, Γ T ) domain of states of the differential system as well as additionally W 1,0 (a, Γ T ) and W 1 (a, Γ T ). The main difference between V 1,0 (a, Γ T ) and the spaces The proven statement is of great practical importance in the formation of economic indicators for the evaluation of network logistics. Indeed, the integrated criterion that reflects the cost of moving goods depends on both the distance traveled and the time. Such dependencies have jumps caused by changes in tariffs for multimodal transportation, the approach to the expiration date of a wide range of consumer goods in the fast-moving consumer goods (FMCGs) segment, the terms of the transport lease, port infrastructure, demurrage, detention, or contractual penalties for short delivery or late shipment, and other factors related to logistics activities.

Mathematical Model of Transfer with Continuously Changing Time
The use of the introduced formalisms for modeling the processes of moving commodities through the logistics network structure allows us to strictly mathematically describe the following algorithm.
On the elements of the spaces W 1,0 (a, Γ T ) and V 1,0 (a, Γ T ), we consider the differential equation It is a system of parabolic differential equations with distributed parameters a(x) and b(x), which characterize the quantitative indicators of the flow state y(x, t) (the transfer rate, transport efficiency, and other criteria) along the edges γ of the graph Γ as a portrait of the logistics network f (x, t) ∈ L 2,1 (Γ T ).
The state function (x, t ∈ Γ T ) (x, t ∈ Γ T ) for system (4) in the network domain Γ T is defined by a generalized solution y(x, t) of system (4) that satisfies the initial and boundary conditions of the form: Above, we have already defined the summability with a square for the established functions a(x) and b(x).
From the condition of belonging to y(x, t) ∈ V 1,0 (a, Γ T ) and summability with the square of the functions a(x), b(x), it follows that the state y(x, t), as a map of the economic criterion on the graph, i.e., the map y : , is a continuous function on the graph, so that the first relation in (5) makes sense and is understood for almost all values of its variable.
For simplicity of further presentation, we use the Dirichlet boundary condition (the second relation in (5)), which in the applications means the absence of a flow (inflow) of a continuous medium in the boundary nodes of the network, i.e., we consider the first initial boundary value problem (4), (5).
We present the main statements and fragments of the proof of the validity of the obtained results in the part that is necessary for the study of the sustainability of solutions. This is primarily applied in the analysis of the adoption of management strategies by managers of logistics services and risk management of logistics business processes. To do this, consider the solution first in the auxiliary space W 1,0 (a, Γ T ), then in the main space V 1,0 (a, Γ T ).

Definition 1.
A generalized solution of the initial boundary value problem (4), (5) is understood as an element (function) y(x, t) ∈ W 1,0 (a, Γ T ) that satisfies the identity of the integral type: for an arbitrary element η(x, t) ∈ W 1 (a, Γ T ) equal to zero for t = T; t (y, η) denotes the bilinear form with respect to the elements u, η: It should be noted that the solvability of problem (4), (5) in W 1,0 (a, Γ T ) (and V 1,0 (a, Γ T )) is established by representing the solution of y(x, t) ∈ V 1,0 (a, Γ T ) as a series using the system of generalized eigenfunctions of the boundary value problem (the problem for an elliptic equation): in the class of functions of the W 1 (a, Γ) space [6,7]. It follows that, under assumption (2), the spectral problem (7) has a set of eigenvalues {λ i } i≥1 forming a countable set. The following properties of the spectral characteristics of problem (7) are valid: 1.
The eigenvalues are real and have finite multiplicity, they are numbered in ascending order of modules (taking into account multiplicities): {λ k } k≥1 , each eigenvalue corresponds to its own real generalized eigenfunction, and the set of generalized eigenfunctions forms the sequence {φ k (x)} k≥1 ; 2.
The numbers λ k are positive except for a finite number of the first ones; for b(x) > 0, all numbers λ k are positive; 3.
For the boundary value problem Λφ = λφ + g, g ∈ L 2 (Γ) generated by the differential expression: is an alternative to Fredholm in the W 1 0 (a, Γ) space. Note 1. In applied problems of practical economics, the coefficient b(x) is a positive function, which means that all eigenvalues of the boundary value problem described in the expression (7) are positive. To find the sustainability conditions of evolutionary systems with an elliptic part defined by the differential operator the a priori inequality b(x) > 0 is the main one for establishing the stability conditions. Theorem 1. Let the following conditions be satisfied: is a non-negative function. The initial boundary value problem (4), (5) in the space W 1 0 (a, Γ) (and the space V 1,0 (a, Γ T )) has a unique generalized solution.
To prove unambiguous weak solvability in the space W 1 0 (a, Γ) (and then in the space V 1,0 (a, Γ T )), we use the Faedo-Galerkin method with a special basis-the set of generalized eigenfunctions {φ k (x)} k≥1 of the class W 1 0 (a, Γ) of the one-dimensional elliptic operator generated by the differential expression: The representation of the solution as a series in the system {φ k (x)} k≥1 and the subsequent analysis of the convergence of this series together with its derivative in the time variable complete the proof. The uniqueness of the solution is established by the classical method using the linearity of the elliptic part of Equation (4).

Mathematical Model of Transport with Discretely Varying Time
Next, we consider a mathematical model of the transfer process in a discrete time change. The mathematical description is based on the Rote method of semi-discretization, which allows us to represent the mathematical model by formalisms of a differentialdifference system of ordinary differential equations on a graph, i.e., a system of equations with a spatial variable changing on the graph. In this case, the theorem of the existence of a solution in the space W 1 0 (a, Γ) is established and the ways of analyzing the stability of this solution are indicated. In the space W 1,0 0 (a, Γ T ) of the states u(x, t)) of system (4), we consider the initial boundary value problem (4), (5).
The functions u(k) k = 1, 2, . . . , M are defined as weak solutions of the elliptic equation (8) satisfying the conditions Thus, for a fixed k relation (8), (9) is the boundary value problem for the elliptic equation (8) with respect to u(k). (8), (9) are analogous to the implicit first-order difference scheme of approximation with respect to the time variable t for the initial boundary value problem (4), (5), given in the space W 1,0 0 (a, Γ T ), with the elliptic operator:
Proof. Note, first of all, that from the condition f (x, t) ∈ L 2,1 (Γ T ) and the representation f τ (k), f τ (k) ∈ L 2 (Γ) follows by virtue of the very definition of the space L 2,1 (Γ). Let k = 1. Based on properties 3 and 4 of the boundary value problem Λφ = λφ, φ| ∂Γ = 0 , and the relation: the statement for the function u(1) follows. The same statement, given the ratio remains true for k = 2, 3, . . . , M. Below, when obtaining an a priori estimate for the functions u(k), the boundary τ 0 for the change τ will be specified. The theorem is proved.

Sustainability of the Control Solutions of the Transfer Processes
The stability of mathematical models of the transfer process (controlling the decision process) with continuously and discretely changing time is considered below. Note that the term "process control solution" in applied problems of practical economics means fixing (selecting) the initial data of the problem in accordance with the current economic need and the capabilities of the economic entity.

Sustainability of Solutions with Continuously Changing Time
We will conduct a mathematically based study of the stability parameters of the found solution of system (4), which describes the dynamics of the logistics process, formalized in the form of mass transfer of goods and cargo through the network. The transport network, in turn, is represented as a mathematical graph.
Assume that 0 ≤ b(x) ≤ β for x ∈ Γ, which guarantees the positivity of the eigenvalues λ i , i ≥ 1. In the cylinder Γ ∞ = Γ 0 × (0, ∞), consider the system (4) defined above. Let us introduce the notation: By the assumption given above, f (x, t) ∈ CL 2,1 (Γ T ) (CL 2,1 (Γ T ) is denoted by the space of functions continuous over the variable t and the norm of the space L 2,1 (Γ T )), and: Let the weak solution for (4) be y (x, t) ∈ V 1,0 (a, Γ t 0 ,∞ ), which is a generalized solution of system (4) in Γ t 0 ,∞ . In this case, we formulate the initial and boundary conditions as follows: We introduce the relation y(x, t) ∈ V 1,0 (a, Γ t 0 ,∞ ) valid for the weak solution of system (4) in the Γ t 0 ,∞ graph under consideration. At the same time, we also formulate the initial and boundary conditions as follows: The weak solution y (x, t) for (4) is denoted by the term 'unperturbed'. On the other hand, for the solution y (x, t), we introduce the definition of 'perturbed'. The previously obtained relations for the generalized solution for the system (4), (5) determine the dependencies y(x, t), y(x, t) in Γ t 0 ,∞ corresponding to the conditions. They also make it possible to determine V 1,0 (a, Γ t 0 ,∞ ) at f (x, t) ∈ CL 2,1 (Γ ∞ ).
We define the stability of its weak solution for (4), (5) in the same way as the Lyapunov stability for the problem under consideration.
It is also possible for the practical needs of management to introduce a definition of the conditions of uniform sustainability. This is necessary for long-term planning of logistics activities. In this case, the unperturbed state of system (4) in the domain Γ t 0 ,∞ under consideration will be determined completely by analogy with Definition 3.
Additionally, for practical applications in long-term planning, the asymptotic and exponential stability can be determined in the framework of (4). Note that the reformulation is permissible due to the linearity of the system (4).
We show that the unperturbed weak solution (state) of system (4) This is a consequence of the linearity of system (4). Let θ(x, t) = y(x, t) − y(x, t), which means θ(x, t) ∈ V 1,0 (a, Γ t 0 ,∞ ), and let θ(x, t) be a generalized solution of the problem for a homogeneous system (4) satisfying the initial and boundary conditions

It follows that problem (4) is uniquely generalized and solvable. Its weak solution
φ i e −λ i t u i (x), φ i = (φ, u i ) and is the limit of the weakly converging sequence θ N N≥1 of its approximations: In this case, the following inequality is performed simultaneously for all N: In this inequality, moving to the limit at N → ∞ , we obtain an upper bound, given e −2λ i t < 1 (i = 1, 2, . . .): for all t ∈ [t 0 , ∞); C * is a constant independent of t. This means that: Fix ε > 0 and take δ = ε C * , then by virtue of the inequality: we obtain the inequality θ(·, t) W 1 (a,Γ) = y(·, t) − y(·, t) W 1 (a,Γ) < ε for arbitrary t > t 0 and, hence, the stability of the unperturbed weak solution (state) of system (4) in Γ T .

Sustainability of Solutions with Discretely Varying Time
First of all, we note that due to the representation (8) of the differential-difference analog of the continuous system (4), the discrete time analogue for the continuously changing time variable t ∈ (0, T] is the set {t = kτ, k = 0, 1, 2, . . . , M}, where τ = T/M. For the differential-difference system (8), (9), stability is defined as the continuous dependence of its solution on the initial data: the functions f τ (k) ∈ L 2 (Γ) (k = 1, 2, . . . , M) and ϕ (x) ∈ L 2 (Γ). Definition 4. A differential-difference system (4), (5) is called stable if the relation is valid for small: for any k = 1, 2, . . . , M; C is a positive constant independent of τ; u(k) 2,Γ is the norm in the space L 2 (Γ); and the norm f τ (k) 2,1,Γ is defined by the formula
For any m = 1, 2, . . . , M: given the evaluation of Theorem 3, we obtain the inequality: Summing up the inequalities for k from 1 to m ≤ M and using the estimates obtained above, we arrive at the inequality: Hence, to the inequality: where the constant C depends only on a * , β, and T.

Connection of Transport Models with Continuously and Discretely Changing Time
Inequalities (14) and (15) make it possible to justify the applicability of the Rote method and to establish a connection between the initial boundary value problem (4), (5) and the differential-difference system (8), (9).
converge uniformly on Γ T to respectively η(x, t), (10) and summing up the obtained identities with respect k and from 1 to M, we arrive at the relation From the obtained relation (17) at M → ∞ , let us pass to the limit with respect to the sequence { u M (x, t)} and obtain identity (6), which means that u(x, t) (the function u(x, t) ∈ W 1,0 0 (a, Γ T )) is a weak solution (state) of system (4) and (5). By virtue of the uniqueness of the solution u(x, t) to problem (4) and (5) on the basic fundamentals and estimate (16), the sequence {u M (x, t)} converges weakly to u(x, t). The theorem is proved.

Note 4.
The results can be used to analyze evolutionary and differential-difference systems with boundary conditions of the second and third types. For this, it is sufficient to set the second or third boundary condition instead of the Dirichlet boundary condition for the boundary value problem (7).

On Computational Aspects of Transport Models with Continuously and Discretely Varying Time
Mathematical models (4), (5) and (8), (9) allow for effective approximation by difference schemes on a two-dimensional grid Γ h,τ T = {(ih, jτ), ih ∈ γ, jτ ∈ [0, T]} and a one-dimensional grid Γ h = {ih} (the points of splitting the edges γ ∈ Γ and the time interval [0, T] are numbered by indices i, j) with approximation errors equal to O(h, τ) and O(h), respectively. In this case, the error in the approximation of systems (4) and (8) on the edges is equal to O(h 2 + τ) and O(h 2 ) and at the nodes is equal to O(h) in both cases; the approximation of boundary conditions (5) and (9) is exact. Numerical analysis of test examples according to schemes for systems (4), (5) and (8), (9) showed almost the same result in the accuracy of calculations, but preference should be given to calculation schemes (8), (9); due to the one-dimensionality of Γ h , the number of arithmetic operations is several orders of magnitude fewer. The restrictions on the use of the presented results are determined by the linearity property of the differential operators of systems (4), (5) and (8), (9).

Discussion
Digital logistics is based on some platform-based business models, so the issues of the implementation of digital business models for sustainable development could be a topic for future research. Companies' information about their socio-environmental actions is important as well as both the rankings of corporate social responsibility [8] and the degree of similarity in international brand valuation rankings that apply to the IT sector [9].
The digital transformation of logistics networks has been driven by different factors, including the approaches of such giant companies as Apple, Alibaba, Facebook, and Google [10] to developing the platform's ecosystem [11] as well as new business models on the basis of the digital twin concept, which was initially implemented in the manufacturing industry [12]. Researchers take into account some different aspects of the digital logistics networks [13][14][15][16] and the features of digital ecosystems [17,18] as well.
The competitiveness of all types of commercial activities using the services of logistics operators, including the concepts of 3PL and 4PL outsourcing, depends on the quality of algorithms for processing data arrays on the movement of goods (taking into account sustainable human resource management [19]). The reader can consider a graph as a prototype of a logistics network. The authors propose the use of the formalisms of initial boundary value problems or differential-difference systems to construct a mathematical model of transfer problems in the class of summable functions, since such systems and the class of functions used help to more accurately take into account the key property of the transferred continuous or discrete media-multiphase media, which involves a lack of continuity in the functions describing the quantitative characteristics of the transferred media. Modern mathematical methods and approaches developed by the authors provide the necessary mathematical tools not only to describe the process on linear network fragments but to obtain adequate mathematical relationships describing the process at the junction nodes of linear network fragments. Additionally, the authors applied the integral functional of the quality of the logistics network [20][21][22][23][24]. The article investigates the boundaries of the stability of an evolutionary parabolic system and a differential-difference system with distributed parameters on a selected graph [25][26][27][28][29]. A computational complexity analysis to demonstrate the computational performance of the proposed solution methodologies could be a topic for future research because of the limitation on the size of the article.

Conclusions
The modern transnational business structure is based on a complex multi-level and constantly evolving transport network. In the middle of the 20th century, simple models for finding optimal routes were developed using linear programming methods. The need to work under conditions of uncertainty and in a dynamic market environment requires us to consider the fundamental difference between discrete and quasi-continuous transportation.
Currently, the theoretical fundamentals of logistics comprise a complex of different methods and models for managing logistics networks. The activity of logistics operators is based on predictive planning and operates with a stream of data on the condition of the cargo and its location in real time. Projecting the movement of material flows onto the logistics network and choosing the optimal supply chains require the creation of effective digital logistics platforms. The core of these systems is formed by an algorithm implemented on mathematical models and the theory of making optimal decisions based on the criteria of economic efficiency. Real logistic activity takes place under the influence of numerous random disturbances. For practical business tasks, it is critically important to analyze the sustainability of the proposed solutions, which not only provides an assessment of investment risks but also allows us to switch to predictive indicators of advanced planning.
In this study, the development of classical approaches to the study of the sustainability of solutions obtained by using mathematical modeling based on systems of differential equations is carried out. The complex topology of logistics networks within which supply chains are built makes it necessary to apply deeper methods in the study of the dynamics of logistics flows. Based on the real conditions of logistics activities, it was assumed that this functionality is not always continuous in terms of movement and time parameters. The complexity of solving the problems of forming effective supply chains has increased under the influence of the trend of the consolidation of business into network structures. The application of the results obtained in this work is not limited to the analysis of the sustainability of solutions. Their use is relevant in assessing investment risks, the sustainability of the business model of outsourcing logistics operators, as well as for gaining competitive advantages.