Published February 5, 2016 | Version 10004053

Modeling Aerosol Formation in an Electrically Heated Tobacco Product

Description

Philip Morris International (PMI) is developing a range
of novel tobacco products with the potential to reduce individual
risk and population harm in comparison to smoking cigarettes.
One of these products is the Tobacco Heating System 2.2 (THS
2.2), (named as the Electrically Heated Tobacco System (EHTS) in
this paper), already commercialized in a number of countries (e.g.,
Japan, Italy, Switzerland, Russia, Portugal and Romania). During use,
the patented EHTS heats a specifically designed tobacco product
(Electrically Heated Tobacco Product (EHTP)) when inserted into
a Holder (heating device). The EHTP contains tobacco material in
the form of a porous plug that undergoes a controlled heating process
to release chemical compounds into vapors, from which an aerosol
is formed during cooling. The aim of this work was to investigate
the aerosol formation characteristics for realistic operating conditions
of the EHTS as well as for relevant gas mixture compositions
measured in the EHTP aerosol consisting mostly of water, glycerol
and nicotine, but also other compounds at much lower concentrations.
The nucleation process taking place in the EHTP during use when
operated in the Holder has therefore been modeled numerically using
an extended Classical Nucleation Theory (CNT) for multicomponent
gas mixtures. Results from the performed simulations demonstrate
that aerosol droplets are formed only in the presence of an aerosol
former being mainly glycerol. Minor compounds in the gas mixture
were not able to reach a supersaturated state alone and therefore
could not generate aerosol droplets from the multicomponent gas
mixture at the operating conditions simulated. For the analytically
characterized aerosol composition and estimated operating conditions
of the EHTS and EHTP, glycerol was shown to be the main aerosol
former triggering the nucleation process in the EHTP. This implies
that according to the CNT, an aerosol former, such as glycerol
needs to be present in the gas mixture for an aerosol to form
under the tested operating conditions. To assess if these conclusions
are sensitive to the initial amount of the minor compounds and to
include and represent the total mass of the aerosol collected during
the analytical aerosol characterization, simulations were carried out
with initial masses of the minor compounds increased by as much
as a factor of 500. Despite this extreme condition, no aerosol
droplets were generated when glycerol, nicotine and water were
treated as inert species and therefore not actively contributing to the
nucleation process. This implies that according to the CNT, an aerosol
cannot be generated without the help of an aerosol former, from
the multicomponent gas mixtures at the compositions and operating
conditions estimated for the EHTP, even if all minor compounds are
released or generated in a single puff.

Files

10004053.pdf

Files (292.4 kB)

Name Size
md5:3bf1998d86feefcebbd30fa1d9eecb8e
292.4 kB Preview Download

Additional details

References

  • H. Vehkam¨aki, Classical nucleation theory in multicomponent systems. Springer-Verlag Berlin Heidelberg, 2006.
  • W. Hinds, Aerosol Technology. New York: Wiley-Interscience, 1999.
  • H. Nguyen, K. Okuyama, T. Mimura, Y. Kousaka, R. Flagan, and J. Seinfeld, "Homogeneous and heterogeneous nucleation in a laminar flow aerosol generator," Journal of Colloid And Interface Science, vol. 119, no. 2, pp. 491–504, 1987.
  • M. Wilck, K. H¨ameri, F. Stratmann, and M. Kulmala, "Experimental and theoretical examination of homogeneous nucleation in a laminar flow reactor (ucpc)," Journal of Aerosol Science, vol. 27, no. SUPPL.1, pp. S587–S588, 1996.
  • M. P. Anisimov, E. G. Fominykh, S. V. Akimov, and P. K. Hopke, "Vapor-gas/liquid nucleation experiments: A review of the challenges," Journal of Aerosol Science, vol. 40, no. 9, pp. 733 – 746, 2009.
  • C. Winkelmann, M. Nordlund, A. Kuczaj, S. Stolz, and B. Geurts, "Efficient second-order time-integration for single-species aerosol formation and evolution," International Journal for Numerical Methods in Fluids, vol. 74 (5), pp. 313–334, 2014.
  • E. Frederix, M. Stanic, A. Kuczaj, M. Nordlund, and B. Geurts, "Extension of the compressible PISO algorithm to single-species aerosol formation and transport," International Journal of Multiphase Flow, vol. 74, pp. 184–194, 2015.
  • H. Arstila, P. Korhonen, and M. Kulmala, "Ternary nucleation: Kinetics and application to water-ammonia-hydrochloric acid system," Journal of Aerosol Science, vol. 30, pp. 131–138, 1999.
  • G. Wilemski, "Composition of the critical nucleus in multicomponent vapor nucleation," The Journal of Chemical Physics, vol. 80, no. 3, pp. 1370–1372, 1984. [10] D. S. van Putten, "Efficient methods for N-component condensation," Ph.D. dissertation, University of Twente, Enschede, The Netherlands, 2011. [11] R. Becker and W. Dring, "Kinetische behandlung der keimbildung in bersttigten dmpfen," Annalen der Physik, vol. 416, no. 8, pp. 719–752, 1935. [12] "The OpenFOAM Foundation," http://www.openfoam.com, 2015. [13] H. Jasak, "Error analysis and estimation for the finite volume method with applications to fluid flows," Ph.D. dissertation, Imperial College, University of London, 1996. [14] M. Nordlund and A. Kuczaj, "Aerosol dosimetry modeling using computational fluid dynamics," in Computational Systems Toxicology, ser. Methods in Pharmacology and Toxicology, J. Hoeng and M. Peitsch, Eds. Humana Press: Springer, 2015. [15] R. Bird, W. Stewart, and E. Lightfoot, Transport Phenomena, revised second ed. New York: John Wiley & Sons Inc., 2007. [16] D. Kashchiev, Nucleation: Basic Theory with Applications. Butterworth-Heinemann: Elsevier Science, 2000. [17] G. Wilemski and B. E. Wyslouzil, "Binary nucleation kinetics. I. Self-consistent size distribution," The Journal of Chemical Physics, vol. 103, no. 3, pp. 1127–1136, 1995. [18] D. S. van Putten, S. P. Glazenborg, R. Hagmeijer, and C. H. Venner, "A multigrid method for N-component nucleation," The Journal of Chemical Physics, vol. 135, no. 1, 2011. [19] R. Issa, B. Ahmadi-Befrui, K. Beshay, and A. Gosman, "Solution of the implicitly discretised reacting flow equations by operator-splitting," Journal of Computational Physics, vol. 93, pp. 388–410, 1991. [20] R. Baker, "The development and significance of standards for smoking-machine methodology," Beitrge zur Tabakforschung / Contributions to Tobacco Research, vol. 20, no. 1, pp. 23–41, 2002. [21] "Philip Morris International," http://www.pmiscience.com, 2016. [22] F. Barontini, A. Tugnoli, V. Cozzani, J. Tetteh, M. Jarriault, and I. Zinovik, "Volatile products formed in the thermal decomposition of a tobacco substrate," Industrial & Engineering Chemistry Research, vol. 52, no. 42, pp. 14 984–14 997, 2013. [23] V.-G. V. und Chemieingenieurwesen and V. Gesellschaft, VDI Heat Atlas, ser. Springer reference. Springer, 2010. [24] "CAMEO Chemicals database," http://cameochemicals.noaa.gov, 2015. [25] T. Daubert, R. Danner, and N. Design Institute for Physical Property Data (New York, Data compilation tables of properties of pure compounds. American Institute of Chemical Engineers, 1985, no. Teil 2. [26] B. Poling, J. Prausnitz, and J. O'Connell, The properties of gases and liquids. McGraw-Hill, 2001. [27] R. Warming and R. Beam, "Upwind second-order difference schemes and applications in aerodynamic flows," AIAA Journal, vol. 14, pp. 1241–1249, 1976.