SciELO - Scientific Electronic Library Online

 
vol.88 número219Modificación superficial de una fundición nodular ferrítica a través de una nitruración iónicaPérdidas no técnicas en un sistema de distribución eléctrica. La eficiencia energética como proceso continuo índice de autoresíndice de materiabúsqueda de artículos
Home Pagelista alfabética de revistas  

Servicios Personalizados

Revista

Articulo

Indicadores

Links relacionados

  • En proceso de indezaciónCitado por Google
  • No hay articulos similaresSimilares en SciELO
  • En proceso de indezaciónSimilares en Google

Compartir


DYNA

versión impresa ISSN 0012-7353versión On-line ISSN 2346-2183

Dyna rev.fac.nac.minas vol.88 no.219 Medellín oct./dic. 2021  Epub 15-Mar-2022

https://doi.org/10.15446/dyna.v88n219.93809 

Articles

Numerical analysis of aircraft flow exhaust gases to design a curved deflector for Mexico City international airport

Análisis numérico del flujo de gases de escape de aviones para diseñar un deflector curvo para el aeropuerto internacional de la Ciudad de México

Jorge Luis Garrido-Télleza 
http://orcid.org/0000-0002-3126-7467

Eusebio Eduardo Hernández-Martíneza 
http://orcid.org/0000-0001-7020-7929

José Carlos Jiménez-Escalonaa 
http://orcid.org/0000-0001-9309-5245

aESIME Ticomán, Instituto Politecnico Nacional, Mexico City, Mexico. jgarrido@ipn.mx, euhernandezm@ipn.mx, jjimeneze@ipn.mx


Abstract

The numerical simulation results of the flow exhaust gases emitted by different aircraft engines are presented in this work. These results permitted the design of the most appropriate circular arc for a curved deflector to correctly direct the flow at the time planes are taking-off on the left runway number five of Mexico City International Airport (AICM) and not to affect the operability and route of the inter-terminal train, called as “Aerotren”. This train has four carriages of 25 passengers each and transports users of Mexico City International Airport from terminal T-1 to terminal T-2 and vice versa. For the numerical simulation, several engine models of different aircraft were used, including the Boeing 747-400, the Boeing 777-200LR/-300ER, the Airbus A340-600, and the McDonnell Douglas MD11. The results presented here are those from the McDonnell Douglas MD11 airplane, which is the most critical case since it has an engine in the rear part of the fuselage on the vertical empennage, and this causes the flow of exhaust gases to arrive in a direction from top to bottom on the upper part of the Aerotren, which for safety should not exceed 15.0 m/s (54 km/hr). The impact of the explosion of these engine models was used to choose the most appropriate deflector curvature angle in the take-off area. The numerical simulation was carried out through the Fluent CFD program, solving the Navier-Stokes equations, the standard (-( model, turbulence, and atmospheric conditions in Mexico City. An experimental stage and the visualization of the behavior of the flow of exhaust gases within a wind tunnel are also presented. A 1: 125 scale model was used for both the deflectors with the radius of curvature obtained and the train car, in an area of one square meter and a grid background screen.

Keywords: numerical analysis; aircraft gas turbine engine flow; computational fluid dynamics; curved deflector; Mexico City international airport; wind tunnel test; visualization

Resumen

En este trabajo se presentan los resultados obtenidos de la simulación numérica del flujo sobre un deflector curvo, con la finalidad de obtener el arco circular más apropiado para direccionar el flujo de los gases de escape en el momento de despegue de los aviones en la pista 05 izquierda del Aeropuerto Internacional de la Ciudad de México (AICM) y no afectar la operatividad y vialidad del tren inter-terminales, también llamado como Aerotrén. Este tren cuenta con cuatro vagones de 25 pasajeros cada uno y tiene la finalidad de transportar a los usuarios del aeropuerto de la terminal T-1 Internacional a la terminal T-2 Nacional y viceversa. Para la simulación numérica se utilizaron varios modelos de motores de diferentes aviones entre ellos: el Boeing 747-400, el Boeing 777-200LR/-300ER, el Airbus A340-600 y el McDonnell Douglas MD11. Los resultados que se presentan en este trabajo son del avión McDonnell Douglas MD11 por ser el caso crítico ya que cuenta con un motor en la parte posterior del fuselaje sobre el empenaje vertical y esto provoca que el flujo de los gases de escape llegue en una dirección desde arriba hacia abajo sobre la parte superior del tren, el cual por seguridad no deben sobrepasar la velocidad de 15.0 m/s (54 Km/hr). También se presentan los resultados experimentales de la visualización del comportamiento del flujo dentro de un túnel de viento, para esto se utilizó un modelo a escala de 1:125 tanto del deflector como del vagón del tren, en un área de 1 m2 y una pantalla de fondo cuadriculada.

Palabras claves: análisis numérico; flujo de motores de turbina de gas en aviones; dinámica computacional de fluidos; deflector curvo; aeropuerto internacional de la Ciudad de México; prueba en túnel de viento; visualización

1. Introduction

An aircraft in the take-off process can generate a flow in the form of a torrent, being able to present flow speeds up to 95 km/h for a distance of 300 m from the position of the engine departure [1]. These high-speed gases can have negative effects on the airport infrastructure or on other vehicles that are at an influence distance from the jet engine output [2]. Due to this effect, since the beginning of the 50's, mechanical elements have been designed that seek to divert the flow of the output of these engines in airport facilities (see [3-5]). More recently, the use of these mechanical deflectors in aircraft carriers has been studied for take-off operations that involve various problems that must be solved for the safety of the aircraft, personnel, and facilities on the take-off deck, such as high pressure and temperature of the exhaust gas flow, the influence of the wind and the correct location and angle of the mechanical elements to deviate the flow. In [6-8] the results of the analysis of the impact of the flow on a portable deflector adapted in the take-off area of aircraft carriers are presented. Based on a CFD simulation, safe conditions and locations for the baffle-mount operation are suggested, which are very useful for the minimum safety of take-off operations. In general, the main difference that can be highlighted is that the standard (-( equations and the Navier-Stokes equations in three dimensions were used in the first two works, while the Reynolds model was used in [8,9], with Navier-Stokes averaged with different turbulence models. On the other hand, Mexico City International Airport is the main airport in Mexico, and in 2019 it reported more than 458,700 landing and takeoff operations. The physical characteristics of the runways classify it as a category 4E airport according to the classification of the International Civil Aviation Organization (ICAO). The aircraft that use left runway 05 affect the area where the Aerotren rail passes, also known as the Interterminal Train, which transports passengers between terminals 1 and 2 of the AICM. The area of the runway used for the aircraft to begin their takeoff run is relatively close to the periphery of the airport and in the area where the Aerotren tracks are installed, as shown in Fig. 1.

Source: Obtained from Google Earth.

Figure 1 Location of the five-left header area of the Mexico City International Airport runways. 

1.1 Contribution

The purpose of this work is to carry out a numerical simulation of the flow of exhaust gases from different engines, for the design of a curved deflector, in order to obtain the most appropriate circular arc to direct the flow of exhaust gases in the take-off process of aircraft on the left runway number five of the AICM so as not to affect the operation and roads of the Aerotren, which is the system for transporting passengers between the two airport terminals. For the numerical simulation process, the impact of the explosion of several aircraft engine models is carried out to obtain the best angle of curvature of the deflector in the take-off zone and thus be able to avoid any inconvenience and not affect the operation of the Aerotren. In the simulation, CFD Fluent software is used, solving the Navier-Stokes equations, the standard (-( model, for turbulence, as well as the atmospheric conditions in Mexico City.

To obtain the most appropriate arc of curvature of the deflector, several analyses were carried out with different geometries and radius ranging from 3.0, 3.5, 4.0 to 6.0 m, obtaining the best results with a deflector with a curvature radius of 4.0 m and dimensions of 4.0 m height and 3.0 m base, implementing a numerical simulation in all cases with this deflector. The numerical simulation of the behavior of the flow with several engines on a circular arc deflector makes it possible to observe the behavior of the exhaust gases in such a way that they do not affect the operation of the Aerotren, which is the exclusive means of transport for connecting passengers who need to transfer between terminals T1 and T2 of the International Airport “Benito Juárez” in Mexico City. In its first stage, the Aerotren consists of four carriages that have an installed capacity to serve up to 100 passengers in a critical time; the 3 km journey in a single trip takes less than 6 minutes traveling at an average speed of 45 km/hr.

2. Materials and methods

2.1 Computational model

The numerical simulation was carried out using the commercial FLUENT® software [10]. A computer program under the computational fluid dynamics (CFD) scheme, this code is structured by the finite volumes’ method. For the generation of the geometric model, the computational model and the boundary conditions, the GAMBIT module was used, which is the preprocessor of FLUENT. The geometric model is a two-dimensional control volume that defines the space where the numerical simulation will be carried out. The computational model is the discretization of the geometric model in small squared geometry control volumes, in each element of which a finite value of a physical property is obtained, which can be speed, pressure, temperature, etc., and the boundary conditions, which are restrictions that are assigned to the control volume and which can be entrance, exit, wall, symmetry, etc. The dimensions of the geometric model and boundary conditions were defined as follows: distance from the front of the outer wing engine of the Boeing 747-400, Boeing 777-200LR/-300ER, Airbus A 340-600 aircraft, distance to the rear of the deflector location 56.0 m; distance from the rear of the deflector to the train right of way 7.80 m; railroad rights-of-way and height of train, 4.0 m by 5.0 m, respectively; deflector dimensions 4.0 m high by 3.0 m base and 4.0 m radius of curvature. For the McDonnell Douglas MD11 aircraft, the distance from the rear of the tail engine to the rear of the deflector is 48.0 m. The heights from the runway floor to the base of the aircraft engines and the values of the speeds of the engine gases flows were taken from the manuals provided by the AICM authorities [11].

Figs. 2-5 present the models of the airplanes and the parameters used for the numerical simulation. Specifically, Fig. 2 presents the model for the Boeing 747-400, Fig. 3 the Airbus A 340-600, Fig. 4 the Boeing 777- 200LR/300ER, and Fig. 5 the McDonnell Douglas MD11.

Source: Authors.

Figure 2 Schematic of the geometric modeling for the Boeing 747-400 airplane. 

Source: Authors.

Figure 3 Schematic of the geometric modeling for the Airbus A 340-600 aircraft. 

Source: Authors.

Figure 4 Schematic of the geometric modeling for the Boeing 777-200LR/300ER airplane. 

Source: Authors.

Figure 5 Schematic of the geometric modeling for the McDonnell Douglas MD11 aircraft. 

2.2 Numerical formulation

The steady-state 2D numerical simulation of each of the aforementioned geometries was performed using the commercial software Fluent. This program solves the Navier-Stokes equations on a structured mesh, for incompressible flow with an implicit formulation. Turbulence is simulated with the (-( standard model (two equations) [12]. The pressure-velocity coupling is calculated through the simple algorithm, with a second order discretization for the momentum and the turbulent kinetic energy [13]. Traditionally, CFD theory has been used to numerically simulate the high temperature and high outflow speeds of aircraft [7], both in aircraft carriers and in ground facilities. Other approaches, such as empirical approximation formulations, the wind tunnel test method, and actual measurement methods, have also been applied for this same purpose. The most precise methods are the last two approaches [6]; however, they cannot be addressed in the initial stages of the analysis.

2.3The K-ε standard model

The K-ε standard model is a semi-empirical model based on the transport equations for turbulent kinetic energy (K) and for turbulent kinetic energy dissipation (ε), and was adopted in this study for simulation of the turbulence. When the model is derived, the flow is considered to be totally turbulent and the effects of molecular viscosity are negligible. Therefore, the K-ε model is valid only for turbulent flows [14].

Turbulent kinetic energy (K):

The dissipation of turbulent kinetic energy (ε) is considered as:

From Eq. (1,2), P k represents the generation of turbulent kinetic energy due to the velocity gradient:

and G k is the turbulent kinetic energy generation due to buoyancy forces:

Y M represents the contribution of dilation fluctuations in compressible turbulence due to the dissipation rate and is determined by:

M t is the turbulent Mach number

where a is defined as

2.4 Model constants

The model constants have the following values: Cμ = 0.09, C = 1.44, C = 1.92 and the Prandtl numbers for the equations of ( and ( are (k = 1.0 and σε = 1.3, respectively. The degree to which ε is affected by buoyancy is determined by the constant C and is calculated according to the following relationship:

3. Analysis of study cases

3.1 Conditions for numerical simulation

Figs. 6-9 show the geometric models of each aircraft used for this analysis, representing the dimensions of the control volume, the distances from the aircraft engine to the front base of the deflector location, as well as the railroad rights-of-way, and the height and width of the air-train. Figs. 10-13 present the computational models of each of the engines of these aircrafts with the number of squared cells for each case.

Source: Authors.

Figure 6 Diagram, conditions and control volume for the Boeing 747-400. 

Source: Authors.

Figure 7 Diagram, conditions and control volume for the Boeing 777-200ER. 

Source: Authors.

Figure 8 Diagram, conditions and control volume for the Airbus A 340-600. 

Source: Authors.

Figure 9 Diagram, conditions and control volume for the McDonnell Douglas MD11. 

Source: Authors.

Figure 10 Computational model and conditions for the evaluation of the Boeing 747-400 airplane with 364,940 cells. 

Source: Authors.

Figure 11 Computational model and conditions for the evaluation of the Boeing 777-200LR airplane with 324,798 cells. 

Source: Authors.

Figure 12 Computational model and conditions for the evaluation of the Airbus A 340-600 aircraft with 331,904 cells. 

Source: Authors.

Figure 13 Computational model and conditions for the evaluation of the McDonnell Douglas MD11 airplane with 338,530 cells. 

4. Results and discussion

As mentioned at the beginning, the results presented are only from the McDonnell Douglas MD11 aircraft, which is the critical case. The results of the numerical evaluation are shown in Figs. 14-16. In particular, in Fig. 14 it is possible to observe in the form of contours, the speed of the flow of the exhaust gases from the aircraft engine, over the deflector and the train. The exhaust gases that are emitted by the aircraft engine from top to bottom impact the top of the deflector and are deflected onto the upper surface of the Aerotren. In Fig. 15 one can see a close-up only of the deflector and the train, and in Fig. 16 a velocity field is presented, where it can be seen that the velocity of the flow of the exhaust gases over the upper part of the deflector is deviated towards the upper part of the train and this allows a damping of the flow speed; we can also observe that the speed reached in the upper left part of the train is of the order of 15 m/s, a condition that allows satisfactory operation of the Aerotren.

Source: Authors.

Figure 14 Results for the critical case study: speed magnitude contours for the aircraft engine (m/s), deflector and Aerotren. 

Source: Authors.

Figure 15 Results for the critical case study: speed magnitude contours (m/s). 

Source: Authors.

Figure 16 Results for the critical case study: velocity vector field (m/s). 

4.1 Design of the curved deflector

From the numerical results of the critical case study, the final considerations for the design proposal of the curved deflector were obtained. In Fig. 17 the graph of the behavior of the static pressure with respect to the length of the deflector chord is presented. On the other hand, Fig. 18 shows a diagram with the values of the force components with respect to the X and Y axes, as well as the respective moment, with a magnitude in the clockwise direction of -5,615.71 mN.

The Fluent software has the option to integrate the pressure distribution and extract moments at any point.

Source: Authors.

Figure 17 Design proposal for the curved deflector: diagram of static pressure versus length of the curve. 

Source: Authors.

Figure 18 Design proposal for the curved deflector: schematic of the force components with respect to the X and Y axes, as well as the total moment. 

4.2 Visualization

In addition to the numerical evaluation and based on the proposed design for the curved deflector, a scaled physical model was built for the purpose of its experimental evaluation within a wind tunnel. Fig. 19 presents the model of the deflector and the Aerotren at 1:125 scale, which were used for the visualization of the flow. Fig. 20 presents the model with a background gridded screen representing the space inside the wind tunnel. It can be seen how the air flow is deviated by the deflector towards the upper part of the train, allowing the air speed not to impact directly on the upper left side of the train, which was one of the critical points for this analysis.

Source: Authors.

Figure 19 Wind tunnel tests for the proposed deflector: model of the deflector and the Aerotren at 1: 125 scale. 

Source: Authors.

Figure 20 Wind tunnel tests for the proposed deflector: visualization of the flow on the proposed model. 

Based on the results of this work, the calculated curved deflector was built on the left runway number five of the AICM. In particular, Fig. 21 shows the curved deflector that was built as a result of this analysis for the head of that runway, while Fig. 22 shows the Aerotren successfully operating and transporting passengers from terminal T1 to T2 and vice versa.

Source: Captured by authors on site.

Figure 21 Curved deflector calculated and installed on five-left runway of the AICM. 

Source: Captured by authors on site.

Figure 22 Image of the operation of the Aerotren transporting passengers from T1 to T2 and vice versa. 

5. Conclusions

Based on a CFD analysis, this article presents the numerical simulation results of the flow exhaust gases emitted by different aircraft engines, intended to design a curved deflector to correctly direct the flow at the time planes are taking-off on the left runway number five of AICM and not to affect the operability and route of the inter-terminal train. The curved deflector with 4.0 m height, 3.0 base, and 4.0 m radius of curvature presented the best performance according to the numerical analysis, in this way the curved deflector was built with these dimensions. It was identified that McDonnell Douglas MD11 airplane is the critical case since it has an engine in the rear part of the fuselage on the vertical empennage, and this causes the flow of exhaust gases to arrive in a direction from top to bottom on the upper part of the Aerotren. However, it was observed the upper-left speed on the train does not exceed the limit of safety speed defined by authorities. In addition, visualization results of the flow behavior of exhaust gases in a wind tunnel permitted to observe that curved deflector deviate in a proper way the flow on the upper part of Aerotren, without affecting its operability, and transporting users from terminal T1 to T2 and vice versa.

References

[1] Nadal-Mora, V., Pezzotti, S., Di Barnardi, A., Di Gregorio, P. and Pesarini, A., Theoretical analysis and experimental validation of the interaction between the exhaust flow of aircraft engines in take-off operations, an anti-jet barrier and a vehicular route, Proceedings of the First Congress of the Ibero-American Network for Research in Air Transport (RIDITA), 2007, pp. 339-358. [ Links ]

[2] Morrison, R., Ground jet blast hazard, ASRS Directline, Issue 6, August 1993. [ Links ]

[3] Hayden, H.J., Jet engine exhaust deflector. U.S. Patent 2,826,382 issued March 11, 1958. [ Links ]

[4] Stanley, B.L. and Dunes, P., Split exhaust jet blast deflector fence. U.S. Patent 5,429,324 issued July 4, 1995. [ Links ]

[5] Campion, G.P. Blast deflector. U.S. Patent 6,802,477 issued October 12, 2004. [ Links ]

[6] Kuishi, Y., Liangliang, C., Hu, L. and Yunliang, W., Analysis of jet blast impact of embarked aircraft on deck takeoff zone. Aerospace Science and Technology, 45, pp. 60-66, 2015. DOI: 10.1016/j.ast.2015.04.010 [ Links ]

[7] Kuishi, Y., Yicheng, S., Hu, L. and Weigang, G., Analysis of the flow of carrier-based aircraft exhaust jets impact on the flight deck. International Journal of Aeronautical and Space Sciences, 16(1), pp. 1-7, 2015. DOI: 10.5139/IJASS.2015.16.1.1 [ Links ]

[8] Fu-Dong, G., De-Xin, W., Hai-Dong, W. and Ming-Ming, J., Numerical analysis and verification of the gas jet from aircraft engines impacting a jet blast deflector. Chinese Journal of Mechanical Engineering, 31(86), pp. 1-11, 2018. DOI: 10.1186/s10033-018-0285-7 [ Links ]

[9] Song, M., Jianguo, T., Xiankai, L. and Hao, J., The effect analysis of an engine jet on an aircraft blast deflector. Transactions of the Institute of Measurement and Control Journal of Mechanical Engineering, 41(4), pp. 990-1001, 2018. DOI: 10.1177/0142331218755892 [ Links ]

[10] Manual de usuario de Fluent, s.a. [ Links ]

[11] Airplane characteristics for airport planning manuals: Boeing 747-400, Boeing 777-200LR/-300ER, Airbus A340-600, McDonnell Douglas MD11, s.a. [ Links ]

[12] Launder, B.E. and Spalding, D.B., Lectures in mathematical models of turbulence. Academic Press, London, U.K., 1972. [ Links ]

[13] Sarkar, S. and Balakrishnan, L., Application of a Reynolds-stress turbulence model to the compressible shear layer, ICASE Report 90-18, NASA CR 182002, 1990. [ Links ]

[14] Baxter, L.L., Turbulent transport of particles, PhD Thesis, Brigham Young University, Provo, Utah, USA,1989. [ Links ]

J.L. Garrido-Téllez, is a BSc. Eng in Mechanical Engineering, in 1982, MSc. in Mechanical Engineering in 2002 and Dr. of Science in Mechanical Engineering in 2018, Diplomat in Computer Engineering in 1996, Diplomat in Thermal Engineering in 2002, all of them form the National Polytechnic Institute, Mexico City. He is currently professor-researcher at ESIME Ticoman of the National Polytechnic Institute, in Mexico City, Mexico. ORCID: 0000-0002-3126-7467

E.E. Hernández-Martínez, is a BSc. Eng in Electronics Engineer from the Technological Institute of Celaya, Mexico. In 2004, he received the MSc. in Instrumentation and Automatic Control from the Autonomous University of Queretaro, Mexico. In 2010, he obtained a Dr. in Mechanical Engineering from the University of Cassino, Italy and the Autonomous University of Queretaro, Mexico. He is with ESIME Ticoman of the National Polytechnic Institute, in Mexico City. His research interests include the simulation and mechatronic design of systems, with a special interest in dynamic and kinematic modeling, optimized design and error modeling in parallel robots. He has authored or co-authored approximately 30 research papers published in international journals and / or extensive conference publications. ORCID: 0000-0001-7020-7929

J.C. Jiménez-Escalona, is a BSc. Eng in Aeronautical Engineer in 1988 from the National Polytechnic Institute, Mexico City. MSc. in Physics of the Atmosphere in 2000 and Dr. of Science in Physics of the Atmosphere in 2005, all of them form the National Autonomous University of Mexico. He is currently professor-researcher at ESIME Ticoman of the National Polytechnic Institute, in Mexico City. His research interests include risk mitigation in aviation due to natural phenomena in the atmosphere. ORCID: 0000-0001-9309-5245

How to cite: Garrido-Téllez, J.L., Hernández-Martínez, E.E. and Jiménez-Escalona, J.C., Numerical analysis of aircraft flow exhaust gases to design a curved deflector for Mexico City international airport.. DYNA, 88(219), pp. 210-217, October - December, 2021.

Received: February 22, 2021; Received: October 21, 2021; Accepted: November 10, 2021

Creative Commons License The author; licensee Universidad Nacional de Colombia.