International Journal of Mechanical and Production

Engineering Research and Development (IJMPERD)

ISSN (P): 2249-6890; ISSN (E): 2249-8001
Vol. 9, Issue 2, Apr 2019, 679-690
© TJPRC Pvt. Ltd.

EXPERIMENTAL FLOW ANALYSIS OF AEROSPACE DUCT

OF EQUIPMENT COOLING SYSTEM IN AN AIRCRAFT

NEERAJ CHAURASIA1 & S. P. S. RAJPUT2

1
Senior Scientific Officer-I, Regional Director Aeronautical Quality Assurance (Guided Weapon & Materials),
Defence Research Development Laboratory, Hyderabad, Telangana, India
2
Professor & Head, Department of Mechanical Engineering, Maulana Azad National Institute of Technology,
Bhopal, Madhya Pradesh, India
ABSTRACT

This paper discusses about the Flow analysis of Aerospace ducts of Equipment cooling system of an aircraft.
Equipment Cooling is provided to maintain assigned temperature conditions in equipment compartment for satisfactory
working of Line Replaceable units (LRUs), which are installed inside equipment compartment. Forced air cooling is
provided to these LRUs using Air Conditioning system (ACS) air. The life span of aircraft is approximately 30-40 years,
and various up-gradations are required to be carried out to meet the operational &mission requirements. Up-gradation

Original Article
programmes of aircraft required installation of new LRUs, thus a requirement of redistribution of flow through complex
duct geometries to meet their additional cooling requirements is mandatory. In the present practice, designers carryout
the redistribution of flow of air through these complex ducts using Continuity & Bernoulli’s equation. It was
experienced that new LRUs at so many times not function properly due to inadequate cooling condition resulting
malfunctioning of LRUs. This Paper aims to identify the reason for not achieving the designed cooling condition at
designated LRUS. To analyse the flow through Aerospace ducts in transient condition, an experimental Environmental
Control System (ECS) test rig has been newly designed to simulate the actual layout of installation of Aerospace ducts
into aircraft. Theoretical analysis correlated by experimental trials analysis was carried out on Environmental Control
System Test Rig. The results analysed and concluded the reason for not achieving the required cooling to LRUs.

KEYWORDS: Equipment Cooling System, Environmental Control System, Flow analysis & Line Replaceable Unit
(LRU)

Received: Jan 05, 2019; Accepted: Jan 25, 2019; Published: Mar 23, 2019; Paper Id.: IJMPERDAPR201968

Notations

ρ = Density of air

µ = Dynamic Viscosity

ν = Kinematic Viscosity

m* =Mass flow rate

m13 = air mass flow rate in branch 13

V= Initial Velocity of the air

f = Darcy friction factor

www.tjprc.org SCOPUS Indexed Journal [email protected]

680 Neeraj Chaurasia & S. P. S. Rajput

Rh = Hydraulic radius

Dh = Hydraulic Diameter

Re = Reynolds no of flow

(hf)01 = Head loss due To Friction in section 0 to1

hf)13 = Head loss due To Friction in section 1 to3

V01 = Mean velocity between section 0 & 1

V13 = Mean velocity between section 1 & 3.

E0 = Energy at the section 0

E1 = Energy at the section 1

E3 = Energy at the section 3

P0 = Pressure drop at section 0

P1 = Pressure drop at section 1

P3 = Pressure drop at section 3

D01 = Diameter of sections 0 to 1

D13 = Diameter of sections 1 to 3

A01 = Area of sections 0 to 1

A13 = Area of sections 1 to 3

L01 = Length of sections 0 to 1

L13 = Length of sections 1 to 3

ε = Surface roughness

1. INTRODUCTION

This research paper discusses about the Theoretical & Experimental analysis of mass flow rate to maintain
assigned temperature conditions in equipment compartment for satisfactory working of Line Replaceable units (LRUs),
which are installed inside equipment compartment. Force air cooling is provided to these LRUs using Air Conditioning
system (ACS) air. The life span of the aircraft is approximately 30-40 years and during this life spans various up gradation
to be carried out to meet the operational & Mission requirements. During up gradation programmes of aircraftdue to
introduction of new LRUs, there remains a requirement of redistribution of flow through complex duct geometries to meet
their additional cooling requirements. In the present practice, designers carryout the redistribution of flow of air through
these complex ducts using Continuity & Bernoulli’s equation. But in order to avoid the malfunctioning of the LRUs, it
needs to analyse the flow experimentally under transient condition. A new Environmental Control Test Rig is designed and
experiment performed. It is also to be noted that to increase or decrease the speed of aircraft, pilots has to change the
engine throttle rating. As engine rating changes, the bleed air pressure fed to Air Conditioning System LRUs varies. This

Impact Factor (JCC): 7.6197 SCOPUS Indexed Journal NAAS Rating: 3.11
Experimental Flow Analysis of Aerospace Duct of 681
Equipment Cooling System in an Aircraft

leads to increment or decrement of turbo cooler turbine inlet pressure. Subsequently, based on differential pressure across
the turbine of turbo cooler, the mass flow rate varies. The sudden variation of this engine rating leads to variation in flow
through ducts and variation of pressure drop frequently due to variable Reynolds’s no. Hence, there is a need to study the
effects of variation of engine throttle on sudden increment and decrement of flow? Generally, it is very difficult for the
designer’s to predict the effect of variation of engine throttle on flow through complex duct geometries due to change in
transient condition.

2. LITERATURE REVIEW

To start a research project, it is very important to know the current methodologies used in industry and explore
simple and accurate methodologies for analysis. The first step for flow analysis is through these complex ducts using
Continuity & Bernoulli’s equation. Initial velocity of air is to be found out for given diameter of pipe by continuity
equation. The redistribution of flow is carried out by Bernoulli’s equation, for which Reynolds number to be found out.
The head loss is important factor for redistribution of flow. Total head loss in Pipe flow is calculated by Darcy-Weisbach

formula hf = , For head loss, friction factor to be calculated. The Darcy friction factor formulae are equations based on

experimental data and theory. Darcy-Weisbach friction factor, resistance coefficient or simply —friction factor and is four
times larger than the Fanning friction factor [1]. For laminar flow, The Darcy friction factor for laminar flow in a circular
pipe (Reynolds number less than 2320) is given by the formula: f = Transition flow(neither fully laminar nor fully

turbulent) flow occurs in the range of Reynolds numbers between 2300 and 4000. The value of the Turbulent flow in
smooth conduits: The Blasius correlation is the simplest equation for computing the Darcy friction factor. Turbulent flow
in rough conduits: The Darcy friction factor for fully turbulent flow (Reynolds number greater than 4000) in rough
conduits is given by the Colebrook equation. Colebrook–White equation: Colebrook–White equation (or Colebrook
equation) expresses the Darcy friction factor(f) as a function of Reynolds number Re and pipe relative roughness ε / Dh,
fitting the data of experimental studies of turbulent flow in smooth and rough pipes.[2][3] The equation can be used to
(iteratively) solve for the Darcy–Weisbach friction factor f. For a conduit flowing completely full of fluid at Reynolds
numbers greater than 4000, it is expressed as:

ε .
= -2*log10*( + )
√ . √

Some sources use a constant of 3.71 in the denominator for the roughness term in the first equation above. [4]
Recently, the Lambert W function has been employed to obtain explicit reformulation of the Colebrook equation. [5]

Colebrook equation can be solved by iteration using the Newton–Raphson method. Approximation of Cole brook equation
is provided by Halland eqaution, Swamee Jain Equation, Serghides’s solution and Goudar–Sonnad equation. The Haaland
equation was proposed by Norwegian Institute of Technology professor Haaland in 1984. The Haaland equation is defined
ε .
as: = -1.8*log10 [ + . )]-2.
√ .

The Swamee–Jain equation [6] is used to solve directly for the Darcy–Weisbach friction factor f for a full-flowing
circular pipe. It is an approximation of the implicit Colebrook–White equation. Swamee Jain equation is given by f=0.25
ε .
[log10 ( ( ) +( . ))]-2. Goudar equation [8] is the most accurate approximation to solve directly for the Darcy–Weisbach
.

friction factor for a full-flowing circular pipe. It is an approximation of the implicit Colebrook–White equation. Brkić
shows one approximation of the Colebrook equation based on the Lambert W-function [9]. Early approximations by Paul

www.tjprc.org SCOPUS Indexed Journal [email protected]

682 Neeraj Chaurasia & S. P. S. Rajput

Richard Heinrich Blasius in terms of the Moody friction factor [10] f = 0.316Re0.25. Johann Nikuradse in 1932 proposed that
this corresponds to a power law correlation for the fluid velocity profile. Mishra and Gupta in 1979 proposed a correction
for curved or helically coiled tubes, taking into account the equivalent curve radius. From literature survey, it is concluded
that Swamee Jain equation is latest refined available equation considering all the relevant factors. Hence, Swamee Jain
equation was selected for our application.

2.1 Equipment Cooling System during up Gradation Programme

During up gradation of aircraft, old avionics LRUs have been replaced by the new one like: UNIT 1, UNIT 2,
UNIT 3, UNIT 4, UNIT 5, UNIT 6, UNIT 7 & UNIT 8. The cooling air requirement of these LRUs along with their
positioning on aircraft is placed at Table 1.

Table 1: Cooling Air Requirement of these LRUs

along with their Positioning on Aircraft
Required Air Mass
LRUs Name
Flow Rate Through LRU
(Newly Introduced)
(Kg/hr)
UNIT 1 20
UNIT 2 20
UNIT 3 55
UNIT 4 20
UNIT 5 65
UNIT 6 164
UNIT 7 14
UNIT 8 15

Configuration of placement of New LRUs at Nose Compartment of aircraft is shown at Figure 1.

Figure 1: Configuration of Equipment Cooling System

Detailed flow analysis regarding distribution of air through various LRUs of aircraft is carried out with the help of
Continuity & Bernoulli’s equation, to ensure appropriate flow distribution in respective LRUs. Head losses due to friction
.
& bends during flow were calculated using friction factor f = 0.0032 + . .

A sample theoretical calculation based on available mass flow rate is shown below:

Air mass Flow Rate outlet of turbo cooler = 400kg/hr

Mass flow rate of cooling air supplied to Compartment 5 = 18 kg/hr

Mass flow rate of cooling air supplied to newly introduced LRUs = 400-18 = 382 Kg/hr

Impact Factor (JCC): 7.6197 SCOPUS Indexed Journal NAAS Rating: 3.11
Experimental Flow Analysis of Aerospace Duct of 683
Equipment Cooling System in an Aircraft

Absolute outlet Pressure from Turbo cooler = 1.2 kgf/cm2

Temperature of the air outlet from turbo cooler ≈150 C

Properties of air at initial temperature are taken from thermodynamic table.

Density of air ρ = 1.2kg/m3

Dynamic Viscosity µ = 1.81 x 10-5N-s/m2

Kinematic Viscosity ν = 1.51 x 10-5m2/sec

Mass flow rate m* = ρ*A*V (Applying Continuity Equation)

∏
382/3600 = 1.2* *D2*V

Taking the initial Diameter of the pipe as 36mm= 0.036m

Initial Velocity of the air V = 86.91 m/s

Assumed that there is no pressure drop through pipelines during the flow due to less pressure gradient across
turbo cooler outlet & newly introduced LRUs. The flow is assumed to be due to dynamic pressure gradient (i. e kinetic
head gradient).

To find the head loss across the pipe lines, we will calculate the friction Factor ‘f’

To calculate f, it is required calculate Reynolds no of flow:

Re = ρ*D*V/ µ

Re = 1.2*0.036*86.91/1.81* 10-5

Re = 207431.6

As Re> 4000, the flow is Turbulent.

For 105< Re < 4 x 107 using the empirical relation for friction factor f

f = 0.0032 + 0.210/ Re 0.257

f = 0.0032 + 0.210/ (207431.6) 0.257

f = 0.01223

3. CALCULATIONS FOR SECTION 0 TO 1

Energy balance equation between section 0 and 1

E0 =E1 + (hf)01

Using Bernoulli’s equation between section 0 & 1

www.tjprc.org SCOPUS Indexed Journal [email protected]

684 Neeraj Chaurasia & S. P. S. Rajput

# #
+ + " = + + " + head Losses
ρ ! ρ !

Due to less pressure gradient across turbo cooler outlet & New LRUs, itis assumed that pressure drop between
sections 0 to 1 is negligible. Thus,

= #, & Z0 = Z1
ρ ρ

L01 = 0.085 m, D01 = 0.036m, V01 =86.91 m/s

V01: mean velocity between section 0 & 1.

#
- head Losses= =E1
! !

Head loss due To Friction in section 0 to1,

(hf)01 = f * L01* V 201 / 2gD01 n

(hf)01 =0.01223*0.085*86.912 / 2*9.81*0.36

(hf)01 = 1.112 m of air

Velocity Head Available at the inlet = V201 /2g= 86.912 / 2 * 9.81= 384.98

Hence, Energy at the section 1

E1 = E0 – ( hf)01

E1 = 384.98 –1.112

E1 = 383.868 m of Air

In a similar manner, calculations for remaining sections have been undertaken and results are tabulated in
Table 2.

Table 2: Calculated Diameters of Pipelines at Respective Section Supplying ECS Air to LRUs
Energy at the Energy at
Frictional Iterative Calculated Fabricated
Mass Flow Length of Inlet of Section Outlet of
Section Loss Velocity Diameter Pipe
Rate Kg/hr Pipe (m) (Einlet) Section (E1)
m of air V (m/s) (mm) (mm)
m of air = Einlet- FL
0-1 382 0.085 385.4 1.112 383 86.91 35.9 36
1 to 3 332 0.305 383 37.98 345.01 82.2 34.5 36
1 to 2 50 0.385 383 97.59 285.4 74.8 14.03 14
2 to A 25 0.23 285.40 60.98 224.42 66.3 10.54 12
2 to B 25 0.062 285.40 20.27 265.13 72.1 10.10 12
3 to 4 291.8 0.175 345.01 20.60 324.40 79.7 34.5 34
3 to C 40.2 0.41 345.01 76.42 268.58 58.8 14.19 16
4 to 5 266.8 0.46 324.4 47.12 277.28 72.6 32.9 32
4 to D 25 0.375 324.4 67.30 257.09 56.6 11.41 10
5 to 6 232 0.14 277.3 14.33 262.94 70.6 31.12 32
6 to E 55 0.195 262.9 21.69 241.25 54.71 17.21 18
6 to F 177 0.4 262.9 24.42 238.52 54.4 30.97 32
5 to 7 34.8 0.46 277.2 81.92 195.35 60.3 13.04 12
7 to G 17 0.32 195.3 36.50 158.85 43.7 10.7 12
7 to H 17 1.045 195.3 99.19 96.16 40.6 11.1 12

Impact Factor (JCC): 7.6197 SCOPUS Indexed Journal NAAS Rating: 3.11
Experimental Flow Analysis of Aerospace Duct of 685
Equipment Cooling System in an Aircraft

3.1 Experimental Determination of Velocity of air flow at Discrete Air Outlet Points of Air
Conditioning System and Mass Flow Rate of Cooled Air Supplied to Avionics LRUs

Velocity of Air flow at discrete air outlet points of Air Conditioning System i. e. Inlet to LRUs has been
experimentally determined by designing an Environmental Control System (ECS) test rig as described below:

3.1.1 Experimental Setup of ECS Test Rig

Aerospace ducts based on theoretical calculations mentioned at Table 2 have been manufactured. These ducts
were made of Aluminum alloy. An experimental Environmental Control System (ECS) test rig has been designed to
simulate the actual layout of installation of Aerospace ducts into aircraft as shown in Figure 2.

AIR CONDITIONER
1.5 Tonne

Figure 2: ECS Test Rig Layout

Aircraft turbo cooler is used on ECS test rig to simulate the system parameters based on variable input pressure
and temperature. Thermal flow meters were installed at the duct to measure the Air flow rate supplied to particular LRU.
To record and online recording of test rig parameters, PLC based control system has been designed. Physical components
of ECS test rig are shown in Figure 3. Placement of thermal flow meters at the inlet of respective LRUs is shown in
Figure 4.

www.tjprc.org SCOPUS Indexed Journal [email protected]

686 Neeraj Chaurasia & S. P. S. Rajput

Figure 3: Components of ECS Test Rig

Figure 4: Installation of Flow Meter at the Inlet of ECS ducts

3.2 Experimentally Determination of Mass Flow Rate of Cooled Air Supplied to Avionics LRUs

Mass flow rates supplied to particular LRUs have been measured using thermal flow meters on ECS test rig
shown at Figure 2. The results are placed in Table 3

Table 3: Experimentally Measured Mass Flow Rates Supplied to Respective LRUs

Required Air Mass Theoretical Calculated Experimentally
LRUs
Flow Rate Through LRU Mass Flow Rate Through Measured Mass
Name
(Kg/hr) LRU(Kg/hr) Flow Rate(Kg/hr)
UNIT 1 20 21 32
UNIT 2 20 21 28
UNIT 3 55 57 43
UNIT 4 20 25 31
UNIT 5 65 66 55.5
UNIT 6 164 166 180
UNIT 7 14 17 16.5
UNIT 8 15 17 16

Impact Factor (JCC): 7.6197 SCOPUS Indexed Journal NAAS Rating: 3.11
Experimental Flow Analysis of Aerospace Duct of 687
Equipment Cooling System in an Aircraft

Table 4: The Comparative Deviation of Mass Flow Rates Based on Theoretical

Deviation, %
(Experimental Mass Flow Rate) Kg/hr) Theoretical w. r. t
Section (Theoretical Mass
Average Values Taken after 05 Times Experimental Mass
Flow Rate) Kg/hr
Experiment Run Flow Rate
0-1 382 381.5 0.131062
1 TO 3 332 334 -0.5988
1 TO 2 50 44 13.63636
2 TO A 25 32 -21.875
2 TO B 25 28 -10.7143
3 TO 4 291.8 296 -1.41892
3 TO C 40.2 43 -6.51163
4 TO 5 266.8 261 2.222222
4TO D 25 31 -19.3548
5 TO 6 232 219 5.936073
6 TO E 55 55.5 -0.9009
6 TO F 177 180 -1.66667
5 TO 7 34.8 32 8.75
7 TO G 17 16 6.25
7 TO H 17 16 6.25

4. RESULTS AND DISCUSSIONS

Figure 5: Diameter Vs Deviation in Mass Flow Rate

Figure 6: % Deviation of Experimental Mass Flow Rate

at Unit with Respect to Required Rate

www.tjprc.org SCOPUS Indexed Journal [email protected]

688 Neeraj Chaurasia & S. P. S. Rajput

Figure 7: Diameter Vs Deviation in Mass Flow Rate

Figure 8: % Deviation of Experimental Mass Flow Rate at

Sections with Respect to Calculated Rate

Engineering practices for pipe systems required calculation of mass flow rate and head loss. In this experimental
setup, we have calculated mass flow rate to each LRUs and complex duct geometry. The experimental mass flow rate is
compared with the theoretical and required one. It is observed that the mass flow rate is more in case of unit 1, 2, 4, and 6.
In case of unit 3 and 5, the amount of mass flow rate is less than the required one.

Figure 6 represents, the graph between % deviation of experimental mass flow rate at Unit with respect to
required mass flow rate and observed significant variations from -21% to 60%. This % deviation in mass flow rate in
various units is due to the abrupt change in diameter of ducts and the position of the unit in the complex duct geometry.
Figure 8 represents the graph between % deviation of experimental mass flow rate at sections with respect to calculated
rate and observed that the deviation are due to the abrupt change in diameter of duct. The graph between diameters of the
complex duct geometry and percentage deviation in mass flow rate has plotted and observed that the percentage deviation
is more in case of abrupt changes in pipe diameter from large diameter pipes to the smaller diameter.

5. CONCLUSIONS

This paper presents the experimental results of mass flow rate to the LRU’s and the complex duct geometry of
aircraft. From the results, it has been concluded that there is variation in mass flow rate required by the LRU’s as suggested
by original unit manufacture. The variation in mass flow rate in complex duct geometries is observed, which leads to

Impact Factor (JCC): 7.6197 SCOPUS Indexed Journal NAAS Rating: 3.11
Experimental Flow Analysis of Aerospace Duct of 689
Equipment Cooling System in an Aircraft

improper distribution of flow to the LRU’s, cause the malfunctioning of LRU’s. For the safety and functioning of LRU’s,
the assumptions taken by the designers for redistribution of flow does not hold good. This again needs to be analyzed.

REFERENCES

1. Manning, Francis S.; Thompson. Richard E. (1991). Oilfield Processing of Petroleum. Vol. I: Natural Gas. PennWell Books.
ISBN 0-87814-343-2.

2. Colebrook, C. F. and White, C. M. (1937). "Experiments with Fluid Friction in Roughened Pipes". Proceedings of the Royal
Society of London. Series A, Mathematical and Physical Sciences 161 (906): 367-381. Bibcodc: 1937RSPSA.161... 367C. doi:
I 0.1 098/rspa.1937.0 150.

3. Nikam, HARSHAL., Mishra, PREM., & Bharambe, SAYALI. (2014). Design and Analysis of Brake Rotor with Parameter
Optimization. International Journal of Automobile Engineering Research and Development (IJAuERD), 4(4), 21-30.

4. Colebrook, C. F. (February 1939). "Turbulent flow in pipes, with particular reference to the transition region between smooth
and rough pipe laws". Journal of the institution of Civil Engineers (London).

5. VDI Heat Atlas second edition page 1058 (TSBN 978-3- 540-77876-9)

6. More, A. A. (2006). "Analytical solutions for the Colebrook and White equation and for pressure drop in ideal gas flow in
pipes". Chemical Engineering Science 61 (16): 5515-5519. Doi:10.1016/j.ces.2006.04.003.

7. P. K. Swamee, A. K. Jaine, Explicit equations for pipe flow problems, J Hydraulics Div, ProcASCE (1976), pp. 657- 664 May

8. Serghides, T. K (1984). "Estimate friction factor accurately". Chemical Engineering Journal 91(5): 63-64.

9. Mitra, N. S., Doloi, B., & Bhattacharyya, B. (2014). Analysis of traveling wire electrochemical discharge machining of Hylam
based composites by Taguchi method. International Journal of Research in Engineering and Technology, 2(2), 223-236.

10. Goudar, C. T., Sonnad, J. R. (August 2008). "Comparison of the iterative approximations of the Colebrook-White equation".
Hydrocarbon Processing Fluid Flow and Rotating Equipment Special Report(August 2008): 79- 83.

11. Brkic, Dejan (2011). "An Explicit Approximation of Colebrook's equation for fluid flow friction factor". Petroleum Science
and Technology 29 (15): 1596-1602. doi: 10.1080/10916461003620453.

12. Mahrous, A., Mahmoud, S., Al-Dadah, R., & Matrawy, K. (2015). Cooling Via Human Forearm for Multiple Sclerosis.

13. Adrian Bejan, Allan D. Kraus, Heat transfer handbook, John Wiley & Sons, 2003.

www.tjprc.org SCOPUS Indexed Journal [email protected]