Numerical analysis of single and multiple particles of Belchatow lignite dried in superheated steam

Low production costs have contributed to the important role of lignite in the energy mixes of numerous countries worldwide. High moisture content, though, diminishes the applicability of lignite in power generation. Superheated steam drying is a prospective method of raising the calorific value of this fuel. This study describes the numerical model of superheated steam drying of lignite from the Belchatow mine in Poland in two aspects: single and multi-particle. The experimental investigation preceded the numerical analysis and provided the necessary data for the preparation and verification of the model. Spheres of 2.5 to 30 mm in diameter were exposed to the drying medium at the temperature range of 110 to 170 °C. The drying kinetics were described in the form of moisture content, drying rate and temperature profile curves against time. Basic coal properties, such as density or specific heat, as well as the mechanisms of heat and mass transfer in the particular stages of the process laid the foundations for the model construction. The model illustrated the drying behavior of a single particle in the entire range of steam temperature as well as the sample diameter. Furthermore, the numerical analyses of coal batches containing particles of various sizes were conducted to reflect the operating conditions of the dryer. They were followed by deliberation on the calorific value improvement achieved by drying, in terms of coal ingredients, power plant efficiency and dryer input composition. The initial period of drying was found crucial for upgrading the quality of coal. The accuracy of the model is capable of further improvement regarding the process parameters.


Introduction
Lignite (also known as brown coal) has a vital position in a number of power generation systems, particularly in Europe. Shallow deposits of this widely distributed resource are cheap in extraction, therefore the total cost of electricity generation in lignite-fired power plants is relatively low [1]. On the other hand, coalification processes did not develop as much as in the case of bituminous coal, what negatively influenced the elemental carbon share in brown coal. Unfortunately, in terms of power generation applicability, water content often prevails over combustible matter, what undermines the calorific value of this energy carrier. The occurrence of water also makes transportation and storage economically unjustified. Yet, the aforementioned low production costs combined with significant volatile matter share, which boosts gasification potential, makes upgrading the quality of lignite by drying worthy of consideration.
The key aspect of drying in any technological process is maximization of the ratio between the achieved increase in material quality (measured by the expiration date, water percentage, calorific value etc. [2][3][4]) and the heat consumed by the drying system. Superheated steam drying is one of the prospective methods, promising the possibility of latent heat recovery, what might considerably reduce the heat demand of fuel preparation devices.
This drying medium has been discussed since the early twentieth century [5] and its application as an alternative to hot air for dewatering granular objects was proposed in 1951 [6]. The advantages of superheated steam drying include the reduction of spontaneous ignition hazard and, above a certain temperature, higher drying efficiency than hot air [7]. The key benefit, though, is the capability for latent heat recovery. One of the applications realizing such an idea is the superheated steam fluidized bed dryer [8,9]. In this concept, the drying system operating in a self-heat recuperation configuration enables the collection of steam from the moisture in the raw fuel, followed by the recirculation of its latent and sensible heat for the drying of another fuel portion [10,11]. By that means, the heat of the phase transition is not given off to the ambience, what results in the decline of energy consumption, thus elevating the thermal efficiency of the system.
The reliable approach to the design and construction of a superheated steam drying (SSD) system coupled with a power generation unit, entails a thorough study on fuel taken into consideration. This kind of works often include numerical analysis as a cheap and precise manner of describing heat and mass transfer phenomena which occur in the drying process. The food industry has commonly adapted such a method of drying optimization. The model of fixed bed dryer of brewers' spent grain using superheated steam one of the examples [12]. Applying the finite difference method, changes in moisture content and temperature in the slice of pork were simulated for SSD [2]. Power engineering also reached for the superheated steam drying models, what resulted in the numerical analysis of the combined heat and power plant fueled by corn ethanol [4] as well as the self-heat recuperative fluidized bed dryer of biomass, which consumed 95% less energy than conventional systems [3]. Coal became the object of drying simulations with regard to its granular form [13]. Another study [14] came out with a single particle SSD model prepared for Australian lignite from the Loy Yang mine. The work was followed by a simulation of multi-object drying [15]. The experimental methodology of the two last mentioned works was analogical as in this study.
A previously verified [16] numerical model of superheated steam drying of Belchatow lignite was applied to this study, relying on the experimental research conducted beforehand on the drying characteristics of Polish lignite excavated from the Belchatow deposit. Outcomes acquired from the drying of variously sized objects is applied in this study to numerically analyze multi-particle SSD, which can reflect the actual operation mode of a dryer more precisely. The simulation is utilized to evaluate the potential of calorific value enhancement, including diverse scenarios of input composition and its influence on the performance of the superheated steam dryer.

Samples
Spherical objects of four diameters were included in the examination: 30, 10, 5 and 2.5 mm. The investigated material originated in the Belchatow lignite mine and was excavated in a batch denominated B2013. Each sample type was prepared by rolling a piece of lignite on a plate with a series of punctuated holes, which diameter declined stepwise to the demanded size. The particles of various sizes differed in methods of suspension inside the test section and the temperature measurement. Figure 1 presents the lignite samples together with the installation methods. Particles of 30 (Fig. 1a) and 10 mm (Fig. 1b) were equipped with two thermocouples per each, measuring temperatures in the half of the radius (later called Bmidpoint^) and the center. The other function of the thermocouples was the suspension of the sphere in the experimental vessel, however, the largest sample was additionally stabilized with polyester thread passed through its interior. One thermocouple, measuring temperature in the center, was placed in the object of 5 mm (Fig. 1c). In the case of 2.5 mm (Fig. 1d) no thermocouple was applied and a quadruple glass hanger was used to suspend the four samples at a time.

Setup and procedure
The scheme of the experimental apparatus is presented in Fig. 2. Before the tests, the examined sphere (or spheres) was inserted into the test cylinder, which was heated and insulated in order to assure the thermal stability. The distilled water was used to generate superheated steam. It was degassed, pumped to an evaporator, turned to steam and then superheated. The gaseous medium filled the test vessel through a baffle plate, used for the uniform distribution of steam within the chamber. The input rate of the water into pump and degasser was 8 cm 3 min −1 .
The temperature of the steam varied between 110 and 170°C, depending on the test case, while its pressure was kept at 1 atmosphere. The assumptions of water condensation and coal volatilization limited the temperature range from the bottom and top, respectively. The decomposition of lignite structure was not observed [17] below 180°C, hence the water removal is treated as an exclusive reason for the descent of the sample weight.
The temperature measurement of the coal interior was realized depending on the particle diameter, as mentioned in the previous subsection. Moreover, an independent measurement of coal surface temperature was conducted by means of thermography for objects of all sizes. The electronic balance with a metal rod, on which the lignite was suspended, enabled the weight measurement. The completion of steam drying was followed by the substitution of medium to nitrogen, which was applied to remove the residual water within the coal. More details of the experimental procedure, with respect to each type of examined object, are given in the previous works [18][19][20].

Drying indicators
Relying on the mass of the sample after the completion of nitrogen drying, two parameters, applied for quantitative description of the drying process, were calculated. They are introduced below.
Moisture content X, is expressed as a ratio of water mass in coal (variable in time) to the mass of the dry coal remaining after drying. The formula is given in Eq. 1: Water percentage WP, represents the share of the water in the lignite sample, related to the initial mass of the sample. It is formulated as mass percentage and shown in Eq. 2: 3 Model of superheated steam drying

Physical model
The physical model of drying was assumed to be a sphere, comprising of three phases of different thermophysical properties: dry coal matter, free water and bound water. Dry coal matter represents the solid mixture of carbon, carbohydrates and ash that are not subject to evaporation and therefore their total mass does not vary through the drying process. Free water is assumed to evaporate at saturation temperature, which is 100°C under atmospheric pressure, while bound water turns to steam above this temperature. The model was constructed on the assumption that the object is a perfect sphere of isotropic properties. The parameters of dry coal, such as density ρ c , specific heat c c or thermal conductivity k c , are constant during drying. On the other hand, the analogical properties of water are significantly dependent The significantly important parameters for conducting this study are given in Table 1. For the purpose of simulation, some physical and chemical properties, due to the lack of proper data on Belchatow coal, were adapted from other lignite types. A few of them were derived directly for Belchatow, based on empirical data, found in an auxiliary investigation or acquired from commissioned research.

Mathematical model
Due to the above-mentioned isotropy of the properties, a onedimensional model of transient heat transfer, along the sample radius, might be applied. The formula describing the temperature field under those conditions, is given by Eq. 3: where q̇V stands for the volumetric heat sink represents heat loss by the evaporation of free and bound water (Section 3.3.3.). The first component of the right-hand side of Eq. 3 can also be expressed in that case as: The initial condition is: The boundary conditions reflect the symmetry condition and the transfer of heat from the steam to the surface of a sample, respectively: Heat transfer coefficient h in Eq. 7 is a function of the steam temperature (the surrounding space around the sample is fulfilled with steam) and sample diameter. It was previously derived relying on thermodynamic deliberation on superheated steam drying, performed for the same lignite type [18]. The used formula is given in Table 1.
Whether the moisture content X n in layer n is higher than equilibrium moisture content X eq in the current temperature (see Table 1), the diffusion of water to the layer n-1 occurs. The mass flow is governed by the Fick's law: where ρ b,w is a bulk density or mass concentration of water within lignite and D, which value is given in Table 1, stands for diffusion coefficient of free water [14].

Discretization scheme
In order to transform the differential equations into algebraic formulas, the control-volume method was applied linear shrinkage 1-l/l ini adapted from [14] diffusion coefficient of free water D LY 3.0 × 10 −9 m s −2 adapted from [14] [ 21,22]. Assuming element n within a sphere, which volume is V n and area equals to A n , Eq. 3 can be integrated on this element as: The first component of right-hand side in Eq. 9 is equal to: where n is a normal vector. Density ρ n , thermal conductivity λ n and specific heat c n are assumed to be constant over the entire element n. Therefore, the apparent values dependent of volume fraction of coal and steam in each time step are used in the calculation domains, so Eqs. 9 and 10 might be integrated to: and: where the right-side of Eq. 12 represent heat fluxes by the conduction driven into and out of layer n.
As to the mass diffusion between the layers, integration of Eq. 8 from r n to r n + 1 gives:

Particular discretization assumptions
Besides those introduced in the previous section, the following assumptions related to the discretization need to be listed: & The uniform spatial distribution of the elements along the radius was applied, forming a series of concentric layers characterized with an initial width equal to 2b. Only the initial width of extreme elements, number 1 and N + 1, was equal to b in order to avoid overlapping of the control-volume balances (Fig. 3). & The sample division N was chosen to be 50, what ensures the accuracy of the numerical prediction independent from the grid size of the computational domain.
& The simulation was performed using the explicit method with a specific value of the time step width Δt adjusted to the conditions: By default, Δt was equal to 10 −3 s, nonetheless it was set at 10 −4 s for samples smaller than 2.5 mm in diameter and for objects of diameter smaller than 1 mm it was 10 −5 s.

Specificity of particular drying stages
During the drying process, the mechanisms of heat and mass transfer vary, with the occurrence of water and temperature level as the main criteria. The numerical model, prepared with regard to the assumptions pointed out in sections 3.1-3.3 was structurally segmented with respect to the mechanism of dewatering. Consecutive drying stages, characterized with various sets of numerical equations, are described in the following subsections. However, the heat conduction between adjacent layers was a common mechanism for all stages and is presented by Eq. 15.
Condensation of water on the surface The temperature of the sample in the initial stage of drying is lower than the saturation temperature. That is why, a condensation of water occurs on the coal surface. The total heat input to the dried object is covered by this mechanism, what is reflected by Eq. 16: The heat consumed at that time causes the increase in the temperature of the dry coal and moisture within the sphere.
Evaporation of water from the surface As the temperature of the surface of the coal sample reaches 100°C, heat transfer to the sample becomes proportional to the temperature difference between the superheated steam and saturation The heat transfer coefficient applied to this period was evaluated relying on the mass evaporated during the constant drying rate period, as shown in Eq. 7. Then, it was linearized in the form of function, of which the variables are the temperature of the superheated steam and the reciprocal of the sample diameter, as seen in Table 1. In this stage of drying, the water film formed by condensation and exudation, remained on the surface. In the case of the larger samples it tended to exceed the critical mass and a part of it fell down in the form of droplets. The assumptions beyond the droplet formation were adapted from another study [15].
Evaporation of free water The mechanism beyond heat transfer to the sample surface does not change after the removal of water from the surface, however, heat input becomes proportional to the difference between the gas temperature and current temperature, which exceeds 100°C and gradually increases (Eq. 19).
Free water in coal, similarly to water on the surface, evaporates at saturation temperature. The relation between heat consumed by the layer and the mass of water evaporated from it can be formulated as: Due to evaporation of water, the moisture content gradient occurs between the layers. The gradient induces water transport between the layers. The bulk density of the water may be expressed as: Taking formula 21 into account, the transfer mechanism of free water, is illustrated by Eq. 22.
Evaporation of bound water When the temperature of a particular layer exceeds 100°C, it signalizes that free water evaporation is over and bound water only remains in the coal sample. Bound water is held firmly within the lignite structure, therefore its evaporation poses a greater challenge than the removal of free water. The evaporation rate at a certain moment is limited, therefore a surplus of heat input is allocated on raising the temperature. Instead of the simple kinetic of latent heat consumed in the evaporation mechanism, an enthalpy change of bound water evaporation which can be described as the function of the layer temperature, applies at that stage. The heat consumption during bound water evaporation is presented in Eq. 23.
The observation of the sample appearance indicated a considerable shrinkage phenomenon below the moisture content of around 0.6 [23]. Such a value is characteristic for the stages of drying when the surface temperature is over 100°C, what means bound water is being evaporated from the shallow parts of the coal. Hence, the shrinkage equation, dependent on the moisture content, was implemented to modify the particle size during simulation, as presented in Table 1.

Results and discussion
4.1 Simulation of single particle drying in superheated steam

Comparison of empirical and numerical results
The simulation results were put together with the experimental data in Fig. 4. The changes in the moisture content, the drying rate (which is a time derivative of the moisture content) and the temperature profiles (surface and center) were included in the drying characteristics. Midpoint temperatures were neglected on purpose due to the clarity of the presented graphs. The comparison was made for 10 mm samples, exposed to the steam at the entire temperature range. The empirical results were appropriated from the previous study [16]. The simulation was conducted for the initial water percentage agreed in each case with the experimental value in order to check the coherence of the results. The other independent variables, in particular those indicated in Table 1 remained unaffected for all instances.
Concluding from the charts, a high temperature of the superheated steam entailed the elevated speed of drying. For that reason, water from lignite was discarded sooner and the sample reached thermal equilibrium with the steam, reducing the time required for the completion of the process.
The rise of the temperature, especially in the case of the temperature at the center of the sample, was more gradual for the experimental results and it occurred sooner than the time predicted by the simulation. Moreover, the apparent distinction between the surface and center temperature predicted numerically, does not occur in the experiments with the exception for a) and b) instances. The model adapted in this study assumed the increase of the temperature over 100°C in the particular layer only after the complete removal of the free water in this layer. Uniformity of structure within the lignite and deformations occurring in the drying process may influence the boundary between the free and bound water moved at various paces within the different cross-sections, and made the temperature profile inconsistent with the calculated prediction as seen in Fig. 4.
Similarly to the set of results presented for the particular diameter, another one was prepared for the temperature fixed at 150°C and a wide range of sizes. It is illustrated in Fig. 5. Note that due to a lack of interior temperature measurement, only one experimental curve representing the surface is given in the temperature profile of the 2.5 mm sample chart.
The major disparity is visible in the time and speed of drying, however, the trends followed by the corresponding curves are alike. The oscillations on the drying rate experimentally obtained, are enhanced with the decrease in diameter and are resulted from the precision of the electronic balance (0.1 mg), which influences the smaller objects to a greater degree. In the case of Fig. 5a) the predicted drying rate slightly overshoots the empirical level, unlike the other cases, but the agreement between the analogical parts of the temperature profile seems the most precise. The same observations apply to the instance of 170°C in Fig. 4a). The issue requires further improvements, including deliberation on the empirically obtained heat transfer coefficient and the parameters related to the thermal conduction within the sphere.

Prediction of temperature dependence on drying performance
In order to quantify the particular features of drying kinetics, another computation was conducted at the same temperature range with the interval of 10°C. The initial water percentage was assumed to be 51.3%, a value averaged over multiple experimental attempts on B2013 coal. The selected parameters are presented in Table 2.
The time of drying t dry predicted by the simulation inclines sharply with the decrease of steam temperature, especially in the vicinity of 100°C. For instance, temperature drop from 120 to 110°C induces the extension of the drying process by 76%. The drying rate (dX/dt) max achieves peak value at the temperature of 170°C. This value was simulated during the period, when dewatering is the most intensethe constant drying rate period (CDRP). The duration of this stage, t CDRP , was expected to last from 1.7% (170°C) to 8.9% (110°C) of the entire drying time. This relatively short period of the process has a decisive impact on the final effect of drying, what is displayed in section 4.2.2. The increase in the center temperature of over 100°C, triggered by the complete evaporation of the free water in the sample, happens around half of the total time for higher steam temperatures and reduced to 1/3 when approaching the saturation temperature. It is marked as t free . The largest gap between the temperatures of the surface and center (ΔT 1,51,max ) was computed at 32.3°C, while in the opposite case it was as low as 1.8°C. Therefore, the objects dried at high temperatures, what results in considerable thermal stresses, are more exposed to the hazard of structure deformation, what was illustrated in the video analysis included in the previous study [18]. An analogical juxtaposition of the parameters for the drying of the 5 mm sample is gathered in Table 3. The time required for drying is around 3 times shorter, regardless of the gas temperature. The constant drying rate period occupies a similar part of the entire process (3.0% to 8.2%), however, the peak speed in dewatering is 2.5 times higher than of a larger object. The evaporation of the free water in the case of the 5 mm sample at 120°C takes as much time as the analogical process conducted at 170°C for the 10 mm object. Due to the smaller dimensions of the 5 mm sample, heat is transported through the sample much faster, what reduces the maximal level of the temperature gradient and limits the destructive impact of thermal stress on the coal structure.

Quantitative verification of numerical computations
In order to comprehensively analyse the coherence of the computational and empirical results, a summary of the values attained at the end of drying, total time of drying and residual water percentage, was prepared for the 10 mm sample in Fig. 6. The simulated results (square marker) are given for the entire temperature range with an interval of 10°C, while three representative experimental cases and the averaged value those from Ref. [16] are presented only for the conditions applied in the tests.
The experimental attempts served to formulate the equation, which was applied in the model for the estimation of the residual amount of water with regard to the drying temperature. Hence, a good agreement was achieved. In the case of the drying time, numerical values oscillated around the empirical results as well as the dashed line, which represents the approximation equation (Eq. 24) for the time of drying proposed in the previous work of the authors [18]. The empirical results of the drying time t exp from Fig. 6 are shown altogether with their standard deviations s(t exp ) in Table 4. With reference to the values of t exp , the root-meansquare error related to both approximation R(t dry,app ) and model R(t dry ) was evaluated. The coherence of the empirical results increases with the rise of the drying temperature. Similarly, the approximated and numerical results evince greater consistence with the experiment, however, the simulation offers precise prediction only for higher temperatures. The case of 110°C was overshot by the model, whereas drying in 130°C lasted longer than the one predicted from the simulation.

Drying of various coal assortments
The coal input to the grinder, which supplies the dryer, comes in various size assortment. As demonstrated in Fig. 5, drying time significantly decreases with the decline in object diameter. Hence, the energy consumed by the grinder on granulating  the fuel entails reduction of time, and thus energy, required for the operation of the dryer. The present study investigates the latter component of this interdependence. The diameters of lignite spheres were assumed to fit the Rosin-Rammler equation, applied frequently to describe size distribution in fragmented minerals [24,25]: In the formula above, f(d) stands for a mass percent ratio of particles of diameter higher than d to total mass of the lignite group, d e is a characteristic particle size for the particular group and B is a uniformity constant equal to 1 [15]. Due to the asymptotical character of Eq. 25 which tends to 0 as the diameter tends to infinity, the d max value is introduced to represent the size over which the particles account for 1% of the mass distribution. Hence, to simplify the calculation, the mass percent ratio is divided into 10 subgroups. The derivative of f(d) is presented for the case of d max = 80 mm in Fig. 7. Particles larger than 80 mm are not taken into account, subgroup 1 represents 9% of the total mass distribution, while subgroups 2-10 stand for 10% each. The square markers in the bottom axis illustrate the mean particle diameters in the particular subgroups.
Three scenarios of size distribution in the lignite batch exposed to the superheated steam were taken into consideration. The maximum diameter values of Bfine^, Bmedium^and Bthick^lignite groups were 40, 80 and 300 mm, respectively, what corresponds to the industrial practice of the Belchatow coal producer [26,27]. The representative values for corresponding subgroups are shown in Table 5.
In Fig. 8, for the first 60 min of the process, the moisture content and drying rate of the multi-particle simulation for medium coal is plotted alongside the moisture content curves for the particular subgroups. The stepwise decline of the drying rate corresponds to the values attained by the samples of particular size at the subsequent stages of the process.
The visualized drying process of the particle grouped defined in the size quality, illustrated as the curves of the drying rate along the moisture apparent content, is presented in Fig. 9. The composition of the coal batch strongly influence the drying rate [28]. The dominant content of the smaller samples boosts the peak speed of the dewatering, determining the shape of the correlation curves. The computational graphs display irregularities, which were caused by a series of falling droplets formed from the water on the surface. Such a phenomenon, specific for the objects of around 30 mm and higher [19], occurred rarely in the drying of other samples [18]. The

The influence of drying on the efficiency of power generation
Two parameters were used to describe the drying impact on the optimization of power generation: thermal efficiency and carbon dioxide emission per generated unit of electricity. For the sake of the optimization analysis, a variable of water share, is introduced in Eq. 26. It expresses the absolute share of water within coal in the particular moment, unlike previously discussed water percentage relating to the initial mass of coal before drying.
In order to quantify those magnitudes, a brief insight into the coal composition is required. Table 6 contains the major components of B2013 lignite, relevant to the present analysis. The share of carbon and hydrogen was reported by a commissioned investigation with respect to coal deprived of water content (Bdry coal^column) and recalculated to as-received state (Braw coal^). The structural division into the dry and liquid part of lignite is shown in Fig. 10.
The calculation of carbon dioxide emission relates to the stoichiometry of carbon combustion. Every 12 kg of carbon binds with oxygen to form 44 kg of CO 2 .
Basing on the assumption that the carbon in the lignite is subject to complete combustion, the emission of CO 2 per kilogram of lignite is formulated as: where WS represents a variable of water share that decreases during drying. Emission of carbon dioxide needs to be related to energy generation. Two indicators are commonly used to describe the thermal effect of fuel combustion. The higher heating value (HHV) represents the total amount of heat produced from the mass unit fuel with the assumption of no heat losses on the evaporation of water present in the fuel or produced parallel to combustion. Basing on the HHV value at the air-dried state, the parameter related to the dry state is: The HHV and water share for the air-dried B2013 coal are 19.08 MJ kg −1 and 14.6 mass%, respectively, what results in HHVat 22.34 MJ kg −1 for the dry coal. Therefore, the formula for the higher heating value of lignite during the drying process is shown in Eq. 29:  The lower heating value (LHV), on the other hand, might be expressed as HHV less the heat needed for inherent moisture evaporation and water synthesized from hydrogen within dry coal. Given that 4 kg of hydrogen takes part in the synthesis of 36 kg of water, LHV with regard to coal in the drying process might be expressed as shown in Eq. 30.
Electricity generation per mass unit on fuel with certain LHV is directly proportional to thermal efficiency.
Subsequently, the Eq. 32 formula presents carbon emission per generated kilowatt-hour: For the purposes of this study, it should be recognized that power generation companies often refer to LHVof applied fuel and calculate thermal efficiency based on this parameter, while in the research publications [29,30], correlation between emission ratio and thermal efficiency rely on HHV. Therefore, thermal efficiency on HHV with regard to the one on LHV might be expressed as: The changes in higher and lower heating value against declining water percentage are presented in Fig. 11. Data points are distributed every 5% of WS, with the exception of the one, which marks 3.94%, a residual water share (equivalent to the residual water percentage of 1.98%) for the superheated steam drying at 150°C. The initial gap between HHV (10.21 MJ kg −1 ) and LHV (8.56 MJ kg −1 ) decreases as the drying progresses, reaching a minimum value of 0.97 MJ kg −1 at the end of drying. This final difference is attributed mostly to hydrogen content in dry coal part, as the theoretical margin (marked in red) for the completely dewatered lignite becomes 0.92 MJ kg −1 .
Enhancement in calorific value by means of drying involve economic and ecological benefits. Efficiency data provided by the electricity producer [31] relates to the LHV of raw lignite, what might be supported by information from the Belchatow power plant boilers manufacturer [32], providing data on the efficiency of the boiler accompanied by the LHV of supplying coal at the level of 8 MJ kg −1 . On this assumption, a chart in Fig. 12 was prepared describing the reduction of the CO 2 unit    [29,34] emission caused by drying. The horizontal axis represents thermal efficiency recalculated to HHV according to Eq. 33. Three cases, gathered in Table 7 are included with regard to the net thermal efficiency of the power generation unit. Note that Cases 1 and 2 refer to lignite power plants, while Case 3 to the bituminous power plant.
Particular data points in each case are placed in the same water share interval as in Fig. 11: beginning with 55% on the left-hand side and descending to 5% and 3.94% on the right. The level of carbon dioxide emission for raw coal in Case 1 is 1.054 kg kWh −1 , a value compatible with data provided by the Belchatow power plant (1.069 kg kWh −1 on average in 2015) [35]. Analogical values for Cases 2 and 3 are 0.954 and 0.871 kg kWh −1 . The total decline in the emission rate is 12.1% below the initial level, regardless the case. However, the distribution of data becomes denser as the water content drops, so it can be found out that the greatest profit for environmental efficiency is achieved in the early period of drying. The reference curves represent the carbon intensity along thermal efficiency for bituminous coal, oil and natural gas. Altogether with the correlation modelled hereby for lignite, they fit to the tendency that lower CO 2 emission levels are characteristic for fuels of higher hydrogen-to-carbon ratio (regarding dry part). Figure 13 illustrates the simulated progress of drying for 3 batch compositions introduced in section 4.2.1, related to the boost in lignite lower heating value. The biggest enhance in the quality of this fuel is achieved in the initial period of drying. Given 8.56 MJ kg −1 as raw coal's LHV, doubling the calorific value of Mixtures 2-4 (mass share of 30 mm samples at 10% or less) takes around 30 min for fine, 60 min for medium and over 400 min for thick coal. The differences extend as the drying continues, with regard to size distribution. For instance, whether the objective is to reach 18 MJ kg −1 (WS≈14%, WP≈8%), a time interval between fine and medium coal is 47 min. However, the time required for raising LHV to 20 MJ kg −1 (WS≈6%, WP≈3%), is 80 min in the case of the fine coal and almost 200 min for the medium one.
In conclusion, improving the quality of lignite and applied power generation technology might be beneficial in terms of ecology, nevertheless, the thermodynamic and financial costs of such efforts should also be taken into account. Therefore, an optimization of power generation unit (in particular fuel preparation devices: grinder and dryer), the general idea of which was proposed in the previous study [36], needs to be performed with regard to available industrial appliances.

Conclusions
Numerical modelling was used to simulate superheated steam drying for single and multiple samples made of Belchatow lignite, as an extension of the earlier work, concerning singular particles only [16]. The results of the investigation were presented in the form of drying charts containing experimental and computational curves as well as quantitative comparisons between those two data acquisition methods.
The application of multi-particle modelling was adapted in order to reflect a real operating mode of the drying system. By such means, a computational device is offered to calculate the time of drying for the demanded heating value of fuel. As a consequence, knowledge about dryer operational characteristics, may lead to heat input estimation. Taking into consideration several cases of fuel input composition, reflecting the actual assortments utilized in the Polish lignite industry, the potential enhancement of coal utilization effectivity was evaluated. The initial period of drying was found to be decisive in upgrading the quality of the fuel.
Further studies, complemented with the analysis of grinder energy demand as a function of coal fragmentation, may apply for the optimization of the drying system. Such works could also include intermediate scenarios, where different assortments are mixed prior to input or the effect of sieving the most thick coal subgroups.