
==== Front
Sci Rep
Sci Rep
Scientific Reports
2045-2322
Nature Publishing Group UK London

38605156
58758
10.1038/s41598-024-58758-5
Article
Effects of decongestion on nasal cavity air conditioning efficiency: a CFD cohort study
Xiao Qiwei 12
Bates Alister J. 123
Doorly Denis J. d.doorly@imperial.ac.uk

4
1 https://ror.org/01hcyya48 grid.239573.9 0000 0000 9025 8099 Center for Pulmonary Imaging Research, Cincinnati Children’s Hospital Medical Center, Cincinnati, OH USA
2 https://ror.org/01hcyya48 grid.239573.9 0000 0000 9025 8099 Division of Pulmonary Medicine, Cincinnati Children’s Hospital Medical Center, Cincinnati, OH USA
3 https://ror.org/01e3m7079 grid.24827.3b 0000 0001 2179 9593 Department of Pediatrics, University of Cincinnati, Cincinnati, OH USA
4 https://ror.org/041kmwe10 grid.7445.2 0000 0001 2113 8111 Department of Aeronautics, Imperial College London, South Kensington Campus, London, UK SW7 2AZ
11 4 2024
11 4 2024
2024
14 848216 2 2024
2 4 2024
© The Author(s) 2024
2024
https://creativecommons.org/licenses/by/4.0/ Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
Decongestion reduces blood flow in the nasal turbinates, enlarging the airway lumen. Although the enlarged airspace reduces the trans-nasal inspiratory pressure drop, symptoms of nasal obstruction may relate to nasal cavity air-conditioning. Thus, it is necessary to quantify the efficiency of nasal cavity conditioning of the inhaled air. This study quantifies both overall and regional nasal air-conditioning in a cohort of 10 healthy subjects using computational fluid dynamics simulations before and after nasal decongestion. The 3D virtual geometry model was segmented from magnetic resonance images (MRI). Each subject was under two MRI acquisitions before and after the decongestion condition. The effects of decongestion on nasal cavity air conditioning efficiency were modelled at two inspiratory flowrates: 15 and 30 L min−1 to represent restful and light exercise conditions. Results show inhaled air was both heated and humidified up to 90% of alveolar conditions at the posterior septum. The air-conditioning efficiency of the nasal cavity remained nearly constant between nostril and posterior septum but dropped significantly after posterior septum. In summary, nasal cavity decongestion not only reduces inhaled air added heat by 23% and added moisture content by 19%, but also reduces the air-conditioning efficiency by 35% on average.

Subject terms

Translational research
Computational biophysics
Engineering
Biomedical engineering
http://dx.doi.org/10.13039/501100000266 Engineering and Physical Sciences Research Council EP/M506345/1 Bates Alister J. issue-copyright-statement© Springer Nature Limited 2024
==== Body
pmcIntroduction

Nasal cavity air-conditioning protects our delicate lung tissues from harm by heating and humidifying inhaled air close to alveolar conditions. The inner surface of the nasal cavity is covered with mucosa and these tissues serve as the heat and water source, typically supplying approximately 360 kcal of heat and 250–400 ml of water per day during the air-conditioning process1,2. A recent study has shown that heat and water transfer at the nasal surface is responsible for the perception of nasal airflow and when air conditioning is limited due to pathology, patients report the sensation of nasal obstruction, regardless of airflow3.

Investigations into nasal cavity air-conditioning have been conducted through both in-vivo and computational methods. Notably, Keck et al. mapped the spatial distribution of air conditioning within the nose using in-vivo measurements over a cohort of 50 healthy volunteers and found the inhaled air temperature significantly increased within the anterior nasal cavities and the humidification process was almost synchronous with the heating process4,5. Computational Fluid Dynamics (CFD) simulations revealed that even in extreme external temperatures, the nose is capable of heating inhaled air to near body temperature6. However, earlier CFD studies of nasal cavity air conditioning have used simplified models; others such as Lindemann et al. and Pless et al., did not incorporate water vapor transport8,9. Garcia et al. investigated the nasal cavity heat and water vapor transport between control and pathological subjects by assuming both constant cavity surface temperature and relative humidity7. A two-film theory was applied in simulating the heating of nasal cavity air by Kumahata et al. and Hanida et al. that evaluated the effects of latent heat, while the water transfer model assumed fully saturated conditions at the wall of the nasal cavity8. More recently, a study by Inthavong et al. applied a two-film theory to both the heat and water model in simulating the nasal cavity air-conditioning over a limited number of subjects9.

The nasal cycle, a natural physiological process, alternately decongests one side of the nasal cavity approximately every two hours, leading to anatomical variations that can affect air-conditioning efficiency10. Prior studies often utilized single geometric models per subject, based on the nasal anatomy at the time of imaging, although it has been shown considering the airway dynamics can make big difference in CFD predicted resutls9,11–19.

In this study, we aim to quantify the influence of anatomical variation on the efficiency of nasal air conditioning. By applying a nasal decongestant to both sides of the nasal cavity, we mimicked not only an extreme state of the nasal cycle using decongestant, representing the most decongested state possible, but also the use of decongestants is relevant for their practical application in temporarily alleviating nasal obstruction. This approach allowed us to measure the impact of nasal anatomy on air conditioning efficiency under both normal and decongested conditions in a paired sample study.

Materials and methods

Subjects

A cohort of 10 healthy volunteers aged between 21 and 38 were recruited. Ethical approval was granted by National Health Service Health Research Authority NRES Committee Southeast Coast-Surrey with reference number 06/Q0602/18. The image acquisition protocol and the translation of magnetic resonance images (MRI) to 3D virtual airway surfaces were the same as in previously published work20.

Informed consent was obtained from each of the subjects included in this study. All methods were carried out in accordance with relevant guidelines and regulations.

CFD simulations and boundary conditions

CFD simulations

The CFD software STAR-CCM+ (Siemens PLM Software, Plano, TX) was used to model the heat and water transported by the airflow through the nasal cavity. The flow computations were performed for inspiratory flowrates at 15 and 30 L min−1 using the methodology as described more fully in previous work. Briefly, the airflow through the nasal cavity was consider incompressible, so the mass conservation formula over a control volume can be written in the following form:∇·ρu→=0

where ∇ is divergence operator, u→ (m/s) means the velocity vector, ρ(kgm-3) is flow density. The momentum conservation can be expressed as follows:ρ∂u→∂t+u→·∇u→=-∇p+μ∇2u→

where μ (Pa s) is dynamic viscosity and ∇2 is Laplace operator. Both flowrates were simulated using a laminar model, because the flow characteristic in the majority of the nasal cavity is mainly laminar as shown by multiple previous studies20–24. The effects of gravity and buoyancy were neglected. The numerical schemes used for the temporal term is second order backward differentiation which involves the current and previous two time levels, and for the convection term the second-order upwind scheme was applied.

A mass flow inlet was specified at a far-removed external domain inlet and a pressure outlet condition (0 Pa gauge pressure) was specified at an extended nasopharyngeal outlet20. This setup of the boundary conditions is suitable because the pressure difference between locations, not the absolute pressure is the metric of interest in this study given the known inlet flowrate. There are 4 million total CFD polyhedral volume cells for each CFD simulation, and each simulation has 7 prism layers near the airway wall to account for both high flow and temperature gradient. The total thickness of prism layer is 0.3 mm, and the first layer is 0.01 mm with a stretch factor of 1.7 to allow a smooth transition of cells in between. Figure 1A demonstrates the mesh details at plane 1V (defined in Fig. 2) that was used for all the simulations in this study. More detailed mesh scenes can be found in the previous study on nasal airflow distribution in this cohort, as well as in the appendix section on mesh independence results20. Results using our model for heat and water transfer are compared with previous in-vivo measurements in the results section.Figure 1 (A) A visual representation of polyhedral volume mesh created using STAR-CCM + on plane 1V (defined in Fig. 2) with the setting that gives 4 million volume meshes in total. (B) Diagram of the thermal and moisture conditions at the nasal cavity wall. TO represents the temperature at the surface between nasal tissue and mucosa, which is assumed at constant temperature of 34 °C. Ts represents the temperature at the surface between the mucosa and the air. Ta is the temperature of the air at the center of the first CFD volume mesh cell adjacent to the wall. Membrane represents a membrane layer which has the same properties as water.

Figure 2 Left side shows the perspective view of a half meter sphere that was attached to the face, and extrusion at the end of nasopharynx with length of 50 times diameter. The right side shows the anatomy of a subject specific nasal cavity with cross-sectional planes as gray line overlay. There are several anatomical landmarks labelled by blue lines: Nostril, nasal valve (NV), anterior septum (AS), posterior septum (PS), and nasopharynx (Naso). 1V, 2V, 3V are three equal spaced landmarks between AS and PS. The overall centerline length was indicated by red dashed line. The normalized length of red dashed line scales between 0 and 1.

Heat and water vapor boundary conditions

The transport of heat and water vapor by the airflow may be described respectively by:1 ∂T∂t+u→·∇T=kρCp∇2T

2 ∂Yωv∂t+u→·∇Yωv=D∇2Yωv

where T(K) is temperature, u→(m/s) is the flow velocity, k(Wm-1K-1) thermal conductivity, ρ(kgm-3) density, CP(Jkg-1K) specific heat capacity, ∇2 is the Laplace operator, Yωv the water vapor fraction and D(cm-2s-1) the diffusion coefficient of water vapor in air. At the mucosal wall, a two-film model was employed to describe the exchange of heat and water with the mucosal surface.

Figure 1B shows schematically the modelling of heat and water vapor transfer from the nasal cavity wall into the inhaled air. The layer labeled “Organ” represents the nasal capillary bed, “Membrane” represents the layer of liquid mucus and the upper portion of the water supplying mucosal tissue on the nasal wall; the first wall layers represent the airway lumen as discretized in the CFD mesh (a detailed CFD mesh scene is shown in the Fig. 1A). The organ layer is responsible for providing the heat source in the model. This heat must then flow through the membrane to the air. The membrane layer is assumed to have the same specific heat capacity as water25. The temperature at the interface (Ts) between the exposed surface of the nasal wall and the cavity air is determined by the resultant heat flux. The inhaled air temperature is represented by Ta, whilst To refers to the temperature of the capillary bed, which is assumed at 34 °C based on previous studies26–28. The species boundary condition at the interface is assumed to be fully saturated at the local surface temperature.

The total amount of heat flux into the air (q˙air) is the sum of the flux from the membrane (q˙membrane) and the latent heat flux (q˙latent) carried by the water during evaporation (or condensation):q˙air=q˙latent+q˙membrane

The fluxes q˙air, q˙latent and q˙memb are calculated using the equations below:q˙air=-kairTa-Tsδwall

q˙latent=LH2O·ω˙

q˙membrane=-kmembraneTs-Toδmembrane

where ω˙(kgm-2s-1) is the water mass flux between the fully saturated wall and the cavity air, and δwall(m) is the half thickness of the first wall layer of the CFD mesh as labelled in Fig. 1B. The latent heat of water LH2O(Jg-1) is defined by the formula29:LH2O=2500.8-0.00006·Ts3+0.0016·Ts2-2.36·Ts,

where Ts is in Celsius. Using Fick’s law of diffusion, ω˙ is calculated as:ω˙=-ρfluid·D·∂Ywv∂δwall

where ρfluid is the density of the mixture of air and water vapor. The mass fraction of water vapor under fully saturated conditions was quantified from an empirical fit of psychrometric data30 by the formula below.Ywv=2·10-5Ts2+0.0003Ts+0.0025

Heat flux and water flux into the air were determined by the surface temperature Ts, and were iteratively updated so that heat and water fluxes at the membrane and at the air were balanced at each time step. Combining the above equations, the surface temperate Ts can be calculated as:Ts=kairTaδmembrane+kmembraneToδwall-LH2Oω˙δwallδmembranekmembraneδwall+kairδmembrane

For this cohort study, computations were performed with Ta set to 25 °C, 50% relative humidity (corresponding to a mass fraction of 0.01125 and To = 34 °C). The rendering of both the applied temperature and moisture boundary conditions can be found in appendix.

Geometric definition

The segmented 3D surface from the MRI scan of each subject includes anatomy from the face to the end of the nasopharynx. In addition, for each subject, an external half sphere (diameter = 0.5m) was attached to the face to ensure natural inflow profile and an extrusion (order of 50 diameter) was added from the end of nasopharynx to prevent reverse flow. The velocity of air approaching the nose during inhalation reduces rapidly, and an extension of this size is more than sufficient to ensure a natural flow profile develops at the entrance to the nose21. Although we have not investigated the minimum diameter needed, a smaller extension of order of 50mm diameter appears more than sufficient for modelling inhalation21. However, the extended domain facilities other expiratory flow studies and requires few additional elements. Extruding the outlet allows flow to stabilize and reduces the impact of artificial outlet boundary condition on the flow field. Figure 2 illustrates the reconstructed 3D surface of one subject with the anatomical landmarks used, as defined in our previous flow only study20.

Measurements

Air temperature

The mean air temperature was calculated over each cross-sectional plane (Fig. 2 right) from anterior to posterior cavity in each subject; in addition, the cohort mean air temperature at each corresponding plane was also calculated.

Relative humidity and moisture content

Relative humidity (RH) is the ratio of the mass fraction of water vapor mixture to the mass fraction of water vapor if fully saturated, at the corresponding temperature and pressure. Lindemann et al., reported that RH can be 90% at 32 °C in the nasopharynx, however, RH is a temperature and pressure dependent metric, and it is not ideal for inter-subject comparisons, given the absolute amount of water vapor is the metric of interest. Therefore, moisture concentration was compared to alveolar conditions (100% RH at 37 °C), rather than local conditions, using a new quantity defined as the ratio between mass fraction of H2O at local conditions and the mass fraction of H2O at alveolar conditions, which was named moisture content (MC). The equation below shows the calculation of MC where Yt is the local mass fraction of water vapor, while Y37∘C refers to the mass fraction of water vapor at alveolar conditions.MC=100·YtY37∘C

Thus, MC is suitable to quantify how much the water vapor is in the air compared to alveolar conditions, irrespective of local conditions.

Intra and inter-subject variations in MC were plotted in a similar manner to the temperature variations described above.

Film coefficients

To quantify the efficiency of heat and water transfer between the nasal cavity wall and inhaled air, heat and mass film coefficients were introduced. The equation below shows the calculation of the heat film coefficient. The left side of the equation Q(W) represents the heat flux while the right side of the equation is the product of heat film coefficient h(Wm2K), nasal cavity surface area A(m2) and temperature gradient at the wall (ΔT):Q=h·A·ΔT

Analogously to heat transfer, the mass film coefficient was calculated using the formula below, indicating water vapor flux is a product of mass film coefficient hc(ms-1), nasal cavity surface area A(m2) and gradient of mass fraction of water vapor at the wall (ΔcA):ΩA=hcAΔcA

Data analysis

The Wilcoxon signed-rank test was performed to evaluate the statistical significance of decongestion-associated changes in nasal cavity air-conditioning efficiency31. Following convention, a 95% confidence level was chosen to distinguish significant from non-significant results32. Box and whisker plots were used to illustrate the results.

Results

Comparison with in-vivo measurements

Keck et al. measured the temperature and humidity profile at multiple positions within nasal cavities of 50 subjects at 25 °C room temperature and 35% RH under restful breathing condition4. They reported the mean temperature and RH at the nasal valve, anterior turbinate and nasopharynx as shown in Table 1. To compare with the results of Keck et al., we conducted simulations under the same boundary temperature and humidity conditions and using 15 L min−1 as a representative restful breathing flow rate. Virtual probes were created at the anatomical similar locations in the CFD simulation. Given uncertainty in matching locations and the associated error ranges, results show reasonable agreement.Table 1 The top 3 rows of the table show both mean temperature and relative humidity results over the cohort of 10 subjects at normal condition at three distinct anatomical locations in nasal cavity.

Exemplar subject	Nasal valve	Anterior turbinate	Nasopharynx	
Temperature (°C)	27.9	28.9	32.4	
Relative humidity (%)	68	83	98	
Keck et al	Nasal valve	Anterior turbinate	nasopharynx	
Temperature (°C)	28.9 ± 2.3	30.3 ± 1.6	32.6 ± 1.5	
Relative humidity (%)	69.0 ± 6.5	78.7 ± 7.2	90.3 ± 5.3	
The bottom 3 rows show the same information but from experimental measurements by Keck et al.5.

Temperature and moisture content predictions

Figure 3a–d, shows results from CFD simulations of inhaled air temperature increase from anterior to posterior nasal cavity averaged over the cohort at both normal and decongested conditions. All simulations shared the same thermal and humidity inflow boundary conditions (25C, 50% RH). The solid line represents the mean value over the cohort and shaded area indicates the variations of one standard deviation.Figure 3 (a)–(d) The mean inhaled air temperature increase from nostril to nasopharynx over the cohort at both normal and decongested conditions with different flowrates as labelled. Panel e and f show the Wilcoxon-signed rank test result of temperature differences between nostril and nasopharynx at normal and decongested conditions at 15 and 30 L min−1 flowrates.

At 15 L min−1 flowrate (Fig. 3a,b), decongestion reduced the amount of heat added to the inhaled air by 24% compared to normal condition on average for the cohort, with the final temperature at nasopharynx reduced from 31.8 to 30.5 °C After doubling the flowrate, decongestion reduced the same measurement of added heat by 22%, with temperature at nasopharynx reduced from 30.3 to 29.4 °C . Figure 3e,f showed the reduction in temperature to be statistically significant with a median reduction by 1.3 °C at flowrate of 15 L min−1, and 0.9 °C at flowrate of 30 L min−1 accordingly.

Similarly, in Fig. 4 from a–d the solid line indicates the mean value of the computed MC distribution in the nasal cavity over the cohort, with the shaded area representing variation of one standard deviation. Results show the MC added to the inhaled air was reduced by decongestion, respectively by 18% at 15 L min−1 flowrate and by 20% for 30 L min−1 flowrate at nasopharynx plane. Applying the Wilcoxon signed rank test confirmed the reductions in added MC to be statistically significant, with median reduction by 20% and 23% respectively. More detailed temperature and MC rendering in the nasal cavities are shown in the appendix.Figure 4 (a)–(d) The mean inhaled air moisture content increase from nostril to nasopharynx over the cohort at both normal and decongested conditions with different flowrates as labelled. Panel e and f show the Wilcoxon-signed rank test result of moisture content differences between nostril and nasopharynx at normal and decongested conditions at 15 and 30 L min−1 flowrates.

Heat and mass film coefficient

The heat film coefficient, h, is a measure of heating efficiency of the nose. Since the vast majority of heating occurs from the nostril to the posterior septum, the heat film coefficient was calculated over this region. Table 2 shows that on average across all subjects, the heat film coefficient is reduced by decongestant, as might be expected as a consequence of mean airflow speed reduction concomitant with increased cavity lumen area.Table 2 Comparison of normal and decongested cohort-mean heat and mass film coefficients at two inspiratory flowrates.

Heat film coefficient h (W m−2 K−1), mean value of region between nostril and posterior septum	
Flowrate (L min−1)	Normal	Decongested	Change in h due to decongestion	
15	125.4	78.0	–38%	
30	174.0	112.1	–36%	
Change in h due to increase in flowrate	40%	40%		
Mass film coefficient hc (m s−1), mean value of region between nostril and posterior septum	
Flowrate (L min−1)	Normal	Decongested	Change in hc due to decongestion	
15	0.061	0.039	–36%	
30	0.073	0.052	–29%	
Change in hc due to increase in flowrate	20%	33%		

Similarly, the mass film coefficient, hc, is a measure of the humidification efficiency of the nose. The mass film coefficient follows the same trends as the heat film coefficient in terms of changes with decongestion and flowrate. However, the changes in hc are less than those for h, particularly those due to flowrate increases.

Discussion

CFD was used to investigate the air conditioning capabilities of the nasal cavity in a cohort of 10 subjects, using in-vivo acquired data pre- and post-decongestion. The air-conditioning performance of the nasal cavity varies with the nasal cycle, and it is of both physiological and clinical interest to investigate the magnitude of the performance decrement caused by decongestion to further our understanding of nasal function.

Results indicated considerable variation in performance across the cohort, in both normal and decongested states. The distribution of modelled mucosal surface temperature (Fig. 3a–d) from anterior to posterior cavity show an initial rapid increase followed by a levelling off after the posterior septum as the surface area for heat exchange reduces. Similar trend is likewise found for the moisture content (MC). The computational modelling predicted a decongestion-associated reduction of added heat to inhaled air by 24% (15 L min−1) and 22% (30 L min−1) respectively. Similarly, the added MC to inhaled air were reduced by 18% (15 L min−1) and 20% (30 L min−1). All the decongestion-associated changes observed in this paired-sample study were found to be statistically significant.

The film coefficient is a parameter to quantify the efficiency of heat or mass transfer from nasal cavity inner wall to inhaled air. The results show decongestion can reduce the heat transfer efficiency by 37% and mass transfer efficiency by 33% for both lower and higher flowrates. In addition, doubling the flowrate can increase the heat (mass) transfer efficiency by 40% (25%) both before and after decongestion, which very likely due to relatively higher temperature gradient caused by fast moving flow near the cavity wall.

Limitations of the present study are firstly, due to the limitation of available computational power, only constant flowrates were simulated in this study. Future study should use subject specific breathing profile to increase the fidelity of the simulations. In addition, the simplified heat and water boundary conditions also can introduce errors especially by assuming fully saturated air at the nasal cavity wall and flow interface. Further work should focus on improving the boundary conditions as well as considering a far larger cohort than the present 10 subjects. It would also be of interest to extend the methodology reported here to consider pathological nasal anatomies. Lastly, the air conditioning performance of the nose depends on the external environment. An indication of the scale of this dependence is shown in the appendix for increasingly demanding external regimes and at an elevated flow rate of 30 L min−1. In the present work, we have studied 10 pairs of different nasal geometries at two flowrates to determine ranges of variation in air conditioning capability under the same external conditions. Extending such cohort studies to cover such broad range of conditions as in the appendix is a goal of further research, but brings with it a much enlarged burden of associated computational and data reduction requirements.

Conclusion

In summary, the majority of the air conditioning (90%) process was completed before nasal posterior septum. The effect of decongestion can significantly reduce the final state of conditioned air at the level of the nasopharynx. On average, doubling the flowrate in the normal nose also lowers the air conditioning capability to similar levels as observed with decongestion. Specifically, the air conditioning efficiency of the nasal cavity can be reduced by 35% on average (including heating and humidifying) after decongestion, whilst doubling the flowrate can reduce the efficiency by 33% on average. This study has quantified how the changes in the nasal anatomy up to the nasal septum, which can be caused by decongestion, affects the conditioning of inhaled air.

Supplementary Information

Supplementary Information.

Supplementary Information

The online version contains supplementary material available at 10.1038/s41598-024-58758-5.

Author contributions

Q.X. conducted the computational fluid dynamics simulations and wrote the manuscript. Q.X. D.D. and A.B. analyzed the data. A.B and D.D. revised the manuscript. All authors reviewed the manuscript.

Funding

Funding was provided by the Engineering and Physical Sciences Research Council (Grant No. EP/M506345/1). The Imperial College Open Access Fund supported publication costs.

Data availability

The datasets generated during and/or analyzed during the current study are available from the corresponding author on reasonable request.

Competing interests

The authors declare no competing interests.

Publisher's note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
==== Refs
References

1. Walker JEC Wells RE Heat and water exchange in the respiratory tract Am. J. Med. 1961 30 259 267 10.1016/0002-9343(61)90097-3 13782729
Walker, J. E. C. & Wells, R. E. Heat and water exchange in the respiratory tract. Am. J. Med. 30, 259–267 (1961).13782729 10.1016/0002-9343(61)90097-3
2. Cole P Further observations on the conditioning of respiratory air J. Laryngol. Otol. 1953 67 669 681 10.1017/S0022215100049161 13109356
Cole, P. Further observations on the conditioning of respiratory air. J. Laryngol. Otol. 67, 669–681 (1953).13109356 10.1017/S0022215100049161
3. Zhao K Regional peak mucosal cooling predicts the perception of nasal patency Laryngoscope 2014 124 589 10.1002/lary.24265 23775640
Zhao, K. et al. Regional peak mucosal cooling predicts the perception of nasal patency. Laryngoscope 124, 589 (2014).23775640 10.1002/lary.24265
4. Keck T Leiacker R Heinrich A Kuhnemann S Rettinger G Humidity and temperature profile in the nasal cavity Rhinology 2000 38 167 171 11190750
Keck, T., Leiacker, R., Heinrich, A., Kuhnemann, S. & Rettinger, G. Humidity and temperature profile in the nasal cavity. Rhinology 38, 167–171 (2000).11190750
5. Keck T Leiacker R Riechelmann H Rettinger G Temperature profile in the nasal cavity Laryngoscope 2000 110 651 654 10.1097/00005537-200004000-00021 10764013
Keck, T., Leiacker, R., Riechelmann, H. & Rettinger, G. Temperature profile in the nasal cavity. Laryngoscope 110, 651–654 (2000).10764013 10.1097/00005537-200004000-00021
6. Naftali S Schroter RC Shiner RJ Elad D Transport phenomena in the human nasal cavity: A computational model Ann. Biomed. Eng. 1998 26 831 839 10.1114/1.108 9779956
Naftali, S., Schroter, R. C., Shiner, R. J. & Elad, D. Transport phenomena in the human nasal cavity: A computational model. Ann. Biomed. Eng. 26, 831–839 (1998).9779956 10.1114/1.108
7. Garcia GJM Atrophic rhinitis: A CFD study of air conditioning in the nasal cavity PLoS ONE 2007 27709 1082 1092
Garcia, G. J. M. et al. Atrophic rhinitis: A CFD study of air conditioning in the nasal cavity. PLoS ONE 27709, 1082–1092 (2007).
8. Hanida S Influence of latent heat in the nasal cavity* J. Biomech. Sci. Eng. 2013 8 209 224 10.1299/jbse.8.209
Hanida, S. et al. Influence of latent heat in the nasal cavity*. J. Biomech. Sci. Eng. 8, 209–224 (2013).10.1299/jbse.8.209
9. Inthavong K Fletcher DF Khamooshi M Vahaji S Salati H Wet surface wall model for latent heat exchange during evaporation Int. J. Numer. Method Biomed. Eng. 2022 38 e3581 10.1002/cnm.3581 35142094
Inthavong, K., Fletcher, D. F., Khamooshi, M., Vahaji, S. & Salati, H. Wet surface wall model for latent heat exchange during evaporation. Int. J. Numer. Method Biomed. Eng. 38, e3581 (2022).35142094 10.1002/cnm.3581
10. Pendolino AL Scarpa B Ottaviano G Relationship between nasal cycle, nasal symptoms and nasal cytology Am. J. Rhinol. Allergy 2019 33 644 10.1177/1945892419858582 31219310
Pendolino, A. L., Scarpa, B. & Ottaviano, G. Relationship between nasal cycle, nasal symptoms and nasal cytology. Am. J. Rhinol. Allergy 33, 644 (2019).31219310 10.1177/1945892419858582
11. Issakhov A Zhandaulet Y Abylkassymova A Issakhov A A numerical simulation of air flow in the human respiratory system for various environmental conditions Theor. Biol. Med. Model 2021 18 2 10.1186/s12976-020-00133-8 33407610
Issakhov, A., Zhandaulet, Y., Abylkassymova, A. & Issakhov, A. A numerical simulation of air flow in the human respiratory system for various environmental conditions. Theor. Biol. Med. Model 18, 2 (2021).33407610 10.1186/s12976-020-00133-8
12. Kumahata K Mori F Ishikawa S Matsuzawa T Nasal flow simulation using heat and humidity models* J. Biomech. Sci. Eng. 2010 5 565 577 10.1299/jbse.5.565
Kumahata, K., Mori, F., Ishikawa, S. & Matsuzawa, T. Nasal flow simulation using heat and humidity models*. J. Biomech. Sci. Eng. 5, 565–577 (2010).10.1299/jbse.5.565
13. Gunatilaka CC Schuh A Higano NS Woods JC Bates AJ The effect of airway motion and breathing phase during imaging on CFD simulations of respiratory airflow Comput. Biol. Med. 2020 127 104099 10.1016/j.compbiomed.2020.104099 33152667
Gunatilaka, C. C., Schuh, A., Higano, N. S., Woods, J. C. & Bates, A. J. The effect of airway motion and breathing phase during imaging on CFD simulations of respiratory airflow. Comput. Biol. Med. 127, 104099 (2020).33152667 10.1016/j.compbiomed.2020.104099
14. Xiao Q Assessing changes in airflow and energy loss in a progressive tracheal compression before and after surgical correction Ann. Biomed. Eng. 2020 48 822 833 10.1007/s10439-019-02410-1 31792705
Xiao, Q. et al. Assessing changes in airflow and energy loss in a progressive tracheal compression before and after surgical correction. Ann. Biomed. Eng. 48, 822–833 (2020).31792705 10.1007/s10439-019-02410-1
15. Xiao Q The interaction between neuromuscular forces, aerodynamic forces, and anatomical motion in the upper airway predicts the severity of pediatric OSA J. Appl. Physiol. 2024 136 70 78 10.1152/japplphysiol.00071.2023 37942529
Xiao, Q. et al. The interaction between neuromuscular forces, aerodynamic forces, and anatomical motion in the upper airway predicts the severity of pediatric OSA. J. Appl. Physiol. 136, 70–78 (2024).37942529 10.1152/japplphysiol.00071.2023
16. Bates AJ Assessing the relationship between movement and airflow in the upper airway using computational fluid dynamics with motion determined from magnetic resonance imaging Clin. Biomech. (Bristol, Avon) 2019 66 88 96 10.1016/j.clinbiomech.2017.10.011 29079097
Bates, A. J. et al. Assessing the relationship between movement and airflow in the upper airway using computational fluid dynamics with motion determined from magnetic resonance imaging. Clin. Biomech. (Bristol, Avon) 66, 88–96 (2019).29079097 10.1016/j.clinbiomech.2017.10.011
17. Yu, S., Sun, X. Z. & Liu, Y. X. Numerical analysis of the relationship between nasal structure and its function. Sci. World J. 2014 (2014).
18. Gunatilaka CC Tracheomalacia reduces aerosolized drug delivery to the lung J. Aerosol. Med. Pulm. Drug Deliv. 2024 37 19 29 10.1089/jamp.2023.0023 38064481
Gunatilaka, C. C. et al. Tracheomalacia reduces aerosolized drug delivery to the lung. J. Aerosol. Med. Pulm. Drug Deliv. 37, 19–29 (2024).38064481 10.1089/jamp.2023.0023
19. Inthavong K Wen J Tu J Tian Z From CT scans to CFD modelling—fluid and heat transfer in a realistic human nasal cavity Eng. Appl. Comput. Fluid Mech. 2009 3 321 335
Inthavong, K., Wen, J., Tu, J. & Tian, Z. From CT scans to CFD modelling—fluid and heat transfer in a realistic human nasal cavity. Eng. Appl. Comput. Fluid Mech. 3, 321–335 (2009).
20. Xiao Q Bates AJ Cetto R Doorly DJ The effect of decongestion on nasal airway patency and airflow Sci. Rep. 2021 11 1 13 10.1038/s41598-021-93769-6 33414495
Xiao, Q., Bates, A. J., Cetto, R. & Doorly, D. J. The effect of decongestion on nasal airway patency and airflow. Sci. Rep. 11, 1–13 (2021).33414495 10.1038/s41598-021-93769-6
21. Doorly DJ Taylor DJ Gambaruto AM Schroter RC Tolley N Nasal architecture: Form and flow Philos. Trans. R. Soc. A Math. Phys. Eng. Sci. 2008 366 3225 3246 10.1098/rsta.2008.0083
Doorly, D. J., Taylor, D. J., Gambaruto, A. M., Schroter, R. C. & Tolley, N. Nasal architecture: Form and flow. Philos. Trans. R. Soc. A Math. Phys. Eng. Sci. 366, 3225–3246 (2008).10.1098/rsta.2008.0083
22. Doorly DJ Taylor DJ Schroter RC Mechanics of airflow in the human nasal airways Respir. Physiol. Neurobiol. 2008 163 100 110 10.1016/j.resp.2008.07.027 18786659
Doorly, D. J., Taylor, D. J. & Schroter, R. C. Mechanics of airflow in the human nasal airways. Respir. Physiol. Neurobiol. 163, 100–110 (2008).18786659 10.1016/j.resp.2008.07.027
23. Garcia GJM Dosimetry of nasal uptake of water-soluble and reactive gases: A first study of interhuman variability Inhal. Toxicol. 2009 21 607 618 10.1080/08958370802320186 19459775
Garcia, G. J. M. et al. Dosimetry of nasal uptake of water-soluble and reactive gases: A first study of interhuman variability. Inhal. Toxicol. 21, 607–618 (2009).19459775 10.1080/08958370802320186
24. Zhang Y Computational fluid dynamics (CFD) investigation of aerodynamic characters inside nasal cavity towards surgical treatments for secondary atrophic rhinitis Math. Probl. Eng. 2019 10.1155/2019/6240320
Zhang, Y. et al. Computational fluid dynamics (CFD) investigation of aerodynamic characters inside nasal cavity towards surgical treatments for secondary atrophic rhinitis. Math. Probl. Eng.10.1155/2019/6240320 (2019).10.1155/2019/6240320
25. Investigations of the thermal properties of human and animal tissues—Enlighten: Theses. https://theses.gla.ac.uk/1019/.
26. McFadden ER Direct recordings of the temperatures in the tracheobronchial tree in normal man J. Clin. Invest. 1982 69 700 705 10.1172/JCI110498 7061708
McFadden, E. R. et al. Direct recordings of the temperatures in the tracheobronchial tree in normal man. J. Clin. Invest. 69, 700–705 (1982).7061708 10.1172/JCI110498
27. McFadden ER Thermal mapping of the airways in humans J. Appl. Physiol. 1985 1985 58 564 570 10.1152/jappl.1985.58.2.564
McFadden, E. R. et al. Thermal mapping of the airways in humans. J. Appl. Physiol. 1985(58), 564–570 (1985).10.1152/jappl.1985.58.2.564
28. Hamdan AT Cherobin GB Voegels RL Rhee JS Garcia GJM Effects of mucosal decongestion on nasal aerodynamics: A pilot study Otolaryngol. Head Neck Surg. 2024 10.1002/OHN.713 38461407
Hamdan, A. T., Cherobin, G. B., Voegels, R. L., Rhee, J. S. & Garcia, G. J. M. Effects of mucosal decongestion on nasal aerodynamics: A pilot study. Otolaryngol. Head Neck Surg.10.1002/OHN.713 (2024).38461407 10.1002/OHN.713
29. Federer B A short course in cloud physics Pure Appl. Geophys. 1976 114 495 495 10.1007/BF00876948
Federer, B. A short course in cloud physics. Pure Appl. Geophys. 114, 495–495 (1976).10.1007/BF00876948
30. ASHRAE Psychrometric Chart no. 3, American Society of Heating, Refrigerating and Air-conditioning Engineers, 1992.
31. Rey, D. & Neuhäuser, M. Wilcoxon-signed-rank test. In International Encyclopedia of Statistical Science 1658–1659 (2011). 10.1007/978-3-642-04898-2_616
32. Hazra A Using the confidence interval confidently J. Thorac. Dis. 2017 9 4125 10.21037/jtd.2017.09.14 29268424
Hazra, A. Using the confidence interval confidently. J. Thorac. Dis. 9, 4125 (2017).29268424 10.21037/jtd.2017.09.14
