Sample size estimation
Maximum midsagittal velocity of the airflow in the upper airway is considered to be the primary outcome variable based on a previous study by Feng et al [27]. A sample size of 30 will be needed to ensure a 80% power to reject the null hypothesis at a significance level of 5%, assuming differences in maximum midsagittal velocity and its standard deviation is 0.7 m/s and 0.5 m/s between the cases with an AN ratio less than 0.6 and equal or more than 0.6, and the ratio between the two groups is 2:1.
Samples collection
This cross-sectional study is a subset of a longitudinal prospective study that was performed at Dalian Stomatological Hospital between 2015 - 2017, in which 2D and 3D image modalities were compared for tracing anatomic landmarks before and after orthodontic treatment. The study was approved by the regional ethics review boards in Dalian, China (DLKQLL201604) and Bergen, Norway (2018/1547 REK Vest). Informed consent was obtained from all patients or their legal guardians. The baseline images from 2015 were retrospectively collected and employed in the current study. The inclusion criteria were individuals aged 9 to 15 years who had had both a lateral cephalogram and CBCT scan examined within one week. For CBCT images, the field of view was required to cover the UA regions, including the nasal cavity, nasopharynx, and oropharynx. The exclusion criteria were severe maxillofacial abnormalities and previous surgery on skeletal and soft tissue related to respiration. In the present study, ninety-two cases were initially included. X.F previewed all the CBCT scans and lateral cephalograms. Fifty-seven cases were excluded, of which 53 did not cover the UA, 3 scans had motion artefact and 1 showed suboptimal patient positioning. Eventually, thirty-five cases were recruited. All the cases were divided into two groups: Group 1 with AN ratios < 0.6 (n=25) and Group 2 with AN ratios ≥ 0.6 (n=10).
Lateral cephalogram
The AN ratios were measured and calculated on the lateral cephalograms captured by a digital pan/ceph system (ORTHOPHOS XG 5; Sirona Dental Systems, Bensheim, Germany) at 73 kVp and 15 mA with exposure times of 9.4 seconds and a contrast resolution of 16-bit depth. A is a perpendicular distance between the point of maximal convexity of the adenoid to the anterior margin of the basiocciput. N is the distance between the posterosuperior edge of the hard palate and the anteroinferior edge of the spheno-occipital synchondrosis [9] (Figure 1).
CBCT scans
All CBCT scans were obtained by 3D eXam (KaVo, Biberach an der Riss, Germany). The following parameters were used: a field of view (FOV) of 16 × 13 cm, tube voltage of 120 kV, tube current of 5 mA, scanning time of 14.7 seconds, voxel size of 0.2 mm, and contrast resolution of 14-bit depth.
CFD simulation
The CBCT images were imported in the digital imaging and communications in medicine (DICOM) format to MIMICS software (MIMICS, Materialise, Belgium) for later analysis. 3D renderings of the CBCT scans were oriented with axial planes parallel to the Frankfurt horizontal plane; the midsagittal planes intersected the nasion and anterior nasal spine; and the coronal plane was adjusted to the level of the porions. For each case, a mask was reconstructed, making sure the integrity of UA displayed correctly. CFD simulation was then conducted on the 3D model within the mask region. The superior boundary of the studied UA was defined as a vertical plane, in the nasal cavity, passing through the most posterior point of the middle turbinate, whereas the inferior boundary was a horizontal plane, in the pharynx, in line with the most anterior-inferior point of cervical vertebra 4. Each end of the boundary was extended by 20mm to avoid flow reversing during the simulating process. The inlet and outlet of UA were set on the extended planes. A surface model was then created according to the extended 3D model for mesh generation. We chose tetrahedral and prismatic cells to construct the main body and boundary layer of the UA mesh (ANSYS, Inc., Canonsburg, Pennsylvania). The SST κ-ω model was used to calculate the aerodynamic characteristics of UA by applying ANSYS Fluent (ANSYS, Inc., Canonsburg, Pennsylvania). The wall of the UA was defined as no-slip, stationary, and rigid. The temperature and density of air were set as fixed. At inspiration, the inlet was set with the pressure 0 Pa and the outlet at a flow rate of - 200 mL/s [28]. The corresponding values were - 200 mL/s and 0 Pa at inlet and outlet at expiration.
Data analyses
The aerodynamic parameters applied and computed throughout the CFD simulation are listed in Table 1. The pressure drop (△P) refers to the pressure difference between a vertical plane through the most posterior point of the middle turbinate and a horizontal plane through the tip of the epiglottis. Data were processed using IBM-SPSS, version 25.0 (IBM, New York, NY, USA). Significance was set at p-values less than 0.05. The assumption of normal distribution for all variables were tested. An independent-samples T test or Mann-Whitney U test was used for comparing each variable according to genders as well as groups. The Intraclass Correlation Coefficient (ICC) was applied for testing the intra- and inter-observer reliability on the selected ten cases using a random number generator.
Table 1
Description of the aerodynamic parameters evaluated applying the CFD simulation
Name
|
|
Unit
|
Definition
|
Maximum Vms
|
m/s
|
The maximum velocity on the midsagittal plane
|
|
△P
|
Pa
|
The pressure drop of airflow between the defined two planes
|
Maximum Pws
|
Pa
|
The maximum lateral pressure of airflow acting on the UA wall
|
Minimum Pw
|
Pa
|
The minimum vertical pressure of airflow acting on the UA wall
|