Active Disturbance Rejection Control of a Longitudinal Tunnel Ventilation System

This paper proposes an innovative approach for controlling pollutant release in a long-distance tunnel via longitudinal ventilation. Enhanced by an active disturbance rejection control (ADRC) method, a ventilation controller is developed to regulate the forced air ventilation in a road tunnel. As a result, the pollutants (particulate matter and carbon monoxide) are reduced by actively regulating the air flow rate through the tunnel. The key contribution of this study lies in the development of an extended state observer that can track the system disturbance and provide the system with compensation via a nonlinear state feedback controller equipped by the ADRC. The proposed method enhances the disturbance attenuation capability in the ventilation system and keeps the pollutant concentration within the legitimate limit in the tunnel. In addition to providing a safe and clean environment for passengers, the improved tunnel ventilation can also achieve better energy saving as the air flow rate is optimized.


Introduction
Longitudinal ventilation systems are widely used in long-distance tunnels for air ventilation. Generally, moving vehicles in tunnels bring about air circulation from the entrance to exit. However, when the traffic is heavy and the meteorological pressure differences are low, there is a need to operate tunnel ventilators to keep pollutant concentrations below the relevant limit. Thus, there is a balance to strike between the minimal use of electricity for fan ventilators and high air quality in the tunnel because of the existence random and chaotic factors. The pollutants in the air are randomly ejected in the tunnel by incoming vehicles whilst the wind load varies in time and direction. Uncertainty and nonlinear factors make it difficult to accurately model the aerodynamics of the air and the diffusion of the pollutants within the tunnel. Therefore, a model-based control scheme cannot perform effectively in these cases.
In the literature, fuzzy logic [1] is introduced to deal with the nonlinearity and complexity of the ventilation system, and applied in the actual tunnel. The results show that energy saving is achieved whilst improving the tunnel air quality. Fuzzy control has become a popular option in ventilation control schemes. In order to improve the performance of fuzzy logic control, many scholars have developed a number of technologies. Koyama et al. [2] used the minimum power consumption of the ventilation equipment as the optimization target and solved the problems using nonlinear programming with fuzzy control. As a result, an optimal scheme to determine the number of running ventilators was developed [2]. Chu et al. developed a fuzzy controller enhanced by a genetic algorithm to maintain the concentration of pollutants at a certain level and to minimize power consumption [3]. Bogdan et al. combined model prediction with fuzzy control [4]. By measuring traffic intensity, weather conditions Parameter uncertainties or model errors in line with these two subsystem models are corrected by a nonlinear controller equipped with ADRC. By tracking the reference, an extended state observer (ESO) compensates for the inaccuracies and improves the control performance. The configuration of the ventilation control system is depicted in Figure 1. It consists of a plant (a tunnel system), an air flow velocity controller, two pollutant concentration observers, and a number of measurement devices.

Dynamical Model of a General Tunnel System with Nonlinear Method
A general tunnel system is very complex and no simple model can be designed for it. To handle the system more easily, Kurka et al. and Ferkl and Meinsma, performed a functional decomposition [19,20]; therefore, the tunnel model comprises three main functional parts (or subsystems)-ventilation, traffic and exhaust.

Ventilation model based on generalized Bernoulli equation
Consider a one-dimensional model with lumped parameters, meaning the airflow velocity v a does not depend on the y and z coordinates, i.e., ∂v a ∂y = ∂v a ∂z = 0. For simplification purposes, consider the incompressible flow, meaning the density of air ρ is constant along the entire length of the tunnel and does not depend on time.
A road tunnel can be divided into sections with constant geometry, cross-section area, hydraulic diameter, etc. Based on the assumptions stated above, the airflow velocity is constant within the whole section. The airflow dynamics in each section is described by the generalized Bernoulli equation.
The generalized Bernoulli equation for incompressible and unsteady flow can be expressed in different ways. For this purpose, the following equation is used [21,22]: where dv a dt represents the acceleration of airflow, L is the length of the given section of the tunnel, ρ is the density of air, v a is the airflow velocity, and ∆P denotes the total pressure change in the given section of the tunnel.
The total pressure change ∆P in the given section of the tunnel can be divided into pressure losses caused by air friction ∆P f ric , changes in tunnel geometry ∆P area , pressure change due to running jet fans ∆P JF , piston effect force of vehicles ∆P pist , and disturbance effects ∆P wind caused by natural wind: The pressure drop ∆P area caused by local losses depends on a loss coefficient ζ, which represents losses due to the cross section changes, the shape of transition, the direction of flow, etc., as in Equation (3): where ∆P area mainly undergoes two types of losses caused by cross section changes. One kind is the contraction loss at the entrance ∆P in = ρ 2 ·ζ in ·v 2 a ; the other is expansion loss at the outlet ∆P out = ρ 2 ·ζ out ·v 2 a ·ζ in denotes the contraction head loss coefficient at the tunnel entrance, and ζ out describes the expansion head loss coefficient at the outlet of the tunnel.
The pressure drop ∆P f ric caused by friction losses depends on wall roughness, traffic signs on the walls, etc. ∆P f ric can be calculated as: where λ represents the coefficient friction of the tunnel walls, and D H is the hydraulic diameter of the tunnel. The pressure drop ∆P pist caused by the piston effect of the vehicles represents the influence of the vehicles passing through the tunnel, shown by Equation (5): where v veh represents the average speed of the vehicles, A Tun describes the tunnel cross section area, S i is the average head surface area of the i th vehicle type, n i is the number of vehicles of the i th type, and γ conv defines the convoy polynomial [23]. The pressure drop ∆P JF describes the effect of the N JV jet fans on the flow velocity, shown by Equation (6): The pressure drop ∆P JF_i indexed by i is characterized by the average outlet velocity v JV_i and the factor k JV_i of each jet fan: where A JV_i is the cross sectional area of the i th ventilator, and k e_i is the correction factor taking account of the installation situation for the i th ventilator [5]. Additionally, the wind acting on the tunnel port affects the air flow in the tunnel dependent on the wind speed v wind by producing a dynamic pressure, ∆P wind , which is calculated as: The dynamic behavior of the air flow is described by the non-linear ordinary differential equation (Equation (3)) of first order. ∆P wind is commonly considered as the resistance pressure. ∆P JF and ∆P pist provide pressure rises for the air flow. ∆P pist is usually uncontrollable unless traffic control is implemented in an emergency. Therefore, what can be controlled is ∆P JF ; by adjusting the number of jet fans running, the velocity of air flow in the tunnel can be changed.

Exhaust distribution
Vehicles passing through a tunnel produce various types of poisonous gasses and small, visibility-reduced particulate matter (PM), especially in the case of heavy vehicles with diesel engines. At present, the main dilution objects of tunnel ventilation are carbon monoxide (CO), nitrogen oxides (NO x ), and PM.
where c(x, t) is the exhaust concentration, D is the diffusion coefficient, and R is exhaust produced by traffic traveling in the tunnel.
According to the paper [23], R is computed as: where R is emissions of CO and NO x (m 3 /s), and particle matter (m 2 /s); q ex is the base emission factor for HGVs (heavy goods vehicles and buses) with diesel engines, depending on average speed and road gradient for the base year 2000; f a is the influence factor for highway grade; f d is the correction factor for average vehicle density; f h is the altitude factor; f iv is the correction factor for gradient and speed; N m is the number of each vehicle type; and f m is the influence factor for vehicle gross masses.
Assuming that δt = 1s and δx = 2m, using the partial differential equation function package provided by MATLAB software, the concentration distribution of different pollutants in tunnels can be obtained.

Traffic
The influence of traffic on the phenomena within a tunnel is dual: it produces pollution and, on the other hand, reduces concentration of pollutants by the piston effect. The characteristics of traffic flow are described by Greenshield's macroscopic stream model as follows: where k tra f f ic is the vehicle density (vehicles/kilometer); k j_tra f f ic is the jam density (vehicles/kilometer); v tra f f ic is the average traffic speed (kilometer/hour); and v f _tra f f ic is the free flow speed (kilometer/hour). The system takes as input the traffic speed v tra f f ic , that is, the average speed of vehicles per hour that enter each of the tunnel entrances, and computes the distributions of the vehicles k tra f f ic in the tunnel.

Principle of the ADRC
The ADRC offers a new and inherently robust controller that requires very little information of the plant, as shown in Figure 2. It is composed of a tracking differentiator (TD), an extended state observer, and a nonlinear PD controller. For the sake of simplicity, the plant described by the first-order differential equation of Equation (13) is used as the theme problem to illustrate the ADRC-based control design.
where y is the output, u is the input signal, x is the state variable, w denotes external disturbances, b is the system parameter, and f (x, w, t) is a multivariable function of the state, external disturbance, and time.
where is the output, is the input signal, is the state variable, denotes external disturbances, is the system parameter, and ( , , ) is a multivariable function of the state, external disturbance, and time.

Tracking differentiator (TD)
The tracking differentiator (TD) can quickly track the input signal without set point jump by constructing a transient profile that the output of the plant can reasonably follow, and generate an approximate differential signal of the input signal. The discrete equation is: where r(t) is the reference input signal of the ADRC. x 1 (t) tracks the input signal r(t) and x 2 (t) is the first-order derivative of x 1 (t). T s is the discrete control period. The dynamic performance of the TD can be tuned by r and h. h is the filtering parameter. A larger h means a better filtering effect on noise. r determines the tracking effect. The larger the value of r, the faster the tracking speed and the greater the overshooting. Therefore, both parameters need to be designed according to the output profile of the TD. The optimal control function of the discrete time system is as follows:

Extended State Observer (ESO)
Considering w in Equation (13) as a disturbance, a second-order extended state observer can be proposed as: where x r (t) is the actual feedback signal of the state variable x(t); z 1 (t) is the estimate of the state variable; e(t) is the estimation error of the observer; z 2 (t) is the estimation of the total disturbance of the system; and a 1 , δ 1 , β 01 , β 02 are four adjustable parameters in the ESO. Taking 0 < a 1 < 1 for instance, the smaller the value of a 1 , the stronger the nonlinearity of function f al(e(t), a 1 , δ 1 ). This means that the ESO has better resistance to system disturbances. δ 1 is the filter parameter of the ESO. β 01 and β 02 greatly affect the dynamic performance of the system. The estimation of state variables is mainly influenced by β 01 . The estimation of system disturbance is mainly affected by β 02 . The larger the values of β 01 and β 02 , the faster the convergence is estimated. However, if they are too large, the output of the ESO might produce an oscillating divergence and a high frequency noise signal. With a well behaved observer, z 1 (t) will closely track the state variable of the system x r (t), and z 2 (t) can estimate the total disturbances f (x, w, t) accurately enough. A nonlinear combination function f al(e(t), a 1 , δ 1 ) is:

Nonlinear State Error Feedback Control Laws (NLSEF)
According to the output of the TD and ESO, the first-order nonlinear state feedback control law is constructed as: where u(t) is the control variable of the ADRC. a 01 , δ 01 , β 1 are adjustable parameters. Their meaning is similar to the adjustable parameters in the ESO. The ADRC consists of TD, ESO, and NLSEF. According to different application specifications, these three parts can be used separately or in combination.

Tunnel Ventilation Control System Based on the Active Disturbance Rejection Technique
According to the basic principle of the ADRC, this paper presents a cascade ventilation control system consisting of the pollution control loop and the airflow control loop (Figure 1). The pollution control loop Σ 1 aims to maintain the pollutant concentration below a specified limit. The control variable of airflow control loop Σ 2 is the running number of jet fans. These two control subsystems will be introduced in the following subsections.

Pollutant Concentration Control Subsystem Based on ADRC
Consider the mass balance equation (Equation (9)) mentioned above; c(x, t) is a multivariable function of both time and position, and is difficult to mathematically analyze with accuracy. Here, c(x, t) does not need to be expressively known and can be simplified by the form: where b 0_c is a certain value, and the remaining uncertainties b 0_c − ∂c ∂x ·v a are attributed to the total disturbance of the system. This is clearly an ADRC problem and should be treated with Equation (23) taken into account as a total disturbance: By combining various known and unknown quantities for an overall disturbance, a complex control problem can be converted into a simple one that the ADRC can solve. Since an ADRC controller consists of TD, ESO and NLSEF, the key calculations in line with these can be presented by Equations (24)-(26) accordingly.
Energies 2020, 13, 1871 where r c (t) is the upper limit of the pollutant concentration. c 1 (t) is the concentration estimation and c 2 (t) is the first-order derivative of c 1 (t).
where v a (t) is the control variable.

Airflow Velocity Control Subsystem Based on ADRC
Dynamic behavior of the flow velocity v a (t) described by the incompressible non-stationary Bernoulli equation (Equation (2)) can be rewritten by this form: ∆P JF_i can be represented by the following equation: where v r is the maximum airflow velocity. ∆P JF_i_1 is the known pressure drop and ∆P JF_i_2 refers to the uncertain part of the pressure drop. Similar to the calculation of ∆P JF_i , ∆P pist , ∆P f ric , and ∆P area can be written as the equations below: Energies 2020, 13, 1871 9 of 16 Therefore, Equation (24) can be rewritten by: where b 1_air and b 0_air refer to the model errors which can be regarded as part of the total disturbance. N JV is the control variable. The key parameters in the ESO and NLSEF can be calculated by Equations (36) and (37): where v ar (t) is the real airflow velocity measured by anemometers, z 1_a (t) refers to the state variable, and z 2_a (t) refers to the estimate of the total disturbances.
where v a (t) is the value of the control variable in Equation (26), which is the reference for this subsystem. By tracking the pollutant concentration curve, the number of running jet fans N JV (t) can be obtained, and the airflow velocity changed accordingly in a longitudinal ventilation tunnel.

Test Results
In this work, a range of tests were performed to optimize the control algorithms and energy consumption.

The Status of the Experimental Tunnel
A ventilation control system based on cascaded ADRC controllers was tested in the simulation model of the right tube of the Zhongliangshan tunnel on the Chengdu-Chongqing expressway in China [24]. Figure 3 shows the layout and installation scheme of the right tube. The parameters of the right tube and vehicle proportion statistics of the peak hours are listed in Tables 1 and 2,  where ( ) is the value of the control variable in Equation (26), which is the reference for this subsystem. By tracking the pollutant concentration curve, the number of running jet fans ( ) can be obtained, and the airflow velocity changed accordingly in a longitudinal ventilation tunnel.

Test Results
In this work, a range of tests were performed to optimize the control algorithms and energy consumption.

The Status of the Experimental Tunnel
A ventilation control system based on cascaded ADRC controllers was tested in the simulation model of the right tube of the Zhongliangshan tunnel on the Chengdu-Chongqing expressway in China [24]. Figure 3 shows the layout and installation scheme of the right tube. The parameters of the right tube and vehicle proportion statistics of the peak hours are listed in Tables 1 and 2, respectively.
The test results shown in this section are obtained based on Tables 1 and 2. If this ADRC scheme is applied to other tunnels or tubes, the parameters will be required to be optimized from measurements to obtain a better model accuracy.     The test results shown in this section are obtained based on Tables 1 and 2. If this ADRC scheme is applied to other tunnels or tubes, the parameters will be required to be optimized from measurements to obtain a better model accuracy.

Effectiveness of the Tunnel Model
The right tube of the experimental tunnel uses a number of jet fans to form a longitudinal ventilation system. Under these circumstances, the concentration of pollutants (CO, PM) increases from the entrance to the exit, and concentration reaches a maximum near the exit. Figure 4 shows that the distribution for CO and PM complies with this change law. These results confirm the effectiveness of the developed model. Moreover, the control result based on this model is also effective, tracking the corresponding changes accordingly.
Energies 2020, 13, x FOR PEER REVIEW 10 of 16 from the entrance to the exit, and concentration reaches a maximum near the exit. Figure 4 shows that the distribution for CO and PM complies with this change law. These results confirm the effectiveness of the developed model. Moreover, the control result based on this model is also effective, tracking the corresponding changes accordingly.

The Performance ADRC on Longitudinal Ventilation System
Among the ADRC components, the TD, ESO, and NLSEF can be used separately or together. In order to make full use of the ESO, they start to work by close monitoring CO, PM, and airflow velocity once the system launches. The TD and NLSEF are not activated until the mean concentration of CO or the maximum concentration of PM exceeds the upper limits. A flow chart of operating the system is shown in Figure 5.
The scheme starts with an initialization. Next, the traffic speed is measured and used to define the pollutant concentration limit. Three important elements, namely, the exhaust concentration, pollutant distribution, and vehicle density, are calculated according to Equations (10)- (12). Then, the pollutant concentration and airflow velocity are estimated according to Equations (25) and (36). Four cases of the combination of CO and PM concentration are considered. In each case, the desired

The Performance ADRC on Longitudinal Ventilation System
Among the ADRC components, the TD, ESO, and NLSEF can be used separately or together. In order to make full use of the ESO, they start to work by close monitoring CO, PM, and airflow velocity once the system launches. The TD and NLSEF are not activated until the mean concentration of CO or the maximum concentration of PM exceeds the upper limits. A flow chart of operating the system is shown in Figure 5.  Figure 6 shows that ESO performance is in line with CO concentration in a control period. Figure 6a illustrates the performance of CO concentration observations. The dotted line (estimated value) accurately traces the solid line (actual value) at the value of 20 ppm. Meanwhile, Figure 6b, the total disturbance of the system can be observed ( 2 = 229 ppm). The ESOs for PM concentration and airflow velocity have the same excellent performance.

Start
Calculating the numbers of fan running based on the desired ( ) Maintaining the previous airflow velocity as the current one The scheme starts with an initialization. Next, the traffic speed is measured and used to define the pollutant concentration limit. Three important elements, namely, the exhaust concentration, pollutant distribution, and vehicle density, are calculated according to Equations (10)- (12). Then, the pollutant concentration and airflow velocity are estimated according to Equations (25) and (36). Four cases of the combination of CO and PM concentration are considered. In each case, the desired airflow is defined accordingly. If CO exceeds the upper limit and PM is lower than the threshold, a desired velocity can be determined by CO emission control by ADRC, and then used. In the following step, the number of running fans is defined as per the desired airflow velocity. Then, the airflow speed is updated. This is a whole loop of the operation.

Performance Analysis of the ESO
Concentration of CO, concentration of PM, and airflow velocity are measured every 15 minutes, and sent to their corresponding ESO. Figure 6 shows that ESO performance is in line with CO concentration in a control period. Figure 6a illustrates the performance of CO concentration observations. The dotted line (estimated value) accurately traces the solid line (actual value) at the value of 20 ppm. Meanwhile, Figure 6b, the total disturbance of the system can be observed (z 2 = 229 ppm). The ESOs for PM concentration and airflow velocity have the same excellent performance.

Steady tracking performance
Continuous congestion is used to test the steady tracking performance of the designed control system. When the average speed of the traffic flow drops to 35 km/h, continuous congestion occurs. At this time, the concentrations of CO and PM both exceed the upper limit. The upper limits of the CO and PM concentrations in a congested state are 70 ppm and 0.0075 m −1 , respectively. The ADRC controller of CO concentration and PM concentration will give the required airflow velocity. Because the deviation of PM concentration from the upper limit is larger than the deviation of CO concentration, the required wind speed to dilute PM concentration is larger than that required to dilute CO concentration. A larger airflow velocity value than the output reference is selected. The ADRC controller of the air flow velocity will receive the required number of jet fans and then act on it. In order to save the power, it is sufficient that the number of jet fans in operation can reduce the concentration of pollutants to near the upper limits. Therefore, a suitable concentration of PM tracking curve (solid line) is designed in Figure 7a-c, which presents the CO concentration and air flow velocity under control by the proposed method. Figure 7d demonstrates the number of jet fans required.

Steady tracking performance
Continuous congestion is used to test the steady tracking performance of the designed control system. When the average speed of the traffic flow drops to 35 km/h, continuous congestion occurs. At this time, the concentrations of CO and PM both exceed the upper limit. The upper limits of the CO and PM concentrations in a congested state are 70 ppm and 0.0075 m −1 , respectively. The ADRC controller of CO concentration and PM concentration will give the required airflow velocity. Because the deviation of PM concentration from the upper limit is larger than the deviation of CO concentration, the required wind speed to dilute PM concentration is larger than that required to dilute CO concentration. A larger airflow velocity value than the output reference is selected. The ADRC controller of the air flow velocity will receive the required number of jet fans and then act on it. In order to save the power, it is sufficient that the number of jet fans in operation can reduce the concentration of pollutants to near the upper limits. Therefore, a suitable concentration of PM tracking curve (solid line) is designed in Figure 7a-c, which presents the CO concentration and air flow velocity under control by the proposed method. Figure 7d demonstrates the number of jet fans required.
it. In order to save the power, it is sufficient that the number of jet fans in operation can reduce the concentration of pollutants to near the upper limits. Therefore, a suitable concentration of PM tracking curve (solid line) is designed in Figure 7a-c, which presents the CO concentration and air flow velocity under control by the proposed method. Figure 7d demonstrates the number of jet fans required.

Dynamic Tracking Performance
When the state of traffic flow develops from congestion to stagnation, the concentration of pollutants continues to rise. This process can be used to test the dynamic performance of the designed control system. The upper limits of the CO and PM concentrations in a stagnation state are 70 ppm and 0.009 m −1 , respectively. Similar to the congestion state, PM pollutes the air more severely than CO, and is the main control target. This is related to the composition of vehicles in the traffic flow. As can be seen from Table 2, the higher the proportion of diesel heavy vehicles, the larger the amount of PM. Therefore, PM was the main dilution target. The tracking curve is designed as the solid line in Figure 8a. The actual PM concentration curve is shown as the dotted line in Figure 8a. As can be seen, the actual curve is able to track the reference well. The absolute error between the two is less than 0.00079 m −1 . Figure 8b reflects the change in CO concentration, which is below the upper limit. Figure 8c,d shows the air speed in the tunnel and the required number of fans running. It can be seen that 10 fans running is needed in ventilation when the traffic is at standstill.

Dynamic Tracking Performance
When the state of traffic flow develops from congestion to stagnation, the concentration of pollutants continues to rise. This process can be used to test the dynamic performance of the designed control system. The upper limits of the CO and PM concentrations in a stagnation state are 70 ppm and 0.009 m −1 , respectively. Similar to the congestion state, PM pollutes the air more severely than CO, and is the main control target. This is related to the composition of vehicles in the traffic flow. As can be seen from Table 2, the higher the proportion of diesel heavy vehicles, the larger the amount of PM. Therefore, PM was the main dilution target. The tracking curve is designed as the solid line in Figure 8a. The actual PM concentration curve is shown as the dotted line in Figure 8a. As can be seen, the actual curve is able to track the reference well. The absolute error between the two is less than 0.00079 m −1 . Figure 8b reflects the change in CO concentration, which is below the upper limit. Figure 8c,d shows the air speed in the tunnel and the required number of fans running. It can be seen that 10 fans running is needed in ventilation when the traffic is at standstill.  Figure 8a. The actual PM concentration curve is shown as the dotted line in Figure 8a. As can be seen, the actual curve is able to track the reference well. The absolute error between the two is less than 0.00079 m −1 . Figure 8b reflects the change in CO concentration, which is below the upper limit. Figure 8c,d shows the air speed in the tunnel and the required number of fans running. It can be seen that 10 fans running is needed in ventilation when the traffic is at standstill.

Energy-Saving Performance Analysis
Step control (on-off control) is usually used in ventilation systems [24]. Table 3 shows a comparison of the step control with the proposed method. It can be seen that although step control is simpler and easier to implement, the number of running jet fans is much higher than the cascaded ADRC strategy proposed in this paper. Therefore, the latter can significantly reduce electricity consumption.

Conclusions
In this paper, a cascaded longitudinal ventilation control approach based on the ADRC was presented for a road tunnel. This control system consists of a CO and PM concentration control subsystem and an air flow velocity control subsystem. This cascade control method is designed to achieve multiple control objectives. In order to reduce the concentration of pollutants in the tunnel, the forced air flow rate in the tunnel needs to be dynamically controlled. Furthermore, the number of jet fans in operation needs to be optimized to save energy. Therefore, by adopting corresponding control measures for each link, more accurate control schemes can be obtained.
Conventionally, step control is widely used in tunnel ventilation systems. However, this

Energy-Saving Performance Analysis
Step control (on-off control) is usually used in ventilation systems [24]. Table 3 shows a comparison of the step control with the proposed method. It can be seen that although step control is simpler and easier to implement, the number of running jet fans is much higher than the cascaded ADRC strategy proposed in this paper. Therefore, the latter can significantly reduce electricity consumption.

Conclusions
In this paper, a cascaded longitudinal ventilation control approach based on the ADRC was presented for a road tunnel. This control system consists of a CO and PM concentration control subsystem and an air flow velocity control subsystem. This cascade control method is designed to achieve multiple control objectives. In order to reduce the concentration of pollutants in the tunnel, the forced air flow rate in the tunnel needs to be dynamically controlled. Furthermore, the number Energies 2020, 13, 1871 15 of 16 of jet fans in operation needs to be optimized to save energy. Therefore, by adopting corresponding control measures for each link, more accurate control schemes can be obtained.
Conventionally, step control is widely used in tunnel ventilation systems. However, this method might cause large fluctuations in the air flow speed over the on/off operation to jet fans, which leads to undesirable dramatic changes in pollutant concentration. By contrast, the proposed ADRC scheme can maintain the pollutant concentration level within a permitted range with minimal energy consumption. The ADRC scheme uses a state observer to capture the external disturbance and compensates for the model uncertainty at the subsystem level to improve the accuracy of the nonlinear control scheme. In this way, the pollutant concentration can be controlled close to the upper limit while running the minimal number of fans. The test results show the effectiveness of the proposed method where the ADRC is used for pollutant control in a tunnel's longitudinal ventilation in an efficient manner.
However, the ADRC algorithms also have disadvantages. The control algorithms are more complex than step control. Moreover, they include several control parameters in the control loop that are determined by empirical methods based on historical data and trial-and-error. In further work, control parameters will be developed and optimized to improve the system performance. A data-driven method will be utilized to enhance the accuracy of NLSEF for nonlinear problems.