Integration of Safety Aspects in Modeling of Superheated Steam Flash Drying of Tobacco

: Knowledge of the drying properties of tobacco in high temperatures above 100 ◦ C and its dust are crucial in the design of dryers, both in the optimization of the superheated-steam-drying process and in the correct selection of innovative explosion protection and mitigation systems. In this study, tobacco properties were determined and incorporated into the proposed model of an expanding superheated steam ﬂash dryer. The results obtained from the proposed model were validated by using experimental data yielded during test runs of an industrial scale of a closed-loop expansion dryer on lamina cut tobacco. Moreover, the explosion and ﬁre properties of tobacco dust before and after the superheated steam-drying process at 160, 170, 180, and 190 ◦ C were experimentally investigated, using a 20 L spherical explosion chamber, a hot plate apparatus, a Hartmann tube apparatus, and a Godbert–Greenwald furnace apparatus. The results indicate that the higher the drying temperature, the more likely the ignition of the dust tobacco cloud, the faster the explosion ﬂame propagation, and the greater the explosion severity. Tobacco dust is of weak explosion class. Dust obtained by drying with superheated steam at 190 ◦ C is characterized by the highest value of explosion index amounting to 109 ± 14 m · bar · s − 1 , the highest explosion pressure rate (405 ± 32 bar/s), and the maximum explosion pressure (6.7 ± 0.3 bar). The prevention of tobacco-dust accumulation and its removal from the outer surfaces of machinery and equipment used in the superheated steam-drying process are highly desirable.


Introduction
Drying is a necessary process and the most energy-consuming step in tobacco production.The selection of drying method affects the tobacco's chemical and physical properties and the quality of cigarette smoke [1][2][3].Tobacco undergoes drying at least three times in its life before it leaves the cigarette factory.First, the tobacco is cured and dried just after harvesting.This stage is necessary to lower the water activity of green tissue below a lever, ensuring that all microbial activity in the product will be terminated and the product can be handled safely.Thus, the green leaf is dried from ca. 85% moisture content to 13 ÷ 18% in atmospheric air (air-curing) or in heated air (flue-curing).Curing that takes place simultaneously with drying is responsible for developing the necessary sensory characteristics of the product [4].Second, leaves from farms are processed by tobacco wholesalers.Tobacco leaves are conditioned to 20 ÷ 22% and threshed (leaf stem is separated from lamina).Tobacco is then re-dried down to 7%, followed by cooling and remoistening to 11%.In this form, tobacco is baled and ready for export.Third, baled tobacco is rewetted at a cigarette Energies 2021, 14, 5927 2 of 22 factory to restore the plant tissue's elasticity, so that the product can be handled without cracking and generating excessive dust.After resting awhile, the tobacco is shredded, and the cut tobacco is dried again to its final moisture content.In this drying process, tobacco can be simultaneously expanded to increase its filling power (and so reduce tissue density).Most often, an ensemble of expansion tunnel and rotary dryers or a one-stage superheated steam expansion dryer (closed-loop expansion dryer (CLED)) are used [4].After this stage, the tobacco is conditioned again by adding some additional components (called sauces), which improve flavor and taste.However, the tobacco will no longer be dried in a dryer, but will only lose some moisture (1-2%) during the cigarette-making process.This paper concerns modeling the drying process of tobacco that takes place in the expanding superheated steam dryer in order to predict the properties of tobacco in high temperatures above 100 • C and the efficiency of drying systems.
During cigarette production, tobacco dust from processed leaves is generated on equipment near the superheated-steam-drying facility and can result in fire and explosion hazards that can lead to human fatalities, serious injuries, and financial loss.Tobacco dust poses a risk of explosion due to the formation of dust-air mixtures from the sources of emissions and the formation of a dust layer, which may be ignited by self-heating [5].The ignited layer of dust may then play the role of an ignition source for subsequent explosions.Very often, designers of technological processes do not have sufficient knowledge concerning the basic physical and chemical parameters of materials [6] and their dusts occurring under specific process conditions.Data related to cut tobacco in high temperatures above 100 • C are scarce and not sufficiently available in the literature [7].Furthermore, there is no information about the influence of the superheated-steam-drying process on tobacco dust flammability and explosiveness parameters.
The present research is intended to develop a model of an expanding superheated steam flash dryer for cut tobacco incorporating determined tobacco properties (sorption isotherms, sorption isobars, particle morphology, drag coefficient, and heat capacity) and validate the results obtained from the proposed model, using experimental data obtained during test runs of an industrial scale CLED on lamina cut tobacco.Additional explosion and fire properties of dust tobacco, such as maximum explosion pressure, maximum rate of pressure rise, explosion index, minimum explosion concentration, minimum ignition energy, minimum ignition temperature, and limiting oxygen concentration, were determined that are critical for the risk assessment and safety management of dust explosions [8].Moreover, the influence of the superheated-steam-drying process on tobacco dust flammability and explosiveness parameters was studied.The results obtained are crucial in the design of dryers, including the optimization of the superheated-steam-drying process, the correct selection of innovative explosion protection, and a mitigation system in order to assure the safety of working people.

Materials
Virginia tobacco was used as the raw material in this study.Samples of tobacco were collected in 2019/2020 from the Poland Tobacco Industrial Corporation and sent to laboratories for physicochemical property analysis.The ultimate analysis of tobacco was performed by using a Vario EL Cube elemental analyzer (Elementar Analysensysteme GmbH, Langenselbold, Germany) and fractional analysis, using a thermogravimetric analyzer (TGA/SDTA851 LF, Mettler-Toledo, Greifensee, Switzerland).The results are presented in Table 1.

Properties of Tobacco
Knowledge of the drying properties of tobacco and its dust is critical in designing a superheated steam dryer and optimizing production costs [9,10].In the present work, tobacco properties were obtained and incorporated in the proposed model of an expanding Energies 2021, 14, 5927 3 of 22 superheated steam flash dryer.Furthermore, the explosion and fire properties of dust tobacco were determined by using a 20 L spherical explosion chamber, hot plate apparatus, Hartmann tube apparatus, and Godbert-Greenwald furnace apparatus.

Sorption Isotherms and Isobars
Water in plant tissue may exist in two forms: bound water and unbound water.Water is considered bound when its partial pressure is lower than the saturated pressure of pure water.The ratio of partial pressure of water over tobacco (p) to saturated water pressure (p s ) at the same temperature (approximately) is called water activity, a w , and is described in Equation (1).
The dependence of water activity on moisture content b (in %) or X (in kg of water/kg dry mass) at constant temperature is called a sorption isotherm.A sorption isotherm represents the thermodynamic equilibrium between the wet solid and humid air and is indispensable for calculating rehumidification, storage, drying, and other processes involving moist air and tobacco.The sorption isotherm can be obtained experimentally by several methods, of which conditioning tobacco samples' oversaturated solutions of various salts is the simplest and cheapest of all.By changing salt solutions, the partial pressure of water in the air is controlled at a constant temperature.The experimental data are obtained either by starting with a wet solid and gradually lowering the moisture content (desorption isotherms) or by measuring moisture uptake of a previously dried solid (sorption isotherm).The two may not be identical, and certain sorption hysteresis almost always exists.In general, sorption isotherms have the shape of a learning curve (S-shape) and can be described by one of several dozens of empirical or semi-empirical equations.Ono and Yoshida [11] obtained the following equations for sorption isotherms of tobacco.
For flue-cured tobacco: For air-cured tobacco: These equations are valid for temperatures up to 80 • C. For green Virginia-type tobacco, Kiranoudis et al. [12] recommended Equation (4).
Legros et al. [13] derived Equation ( 5) based on cut lamina obtained from the BAT company for blend 72 in the range of 20 to 60 • C.
Energies 2021, 14, 5927 The above equations are plotted in Figure 1, and they are compared with experimental results obtained by Miyauchi et al. [14] for bright yellow and Burley tobaccos.The above equations are plotted in Figure 1, and they are compared with experimental results obtained by Miyauchi et al. [14] for bright yellow and Burley tobaccos.In superheated steam-drying, there is no inert gas (air) present in the system, so water activity is defined as the ratio of actual water pressure (ambient pressure) to the saturated pressure of pure water at a given temperature.Since the measurements are made at constant ambient pressure, the curves obtained are called sorption isobars.By changing temperature, the saturated pressure of water is controlled at constant pressure.The measurements of sorption isobars require a more sophisticated setup, especially when operating pressure is different than atmospheric pressure.The sorption isobars of tobacco were investigated in a setup built at the Lodz University of Technology (Poland), as shown in Figure 2. A detailed description of the apparatus and test procedure was presented in our previous publication [15].In superheated steam-drying, there is no inert gas (air) present in the system, so water activity is defined as the ratio of actual water pressure (ambient pressure) to the saturated pressure of pure water at a given temperature.Since the measurements are made at constant ambient pressure, the curves obtained are called sorption isobars.By changing temperature, the saturated pressure of water is controlled at constant pressure.The measurements of sorption isobars require a more sophisticated setup, especially when operating pressure is different than atmospheric pressure.The sorption isobars of tobacco were investigated in a setup built at the Lodz University of Technology (Poland), as shown in Figure 2. A detailed description of the apparatus and test procedure was presented in our previous publication [15].Sorption isobars are traditionally presented to mean moisture content dependent on temperature.Exemplary sorption isobars for tobacco lamina and stem measured in this setup are shown in Figure 3.   Sorption isobars are traditionally presented to mean moisture content dependent on temperature.Exemplary sorption isobars for tobacco lamina and stem measured in this setup are shown in Figure 3. Again, sorption isobars may be approximated by a selected empirical equation.However, one may consider replotting Figure 3 as aw = f(X), defining aw by the following equation.The curve obtained is shown in Figure 4 for lamina only.Again, sorption isobars may be approximated by a selected empirical equation.However, one may consider replotting Figure 3 as a w = f(X), defining a w by the following equation.
The curve obtained is shown in Figure 4 for lamina only.The curve that can be fitted to the points of Figure 4 is not an isotherm, since the points are obtained at varying temperatures, but an isobar.Adding a third axis, i.e., temperature, would allow the data of sorption isotherms and sorption isobars to be shown simultaneously.Such a 3D chart can be thought of as a generalization of sorptional equilibrium over a wide temperature and moisture-content range [16] that can be represented as a 3D surface that can be fitted with a suitable empirical equation.Such an empirical equation by Chen and Clayton [17] was fitted to scarce data available for tobacco, and the result is shown in Figure 5.The curve that can be fitted to the points of Figure 4 is not an isotherm, since the points are obtained at varying temperatures, but an isobar.Adding a third axis, i.e., temperature, would allow the data of sorption isotherms and sorption isobars to be shown simultaneously.Such a 3D chart can be thought of as a generalization of sorptional equilibrium over a wide temperature and moisture-content range [16] that can be represented as a 3D surface that can be fitted with a suitable empirical equation.Such an empirical equation by Chen and Clayton [17] was fitted to scarce data available for tobacco, and the result is shown in Figure 5.The curve that can be fitted to the points of Figure 4 is not an isotherm, since the points are obtained at varying temperatures, but an isobar.Adding a third axis, i.e., temperature, would allow the data of sorption isotherms and sorption isobars to be shown simultaneously.Such a 3D chart can be thought of as a generalization of sorptional equilibrium over a wide temperature and moisture-content range [16] that can be represented as a 3D surface that can be fitted with a suitable empirical equation.Such an empirical equation by Chen and Clayton [17] was fitted to scarce data available for tobacco, and the result is shown in Figure 5.The equation has the following form: where a 1 = −10.593,a 2 = −0.54254, a 3 = −2.9855,and a 4 = 0.3993.It is obvious that producing a more reliable equation would require a substantial experimental effort.It may be added that having in mind the three types of drying processes mentioned in the introduction, each one would require its sorption equilibrium to be measured individually.It is well-known that green tissue has different sorptional properties than rehydrated tissue.Similarly, adding sauces that contain soluble components can significantly alter the sorptional properties making tobacco more hygroscopic.Differences between tobacco cultivars in terms of sorptional equilibrium are less pronounced than differences between lamina and stem.Apparently, it is the cell morphology that decides water activity.Further studies to support this thesis are needed.

•
Heat capacity of tobacco: The heat capacity of tobacco is necessary for the calculation of solid temperature in the falling drying rate period.It can be measured in a scanning calorimeter as a function of temperature.The specific heat of dry tobacco was measured by a scanning calorimeter in a range of 25-140 • C by Pakowski et al. [4] and an average value of 1.467 kJ/kgK was adopted.Kerekes et al. [18] quotes a slightly smaller value of 1.382 kJ/kg K obtained during air-curing of a Hungarian tobacco variety.

•
Particle morphology and drag coefficient: Cut tobacco has the form of thin ribbons, which are extremely difficult to describe.In pneumatic transport, the force of resistance to air flow of a single particle can be described in Equation (8).
where C D is the drag coefficient, u r = u g − u m is the velocity difference of solid and air, and S p is the cross-section area perpendicular to flow.Formulas for the drag coefficient depend on particle morphology.For isometric particles (those having all three basic dimensions of the same order of magnitude), the following formula of Haider and Levenspiel [19] is recommended. where ψ = sur f ace area o f equivalent sphere particle area In a single-particle model of pneumatic transport, it is necessary to represent the solid phase by individual monodisperse particles.Despite inherent difficulties in describing the morphology of cut tobacco [20], we found out that for the purpose of our model the combination of equivalent particle diameter, d p , of 1.2 mm and sphericity ψ of 0.15 together with a particle density of 800 kg/m 3 works well and gives good simulation results.

•
Explosion and fire properties of dust tobacco: Explosion and fire properties of tobacco after drying with superheated steam, such as maximum explosion pressure, P max ; maximum rate of pressure rise, (dP/dt) max ; explosion index, K St ; minimum explosion concentration, MEC; minimum ignition energy, MIE; minimum ignition temperature, MIT; and limiting oxygen concentration, LOC [6,7], were experimentally determined according to standards shown in Table 2.
Samples of tobacco were milled with a vibration laboratory mill (LMW, Testchem, Pszów, Poland), and powder with a grain size lower than 500 µm was chosen for the experiments.An IPS particle size analyzer (Kamika Instruments, Warsaw, Poland) was used to evaluate the particle size distribution.The Sauter mean diameters of the tobacco dust particles were 86.7 µm.Such sizes of tobacco dust can deposit on device surfaces and pose potential fire and explosion hazards.
A standard 20 L dust-explosion vessel (SPD-2.2,ANKO, Warsaw, Poland) with a dispersion, vacuum, pressure sensor, ignition, control, and data-acquisition system, as shown in Figure 6, was applied to estimate P max , (dP/dt) max , K St , LOC, and MEC.The tobacco dust was loaded in a 0.6 L reservoir, premixed with pressurized air (20 bar), and then dispersed into the vessel through the fast-acting valve and ignited after 60 ms.Two pyrotechnical igniters (Fr.Sobbe GmbH, Dortmund, Germany) with ignition energy each of 1 kJ were used to measure MEC and LOC, and of 5 kJ to evaluate P max and (dP/dt) max .Tests in the 20 L dust-explosion vessel were conducted by systematically increasing the dust concentration from 125 to 1500 g/m 3 , until P ex and (dP/dt) ex results reached a maximum value (P max , (dP/dt) max ).In the case of MEC, three measurement series were made, starting with a 500 g/m 3 concentration, followed by 250, 125, and 60 g/m 3 .dispersion, vacuum, pressure sensor, ignition, control, and data-acquisition system, as shown in Figure 6, was applied to estimate Pmax, (dP/dt)max, KSt, LOC, and MEC.The tobacco dust was loaded in a 0.6 L reservoir, premixed with pressurized air (20 bar), and then dispersed into the vessel through the fast-acting valve and ignited after 60 ms.Two pyrotechnical igniters (Fr.Sobbe GmbH, Dortmund, Germany) with ignition energy each of 1 kJ were used to measure MEC and LOC, and of 5 kJ to evaluate Pmax and (dP/dt)max.Tests in the 20 L dust-explosion vessel were conducted by systematically increasing the dust concentration from 125 to 1500 g/m 3 , until Pex and (dP/dt)ex results reached a maximum value (Pmax, (dP/dt)max).In the case of MEC, three measurement series were made, starting with a 500 g/m 3 concentration, followed by 250, 125, and 60 g/m 3 .The minimum ignition temperature of the tobacco dust cloud was determined by using the equipment presented in Figure 7.The central part of the Godbert-Greenwald furnace (MAIT/MIT, Anko, Warsaw, Poland) was an open-bottom silica tube with a diameter of 3.9 cm and a height of 23 cm.This tube was heated to the desired temperature with a device that controls the temperature of the furnace [21].The tests confirmed that the temperature distribution of the furnace tube was homogeneous.During each experiment, compressed air was introduced into the tube at pressures of 10, 20, 30, and 50 kPa to create a dust cloud.Tobacco dust loads in measurements of 0.1, 0.2, and 0.3 g, and 0.5 g were used.The first measurement was made at 500 • C, a weight of 0.3 g, and pressure of blown air of 30 kPa.The dust mass of tobacco at which the highest flame intensity was observed was determined.The minimum ignition temperature of the tobacco dust cloud was determined by using the equipment presented in Figure 7.The central part of the Godbert-Greenwald furnace (MAIT/MIT, Anko, Warsaw, Poland) was an open-bottom silica tube with a diameter of 3.9 cm and a height of 23 cm.This tube was heated to the desired temperature with a device that controls the temperature of the furnace [21].The tests confirmed that the temperature distribution of the furnace tube was homogeneous.During each experiment, compressed air was introduced into the tube at pressures of 10, 20, 30, and 50 kPa to create a dust cloud.Tobacco dust loads in measurements of 0.1, 0.2, and 0.3 g, and 0.5 g were used.The first measurement was made at 500 °C, a weight of 0.3 g, and pressure of blown air of 30 kPa.The dust mass of tobacco at which the highest flame intensity was observed was determined.The Hartmann tube equipment (MINOR 2, ANKO, Warsaw, Poland) depicted in Figure 8 was applied to determine the minimum ignition energy of the tobacco dust-air mixture.The equipment was composed of a spark-generating circuit (1 mH and 1μH), a glass tube with a pneumatic dust dispersion system, and a control panel that allowed the The Hartmann tube equipment (MINOR 2, ANKO, Warsaw, Poland) depicted in Figure 8 was applied to determine the minimum ignition energy of the tobacco dust-air mixture.The equipment was composed of a spark-generating circuit (1 mH and 1 µH), a glass tube with a pneumatic dust dispersion system, and a control panel that allowed the delay time to be set from 0 to 10,000 ms (±10 ms) and regulated the discharge spark energy from 1 to 1000 mJ [22].In MIE, eleven measurement series were performed, starting with a 125 g/m 3 concentration and ending with 3000 g/m 3 .
Energies 2021, 14, x FOR PEER REVIEW 10 of 23 The measuring station of the hot plate (LIT-3, ANKO TRADING, Warsaw, Poland) presented in Figure 9 was used to determine the minimum ignition temperatures at which a layer of tobacco dust with the thickness of 5 and 50 mm deposited on a hot plate at a constant temperature and ignited and/or thermally decomposed [23].To indicate the minimum ignition temperatures of the tobacco dust layer, two methods were used: the constant temperature of the heated plate (CTHP) and the constant rate of heat generation of 3 °C/min (CRHG).

Principles of Superheated Steam-Drying
When the steam is superheated, it can take in moisture from the wet solid until it becomes saturated.One part of the exiting steam can be superheated again and recycled into the dryer, while the other part is bled from the dryer and its heat can be recovered by condensation.This is, in short, a principle of superheated steam-drying.The idea was developed in Germany and was first mentioned in the literature in 1898 [20].Although technically more difficult to accomplish than air-drying, superheated steam-drying has several advantages, which are shown below.

•
Heat used for drying can be recovered thus lowering the net heat consumption even below the heat of vaporization.The measuring station of the hot plate (LIT-3, ANKO TRADING, Warsaw, Poland) presented in Figure 9 was used to determine the minimum ignition temperatures at which a layer of tobacco dust with the thickness of 5 and 50 mm deposited on a hot plate at a constant temperature and ignited and/or thermally decomposed [23].To indicate the minimum ignition temperatures of the tobacco dust layer, two methods were used: the constant temperature of the heated plate (CTHP) and the constant rate of heat generation of 3 • C/min (CRHG).The measuring station of the hot plate (LIT-3, ANKO TRADING, Warsaw, Poland) presented in Figure 9 was used to determine the minimum ignition temperatures at which a layer of tobacco dust with the thickness of 5 and 50 mm deposited on a hot plate at a constant temperature and ignited and/or thermally decomposed [23].To indicate the minimum ignition temperatures of the tobacco dust layer, two methods were used: the constant temperature of the heated plate (CTHP) and the constant rate of heat generation of 3 °C/min (CRHG).

Principles of Superheated Steam-Drying
When the steam is superheated, it can take in moisture from the wet solid until it becomes saturated.One part of the exiting steam can be superheated again and recycled into the dryer, while the other part is bled from the dryer and its heat can be recovered by condensation.This is, in short, a principle of superheated steam-drying.The idea was developed in Germany and was first mentioned in the literature in 1898 [20].Although technically more difficult to accomplish than air-drying, superheated steam-drying has several advantages, which are shown below.

•
Heat used for drying can be recovered thus lowering the net heat consumption even

Principles of Superheated Steam-Drying
When the steam is superheated, it can take in moisture from the wet solid until it becomes saturated.One part of the exiting steam can be superheated again and recycled into the dryer, while the other part is bled from the dryer and its heat can be recovered by condensation.This is, in short, a principle of superheated steam-drying.The idea was developed in Germany and was first mentioned in the literature in 1898 [20].Although technically more difficult to accomplish than air-drying, superheated steam-drying has several advantages, which are shown below.

•
Heat used for drying can be recovered thus lowering the net heat consumption even below the heat of vaporization.

•
Solid temperature exceeds boiling point temperature: This effect causes the expansion of tobacco and sterilizes the product biologically.

•
The atmosphere is virtually oxygen-free and thus oxidation reactions and danger of ignition are eliminated.

•
Steam-drying is more effective than air-drying, provided that steam temperature is higher than the so-called inversion temperature (at atmospheric pressure the inversion temperature is ca.175 • C).
The main disadvantages of superheated steam-drying are the following: materials undergoing drying must be temperature resistant, superheated steam is corrosive and requires better construction materials than air-drying, and thermal insulation must be good enough to prevent condensation on the walls [4,7].
In convective air-drying, when evaporating free moisture, heat must be transferred through the boundary layer to the material and the resulting vapor must diffuse back through the boundary layer into the core of the air flow.Both processes limit the rate of drying.In superheated steam-drying, there is no diffusion of vapor through the gas phase, because the gas phase is pure vapor.Therefore, the drying rate of free moisture is entirely controlled by heat transfer.
When a wet solid is immersed into superheated steam, the following three periods can be observed.
Condensation period: When the solid temperature is lower than the saturation temperature, steam condenses on the solid, and its temperature and moisture content increase.The drying rate in this period is negative and can be described as follows: Constant drying rate period: After reaching the saturation temperature, free moisture is evaporated, and all heat transferred to the solid is used for moisture evaporation.This holds the solid temperature and the drying rate constant.The drying rate can be described by the following equation.
Falling drying rate period: When all free moisture is evaporated, the drying rate decreases.Now moisture must diffuse to the surface from the depths of the solid phase, and this slows the process down.In effect, the net heat transferred to the solid is only partly used for evaporation of moisture, and thus the temperature of the solid increases.In order to describe the drying rate in this period, a detailed model of water diffusion in the solid would be needed.Such models are complex and difficult to use in practice.Instead, a simple empirical model can be used where the drying rate is proportional to the portion of moisture content remaining in the solid according to the following equation.
In this equation, n is an empirical parameter that can be evaluated on the basis of experimental data.If no information is available, n = 1 is a good initial guess.Both critical moisture content, X c , and equilibrium moisture content, X eq , can be obtained from the sorption isobar.
A superheated steam flash dryer (CLED) is shown schematically in Figure 10.The drying pipe is the section between control points 0 and 1. Wet tobacco is introduced through a rotating airlock into the pipe when it is carried away by steam flowing with a velocity of ca.27-30 m/s.During its short residence in the pipe, the tobacco loses moisture and heats up.It is then separated from the steam in the cyclone and leaves the dryer through another airlock.The steam leaving the cyclone is partly vented and partly circulated by a circulating fan.It is then indirectly reheated in the superheater and returned to the dryer pipe.

Earlier Work
Models for superheated-steam flash drying available vary in complexity from relatively simple (e.g., Reference [24]) to complex ones (e.g., Reference [25]).Generally, the solid is represented by monodisperse or polydisperse isotropic particles characterized by certain equivalent diameters (diameter of a sphere of equal volume) and sphericity (Equation ( 11)).In the case of polydisperse solids [24], the particle size range needs to be divided into seven fractions each, calculated separately.In this case, small particles may be already dry, while large particles are still in the heat-up period: it is also recommended that the heat transfer coefficient be adjusted for simultaneous mass transfer, using the correction shown below: where However, for temperatures of steam up to 300 °C, at atmospheric pressure, this correction does not exceed 94% of the value for pure heat transfer, so it can be neglected in engineering calculations.A two-dimensional model of pneumatic transport used by Reference [25] is recommended for more accurate calculations: the degree of complexity associated with it makes the calculations too expensive for ordinary applications.Furthermore, the improved accuracy it brings may be offset by the inaccuracy of empirical coefficients used in the model. Assumptions: In flash-drying, dispersed tobacco particles are fed into superheated steam, flowing at a velocity of ca. 30 m/s, which is much higher than the entrainment velocity.The particles are carried away by steam in the dryer pipe when they undergo drying and expansion.Dry particles are separated in a cyclone.The process is complex, and its math-

Proposed Model of Superheated Steam Flash Drying of Cut Tobacco 4.1. Earlier Work
Models for superheated-steam flash drying available vary in complexity from relatively simple (e.g., Reference [24]) to complex ones (e.g., Reference [25]).Generally, the solid is represented by monodisperse or polydisperse isotropic particles characterized by certain equivalent diameters (diameter of a sphere of equal volume) and sphericity (Equation ( 11)).In the case of polydisperse solids [24], the particle size range needs to be divided into seven fractions each, calculated separately.In this case, small particles may be already dry, while large particles are still in the heat-up period: it is also recommended that the heat transfer coefficient be adjusted for simultaneous mass transfer, using the correction shown below: where However, for temperatures of steam up to 300 • C, at atmospheric pressure, this correction does not exceed 94% of the value for pure heat transfer, so it can be neglected in engineering calculations.A two-dimensional model of pneumatic transport used by Reference [25] is recommended for more accurate calculations: the degree of complexity associated with it makes the calculations too expensive for ordinary applications.Furthermore, the improved accuracy it brings may be offset by the inaccuracy of empirical coefficients used in the model.

Assumptions:
In flash-drying, dispersed tobacco particles are fed into superheated steam, flowing at a velocity of ca. 30 m/s, which is much higher than the entrainment velocity.The particles are carried away by steam in the dryer pipe when they undergo drying and expansion.Dry particles are separated in a cyclone.The process is complex, and its mathematical description must be considerably simplified in order to construct a model that can be solved.The following assumptions are made:

•
The solid is represented by single isometric particles of the same size;

•
The flow of the suspension is one-dimensional.That is to say, no radial distribution of velocity, temperature, and/or concentration in the pipe are allowed; • A steady-state plant operation is assumed; • Only the process of drying in the pipe is considered: drying in a cyclone would require another model.

Model Formulation
The steady-state one dimensional distributed parameter model of the flash dryer can be constructed by establishing the momentum, heat, and mass balances over a differential volume element of the pipe shown in Figure 11.The variables used in Section 4.2 are defined in the nomenclature section.
Energies 2021, 14, x FOR PEER REVIEW 13 of 23 ematical description must be considerably simplified in order to construct a model that can be solved.The following assumptions are made: The solid is represented by single isometric particles of the same size; The flow of the suspension is one-dimensional.That is to say, no radial distribution of velocity, temperature, and/or concentration in the pipe are allowed; • A steady-state plant operation is assumed; • Only the process of drying in the pipe is considered: drying in a cyclone would require another model.

Model Formulation
The steady-state one dimensional distributed parameter model of the flash dryer can be constructed by establishing the momentum, heat, and mass balances over a differential volume element of the pipe shown in Figure 11.The variables used in Section 4.2 are defined in the nomenclature section.The following balances can be written.

•
Mass balance for the solid phase: from which we have the following: • Enthalpy balance for the solid phase: where, by using enthalpy definitions, we have the following: we have the following: Then, by using ( 20)-( 21) and rearranging, we finally have the following: The following balances can be written.

•
Mass balance for the solid phase: from which we have the following: • Enthalpy balance for the solid phase: where, by using enthalpy definitions, we have the following: we have the following: Then, by using ( 20)-( 21) and rearranging, we finally have the following: • Mass balance for the gas phase: since Then, by combining ( 24) and ( 25), we have the following: and Introducing the ideal gas equation for calculation of gas density change with temperature, we finally have the following: • Enthalpy balance for gas phase: Since the gas phase is pure steam, by introducing into (23) and rearranging, we have the following: Now, to be able to calculate changes of solid velocity and voidage with height, two additional equations must be introduced.

•
The continuity equation for the solid phase: Then, after differentiation and separation of terms, we have the following: • Kinetic energy balance for the solid phase: Now, by differentiation, we have the following: Since the first term is zero, we have the following: The drag force can be calculated from numerous theories.One often used reads as follows: Finally, the kinetic energy balance takes the following form: The momentum balance equation for the gas phase can be used for the calculation of pressure drop along the dryer pipe: The above-framed relationships form a one-dimensional model of flash drying with superheated steam.For the calculation of several variables involved in this formulation, constitutive-and property-defining equations are necessary.The following relationships are used.

•
Sensible heat flux: The sensible heat flux, q, is calculated as follows: where α is the heat transfer coefficient.Numerous correlations exist for the calculation of heat transfer coefficients in pneumatic transport.For the purpose of this model, the following equations of Bandrowski and Kaczmarzyk [17] were selected.For u g − u m < 1.2u t , we have the following: For u g − u m ≥ 1.2u t , we have the following: where • Drying rate: Drying rate equations were defined earlier as Equations ( 7) through (9).

• Transport velocity:
The following formula of Goroshko et al. [19] was used.
Re t = Ar where • Voidage: • Friction factors: The gas-wall friction factor can be calculated from the following formula by Reference [26].
where Re D is the Reynolds number based on ductwork diameter.The particle-wall friction factor in pneumatic transport was studied by Reference [27], and the following formulas were obtained.
For the acceleration zone, we have the following: Outside the acceleration zone, we have the following: • Specific area of contact:

Model Simulation
The proposed model of superheated steam flash drying of cut tobacco is composed of seven ordinary differential equations (ODEs), 18, 23, 29, 32, 35, 40, and 41.Equations (42) through (53) provide necessary auxiliary variables.All necessary steam properties are calculated from suitable relationships [28,29].The set of ODEs can be solved, provided that initial conditions are available.Since the dryer operates in a closed loop, the external heat and mass balances for other components of the loop must be calculated until the parameters of steam at control point 5 (Figure 10) are identical to those of point 0. This requires several iterations.In our study, the model was solved by using Mathcad mathematical software, assuming the general capacity of tobacco drying in a flash dryer of 3 tons/hour.Typical results of the basic parameters of the tobacco drying in the function of the apparatus length are shown in Figures 12-14, while Table 3 gives the results of the simulation using a superheated steam flash pneumatic dryer.0. This requires several iterations.In our study, the model was solved by using Math mathematical software, assuming the general capacity of tobacco drying in a flash d of 3 tons/hour.Typical results of the basic parameters of the tobacco drying in the f tion of the apparatus length are shown in Figures 12-14, while Table 3 gives the resul the simulation using a superheated steam flash pneumatic dryer.mathematical software, assuming the general capacity of tobacco drying in a flash dr of 3 tons/hour.Typical results of the basic parameters of the tobacco drying in the fu tion of the apparatus length are shown in Figures 12-14, while Table 3 gives the result the simulation using a superheated steam flash pneumatic dryer.Table 3. Results for a superheated steam pneumatic dryer of tobacco.

Dimension and Parameters of Superheated Steam Pneumatic Dryer Results
Length, m 8.02 Diameter, m 0.5 Dryer capacity, tons/h 3 Inlet steam temperature, °C 300 Steam velocity, m/s 30 Coefficient of heating medium excess, kg/kg 9 Material dwell time, s 0.4

Fire and Explosion Properties of Tobacco
The explosion and fire properties of tobacco before and after using the superheated steam-drying process at 160, 170, 180, and 190 °C were experimentally investigated.The influence of the impact of drying conditions on the explosion sensitivity and severity parameters was noted, and it is shown in Table 4 and Figure 15.The results indicate that the higher the drying temperature, the more likely is the dust tobacco cloud to be ignited, the faster the explosion flame propagation, and the greater the explosion severity.

Fire and Explosion Properties of Tobacco
The explosion and fire properties of tobacco before and after using the superheated steam-drying process at 160, 170, 180, and 190 • C were experimentally investigated.The influence of the impact of drying conditions on the explosion sensitivity and severity parameters was noted, and it is shown in Table 4 and Figure 15.The results indicate that the higher the drying temperature, the more likely is the dust tobacco cloud to be ignited, the faster the explosion flame propagation, and the greater the explosion severity.Tobacco dust is of weak explosion class (St 1) [7].Dust obtained by drying with superheated steam at 190 • C is characterized by the highest value of K stmax , amounting to 109 ± 14 m•bar•s −1 , and the highest explosion pressure jump (405 ± 32 bar/s), as well as the maximum explosion pressure (6.7 ± 0.3 bar), which may cause the most significant damage during a potential event in premises where tobacco dust is dried with superheated steam, especially inside enclosed spaces.The value of the maximum pressure increase rate over time, K stmax , was observed at the same dust concentration level at which the maximum explosion pressure values were estimated (Figure 15).± 14 m•bar•s −1 , and the highest explosion pressure jump (405 ± 32 bar/s), as well as the maximum explosion pressure (6.7 ± 0.3 bar), which may cause the most significant damage during a potential event in premises where tobacco dust is dried with superheated steam, especially inside enclosed spaces.The value of the maximum pressure increase rate over time, Kstmax, was observed at the same dust concentration level at which the maximum explosion pressure values were estimated (Figure 15).Studies on the minimum explosive concentration revealed that the lowest concentration value for which tobacco dust before being dried with superheated steam may form an explosive atmosphere amounts to 125 g/m 3 .The MEC values of the dust after drying at 160, 170, 180, and 190 °C are lower, and the tobacco falls within the explosive range as soon as a concentration level of 60 g/m 3 is reached.
The MIE value of a tobacco cloud is indispensable for the assessment of fire and explosion hazards in a superheated-steam-drying system.According to the results summarized in Table 4, tobacco dust before and after drying is susceptible to ignition from the initials over 1000 mJ.The MIE value identified in laboratory conditions can be lower in process conditions, since the air temperature during drying is higher.The probability of tobacco dust cloud ignition even in optimum conditions is low, as the ignition device energies used in experimental tests are much higher than the energy of authentic ignition sources.Dust explosion seems highly unlikely in factories where dust distribution is generally non-uniform, contrary to the test conditions.Despite many investigations, no single case of an explosion in tobacco factories has been discovered.
The minimum ignition temperature of tobacco dust, Tcl, after superheated steam-drying at 190 °C is 90 °C higher than the minimum ignition temperature of dust, which can be generated upon contact with non-thermally treated tobacco.It is caused by a higher moisture content in the material and reduced dust dispersion by increasing interparticle cohesion.The influence of the temperature of tobacco drying with superheated steam on the minimum ignition temperature of the tobacco-dust layer in contact with the surface is not so evident.For instance, the minimum ignition temperature value of dry tobacco dust, subject to tests, amounted to MIT5mm 290 °C ± 2.9, while for tobacco dust before drying, the MIT5mm at the heating-plate constant temperature amounted to 320 °C ± 2.9.When analyzing the minimum ignition temperature, one should note that the MIT depends on the deposited-dust-layer thickness.Whenever the temperature was Studies on the minimum explosive concentration revealed that the lowest concentration value for which tobacco dust before being dried with superheated steam may form an explosive atmosphere amounts to 125 g/m 3 .The MEC values of the dust after drying at 160, 170, 180, and 190 • C are lower, and the tobacco falls within the explosive range as soon as a concentration level of 60 g/m 3 is reached.
The MIE value of a tobacco cloud is indispensable for the assessment of fire and explosion hazards in a superheated-steam-drying system.According to the results summarized in Table 4, tobacco dust before and after drying is susceptible to ignition from the initials over 1000 mJ.The MIE value identified in laboratory conditions can be lower in process conditions, since the air temperature during drying is higher.The probability of tobacco dust cloud ignition even in optimum conditions is low, as the ignition device energies used in experimental tests are much higher than the energy of authentic ignition sources.Dust explosion seems highly unlikely in factories where dust distribution is generally non-uniform, contrary to the test conditions.Despite many investigations, no single case of an explosion in tobacco factories has been discovered.
The minimum ignition temperature of tobacco dust, T cl , after superheated steamdrying at 190 • C is 90 • C higher than the minimum ignition temperature of dust, which can be generated upon contact with non-thermally treated tobacco.It is caused by a higher moisture content in the material and reduced dust dispersion by increasing interparticle cohesion.The influence of the temperature of tobacco drying with superheated steam on the minimum ignition temperature of the tobacco-dust layer in contact with the surface is not so evident.For instance, the minimum ignition temperature value of dry tobacco dust, subject to tests, amounted to MIT 5mm 290 • C ± 2.9, while for tobacco dust before drying, the MIT 5mm at the heating-plate constant temperature amounted to 320 • C ± 2.9.When analyzing the minimum ignition temperature, one should note that the MIT depends on the deposited-dust-layer thickness.Whenever the temperature was measured for a 50 mm-thick layer, the MIT was on average 10-30 • C lower than the temperature for a 5 mm-thick dust layer.The identification was also carried out at a constant rate of heat generation (CRHG method).In the standard CTHP method, the ignition is assumed to occur if the temperature in the dust layer exceeds the heating-plate temperature by at least 250 • C. In the CRHG method, the heating-plate temperature changes during the measurements.In the reference case, the momentary temperature of the heating plate was taken into consideration by controlling the tobacco dust and heating-plate-temperature difference.For a 50 mm-thick dust layer, the heating-plate temperature increased for the initial 40 min to become stable afterwards.Dust ignition in the thin layer occurred a long time before the heating-plate temperature became stable.A thicker dust layer reduces the MIT and extends the ignition time in the layer.The minimum ignition temperature value for tobacco dust after drying at 190 • C, tested at a constant increase of the heating-plate temperature, amounted to MIT 5mm 260 • C ± 3.2, while for tobacco dust before drying, it was 280 • C ± 2.9.The differences become less significant as the dust layer thickness increases.For 50 mm-thick layers, the differences stabilize at 10 • C, which corresponds to the heating-plate-temperature rise when the CRHG method is used for its determination.The results reveal that in some cases, the standard method of a dust layer minimum ignition temperature determination may lead to hazard misestimation, whereby it may apply to the cases of the parameter determination for thick and thin dust layers.The maximum temperature for superheated-steam-drying equipment coated with an up to 5 mm-thick tobacco layer should not exceed a temperature of 75 • C lower than the dust layer minimum ignition temperature, amounting to 215 • C in the reference case, according to the standard method, or 195 • C based on the results according to the CRHG method.Based on experimental test results and previous experiments, it is recommended that general fire prevention measures be used.The drying process should be carried out in such a way that dustiness is avoided in the buildings.In order to reduce dust deposition, a cleaning procedure including cleaning at regular intervals should be followed including wet cleaning of equipment, or the equipment should be vacuumed with EX design devices meeting ATEX requirements for dust explosion hazard zones to avoid airborne dust [30,31].Effective equipment maintenance systems, regular staff training on safety rules and fire regulations, observance of the regulations concerning hot work, and hot-work permits should be implemented.

Validation of Results Obtained from the Proposed Model
The superheated-steam-drying results of tobacco obtained from the proposed model were compared to experimental data summarized in Reference [4] that were obtained during test runs made by ITM on CLED installed in Minsk, Belarus, and performed in January and April 2001.Two types of cut tobacco blends were analyzed, T1 and T2, differing in their recipe.Twelve runs were performed in which the following parameters were measured: tobacco inlet and moisture content, steam inlet and outlet temperature, inlet tobacco flow rate, inlet steam flow rate, tobacco size distribution, and tobacco filling power before and after CLED.The inlet steam velocity occurred in a range of 27 to 30 m/s.Pipe diameter was 0.5 m, and total pipe length, including horizontal section, vertical section, and two elbows, was 8.02 m.In simulations, the value of empirical parameter m in the kinetic equation for drying rate was taken as 0.4.The experiments described above allow only for comparison of two parameters measured and simulated by the model, tobacco moisture content after dryer, and steam temperature at the pipe exit.The tobacco moisture content is best measured after the dryer pipe.Here, however, it was measured after the cyclone.The comparison of measured and simulated values is shown in Figures 16 and 17.The error-of-fit does not exceed several percents on average.

Conclusions
A simulation is a convenient tool for developing and scaling up superhea flash dryers for tobacco.Even a simple one-dimensional model of isometric perse particle drying was proven useful in predicting the drying process in a C proposed model is suitable for a full-size CLED where the drying pipe contribu overall drying effect is substantial.In a "baby CLED", the share of the drying tribution is smaller and probably comparable to that of the cyclone.
The experimental studies indicate that the higher the drying temperature bacco, the more likely the ignition of the dust cloud, the greater the explosion and the faster the propagation of the explosion flame.Lack of oxygen in the su steam dryer led us to us assume that this process is safe.However, there is a p of occurrence of fire or explosion due to the generation of tobacco dust.The m temperature for superheated steam-drying equipment coated with up to 5 tobacco layer should not exceed a temperature of 75 °C lower than the dust-la mum ignition temperature, amounting to 215 °C in the reference case, accord standard method, or 195 °C based on the results according to the CRHG me results allow the fire and explosion hazards to be minimized and justify the ne ensuring comprehensive and frequent removal of accumulated tobacco dust, t of drying technological regimes, and the required frequency of maintenance an tion activities of equipment.
The above work illustrates present-day possibilities in a simulation and stages: measure and collect property data, build a model and solve it with rea

Conclusions
A simulation is a convenient tool for developing and scaling up superheated steam flash dryers for tobacco.Even a simple one-dimensional model of isometric monodisperse particle drying was proven useful in predicting the drying process in a CLED.The proposed model is suitable for a full-size CLED where the drying pipe contribution to the overall drying effect is substantial.In a "baby CLED", the share of the drying pipe contribution is smaller and probably comparable to that of the cyclone.
The experimental studies indicate that the higher the drying temperature of the tobacco, the more likely the ignition of the dust cloud, the greater the explosion severity, and the faster the propagation of the explosion flame.Lack of oxygen in the superheated steam dryer led us to us assume that this process is safe.However, there is a probability of occurrence of fire or explosion due to the generation of tobacco dust.The maximum temperature for superheated steam-drying equipment coated with up to 5 mm-thick tobacco layer should not exceed a temperature of 75 • C lower than the dust-layer minimum ignition temperature, amounting to 215 • C in the reference case, according to the standard method, or 195 • C based on the results according to the CRHG method.The results allow the fire and explosion hazards to be minimized and justify the necessity of ensuring comprehensive and frequent removal of accumulated tobacco dust, the respect of drying technological regimes, and the required frequency of maintenance and inspection activities of equipment.
The above work illustrates present-day possibilities in a simulation and its three stages: measure and collect property data, build a model and solve it with readily available mathematical software, and validate the results against industrial process data.The first three stages can be performed by academia; the third is required from industry.As a result of such cooperation, useful simulation models can be created for almost all industrial processes.
Further work is required.First, better sorption isotherms and isobars must be obtained through extensive experiments.Second, a better representation of particles in the model is necessary.For this purpose, an experimental setup for the determination of drag coefficient for ribbon-like particles is under construction in our laboratory.Introducing polydispersity into the model will then be possible.Finally, the drying-rate equation must be experimentally verified.Modification of the setup shown in Figure 10 will allow for the direct measurement of the drying kinetics in controlled conditions.Finally, the model of the expansion process can be introduced and solved simultaneously with the present model.
.[13] derived Equation (5) based on cut lamina obtained from the BAT company for blend 72 in the range of 20 to 60 °C.

Figure 1 .
Figure 1.Comparison of experimental data obtained by References[11,12] and several empirical sorption isotherm equations for tobacco.

Figure 1 .
Figure 1.Comparison of experimental data obtained by References[11,12] and several empirical sorption isotherm equations for tobacco.
traditionally presented to mean moisture content dependent on temperature.Exemplary sorption isobars for tobacco lamina and stem measured in this setup are shown in Figure3.

Figure 4 .
Figure 4. Data from Figure 3 replotted in a w vs. X coordinates (lamina only).

Figure 5 .
Figure 5. Three-dimensional sorption-equilibrium equation fitted to our experimental results (blue dots).The equation has the following form:

Energies 2021 ,
14, x FOR PEER REVIEW 12 of 23 through another airlock.The steam leaving the cyclone is partly vented and partly circulated by a circulating fan.It is then indirectly reheated in the superheater and returned to the dryer pipe.

Figure 10 .
Figure 10.Schematic of a superheated steam flash dryer.

Figure 10 .
Figure 10.Schematic of a superheated steam flash dryer.

Figure 11 .
Figure 11.Schematic of the momentum, heat, and mass balances over a differential volume element of the dryer pipe.

Figure 11 .
Figure 11.Schematic of the momentum, heat, and mass balances over a differential volume element of the dryer pipe.

Figure 12 .
Figure 12.Moisture-content change in the function of the apparatus length.

Figure 13 .
Figure 13.Steam temperature (dotted line) and tobacco temperature (solid line) in the functi the apparatus length.

Figure 12 .
Figure 12.Moisture-content change in the function of the apparatus length.

Figure 12 .
Figure 12.Moisture-content change in the function of the apparatus length.

Figure 13 .
Figure 13.Steam temperature (dotted line) and tobacco temperature (solid line) in the functio the apparatus length.

Figure 13 .
Figure 13.Steam temperature (dotted line) and tobacco temperature (solid line) in the function of the apparatus length.

Table 3 . 4 Figure 14 .
Figure 14.Steam velocity (dotted line) and tobacco velocity (solid line) in the function of the apparatus length.

Figure 14 .
Figure 14.Steam velocity (dotted line) and tobacco velocity (solid line) in the function of the apparatus length.

Figure 15 .
Figure 15.Impact of tobacco dust concentration on values of (a) P max , (b) dP/dt ex , and (c) K st .

Energies 2021 , 23 Figure 16 .
Figure 16.Results of experimental and modeling final tobacco moisture content.

Figure 16 . 22 Figure 16 .
Figure 16.Results of experimental and modeling final tobacco moisture content.

Figure 17 .
Figure 17.Results of experimental and modeling final steam temperature.

Figure 17 .
Figure 17.Results of experimental and modeling final steam temperature.

Table 1 .
Chemical and physical composition of the tobacco (wt.-%, dry basis).
* Calculated by difference.

Table 4 .
Results of explosion and fire properties of tobacco.

Table 4 .
Results of explosion and fire properties of tobacco.