# Enhancing the Supersonic Wind Tunnel Performance Based on Plenum Temperature Control

Authors: A. Nazarian Shahrbabaki; M. Bazazzadeh; A. Shahriari; M. Dehghan Manshadi
Journal: ISRN Aerospace Engineering
 (2014)

---

## Abstract

The application of fuzzy logic controllers (FLCs) to the control of nonlinear processes, typically controlled by a human operator, is a topic of much study. In this paper, the design and application of a FLC is discussed to control the plenum chamber temperature for a blowdown supersonic wind tunnel (BSWT) with the aim of achieving the accurate and desired results. In this regard, first, a nonlinear mathematical model of special BSWT is developed in Matlab/Simulink software environment. Next, an artificially intelligent controller is designed using fuzzy logic approach. For this purpose, a proportional-derivative FLC (PD-FLC) system is developed in the Simulink toolbox to control the plenum stagnation temperature using a heater upstream of the plenum chamber. Finally, the system simulation results inside of the temperature and pressure controllers in comparison with the experimental run are presented. The results for Mach 2.5 blowdown run show the great performance of the Wind Tunnel Simulator Model and its temperature controller system.


---

## Body

## 1. Introduction

Blowdown supersonic wind tunnels (BSWTs) deliver flow at constant stagnation temperature and pressure. The stagnation temperature is generally regarded to be equal to the plenum temperature which is controlled by heater upstream of the plenum chamber. During a blowdown run, the storage tank temperature and pressure that supplies plenum chamber flow decrease continuously. Thus, to maintain a constant plenum pressure as close as possible to a setpoint pressure signal, the regulator or control valve must open progressively. Besides, to maintain a constant stagnation temperature in the plenum chamber, a heater must operate continuously during a supersonic run [1–3]. The controller must operate at different stagnation pressures and Mach numbers and has to be robust to accommodate the varying pressure and mass flow requirements safely. New concepts for control are under implementation with the goal of reducing transition phase and overall loads on the models [4, 5].The block diagram of the control systems in the proposed BSWT is illustrated in Figure1.Figure 1
Schematic diagram of the pressure and temperature control systems in the supersonic wind tunnel.Fuzzy logic has been the area of heated debate and much controversy during the last decades. The first paper in fuzzy set theory, which is now considered to be the seminar paper of the subject, was written by Zadeh et al. [6], who is considered the founding father of the field. In that work, Zadeh was implicitly advancing the concept of approximate human reasoning to make effective decisions on the basis of the available imprecise, linguistic information.In the 1970s, King and Mamdani [7] studied the application of FLCs to the control of nonlinear industrial processes that typically can only be controlled successfully by a human operator. The idea of FLCs has become a common solution in recent years, with applications ranging from automation of industrial processes to control of electronic devices in consumer products. The design and application of a FLC for the control of plenum chamber temperature in several wind tunnels at NASA’s Langley Research Center (LaRC) in Hampton, VA, is described [8].The intelligent controlling approaches like fuzzy logic (FL) will provide the required scope for wind tunnels to be more efficient, safe, and economic. The approaches will help to enable a level of performance that far exceeds that of today’s wind tunnel in terms of reduction of harmful emissions, maximization of run time, and minimization of noise, while improving system affordability and safety.In this paper, a nonlinear mathematical model of a special BSWT is developed in Simulink software simulator. Next, a PD-FLC system is developed in the Simulink toolbox to control the plenum stagnation temperature. Finally, the simulation results for a Mach 2.5 blowdown run in comparison with the experimental results are presented along.


## 2. Wind Tunnel Mathematical Model

The schematic diagram of the proposed BSWT is shown in Figure2. In the figure, the major subsystems are a storage tank, a gate valve, a pressure regulator valve or control valve, a heater, a plenum chamber, a supersonic nozzle, a test section, a diffuser, and a silencer.Figure 2
Schematic diagram of the blowdown supersonic wind tunnel.


### 2.1. Heater

There are two ways in which BSWTs are customary operated: (1) operation with constant pressure and (2) operation with constant-mass flow rate [8].For constant pressure operation, the only control necessary is a pressure control valve that holds the stagnation pressure constant. By the way, for constant-mass runs, the stagnation temperature must be held constant and either a heater or a thermal mass external to the storage tank is required. Even more important, the constant temperature of the constant-mass run keeps the Reynolds number constant in the test section. This drop in stagnation temperature affects the Mach number in the test section only secondarily through a change in boundary layer thickness, but it can affect balance strain gage readings significantly, and it does change the Reynolds number during a run.Obviously a heater downstream of the storage tank could be designed to yield essentially zero temperature drops, and such heaters are sometimes employed. In this regard, the controlled flow stagnation temperature in the heater output is presented as [9]
(1)
T
H
=
T
in
+
[
W
th
q
v
C
p
]
,
where T
H is the flow stagnation temperature in the heater output (
K
), T
in is the air initial stagnation temperature (
K
), W
th is the input thermal power (
W
), q
v is the mass flow rate through the valve (
kg
·
s
-
1
), and C
p is the gas constant pressure specific heat (
J
·
k
g
-
1
·
K
-
1
).


### 2.2. Control Valve

BSWTs are generally operated with a constant stagnation pressure in the plenum chamber, with control usually provided by one or more pressure control valves. In this regard, a main step of the control valve simulation is calculating the mass flow rate through the valve. This parameter is determined using the standard valve sizing relation as [10, 11]
(2)
q
V
=
C
v
N
8
F
P
Y
P
T
(
X
M
a
)
(
T
T
Z
)
.In (2), C
v is the valve sizing coefficient, N
8 is a constant for units, F
P is a correction factor that accounts for pressure losses due to piping fittings and elbows, Y is expansion factor, P
T is the tank stagnation pressure, X is the valve pressure drop ratio, M
a is the air molecular weight, T
T is the tank stagnation temperature, and Z is the air compressibility factor.It should be noted thatC
v is the functions of control valve opening and should be provided by the manufacturer of the control valve. In this paper, for the valve opening simulation during the tunnel run, some lookup tables were set to the valve manufacturer’s data.


## 2.1. Heater

There are two ways in which BSWTs are customary operated: (1) operation with constant pressure and (2) operation with constant-mass flow rate [8].For constant pressure operation, the only control necessary is a pressure control valve that holds the stagnation pressure constant. By the way, for constant-mass runs, the stagnation temperature must be held constant and either a heater or a thermal mass external to the storage tank is required. Even more important, the constant temperature of the constant-mass run keeps the Reynolds number constant in the test section. This drop in stagnation temperature affects the Mach number in the test section only secondarily through a change in boundary layer thickness, but it can affect balance strain gage readings significantly, and it does change the Reynolds number during a run.Obviously a heater downstream of the storage tank could be designed to yield essentially zero temperature drops, and such heaters are sometimes employed. In this regard, the controlled flow stagnation temperature in the heater output is presented as [9]
(1)
T
H
=
T
in
+
[
W
th
q
v
C
p
]
,
where T
H is the flow stagnation temperature in the heater output (
K
), T
in is the air initial stagnation temperature (
K
), W
th is the input thermal power (
W
), q
v is the mass flow rate through the valve (
kg
·
s
-
1
), and C
p is the gas constant pressure specific heat (
J
·
k
g
-
1
·
K
-
1
).


## 2.2. Control Valve

BSWTs are generally operated with a constant stagnation pressure in the plenum chamber, with control usually provided by one or more pressure control valves. In this regard, a main step of the control valve simulation is calculating the mass flow rate through the valve. This parameter is determined using the standard valve sizing relation as [10, 11]
(2)
q
V
=
C
v
N
8
F
P
Y
P
T
(
X
M
a
)
(
T
T
Z
)
.In (2), C
v is the valve sizing coefficient, N
8 is a constant for units, F
P is a correction factor that accounts for pressure losses due to piping fittings and elbows, Y is expansion factor, P
T is the tank stagnation pressure, X is the valve pressure drop ratio, M
a is the air molecular weight, T
T is the tank stagnation temperature, and Z is the air compressibility factor.It should be noted thatC
v is the functions of control valve opening and should be provided by the manufacturer of the control valve. In this paper, for the valve opening simulation during the tunnel run, some lookup tables were set to the valve manufacturer’s data.


## 3. Wind Tunnel Simulator Model (WTSM)

### 3.1. Facility and Requirements

The tunnel has a Mach number range of 1.5 to 3, a Reynolds number range of 25~100 × 106 per meter, and up to 70 seconds of usable run time. Also the tunnel is operated using a 12-inch diameter, pneumatically driven automatic butterfly valve (Bray control valve). The storage tank has a volume of 220 m3 and it is filled with compressed, dry air at 1.2 MPa before the run. The heaters for the BSWT system are electric heaters requiring three-phase AC voltage input which is controlled by a silicon controlled rectifier (SCR). The air flowing into the heaters passes through a bundle of electrical resistance heated tubes to maximize air contact with heated metal.


### 3.2. WTSM Development

Mathematical models of the wind tunnel process are very complicated because they involve viscous effects and distributed characteristics. In this paper, a nonlinear mathematical model of a BSWT based on isentropic relation was developed in Simulink software environment as shown in Figure3.Figure 3
Wind tunnel simulator model and pressure controller development in Matlab/Simulink software environment.The wind tunnel model is constructed with a component approach for ease of modification and replacement with different tunnel components. The WTSM includes some component modules such as storage tank, control valve, heater, plenum pressure, nozzle, and test section.


## 3.1. Facility and Requirements

The tunnel has a Mach number range of 1.5 to 3, a Reynolds number range of 25~100 × 106 per meter, and up to 70 seconds of usable run time. Also the tunnel is operated using a 12-inch diameter, pneumatically driven automatic butterfly valve (Bray control valve). The storage tank has a volume of 220 m3 and it is filled with compressed, dry air at 1.2 MPa before the run. The heaters for the BSWT system are electric heaters requiring three-phase AC voltage input which is controlled by a silicon controlled rectifier (SCR). The air flowing into the heaters passes through a bundle of electrical resistance heated tubes to maximize air contact with heated metal.


## 3.2. WTSM Development

Mathematical models of the wind tunnel process are very complicated because they involve viscous effects and distributed characteristics. In this paper, a nonlinear mathematical model of a BSWT based on isentropic relation was developed in Simulink software environment as shown in Figure3.Figure 3
Wind tunnel simulator model and pressure controller development in Matlab/Simulink software environment.The wind tunnel model is constructed with a component approach for ease of modification and replacement with different tunnel components. The WTSM includes some component modules such as storage tank, control valve, heater, plenum pressure, nozzle, and test section.


## 4. Fuzzy Logic Controller for the Plenum Temperature

The control objective is to provide the plenum chamber temperature desired by the researcher performing tests in the wind tunnel. The application of proportional-integral-derivative (PID) control loops to the temperature processes in the system results in unacceptable control with respect to operation within system temperature constraints and regulation of the controlled process at the desired temperature. A PID control loop is unable to reasonably handle the long delay time in temperature response associated with air flow through long lengths of metal piping which absorb thermal energy. The recent application of a FLC provides appropriate control of the heaters to obtain the desired plenum chamber temperature.A FLC utilizes FL to convert linguistic information based on expert knowledge into an automatic control strategy. In order to use the FL for control purposes, a front-end “fuzzifier” and a rear-end “defuzzifier” are added to the usual input-output data set. A FLC commonly consists of four sections: rules, fuzzifier, inference engine, and defuzzifier. Once the rule has been established, the controller can be considered as a nonlinear mapping from the input to the output [12]. The block diagram of the generalized indistinct controller is presented in Figure 4.Figure 4
The block diagram of the fuzzy logic controller.The designed FLC would determine the amount of heater thermal power over its transient operation. The FLC and the available WTSM constitute a closed loop which is shown in Figure5.Figure 5
Layout of the wind tunnel model and the temperature controller.The first input of the fuzzy controller is the temperature error (Error). The second input is the first input variations over time (Error Rate). The defined membership functions (MFs) for the inputs are shown in Figure6. The FLC output is heater power over time. The defined fuzzy MFs for the FLC output are shown in Figure 7.The defined membership functions for (a) Error and (b) Error Rate.
(a)
(b)Figure 7
Heater power rate membership functions (FLC output).Conventionally, fuzzy rules are established by a combination of knowledge, experience, and observation and may thus not be optimal. Additionally, in spite of efforts to formalize a development standard for fuzzy controllers, fine tuning its performance is still a matter of trial and error [13, 14].The fuzzy rules list is presented in Table1. For instance, one of the table entries is the equivalent of “if Error is Positive-Large and Error Rate is Zero then Heater Power Rate is Positive Large.”Table 1
The fuzzy rules list.
Heater power rate
Error Rate
NeL
NeS
Zero
PoS
PoL
Error
NeL
Zero
NeM
NeL
NeL
NeL
NeM
PoL
Zero
NeS
NeL
NeL
NeS
PoM
PoS
Zero
NeM
NeL
Zero
PoM
PoS
Zero
NeS
NeM
PoS
PoL
PoM
Zero
NeS
NeM
PoM
PoL
PoL
PoS
Zero
NeL
PoL
PoL
PoL
PoL
PoM
ZeroThe defuzzification process takes place after the generation of the fuzzy control signals is completed using the inference mechanism. The resulting fuzzy set must be converted to a quantity which would be sent to the process regulating valve as a control signal. In this part, the inference results of all activated logic rules are synthesized into crisp output for making a decision. In this research study, the logic AND has been implemented with the minimum operator, and the defuzzification method is based on bisector area. The variation of heater power rate versus the two controller inputs is depicted in Figure8.Figure 8
Control surface output of the FLC.By designing the controller, the influence of controller parameters on wind tunnel should be examined in more detail. So it would be essential to simulate the wind tunnel and controller simultaneously.


## 5. Simulation Results

This paper is focused on the tunnel performance for Mach 2.5 blowdown run. The tunnel operational map (Figure9), which is based on the tests data, is used to select the proper pressure controller setpoint for the proposed Mach number in the test section. It is previously noted that the stagnation temperature drops according to the polytropic process in the storage tank. In this regard, Figure 10 shows the temperature drops during the wind tunnel run without heater for a Mach 2.5 run [1].Figure 9
The tests range for the proposed BSWT.Figure 10
Temperature variation for a Mach 2.5 run.The wind tunnel simulation results with temperature controller for Mach 2.5 are presented in Figure11. It should be noted that the reference temperature for heater controller is considered to be equal to 300 K. Also experimental verification has been carried out for supersonic conditions to show the accuracy of the fuzzy logic applied to the BSWT.Simulation results for a Mach 2.5 run: (a) plenum and storage tank pressures with experimental results and (b) Reynolds number in the test section.
(a)
(b)Figure11(a) shows plenum and tank stagnation pressures for a Mach 2.5 run. Also the figure shows the pressure variation in comparison with experimental run. By the way, Figure 11(b) shows that using heater is caused to maintain the Reynolds number in the test section during the tunnel operation.Also the valve opening and flow rate variations for a Mach 2.5 blowdown run are presented in Figures12(a) and 12(b). According to Figure 12(a), after 54 seconds of tunnel operation, the control valve opening reaches its maximum opening condition. In this situation, the plenum pressure and Reynolds number are dropped suddenly.Simulation results for a Mach 2.5 run: (a) valve opening angle (VOA) function, (b) valve flow rate, (c) plenum stagnation temperature, and (d) heater power variation during the tunnel run.
(a)
(b)
(c)
(d)Plenum temperature and heater power variations are presented in Figures12(c) and 12(d). According to the figures, plenum temperature is kept as close as the reference temperature (300 K) with the lowest overshoot, settling time, and steady state error.


## 6. Conclusion

In this paper, first, a nonlinear mathematical model of special supersonic wind tunnel is developed in Simulink software environment. Then, a fuzzy logic controller system is developed to control the plenum stagnation temperature. Finally, the simulation results for Mach 2.5 blowdown run in comparison with the experimental results are presented. The stagnation temperature affects the Reynolds number and the Mach number of a blowdown supersonic test. For example, in the proposed wind tunnel, the Reynolds number increases 13.16 percent from30.8
×
10
6
/
m at Mach 2.5 when the stagnation temperature drops from 290 to 260 K. In this paper, the temperature drop is suppressed by installing heater and temperature controller upstream of the plenum chamber to maintain a constant Reynolds number with keeping the stagnation pressure constant. In this way, the Error in the heater performance during the run is minimized with the control program for the wind tunnel facility characteristics. The new intelligent control system leads to a control response with the lowest overshoot, settling time, and steady state error.



---

*Source: 317049-2014-04-02.xml*