Numerical Analysis of Airflow in Trachea Affected by Thyroid Cancer Using Finite Volume Method

Arif Fatahillah, Miftahul Jannah, Susi Setiawani, Toto Bara Setiawan, Arika Indah Kristiana


The study aims to analyze the airflow in the trachea due to thyroid cancer. Thyroid cancer is a cancer that attacks the thyroid gland. The thyroid gland is located in the neck below the Adam's apple. In general, thyroid cancer causes stenosis of the trachea due to cancer compression of the tracheal wall. This problem is solved using the finite volume method with SIMPLE (Semi Implicit Method for Pressure Linked Equation) discretization. This method is well suited for the discretization of irregularly shaped fluid flows. The software used is Matlab and Ansys Fluent. Matlab is used to compute numerical solutions and display graphs, while Fluent is used to visualize the air flow. The result of this study is that the greater the percentage of stenosis, the velocity of the airflow is faster and the pressure of the airflow is lower. The most critical condition for respiratory power loss is when the airflow velocity exceeds the normal limit, i.e. when the narrowing is more than 50%, the initial velocity is more than 0,1 m/s.


Thyroid cancer, Tracheal stenosis, Finite Volume Method, SIMPLE

Full Text:



W. H. Man, W. A. A. de S. Piters, and D. Bogaert, “The Microbiota of the Respiratory Tract: Gatekeeper to Respiratory Health,” Nature Reviews Microbiology, vol. 15, no. 5. Nature Publishing Group, pp. 259–270, Apr. 10, 2017. doi: 10.1038/nrmicro.2017.14.

B. E. M. Brand-Saberi and T. Schäfer, “Trachea: Anatomy and physiology,” Thoracic Surgery Clinics, vol. 24, no. 1. pp. 1–5, Feb. 2014. doi: 10.1016/j.thorsurg.2013.09.004.

P. W. Furlow and D. J. Mathisen, “Surgical anatomy of the trachea,” Ann Cardiothorac Surg, vol. 7, no. 2, pp. 255–260, Mar. 2018, doi: 10.21037/acs.2018.03.01.

S. Haykal, M. Salna, T. K. Waddell, and S. O. Hofer, “Advances in Tracheal Reconstruction,” Plastic and Reconstructive Surgery, vol. 2, no. 7. Lippincott Williams and Wilkins, Jul. 09, 2014. doi: 10.1097/GOX.0000000000000097.

Y. Khan and A. Farhan, Histology, Thyroid Gland. Treasure Island (FL): Stat Pearls, 2019.

K. Arrumugam, N. P. Ekawati, and I. M. Gotra, “Characteristic of Thyroid Carcinoma based on Age, Sex and Histology Type at Sanglah General Hospital, Bali, Indonesia, 2016,” DiscoverSys, vol. 9, no. 3, pp. 124–126, 2018, doi: 10.1556/ism.v9i3.285.

“Cancer Facts & Figures 2015,” Atlanta, 2015.

M. E. Cabanillas, D. G. McFadden, and C. Durante, “Thyroid cancer,” The Lancet, vol. 388, no. 10061. Lancet Publishing Group, pp. 2783–2795, May 2016. doi: 10.1016/S0140-6736(16)30172-6.

W. J. Poynot, K. A. Gonthier, M. E. Dunham, and T. W. Crosby, “Classification of Tracheal Stenosis in Children based on Computational Aerodynamics,” J Biomech, vol. 104, pp. 1–11, May 2020, doi: 10.1016/j.jbiomech.2020.109752.

Q. T. Nguyen et al., “Diagnosis and Treatment of Patients with Thyroid Cancer,” vol. 8, no. 1, pp. 30–40, 2015, [Online]. Available:

G. Pellegriti, F. Frasca, C. Regalbuto, S. Squatrito, and R. Vigneri, “Worldwide increasing incidence of thyroid cancer: Update on epidemiology and risk factors,” Journal of Cancer Epidemiology, vol. 2013. 2013. doi: 10.1155/2013/965212.

T. Zobaer and A. Sutradhar, “Modeling the effect of tumor compression on airflow dynamics in trachea using contact simulation and CFD analysis,” Comput Biol Med, vol. 135, Aug. 2021, doi: 10.1016/j.compbiomed.2021.104574.

A. Issakhov, A. Mardieyeva, Y. Zhandaulet, and A. Abylkassymova, “Numerical study of air flow in the human respiratory system with rhinitis,” Case Studies in Thermal Engineering, vol. 26, Aug. 2021, doi: 10.1016/j.csite.2021.101079.

H. Shamohammadi, S. Mehrabi, S. Sadrizadeh, M. Yaghoubi, and O. Abouali, “3D numerical simulation of hot airflow in the human nasal cavity and trachea,” Comput Biol Med, vol. 147, Aug. 2022, doi: 10.1016/j.compbiomed.2022.105702.

F. Moukalled, L. Mangani, and M. Darwish, “Fluid Mechanics and Its Applications The Finite Volume Method in Computational Fluid Dynamics.” [Online]. Available:

H. Khawaja and M. Moatamedi, “Semi-Implicit Method for Pressure-Linked Equations (SIMPLE) - solution in MATLAB®,” 2018.

M. M. Bhatti, M. Marin, A. Zeeshan, and S. I. Abdelsalam, “Editorial: Recent Trends in Computational Fluid Dynamics, 2nd Edition,” Mathematical and Statistical Physics, 2020, doi: 10.3389/978-2-88966-900-4.

Y. Zhao and X. Su, Computational Fluid-Structure Interaction: Methods, Models, and Applications. Chennai: Matthew Deans, 2019. doi:

P. D. Dwidmuthe, G. G. Dastane, C. S. Mathpati, and J. B. Joshi, “Study of blood flow in stenosed artery model using computational fluid dynamics and response surface methodology,” Can J Chem Eng, vol. 99, no. S1, pp. S820–S837, Oct. 2020, doi: 10.1002/cjce.23991.

J. H. Kao and P. Y. Tseng, “Application of Computational Fluid Dynamics (CFD) Simulation in a Vertical Axis Wind Turbine (VAWT) System,” in IOP Conference Series: Earth and Environmental Science, Institute of Physics Publishing, Feb. 2018. doi: 10.1088/1755-1315/114/1/012002.

J. M. Durán, “What is a Simulation Model?,” Minds Mach (Dordr), vol. 30, no. 3, pp. 301–323, Sep. 2020, doi: 10.1007/s11023-020-09520-z.

J. Wang, Y. Wang, and J. Shi, “On Efficiency and Effectiveness of Finite Volume Method for Thermal Analysis of Selective Laser Melting,” Engineering Computations (Swansea, Wales), vol. 37, no. 6, pp. 2155–2175, May 2020, doi: 10.1108/EC-03-2019-0106.

Y. Zhao, Y. Kim, M. S. Hossain, M. Nazem, J. Liu, and Y. Hu, “Numerical advancements on the analysis of dynamically installed anchors,” Ocean Engineering, vol. 182. Elsevier Ltd, pp. 343–359, Jun. 15, 2019. doi: 10.1016/j.oceaneng.2019.04.079.

H. K. Versteeg and W. Malalasekera, An Introduction to Computational Fluid Dynamics Second Edition, 2nd ed. 2007. [Online]. Available:



  • There are currently no refbacks.

Copyright (c) 2024 Arif Fatahillah

Creative Commons License
This work is licensed under a Creative Commons Attribution-ShareAlike 4.0 International License.

Editorial Office
Mathematics Department,
Universitas Islam Negeri Maulana Malik Ibrahim Malang
Gajayana Street 50 Malang, East Java, Indonesia 65144
Faximile (+62) 341 558933

Creative Commons License
CAUCHY: Jurnal Matematika Murni dan Aplikasi is licensed under a Creative Commons Attribution-ShareAlike 4.0 International License.