Test Data AC7-01: Difference between revisions
(7 intermediate revisions by the same user not shown) | |||
Line 141: | Line 141: | ||
of voxels analysed, and the volume of the region. The software was tested and compared | of voxels analysed, and the volume of the region. The software was tested and compared | ||
to commercial tools with an excellent agreement as documented by | to commercial tools with an excellent agreement as documented by | ||
[[Best_Practice_Advice_AC7-01# | [[Best_Practice_Advice_AC7-01#nesterov2009|Nesterov ''et al.'' (2009)]] | ||
and [[Best_Practice_Advice_AC7-01# | and [[Best_Practice_Advice_AC7-01#nesterov2009|Harms ''et al.'' (2014)]]. | ||
To facilitate the identification of the hot spots in the Carimas | To facilitate the identification of the hot spots in the Carimas | ||
software, the BGRY colour system was used to display the PET images and the system | software, the BGRY colour system was used to display the PET images and the system | ||
Line 401: | Line 401: | ||
The measured Deposition data reported in [[Test_Data_AC7-01#figure8|figure 8]] can be found | The measured Deposition data reported in [[Test_Data_AC7-01#figure8|figure 8]] can be found | ||
in | in file [[Media:AC7-01_DF_in-vitro.xlsx|DF_in-vitro.xlsx]]. | ||
<div id="figure9"></div> | <div id="figure9"></div> | ||
Line 413: | Line 413: | ||
---- | ---- | ||
{{ACContribs | {{ACContribs | ||
|authors=P. Koullapis, F. Lizal, J. Jedelsky, L. Nicolaou, K. Bauer, O. Sgrott, M. Jicha, M. Sommerfeld, S. C. Kassinos | |authors=P. Koullapis<sup>a</sup>, F. Lizal<sup>b</sup>, J. Jedelsky<sup>b</sup>, L. Nicolaou<sup>c</sup>, K. Bauer<sup>d</sup>, O. Sgrott<sup>e</sup>, M. Jicha<sup>b</sup>, M. Sommerfeld<sup>e</sup>, S. C. Kassinos<sup>a</sup> | ||
|organisation=University of Cyprus | |organisation=<br><sup>a</sup>Department of Mechanical and Manufacturing Engineering, University of Cyprus, Nicosia, Cyprus<br><sup>b</sup>Faculty of Mechanical Engineering, Brno University of Technology, Brno, Czech Republic<br><sup>c</sup>Division of Pulmonary and Critical Care, School of Medicine, Johns Hopkins University, Baltimore, USA<br><sup>d</sup>Institute of Mechanics and Fluid Dynamics, TU Bergakademie Freiberg, Freiberg, Germany<br><sup>e</sup>Institute Process Engineering, Otto von Guericke University, Halle (Saale), Germany | ||
}} | }} | ||
{{ACHeader | {{ACHeader |
Latest revision as of 11:11, 30 September 2021
Aerosol deposition in the human upper airways
Application Challenge AC7-01 © copyright ERCOFTAC 2019
Test Data
Overview of Tests
Major portions of this section were adopted from Lizal et al. (2015). The positron emission tomography (PET) method provides the best spatial resolution (among radiological methods). In addition to local deposition in the various sections, the deposition hot spots can also be evaluated. However, in comparison to the PET methodology, which is routinely applied to clinical examination, using this method in the in vitro design requires major modifications both in the aerosol preparation and, in particular, in the experiment evaluation approach. The method, based on PET and fulfilling the above mentioned criteria, is presented in the following.
The aerosol exposure procedure
It is a common practice to coat the inner surface of the model, especially when using solid particles, to prevent bouncing of the particles hitting the surface. Since we used liquid di-ethylhexyl sebacate (DEHS) particles, we did not need to coat the inner surface of the model. Another reason for the coating is to prevent surface wetting. In our case the exposure time was short (5 to 15 mins) and only small amounts of DEHS deposited on the walls, therefore the possible flooding of the surface was not an issue. Aerosol particles were generated by a TSI 3475 Condensation Monodisperse Aerosol Generator (CMAG) from TSI, Inc., which works on the controlled heterogeneous condensation principle. Vapours of a suitable material, specifically DEHS, condense by a controlled method on small sodium chloride particles serving as the condensation nuclei. The advantage of DEHS is that it is not hydrophilic and does not evaporate, resulting in a constant size of generated particles. In a standard operating mode, the generator can produce particles with aerodynamic diameters within the 0.1 to 8μm range. The density of DEHS used for the experiments was 0.914 g/cm3 at 25°C. Radioactive aerosol particles were needed for the PET measurement of deposition. Therefore, the solution in the atomizer of the generator had to be tagged by a suitable radioactive substance. Fluorine 18 was the logical Choice of the positron emitter, being easily available at the cooperating PET center and possessing a suitable half-life (109 minutes). A solution of fluorine 18 in the form of fluoride ions was prepared by irradiation of H2180 enriched water on an IBA Cyclone 18/9 cyclotron (irradiation time 25 min, integrated irradiation current 11 μAh) at the UJV Rez's PET Centre in Brno. The irradiated water was transferred by a capillary transport system to a shielded dispensing box, where the fluorine 18 ions were captured on an ion—exchange resin (AG1—X8, BioRad) column and subsequently eluted with 300 ml of 10% sodium Chloride solution, followed by 1Â ml of water for injection. The resulting solution was repeatedly diluted with water for injection until the desired initial radioactivity was achieved. The CMAG was modified for the deposition measurement by using PET so that the atomizer vessel was accommodated in a protective lead container to shield of ionizing radiation. The atomizer was filled with a sodium Chloride solution containing 18F at an initial activity of 2.5 GBq. The concentration of sodium Chloride solution was 20 mg/L. The experimental rig is shown in Figure 5.
Figure 5: A scheme of the experimental setup during the PET measurement of aerosol deposition. |
The generated aerosol was fed through a 85Kr based NEKR-10 charge equilibrator
(Eckert & Ziegler Cesio) to a PAM aerosol monitor (TSI 3375) for continuous particle
size and concentration measurement. The operation and precision of the aerosol monitor
was validated using Aerodynamic particle sizer (APS) TSI 3321. The validation was
performed prior to the experiment using the identical setup, apart from adding of the
radioactive substance into the atomizer. The size of the particles generated by CMAG was
adjusted according to APS and the size displayed by the aerosol monitor was recorded.
Subsequently, during the experiment, only the radioactive substance was added to the
atomizer and particles of the same size were produced. Usually only a small correction of
saturator flow was needed at the beginning of the experiments. Only one size of particles
was measured during one day; therefore no further adjustments of the generator were
needed. The aerosol monitor served as an on-line indication of the process of aerosol
generation being stable. It was easily accomplished, as the exposure of the model to
the aerosol lasted only 5 to 15 mins. Filters consisting of Millipore AP40 glass fibers
were attached to the output branches of the respiratory tract model. The entire system
(Figure 6)
was enclosed in a plastic bag which was kept in a vacuum to prevent the active
aerosol from leaking into the laboratory. All the 10 terminal branches with flow meters
for flow rate control were combined into one branch with a protective High efficiency
particulate air (HEPA) filter. The vacuum was generated with a Busch R, 5 PA 0008 C
rotary oil vacuum pump. The flow distribution in each section of the model is provided
in Table 3.
Figure 6: A photograph of the physical model prior to the PET measurement. |
The whole exposure of the model to the aerosol was performed in a shielded laboratory
with an underpressure ventilating system, which would prevent the aerosol from escaping
the room in case of a primary safety system failure. The laboratory personnel were not
present in the laboratory during the exposure, with the exception of the regular instrument
supervision. Whenever they had to enter the lab, they wore half mask respirators. The
experiments were performed in a steady-state inhalation mode with the flow rates of 15,
30, and 60 L/min. Liquid monodisperse particles with mass median aerodynamic diameter
of 2.5 and 4.3 were used.
The standard geometric deviation of size was less than 1.24
for all measured regimes. The models were exposed for 10 to 15 minutes depending on
radioactivity decrease by radionuclide decay. The peak activity in the models was 4 to 60
(depending on the particle size and concentration, measuring mode, and model
used), as measured with an RP-2000 portable contamination meter (VF Zilina, CZ).
Deposited activity evaluation method
The model was transported to a Siemens Biograph 64 Truepoint PET-CT scanner immediately after the radioactive exposure. The transportation took approximately 3 minutes. The PET-CT scanner acquired firstly CT images, which were promptly followed by PET images.
Table 3: Flowrates through the individual sections of the model (section identification numbers are presented in Figure 4). Sections 13 – 22 lead to the 10 outlets. Flowrates in the intermediated sections (4 – 12) are calculated using mass continuity. |
Both the CT and PET images could be attributed unambiguously to the given geometry, owing to the assigned geometrical coordinates. The CT images are important for a precise localization of the edges of the sections, whereas the PET images contain the essential information about aerosol deposition. The CT and PET images were imported into Carimas 2.47 ver. 2.2.44.7002 SW (Turku PET Centre 2012) with the CT images as the main images and the PET images as background images. The software Carimas is a medical image processing tool, which was developed primarily for analysis of PET images at Turku PET Centre in Finland and is available as a freeware for non-commercial use from: http://www.turkupetcentre.fi/carimas/download/. The Carimas was developed in an “interactive data language” (IDL) and runs on IDL Virtual Machine™, a cross-platform utility for running IDL code. It supports multiple input formats (DICOM, ECAT, Analyse, Interfile, Nifti, Interfile, MicroPET) and general bitmap formats (JPG, TIFF, PNG and BMP). Researchers can perform visualization, segmentation, statistical analysis, or modelling of PET data. It is possible to perform the image fusion, i.e. to coregistrate PET and CT or MR, images. Users can define a volume of interest (VOI) and, by using the static image analysis, they can calculate the mean activity in the VOI (in units of the PET image), standard deviation, minimal and maximal values, the number of voxels analysed, and the volume of the region. The software was tested and compared to commercial tools with an excellent agreement as documented by Nesterov et al. (2009) and Harms et al. (2014). To facilitate the identification of the hot spots in the Carimas software, the BGRY colour system was used to display the PET images and the system range was reduced so that the hot spots were clearly seen in the model sections (Figure 7). The sections of the model were marked as volumes of interest (VOI), and the mean volume activity (in Bq/mL) was evaluated using the Carimas software. Each section was enveloped in an independent VOI. A cylinder or a sphere was selected as the starting shape of the VOI, depending on the shape of the sections. Firstly, the starting shape of the VOI was positioned and shaped as a whole to attain a suitable orientation and an adequate size. In the next step, each VOI was shaped in the vertex mode for a local mesh adjustment, so that the VOI surrounded the section with an overreach, while preventing an overlap of different VOIs. If the latter was not avoided, the observed activity would be attributed to both VOIs, resulting in biased results. The degree of mesh density can be modified by reducing the number of nodes or by enlarging or shrinking the mesh.
Figure 7: CT and PET image presentation in the Carimas software. |
With regard to the PET precision, a VOI should optimally overreach the section by
roughly 5mm on all sides to cover all the radiation emitted from the section. However, due
to the complex lung geometry, it was impossible to form all the VOIs with such a large
overreach, and so the uncounted radioactivity was accounted for in a correction VOI. Two
additional VOIs had to be created to determine the magnitude of the correction: a VOI for
the top part of the box, accommodating the model with the hoses, and a correction VOI
for the bottom part of the box, accommodating the filters downstream of the terminal
branches. They encompass the large top space with the model and the bottom space
between the two horizontal partitions respectively in Figure 6.
Two separate corrections
were needed for the model and filters, because the difference in volume radioactivity is
within few orders of magnitude and therefore it is necessary to assign the uncounted
radioactivity to its real source. It is essential to preserve the true ratio of the radiation
deposited in model sections and on filters. The correction factor for the sections was
calculated as the ratio of the total activity obtained from the VOI encompassing the
whole model to the activity obtained by summing up the section VOIs. Similarly, the
correction factor for the filters was calculated as the ratio of the total activity of the large
VOI encompassing all the filters, to the sum of the activities of the filter VOIs. The
corrections for the model sections
() and for the filters () were calculated as
follows:
where and are the activities measured in the model
sections and on the filters respectively, and and
are the activities measured in the correction VOI for the
model and filters respectively. We assume that the distribution of uncounted radioactivity
is proportional to the measured activity in the sections of the model. The higher the
activity measured in the filter, the higher the activity spread around the section.
Therefore, the uncounted activity is to be accounted for by the sections proportionally to their
measured activity. The correction formula then is simply:
The asterisk () indicates that the data has been corrected. The magnitude of the
correction factor can also be regarded as a model tightness indicator. If the model was not
tight enough and some aerosol leaked beyond the model, the amount of aerosol in the
correction VOI would be appreciably higher than in the sections; also, the calculated
correction factor would be higher.
The correction factors for the model sections and for the
filters lie largely within the ranges of 1.0 to 1.2 and 1.0 to 1.1 respectively. Aerosol leaks
should be suspected if these ranges were exceeded. The VOIs created can be saved in
VRML format in the Carimas software and imported to another model. Therefore, the
same set of VOIs was used for the analysis of all models. The corrected activities served to
calculate the deposition fractions (DF) for two aerosol particle sizes (DF is defined as the
ratio of inhaled particles deposited in the lungs section
to the total number of particles
entering the lungs). The total amount of aerosol was calculated as the sum of activities in
the model, in the output sections, in the connection hoses and on the filters. The observed
activity is assumed to be directly proportional to the number of deposited particles. This
is based on the fact that the radioactive atoms are dispersed uniformly in the constantly
stirred atomizer and that each particle generated contains a nucleus of the same size, on
which DEHS has condensed.
where is the total number of all sections including the hoses and filters.
Measurement errors
Sources of error in PET imaging
The quality of the images obtained with PET is affected by the following parameters: the spatial resolution, sensitivity, noise, scattered radiation and contrast. The parameters are interrelated, thus any efforts to reduce the effect of one parameter often result in an increased effect of another parameter (Saha, 2010). The accuracy of determining of the deposition point for a specific particle is related to the attainable spatial resolution of the PET scanner. The resolution is affected by the detector size, positron range, noncolinearity, reconstruction method and localization of the detectors: The actual spatial resolution of a particular PET scanner is measured by the standardized method recommended by the NEMA – National Electrical Manufacturers Association (2001). The full width at half maximum (FWHM) is the quantity serving as the measure of resolution. This quantity can be understood as the magnitude of point source “blurring” in the photograph or as the shortest distance at which two point sources can be discerned as such. The spatial resolutions (based on NEMA (2001)) for the PET Siemens Biograph 64 system are listed in Table 4.
Table 4: Spatial resolution of the PET scanner used for experiments (Siemens 2014) |
The remaining parameters of our PET scanner were: the sensitivity 4.2 cps/kBq, noise equivalent count rate 96 kcps at 35 kBq/cc, and scatter fraction < 38% (Lizal et al., 2015). As noted in the model limitations section the possible electrical charge on the walls of the physical model should be considered. The amount of depositing particles due to the mirror charge is probably low, however, it should be confirmed by an additional experiment.
Inter-person variability
To assess the error connected to the analysis of images, in particular to the definition of VOIs, the whole process, beginning with forming the VOIs, was performed by two technicians independently and their results were compared. The average of the differences between corresponding sections was 3% of the average value of the calculated radioactivity. More experiments should be performed to calculate the actual uncertainty of the whole method, but the above mentioned test confirmed that even if the VOIs were formed independently by two different individuals, their results were in close agreement (Lizal et al., 2015).
Validation of the method
The experimental setup and deposition measuring technique was validated by optical microscopy. Nickel coated particles (M-18-Ni, Cospheric, Santa Barbara, CA, USA) with MMAD 14.1μm were dispersed by the Small Scale Powder Disperser TSI 3433 and deposited in the model. The model was disassembled and particles were rinsed out of each section using ultrasonic bath. The particles were then filtered out and counted using optical microscopy. The results were consistent with the results acquired by PET. However, the presented PET method is capable of achieving a significantly higher spatial resolution. Therefore, more detailed data on local deposition can be obtained, as the size of VOIs can be smaller than the size of the sections of the model and is limited only by the spatial resolution of the PET scanner, in our case 5 mm (Lizal et al., 2015).
Measured data
The proposed method was used to evaluate deposition during steady-state inhalation for the following cases:
- 15 L/min,
- 30 L/min,
- 60 L/min,
The measurement for 60 L/min was made in duplicate. The Deposition Fraction difference between the two measurements was calculated for each section, and the observed difference within the repeated measurements was drawn in the plot in the form of error lines as the estimate of measurement variability.
The plot in fig. 8 demonstrates that the majority of particles was deposited in the oral cavity and in the trachea in all the applied variants. A comparison of the regimes shows that the amount of deposited particles increases with their increasing momentum, which confirms that inertial impaction was the dominant deposition mechanism in our case. The results demonstrate that the in vitro deposition measurement method based on the PET technique enables the DF to be evaluated with adequate precision for distinguishing the effects of the particle size and the breathing regime.
Figure 8: Results obtained by measurements in the realistic model (section identification numbers are presented in Figure 4). Error bars indicate the difference between two measurements performed on two separate, but geometrically identical models. Measurements for 15 and 30 L/min were performed just once. |
In the following step, a comparison with previously published data was performed.
Few studies presenting in vitro measured local aerosol deposition in tracheobronchial
airways have been published in past years. The most suitable for comparison with our data
were the results acquired by
Zhou & Cheng (2005) and
Chan & Lippmann (1980) as they
both used realistic lung casts. However, their models had less generations of branching
(four and six, respectively) compared to our model. The comparison of deposition efficiency (DE)
defined as the ratio of the number of particles deposited in the section to
the number of particles entering the section showed close agreement with both studies
mentioned above (see Figure 9).
The practical calculation of DE was performed as a ratio
of corrected activity in the examined section to the sum of activities in all downstream
sections including filters and hoses, and the activity in the examined section:
where is the corrected activity in the section downstream the i-th section. DE was
plotted as a function of Stokes number defined as:
where is the aerodynamic diameter of a particle,
is the velocity,
is the Viscosity of air and is the characteristic diameter (of the airway), which was
measured in the digital geometry in Rhinoceros 4.0 software (McNeel Seattle, WA7 USA).
The diameter of the entrance airway was used as the characteristic diameter, in accordance
with Zhou & Cheng (2005).
Chan & Lippmann (1980) introduced the following correlation
fit based on their experimental data for the first six generations:
Similar correlation fit was calculated based on the current data for 1st to 4th generation.
Both the empirical equations were plotted in the graph of Figure 9.
The airways in the
4th to 7th generation, which create separate parts of the model (identification numbers 13
to 22), were evaluated as a whole section and plotted in
Figure 9 as well. The deposition
fraction can be calculated in a higher resolution (creating smaller sections) down to the
resolution of the PET scanner (5Â mm). This resolution is significantly better than
the resolution of other commonly used methods (e.g. gravimetry or fluorometry), but
the evaluation of deposition in single airways with diameters smaller than PET scanner
resolution is not possible.
The measured Deposition data reported in figure 8 can be found in file DF_in-vitro.xlsx.
Figure 9: The comparison of the current results with previously published experimental data (symbols) and correlation fits (lines). The filled triangles represent the current data for single bifurcations in the 1st to 4th generation. The filled circles represent the current data for the model sections no. 13 to 22 (Figure 4), which contain multiple bifurcations. |
Contributed by: P. Koullapisa, F. Lizalb, J. Jedelskyb, L. Nicolaouc, K. Bauerd, O. Sgrotte, M. Jichab, M. Sommerfelde, S. C. Kassinosa —
aDepartment of Mechanical and Manufacturing Engineering, University of Cyprus, Nicosia, Cyprus
bFaculty of Mechanical Engineering, Brno University of Technology, Brno, Czech Republic
cDivision of Pulmonary and Critical Care, School of Medicine, Johns Hopkins University, Baltimore, USA
dInstitute of Mechanics and Fluid Dynamics, TU Bergakademie Freiberg, Freiberg, Germany
eInstitute Process Engineering, Otto von Guericke University, Halle (Saale), Germany
© copyright ERCOFTAC 2019