1. Introduction
The customer behavior analysis of queueing systems has received considerable attention due to the widespread application of random service systems. Most references assume that the server is always available, but perfectly reliable servers are virtually nonexistent. Economou and Kanta [Reference Economou and Kanta3] considered the equilibrium behavior in the observable queues with breakdowns and repairs. Subsequently, Li et al. [Reference Li, Wang and Zhang7] extended the research results in [Reference Economou and Kanta3] and derived the joining strategies for customers in unobservable cases. Wang and Zhang [Reference Wang and Zhang10] studied an unreliable M/M/1 queue with delayed repairs in observable cases. Inspired by Wang and Zhang [Reference Wang and Zhang10], Yu et al. [Reference Yu, Liu and Wu15] explored the same settings in unobservable cases. Besides, Yang et al. [Reference Yang, Wang and Zhang14] analyzed the Geo/Geo/1 queue with fault characteristics and determined the Nash equilibrium and the social optimum. Chen and Zhou [Reference Chen and Zhou2] incorporated the setup times into a repairable queue and derived the mixed strategies of customers. Moreover, Li et al. [Reference Li, Wang and Zhang5] obtained the equilibrium balking strategies and socially optimal behavior in an M/M/1 queue with partial breakdowns and immediate repairs, where the systems can continue the service at a lower rate during the breakdowns. Yu et al. [Reference Yu, Liu and Wu16] compensated the game analysis in [Reference Li, Wang and Zhang5] by studying the corresponding partially observable cases. Xu and Xu [Reference Xu and Xu12] further discussed the individual Nash equilibrium in an M/M/1 queue with partial failures and delayed repairs. Considering more complex service systems, Zhang and Xu [Reference Zhang and Xu18] made an economic analysis in a queue with setup times and partial failures. Aghsami and Jolai [Reference Aghsami and Jolai1] studied the equilibrium balking strategies in the single-server Markovian queue with partial breakdowns and interruptible setup policy. The economic analysis of queueing systems with breakdowns has been widely applied in the inventory management, communication networks, and supply-chain systems.
In a complex queueing system, the interarrival time of customers is getting smaller, and the service time is getting shorter. Discrete customers can be considered as the continuous fluid. There is an emerging tendency to study the fluid models strategically due to the limitation and the complexity of discrete queueing systems. Economou and Manou [Reference Economou and Manou4] considered the fluid model with two service modes and derived the individual equilibrium thresholds and socially optimal strategies. Subsequently, Wang and Xu [Reference Wang and Xu11] studied the strategic behavior in a fluid vacation queue. Recently, Xu and Wang [Reference Xu and Wang13] explored an on-off fluid queue and characterized the mixed equilibrium strategies in the unobservable cases. Logothetis et al. [Reference Logothetis, Manou and Economou9] also analyzed a fluid on-off model and studied the equilibrium strategies of customers under conditions of reneging and no-reneging, respectively. However, none of the above literature considers the entrance fee problems, which can eliminate the gap between the individual equilibrium and socially optimal strategies.
In real life, parallel customers who have the same service levels are widespread in public services, such as volunteer clinics for specific populations, psychological counseling, and legal advice. Li et al. [Reference Li, Xu and Ni6] analyzed the service policy in an M/M/1 queue with two types of customers and different waiting costs. Zhang and Xu [Reference Zhang and Xu17] incorporated two types of parallel customers into a fault queueing system and obtained the equilibrium balking strategy in the observable cases. Based on Zhang and Xu [Reference Zhang and Xu17], Liu et al. [Reference Liu, Xu, Wang and Yue8] studied a repairable fluid queue with two types of parallel customers from an economic viewpoint. However, in many real-life situations, the server may not stop working completely in the maintenance period due to the durability and complexity of modern operating systems and system failures provide lower service rates. For example, a running computer system may slow down due to aging hardware or virus attacks. Therefore, incorporating the incomplete fault schedule into the fluid queue is reasonable and valuable, and this work compensates the equilibrium game analysis in [Reference Liu, Xu, Wang and Yue8] by studying the incomplete breakdowns.
The proposed model has broad applications and contributions in many fields, one of which is wireless communication technology. With the rapid development of telecommunication networks, wireless resources are shared among various services, and suppliers divide the information flowing into the processor into data, voice, and video categories to provide communication services. The unreliable router can be interrupted when transmitting information due to hardware failures and needs a random amount of time for recovery. During the maintenance periods, the router reduces the efficiency and can still transmit the information fluid at a lower rate. In this scenario, the transmitted information, the router, and the breakdown correspond to the arriving fluids, the buffer, and the incomplete fault, respectively, in the queueing terminology. The single service and pricing standard in the queueing system cannot satisfy the requirements of the growing variety of scenarios. Therefore, considering the equilibrium behavior of two or even more types of customers in an unreliable fluid queue is significant and necessary.
Motivated and inspired by the above practical cases, this paper makes an economic analysis for two types of the fluid in an unreliable server under the fully observable case. Different from [Reference Liu, Xu, Wang and Yue8], this paper assumes that the service can provide relatively low service rates after failure. To the best of the author’s knowledge, no research has been found in the literature on conducting equilibrium strategies for fluid queues with incomplete fault and parallel customers. The individual equilibrium strategy can be obtained based on the non-cooperative game theory, and an exponential utility function is constructed to compute the expected social benefit. The reasonable entrance fee strategy can dynamically regulate system parameters and optimize the expected social benefits. Besides, the sensitivity analysis between the globally optimal thresholds and the expected social benefit is illustrated by several numerical examples. The research results provide a theoretical basis and valid suggestions for signal transmissions in network systems, which can obtain more stable response performance in the practical production life.
This article is organized as follows. Section 2 gives a detailed description of the fluid model. Sections 3 and 4 are devoted to the individual balking strategy and social equilibrium analysis in the fully observable case, respectively. Furthermore, the effects of the expected social benefit on globally optimal thresholds are illustrated by numerical examples in Section 4. Section 5 presents an entrance fee policy and analyzes it theoretically and numerically. Finally, we briefly conclude the paper in Section 6.
2. Model description
Assume that two types of parallel customers flow into the buffer according to the exponential distribution with rates 
${{\lambda }_{1}}$ and 
${{\lambda }_{2}}$ 
$\left( {{\lambda }_{1}} \lt {{\lambda }_{2}} \right)$, respectively, and the arrival processes are independent of each other. 
$I\left( t \right)$ is the state of the buffer at time t, and 
$I\left( t \right)=0,1$ indicates that the buffer is in maintenance period and normal working period, respectively. The buffer has an exponentially distributed lifetime with parameter 
${\mathrm{\theta}}$ and the output rate of the buffer in the normal working state is 
${{\mu }_{1}}$. When the system breaks down, the buffer immediately enters the maintenance state and provides services for the fluid at a lower rate of 
${{\mu }_{0}}\ (\mu_0 \lt \mu_1$). The maintenance time follows the exponential distribution with parameter 
${\mathrm{\xi}}$. Let 
$X\left( t \right)$ be the fluid level at time t, and the net input rate of the buffer can be expressed as:
\begin{equation*}\frac{{{\rm d}X(t)}}{{{\rm d}t}} = \left\{{\begin{array}{*{20}{l}}
{{\lambda _1} + {\lambda _2} - {\mu_1},}&{I(t) = 1,X(t) \gt 0,}\\
{{\lambda _1} + {\lambda _2} - {\mu_0},}&{I(t) = 0,X(t) \gt 0,}\\
{\max \{{\lambda _1} + {\lambda _2} - {\mu_1},0\} ,}&{I(t) = 1,X(t) = 0.}
\end{array}} \right.\end{equation*} The condition 
${{\lambda }_{j}} \gt {{\mu }_{0}},j=1,2$ ensures that the buffer can reach the steady-state.
Suppose that the unit of the type 
$j(j=1,2)$ fluid receives a reward of 
${{R}_{j}}$ when it flows out of the buffer and costs 
${{C}}_{j}$ per unit time during its sojourn time, and the inequality 
${{R}_{j}} \gt \frac{{{C}_{j}}}{\mu_0}$ is satisfied hereinafter to ensure that the fluid prefers to enter when the buffer is empty.
In this paper, the fluid level X(t) and the buffer’s state I(t) can be observed when the fluid arrives the buffer, and the decisions are irreversible, that is, neither retrials of balking customers or reneging of entering customers are permitted. When the fluid reaches the buffer at time t, we assume that the observed buffer state is 
$\left( X\left( t \right),I\left( t \right) \right)=\left( x,i \right)$.
3. Individual equilibrium strategy
In the financial and insurance industry, the exponential utility function is more practical and significant than the linear utility function. The expected net benefit per unit time for the type j fluid after its service when the buffer is in state 
$\left( x,i \right)$ can be denoted as:
\begin{equation*}{B_j}\left( {x,i} \right) = {R_j} - {C_j}E\left[ {{{\rm e}^{\alpha {S_j}\left( {x,i} \right)}}} \right],\end{equation*} where 
${{S}_{j}}\left( x,i \right), j=1,2$ represents the expected sojourn time for the type j fluid. 
${{{\rm e}^{\alpha {S_j}\left( {x,i} \right)}}}$ is integrable and α is a constant, and 
$0 \lt \alpha \lt \theta, 0 \lt \alpha \lt \xi $.
The pre-condition of the fluid flowing into the buffer is that the expected net benefit 
${{B}_{j}}\left( x,i \right)$ is positive, then we get the following result.
Theorem 3.1. In the fully observable case, there exists a pair of thresholds 
$\left( {{x^*}_{j}}\left( 0 \right),{{x^*}_{j}}\left( 1 \right) \right)$ for type 
$j\,(\,j=1,2)$ customers in the fluid model with two types of parallel customers and incomplete fault, and the equilibrium strategies 
${{x^*}_{j}}\left( 0 \right)$ and 
${{x^*}_{j}}\left( 1 \right)$, respectively, are the unique roots of the equations described below.
\begin{equation*}R_j-C_j\frac{{\left[ {\alpha ({\mu_0} - {\mu_1}) - \theta {\mu_0}} \right]{{\rm e}^{\frac{{\alpha - \xi }}{{{\mu_0}}}x}} + \xi {\mu_1}{{\rm e}^{\frac{{\alpha - \theta }}{{{\mu_1}}}x}}}}{{\alpha ({\mu_0} - {\mu_1}) - \theta {\mu_0} + \xi {\mu_1}}}=0,\end{equation*}
\begin{equation*}R_j-C_j \frac{{\left[ {\alpha ({\mu_0} - {\mu_1}) + \xi {\mu_1}} \right]{{\rm e}^{\frac{{\alpha - \theta }}{{{\mu_1}}}x}} + \theta {\mu_0}{{\rm e}^{\frac{{\alpha - \xi }}{{{\mu_0}}}x}}}}{{\alpha ({\mu_0} - {\mu_1}) - \theta {\mu_0} + \xi {\mu_1}}}=0.\end{equation*} The equilibrium strategy of the fluid is defined as “while arriving at time t, observe the system is in state 
$(X(t),I(t))$, the type j customers enter if 
$X(t)\le {{x^*}_{j}}\left( I(t) \right)$ and balk otherwise.”
Proof. Assume that 
${{T}_{i}}$ is the remaining sojourn time of the buffer in the state 
$i,i=0,1$ when the fluid arrives, then
\begin{equation}
{S_j}(x,i) = \left\{{\begin{array}{*{20}{l}}
{\frac{x}{{{\mu_i}}},}&{{T_i} \ge \frac{x}{{{\mu_i}}},}\\
{{T_i} + {S_j}(x - \mu {T_i},1 - i),}&{{T_i} \lt \frac{x}{{{\mu_i}}}.}
\end{array}} \right.
\end{equation}Using the conditional expectation formula, we have:
\begin{equation*}E\left( {{{\rm e}^{\alpha {S_j}(x,0)}}} \right) = {{\rm e}^{(\alpha - \xi )\frac{x}{{{\mu_0}}}}} + \int_0^{\frac{x}{{{\mu_0}}}} {E\left( {{{\rm e}^{\alpha {S_j}(x - {\mu_0}t,1)}}} \right)\xi {{\rm e}^{(\alpha - \xi )t}}{\rm d}t}, \end{equation*}
\begin{equation*}E\left( {{{\rm e}^{\alpha {S_j}(x,1)}}} \right) = {{\rm e}^{(\alpha - \theta )\frac{x}{{{\mu_1}}}}} + \int_0^{\frac{x}{{{\mu_1}}}} {E\left( {{{\rm e}^{\alpha {S_j}(x - {\mu_1}t,0)}}} \right)\theta {{\rm e}^{(\alpha - \theta )t}}{\rm d}t}. \end{equation*}After manipulating, we obtain the differential equations,
\begin{equation}
\frac{{\xi - \alpha }}{{{\mu_0}}}E\left( {{{\rm e}^{\alpha {S_j}(x,0)}}} \right) + \frac{{{{\rm d}}E({{\rm e}^{\alpha {S_j}(x,0)}})}}{{{{\rm d}}x}} = \frac{\xi }{{{\mu_0}}}E\left( {{{\rm e}^{\alpha {S_j}(x,1)}}} \right),\end{equation}
\begin{equation}\frac{{\theta - \alpha }}{{{\mu_1}}}E\left( {{{\rm e}^{\alpha {S_j}(x,1)}}} \right) + \frac{{{{\rm d}}E({{\rm e}^{\alpha {S_j}(x,1)}})}}{{{{\rm d}}x}} = \frac{\theta }{{{\mu_1}}}E\left( {{{\rm e}^{\alpha {S_j}(x,0)}}} \right),\end{equation} with the boundary condition 
$E\left( {{{\rm e}}^{\alpha {{S}_{j}}(0,i)}} \right)=1$.
From formulas (2) and (3), the exponential form of the expected sojourn time in the buffer can be obtained as:
\begin{equation*}E({{\rm e}^{\alpha {S_j}(x,0)}}) = \frac{{\left[ {\alpha ({\mu_0} - {\mu_1}) - \theta {\mu_0}} \right]{{\rm e}^{\frac{{\alpha - \xi }}{{{\mu_0}}}x}} + \xi {\mu_1}{{\rm e}^{\frac{{\alpha - \theta }}{{{\mu_1}}}x}}}}{{\alpha ({\mu_0} - {\mu_1}) - \theta {\mu_0} + \xi {\mu_1}}}, i=0,\end{equation*}
\begin{equation*}E({{\rm e}^{\alpha {S_j}(x,1)}}) = \frac{{\left[ {\alpha ({\mu_0} - {\mu_1}) + \xi {\mu_1}} \right]{{\rm e}^{\frac{{\alpha - \theta }}{{{\mu_1}}}x}} + \theta {\mu_0}{{\rm e}^{\frac{{\alpha - \xi }}{{{\mu_0}}}x}}}}{{\alpha ({\mu_0} - {\mu_1}) - \theta {\mu_0} + \xi {\mu_1}}}, i=1.\end{equation*} Obviously, 
$E({{\rm e}^{\alpha {S_j}(x,i)}}),i=0,1$ is a monotonically increasing function with respect to x, then the unique roots 
$ {{x^*}_{j}}\left( 0 \right)$ and 
${{x^*}_{j}}\left( 1 \right) $ for type 
$j\,(j=1,2)$ customers are given in Theorem 3.1.
4. Social equilibrium analysis
This section studies the steady-state probability distribution of the fluid level and the socially optimal strategy in a global optimization problem. The fluid considers its own interests first when encountering a joining-balking problem, then the buffer will be over-crowded and cannot achieve the global optimum. Assume that all customers follow the socially optimal threshold strategy 
$\left( {{x}_{ej}}\left( 0 \right),{{x}_{ej}}\left( 1 \right) \right),j=1,2,$ and the fluid flows into the buffer if 
$X\left( t \right) \lt {{x}_{ej}}\left( i \right)$ when the system is in state (x, i), otherwise the fluid balks. Besides, the fluid prefers to enter when the buffer is in a normal working period because the expected sojourn time of the fluid in state 1 is always less than the expected sojourn time in state 0 with the same fluid level. Therefore, the inequalities 
${{x}_{e1}} \left( 0 \right) \lt {{x}_{e1}}\left( 1 \right)$ and 
${{x}_{e2}} \left( 0 \right) \lt {{x}_{e2}}\left( 1 \right)$ can be obtained. Then, we suppose that the thresholds satisfy 
${{x}_{e1}} \left( 0 \right) \lt {{x}_{e2}}\left( 0 \right) \lt {{x}_{e1}}\left( 1 \right) \lt {{x}_{e2}}\left( 1 \right),$ and the others can be calculated based on the similar calculation method.
The overall expected social benefit in unit time is:
\begin{equation*}
\begin{aligned}
{B}\left( {{x_{e1}}\left( 0 \right),{x_{e2}}\left( 0 \right),{x_{e1}}\left( 1 \right),{x_{e2}}\left( 1 \right)} \right) &= \sum\limits_{j = 1}^2 {{B_j}({x_{ej}}(0),{x_{ej}}(1))}= \sum\limits_{j = 1}^2 {\left( {{\lambda _{ej}}{R_j} - {C_j}E\left( {{{\rm e}^{\alpha {X_j}}}} \right)} \right)},\end{aligned}
\end{equation*} where 
${{B}_{j}}({{x}_{ej}}(0),{{x}_{ej}}(1))$ and 
${{\lambda }_{ej}}$ are the expected social benefit per unit time and the effective arrival rate for type j customers, respectively. The effective arrival rate λej can be calculated from the arrival rate and the entrance probability, and 
${{X}_{j}}$ is a random variable that denotes the fluid level of type j customers.
4.1. Steady-state probability distribution
The stationary probability distribution of the fluid level in state i can be defined as:
\begin{equation*}{F_i}(x) = \mathop {\lim }\limits_{t \to \infty }F_i(t,x)= p\{X(t) \le x,I(t) = i\},\quad x \ge 0,i = 0,1.\end{equation*} According to the alternating renewal process, we can get the steady-state distribution of the buffer is 
${{\pi }_{0}}=\frac{\theta }{\theta +\xi },{{\pi }_{1}}=\frac{\xi }{\theta +\xi }$.
Theorem 4.1. In the fully observable fluid queue with two types of parallel customers and incomplete fault, if all customers follow the threshold strategy 
$\left( {{x}_{ej}}\left( 0 \right),{{x}_{ej}}\left( 1 \right) \right)$ and 
${{x}_{e1}} \left( 0 \right) \lt {{x}_{e2}}\left( 0 \right) \lt {{x}_{e1}}\left( 1 \right) \lt {{x}_{e2}}\left( 1 \right)$, the steady-state probability distribution of the fluid level is as follows.
Case I. 
${{\lambda }_{2}} \gt {{\mu }_{1}}$.
\begin{equation}
{F_0}(x) = \left\{{\begin{array}{*{20}{l}}
0,&{x \lt {x_{e2}}(0)},\\
{{\pi _0}\frac{{ \xi (q{\lambda _1} + {\lambda _2} - {\mu_1}){{\rm e}^{- \left( {\frac{\theta }{{q{\lambda _1} + {\lambda _2} - {\mu_1}}} - \frac{\xi }{{{\mu_0}}}} \right)\left( {x - {x_{e2}}(0)} \right)}- \theta{\mu_0}}}}{{\xi (q{\lambda _1} + {\lambda _2} - {\mu_1}){{\rm e}^{- \left( {\frac{\theta }{{q{\lambda _1} + {\lambda _2} - {\mu_1}}} - \frac{\xi }{{{\mu_0}}}} \right)\left( {{x_{e2}}(1) - {x_{e2}}(0)} \right)}} - \theta {\mu_0}}}},&{{x_{e2}}(0) \le x \le {x_{e2}}(1)},\\
{{\pi _0}},&{x \ge {x_{e2}}(1)}.
\end{array}} \right.\end{equation}
\begin{equation}{F_1}(x) = \left\{{\begin{array}{*{20}{l}}
0,&{x \le {x_{e2}}(0)},\\
{{\pi _1}\frac{{\theta {\mu_0}\left[ {{{\rm e}^{- \left( {\frac{\theta }{{q{\lambda _1} + {\lambda _2} - {\mu_1}}} - \frac{\xi }{{{\mu_0}}}} \right)\left( {x - {x_{e2}}(0)} \right)}} - 1} \right]}}{{\xi \left( {q{\lambda _1} + {\lambda _2} - {\mu_1}} \right){{\rm e}^{- \left( {\frac{\theta }{{q{\lambda _1} + {\lambda _2} - {\mu_1}}} - \frac{\xi }{{{\mu_0}}}} \right)\left( {{x_{e2}}(1) - {x_{e2}}(0)} \right)}} - \theta {\mu_0}}}},&{{x_{e2}}(0) \le x \le {x_{e2}}(1)},\\
{{\pi _1}},&{x \gt {x_{e2}}(1)}.
\end{array}} \right.\end{equation} Case II. 
${{\lambda }_{1}}+{{\lambda }_{2}} \gt {{\mu }_{1}}$ and 
${{\lambda }_{2}}\le {{\mu }_{1}}$.
\begin{equation}{F_0}(x) = \left\{{\begin{array}{*{20}{l}}
0,&{x \lt {x_{e2}}(0)},\\
{{\pi _0}\frac{{\xi ({\lambda _1} + {\lambda _2} - {\mu_1}){{\rm e}^{- \left( {\frac{\theta }{{{\lambda _1} + {\lambda _2} - {\mu_1}}} - \frac{\xi }{{{\mu_0}}}} \right)\left( {x - {x_{e2}}(0)} \right)}- \theta{\mu_0}}}}{{\xi ({\lambda _1} + {\lambda _2} - {\mu_1}){{\rm e}^{- \left( {\frac{\theta }{{{\lambda _1} + {\lambda _2} - {\mu_1}}} - \frac{\xi }{{{\mu_0}}}} \right)\left( {{x_{e1}}(1) - {x_{e2}}(0)} \right)}} - \theta {\mu_0}}}},&{{x_{e2}}(0) \le x \le {x_{e1}}(1)},\\
{{\pi _0}},&{x \ge {x_{e1}}(1)}.
\end{array}} \right.\end{equation}
\begin{equation}{F_1}(x) = \left\{{\begin{array}{*{20}{l}}
0,&{x \le {x_{e2}}(0)},\\
{{\pi _1}\frac{{\theta {\mu_0}\left[ {{{\rm e}^{- \left( {\frac{\theta }{{{\lambda _1} + {\lambda _2} - {\mu_1}}} - \frac{\xi }{{{\mu_0}}}} \right)\left( {x - {x_{e2}}(0)} \right)}} - 1} \right]}}{{\xi \left( {{\lambda _1} + {\lambda _2} - {\mu_1}} \right){{\rm e}^{- \left( {\frac{\theta }{{{\lambda _1} + {\lambda _2} - {\mu_1}}} - \frac{\xi }{{{\mu_0}}}} \right)\left( {{x_{e1}}(1) - {x_{e2}}(0)} \right)}} - \theta {\mu_0}}}},&{{x_{e2}}(0) \le x \le {x_{e1}}(1)},\\
{{\pi _1}},&{x \gt {x_{e1}}(1)}.
\end{array}} \right.\end{equation} Case III. 
${\lambda _{\rm{1}}}{\rm{+ }}{\lambda _{\rm{2}}}{\rm{= }}{\mu_{\rm{1}}}$. The fluid level is stabilized at 
${{x}_{e2}}(0)$.
Case IV. 
${\mu_{\rm{0}}} \lt {\lambda _1} + {\lambda _2} \lt {\mu_1}$.
\begin{equation}{F_0}(x) = \left\{{\begin{array}{*{20}{l}}
0,&{x \le 0},\\
{{\pi _0}\frac{{\xi \left( {{\lambda _1} + {\lambda _2} - {\mu_1}} \right)\left[ {1 - {{\rm e}^{- \left( {\frac{\xi }{{p{\lambda _1} + {\lambda _2} - {\mu_0}}} + \frac{\theta }{{{\lambda _1} + {\lambda _2} - {\mu_1}}}} \right)x}}} \right]}}{{\theta \left( {p{\lambda _1} + {\lambda _2} - {\mu_0}} \right){{\rm e}^{- \left( {\frac{\xi }{{p{\lambda _1} + {\lambda _2} - {\mu_0}}} + \frac{\theta }{{{\lambda _1} + {\lambda _2} - {\mu_1}}}} \right){x_{e2}}(0)}} + \xi \left( {{\lambda _1} + {\lambda _2} - {\mu_1}} \right)}}},&{0 \le x \le {x_{e2}}(0)},\\
{{\pi _0}},&{x \gt {x_{e2}}(0)}.
\end{array}} \right.\end{equation}
\begin{equation}{F_1}(x) = \left\{{\begin{array}{*{20}{l}}
0,&{x \lt 0},\\
{{\pi _1}\frac{{\xi ({\lambda _1} + {\lambda _2} - {\mu_1}) + (p{\lambda _1} + {\lambda _2} - {\mu_0})\theta {{\rm e}^{- \left( {\frac{\xi }{{p{\lambda _1} + {\lambda _2} - {\mu_0}}} + \frac{\theta }{{{\lambda _1} + {\lambda _2} - {\mu_1}}}} \right)x}}}}{{\theta (p{\lambda _1} + {\lambda _2} - {\mu_0}){{\rm e}^{- \left( {\frac{\xi }{{p{\lambda _1} + {\lambda _2} - {\mu_0}}} + \frac{\theta }{{{\lambda _1} + {\lambda _2} - {\mu_1}}}} \right){x_{e2}}(0)}} + \xi ({\lambda _1} + {\lambda _2} - {\mu_1})}}},&{0 \le x \le {x_{e2}}(0)},\\
{{\pi _1}},&{x \ge {x_{e2}}(0)}.
\end{array}} \right.\end{equation}Proof. By comparing the arrival rate and the outflow rate, the steady-state probability distribution of the fluid level can be discussed in the following four cases.
Case I. 
${\lambda _2} \gt \mu_{1}$.
When the buffer is in state 1, the fluid in the buffer increases at the rate of 
${{\lambda }_{1}}+{{\lambda}_{2}}-{\mu_{1}}$ until the fluid level reaches the threshold 
${{x}_{e1}}(1)$. When the fluid level exceeds the optimal threshold 
${{x}_{e1}}(1)$, only type 2 customers flow in and the fluid level increases to the threshold 
${{x}_{e2}}(1)$ at the rate of 
${{\lambda }_{2}}-{\mu_{1}}$. If the buffer is still in state 1, the fluid will enter the system with probability 
$\frac{{\mu_{1}}}{{{\lambda}_{2}}}$ to ensure that the fluid level is stabilized at 
${{x}_{e2}}(1)$. When the buffer is in state 0, the fluid decreases at the rate of 
${\mu_{0}}$. The system repeats in this way, and the fluid level fluctuates within 
$\left[{{x}_{e2}}(0),{{x}_{e2}}(1) \right]$.
Considering the fluid level in a tiny time interval, we have
\begin{equation*}
\begin{aligned}
&{F_0}(t + \Delta t,x) = {F_0}(t,x +{\mu_0}\Delta t){{\rm{e}}^{\xi \Delta t}} + {F_1}(t,x - \left( {q{\lambda _1} + {\lambda _2} - {\mu_1}} \right)\Delta t)(1 - {{\rm{e}}^{\theta\Delta t}}),\\
&{F_1}(t + \Delta t,x) = {F_1}(t,x - \left( {q{\lambda _1} + {\lambda _2} - {\mu_1}} \right)\Delta t){{\rm{e}}^{\theta\Delta t}} + {F_0}(t,x +{\mu_0}\Delta t)(1 - {{\rm{e}}^{\xi\Delta t}}),
\end{aligned}
\end{equation*}Substituting the exponential form of the Taylor formulas with Peano terms, we obtain:
\begin{equation*}
\begin{aligned}
&{F_0}(t + \Delta t,x) = {F_0}(t,x +{\mu_0}\Delta t)(1 - \xi\Delta t) + {F_1}(t,x - \left( {q{\lambda _1} + {\lambda _2} - {\mu_1}} \right)\Delta t)\theta\Delta t + o(\Delta t),\\
&{F_1}(t + \Delta t,x) = {F_1}(t,x - \left( {q{\lambda _1} + {\lambda _2} - {\mu_1}} \right)\Delta t)(1 - \theta\Delta t) + {F_0}(t,x + {\mu_0}\Delta t)\xi\Delta t + o(\Delta t),
\end{aligned}
\end{equation*} Dividing both sides by 
$\Delta t$ and making 
$\Delta t \to 0$, then we can get:
\begin{equation*}\begin{array}{l}
\left\{
\begin{aligned}
&\frac{{\partial {F_0}(t,x)}}{{\partial t}} -{\mu_0} \frac{{\partial {F_0}(t,x)}}{{\partial x}} = -\xi{F_0}(t,x) + \theta{F_1}(t,x),\\
&\frac{{\partial {F_1}(t,x)}}{{\partial t}} + \left( {q{\lambda _1} + {\lambda _2} - {\mu_1}} \right)\frac{{\partial {F_1}(t,x)}}{{\partial x}} = - \theta{F_1}(t,x) + \xi{F_0}(t,x),
\end{aligned}
\right.
\end{array}\end{equation*} When the system is stable, there exists 
$\mathop {\lim }\limits_{t \to \infty } \frac{{\partial {F_i}(t,x)}}{{\partial t}} = 0,{i} = 0,1.$
The ordinary differential equations can be constructed as:
\begin{equation}\begin{array}{l}
\left\{
\begin{aligned}
&\ -{\mu_0} \frac{{{\rm d}{F_0}(x)}}{{{\rm d}x}} = -\xi {F_0}(x) + \theta {F_1}(x),\\
&\left( {q{\lambda _1} + {\lambda _2} - {\mu_1}} \right)\frac{{{\rm d}{F_1}(x)}}{{{\rm d}x}} = \xi {F_0}(x) - \theta {F_1}(x),
\end{aligned}
\right.
\end{array}\end{equation}with boundary conditions
\begin{equation*}
{F_0}({x_{e2}}(1)) = {\pi _0},\ {F_1}({x_{e2}}(0)) = 0,
\end{equation*} where 
$q=\frac{{{x}_{e1}}(1)-{{x}_{e2}}(0)}{{{x}_{e2}}(1)-{{x}_{e2}}(0)}, 0\le q\le 1$ is the probability that the fluid level in the system is less than the threshold 
${{x}_{e1}}(1)$ in state 1.
The formulas (4) and (5) can be obtained by solving the differential equation (10). The point masses and the probability density functions in states 0 and 1, respectively, are:
\begin{equation*}{P_0}({x_{e2}}(0))=F_{0}({x_{e2}}(0)),\
{P_1}({x_{e2}}(1))=\pi_{1}-F_{1}({x_{e2}}(1)),\end{equation*}
\begin{equation*}{f_0}(x) = \frac{{{\pi _0}\xi (q{\lambda _1} + {\lambda _2} - {\mu_1})\left( {\frac{\xi }{{{\mu_0}}} - \frac{\theta }{{q{\lambda _1} + {\lambda _2} - {\mu_1}}}} \right){{\rm e}^{- \left( {\frac{\theta }{{q{\lambda _1} + {\lambda _2} - {\mu_1}}} - \frac{\xi }{{{\mu_0}}}} \right)\left( {x - {x_{e2}}(0)} \right)}}}}{{\xi (q{\lambda _1} + {\lambda _2} - {\mu_1}){{\rm e}^{- \left( {\frac{\theta }{{q{\lambda _1} + {\lambda _2} - {\mu_1}}} - \frac{\xi }{{{\mu_0}}}} \right)\left( {{x_{e2}}(1) - {x_{e2}}(0)} \right)}} - \theta {\mu_0}}},\end{equation*}
\begin{equation*}{f_1}(x) = \frac{{{\pi _1}\theta {\mu_0}\left( {\frac{\xi }{{{\mu_0}}} - \frac{\theta }{{q{\lambda _1} + {\lambda _2} - {\mu_1}}}} \right){{\rm e}^{- \left( {\frac{\theta }{{q{\lambda _1} + {\lambda _2} - {\mu_1}}} - \frac{\xi }{{{\mu_0}}}} \right)\left( {x - {x_{e2}}(0)} \right)}}}}{{\xi \left( {q{\lambda _1} + {\lambda _2} - {\mu_1}} \right){{\rm e}^{- \left( {\frac{\theta }{{q{\lambda _1} + {\lambda _2} - {\mu_1}}} - \frac{\xi }{{{\mu_0}}}} \right)\left( {{x_{e2}}(1) - {x_{e2}}(0)} \right)}} - \theta {\mu_0}}}.\end{equation*} Case II. 
${{\lambda }_{1}}+{{\lambda }_{2}} \gt {{\mu }_{1}}$ and 
${{\lambda }_{2}}\le {{\mu }_{1}}$.
When the buffer is in state 1, the fluid in the buffer increases at the rate of 
${{\lambda }_{1}}+{{\lambda }_{2}}-{{\mu }_{1}}$ until the fluid level reaches the threshold 
${{x}_{e1}}(1)$, then the fluid level remains unaltered. When the buffer is in state 0, the fluid decreases at the rate of 
${{\mu }_{0}}$. The system repeats in this way, and the fluid level fluctuates within 
$\left[ {{x}_{e2}}(0),{{x}_{e1}}(1) \right]$.
Considering the fluid level in a tiny time interval, we have:
\begin{equation*}
\begin{aligned}
&{F_0}(t + \Delta t,x) = {F_0}(t,x +{\mu_0}\Delta t){{\rm{e}}^{\xi \Delta t}} + {F_1}(t,x - \left( {{\lambda _1} + {\lambda _2} - {\mu_1}} \right)\Delta t)(1 - {{\rm{e}}^{\theta\Delta t}}),\\
&{F_1}(t + \Delta t,x) = {F_1}(t,x - \left( {{\lambda _1} + {\lambda _2} - {\mu_1}} \right)\Delta t){{\rm{e}}^{\theta\Delta t}} + {F_0}(t,x +{\mu_0}\Delta t)(1 - {{\rm{e}}^{\xi\Delta t}}),
\end{aligned}
\end{equation*}Substituting the exponential form of the Taylor formulas with Peano terms, we obtain:
\begin{equation*}
\begin{aligned}
&{F_0}(t + \Delta t,x) = {F_0}(t,x +{\mu_0}\Delta t)(1 - \xi\Delta t) + {F_1}(t,x - \left( {{\lambda _1} + {\lambda _2} - {\mu_1}} \right)\Delta t)\theta\Delta t + o(\Delta t),\\
&{F_1}(t + \Delta t,x) = {F_1}(t,x - \left( {{\lambda _1} + {\lambda _2} - {\mu_1}} \right)\Delta t)(1 - \theta\Delta t) + {F_0}(t,x + {\mu_0}\Delta t)\xi\Delta t + o(\Delta t),
\end{aligned}
\end{equation*} Dividing both sides by 
$\Delta t$ and making 
$\Delta t \to 0$, then we can get:
\begin{equation*}\begin{array}{l}
\left\{
\begin{aligned}
&\frac{{\partial {F_0}(t,x)}}{{\partial t}} -{\mu_0} \frac{{\partial {F_0}(t,x)}}{{\partial x}} = -\xi{F_0}(t,x) + \theta{F_1}(t,x),\\
&\frac{{\partial {F_1}(t,x)}}{{\partial t}} + \left( {{\lambda _1} + {\lambda _2} - {\mu_1}} \right)\frac{{\partial {F_1}(t,x)}}{{\partial x}} = - \theta{F_1}(t,x) + \xi{F_0}(t,x),
\end{aligned}
\right.
\end{array}\end{equation*} When the system is stable, there exists 
$\mathop {\lim }\limits_{t \to \infty } \frac{{\partial {F_i}(t,x)}}{{\partial t}} = 0,i = 0,1.$
The ordinary differential equations can be constructed as:
\begin{equation}\begin{array}{l}
\left\{
\begin{aligned}
&\ -{\mu_0}\frac{{{\rm d}{F_0}(x)}}{{{\rm d}x}} = -\xi {F_0}(x) +\theta {F_1}(x),\\
&\left( {{\lambda _1} + {\lambda _2} - {\mu_1}} \right)\frac{{{\rm d}{F_1}(x)}}{{{\rm d}x}} = \xi {F_0}(x) - \theta {F_1}(x),
\end{aligned}
\right.
\end{array}\end{equation} with boundary conditions 
${F_0}({x_{e1}}(1)) = {\pi _0}$ and 
$\ {F_1}({x_{e2}}(0)) = 0.$
The formulas (6) and (7) can be obtained by solving the differential equation (11). The point masses and the probability density functions in two states, respectively, are:
\begin{equation*}{P_0}({x_{e2}}(0))=F_{0}({x_{e2}}(0)),\ {P_1}({x_{e1}}(1))=\pi_{1}-F_{1}({x_{e1}}(1)),\end{equation*}
\begin{equation*}{f_0}(x) = \frac{{{\pi _0}\xi ({\lambda _1} + {\lambda _2} - {\mu_1})\left( {\frac{\xi }{{{\mu_0}}} - \frac{\theta }{{{\lambda _1} + {\lambda _2} - {\mu_1}}}} \right){{\rm e}^{- \left( {\frac{\theta }{{{\lambda _1} + {\lambda _2} - {\mu_1}}} - \frac{\xi }{{{\mu_0}}}} \right)\left( {x - {x_{e2}}(0)} \right)}}}}{{\xi ({\lambda _1} + {\lambda _2} - {\mu_1}){{\rm e}^{- \left( {\frac{\theta }{{{\lambda _1} + {\lambda _2} - {\mu_1}}} - \frac{\xi }{{{\mu_0}}}} \right)\left( {{x_{e1}}(1) - {x_{e2}}(0)} \right)}} - \theta {\mu_0}}},\end{equation*}
\begin{equation*}{f_1}(x) = \frac{{{\pi _1}\theta {\mu_0}\left( {\frac{\xi }{{{\mu_0}}} - \frac{\theta }{{{\lambda _1} + {\lambda _2} - {\mu_1}}}} \right){{\rm e}^{- \left( {\frac{\theta }{{{\lambda _1} + {\lambda _2} - {\mu_1}}} - \frac{\xi }{{{\mu_0}}}} \right)\left( {x - {x_{e2}}(0)} \right)}}}}{{\xi \left( {{\lambda _1} + {\lambda _2} - {\mu_1}} \right){{\rm e}^{- \left( {\frac{\theta }{{{\lambda _1} + {\lambda _2} - {\mu_1}}} - \frac{\xi }{{{\mu_0}}}} \right)\left( {{x_{e1}}(1) - {x_{e2}}(0)} \right)}} - \theta {\mu_0}}}.\end{equation*} Case III. 
${\lambda _{\rm{1}}}{\rm{+ }}{\lambda _{\rm{2}}}{\rm{= }}{\mu_{\rm{1}}}$.
When the buffer is in state 0, the net input rate of the fluid is 
${\lambda _{\rm{1}}}{\rm{+ }}{\lambda _{\rm{2}}}-\mu_0$ or 
${\lambda _{\rm{2}}}-\mu_0$ and the fluid level increases to 
${{x}_{e2}}\left( 0 \right)$, then the fluid level remains unaltered. When the buffer is in state 1, the fluid level remains constant. The fluid level eventually stabilizes at 
${{x}_{e2}}\left( 0 \right)$ with the alternation of the server status.
Case IV. 
${\mu_{\rm{0}}} \lt {\lambda _1} + {\lambda _2} \lt {\mu_1}$.
When the buffer is in state 0, the fluid in the system increases at the rate of 
${{\lambda }_{1}}+{{\lambda }_{2}}-{{\mu }_{0}}$ until the fluid level reaches the threshold 
${{x}_{e1}}(0)$. When the fluid level exceeds the optimal threshold 
${{x}_{e1}}(0)$, only type 2 customers flow in and the fluid level increases to 
${{x}_{e2}}(0)$ at the rate of 
${{\lambda }_{2}}-{{\mu }_{0}}$. If the buffer is still in state 0, the fluid will enter the system with probability 
$\frac{{{\mu }_{0}}}{{{\lambda }_{2}}}$ to ensure that the fluid level is stabilized at 
${{x}_{e2}}(0)$. When the buffer is in state 1, the fluid decreases at the rate of 
$\lambda_1+\lambda_2-{{\mu }_{1}}$. The system repeats in this way, and the fluid level fluctuates within 
$\left[ 0,{{x}_{e2}}(0) \right]$.
Considering the fluid level in a tiny time interval, we have:
\begin{equation*}
\begin{aligned}
&{F_0}(t + \Delta t,x) = {F_0}(t,x -\left( {p{\lambda _1} + {\lambda _2} - {\mu_0}} \right)\Delta t){{\rm{e}}^{\xi \Delta t}} + {F_1}(t,x - \left( {{\lambda _1} + {\lambda _2} - {\mu_1}} \right)\Delta t)(1 - {{\rm{e}}^{\theta\Delta t}}),\\
&{F_1}(t + \Delta t,x) = {F_1}(t,x - \left( {{\lambda _1} + {\lambda _2} - {\mu_1}} \right)\Delta t){{\rm{e}}^{\theta\Delta t}} + {F_0}(t,x -\left( {p{\lambda _1} + {\lambda _2} - {\mu_0}} \right)\Delta t)(1 - {{\rm{e}}^{\xi\Delta t}}),
\end{aligned}
\end{equation*}Substituting the exponential form of the Taylor formulas with Peano terms, we obtain:
\begin{equation*}
\begin{aligned}
{F_0}(t + \Delta t,x) =& {F_0}(t,x -\left( {p{\lambda _1} + {\lambda _2} - {\mu_0}} \right)\Delta t)(1 - \xi\Delta t) \\
&+ {F_1}(t,x - \left( {{\lambda _1} + {\lambda _2} - {\mu_1}} \right)\Delta t)\theta\Delta t + o(\Delta t),\\
{F_1}(t + \Delta t,x) =& {F_1}(t,x - \left( {{\lambda _1} + {\lambda _2} - {\mu_1}} \right)\Delta t)(1 - \theta\Delta t) \\
&+ {F_0}(t,x -\left( {p{\lambda _1} + {\lambda _2} - {\mu_0}} \right)\Delta t)\xi\Delta t + o(\Delta t),
\end{aligned}
\end{equation*} Dividing both sides by 
$\Delta t$ and making 
$\Delta t \to 0$, then we can get:
\begin{equation*}\begin{array}{l}
\left\{
\begin{aligned}
&\frac{{\partial {F_0}(t,x)}}{{\partial t}} +\left( {p{\lambda _1} + {\lambda _2} - {\mu_0}} \right) \frac{{\partial {F_0}(t,x)}}{{\partial x}} = -\xi{F_0}(t,x) + \theta{F_1}(t,x),\\
&\frac{{\partial {F_1}(t,x)}}{{\partial t}} + \left( {{\lambda _1} + {\lambda _2} - {\mu_1}} \right)\frac{{\partial {F_1}(t,x)}}{{\partial x}} = - \theta{F_1}(t,x) + \xi{F_0}(t,x),
\end{aligned}
\right.
\end{array}\end{equation*} When the system is stable, there exists 
$\mathop {\lim }\limits_{t \to \infty } \frac{{\partial {F_i}(t,x)}}{{\partial t}} = 0,i = 0,1.$
The ordinary differential equations can be constructed as:
\begin{equation}\begin{array}{l}
\left\{
\begin{aligned}
&\left( {p{\lambda _1} + {\lambda _2} - {\mu_0}} \right)\frac{{{\rm d}{F_0}(x)}}{{{\rm d}x}} = \theta {F_1}(x)- \xi {F_0}(x) ,\\
&\left( {\lambda _1} + {\lambda _2} - {\mu_1}\right)\frac{{{\rm d}{F_1}(x)}}{{{\rm d}x}} =\xi {F_0}(x)- \theta {F_1}(x),
\end{aligned}
\right.
\end{array}\end{equation} with boundary conditions 
${F_0}({\rm{0}}) = {\rm{0}}$ and 
${F_1}({x_{e2}}(0)) = {\pi _{\rm{1}}}$, where 
$p=\frac{{{x}_{e1}}(0)}{{{x}_{e2}}(0)}, 0 \le p \le1$ is the probability that the fluid level in the system is less than the threshold 
${{x}_{e1}}(0)$ in state 0.
The formulas (8) and (9) can be obtained by calculating equation (12). The point masses and the probability density functions in two states, respectively, are:
\begin{equation*}{P_0}({x_{e2}}(0))=\pi_{0}-F_{0}({x_{e2}}(0)),\ P_1(0)=F_1(0),\end{equation*}
\begin{equation*}{f_0}(x) = \frac{{{\pi _0}\xi \left( {{\mu_1}-{\lambda _1} -{\lambda _2} } \right)\left( {\frac{\theta }{{{\mu_1}-{\lambda _1} - {\lambda _2} }}-\frac{\xi }{{p{\lambda _1} + {\lambda _2} - {\mu_0}}} } \right){{\rm e}^{- \left( {\frac{\xi }{{p{\lambda _1} + {\lambda _2} - {\mu_0}}} + \frac{\theta }{{{\lambda _1} + {\lambda _2} - {\mu_1}}}} \right)x}}}}{{\theta \left( {p{\lambda _1} + {\lambda _2} - {\mu_0}} \right){{\rm e}^{- \left( {\frac{\xi }{{p{\lambda _1} + {\lambda _2} - {\mu_0}}} + \frac{\theta }{{{\lambda _1} + {\lambda _2} - {\mu_1}}}} \right){x_{e2}}(0)}} + \xi \left( {{\lambda _1} + {\lambda _2} - {\mu_1}} \right)}},\end{equation*}
\begin{equation*}{f_1}(x) = \frac{{ {\pi _1}(p{\lambda _1} + {\lambda _2} - {\mu_0})\left( \frac{\theta }{{{\mu_1}-{\lambda _1} - {\lambda _2}}} -\frac{\xi }{p{\lambda _1} + {\lambda _2} - {\mu_0}} \right)\theta {{\rm e}^{- \left( {\frac{\xi }{{p{\lambda _1} + {\lambda _2} - {\mu_0}}} + \frac{\theta }{{{\lambda _1} + {\lambda _2} - {\mu_1}}}} \right)x}}}}{{\theta (p{\lambda _1} + {\lambda _2} - {\mu_0}){{\rm e}^{- \left( {\frac{\xi }{{p{\lambda _1} + {\lambda _2} - {\mu_0}}} + \frac{\theta }{{{\lambda _1} + {\lambda _2} - {\mu_1}}}} \right){x_{e2}}(0)}} + \xi ({\lambda _1} + {\lambda _2} - {\mu_1})}}.\end{equation*}4.2. Analysis of the expected social benefit
Based on Theorem 4.1, when all customers follow the optimal balking strategies 
$ \left({{x}_{e1}}(0),{{x}_{e2}}(0),{{x}_{e1}}(1),{{x}_{e2}}(1)\right)$, the utility functions of the expected social benefit per unit time are as follows.
Case I. 
${\lambda}_{2} \gt \mu_{1}$.
In this case, the fluid level fluctuates within 
$\left[ {{x}_{e2}}(0),{{x}_{e2}}(1) \right]$. The type 1 customers refuse to flow in when the fluid level reaches the socially optimal thresholds 
${{x}_{e1}}(i), i=0,1.$ The entrance probabilities of type 2 customers are 
$\frac{\mu_0}{\lambda_2}$ and 
$\frac{\mu_1}{\lambda_2}$ when the fluid levels satisfy 
$x={{x}_{e2}}(0)$ and 
$x={{x}_{e2}}(1)$, respectively. Based on Theorem 4.1, the stationary probability distribution 
$F_1(x)$ is continuous when the fluid level reaches the thresholds 
${{x}_{e2}}(0)$ and 
${{x}_{e1}}(1)$, then the transient probabilities are 
$P_1({{x}_{e2}}(0))=0$ and 
$P_1({{x}_{e1}}(1))=0$.
The expected social benefit per unit time can be expressed as:
\begin{equation*}B\left( {{x_{e1}}(0),{x_{e2}}(0),{x_{e1}}(1),{x_{e2}}(1)} \right) = {\lambda _{e1}}{R_1} + {\lambda _{e2}}{R_2} - \left[ {{C_1}E\left( {{{\rm e}^{\alpha {X_1}}}} \right) + {C_2}E\left( {{{\rm e}^{\alpha {X_2}}}} \right)} \right],\end{equation*}where the effective arrival rate and the exponential expected fluid level, respectively, are:
\begin{equation*}
{\lambda _{e1}} = {\lambda _1}\int_{{x_{e2}}(0)}^{{x_{e1}}(1)} {{f_1}(x)}{\rm d}x = {\lambda _1}\frac{{{\pi _1}{\mu_0}\theta {{\rm e}^{- {A_1}\left( {{x_{e1}}(1) - 2{x_{e2}}(0)} \right)}}}}{{\xi \left( {q{\lambda _1} + {\lambda _2} - {\mu_1}} \right){{\rm e}^{- {A_1}\left( {{x_{e2}}(1) - {x_{e2}}(0)} \right)}} - \theta {\mu_0}}},
\end{equation*}
\begin{equation*}
\begin{aligned}
{\lambda _{e2}} &= {\lambda _2}\left[ {{P_0}({x_{e2}}(0))\frac{{{\mu_0}}}{{{\lambda _2}}} + {P_1}({x_{e2}}(1))\frac{{{\mu_1}}}{{{\lambda _2}}} + \int_{{x_{e2}}(0)}^{{x_{e2}}(1)} {{f_1}(x){\rm d}x} } \right]\\
&= {\lambda _2}\frac{{{\pi _1}{\mu_0}\theta {{\rm e}^{- {A_1}\left( {{x_{e2}}(1) - 2{x_{e2}}(0)} \right)}} + {\pi _0}{B_1}\frac{{{\mu_0}}}{{{\lambda _2}}} + {\pi _1}{B_1}{{\rm e}^{- {A_1}\left( {{x_{e2}}(1) - {x_{e2}}(0)} \right)}}\frac{{{\mu_1}}}{{{\lambda _2}}}}}{{\xi \left( {q{\lambda _1} + {\lambda _2} - {\mu_1}} \right){{\rm e}^{- {A_1}\left( {{x_{e2}}(1) - {x_{e2}}(0)} \right)}} - \theta {\mu_0}}},
\end{aligned}
\end{equation*}
\begin{equation*}
\begin{aligned}
E({{\rm e}^{\alpha {X_1}}}) &= \int_{{x_{e2}}(0)}^{{x_{e1}}(1)} {{f_1}(x){{\rm e}^{\alpha \frac{{{\lambda _1}}}{{{\lambda _1} + {\lambda _2}}}x}}{\rm d}x + } \int_{{x_{e2}}(0)}^{{x_{e2}}(1)} {{f_0}(x){{\rm e}^{\alpha \frac{{{\lambda _1}}}{{{\lambda _1} + {\lambda _2}}}x}}{\rm d}x} \\
&= \frac{{ {\pi _1}\theta (q{\lambda _1} + {\lambda _2} - {\mu_1} + {\mu_0}){A_1}\left[ {{{\rm e}^{- {A_1}\left( {{x_{e1}}(1) - 2{x_{e2}}(0)} \right)}} + {{\rm e}^{- {A_1}\left( {{x_{e2}}(1) - 2{x_{e2}}(0)} \right)}}} \right]}}{{\left[ {\xi \left( {q{\lambda _1} + {\lambda _2} - {\mu_1}} \right){{\rm e}^{- {A_1}\left( {{x_{e2}}(1) - {x_{e2}}(0)} \right)}} - \theta {\mu_0}} \right]\left( {{A_1} - {E_1}} \right)}},
\end{aligned}
\end{equation*}
\begin{equation*}
\begin{aligned}
E\left( {{{\rm e}^{\alpha {X_2}}}} \right)=& \int_{{x_{e2}}(0)}^{{x_{e1}}(1)} {{f_1}(x){{\rm e}^{{D_1}x}}{\rm d}x + \int_{{x_{e1}}(1)}^{{x_{e2}}(1)} {{f_1}(x){{\rm e}^{\alpha x}}{\rm d}x + {P_1}({x_{e2}}(1)){{\rm e}^{\alpha \left[ {{x_{e2}}(1) - \frac{{{\lambda _1}}}{{{\lambda _1} + {\lambda _2}}}{x_{e1}}(1)} \right]}}} } \\
&+ \int_{{x_{e2}}(0)}^{{x_{e2}}(1)} {{f_0}(x){{{\rm e}}^{\alpha \frac{{{\lambda _2}}}{{{\lambda _1} + {\lambda _2}}}x}}{\rm d}x + {P_0}({x_{e2}}(0)){{{\rm e}}^{\alpha {x_{e2}}(0)}}} \\
= &\frac{{{\pi _1}\theta {\mu_0}{A_1}{{\rm e}^{{A_1}{x_{e2}}(0)}}\left( {\frac{{{{\rm e}^{\left( {{D_1} - {A_1}} \right)\left( {{x_{e1}}(1) - {x_{e2}}(0)} \right)}}}}{{{A_1} - {D_1}}} + \frac{{{{\rm e}^{\left( {\alpha - {A_1}} \right)\left( {{x_{e2}}(1) - {x_{e1}}(1)} \right)}}}}{{{A_1}-\alpha}}} \right)}}{{\xi \left( {q{\lambda _1} + {\lambda _2} - {\mu_1}} \right){{\rm e}^{- {A_1}\left( {{x_{e2}}(1) - {x_{e2}}(0)} \right)}} - \theta {\mu_0}}}\\
&+ \frac{{{\pi _0}{B_1}{{\rm e}^{\alpha {x_{e2}}(0)}}+{\pi _0}\xi (q{\lambda _1} + {\lambda _2} - {\mu_1}){A_1}\frac{{{{\rm e}^{\left( {{D_1} - {A_1}} \right)\left( {{x_{e2}}(1) - {x_{e2}}(0)} \right) + {A_1}{x_{e2}}(0)}}}}{{{A_1} - {D_1}}} }}{{\xi \left( {q{\lambda _1} + {\lambda _2} - {\mu_1}} \right){{\rm e}^{- {A_1}\left( {{x_{e2}}(1) - {x_{e2}}(0)} \right)}} - \theta {\mu_0}}}\\
&+ \frac{{\pi _1}{B_1}{{\rm e}^{\left[ {\left( {\alpha - {A_1}} \right){x_{e2}}(1) - {E_1}{x_{e1}}(1)} \right]}}{{\rm e}^{{A_1}{x_{e2}}(0)}}}{{\xi \left( {q{\lambda _1} + {\lambda _2} - {\mu_1}} \right){{\rm e}^{- {A_1}\left( {{x_{e2}}(1) - {x_{e2}}(0)} \right)}} - \theta {\mu_0}}},
\end{aligned}
\end{equation*} where 
${A_1} = \frac{\theta }{{q{\lambda _1} + {\lambda _2} - {\mu_1}}} - \frac{\xi }{{{\mu_0}}}, {B_1} = \xi \left( {q{\lambda _1} + {\lambda _2} - {\mu_1}} \right) - \theta {\mu_0}, {D_1} = \alpha \frac{{{\lambda _2}}}{{{\lambda _1} + {\lambda _2}}}, {E_1} = \alpha \frac{{{\lambda _1}}}{{{\lambda _1} + {\lambda _2}}}.$
Case II. 
${\lambda}_{1} + {\lambda}_{2} \gt \mu_{1}$ and 
${\lambda}_{2} \le \mu_{1}$.
In this case, the fluid level fluctuates within 
$\left[ {{x}_{e2}}(0),{{x}_{e1}}(1) \right]$. The type 1 customers refuse to flow in when the fluid level reaches the socially optimal thresholds 
${{x}_{e1}}(0)$ in state 0. The entrance probability of type 1 customers is 
$\frac{\mu_1}{\lambda_1+\lambda_2}$ when the fluid levels satisfy 
$x={{x}_{e1}}(1)$. The entrance probabilities of type 2 customers are 
$\frac{\mu_0}{\lambda_2}$ and 
$\frac{\mu_1}{\lambda_1+\lambda_2}$ when the fluid levels satisfy 
$x={{x}_{e2}}(0)$ and 
$x={{x}_{e1}}(1)$, respectively. Based on Theorem 4.1, the stationary probability distribution 
$F_1(x)$ is continuous when the fluid level reaches the threshold 
${{x}_{e2}}(0)$, then the transient probability is 
$P_1({{x}_{e2}}(0))=0$.
The expected social benefit per unit time can be expressed as:
\begin{equation*}B\left( {{x_{e1}}(0),{x_{e2}}(0),{x_{e1}}(1),{x_{e2}}(1)} \right) = {\lambda _{e1}}{R_1} + {\lambda _{e2}}{R_2} - \left[ {{C_1}E\left( {{{\rm e}^{\alpha {X_1}}}} \right) + {C_2}E\left( {{{\rm e}^{\alpha {X_2}}}} \right)} \right],\end{equation*}where the effective arrival rates for two types of customers and the exponential expected fluid levels in the system, respectively, are:
\begin{equation*}
\begin{aligned}
{\lambda _{e1}}& = {\lambda _1}\left[ {\int_{{x_{e2}}(0)}^{{x_{e1}}(1)} {{f_1}(x){\rm d}x + \frac{{{\mu_1}}}{{{\lambda _1} + {\lambda _2}}}{P_1}({x_{e1}}(1))} } \right]\\
& = {\lambda _1}\frac{{{\pi _1}\theta {\mu_0}{{\rm e}^{- {A_2}\left( {{x_{e1}}(1) - 2{x_{e2}}(0)} \right)}} + {\pi _1}{F_2}\frac{{{\mu_1}}}{{{\lambda _1} + {\lambda _2}}}{{\rm e}^{- {A_2}\left( {{x_{e1}}(1) - {x_{e2}}(0)} \right)}}}}{{\xi \left( {{\lambda _1} + {\lambda _2} - {\mu_1}} \right){{\rm e}^{- {A_2}\left( {{x_{e1}}(1) - {x_{e2}}(0)} \right)}} - \theta {\mu_0}}},
\end{aligned}
\end{equation*}
\begin{equation*}
\begin{aligned}
{\lambda _{e2}}& = {\lambda _2}\left[ {{P_0}({x_{e2}}(0))\frac{{{\mu_0}}}{{{\lambda _2}}} + \int_{{x_{e2}}(0)}^{{x_{e1}}(1)} {{f_1}(x){\rm d}x + {P_1}({x_{e1}}(1))\frac{{{\mu_1}}}{{{\lambda _1} + {\lambda _2}}}} } \right]\\
&= {\lambda _2}\frac{{{\pi _0}{D_2}\frac{{{\mu_0}}}{{{\lambda _2}}} + {\pi _1}\theta {\mu_0}{{\rm e}^{- {A_2}\left( {{x_{e1}}(1) - 2{x_{e2}}(0)} \right)}} + {\pi _1}{F_2}\frac{{{\mu_1}}}{{{\lambda _1} + {\lambda _2}}}{{\rm e}^{- {A_2}\left( {{x_{e1}}(1) - {x_{e2}}(0)} \right)}}}}{{\xi \left( {{\lambda _1} + {\lambda _2} - {\mu_1}} \right){{\rm e}^{- {A_2}\left( {{x_{e1}}(1) - {x_{e2}}(0)} \right)}} - \theta {\mu_0}}},
\end{aligned}
\end{equation*}
\begin{equation*}
\begin{aligned}
E({{\rm e}^{\alpha {X_1}}})=& \int_{{x_{e2}}(0)}^{{x_{e1}}(1)} {\left( {{f_0}(x) + {f_1}(x)} \right)} {{\rm e}^{\alpha \frac{{{\lambda _1}}}{{{\lambda _1} + {\lambda _2}}}x}}{\rm d}x + {P_1}({x_{e1}}(1)){{{\rm e}}^{\alpha \frac{{{\lambda _1}}}{{{\lambda _1} + {\lambda _2}}}{x_{e1}}(1)}}\\
=& \frac{{{\pi _0}\xi ({\lambda _1} + {\lambda _2} + {\mu_0} - {\mu_1}){A_2}\frac{{{{\rm e}^{{A_2}{x_{e2}}(0) + \left( {{E_1} - {A_2}} \right)\left( {{x_{e1}}(1) - {x_{e2}}(0)} \right)}}}}{{{A_2} - {E_1}}}}}{{\xi \left( {{\lambda _1} + {\lambda _2} - {\mu_1}} \right){{\rm e}^{- {A_2}\left( {{x_{e1}}(1) - {x_{e2}}(0)} \right)}} - \theta {\mu_0}}}\\
&+\frac{{\pi _1}{F_2}{{\rm e}^{{E_1}{x_{e1}}(1) - {A_2}\left( {{x_{e1}}(1) - {x_{e2}}(0)} \right)}}}{{\xi \left( {{\lambda _1} + {\lambda _2} - {\mu_1}} \right){{\rm e}^{- {A_2}\left( {{x_{e1}}(1) - {x_{e2}}(0)} \right)}} - \theta {\mu_0}}},
\end{aligned}
\end{equation*}
\begin{equation*}
\begin{aligned}
E\left( {{{\rm e}^{\alpha {X_2}}}} \right)=& \int_{{x_{e2}}(0)}^{{x_{e1}}(1)} {\left( {{f_0}(x) + {f_1}(x)} \right){{\rm e}^{\alpha \frac{{{\lambda _2}}}{{{\lambda _1} + {\lambda _2}}}x}}{\rm d}x + {P_0}} ({x_{e2}}(0)){{\rm e}^{\alpha {x_{e2}}(0)}}\\
&+ {P_1}({x_{e1}}(1)){{\rm e}^{\alpha \frac{{{\lambda _2}}}{{{\lambda _1} + {\lambda _2}}}{x_{e1}}(1)}}\\
=& \frac{{ {\pi _0}\xi ({\lambda _1} + {\lambda _2} + {\mu_0} - {\mu_1}){A_2}\frac{{{{\rm e}^{{A_2}{x_{e2}}(0)}}}}{{{A_2} - {D_1}}}{{\rm e}^{\left( {{D_1} - {A_2}} \right)\left( {{x_{e1}}(1) - {x_{e2}}(0)} \right)}}}}{{\xi \left( {{\lambda _1} + {\lambda _2} - {\mu_1}} \right){{\rm e}^{- {A_2}\left( {{x_{e1}}(1) - {x_{e2}}(0)} \right)}} - \theta {\mu_0}}}\\
&+ \frac{{{\pi _1}{F_2}{{\rm e}^{- {A_2}\left( {{x_{e1}}(1) - {x_{e2}}(0)} \right)}}{{\rm e}^{{D_1}{x_{e1}}(1)}} + {\pi _1}{D_2}{{\rm e}^{\alpha {x_{e2}}(0)}}}}{{\xi \left( {{\lambda _1} + {\lambda _2} - {\mu_1}} \right){{\rm e}^{- {A_2}\left( {{x_{e1}}(1) - {x_{e2}}(0)} \right)}} - \theta {\mu_0}}},
\end{aligned}
\end{equation*} where 
${A_2} = \frac{\theta }{{{\lambda _1} + {\lambda _2} - {\mu_1}}} - \frac{\xi }{{{\mu_0}}},\ {D_2} = \xi ({\lambda _1} + {\lambda _2} - {\mu_1})- {\mu_0}\theta,\ {F_2} = \xi \left( {{\lambda _1} + {\lambda _2} - {\mu_1}} \right) + \theta {\mu_0}.$
Case III. 
${\lambda _{\rm{1}}}{\rm{+ }}{\lambda _{\rm{2}}}{\rm{= }}{\mu_{\rm{1}}}$.
In this case, the fluid level stabilizes at 
${{x}_{e2}}(0)$. The type 1 customers refuse to flow in and the entrance probability of type 2 customers is 
$\frac{\mu_0}{\lambda_2}$ when the system is in state 0. Besides, due to 
${{x}_{e2}}(0) \lt {{x}_{e1}}(1) \lt {{x}_{e2}}(1)$, two types of parallel customers flow into the buffer when the system is in state 1.
The expected social benefit per unit time is:
\begin{equation*}
\begin{aligned}
&B\left( {{x_{e1}}(0),{x_{e2}}(0),{x_{e1}}(1),{x_{e2}}(1)} \right) \\
=&{\lambda _1}{\pi _1}{R_1} + {\lambda _2}\left( {{\pi _1} + {\pi _0}\frac{{{\mu_0}}}{{{\lambda _2}}}} \right){R_2} - \left[ {C_1}{E\left( {{{\rm e}^{\alpha \frac{{{\lambda _1}}}{{{\lambda _1} + {\lambda _2}}}{x_{e2}}(0)}}} \right) + {C_2}E\left( {{{\rm e}^{\alpha \frac{{{\lambda _2}}}{{{\lambda _1} + {\lambda _2}}}{x_{e2}}(0)}}} \right)} \right].\end{aligned}
\end{equation*} Case IV. 
${\mu_{\rm{0}}} \lt {\lambda _1} + {\lambda _2} \lt {\mu_1}$.
In this case, the fluid level fluctuates within 
$\left[ 0,{{x}_{e2}}(0) \right]$. The fluids choose to flow in when the fluid level satisfies 
$x\le {{x}_{ej}}(0),j=1,2$, and the entrance probability of type 2 customers is 
$\frac{\mu_0}{\lambda_2}$ when the fluid level satisfies 
$x={{x}_{e2}}(0)$ in state 0. Besides, due to 
${{x}_{e2}}(0) \lt {{x}_{e1}}(1) \lt {{x}_{e2}}(1)$, two types of parallel customers flow into the buffer when the system is in state 1. Based on Theorem 4.1, the stationary probability distribution 
$F_0(x)$ is continuous when the fluid levels satisfy x = 0 and 
$x=x_{e1}(0)$, then the transient probabilities are 
$P_0(0)=0$ and 
$P_0(x_{e1}(0))=0$. The stationary probability distribution 
$F_1(x)$ is continuous when the fluid level satisfies 
$x=x_{e2}(0)$, then the transient probability is 
$P_1(x_{e2}(0))=0$.
The expected social benefit per unit time is:
\begin{equation*}B\left( {{x_{e1}}(0),{x_{e2}}(0),{x_{e1}}(1),{x_{e2}}(1)} \right) = {\lambda _{e1}}{R_1} + {\lambda _{e2}}{R_2} - \left[ {{C_1}E\left( {{{\rm e}^{\alpha {X_1}}}} \right) + {C_2}E\left( {{{\rm e}^{\alpha {X_2}}}} \right)} \right],\end{equation*}where the effective arrival rates for two types of customers and the exponential expected fluid levels in the system respectively are:
\begin{equation*}
{\lambda _{e1}} = {\lambda _1}\left(\int_0^{{x_{e1}}(0)} {{f_0}(x){\rm d}x + {\pi _1}}\right) = {\lambda _1}\left( {{\pi _1} - \frac{{{\pi _0}\xi {E_3}{{\rm e}^{- {A_3}{x_{e1}}(0)}}}}{{\theta {D_3}{{\rm e}^{- {A_3}{x_{e2}}(0)}} + \xi {E_3}}}} \right),
\end{equation*}
\begin{equation*}\begin{aligned}
{\lambda _{e2}} =& {\lambda _2}\left( {\int_0^{{x_{e2}}(0)} {{f_0}(x){\rm d}x + {\pi _1}} } +{P_0}({x_{e2}}(0))\frac{{{\mu_0}}}{{{\lambda _2}}}\right)\\
{\rm{= }}&{\lambda _2}\frac{{\left[ {{\pi _1}\theta {D_3} + {\pi _0}\left( {(\theta {D_3} + \xi {E_3})\frac{{{\mu_0}}}{{{\lambda _2}}} - \xi {E_3}} \right)} \right]{{\rm{e}}^{- {A_3}{x_{e2}}(0)}} + \xi {E_3}}}{{\theta {D_3}{e^{- {A_3}{x_{e2}}(0)}} + \xi {E_3}}},
\end{aligned}
\end{equation*}
\begin{equation*}
\begin{aligned}
E\left( {{{\rm e}^{\alpha {X_1}}}} \right) &= \int_0^{{x_{e1}}(0)} {{f_0}(x){{\rm e}^{\alpha \frac{{{\lambda _1}}}{{{\lambda _1} + {\lambda _2}}}x}}{\rm d}x} + \int_0^{{x_{e2}}(0)} {{f_1}(x){{\rm e}^{\alpha \frac{{{\lambda _1}}}{{{\lambda _1} + {\lambda _2}}}x}}{\rm d}x} \\
&= \frac{{{\pi _1}\theta {A_3}\left[ {{E_3}{{\rm e}^{\left( {{E_1} - {A_3}} \right){x_{e1}}(0)}} - {D_3}{{\rm e}^{\left( {{E_1} - {A_3}} \right){x_{e2}}(0)}}} \right]}}{{\left[ {\theta {D_3}{{\rm e}^{- {A_3}{x_{e2}}(0)}} + \xi {E_3}} \right]\left( {{E_1} - {A_3}} \right)}},
\end{aligned}
\end{equation*}
\begin{equation*}
\begin{aligned}
E\left( {{{\rm e}^{\alpha {X_2}}}} \right)=& \int_0^{{x_{e1}}(0)} {{f_0}(x){{\rm e}^{{D_1}x}}{\rm d}x + \int_{{x_{e1}}(0)}^{{x_{e2}}(0)} {{f_0}(x){{\rm e}^{\alpha x}}} }{\rm d}x + \int_0^{{x_{e2}}(0)} {{f_1}(x){{\rm e}^{{D_1}x}}{\rm d}x} \\
&+ {P_0}({x_{e2}}(0)){{\rm e}^{\alpha \left[ {{x_{e2}}(0) - \frac{{{\lambda _1}}}{{{\lambda _1} + {\lambda _2}}}{x_{e1}}(0)} \right]}}\\
=& \frac{{{\pi _0}\xi {A_3}\left[ {\frac{{{E_3}{{\rm e}^{\left( {{D_1} - {A_3}} \right){x_{e1}}(0)}} + {D_3}{{\rm e}^{\left( {{D_1} - {A_3}} \right){x_{e2}}(0)}}}}{{{D_1} - {A_3}}} + \frac{{{E_3}{{\rm e}^{\left( {\alpha - {A_3}} \right)\left( {{x_{e2}}(0) - {x_{e1}}(0)} \right)}}}}{{\alpha - {A_3}}}} \right]}}{{\theta {D_3}{{\rm e}^{- {A_3}{x_{e2}}(0)}} + \xi {E_3}}}\\
&+ \frac{{{\pi _0}\left[ {\xi {E_3} + \theta {D_3}} \right]{{\rm e}^{\alpha \left[ {{x_{e2}}(0) - \frac{{{\lambda _1}}}{{{\lambda _1} + {\lambda _2}}}{x_{e1}}(0)} \right] - {A_3}{x_{e2}}(0)}}}}{{\theta {D_3}{{\rm e}^{- {A_3}{x_{e2}}(0)}} + \xi {E_3}}},
\end{aligned}
\end{equation*} where 
${A_3} = \frac{\xi }{{p{\lambda _1} + {\lambda _2} - {\mu_0}}} + \frac{\theta }{{{\lambda _1} + {\lambda _2} - {\mu_1}}},\ {D_3} = p{\lambda _1} + {\lambda _2} - {\mu_0},\ {E_3} = {\lambda _1} + {\lambda _2} - {\mu_1}.$
4.3. Numerical comparisons of the expected social benefit
The optimal balking strategy 
$\left({{x}_{e1}}(0),{{x}_{e2}}(0),{{x}_{e1}}(1),{{x}_{e2}}(1) \right)$ is assumed to obtain the maximal expected social benefit per unit time. However, the optimal balking strategy cannot be accomplished by conventional mathematical methods due to the complexity of the utility function 
$B\left( {{x_{e1}}(0),{x_{e2}}(0),{x_{e1}}(1),{x_{e2}}(1)} \right)$. This subsection continues to investigate the effects of the optimal balking strategy on the expected social benefit by some numerical examples. Then, we assume that 
${\lambda _1}{\rm{= 2}}{\rm{.8,}}\ \lambda_2=3,\ \theta {\rm{= 0}}{\rm{.4,}}\ \xi {\rm{= 0}}{\rm{.3,}}\ {\mu_0}{\rm{= 1,}}\ {{{R}}_1}{\rm{=10,}}\ {{{R}}_2}{\rm{= 15,}}\ {{{C}}_1}{\rm{= 3,}}\ {{{C}}_2}{\rm{= 4,}}\ \alpha {\rm{= 0}}{\rm{.2}}$ in the following cases.
Figure 1 shows the variation of the expected social benefit per unit time with the thresholds 
${{x}_{e2}}(0)$ and 
${{x}_{e1}}(1)$ when 
${{\mu }_{1}}=3$. In this case, the type 2 customers’ arrival rate is equal to the normal service rate, that is 
${{\lambda }_{2}}={{\mu }_{1}} \lt {{\lambda }_{1}}+{{\lambda }_{2}}$. The expected social benefit decreases with the increase of the type 2 customers’ threshold 
${{x}_{e2}}(0)$ when the type 1 customers’ threshold 
${{x}_{e1}}(1)$ is fixed. On the contrary, the expected social benefit per unit time increases first and then decreases with the increase of the threshold 
${{x}_{e1}}(1)$ when the threshold 
${{x}_{e2}}(0)$ is fixed. This indicates that if each fluid flows into the buffer selfishly, the system will be over-congested and cannot achieve the global maximum. Therefore, the system designers could appropriately regulate the joining thresholds for better deployment management. The expected social benefit is concave, and the maximal benefit 
$B_{\rm max}=53.7$ can be obtained at the thresholds 
${{x}_{e2}}(0)=0$ and 
${{x}_{e1}}(1)=6.3$.
Figure 1. Expected social benefit versus thresholds 
${x_{ej}}(i)$ when 
${\mu_1} = 3$.
Figure 2. Expected social benefit versus thresholds 
${x_{ej}}(i)$ when 
${\mu_1} = 3.3$.
Figure 2 describes the sensitivity between the expected social benefit per unit time and the thresholds 
${{x}_{e2}}(0)$ and 
${{x}_{e1}}(1)$ when the outflow rate is 
${{\mu }_{1}}=3.3$. When parameters satisfy 
${{\lambda }_{2}} \lt {{\mu }_{1}} \lt {{\lambda }_{1}}+{{\lambda }_{2}}$, the variable tendency of the expected social benefit is positively correlated with the type 1 customers’ threshold 
${{x}_{e1}}(1)$ while negatively correlated with the type 2 customers’ threshold 
${{x}_{e2}}(0)$. As the optimal threshold 
${{x}_{e2}}(0)$ continues to increase, substantial customers are emerging in the system, and the buffer may be heavily loaded, which inevitably has a negative impact on society under the fully observable case. Evidently, when the threshold 
${{x}_{e1}}(1)$ is sufficiently large, the expected social benefit is more likely to reach its maximum. The expected social benefit achieves a maximum 
$B_{\rm max}=76$ at thresholds 
${{x}_{e2}}(0)=0$ and 
${{x}_{e1}}(1)=7.8$.
Figures 3 and 4 present the sensitivity between the socially optimal thresholds and the expected social benefit per unit time when the sum of the arrival rates is less than the outflow rate under the normal working state, that is 
${{\lambda }_{1}}+{{\lambda }_{2}} \lt {{\mu }_{1}}$. The expected social benefit accelerates its decline significantly when the type 1 customers’ threshold 
${{x}_{e1}}(0)$ increases. Moreover, when the thresholds for two types of customers are almost equal, the expected social benefit per unit time grows slowly and gradually stabilizes. Then, the expected social benefit reaches its maximum when the fluid flows into the buffer at this time.
Figure 3. Expected social benefit versus thresholds 
${x_{ej}}(0)$ when 
${\mu_1} = 6$.
Figure 4. Expected social benefit versus thresholds 
${x_{ej}}(0)$ when 
${\mu_1} = 6.5$.
By comparing Figures 3 and 4, we observe that when the service rate µ 1 increases, an arriving customer who is served in a shorter waiting time tends to have a stronger joining willingness. Concretely, the customers have more chances to be served at a higher service rate µ 1, which can reduce their overall sojourn time and increase the expected social benefit. Besides, the impact of the type 1 customers’ threshold 
${{x}_{e1}}(0)$ on the expected social benefit significantly reduces and the overall benefit tends to be stable when the outflow rate increases.
5. Analysis of the entrance fee strategy
Arrivals are self-interested and consider their own benefits first when encountering a joining/balking problem based on the non-cooperative game theory. The fluid inclines to flow out of the system as quickly as possible, then the buffer is saturated ultimately, which results in the over-utilization of resources and overcrowding of the system. The expected social benefit per unit time may not be globally optimal. Consequently, the decision-makers can impose an entrance fee 
${{p}_{j}(j=1,2)}$ on type j customers, which can make the fluid comply with the globally optimal threshold 
$x_{ej}^{*}(i)$, 
$x_{ej}^{*}(i) \gt x_{ej}(i)$. The entrance fees can relieve the queue congestion for systems and provide management insights for the optimal design of the fluid queues, which has theoretical and practical significance. An example is the entrance tickets for “fast-pass” in the amusement parks, and the separate fees and reservations for popular attractions can limit the queue length. Besides, some patients with life-threatening situations must choose to go to the emergency room and pay an expensive fee in order to receive services immediately and reduce the pressure on normal clinics.
The expected net benefit of the type j customers who observe the buffer at state i and decide to enter the system is yielded as:
\begin{equation}
B({x_j}(i)) = {R_j} - {p_j} - {C_j}\frac{{{{x}_{j}}(i) + 1}}{{{\mu_i}}},\end{equation} where 
${{x}_{j}}(i)$ is the fluid level of the type j customers in state i.
The expected net benefit should be non-negative to ensure that customers are willing to join. When the buffer is in a maintenance period, the entrance fee of the type j customers satisfies:
\begin{equation*}\left\{{\begin{array}{*{20}{c}}
{B(x_{ej}^*(0) - 1) = {R_j} - {C_j}\frac{{x_{ej}^*(0)}}{{{\mu_0}}} - {p_j} \ge 0},\\
{B(x_{ej}^*(0)) = {R_j} - {C_j}\frac{{x_{ej}^*(0) + 1}}{{{\mu_0}}} - {p_j} \lt 0},
\end{array}} \right.\end{equation*} that is, the entrance fee 
${p}_{j}$ satisfies:
\begin{equation*}{R_j} - {C_j}\frac{{x_{ej}^*(0) + 1}}{{{\mu_0}}} \lt p_j \le {R_j} - {C_j}\frac{{x_{ej}^*(0)}}{{{\mu_0}}}.\end{equation*} Then the maximal entrance fee 
$p_j^*$ in state 0 satisfies:
\begin{equation} p_j^* = {R_j} - {C_j}\frac{{x_{ej}^*(0)}}{{{\mu_0}}} \lt {R_j} - {C_j}\frac{{x_{ej}(0)}}{{{\mu_0}}}.\end{equation}On the other hand, when the buffer is in a normal working period, we have:
\begin{equation*}
\left\{{\begin{array}{*{20}{c}}
{B(x_{ej}^*(1) - 1) = {R_j} - {C_j}\frac{{x_{ej}^*(1)}}{{{\mu_1}}} - {p_j} \ge 0},\\
{B(x_{ej}^*(1)) = {R_j} - {C_j}\frac{{x_{ej}^*(1) + 1}}{{{\mu_1}}} - {p_j} \lt 0},
\end{array}} \right.
\end{equation*} that is, the entrance fee 
${p}_{j}$ satisfies:
\begin{equation*}{R_j} - {C_j}\frac{{x_{ej}^*(1) + 1}}{{{\mu_1}}} \lt p_j\le {R_j} - {C_j}\frac{{x_{ej}^*(1)}}{{{\mu_1}}}.\end{equation*} Thus, the maximal entrance fee 
$p_j^*$ in state 1 satisfies:
\begin{equation} p_j^*= {R_j} - {C_j}\frac{{x_{ej}^*(1)}}{{{\mu_1}}} \lt {R_j} - {C_j}\frac{{x_{ej}(1)}}{{{\mu_0}}}.\end{equation}The benefit for a service provider can be expressed as:
\begin{equation*}Z\left( {{x_1}(i),{x_2}(i)} \right) = {\lambda _1}[1 - {P_i}({x_1}(i))]{p_1} + {\lambda _2}[1 - {P_i}({x_2}(i))]{p_2},\end{equation*} where 
${P_i}({x_j}(i))$ is the probability that the type j fluid reaches the threshold 
${x_j}(i)$ in state i.
Theorem 5.1. In the fully observable fluid queue with two types of parallel customers and incomplete fault, if the fluid observes the state 
$\left( {x_j}(i),i \right)$ and decides to join, the maximal entrance fee per unit time can be expressed as
Case I. 
${\lambda}_{2} \gt \mu_{1}$.
\begin{equation*}Z\left( {x_{e1}^*(1),x_{e2}^*(1)} \right) = {\lambda _1}\left( {{R_1} - \frac{{{C_1}x_{e1}^*(1)}}{{{\mu_1}}}} \right) + {\lambda _2}\left[ {1 - {P_1}(x_{e2}^*(1))} \right]\left( {{R_2} - \frac{{{C_2}x_{e2}^*(1)}}{{{\mu_1}}}} \right),\end{equation*}where
\begin{equation*}{P_1}(x_{e2}^*(1)) = \frac{{{\pi _1}\left[ {\xi \left( {q{\lambda _1} + {\lambda _2} - {\mu_1}} \right) - \theta {\mu_0}} \right]{{\rm e}^{- \left( {\frac{\theta }{{q{\lambda _1} + {\lambda _2} - {\mu_1}}} - \frac{\xi }{{{\mu_0}}}} \right)\left( {x_{e2}^*(1) - {x_{e2}}(0)} \right)}}}}{{\xi \left( {q{\lambda _1} + {\lambda _2} - {\mu_1}} \right){{\rm e}^{- \left( {\frac{\theta }{{q{\lambda _1} + {\lambda _2} - {\mu_1}}} - \frac{\xi }{{{\mu_0}}}} \right)\left( {x_{e2}^*(1) - {x_{e2}}(0)} \right)}} - \theta {\mu_0}}}.\end{equation*} Case II. 
${\lambda}_{1} + {\lambda}_{2} \gt \mu_{1}$ and 
${\lambda}_{2} \le \mu_{1}$.
\begin{equation*}Z\left( {x_{e1}^*(1),x_{e2}^*(1)} \right) = {\lambda _1}[1 - {P_1}(x_{e1}^*(1))]\left( {{R_1} - \frac{{{C_1}x_{e1}^*(1)}}{{{\mu_1}}}} \right) + {\lambda _2}\left( {{R_2} - \frac{{{C_2}x_{e2}^*(1)}}{{{\mu_1}}}} \right),\end{equation*}where
\begin{equation*}{P_1}(x_{e1}^*(1)) = \frac{{{\pi _1}\left[ {\xi \left( {{\lambda _1} + {\lambda _2} - {\mu_1}} \right) - \theta {\mu_0}} \right]{{\rm e}^{- \left( {\frac{\theta }{{{\lambda _1} + {\lambda _2} - {\mu_1}}} - \frac{\xi }{{{\mu_0}}}} \right)\left( {x_{e1}^*(1) - {x_{e2}}(0)} \right)}}}}{{\xi \left( {{\lambda _1} + {\lambda _2} - {\mu_1}} \right){{\rm e}^{- \left( {\frac{\theta }{{{\lambda _1} + {\lambda _2} - {\mu_1}}} - \frac{\xi }{{{\mu_0}}}} \right)\left( {x_{e1}^*(1) - {x_{e2}}(0)} \right)}} - \theta {\mu_0}}}.\end{equation*} Case III. 
${\lambda _{\rm{1}}}{\rm{+ }}{\lambda _{\rm{2}}}{\rm{= }}{\mu_{\rm{1}}}$.
\begin{equation*}Z\left( {x_{e1}^*(1),x_{e2}^*(1)} \right) = {\lambda _1}\left( {{R_1} - \frac{{{C_1}x_{e1}^*(1)}}{{{\mu_1}}}} \right) + {\lambda _2}\left( {{R_2} - \frac{{{C_2}x_{e2}^*(1)}}{{{\mu_1}}}} \right).\end{equation*} Case IV. 
${\mu_{\rm{0}}} \lt {\lambda _1} + {\lambda _2} \lt {\mu_1}$.
\begin{equation*}Z\left( {x_{e1}^*(0),x_{e2}^*(0)} \right) = {\lambda _1}\left( {{R_1} - \frac{{{C_1}x_{e1}^*(0)}}{{{\mu_0}}}} \right) + {\lambda _2}[1 - {P_0}(x_{e2}^*(0))]\left( {{R_2} - \frac{{{C_2}x_{e2}^*(0)}}{{{\mu_0}}}} \right),\end{equation*}
\begin{equation*}Z\left( {x_{e1}^*(1),x_{e2}^*(1)} \right) = {\lambda _1}\left( {{R_1} - \frac{{{C_1}x_{e1}^*(1)}}{{{\mu_1}}}} \right) + {\lambda _2}\left( {{R_2} - \frac{{{C_2}x_{e2}^*(1)}}{{{\mu_1}}}} \right),\end{equation*}where
\begin{equation*}{P_0}(x_{e2}^*(0)) = \frac{{{\pi _0}\left[ {\theta \left( {p{\lambda _1} + {\lambda _2} - {\mu_0}} \right) + \xi \left( {{\lambda _1} + {\lambda _2} - {\mu_1}} \right)} \right]{{\rm e}^{- \left( {\frac{\xi }{{p{\lambda _1} + {\lambda _2} - {\mu_0}}} + \frac{\theta }{{{\lambda _1} + {\lambda _2} - {\mu_1}}}} \right)x_{e2}^*(0)}}}}{{\theta \left( {p{\lambda _1} + {\lambda _2} - {\mu_0}} \right){{\rm e}^{- \left( {\frac{\xi }{{p{\lambda _1} + {\lambda _2} - {\mu_0}}} + \frac{\theta }{{{\lambda _1} + {\lambda _2} - {\mu_1}}}} \right)x_{e2}^*(0)}} + \xi \left( {{\lambda _1} + {\lambda _2} - {\mu_1}} \right)}}.\end{equation*}Proof. By comparing the arrival rate and the outflow rate, the entrance fee per unit time can be obtained in the following four cases.
Case I. 
${\lambda}_{2} \gt \mu_{1}$.
When the buffer is in state 0, the fluid level is always greater than the balking thresholds of the customers, then the fluid is reluctant to flow into the system. When the buffer is in state 1, the entrance fee benefit shows a strictly growing trend with the increase of the fluid level. Therefore, the total entrance fees of the service provider can maximize when the fluid level reaches its maximum. According to (14) and (15), the maximal unit price of the entrance fee is 
${{R}_{j}}-{{C}_{j}}\frac{x_{ej}^{*}(1)}{{{\mu }_{1}}}$. The fluid flows into the buffer when the fluid level is less than the threshold strategy 
$x_{ej}^{*}(1),j=1,2$, then the inflow probability per unit time is 
${{\lambda }_{j}}P\left( x \lt x_{ej}^{*}(1) \right)=\lambda_j[1-{{P}_{1}}(x_{ej}^{*}(1))],j=1,2$. According to the steady-state probability distribution in Theorem 4.1, the function 
$F_1(x)$ is continuous when the fluid level reaches the threshold 
$x_{e1}^{*}(1)$, then the transient probability is 
${{P}_{1}}(x_{e1}^{*}(1))=0$. Therefore, the maximal income of service providers per unit time is 
$Z\left( x_{e1}^{*}(1),x_{e2}^{*}(1) \right)={{{\lambda }_{1}}{{p}_{1}}}+{{\lambda }_{2}}[1-{{P}_{1}}(x_{e2}^{*}(1))]{{p}_{2}}$.
Case II. 
${\lambda}_{1} + {\lambda}_{2} \gt \mu_{1}$ and 
${\lambda}_{2} \le \mu_{1}$.
When the buffer is in state 0, the fluid level is greater than the thresholds of two types of customers, and the fluid chooses to balk. When the buffer is in state 1, the maximal unit price of the entrance fee is 
${{R}_{j}}-{{C}_{j}}\frac{x_{ej}^{*}(1)}{{{\mu }_{1}}}$, and the inflow probability of the type 1 customers is 
$1-{{P}_{1}}(x_{e1}^{*}(1))$. The type 2 customers flow into the buffer when the fluid level is less than the threshold 
$x^{*}_{e2}(1)$, then the maximal income of service providers per unit time is 
$Z\left( x_{e1}^{*}(1),x_{e2}^{*}(1) \right)={{{\lambda }_{1}}[1-{{P}_{1}}(x_{e1}^{*}(1))]{{p}_{1}}}+{{\lambda }_{2}}{{p}_{2}}$.
Case III. 
${\lambda}_{1}{+}{\lambda}_{2}=\mu_{1}$.
When the buffer is in state 0, the fluid chooses to balk and the inflow probability is 0. When the buffer is in state 1, the probability that the fluid level reaches the threshold 
$x_{ej}^{*}(1)$ is 
${{P}_{1}}(x_{ej}^{*}(1))=0$, then we get the maximal income 
$Z\left( x_{e1}^{*}(1),x_{e2}^{*}(1) \right)={{\lambda }_{1}}{{p}_{1}}+{{\lambda }_{2}}{{p}_{2}}$.
Case IV. 
$\mu_{0} \lt {\lambda}_{1} + {\lambda}_{2} \lt \mu_{1}$.
When the buffer is in state 0, the maximal unit price of the entrance fee is 
${{R}_{j}}-{{C}_{j}}\frac{x_{ej}^{*}(0)}{{{\mu }_{0}}}$. The fluid flows into the buffer when the fluid level is less than the threshold strategy 
$x_{ej}^{*}(0),j=1,2$, then the inflow probability per unit time is 
${{\lambda }_{j}}P\left( x \lt x_{ej}^{*}(0) \right)=\lambda_j[1-{{P}_{0}}(x_{ej}^{*}(0))],j=1,2$. According to the steady-state probability distribution in Theorem 4.1, the function 
$F_0(x)$ is continuous when the fluid level reaches the threshold 
$x_{e1}^{*}(0)$, then the transient probability is 
${{P}_{0}}(x_{e1}^{*}(0))=0$. Therefore, the maximal income of service providers per unit time is 
$Z\left( x_{e1}^{*}(0),x_{e2}^{*}(0) \right)={\lambda _1}{p_1} + {\lambda _2}[1 - {P_0}(x_{e2}^*(0))]{p_2}$. When the buffer is in state 1, the probability that the fluid level reaches the threshold 
$x_{ej}^{*}(1)$ is 
${{P}_{1}}(x_{ej}^{*}(1))=0$, and the maximal income per unit time is 
$Z\left( x_{e1}^{*}(1),x_{e2}^{*}(1) \right)={{\lambda }_{1}}{{p}_{1}}+{{\lambda }_{2}}{{p}_{2}}$.
Next, we present some numerical examples intuitively with case IV due to the complexity of the maximal income per unit time in Theorem 5.1. We assume that 
${\lambda _1}{\rm{= 2}},\ {\lambda _2}{\rm{= 3}},\ \theta {\rm{= 0}}{\rm{.6,}}\ \xi {\rm{= 0}}{\rm{.4,}}\ {\mu_1}{\rm{= 6,}}\ {{{R}}_1}{\rm{=15,}}\ {{{R}}_2}{\rm{= 25,}}\ {{{C}}_1}{\rm{= 2,}}\ C_2=3, \alpha {\rm{= 0}}{\rm{.2}}$ in the following cases.
Figure 5. Entrance fee income 
$Z\left( x_{e1}^{*}(0),x_{e2}^{*}(0) \right)$ versus thresholds 
$x_{ej}^{*}(0)$ when 
$\mu_0=1.5$.
Figure 6. Entrance fee income 
$Z\left( x_{e1}^{*}(0),x_{e2}^{*}(0) \right)$ versus thresholds 
$x_{ej}^{*}(0)$ when 
$\mu_0=1.9$.
When the buffer is in state 1, 
$Z\left( {x_{e1}^*(1),x_{e2}^*(1)} \right)$ decreases monotonically with the thresholds 
$x_{ej}^{*}(1), j=1,2$ according to the expression of the maximal income in Theorem 5.1. Figure 5 displays the sensitivity between the maximal income and the thresholds 
$x_{ej}^{*}(0),j=1,2$ in state 0 with parameter 
${{\mu }_{0}}=1.5$. 
$Z\left( x_{e1}^{*}(0),x_{e2}^{*}(0) \right)$ first increases and then decreases with the threshold 
$x_{e2}^{*}(0)$, and always decreases with the threshold 
$x_{e1}^{*}(0)$. The income of the entrance fee is concave and the optimal thresholds 
$x_{e1}^{*}(0)$ and 
$x_{e2}^{*}(0)$ can be set to 0 and 2.6 to maximize the entrance fee 
$Z\left( x_{e1}^{*}(0),x_{e2}^{*}(0) \right)=62$. Figure 6 shows the variation between the maximal income and the optimal thresholds 
$x_{ej}^{*}(0)$ in state 0 with parameter 
${{\mu }_{0}}=1.9$. The overall trend of 
$Z\left( x_{e1}^{*}(0),x_{e2}^{*}(0) \right)$ is similar to Figure 5, and the entrance fee varies significantly when the threshold 
$x_{e1}^{*}(0)$ is relatively small.
From Figures 5 and 6, we can observe that the maximal income per unit time increases with faster service rates and more frequent deliveries, which has a positive impact on society under the fully observable case. This suggests that the social planners can appropriately accelerate the service rate and reasonably control the optimal thresholds for the fluid to gain more entrance fees. However, if the fee-collecting organization imposes excessively expensive fees on customers, the arrivals are reluctant to join the system due to the greater waiting costs, which affects society negatively and cannot achieve the global optimum.
6. Conclusion
Based on the queueing theory and non-cooperative game theory, this paper explores the equilibrium strategies in a fluid model with incomplete fault and parallel arrivals. The existence and uniqueness of the strategic behavior are derived, and an exponential utility function is constructed in the fully observable case. The entrance fee strategy for the arrivals and social designers could be imposed on the fluid and regulated the system parameters dynamically, which can maximize the social benefits without compromising individual interests. Finally, the effects of the inflow and outflow rates on the expected social benefit and the entrance fee income are illustrated by several numerical examples. The research results provide a feasible method and valuable insights for computer technologies, digital communication networks, and flexible manufacturing systems. Further extension of this work may explore the equilibrium behavior of the fluid in the unobservable cases. The economic analysis of the fluid queue with two types of priority customers and incomplete fault is also an intriguing and challenging direction.
Funding statement
This research was supported by the National Natural Science Foundation of China under grant no. 62171143.
Competing interest
The authors declare no conflict of interest.
Open access
