- 1Department of Orthodontics, School and Hospital of Stomatology, Cheeloo College of Medicine, Shandong University & Shandong Key Laboratory of Oral Tissue Regeneration & Shandong Engineering Laboratory for Dental Materials and Oral Tissue Regeneration, Jinan, China
- 2Institute of Thermodynamics and Fluid Mechanics, School of Energy and Power Engineering, Shandong University, Jinan, China
Nasal septal deviation (NSD) is one of the most common nasal diseases. Different from common clinical examination methods, computational fluid dynamics (CFD) can provide visual flow information of the nasal cavity. The dimension and volume of the nasal cavity are easily affected by rapid maxillary expansion (RME). The purpose of this study was to use CFD to evaluate the effect of RME on the aerodynamics of the nasal cavity in children with maxillary transverse deficiency and NSD. Computational fluid dynamics was implemented after 3D reconstruction based on the CBCT of 15 children who have completed RME treatment. After treatment, the volume increases in the nasal cavity, nasopharynx, oropharynx, and pharynx were not statistically significant. The wall shear stress of the nasal cavity after RME, 1.749 ± 0.673 Pa, was significantly lower than that before RME, 2.684 ± 0.919 Pa. Meanwhile, the maximal negative pressure in the pharyngeal airway during inspiration was smaller after RME (−31.058 Pa) than before (−48.204 Pa). This study suggests that RME has a beneficial effect on nasal ventilation. The nasal airflow became more symmetrical in the bilateral nasal cavity after RME. Pharyngeal resistance decreased with the reduction in nasal resistance and the increase in the volume of oropharynx after RME.
Introduction
Nasal septal deviation (NSD) (1) is one of the most frequently encountered diseases in the rhinology clinic, with a 39.9% incidence in children (2). It can result in imbalanced airflow in the bilateral nasal cavity, hypertrophy of the turbinate, and increased nasal resistance (3–5). Increased nasal resistance may lead to a chronic mouth-breathing pattern, causing dentofacial deformities, such as bilateral maxillary crossbite and maxillary transverse deficiency, during the growth and development of the craniofacial complex (6). Furthermore, increased nasal resistance may be responsible for obstructive sleep apnea syndrome (OSAS) because the pharyngeal airway tends to collapse with greater negative pressure in the pharyngeal (7). Therefore, it is necessary to evaluate the changes in nasal airflow dynamics after rapid maxillary expansion (RME) in children with maxillary transverse deficiency and NSD.
Rapid maxillary expansion is a well-known orthodontic procedure used in growing patients to correct maxillary transverse deficiency by opening the midpalatal suture and moving maxilla laterally (8, 9). It has been demonstrated that RME can increase nasal volume, reduce nasal resistance, and improve nasal airway ventilation by many studies using two-dimensional (2D) cephalometric measurement (10), CBCT (11), rhinomanometry (12), PSG (13), acoustic rhinometry (14), and computational fluid dynamics (CFD) (15, 16). Some studies have reported that RME increases the length of the septum and corrects NSD during childhood (17, 18). These RME-related studies were designed to explore morphological changes of the nasal structure and the measurement of nasal resistance. At present, studies on NSD in children are mainly divided into three categories: One includes comparative studies on the effect of NSD on craniofacial growth in children (19, 20); another focuses on the effect of septoplasty on nasal function (21–23); and the other includes studies that discuss the effect of maxillary expansion on nasal anatomy (17). Considering that NSD may cause the patient to develop mouth-breathing patterns leading to high vault and maxillary transverse deficiency, RME may be an effective treatment (17). However, there has been no reports related to changes of nasal airflow dynamics after RME in children with NSD and maxillary transverse deficiency. Therefore, a detailed investigation of nasal airflow is beneficial for us to further understand the relationship between nasal structure and nasal aerodynamics, and it could provide an important theoretical basis for early orthopedic treatment.
Currently, CFD has been recognized as an appropriate tool for simulating the dynamics of airflow, as it enables us to visualize three-dimensional (3D) airflow characteristics (24–28). The objective of this study was to evaluate nasal aerodynamics in children with maxillary transverse deficiency and NSD after RME using CFD.
Materials and Methods
Subjects
The study was reviewed and approved by the Shandong University School of Stomatology Research Ethic Board (protocol number 20200802). All written informed consents were received from the patients. The sample size was calculated based on an α of 0.05 and a β of 0.2 to detect the difference of 242.66 Pa in maximum pressure between groups, with a 300.22-Pa estimated standard deviation (29). The power analysis indicated that a sample size of 13 was required. The sample of this study consisted of 15 (6 boys and 9 girls with a median age of 9.57 ± 1.51 years) children with maxillary transverse deficiency (30) and NSD (31) who were treated with RME. The degree of NSD depends on the angle of the line between the crista galli and the premaxilla and the line between the crista galli and the most prominent point of the septal bone or cartilage. Septal body asymmetry is defined as the width measured from the lateral aspect of the septal body to the septal bone or cartilage. The overall degree of septal deviation and septal body asymmetry were 4.56 ± 1.33° and 2.16 ± 0.32 mm, respectively (see Supplementary Table 1 for other anatomical measurement results). The exclusion criteria included (1) history of orthodontic or orthopedic treatment, (2) systemic diseases or any other severe craniofacial anomalies except maxillary transverse deficiency and nasal septum deviation, (3) temporomandibular joint disorders, and (4) no history of tonsillectomy. The expander was a Hyrax device that was bonded to the first deciduous molars or first premolars and the first permanent molars. Parents were instructed to turn the expansion screw 0.5 mm per day until the palatal cusps of the maxillary first molars contacted the buccal cusps of the mandibular first molars after 14 ± 2 days. The opening of the palatal suture was 4.13 ± 0.41 mm.
CBCT Data Acquisition
CBCT scans were taken (which was not routine) with the patient awake in the supine position and the Frankfort horizontal plane vertical to the floor by the same operator using the same CBCT scanner (New Tom 5G, Verona, Italy; scan time: 10 s, slice thickness: 0.3 mm, 110 kV, 5 mA). During scanning, patients were guided to close their mouths with the maximum intercuspation and hold their breath after the end of expiration.
CBCT images were obtained before and after treatment as T1 data and T2 data, respectively, and were stored in the DICOM (Digital Imaging and Communications in Medicine) file format.
Three-Dimensional Virtual Model Reconstruction
Three-dimensional models of the upper airway were constructed from the CBCT scans using the MIMICS 21.0 (Materialism's Interactive Medical Image Control System) software. The upper airway was isolated from other neighboring structures by setting the threshold of gray level between −1,024 Hounsfield Units (HU) and −220 HU. The sinuses were removed from the nasal cavity to simplify the 3D upper airway models in this study, according to Xiong et al.'s research (32). After reconstructing 2D segmented sets, primary smoothing of the models was performed to reduce computational cost and improve computational efficiency without affecting the main flow pattern inside (4). All 3D models were exported as stereolithography (STL) files (Figure 1).
Figure 1. (A) A CBCT slice of the upper airway. (B) 3D model of the upper airway. (C) Mesh generation of the upper airway 3D geometry.
Mesh Generation
Firstly, five inflation layers were set at the wall to generate a finer grid near the wall. Then, the unstructured tetrahedral volumetric grids of the upper airway flow domain were generated (Figure 1). To compromise between the computational cost and accuracy, grid independence tests were performed by repeatedly solving the pressure drops of the nasopharynx, oropharynx, and hypopharynx planes (Figure 2) with four different element size meshes to establish grid independence solutions. The relative difference (RD) is defined as follows:
According to the comparison result, 530,000 grids were selected for subsequent calculations (Supplementary Table 2), resulting in a mesh of the pre- or post-RME model with at least 2.66 million elements. The maximum Skewness of the grids was 0.8. The 3D models were then loaded into the Reynolds Average Navier-Stokes CFD solver FLUENT v15.0 (ANSYS, Inc., Canonsburg, USA) for airflow simulation.
Figure 2. The pressure and velocity profile of the airflow before and after RME. (A) before RME, (B) after RME.
Computational Fluid Dynamics
The inspiration flow in the upper airway was simulated using the Reynolds Average Navier-Stokes (RANS) CFD solver FLUENT v15.0, since the negative pressure resulting from inspiration is regarded as one of the causes of airway collapse in OSAS (33, 34). The air in the calculation domain was supposed to be a Newtonian, homogeneous, and incompressible fluid. The RANS equations were solved to account for the possible existence of turbulence with the standard k-ω model using the SIMPLE (semi-implicit method for pressure linked equation) method (35). For the inspiration process, a pressure condition (gauge pressure was 0 pa) was applied to the both nostrils (the flow outlet), and a negative velocity (V = Q/A, where Q was the volume flow rate of the model, and A was the area of the bottom of the hypopharynx) was applied to the bottom of the hypopharynx (the flow inlet). In this study, the constant inspired flow rate was set as 200 ml/s (16). For simplification, the factors of temperature, humidity and vibrissae were neglected in the simulation, while the gravity factor was considered. The direction of gravity is along the Y-axis as shown in Figure 1C. The calculation residual was set as 10−6. A second-order upwind scheme was used for discretizing the momentum, turbulent kinetic energy, and turbulent dissipation rate equations.
Validation of CFD Model
In this paper, the wall static pressure results of three different CFD numerical models (Sp-Almaras, standard k-ε, standard k-ω, k-ω SST, and LES) were compared using an in vitro experimental system introduced in a previous study (36). Three pressure taps were set on the surface of the 1:1 scaled mechanical upper airway model, which were located in the nasopharynx, oropharynx, and hypopharynx, respectively (Supplementary Figure 1). The inspiratory experiment was carried out at a constant average flow rate of 650 ml/s. Computational fluid dynamics simulations were performed under the same airway geometry and boundary conditions. Regarding the wall pressure, the standard k-ω model yielded better agreement with the experimental data than the other CFD models (Supplementary Figure 2), with a maximum difference of <17%. The reason for this slightly larger RD may be due to the limitation of STL technology, which caused the small resin burrs inside the model to affect the air flow.
CFD Outcome
Three planes were examined: (1) the choanae plane, a plane passing through the posterior nasal spine and the midpoint of the posterior upper edge of the nasal septum; (2) the palatal plane, a plane parallel to the hard palate connecting the anterior nasal spine and the posterior nasal spine; and (3) the superior border of the epiglottis plane, a plane parallel to Frankfort horizontal plane passing through the superior border of the epiglottis. The pharynx was further divided into three parts: nasopharynx, oropharynx, and hypopharynx. The volumes (Vol) of each section were measured.
After the simulation, all of the cross-sections from the nostril to choanae were selected along the Y-axis every 10 mm. Nasal resistance (R = ΔP/Q, the ratio between pressure changes from external nares to the choanae and total flow rate) (37), maximum velocity of the nasal cavity, and pressure drops (ΔP = Pmax-Pmin) were calculated from the outcome of the CFD analysis. Wall shear stress (wss), maximum negative pressure (Pmin), and maximum velocity (Vmax) of the pharynx were measured. Because the flow rate was constant in our study, ΔP reflected the changes in the resistance for each pharyngeal part.
Statistical Analysis
Statistical analysis was performed by the SPSS (version 20.0, IBM, New York, USA) software package. Intra-examiner reliability was determined by performing the measurements for each CBCT image on two separate occasions by one examiner at a 2-week interval. The intraclass correlation coefficients were calculated, and then the mean of the 2 measurements was used in statistical analysis. The error of the method was calculated using the Dahlberg formula: ME= (where d indicates deviations between the two measurements, and n indicates the number of paired objects). All measured variables were described using the mean and the standard deviation. For normal distribution data, a paired t-test was used to compare the difference between T1 and T2. In the case of abnormal distributions, Wilcoxon signed rank test was used for comparison. P < 0.05 indicated that the difference was significant.
Results
post-hoc analysis showed that the average power (1-β) of all calculation results reached 0.84. The intraclass correlation coefficients for all measurements ranged from 0.94 to 0.97, indicating sufficient reliability. The ME ranged from 0.052 to 0.081 cm3 for volume measurements and from 0.42 to 0.95 Pa for maximal negative pressure. According to all repeated analyses, the measurement errors were negligible. Taking a 9-year-old girl as an example, Figures 1–4 illustrate the results of RME.
Volume of the Upper Airway
Table 1 shows that after treatment, the average increments of the nasal cavity, nasopharynx, oropharynx, and hypopharynx between T1 and T2 were 0.81, 0.41, 0.70, and 0.35 cm3, respectively (P < 0.05), indicating that RME did not significantly change the volume of the upper airway.
Velocity of the Nasal Cavity
Figure 3 shows the velocity magnitude contours at the coronal cross-sections in the subject. Before RME, the nasal airflow was obviously asymmetrical between the two sides of the NSD nose. Before RME, the major airflow passed through the middle nasal airway (between the inferior and middle turbinates around the septum) on the convex side, while it passed through the superior part of the nasal airway (between the middle and superior turbinates around the septum) in the concave side. After RME, the main airflow path was similar to that before the treatment, but the flow of the nasal cavity on both sides becomes relatively symmetrical. The maximal airflow velocity after RME, 4.754 m/s, was lower than that before RME, 5.990 m/s.
Figure 3. The velocity magnitude contours at five coronal sections of a typical subject with nasal septal deviation. (A) before RME, (B) after RME.
Nasal Resistance
Nasal resistance after RME, 0.043 ± 0.017 Pa/(ml/s), was significantly smaller than that before RME, 0.073 ± 0.040 Pa/(ml/s) (Table 2; Figure 2). Before RME, the pressure drop was more drastic from the nostril to the choanae (Figure 2A). After RME, the pressure decreased smoothly along the airway from the nostril to the choanae (Figure 2B).
Wall Shear Stress of the Nasal Cavity
Figure 4 shows the side and bottom views of wall shear stress distribution of the nasal cavity. The maximal wall shear stress of the nasal cavity decreased from 2.684 ± 0.919 to 1.749 ± 0.673 Pa after RME. Before RME, a relatively larger shear stress was found in the middle and inferior region of the common meatus in the right nasal cavity (Figure 4A). After RME, the wall shear stress of the nasal cavity decreased significantly, and a relatively larger shear stress was found in the middle region of the common meatus in both sides of the nasal cavity (Figure 4B).
Figure 4. Side and bottom views of the wall shear stress (Pa) of nasal cavity. (A) before RME, (B) after RME.
Pressure in the Pharynx
The minimal negative pressure in the pharynx and nasopharynx drop after RME was significantly smaller than that before RME, whereas there were no significant changes in the oropharynx and hypopharynx pressure drop (Table 2; Figure 2).
Discussion
We investigated the changes of nasal airflow in children with NSD and maxillary transverse deficiency. After RME, nasal airflow became more symmetric and was observed in the middle and inferior region of the common meatus, which is beneficial for heating the inspire airflow, as the inferior turbinates are the main region of heat exchange from mucosa (38). The velocity of the nasal cavity was decreased, showing that the inspired airflow had enough time to contact the nasal mucosa to facilitate heating/cooling and humidification of inspired air.
The wall shear stresses was the vital factor in controlling the defense mechanisms of the upper airway (39). Excessive wall shear stresses related to a greater velocity of airflow along the wall of the nasal cavity may cause irritation of the blood vessels (40) and adversely affect the functioning of epithelial cells (41). The greater wall shear stress in the nasal cavity was decreased after RME, which may play a protective role in the mucous membrane of the nasal cavity.
Before RME, the pressure of pharyngeal decreased sharply, and the pressure drop increased due to increased nasal resistance. After treatment, the maximal negative pressure and the resistance of the pharyngeal during inspiration decreased with decreased nasal resistance and increased oropharyngeal volume after RME. Likewise, it was demonstrated in a study by Iwasaki et al. (37) and Hur et al. (42). The decreased resistance of the pharyngeal (pressure drop) after RME may cause a decrease in pharyngeal compliance and an increase in the capability to prevent pharyngeal collapse. This mechanism may relieve the symptoms of OSAS in children.
This research noted volume increases in the nasal cavity (0.81 cm3) and oropharynx (0.70 cm3). The increase in the oropharyngeal volume may be because the tongue is moved forward with the increase in the oral volume after RME (43). The measurement of the volume of the nasopharynx (0.41 cm3) showed a smaller increase than that of the nasal cavity, as circummaxillary structures, such as the zygomatic bone and sphenoid bone, are more resistant to RME (44). The change in the volume of the hypopharynx (0.35 cm3) may be due to the young growth of the mandible, which can promote the sagittal growth of the mandible and increase the volume of the airway (45, 46). These findings suggested that nasal and pharyngeal airflow characteristics responded positively to RME. Rapid maxillary expansion may be an effective option for children with NSD and maxillary transverse deficiency.
There are still some limitations to this study. Firstly, the sample size was relatively small. Secondly, the interaction between the tissues around the upper airway and the airflow was not taken into account. Thirdly, there was a lack of comparison between traditional RME and bone-borne RME. The latter is related to greater skeletal expansion (47), which may have a greater impact on airflow in the nasal cavity. In the future, we will conduct studies to examine the effects of RME and bone-borne RME on nasal aerodynamics in patients with different locations of septal deviation. We will also investigate different extents of inferior turbinate hypertrophy and further study the fluid–solid interaction model of the upper airway to increase the validity of the results.
Data Availability Statement
The original contributions presented in the article/Supplementary Material, further inquiries can be directed to the corresponding author.
Ethics Statement
The studies involving human participants were reviewed and approved by School of Stomatology Shandong University Research Ethic Board. Written informed consent to participate in this study was provided by the participants' legal guardian/next of kin.
Author Contributions
SC contributed to conception, design, data acquisition, analysis, interpretation, and drafted and critically revised the manuscript. JW contributed to analysis and interpretation and critically revised the manuscript. XX contributed to data acquisition and analysis. YZ and HL contributed to data acquisition. DL contributed to conception, design, interpretation, and drafted and critically revised the manuscript. All authors gave final approval and agree to be accountable for all aspects of the work.
Funding
This work was supported by the Clinical Research Center of Shandong University (Grant Number 2020SDUCRCA005) and Shandong University Postgraduate Education and Teaching Excellence Cultivation Program of China (Grant Number ZY2019004).
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher's Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fped.2021.718735/full#supplementary-material
Supplementary Figure 1. Experimental airway model with pressure taps locations.
Supplementary Figure 2. Comparison of wall static pressure between experiments and three different CFD models.
Supplementary Table 1. Statistical comparisons of morphological measurement.
Supplementary Table 2. Grid independence test.
Supplementary Table 3. Changes of airway volume and aerodynamic parameters before and after rapid maxillary expansion.
References
1. Ahn JC, Kim JW, Lee CH, Rhee CS. Prevalence and risk factors of chronic rhinosinusitus, allergic rhinitis, and nasal septal deviation: results of the Korean National Health and Nutrition Survey 2008–2012. JAMA Otolaryngol Head Neck Surg. (2016) 142:162–7. doi: 10.1001/jamaoto.2015.3142
2. Yildirim I, Okur E. The prevalence of nasal septal deviation in children from Kahramanmaras, Turkey. Int J Pediatr Otorhinolaryngol. (2003) 67:1203–6. doi: 10.1016/j.ijporl.2003.07.005
3. Berger G, Hammel I, Berger R, Avraham S, Ophir D. Histopathology of the inferior turbinate with compensatory hypertrophy in patients with deviated nasal septum. Laryngoscope. (2000) 110:2100–5. doi: 10.1097/00005537-200012000-00024
4. Chen XB, Lee HP, Chong VF, Wang De Y. Assessment of septal deviation effects on nasal air flow: a computational fluid dynamics model. Laryngoscope. (2009) 119:1730–6. doi: 10.1002/lary.20585
5. Kim SK, Heo GE, Seo A, Na Y, Chung SK. Correlation between nasal airflow characteristics and clinical relevance of nasal septal deviation to nasal airway obstruction. Respir Physiol Neurobiol. (2014) 192:95–101. doi: 10.1016/j.resp.2013.12.010
6. Moss ML, Salentijn L. The primary role of functional matrices in facial growth. Am J Orthod. (1969) 55:566–77. doi: 10.1016/0002-9416(69)90034-7
7. Young T, Finn L, Kim H. Nasal obstruction as a risk factor for sleep-disordered breathing. The University of Wisconsin Sleep and Respiratory Research Group. J Allergy Clin Immunol. (1997) 99:S757–62. doi: 10.1016/S0091-6749(97)70124-6
8. Haas AJ. Rapid expansion of the maxillary dental arch and nasal cavity by opening the midpalatal suture. Angle Orthod. (1961) 31:73–90.
9. Lo Giudice A, Leonardi R, Ronsivalle V, Allegrini S, Lagravere M, Marzo G, et al. Evaluation of pulp cavity/chamber changes after tooth-borne and bone-borne rapid maxillary expansions: a CBCT study using surface-based superimposition and deviation analysis. Clin Oral Investig. (2021) 25:2237–47. doi: 10.1007/s00784-020-03539-3
10. Maspero C, Galbiati G, Perillo L, Favero L, Giannini L. Orthopaedic treatment efficiency in skeletal class III malocclusions in young patients: RME-face mask versus TSME. Eur J Paediatr Dent. (2012) 13:225–30.
11. Abdalla Y, Brown L, Sonnesen L. Effects of rapid maxillary expansion on upper airway volume: a three-dimensional cone-beam computed tomography study. Angle Orthod. (2019) 89:917–23. doi: 10.2319/101218-738.1
12. Monini S, Malagola C, Villa MP, Tripodi C, Tarentini S, Malagnino I, et al. Rapid maxillary expansion for the treatment of nasal obstruction in children younger than 12 years. Arch Otolaryngol Head Neck Surg. (2009) 135:22–7. doi: 10.1001/archoto.2008.521
13. Fastuca R, Meneghel M, Zecca PA, Mangano F, Antonello M, Nucera R, et al. Multimodal airway evaluation in growing patients after rapid maxillary expansion. Eur J Paediatr Dent. (2015) 16:129–34.
14. Doruk C, Sokucu O, Sezer H, Canbay EI. Evaluation of nasal airway resistance during rapid maxillary expansion using acoustic rhinometry. Eur J Orthod. (2004) 26:397–401. doi: 10.1093/ejo/26.4.397
15. Ghoneima A, Albarakati S, Jiang F, Kula K, Wasfy T. Computational fluid dynamics analysis of the upper airway after rapid maxillary expansion: a case report. Prog Orthod. (2015) 16:10. doi: 10.1186/s40510-015-0085-x
16. Iwasaki T, Saitoh I, Takemoto Y, Inada E, Kanomi R, Hayasaki H, et al. Improvement of nasal airway ventilation after rapid maxillary expansion evaluated with computational fluid dynamics. Am J Orthod Dentofacial Orthop. (2012) 141:269–78. doi: 10.1016/j.ajodo.2011.08.025
17. Aziz T, Ansari K, Lagravere MO, Major MP, Flores-Mir C. Effect of non-surgical maxillary expansion on the nasal septum deviation: a systematic review. Prog Orthod. (2015) 16:15. doi: 10.1186/s40510-015-0084-y
18. Farronato G, Giannini L, Galbiati G, Maspero C. RME: influences on the nasal septum. Minerva Stomatol. (2012) 61:125–34.
19. D'Ascanio L, Lancione C, Pompa G, Rebuffini E, Mansi N, Manzini M. Craniofacial growth in children with nasal septum deviation: a cephalometric comparative study. Int J Pediatr Otorhinolaryngol. (2010) 74:1180–3. doi: 10.1016/j.ijporl.2010.07.010
20. Holton NE, Yokley TR, Figueroa A. Nasal septal and craniofacial form in European- and African-derived populations. J Anat. (2012) 221:263–74. doi: 10.1111/j.1469-7580.2012.01533.x
21. Cingi C, Muluk NB, Ulusoy S, Lopatin A, Sahin E, Passali D, et al. Septoplasty in children. Am J Rhinol Allergy. (2016) 30:e42–47. doi: 10.2500/ajra.2016.30.4289
22. Sabry O, Dewidar H, Aziz MA, Elemam A, Nassar A. Efficiency of modified Goldman's technique in open pediatric septoplasty. Egypt J Otolaryngol. (2021) 37:16. doi: 10.1186/s43163-021-00074-0
23. Saniasiaya J, Abdullah B. Quality of life in children following nasal septal surgery: a review of its outcome. Pediatr Investig. (2019) 3:180–4. doi: 10.1002/ped4.12145
24. Bailie N, Hanna B, Watterson J, Gallagher G. An overview of numerical modelling of nasal airflow. Rhinology. (2006) 44:53–7.
25. Burgos MA, Sanmiguel-Rojas E, Del Pino C, Sevilla-Garcia MA, Esteban-Ortega F. New CFD tools to evaluate nasal airflow. Eur Arch Otorhinolaryngol. (2017) 274:3121–8. doi: 10.1007/s00405-017-4611-y
26. Li L, Zang H, Han D, Ramanathan M Jr, Carrau RL, London NR Jr. Impact of a concha bullosa on nasal airflow characteristics in the setting of nasal septal deviation: a computational fluid dynamics analysis. Am J Rhinol Allergy. (2020) 34:456–62. doi: 10.1177/1945892420905186
27. Ramanathan M, Ramesh P, Aggarwal N, Parameswaran A, Sailer HF, George AE. Evaluation of airflow characteristics before and after septoplasty in unilateral cleft patients with a deviated nasal septum: a computational fluid dynamics study. Int J Oral Maxillofac Surg. (2021) 50:451–6. doi: 10.1016/j.ijom.2020.07.032
28. Xu C, Sin S, McDonough JM, Udupa JK, Guez A, Arens R, et al. Computational fluid dynamics modeling of the upper airway of children with obstructive sleep apnea syndrome in steady flow. J Biomech. (2006) 39:2043–54. doi: 10.1016/j.jbiomech.2005.06.021
29. Iwasaki T, Yanagisawa-Minami A, Suga H, Shirazawa Y, Tsujii T, Yamamoto Y, et al. Rapid maxillary expansion effects of nasal airway in children with cleft lip and palate using computational fluid dynamics. Orthod Craniofac Res. (2019) 22:201–7. doi: 10.1111/ocr.12311
30. Miner RM, Al Qabandi S, Rigali PH, Will LA. Cone-beam computed tomography transverse analysis. Part I: normative data. Am J Orthod Dentofacial Orthop. (2012) 142:300–7. doi: 10.1016/j.ajodo.2012.04.014
31. Setlur J, Goyal P. Relationship between septal body size and septal deviation. Am J Rhinol Allergy. (2011) 25:397–400. doi: 10.2500/ajra.2011.25.3671
32. Xiong GX, Zhan JM, Jiang HY, Li JF, Rong LW, Xu G. Computational fluid dynamics simulation of airflow in the normal nasal cavity and paranasal sinuses. Am J Rhinol. (2008) 22:477–82. doi: 10.2500/ajr.2008.22.3211
33. Fleck RJ, Ishman SL, Shott SR, Gutmark EJ, McConnell KB, Mahmoud M, et al. Dynamic volume computed tomography imaging of the upper airway in obstructive sleep apnea. J Clin Sleep Med. (2017) 13:189–96. doi: 10.5664/jcsm.6444
34. Tan J, Huang J, Yang J, Wang D, Liu J, Liu J, et al. Numerical simulation for the upper airway flow characteristics of Chinese patients with OSAHS using CFD models. Eur Arch Otorhinolaryngol. (2013) 270:1035–43. doi: 10.1007/s00405-013-2363-x
35. Mylavarapu G, Murugappan S, Mihaescu M, Kalra M, Khosla S, Gutmark E. Validation of computational fluid dynamics methodology used for human upper airway flow simulations. J Biomech. (2009) 42:1553–9. doi: 10.1016/j.jbiomech.2009.03.035
36. Chen S, Wang J, Liu D, Lei L, Wu W, Liu Z, et al. Open oral cavity has little effects on upper airway aerodynamics in children with obstructive sleep apnea syndrome: a computational fluid dynamics study based on patient-specific models. J Biomech. (2021) 121:110383. doi: 10.1016/j.jbiomech.2021.110383
37. Iwasaki T, Takemoto Y, Inada E, Sato H, Suga H, Saitoh I, et al. The effect of rapid maxillary expansion on pharyngeal airway pressure during inspiration evaluated using computational fluid dynamics. Int J Pediatr Otorhinolaryngol. (2014) 78:1258–64. doi: 10.1016/j.ijporl.2014.05.004
38. Lindemann J, Keck T, Wiesmiller K, Sander B, Brambs HJ, Rettinger G, et al. Nasal air temperature and airflow during respiration in numerical simulation based on multislice computed tomography scan. Am J Rhinol. (2006) 20:219–23. doi: 10.1177/194589240602000220
39. Sul B, Wallqvist A, Morris MJ, Reifman J, Rakesh V. A computational study of the respiratory airflow characteristics in normal and obstructed human airways. Comput Biol Med. (2014) 52:130–43. doi: 10.1016/j.compbiomed.2014.06.008
40. Inthavong K, Shang Y, Tu J. Surface mapping for visualization of wall stresses during inhalation in a human nasal cavity. Respir Physiol Neurobiol. (2014) 190:54–61. doi: 10.1016/j.resp.2013.09.004
41. Huh D, Fujioka H, Tung YC, Futai N, Paine R III, Grotberg JB, et al. Acoustically detectable cellular-level lung injury induced by fluid mechanical stresses in microfluidic airway systems. Proc Natl Acad Sci USA. (2007) 104:18886–91. doi: 10.1073/pnas.0610868104
42. Hur JS, Kim HH, Choi JY, Suh SH, Baek SH. Investigation of the effects of miniscrew-assisted rapid palatal expansion on airflow in the upper airway of an adult patient with obstructive sleep apnea syndrome using computational fluid-structure interaction analysis. Korean J Orthod. (2017) 47:353–64. doi: 10.4041/kjod.2017.47.6.353
43. Ribeiro AN, De Paiva JB, Rino-Neto J, Illipronti-Filho E, Trivino T, Fantini SM. Upper airway expansion after rapid maxillary expansion evaluated with cone beam computed tomography. Angle Orthod. (2012) 82:458–63. doi: 10.2319/030411-157.1
44. Gorgulu S, Gokce SM, Olmez H, Sagdic D, Ors F. Nasal cavity volume changes after rapid maxillary expansion in adolescents evaluated with 3-dimensional simulation and modeling programs. Am J Orthod Dentofacial Orthop. (2011) 140:633–40. doi: 10.1016/j.ajodo.2010.12.020
45. Lee B. Timing of Peak Mandibular Growth in Different Facial Growth Patterns and Resultant Mandibular Projection. Master of Science, Graduate Department of Dentistry University of Toronto (2011).
46. Schendel SA, Jacobson R, Khalessi S. Airway growth and development: a computerized 3-dimensional analysis. J Oral Maxillofac Surg. (2012) 70:2174–83. doi: 10.1016/j.joms.2011.10.013
Keywords: maxillary transverse deficiency, nasal septal deviation (NSD), rapid maxillary expansion (RME), nasal resistance, nasal aerodynamics, computational fluid dynamics (CFD)
Citation: Chen S, Wang J, Xi X, Zhao Y, Liu H and Liu D (2022) Rapid Maxillary Expansion Has a Beneficial Effect on the Ventilation in Children With Nasal Septal Deviation: A Computational Fluid Dynamics Study. Front. Pediatr. 9:718735. doi: 10.3389/fped.2021.718735
Received: 01 June 2021; Accepted: 06 December 2021;
Published: 10 February 2022.
Edited by:
Kostas N. Priftis, National and Kapodistrian University of Athens, GreeceReviewed by:
Jingliang Dong, RMIT University, AustraliaAntonino Lo Giudice, University of Catania, Italy
Copyright © 2022 Chen, Wang, Xi, Zhao, Liu and Liu. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Dongxu Liu, liudongxu@sdu.edu.cn