Preprint
Article

Revolutionary Analysis of the Pointwise Stationary Fluid Flow Approximation model of the non-stationary queue with Pointwise Stationary Fluid Flow Approximation applications to the Internet of Things

This version is not peer-reviewed.

Submitted:

18 February 2024

Posted:

20 February 2024

You are already at the latest version

Abstract
This paper provides the first-time ever analytic modelling for Pointwise Stationary Fluid Flow Approximation(PSFFA) model of the non-stationary
Keywords: 
Time Varying  queue; modelling; simulation; Pointwise Fluid Flow Approximation (PSFFA)
Subject: 
Computer Science and Mathematics  -   Probability and Statistics

1. Introduction

The field of transient/non-stationary analysis has limited literature, which can be categorized into simulation, transient analysis, analysis, and applications techniques. These categories encompass various approaches to studying systems that change over time, including simulations, analysing transient behaviour, and exploring non-stationary phenomena. In certain cases, mathematical transformations can be used to obtain a closed form expression for analysing non-stationary queueing systems. However, evaluating these expressions can be computationally complex. As a result, there has been a focus on numerically determining the transient behaviour of such systems instead of deriving closed form expressions.
The current exposition contributes to solving for first time ever, the longstanding unsolved problem of obtaining the state variable of the time varying M / D / 1 queueing system.
The following flowchart shows how this paper is organized.
Preprints 99206 i001

2. PSFFA

Let f i n ( t ) and f o u t ( t ) serve as the temporal flow in, and fow out ,   respectively. Therefore,
d x ( t ) d t = x . t = f o u t t + f i n t ,   x ( t ) a s t h e s t a t e v a r i a b l e
f o u t t links server utilization, ρ ( t ) and the time-dependent mean service rate, μ t by:
f o u t t = μ t ρ t
For an infite queue waiting space:
f i n t = M e a n a r r i v a l r a t e = λ t
Thus (2.1) rewrites to:
x . t = μ t ρ t + λ t   ,   1 > ρ t = λ t   μ ( t ) > 0
The stable subcase of (2.4) (i.e., x . t = 0), implies:
x = G 1 ( ρ )
The numerical invertibility of G 1 ( ρ ) , yields
ρ = G 1 1 x
Hence,
x . t = μ t ( G 1 1 x ( t ) ) + λ ( t )  
The M / D / 1 queueing system is made of Poisson arrival, one exponential (Poisson) server, FIFO (First-In-First-Out). Thus, M / D / 1 queueing system’s - G 1 (c.f., [1]) reads:
G 1 x = ( x + 1 ) ( x 2 + 1 )
Therefore, the time varying M / D / 1 queueing system’s -PSFFA model reads:
x . = μ ( x + 1 ) ( x 2 + 1 ) + λ
Non-staionary queues’ life example [6] is depicted by Figure 1.

3. Solving the non-stationary queueing system’s PSFFA (c,f., (2.9))

Theorem 3.1 Ismail’s ration, β solves (2.9), with a closed form expression to read as:
x β x + 1 1 1 β 2 = γ e [ ( x + 1 ) 1 β μ ( t ) d t ]   , β = λ t μ t                            
Proof
We have
x . = μ ( x + 1 ) ( x 2 + 1 ) + λ
Let, then x . = y .   c o t h y   c s c h y . Setting, β = λ ( t ) μ ( t ) .   Thus, we have
y .   c o t h y   c s c h y = μ ( 1 + c s c h y ) c o t h y + λ = μ [ 1 + c s c h y ) c o t h y β ]
Therefore, we have
c o t h y   c s c h y d y [ 1 + c s c h y ) c o t h y β ] = μ d t
c o s h y d y [ 1 + s i n h y c o s h y + β s i n h y ] s i n h y = μ d t
2 β ( e 3 y + e y ) d y e 2 y e y β + ( 2 β 1 ) ( e 2 y 1 ) = μ d t
e 2 y e y β + 2 β 1 = 0 e y = a , b
where
a = 1 β + ( 1 β 2 8 β + 4 ) 2 ,   b = 1 β ( 1 β 2 8 β + 4 ) 2
Let
2 β ( e 3 y + e y ) e 2 y e y β + ( 2 β 1 ) ( e 2 y 1 ) = A ( e y a ) + B ( e y b ) + C ( e y 1 ) + D ( e y + 1 )
Hence, it is implied that:
A + B + C + D = 2 β   A = 2 β ( B + C + D )
B = 1 a ( C D ) 2 β + a
Thus,
C = 4 β 2 β + a 1 + 2 β + 2 a a b a + b 2 β + a D 2 β + a + a b a + b 2 β + a
  b A + a B + a b C a b D = 0
Therefore, we have
D = 2 b β + a a b b 1 + 2 β 2 β + a + a 1 a 2 β + a [ 4 β 2 β + a 2 β + a + a b a + b 2 β + a ] ( a b + b 1 + 2 β + 2 a 2 β + a + a 1 a 2 β + a ) + a a b b 1 + 2 β 2 β + a + a 1 a 2 β + a [ 1 + 2 β + 2 a a b a + b 2 β + a 2 β + a + a b a + b 2 β + a
Preprints 99206 i002
A = 2 β 1 + 2 β C 2 β + a 1 + 2 β + 2 a D 2 β + a ,   B = 1 a ( C D ) 2 β + a
This finally solves the complicated mathematical computations to obtain:
A e y 1 a e y + B e y 1 b e y + C e y 1 e y + D e y 1 + e y d y = μ d t
Integrating both sides, implies f o r s o m e n o n n e g a t i v e r e a l c o n s t a n t p a r a m e t e r η ,
| 1 a e y | A | 1 b e y | B | 1 e y | C 1 + e y D = η e μ d t
This transforms to the final required closed form solution:
| 1 a e c s c h 1 ( x ) | A | 1 b e c s c h 1 ( x ) | B | 1 e c s c h 1 ( x ) | C 1 + e c s c h 1 ( x ) D = η e μ d t
By mathematical analysis, it is well known that:
c s c h 1 x = l n 1 + 1 + x 2 x
with the domain of real line with zero removed
Thus, one gets
( | 1 a x 1 + 1 + x 2 | ) A | 1 b x 1 + 1 + x 2 | B ( | 1 x 1 + ( 1 + x 2 ) | ) C 1 + x 1 + ( 1 + x 2 ) D = η e μ d t
Corollary 2:
As x ( t ) 0 , we have
η e μ d t = l i m x ( t ) 0 | 1 a x 1 + 1 + x 2 | A | 1 b x 1 + ( 1 + x 2 ) | B | 1 x 1 + 1 + x 2 | C 1 + x 1 + ( 1 + x 2 ) D = 1
Implying
μ d t = l n γ
Corollary 3
As x ( t ) , we have
η e μ d t = l i m x ( t ) | 1 a x 1 + 1 + x 2 | A | 1 b x 1 + ( 1 + x 2 ) | B | 1 x 1 + 1 + x 2 | C 1 + x 1 + ( 1 + x 2 ) D = l i m x ( t ) | 1 a 1 x + 1 + 1 x 2 | A | 1 b 1 x + ( 1 + ( 1 x ) 2 ) | B | 1 1 1 x + 1 + 1 x 2 | C 1 + 1 1 x + ( 1 + ( 1 x ) 2 ) D = | 1 a | A | 1 b | B | 1 1 1 x + 1 + 1 x 2 | C 1 + 1 1 x + ( 1 + ( 1 x ) 2 ) D
Numerical experiments
Numerical experiment One
Let β = 2   ,   then a = 0.5   ,   b = 0   ,   η = 1 ,   μ t = t , A = 7.923076923 , B = 1.846153846 , C = 6 , D = 0.4615384615 We have
2 [ ( | 1 0.5 x 1 + 1 + x 2 | ) 7.923076923 ( | 1 x 1 + ( 1 + x 2 ) | ) 6 1 + x 1 + ( 1 + x 2 ) 0.4615384615 ] = t
Figure 2 shows a new phenomenon to queuing theorist. The possibility that time will converge to a certain value for sufficiently large number in the time varying M / D / 1 queuing system. This is for an increasing temporal mean service rate.
This shows that as the time varying M / D / 1 queuing system’s state variable becomes sufficiently large, time vanishes.
Numerical experiment Two
Let β = 2 , then a = 0.5 , b = 0 , η = 1 , μ t = 1 t
A = 7.923076923 , B = 1.846153846 , C = 6 , D = 0.4615384615
We have
[ ( | 1 0.5 x 1 + 1 + x 2 | ) 7.923076923 1 x 1 + 1 + x 2 6 1 + x 1 + ( 1 + x 2 ) 0.4615384615 ] = t
Figure 3 visualizes a new phenomenon to queuing theorist. The possibility that time will converge to a certain value for sufficiently large number in the time varying M / D / 1 queuing system. This is for a decreasing temporal mean service rate.

4. Some PSFFA applications to IoT

By describing how vehicles in a platoon use 802.11p communication to exchange messages and change their movement characteristics at intersections, a time-dependent model for assessing the platooning communications’ effectiveness at intersections phase was thoroughly investigated[7]. The model evaluates the effectiveness of platooning communications and addresses potential safety concerns by considering variables including vehicle behaviours, traffic signals, and the changing connectivity among vehicles.
The authors [7] used PSFFA to describe the transmission queue’s dynamic behaviour in platooning communications. They also create models that characterise the continuous backoff freeze and four access categories (ACs) of 802.11p as they relate to the time-dependent access procedure.
For 802.11p communication in platooning situations at junctions, the authors created models. [7] consider continual backoff freezing and use the PSFFA to describe the gearbox queue’s dynamic behaviour. The access process with its four Access Categories (ACs), they also use a z-domain linear model was demonstrated by [7].
The display of time-dependent packet transmission delay and packet delivery ratio of four Acs can be visualized by Figure 4a, Figure 4b, Figure 4c, Figure 5a, Figure 5b and Figure 5c, respectively.
The findings of a study[8] into the variables influencing queue utilisation dynamics on routers in telecommunication networks. The investigation shows that the average queue length reaches a steady state value following a transient process lasting from a few to tens of seconds when utilising a dynamic model, specifically PSFFA. It is advised to calculate the average queue length using steady state estimations only after the transient process has finished and a more precise differential model may be used.
To accurately predict the average queue length while analysing the average queue length and Quality of Service[8] in a network, a dynamic model with a nonlinear differential equation must be used. Only when the transient process has ended, and the length of the transient process is controlled by variables like flow rate, router interface capacity, and service discipline, are steady-state estimations useful for determining the average queue length. A better choice of queuing models and smaller packet sizes can further hasten the average queue length's convergence.

5. Closing remarks with next phase of research

In this work, a challenging topic in queueing theory is examined; more precisely, the underlying queue’s state variable is determined. The article offers a solution to this issue by utilizing a pointwise stationary fluid flow approximation (PSFFA) technique to formulate the non-stationary M / D / 1 queueing system. Future work will concentrate on resolving open research issues and investigating applications of non-stationary queues in other scientific areas. The study also examines the effects of time, combined with and queueing parameters on the underlying queue’s stability dynamics. More fundamentally, some interesting PSFFA applications to IoT are provided. Future work involves further investigation of the impact of   σ t ,   1 > σ t > 0 on the stability of G / M / 1 PSFFA model of t h e n o n s t a t i o n a r y   E k / M /   1   a n d   I P P / M / 1 queueing systems.

References

  1. Zhao, x.; et al. A Queuing Network Model of a Multi-Airport System Based on Point-Wise Stationary Approximation. Aerospace 2022, 9. [Google Scholar] [CrossRef]
  2. Dr Ismail A Mageed. Solving the Open Problem of Finding the Exact Pointwise Stable Fluid Flow Approximation (PSFFA) State Variable of a Non-Stationary M/M/1 Queue with Potential Real-Life PSFFA Applications to Computer Engineering, 28 January 2024, PREPRINT (Version 1) available at Research Square. [CrossRef]
  3. Mageed, I.A.; Zhang, Q. Solving the open problem for 𝐺𝐼/𝑀/1 pointwise stationary fluid flow approximation model (PSFFA) of the non-stationary 𝐷/𝑀/1 queueing system. Electron. J. Comput. Sci. Inf. Technol. 2023, 1, 1–6. [Google Scholar]
  4. Mageed, D.I.A. Upper and Lower Bounds of the State Variable of 𝑀/𝐺/1 PSFFA Model of the Non-Stationary M/Ek/1 Queueing System. Preprints 2024, 2024012243. [Google Scholar] [CrossRef]
  5. Hu, L.; Zhao, B.; Zhu, J.; Jiang, Y. Two time-varying and state-dependent fluid queuing models for traffic circulation systems. Eur. J. Oper. Res. 2019, 275, 997–1019. [Google Scholar] [CrossRef]
  6. Mageed, D.I.A. Ismail’s Ratio Conquers New Horizons the Non-stationary M/G/1 Queue’s State Variable Closed Form Expression. Preprints 2024, 2024020548. [Google Scholar] [CrossRef]
  7. Wu, Q.; Zhao, Y.; Fan, Q. Time-dependent performance modeling for platooning communications at intersection. IEEE Internet Things J. 2022, 9, 18500–13. [Google Scholar] [CrossRef]
  8. Yeremenko, O.; Lebedenko, T.; Vavenko, T.; Semenyaka, M. Investigation of queue utilization on network routers by the use of dynamic models. In 2015 Second International Scientific-Practical Conference Problems of Info communications Science and Technology (PIC S&T.
Figure 1.
Figure 1.
Preprints 99206 g001
Figure 2.
Figure 2.
Preprints 99206 g002
Figure 3.
Figure 3.
Preprints 99206 g003
Figure 4. The packet transmission delay that is time dependent. A left turn, a straight shot, or a right turn, respectively[7].
Figure 4. The packet transmission delay that is time dependent. A left turn, a straight shot, or a right turn, respectively[7].
Preprints 99206 g004
Figure 5. The proportion of time-dependent packet deliveries. A left turn, a straight shot, or a right turn, respectively[7].
Figure 5. The proportion of time-dependent packet deliveries. A left turn, a straight shot, or a right turn, respectively[7].
Preprints 99206 g005
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.
Copyright: This open access article is published under a Creative Commons CC BY 4.0 license, which permit the free download, distribution, and reuse, provided that the author and preprint are cited in any reuse.

Downloads

89

Views

26

Comments

0

Subscription

Notify me about updates to this article or when a peer-reviewed version is published.

Email

Prerpints.org logo

Preprints.org is a free preprint server supported by MDPI in Basel, Switzerland.

Subscribe

© 2025 MDPI (Basel, Switzerland) unless otherwise stated