The properties of the fabric samples for drying tests
1. Introduction
Hot air impingement is one of the most widely used methods for material drying. It is also the most traditional drying approach used in industrial process for various kinds of material, such as wood, paper, food, medicine and construction materials. Many research studies have been carried out to see how it is effectively used to process different types of material, and how it is implemented into the design of heat setting machines, such as spray dryers, conveyor dryers, tunnel dryers, fluidized bed dryers and drum dryers.
To study hot air impingement in textile and clothing industries, modeling of porous type fabric drying process would be the key study area. The heat and mass transfer principles are used as tools to assist with the investigation of the hot air impingement mechanism. The mechanism is usually treated as a mass transfer process of the moisture content from the porous material to the impinging air. The transfer of moisture content from the fabric material to the hot air stream is due to a heat transfer process under an in-equilibrium condition. The change of water phases is traditionally described by linear heat transfer equations. As a matter of fact, the driving force in the internal structure of porous materials is not a simple direct proportional relationship between energy exchange and the phase change of the interacting substances, i.e. air and water. Therefore non-linear analytical models based on the physical properties of fabrics will be proposed in this study to provide better simulation results. In the models, the parameters for modeling will be empirically determined and used to describe the drying phenomenon down to microscopic levels. The descriptions will involve the physical and mechanical properties of the drying materials such as mass density, flow viscosity, thermal conductivity, diffusion properties, cohesive properties and flow kinetics. Ip and Wan (2011) have suggested the strategies of using analytical techniques to determine the modeling parameters, and these methods will be investigated in greater depth in this research.
2. Three periods of a fabric drying cycle
Fabric is usually dried up for the purposes of storing or setting. Using thermal energy to dry up and perform setting has been the most traditional and effective method. In this study, heated air is used as a processing agent. Its physical properties will be changed in the gaining of moisture and the loss of thermal energy. The moisture in the fabric will change to vapor after gaining energy from air to create a mass transfer process. The reduction of moisture and increase of fabric temperature is a complicated heat/mass transfer process. Merely using linear conductive and convective heat transfer equations to model the process seems to be inadequate. Diffusion theories are therefore suggested to present the details of the drying process.
“Preheating”, “Constant drying” and “Falling drying” are the three periods of a fabric drying cycle as shown in Fig. 1. In the preheating period, most thermal energy is absorbed by water on the fabric surface because air is a poor thermal conductor. The mass transfer rate of water is not high in this period. When more thermal energy is absorbed, water on the fabric surface will change to vapor by evaporation at a rapid mass transfer rate. The water loss rate will keep constant depending upon the air temperature, velocity and atmospheric pressure. As the mass transfer rate of water is constant in this period, it is labeled as constant drying period. While the moisture content in the fabric is going down from the initial
Kowalski et al. (2007) has employed partial differentiation and numerical analysis tools to present thermo-mechanical properties of a drying process for porous materials. However, the modeling process is computational intensive and time consuming. It seems impractical to be used in the industrial drying process because a quick response is always needed to manage numerous varying conditions. The concept of using analytical approaches to model drying processes for fruits and sugar has been proved (Khazaei et al., 2008). However, its performance heavily relies upon the sample data and the reliability of the assumptions. The scopes of the research are therefore having rooms to improve the inadequacy. Objectives of the research are to explore and develop robust analytical models that can effectively simulate the characteristics of hot air impingement process for porous fabrics of different textile properties.
2.1. Research objectives
The research objectives in the Chapter are:
to investigate how to present the drying characteristics of porous type fabric using non-linear analytical models,
to evaluate the performance of the models in the simulation of drying process, and
to comment their accuracies in the modeling of different fabric types under various air setting conditions.
In this research, four non-linear analytical models will be studied. The modeling parameters for the models will be empirically determined. The performance of the model will be examined through a careful comparison of testing results. A drying test will be set up to assist the determination of the modeling parameters, and evaluate the performance of the developed models.
2.2. Equations for moisture flow through control volume
Using the approach of control volume to the describe characteristics of moisture flow in porous materials would be close to the phenomenon of fabric drying (Kowalski, 2003). The total mass of the constituents within the control volume will remain unchanged when the porous fabric volume shrinks or otherwise. A set of mass balance equations for each individual constituent, i.e. water, vapor and air in a drying process is given as:
where
2.3. Boundary conditions for the period of constant drying
The water evaporated inside fabric material is much less than that on the boundary surface in the constant drying period. Assuming that the mass flux
Wv = Wa = 0
Then, the mass balance Equations (1) – (3) can be rewritten for the calculation of moisture content in the period of constant drying and give:
Equation (5) shows the phase transition of water inside the fabric. Equations (6) and (7) show the phase transition of vapor and air respectively inside the fabric.
2.4. Boundary conditions for the period of falling drying
When the fabric moisture content falls to the critical point at
where
2.5. Calculation of mass fluxes
The mass flux of water in Equation (5) at the constant drying period is given as:
In Equation (11),
The generation of moisture is due to phase transition of water into vapor, in which, the efflux of vapor is significant. The coefficient
3. Fabric drying tests
A series of experiments were conducted to measure the drying characteristics of a group of fabric samples. Cotton is the major studying material in the tests as it is used most widely in clothing industry. The objectives of the experimental tests are to examine the drying characteristics of the fabrics under different boundary conditions, such as fabric texture, density, thickness, air temperature and impingement velocity.
3.1. Drying test set up
Six cotton fabric samples were examined. The samples were labeled from A to F, and their properties are listed in Table 1.
The set up as shown in Fig. 2 is an air heater providing hot air stream for each drying test. The temperature and speed of the impinging air are adjustable to provide different boundary conditions for the study. Disc-shaped fabric samples of 100 cm2 in area are mounted on a polystyrene backing plate with wire gauze facing the impinging hot air.
Plain knitted | 20 s/2 | 224.4 | 0.6594 | |
Plain knitted | 32 s/1 | 147.7 | 0.4363 | |
Plain knitted | 20 s/2 | 271.0 | 0.7769 | |
Plain weaved | - | 182.0 | 0.5638 | |
Plain knitted | 20 s/1 | 193.0 | 0.5025 | |
Plain knitted | 32 s/2 | 200.0 | 0.6188 |
In the tests, each fabric sample was being dried under eight conditions as listed in Table 2. The fabric weight was measured by an electronic balance every 30 seconds of drying under different setting conditions. The repeated drying and weight measurement procedures were conducted until all the moisture in the fabric samples was removed. The same testing procedures were repeated for all the fabric samples as given in Table 1.
80.0 | 1.48 | |
81.5 | 1.45 | |
86.5 | 1.43 | |
54.0 | 1.10 | |
55.5 | 1.15 | |
54.0 | 1.02 | |
57.0 | 1.41 | |
58.0 | 1.46 |
3.2. Results and discussions of the drying tests
Fig. 3 shows testing results from the six fabric samples under air setting condition 1 as listed in Table 2. The normalized water contents instead of the absolute values recorded from the tests were used in order to compensate the variation of fabric weight among the six samples.
The testing results as illustrated in Fig. 3 for the tested fabrics have shown the relationships among water content, fabric density, texture, air temperature and impingement velocity. Their relationships are discussed in the following sections.
3.2.1. Drying rate versus fabric texture
Table 3 lists the results of four fabric samples tested under the air setting conditions of 1, 4, 6 and 8 as listed in Table 2. Among the tested samples, sample D is weaved fabric and the others are knitted fabrics. The drying rate of sample D is the highest among the others at the constant drying period.
0.083 | 0.439 | 0.100 | 0.462 | |||
0.080 | 0.383 | 0.100 | 0.395 | |||
0.095 | 0.392 | 0.083 | 0.398 | |||
0.065 | 0.281 | 0.065 | 0.287 | |||
Average drying rate = 0.374 g/min | Average drying rate = 0.386 g/min | |||||
0.083 | 0.431 | 0.077 | 0.426 | |||
0.067 | 0.304 | 0.057 | 0.323 | |||
0.063 | 0.349 | 0.065 | 0.385 | |||
0.061 | 0.290 | 0.056 | 0.302 | |||
Average drying rate = 0.344 g/min | Average drying rate = 0.359 g/min |
3.2.2. Drying rate versus fabric density and thickness
Table 4 lists the testing results of two selected fabric samples A and C with same texture, yarn structure and different density and thickness under the air setting conditions of 1, 4, 6 and 8.
0.071 | 0.416 | 0.067 | 0.443 | |
0.063 | 0.407 | 0.063 | 0.391 | |
0.065 | 0.378 | 0.053 | 0.396 | |
0.047 | 0.281 | 0.044 | 0.289 | |
Average drying rate = 0.371 g/min | Average drying rate = 0.380 g/min |
The results listed in Table 4 show a similar result at the constant drying for the fabrics with different density and thickness, but same texture and yarn structure.
3.2.3. Drying rate versus air temperature
Table 5 lists the results of the drying rate for all fabric samples tested under air setting conditions 1 and 8 with similar impingement velocity at 1.47 m/s and different temperature.
Temp. (°C) | 1/time (min-1) | Drying rate (g/min) | 1/time (min-1) | Drying rate (g/min) | 1/time (min-1) | Drying rate (g/min) |
0.047 | 0.281 | 0.065 | 0.281 | 0.044 | 0.289 | |
0.071 | 0.416 | 0.083 | 0.439 | 0.067 | 0.443 | |
Temp. (°C) | 1/time (min-1) | Drying rate (g/min) | 1/time (min-1) | Drying rate (g/min) | 1/time (min-1) | Drying rate (g/min) |
0.065 | 0.287 | 0.061 | 0.290 | 0.056 | 0.302 | |
0.100 | 0.462 | 0.083 | 0.431 | 0.077 | 0.426 |
The drying rate at the constant drying period increases with the rise of air temperature for all fabric samples.
3.2.4. Summary of the experimental findings
The period of constant drying as illustrated in Figs. 1 and 3 has constituted a large portion of the drying cycle. The moisture reduction rate at the period could be used as an indicator to show the properties of the fabric, and conditions of the impinging air. The experimental findings in Tables 3 – 5 have shown the performance of the drying process against the boundary conditions including fabric texture, density, thickness, air temperature and impinging velocity. It has been observed that the increase of air temperature and velocity will speed up the drying rate. The fabric properties could also affect the drying rate but not as much as the air properties. These findings could be useful in the setting up of analytical models to simulate each period of a fabric drying cycle.
4. Development of non-linear analytical models to simulate the drying of porous type fabrics
As mentioned in Section 2, the drying rate of porous type fabrics has a non-linear relationship with time at the falling drying period. Some inaccurate results will be found if the traditional linear heat transfer equations are applied because the heat transfer coefficient changes with the change of the moisture contents at the falling drying period. To ensure an accurate modeling of the drying process, the heat transfer coefficient should be adjustable corresponding to the diffusion properties in the forming of dry/wet regions as mentioned in Section 2.4. A non-linear model is therefore used to describe the process characteristics (Haghi, 2006; Moropoulou, 2005). The knitted and weaved fabrics studied in this research are considered porous type materials because they contain unidirectional pores. The randomly distributed pores give an environment to establish a diffusion process when portions of water dry up. It is clear that the moisture diffusion rate has a close relationship to the size and number of fabric pores. The studied models given in Section 4.6 will address these essential modeling parameters. The other fabric parameters including texture, density and thickness that correlate to the drying rate will also be modeled by different modeling principles, and given in Sections 4.3 – 4.5.
4.1. Determination of the critical moisture content
The two periods of a fabric drying cycle as illustrated in Fig. 1 should be modeled separately. Traditional heat transfer equations could be used as modeling tools for the constant drying period, whilst, non-linear modeling equations should be considered when the moisture reduction rate varies with time in the falling drying period. The critical moisture content
The moisture reduction rate
The critical moisture content
Critical moisture content (g/g dry fabric weight) | |
0.8 | |
0.7 | |
0.8 | |
0.8 | |
0.7 | |
0.7 |
4.2. Using diffusion theories to model a fabric drying process
The boundary conditions for heat/mass transfer in the porous fabric have been discussed in Sections 2.2 to 2.5. However, they are not good enough to estimate the fabric moisture content during the drying process. It is necessary to have a further investigation to estimate the moisture content in individual period of drying. The authors have set-up a group of models based upon diffusion theories and Kowalski’s (2003) boundary equations to simulate the moisture changing rate under various boundary conditions (Ip and Wan, 2011). The investigated drying models will be presented in differential forms to address the movement of moisture contents in fabric. The models are based on the principles of chemical diffusion mechanism to calculate the rate of moisture change (
4.3. First order kinetics model
Roberts and Tong (2003) have shown a successful result in the modeling of bread drying process using first order exponential equations. In their research, microwave was used as the drying agent, and the process has been assumed as isothermal. Unfortunately, the experiential setup is quite different from convective drying using impinging air in this study. It is therefore necessary to develop new modeling equations for porous type fabrics. Schlunder (2004) has stressed that the falling drying period should be considered as an isothermal process. First order exponential equations might be appropriate to describe the process, thus, the first model developed in this study is labeled as “First order kinetics model”.
In the First order kinetics model, there is an assumption that the vaporization of water inside fabric can be described as a kinetic reaction motion of water molecules. The reaction rate is treated as the moisture reduction rate at the falling drying period. Thus, the water evaporation rate will correlate with the moisture content. The equation of the kinetic model is given as:
In Equation (17),
where
The testing results of moisture content as shown in Fig. 4 are further plotted in terms of drying cycle time
The kinetic coefficient
A graphical method to determine the kinetic coefficient
where
If the Arrhenius relationship is applied to describe a fabric drying process, a plotting of ln
353.0 | 1.48 | 0.5295 | 0.5198 | |
354.5 | 1.45 | 0.5818 | 0.5230 | |
359.5 | 1.43 | 0.5076 | 0.5336 | |
327.0 | 1.10 | 0.5297 | 0.4634 | |
328.5 | 1.15 | 0.5494 | 0.4667 | |
327.0 | 1.02 | 0.5537 | 0.4634 | |
330.0 | 1.41 | 0.3589 | 0.4699 | |
331.0 | 1.46 | 0.3568 | 0.4722 |
-0.6358 | 0.0028 | 0.3920 | |
-0.5416 | 0.0028 | 0.3716 | |
-0.6781 | 0.0028 | 0.3577 | |
-0.6354 | 0.0031 | 0.0953 | |
-0.5989 | 0.0030 | 0.1398 | |
-0.5911 | 0.0031 | 0.0198 | |
-1.0247 | 0.0030 | 0.3436 | |
-1.0306 | 0.0030 | 0.3784 |
Fig. 7 illustrates the plotting of ln
The calculated values for
Experimental results listed in Table 8 can be further used to determine the coefficients of
Using results from the regression table, the model equation is given as:
5.535 | -1974.7 | -1.521 | |
5.740 | -1993.9 | -1.210 | |
3.714 | -1476.9 | -0.911 | |
4.430 | -1584.6 | -1.276 | |
5.580 | -2054.5 | -0.494 | |
6.087 | -2223.7 | -0.280 |
Table 9 lists the regression results of all the fabric samples. A comparison of the differences of
80.0 | 1.48 | 0.5198 | 0.5195 | 1.83 | 1.90 | |
81.5 | 1.45 | 0.5230 | 0.5487 | 10.10 | 5.68 | |
86.5 | 1.43 | 0.5336 | 0.6056 | 5.12 | 19.30 | |
54.0 | 1.10 | 0.4634 | 0.5227 | 12.52 | 1.32 | |
55.5 | 1.15 | 0.4667 | 0.5022 | 15.06 | 8.59 | |
54.0 | 1.02 | 0.4634 | 0.5863 | 16.32 | 5.89 | |
57.0 | 1.41 | 0.4699 | 0.3786 | 30.95 | 5.50 | |
58.0 | 1.46 | 0.4722 | 0.3656 | 32.34 | 2.48 | |
0.4890 | 0.5037 | 15.53 | 6.333 |
The deviations of
The red curve illustrated in Fig. 9 is the modeled drying cycle obtained from the regression model for fabric sample A in the falling drying period with a kinetic coefficient
In conclusion, the falling drying period in a fabric drying cycle can be modeled by the First order kinetics model using an exponential function. A coefficient
224 | 0.6594 | 0.5037 | |
148 | 0.4363 | 0.6328 | |
271 | 0.7769 | 0.4136 | |
182 | 0.5638 | 0.5607 | |
193 | 0.5025 | 0.5469 | |
200 | 0.6188 | 0.5869 |
4.4. Diffusion model
Most of the Diffusion models presenting the change of moisture content have been based upon Fick’s law (Ramaswamy and Nieuwenhuijzen, 2002). The Fick's first law states the diffusion flux flowing from the regions of higher concentration to lower concentration obeying a magnitude proportional relationship to the concentration gradient. The one dimensional Fick’s first law in differential form is given as:
where
If the Fick’s second law is applied to model the process of drying porous fabric, the fabric will be considered as an infinite thin slab and dried from one direction. If heat transfer from the surrounding to the fabric is negligible, the integration result from Equation (25) will give the Diffusion model. The model equation given in Equation (26) is in terms of moisture content
The Diffusion equation is similar to Equation (18) of the First order kinetics model. The only difference between the two model equations is the fabric thickness
The slope of the fitted straight line in Fig. 6 is -0.5076 for fabric sample A with a thickness of 0.6594 mm.
Using information in Table 8 to determine the constants of
The regression results for all the fabric samples are listed in Table 12.
20.88 | -12070.75 | 2.173 | |
10.59 | -2006.16 | 1.215 | |
11.43 | -1495.8 | 0.9718 | |
11.374 | -1601.97 | 1.277 | |
10.578 | -1895.1 | 0.526 | |
10.08 | -2201.5 | 0.242 |
The effective diffusion coefficient
80.0 | 1.48 | 7.03023x10-7 | |
81.5 | 1.45 | 8.49438x10-7 | |
86.5 | 1.43 | 1.40575x10-6 | |
54.0 | 1.10 | 8.83487x10-8 | |
55.5 | 1.15 | 9.49414x10-8 | |
54.0 | 1.02 | 1.04100x10-7 | |
57.0 | 1.41 | 7.20528x10-8 | |
58.0 | 1.46 | 7.46037x10-8 | |
4.24032x10-7 |
Fig. 11 illustrates a comparison between the experiential records and the modeling results from regression model using
The Diffusion model has been applied to model each of the drying cycle for all fabric samples, the discrepancies between the modeling results and records from experiments are listed in Table 14.
76.63 | 42.86 | 49.03 | 44.44 | 48.35 | 43.34 |
4.5. Kinetics model based on the solutions of diffusion equations
Fick’s second law for diffusion applications is commonly used to simulate mass transfer process in convective drying. However, the exponential term in Equation (26) causes a restriction to the Diffusion model be applied in the falling drying period. A separate modeling process is needed to describe the constant drying period for completed modeling of a drying cycle. Efremov (1998, 2002) has proposed a mathematical solution to solve the Frick’s law using integral error functions:
where
where
The first and second terms in the right-hand-side of Equation (32) represent the characteristics in the constant drying period, and the third term represents the falling drying period. The new kinetics equation consists of two modeling sections to describe the linear and non-linear parts of a drying process. The drying rate at constant drying period
0.0069 | 0.0073 | 0.0074 | 0.0077 | 0.0072 | 0.0071 | |
0.0069 | 0.0079 | 0.0083 | 0.0080 | 0.0084 | 0.0071 | |
0.0065 | 0.0067 | 0.0069 | 0.0069 | 0.0070 | 0.0065 | |
0.0068 | 0.0064 | 0.0065 | 0.0066 | 0.0051 | 0.0054 | |
0.0066 | 0.0065 | 0.0072 | 0.0066 | 0.0074 | 0.0074 | |
0.0063 | 0.0065 | 0.0066 | 0.0066 | 0.0058 | 0.0064 | |
0.0047 | 0.0046 | 0.0051 | 0.0047 | 0.0048 | 0.0051 | |
0.0047 | 0.0047 | 0.0048 | 0.0048 | 0.0048 | 0.0050 | |
0.00618 | 0.00633 | 0.00660 | 0.00649 | 0.00631 | 0.00625 |
The drying rate at constant drying period
Tables 16 and 17 list the discrepancies of the two modeling results from Figs. 12and 13.
Discrepancy (%) | 41.17 | 30.25 | 40.53 | 30.21 | 36.63 | 38.06 |
Discrepancy (%) | 8.70 | 7.92 | 10.12 | 11.69 | 8.87 | 7.15 |
4.6. Wet surface model
The fourth analytical model, Wet surface, was proposed by Schlunder (1988, 2004). He has proved that the remaining moisture in porous materials would be 20 to 30 % of the saturated moisture content at the initial stage when a drying process reaches the critical moisture stage. The material surface is unlikely to be fully wetted at this low moisture content condition. The Wet surface model is therefore designed to address the characteristics of the partially wet surface. The drying rate
The results of converting Equation (36) from the critical moisture content
In Equation (37), the parameters for fabric drying modeling are the initial drying rate
where
No (g/s) | Th (K) | TSL (K) | (J/g) | U (W/m2K) | c (W/mK) | (mm) | |
0.0069 | 353.0 | 301 | 2310 | 30.79 | 0.028 | 0.91 | |
0.0077 | 354.5 | 301 | 2304 | 33.29 | 0.028 | 0.84 | |
0.0065 | 359.5 | 302 | 2291 | 25.70 | 0.028 | 1.09 | |
0.0068 | 327.0 | 293 | 2372 | 47.35 | 0.028 | 0.59 | |
0.0066 | 328.5 | 294 | 2372 | 45.32 | 0.028 | 0.62 | |
0.0063 | 327.0 | 293 | 2372 | 43.89 | 0.028 | 0.64 | |
0.0047 | 330.0 | 294 | 2366 | 31.04 | 0.028 | 0.90 | |
0.0047 | 331.0 | 295 | 2366 | 30.75 | 0.028 | 0.91 | |
0.8125 |
The final parameter to be determined for Equation (37) is the fabric pore size
The fitted line in Fig. 14(a) shows the results from an assigned pore size of 0.5 mm. Results from the plotting have shown that the slope of the fitted line is not a unity, and the y-intercept does not meet the origin. As a result, the calculated viscous sub-layer thickness
Table 19 lists the determined pore sizes for other fabric samples. The calculated pore sizes for each fabric sample are then substituted into Equation (37) to determine the drying rate
0.40 | 0.61 | 0.19 | 0.31 | 0.46 | 0.36 | |
0.40 | 0.55 | 0.35 | 0.70 | 0.40 | 0.45 | |
0.35 | 0.40 | 0.55 | 0.31 | 0.21 | 1.55 | |
0.16 | 0.23 | 0.22 | 0.25 | 0.20 | 0.50 | |
0.19 | 0.31 | 0.22 | 0.40 | 0.37 | 0.35 | |
0.33 | 0.07 | 0.49 | 0.40 | 0.37 | 0.44 | |
0.40 | 1.10 | 0.60 | 0.31 | 0.43 | 0.55 | |
0.55 | 0.53 | 0.49 | 1.50 | 1.10 | 1.20 | |
0.3478 | 0.4750 | 0.3888 | 0.5225 | 0.4425 | 0.6750 |
The findings given in Fig. 15 have shown that the Wet surface model could not produce an accurate modeling result as the measured discrepancy is 77.53 % based upon the experimental records.
5. A performance evaluation of the studied models
The performance of the studied analytical models for the modeling of porous fabric drying process should be reviewed. The percentage of discrepancies from each modeling results in comparison with the experimental records are summarized in Table 20. It is made clear that the Kinetics model based on the solutions of diffusion equations has produced the best performance. The First order kinetics model has provided a better performance than the Diffusion model, and the Wet surface model has given the largest discrepancy from the statistical records listed in Table 20.
42.39 | 34.69 | 38.74 | 34.21 | 42.25 | 37.11 | 38.23 | |
76.63 | 42.86 | 49.03 | 44.44 | 48.35 | 43.34 | 50.77 | |
8.70 | 7.92 | 10.12 | 11.69 | 8.87 | 7.15 | 9.08 | |
77.53 | 37.47 | 93.50 | 56.00 | 69.87 | 64.02 | 66.40 |
The findings have shown that the Kinetics model based on the solutions of diffusion equations could be the best one in the simulation of a porous fabric drying process among the others. The required condition for the model is to have a predictable drying cycle time
6. Conclusion
The principles of water mass movement due to phase change in the drying of porous fabrics have been studied. The boundary equations for mass transfer between water, vapor and air were used to support the establishing a new set of drying models using diffusion theories. Experiments were done to find information for the determination of the modeling parameters. The performance of the developed models has been evaluated. Among the four models, the Kinetics model based on the solutions of diffusion equations has produced the best performance. In the real life applications, they could act as a mathematical tool to assist a precise estimation of the moisture content in fabric drying or heat setting process under various processing conditions. Further work has been started to apply the developed drying models in the design of garment setting machines for clothing industry.
References
- 1.
Cohen J. Cohen P. West S. G. Aiken L. S. 2003 (2nd ed.): Lawrence Erlbaum Associates, Hillsdale, NJ - 2.
Efremov G. I. 1998 Kinetics of Convective Drying of Fibre Materials Based on Solution of a Diffusion Equation. ,30 6 417 422 - 3.
Efremov G. I. 2002 Drying Kinetics Derived from Diffusion Equation with Flux-Type Boundary Conditions. ,20 1 55 66 - 4.
Haghi A. K. 2006 Transport Phenomena in Porous Media: A Review. .40 1 14 26 - 5.
Ip R. W. L. Wan I. C. 2011 New Use Heat Transfer Theories for the Design of Heat Setting Machines for Precise Post-Treatment of Dyed Fabrics. , Vols. 312-315,748 751 - 6.
Khazaei J. Chegini G. R. Bakhshiani M. 2008 A Novel Alternative Method for Modeling the Effects of Air Temperature and Slice Thickness on Quality and Drying Kinetics of Tomato Slices: Superposition Technique. ,26 6 759 775 - 7.
Kowalski S. J. 2000 Toward a Thermodynamics and Mechanics of Drying Processes. ,55 7 1289 1304 . - 8.
Kowalski S. J. 2003 , Springer,3-54000-412-2 Germany - 9.
Kowalski S. J. Musielak G. Banaszak J. 2007 Experimental Validation of the Heat and Mass Transfer Model for Convective Drying. ,25 1-3 ,107 121 - 10.
Moropoulou A. 2005 Drying Kinetics of Some Building Materials. ,22 2 203 208 - 11.
Ramaswamy H. S. van Nieuwenhuijzen N. H. 2002 Evaluation and Modeling of Two-Stage Osmo-Convective Drying of Apple Slices. ,20 3 651 667 - 12.
Roberts J. S. Srikiatden J. 2005 Moisture Loss Kinetics of Apple During Convective Hot Air and Isothermal Drying. ,8 3 493 512 - 13.
Roberts J. S. Tong C. H. 2003 Drying Kinetics of Hygroscopic Porous Materials under Isothermal Conditions and the Use of a First-Order Reaction Kinetic Model for Predicting Drying.6 3 355 367 . - 14.
Schlunder E. U. 1988 On the Mechanism of the Constant Drying Rate Period and Its Relevance to Diffusion Controlled Catalytic Gas-Phase Reactions. ,43 10 2685 2688 - 15.
Schlunder E. U. 2004 Drying of Porous Material during the Constant and the Falling Rate Period: A Critical Review of Existing Hypotheses. ,22 6 1517 1532