Can CFD establish a connection to a milder COVID-19 disease in younger people? Aerosol deposition in lungs of different age groups based on Lagrangian particle tracking in turbulent flow

To respond to the ongoing pandemic of SARS-CoV-2, this contribution deals with recently highlighted COVID-19 transmission through respiratory droplets in form of aerosols. Unlike other recent studies that focused on airborne transmission routes, this work addresses aerosol transport and deposition in a human respiratory tract. The contribution therefore conducts a computational study of aerosol deposition in digital replicas of human airways, which include the oral cavity, larynx and tracheobronchial airways down to the 12th generation of branching. Breathing through the oral cavity allows the air with aerosols to directly impact the larynx and tracheobronchial airways and can be viewed as one of the worst cases in terms of inhalation rate and aerosol load. The implemented computational model is based on Lagrangian particle tracking in Reynolds-Averaged Navier–Stokes resolved turbulent flow. Within this framework, the effects of different flow rates, particle diameters and lung sizes are investigated to enable new insights into local particle deposition behavior and therefore virus loads among selected age groups. We identify a signicant increase of aerosol deposition in the upper airways and thus a strong reduction of virus load in the lower airways for younger individuals. Based on our findings, we propose a possible relation between the younger age related fluid mechanical protection of the lower lung regions due to the airway size and a reduced risk of developing a severe respiratory illness originating from COVID-19 airborne transmission.


Introduction
Flows with dispersed particles are of great interest, since they can be found in numerous fields of engineering and medical science. In this context, the ongoing COVID-19 pandemic can be considered as a recent topic. The causative agent of the disease is known as SARS-CoV-2 that primarily targets the respiratory tract and can manifest mostly as pneumonia when affecting the lower respiratory tract and less likely as acute respiratory distress syndrome (ARDS) [13,32]. The virus is known to be transmissible via contacts and droplets as well as aerosols [32,38]. Duguid [10] showed that by a single sneeze, infected hosts can easily generate up to a few million contaminated droplets and aerosols. In many studies aerosols are defined as particles with a diameter d p ≤ 5 µ m and larger particles as droplets [11,36]. However, there have been some suggestions that widen that range and postulate particles up to an aerodynamic diameter of 10−20 µ m as aerosols, due to their ability to linger prolongedly in the air and reach deeper in the lung, [7]. Thomas et al. [33] stated that larger particles mostly deposited in the upper airways whereas small aerosols are more prone to bypass the mechanical defense mechanisms and penetrate deeper into the lung, causing typical disease profiles related to the alveolar region. After the host released the contaminated aerosols, the SARS-CoV-2 virus has been found to remain viable for 3 h, increasing the risk for airborne spread of COVID-19 [7]. On November 8th, 2020, approximately one year after the outbreak, COVID-19 affects a significant proportion of people with 49.7 million reported cases and 1.2 million deaths worldwide [37]. However, the SARS-COV-2 virus is not yet fully understood and there are still high uncertainties remaining. This involves the size distribution of droplets and aerosols, the aerosol viral load and the minimum number of inhaled SARS-CoV-2 viruses that are required to infect an individual [26]. Moreover, it is still unknown why there is an unequal distribution and course of infections among the population as it is less frequently diagnosed among children [6]. As mentioned by Thomas et al. [33], the biometry and therefore deposition profiles are affected by factors like age, body weight, breathing mode, gender and health state. This leads to the question, whether pure mechanical effects, like biometry and breathing mode of children, render a decisive difference in aerosol deposition behavior and consequently virus load compared to grown ups.
Given the fact that in vivo and in vitro experiments are mainly limited due to human safety and image resolution [19], computational fluid dynamics (CFD) can play an important role to enable new insights in this field. In this paper we first validate our numerical lung model setup by comparison with the in vitro and in silico benchmark case of Koullapsis et al. [18]. Moreover, we discuss the effects of age-related lung sizes on the deposition of cough generated aerosols and therefore virus doses in human airways. This is chosen to investigate the correlation of age and regional aerosol deposition accessing the SARS-CoV-2 pandemic from a fluid mechanical point of view. The general purpose of this paper is to provide new insights in age-related aerosol deposition.
The paper is organized as follows: In Sect. 2, the governing equations are reviewed. Moreover, the lung model is validated with respect to benchmark models for the case of a human lung replica in Sect. 3. Additionally, Sect. 4 contains the setup of the present lung sizes to model aerosol deposition in various human age-groups as well as results and discussion of local aerosol load in different airway sizes. Finally, Sect. 5 summarizes the paper and presents the main conclusions.

Airway geometry
There is a variety of airway geometries that is used in CFD. Early models employed simplified artificial geometries mainly basing on the symmetric model of Weibel et al. [34] In recent years, medical imaging enabled a detailed view of the human airway and provided more realistic replicas [19]. The geometry considered in this paper, see Fig. 1b, is the same as in the benchmark case of Koullapsis et al. [18] and is originally used in [3,4,15,23]. The benchmark case consists of a simplified lung model that was adopted to measure the regional deposition ratios. The realistic airway geometry and the simplified airway model used in the benchmark case are shown in Fig. 1. More details related to the lung model generation can be reviewed in [18].

Flow field
We employ the Reynolds Averaged Navier Stokes (RANS) equations using the k-ω-SST model in order to study the steady state flow in the replica of the human lung. The governing RANS equations for an incompressible fluid, which are solved using OpenFOAM [28,35], are given by with representing the mean viscous stress tensor components and denoting the Reynolds stresses [12]. Hereū andp describe the Reynolds-averaged fluid velocity components and pressure and ρ denotes the fluid density. Moreover, the -sign represents fluctuations. Additional body forces can be found inf D . The closure problem to constitutively express τ RANS can be solved by employing approximate turbulence models like the k-ω-SST model [12]. OpenFOAM uses the finite volume method to discretise the above equations.

Particles
Dispersed flows can be captured by the Lagrangian-Eulerian approach. Hence, the Navier-Stokes equations are solved for the continuous phase (air) as mentioned in Sect. 2.2 and Fig. 1 a Original airway, b benchmark case, c segmented model [18] the motion of particles (aerosols) in fluid is described in a Lagrangian manner [14]. To this end, a set of ordinary differential equations is evaluated along the particle trajectory to obtain the change of particle location and motion. These equations are given by Newton's second law and render for a spherical particle point mass: where x p is the particle's position vector, u p the particle's velocity, m p = ρ p d 3 p π/6 the mass of the spherical particle and F i represents the sum of forces acting on the particle. [14] In this study, the aerosol dimensions are estimated to be at the scale of 1 ≤ d p ≤ 10 µ m. Hence, the forces containing the major influence on the particle trajectory are the drag force F D , the buoyancy force F B and the gravitational force F G . Other forces like Brownian motion, added mass, and Basset history force are considered as negligible [14]. Thus, the force balance equation simplifies to: where m p , V p , ρ p are the mass, volume and density of the particle, respectively, ρ f denotes the fluid density and g is the gravitational acceleration. The drag force for spherical particles in OpenFOAM is implemented as follows: [14] where d p is the particle diameter and C D is the drag coefficient. OpenFOAM uses the following empirical relation for the drag coefficient [14]: where Re r := ρ f d p |u p −u|/μ denotes the particle Reynolds number based on the relative velocity [8]. The buoyancy and gravitational force are usually combined and computed jointly as [17]: The equations are solved with the icoUncoupledKinematic-ParcelFoam 1 solver of OpenFOAM .

Turbulent dispersion model
To account for the interaction of the particles with the turbulent eddies, the instantaneous velocity u =ū + u of the fluid is required. As this field is not accessible from the RANS equations and only the averaged velocityū is available, additional models are required to properly estimate the fluctuation velocity. Here we employ the OpenFOAM model Stochas-ticDispersionRAS [16] where u is computed to disturb the velocity field in a random direction, with a Gaussian distribution of zero mean and variance σ [18]. It relates the velocity fluctuations to the turbulent kinetic energy k as follows  where d is an additional random vector, ξ denotes random numbers with zero mean and unit variance of Gaussian distribution and k is the turbulent kinetic energy [14]. The model assumes isotropic turbulence, rendering the standard deviation σ as with u 1 , u 2 , u 3 describing the velocity coefficients in Cartesian coordinates [14,16].

Limitations
To model the flow of aerosols, we restricted ourselves in the present study to the following conditions: -dilute flow allowing for one-way coupling of particles and fluid, -assumption of isotropic turbulence, -steady state flow field.
In addition, the size of the aerosols under investigation (1 ≤ d p ≤ 10 µ m) is sufficiently small, so their surface tension is strong enough to solely behave like small spherical rigid particles [2].

Numerical setup
The simplified airway model used in the benchmark case is shown in Fig. 1b. The benchmark analysis consists of an in vitro experiment and five different in silico approaches. To evaluate the results of the present model, the benchmark LES1 simulation was selected as the main reference as it showed good agreement with the in vitro results [18]. Furthermore the benchmark study demonstrated that the RANS1 result was the most accurate RANS simulation compared to RANS2 and RANS3. Hence, both LES1 and RANS1 are used for comparison with the present model. Their numerical setup is provided in Table 1. In addition, the mesh statistics are presented in Table 2.

Flow field
In Fig. 2 the contours of mean velocity magnitude and turbulent kinetic energy are compared in the central sagittal plane of the airways mouth-throat region. It is evident that a good qualitative agreement of the mean velocity distribution is achieved for the present mode (Mesh 2) compared to the LES1. The Y-shaped velocity pattern of the LES1 simulation can be resolved with the present model (Mesh 2) approach, whereas it is absent in the RANS1 velocity field. The level of turbulent kinetic energy for both RANS cases is lower in comparison to the LES1 results. However, the present model (Mesh 2) is able to capture the main characteristics e.g. maximum of the turbulent kinetic energy in the throat region close to the (D1-D2) cross-section. The discrepancies between RANS methods are likely due to the lower grid resolution of the RANS1 model and the absence of boundary layers.
In the next step a more detailed analysis of the velocity and kinematic turbulent energy fields is conducted by comparing the respective profiles at selected cross-sections. The results are displayed in Fig. 3. Here also the PRESENT (Mesh 1) cross-sections are included, which were evaluated on a significantly coarser grid. The location of the cross-sections is highlighted in Fig. 2, with exact locations estimated from Koullapsis et al. [18]. In the (A1-A2) section, which represents the inlet region, the turbulent velocity profile can be captured by all numerical setups, except the PRESENT (Mesh1) setting that employed a fixed inlet flow rate. More-over, the typically low velocity profiles in the low mouth and pharynx region, depicted in (B1-B2) and (C1-C2) respectively, can be reproduced by all methods. In addition, the present models are able to capture the effects in the acceleration region (D1-D2) of the pharynx. Besides, the (E1-E2) (F1-F2) cross-sections of the present model achieve good agreement with the LES1 and RANS1 reference results. Even though the exact locations of the reference profiles were not provided in [18], the mean velocity profiles achieved low discrepancies compared to the LES1 and RANS1 cross-sections. In the case of turbulent kinetic energy, the deviation between the present models and the LES1 results slightly increases and the RANS results lead, mainly in the upper mouththroat region (B1-B2) (C1-C2), to an under-prediction. As observed by Koullapsis et al. [18], larger discrepancies in the profiles of the turbulent kinetic energy occur in the upper regions of the airway across simulations. This deviation can partially be linked to different inlet conditions as well as mesh resolution among LES1 and RANS models. A further cause is the pronounced anisotropic nature of the flow in this region, a consequence of some significant changes in the airway's geometry. Capturing the anisotropy effects is a known deficiency of the RANS based model, which can significantly be improved by applying LES models, however at a much higher computational cost. Since the k-ω-SST model still performed relatively well compared to the highend LES validation case, capturing the main characteristics of the turbulent kinetic energy, its reasonable computational cost led to its choice in the framework of the targeted parametric analysis. Additionally, it is indicated that the impact of the different inlet conditions of PRESENT (Mesh 1) on the turbulence kinetic energy as well as the velocity, is low downstream of the mouth region.

Verification of particle deposition
To verify the particle deposition behavior, the present lung model is compared to the LES1 and RANS1 deposition frac-tions of Koullapsis et al. [18] for three different flow rates and a diameter range of 1 ≤ d p ≤ 10 µ m. According to Koullapsis et al. [18], the particles are considered to be diethylhexyl sebacate particles in ambient air temperature (ρ p = 914 kg/m 3 ) and are distributed uniformly across the inlet. Details of the particle-insertion are presented in Table 3. A thorough description of the reference simulations that are compared in Fig. 4 is provided in [18]. In the following, the human airway is subdivided in three different regions of interest. These are the mouth-throat, tracheobronchial tree as well as a combination of both, which is referred to as the overall region. The mouth-throat region includes the oral-cavity and trachea. The tracheobronchial tree describes the human lung more downstream, excluding the collectors, which represent the lower airway regions. Figure 4 displays the deposition behavior of the selected particle range with an inhalation of 60 l/min in the mouth-throat, tracheobronchial tree and overall lung region. A wide range of results is generated by employing the different methods (LES1-3, RANS1-3) in the reference [18]. It is obvious that the deviation between the present model  (PRESENT (Mesh 2)) and the benchmark LES1 results in the overall airway geometry is small, compared to the addition provided method. Moreover, it is shown that the numerical results in the mouth-throat as well as the tracheobronchial tree are respectively in a good agreement, proving the validity of the present model. However, in the tracheobronchial region, the deviation increases by decreasing the particle diameter. In the mouth-throat region, in contrary, the difference between the present model and the reference vanishes towards smaller particles. In Fig. 5

Aerosol deposition in lungs of different age groups
In this section the verified lung model is further employed to account for aerosol deposition in different lung sizes. The main interest is to evaluate possible variations in aerosol deposition of different particle sizes in human airways

Model generation
Due to the limited availability for sufficient lung models of different age groups, the airway geometries of children are generated by scaling the adult lung volume of [18]. The scaling factor is set to represent the ratio of the total volume capacity (TLC) of the considered age group compared to the TLC of a grown-up. The TLC describes the 'maximum volume of air that the lungs can hold after a maximum inspiration' [1]. For adults, the TLC is taken to be 5.7 l as mentioned by Aung et al. [1]. To estimate the TLC of children, the relationship proposed by Lyons et al. [24] is considered: where H is the height in cm and W is the mass in kg. The flow rates are estimated by taking the average breathing frequency of the considered age groups times the tidal volume (TV). The TV describes the 'volume of air that is breathed in and out in a single quiet breath' [1]. The average TV of a resting male adult is estimated as 630 ml [30]. To obtain the TV of children the following estimation for mechanical ventilation is employed [9]: The studied characteristic children, which represent the average weight and height of the given ages, are provided in Table 4. In the present study, the set-up of the simulation is the same as in Sect. 3.1. In Addition, the computational domain is scaled by the TLC ratio. For the youngest child (Age 1) we estimate a mean inlet velocity ofŪ inlet = 0.74 m/s whereasŪ inlet = 0.33 m/s is obtained for the male adult. Considering the different inlet areas, we achieve an average (inlet) Reynolds number of Re = 392 for the child (Age 1) and Re = 436 for the adult. Therefore, the flow characteristics are approximately similar across the simulations. The flow field and turbulent kinetic energy in the central sagittal plane are displayed among various age-groups in Fig. 6. Observe the different scaling of these quantities in Fig. 6 despite otherwise similar patterns of U and k.

Modeling of aerosol
The particles are considered to be cough-generated aerosols with a density of ρ p = 1704 kg/m 3 as proposed by [21]. Moreover, a quantity of 100, 000 aerosols are distributed randomly at the inlet and released over a time span of 5 µ s. Furthermore, the initial parcel/particle velocity (U 0 ) is set to zero. Besides, the considered aerosol range is 1 − 10 µ m as discussed in Sect. 1. In addition, the aerosols are assumed to stick to the airway wall once a particle comes into contact with this boundary, to mimic the mucus layer on the inner walls of the airways [18]. Furthermore, the particle tracking time step was set to 10 µ s for the male adult and children (Age 13-5) as well as 5 µ s for child models (Age 3-1), ensuring a maximal particle Courant number of Co p ≤ 1.0.

Influence of lung size on aerosol deposition
Firstly, after 2.5 s the deposition fractions of the mouththroat, the tracheobronchial tree and the cumulative results (overall) were investigated. Figure 7a-c displays the aerosol deposition fractions over particle size for seven different lung volumes. The cumulative results, displayed in Fig. 7a, show a clear trend for all considered lung sizes. The deposition fraction of all artificial individuals grows with increasing aerosol diameter. However, Fig. 7b indicates a more moderate increase in the mouth-throat region. Therefore, the main impact of the particle diameter on the aerosol deposition occurs in the tracheobronchial tree, which is presented in Fig. 7c. The effect of particle size influence coincides with Sect. 3.3. Additionally, aerosols with a diameter of d p = 1 µ m behave similarly for all lung dimensions and reach an identical deposition fraction in the mouth-throat as well as in the tracheobronchial tree. Consequently, these aerosols experience the lowest deposition rate, indicating that aerosols with d p = 1 µ m are more prone to travel deeper into the lung. However, Fig. 7a, c, indicate that the aerosol deposition frac- tion varies significantly among the considered age groups for particles with d p > 3 µ m. In the mouth-throat region, presented in Fig. 7b, the aerosol deposition is mainly unaffected by the airway size, as the deviations between the individuals are comparatively small. In contrary, for the younger subjects, a notable higher amount of aerosols deposits in the tracheobronchial tree, see Fig. 7c, achieving an overall higher deposition rate. We relate a higher deposition in the upper airways (mouth-throat region and tracheobronchial tree) to a lower amount of aerosols that reach the collectors that represent the alveolar region. Therefore there are less contagious particles that could reach the alveolar region of the lung, as highlighted in Fig. 7d. This trend continues for d p = 10 µ m aerosols, where the highest aerosol deposition is observed. For this diameter the youngest individual (Age 1) deposits approximately 90% of inserted aerosols, resulting in the lowest amount of aerosols and thus virus load that travel further into the collectors, which represent the smaller airway regions. For all considered age groups and diameters, the lung size has a major impact on the deposition effect in the tracheobronchial tree, whereas the mouth-throat region remains mainly unaffected. As stated in [5] there is a close correlation of SARS-CoV-2 viral load in the lower airways and severity in COVID-19 ARDS. In this context, we conjecture that the lower aerosol and consequently lower virus dose in the alveolar region among the children models is related to an increase of the probability of a mild infection of COVID-19. Furthermore, by increasing the age of the subject and thus the lung volume, the deposition fractions approach the results of the adult reference case.
A more detailed view of the local deposition is provided in Fig. 8, where the aerosol deposition fraction is displayed over each lung model component. The allocation of these segment identifiers (segment-ID) to the lung model is provided in Fig. 1c. The component-wise deposition fractions in Fig. 8 underline the effect that the particle deposition fractions of different diameters in the mouth-throat region are nearly identical for all age groups. In addition, Fig. 8b-d indicates the correlation between an increasing deposition effect of the upper airways and a reduced airway dimension. The segmental deposition is nearly identical for all considered age groups, for particles with a diameter of d p = 1 µ m. The aerosols with a diameter d p = 10 µ m are more prone to deposit in the upper airway regions. For the child (Age 1) that is inhaling aerosols with a diameter of d p = 10 µ m most of the particles deposit in the first 13 segments. This effect is further visualized in Fig. 9 which presents the cumulative aerosol deposition fraction along each airway branch for d p = 10 µ m across six age-groups. The cumulative aerosol deposition fraction along each airway branch was calculated by summing fractional deposition in the current and all preceding segments. It is highlighted that the younger the considered individual the higher the cumulative deposition fraction in upper airway generations. Due to this increase, the deposition fraction in the lower airway regions is significantly reduced and less aerosols reach the collectors. This effect is mitigated if smaller particles are considered as presented in Fig. 8. Furthermore Fig. 8 displays that the deposition increases with larger aerosol diameter for all individuals. In all cases, the mouth-throat region has a high deposition rate, however, a difference between different age groups becomes clear in the tracheobronchial tree for the particles with d ≥ 5 µ m, as with the decreasing of the age of the individual and consequently the lung size, the particle deposition significantly increases.
In our study we relate the volume of air that is breathed in and out in a single breath (tidal volume T V ) in resting condition to the estimated weight of the considered child, see Eq. 14. This leads to a T V of 72 ml for the child (Age 1) compared to T V = 630 ml for adults. The breathing frequency of these age-groups are estimated as f Age1 = 32 breaths/min and f Adult = 12 breaths/min respectively. Compared to the adult a child has therefore a higher breathing frequency as well as a reduced tidal volume. We observe, that the higher the frequency-volume ratios the higher the deposition in the upper airways. This consequently leads to a significant reduction of aerosols that can penetrate deeper into the lung. These findings relate exclusively to person's normal activity, excluding exercise activity. More details on the average T V and breathing frequencies during the active phase as well as the regeneration phase of physical activities would be needed to assess the safety of exercising.

Time dependent deposition
The time dependent deposition of aerosol particles with d p = 1, 5, 10 µ m is compared across seven different airway dimensions. The aim is to estimate the percentage of deposited particles after one inhalation. To this end, the mean inhalation time is estimated as half of the inverse of the breathing frequency: The average breathing frequency f avg and the estimated mean inhalation time t inh are provided in Table 5. Figure 10 displays the resulting deposition fractions over a timespan of 2.5 s. Figure 10a-c represents the deposition behavior for d p = 1 µ m aerosols. It underlines the findings of Sect. 4.3.1, which state that the considered age groups reach a similar stationary deposition fraction in the mouththroat as well as the tracheobronchial tree for small aerosols with d p = 1 µ m. Moreover, Fig. 10a-c indicates a different time dependent deposition. Due to the higher flow rates inside the smaller lungs, the particles tend to reach their final depo-  sition state earlier. Figure 10d-f presents the time dependent deposition for d p = 5 µ m aerosols. For this aerosol dimension, a significant increase of the deposition rate is observed for the mouth-throat as well as the tracheobronchial tree, reaching an overall filtering of approximately 20−30 %. It is shown that the lung size impacts the deposition rate as well as the time dependent deposition. In Fig. 10a-i it is visible that all particles are already deposited after the time needed for one inhalation (Table 5), this indicates, that all non-deposited particles reach the collectors.

15 min: Inhalation
In the final step, the deposited aerosols after 15 mininhalation are investigated. This time span was chosen, according to the recommendations of the German Health Departments, which predict a high risk, if an individual had distinct cumulative face-to-face contact with a host for at least 15 min [29].
These results are obtained by taking the particle deposition fractions presented in Sect. 4.3.1 times the inhaled aerosols for each subject after 15 min. To estimate the aerosol load after 15 min, we assumed that an infected individual released contagious aerosol by a single cough in the inhalation region (1 m × 1 m × 1 m) of the subject. Lindsley et al. [22] measured 900-302,200 particles/cough while subjects had influenza. Therefore, an average aerosol load in a cubic room of l aerosols = 150, 000/m 3 is considered. The volume that the specific subjects inhales in the critical time Q 15min is provided in Table 6. The equation for the inhalation load of aerosols n 15min after 15 min renders n 15min = l aerosols × Q 15min .
(16) Table 6 shows that n 15min , the total amount of inhaled aerosols after 15 min, varies significantly across the considered age-stages. The inhalation dose is strongly increasing with age and leads for the adult to approximately one third more aerosols compared to the adolescent (Age 13) and   Figure 11a-d presents the deposited particles over the particle diameter for the seven different lung sizes. Moreover, Fig. 11d shows that the trend of Sect. 4.3.1 is more pronounced. This is due to the effect of increased deposition in the upper airway for younger individual and the likewise reduction of inhalation aerosol load. For the child (Age 1) the aerosol load that reached into the collectors and therefore penetrate deeper into the lung is less than one third for aerosols with d p = 1 µ m and less than one tenth for d p = 10 µ m. For older children models the difference is mitigated, but remains notable.

Conclusions
The simulation of regional aerosol deposition in human airways of different sizes, corresponding to different age groups, is an important step to gather knowledge about the unequal distribution of COVID-19 infection, especially among younger age groups, also from a fluid mechanical point of view. In order to conduct this research, one of the current limitations was the lack of detailed lung models for different age groups. To solve the problem the considered lung geometries were generated by a dedicated scaling of the adult lung of Koullapsis et al. [18], and then used in computational models to study general trends in deposition of aerosols. The generated lung models covered artificial models of children aged 1−13 years as well as the adult model of Koullapsis et al. [18]. The computational model is based on RANS equations with k-ω-SST model employed to account for the turbulent flow in the airways.
With the Lagrangian particle tracking the StochasticRan-domRAS model of OpenFOAM is used to account for the impact of turbulent eddies on the particle trajectories. The benchmark LES1 case of Koullapsis et al. [18] served for verification of the implemented computational model. The present model leads to computational results with reasonably small differences to the reference, rendering it as suitable for investigations of the impact of different lung sizes on aerosol deposition. Aerosol deposition results are obtained and presented for different lung sizes and particle dimensions. Significant variability in regional aerosol deposition is observed across the considered airway models. However, the differences mitigated towards the smaller aerosols. In this context, deposition was found to be particularly sensitive in the tracheobronchial tree and less impacted in the mouththroat region among the considered age groups. In addition,

(c) (d)
a general trend is observed, which indicated a higher deposition of aerosols in the upper airways for younger individuals and therefore a reduced deposition in the lower airways. A higher alveolar virus load in the lower airways is conjectured to trigger lower respiratory tract symptoms, like pneumonia or acute respiratory distress syndrome, which are also known to increase morbidity of COVID-19 patients [5,33]. We conclude that a higher aerosol deposition in the upper airways of children, mainly in the tracheobronchial tree, leads to a significant reduction of virus load in the lower airways. We connect this effect to a higher chance of developing mild to moderate respiratory illness. Regarding our results, we propose a possible relation between the age related mechanical protection of the lung, due to airway size, and the risk of severe respiratory illness originating from COVID-19 airborne transmission. The presented computational study is applicable also to a more complete upper airway geometry case, that would also include the nasal cavity, with possibility to study also other inhalation regimes.
Funding Open Access funding enabled and organized by Projekt DEAL.

Conflict of interest None.
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://creativecomm ons.org/licenses/by/4.0/.