Next Article in Journal
PMLPNet: Classifying Multi-Class Pests in Wild Environment via a Novel Convolutional Neural Network
Previous Article in Journal
Inducing Aquilaria sinensis (Lour.) Spreng “Qinan” Agarwood Formation with Chemical Reagents
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

The Analytical Theory of Water Activity in the Paddy Rice Drying Process

1
College of Engineering, South China Agricultural University, Guangzhou 510642, China
2
School of Intelligent Engineering, Shaoguan University, Shaoguan 512099, China
3
College of Agricultural Engineering, Al-Azhar University, Cairo 11751, Egypt
*
Author to whom correspondence should be addressed.
Agronomy 2024, 14(8), 1728; https://doi.org/10.3390/agronomy14081728
Submission received: 15 May 2024 / Revised: 12 July 2024 / Accepted: 5 August 2024 / Published: 6 August 2024
(This article belongs to the Section Agricultural Biosystem and Biological Engineering)

Abstract

:
Drying involves the evaporation of moisture, accompanied by simultaneous heat transfer, mass transfer, and momentum transfer. While the diffusion law is considered an applicable model for explaining the drying phenomenon, the actual drying process cannot be accurately predicted using an analytical solution with a constant diffusion coefficient. Energy efficiency in the drying process is low due to an insufficient understanding of the mechanisms governing moisture migration from solids to air. The development of drying theory has stalled due to an unsolvable discrepancy between experimental results and analytical results. This study analyzes the effect of the binding energy of moisture in paddy rice on the diffusion coefficient. The theoretical relationship between water activity and drying rate in paddy rice was investigated, and the drying process was successfully explained by analyzing free energy transfer and transition theory. The mechanism of rice grain drying was described using a new theoretical solution for the drying process. These results provide new insights into the development of a scientific evaluation standard for assessing the efficiency of actual drying processes.

1. Introduction

The transfer laws of heat, mass, and momentum regulate the drying phenomenon. However, each of these transfer laws describes a specific behavior within the phenomenon. When applied, these conditions must satisfy not only their respective singularity criteria but also geometric, physical, temporal, and boundary conditions. The simultaneous transfer of heat, mass, and momentum during the drying process makes it challenging to derive a mathematically precise solution for the transfer coefficients [1]. The diffusion models yield significant errors when the diffusion coefficient is treated as a constant, resulting in discrepancies compared to actual conditions. Furthermore, drying mechanism models based on diffusion kinetics face challenges in accurately quantifying the activation energy, pre-exponential factor, and process exponent [2,3].
Researchers have modeled moisture migration as continuous medium flow, using time and space continuous functions and applying Fick’s law and the Fourier heat conduction equation. They have employed the Eulerian method, Lagrangian method, and volume averaging method [4,5], establishing numerous drying models tailored to specific material systems [6,7,8,9]. Researchers have differentiated mass transfer within materials into capillary flow and vapor diffusion motion in accordance with the laws of mass and energy conservation and the principles of irreversible thermodynamics. Various drying differential equation sets have been established based on temperature, concentration, and pressure potential fields [9,10,11,12,13]. Using computer simulation, researchers developed pore network models to analyze the drying process and investigated calculation parameters for a series of models [14,15]. Researchers have observed that materials undergo a slowed drying process under conditions where the evaporation surface recedes inward from the outer surface and the area gradually decreases [16,17,18]. Materials are simplified by considering the vaporization interface as the boundary, with a wet zone inside the interface and a dry zone outside, assumed to be in complete equilibrium. Various diffusion models have been developed to analyze specific material drying processes. As the drying progresses, real-time changes in the material’s physical properties lead to complex physicochemical and biochemical reactions, volume shrinkage, texture, and configuration alterations, which in turn affect pore structure, tortuosity, and moisture transport pathways, influencing overall intensity. This results in discrepancies between the actual process and the idealized diffusion at the interface [19,20,21,22,23]. The drying process does not conform strictly to coefficient diffusion, revealing deficiencies in existing drying models and incomplete theoretical descriptions [24].
During actual drying processes, moisture and temperature vary between the inside and outside of rice grains. The state changes in moisture vaporization and migration processes within rice grains demonstrate significant inertia and nonlinear characteristics, differing from ideal gas behavior [25,26]. This fundamental difference distinguishes rice drying from the quantitative material systems and ideal processes described in general thermodynamics [27,28].
Moisture vaporization and migration are spontaneous processes, making water activity an intrinsic property of rice. Experimental data under steady-state conditions allow for an accurate correlation between water activity and external condition parameters. Moisture migration is influenced by free energy [28,29], which represents the energy available in a material system to perform technical work in the external environment. This phenomenon can be observed through molecular or microscopic motion or described through heat and mass transfer phenomena during the process [29,30,31,32]. Following this research approach, researchers have derived the theoretical expression for grain moisture binding energy [27] and provided mathematical solutions for the vaporization latent heat coefficient and drying characteristic function of grain moisture. By establishing the theoretical link between rice water activity and drying rate and integrating the transfer and conversion of free energy as a unified metric [28,33], this study explores the binding energy of rice moisture. It presents a theoretical framework for understanding the rice drying process and outlines practical applications.
This study aims to establish a rigorous mathematical relationship between the drying rate of rice grains and the consumption of free energy during moisture evaporation while also revealing the theoretical solution for the water-binding energy of rice grains. Therefore, the theoretical heat consumption of rice grain drying can be evaluated under various environmental conditions. In practical applications, this study not only analyzes the drying process of rice grains at specific phases but also establishes a theoretical framework for predicting the drying process under a variety of external conditions through these analyses. These findings provide a theoretical basis for objectively evaluating the drying efficiency of rice grains in any actual drying process and for developing scientific standards to evaluate energy utilization efficiency in drying equipment systems.

2. The Paddy Rice Drying Theory

2.1. The Moisture Content in Paddy Rice

In the field of drying technology, moisture-binding in materials is categorized into chemically bound water, physicochemically bound water, and physical–mechanically bound water. Chemically bound water showed strict quantitative relationships, whereas physicochemical-bound water, which lacks such strict quantification, can be further categorized into (1) adsorption-bound water, a liquid confined by force fields on external and internal surfaces of “capsules”, and (2) osmotically maintained water, encompassing expansion water and structural water, sealed within cells. This category includes water adsorbed by composite capsules through osmotic absorption and fixed structural water. Free water, due to its low binding energy, is classified separately from physical–mechanically bound water, which can vary in quantity and is present in both large and small capillaries within materials. In drying engineering, regardless of the form in which moisture exists in the material, the primary focus is on the amount of moisture that can be removed from the wet material and the associated energy consumption. Based on equilibrium characteristics, moisture can be categorized as bound moisture or free moisture, depending on the relative humidity of the air and the moisture content of the material, as illustrated in Figure 1a,b.
The activity of non-bound water is equal to 1, and its evaporation was unaffected by the material, similar to the evaporation of a free liquid surface. The rate of vaporization remained solely dependent on external conditions. Bound moisture refers to the moisture contained within rice grains, where, under steady-state conditions, evaporation characteristics are primarily governed by the physical properties of the rice. The evaporation rate corresponds to the prevailing environmental conditions and is influenced by water activity [34].

2.2. Two-Staged Characteristics of Paddy Drying

Paddy rice comprises two parts: the husk and the grain, which differ significantly in texture. The husk, with its loose texture, surrounds the grain and facilitates easier drying. Therefore, the drying process of rice is characterized by a two-stage decrease in drying rate [35]. During the initial stage of drying, moisture distribution inside and outside the rice was uniform and consistent, with moisture exhibiting zero binding energy. At this stage, the drying rate was at its maximum, while with the progress of the drying process, the husk first encountered the external medium, gradually disrupting the thermal equilibrium between the husk and the grain. This phenomenon led to a significant decrease in the drying rate during the initial stage (Stage I) of drying. As drying progresses and the moisture content of the husk decreases rapidly, the drying rate subsequently slows down, primarily influenced by internal moisture diffusion within the grains, leading to a two-stage decrease in drying rate, as illustrated in Figure 2.

3. Theory Analytic Methods

3.1. The Free Energy Consumption in the Moisture Evaporation

The drying process can be quantitatively explained by considering both the amount of moisture migration and its migration potential, determined by the material system’s state parameters [36]. Furthermore, all processes dependent on thermodynamic state parameters can be attributed to enthalpy. Understanding the energy input into the drying system allows us to describe changes in both the system and its materials using state functions, regardless of the form of water present or alterations in the drying mechanism parameter. Therefore, the enthalpy of any given drying state point may be used to determine the thermal energy consumption at that state point, regardless of the conditions. This calculation facilitates the derivation of the expression for the material’s internal moisture binding energy, represented as Equation (1) [27].
μ = p s p v d p = R p n T ln φ
where ϕ is the water activity, ϕ = p p s , P represents the vapor fraction pressure, Pa; T is the vaporization temperature of water, K; p s is the saturated vapor pressure, Pa; R p n is the interval characteristic constant R p n = p s p c p d ( ln p ) p p s ; c p is the specific heat of water vapor at constant pressure, kJ/kg, with a fixed value at each state point. In [ p s , p ] interval, isothermal test condition, R p n = i = 1 n c p i n , kJ/kg·K−1; c p i is the specific heat value, kJ/kg, corresponding to different pressures under isothermal conditions.
Similarly, it is feasible to experimentally ascertain the specific heat at constant pressure at varying temperatures under equilibrium pressure conditions. The physical significance of Rpn was almost similar to the average specific heat at constant pressure throughout the state changes occurring during the process of moisture vaporization.
The free energy expended during moisture vaporization reflects the ability of water vapor to perform external work. Letting f denote the specific free energy consumption of water vaporization in its free state (measured in kJ/kg), the rate of change in the system’s free energy is given by d f = p d v .
Because the specific volume of water vapor is much larger than that of liquid water at the same temperature, the specific volume of water can be considered negligible in comparison. Thus, when the water activity equals 1, the integral expression for the free energy consumption required to vaporize 1 kg of moisture in the material system can be expressed as Equation (2).
0 f d f = 0 v p d v
Due to d v = c p p 2 T d p , At the same temperature, p corresponds to v one by one. Thus, expression (2) can be rewritten as: 0 f d f = 0 p c p p T d p . The specific free energy for vaporizing 1 kg of moisture was obtained through an integral calculation, which can be written as.
f = R p n T
The free energy required per 1 kg of vaporized water is the sum of the terms presented in Formulas (1) and (3) upon moisture vaporization from the material. The expression for the free energy consumed per 1 kg of water vaporized from the material is represented as Equation (4), using the symbol fg.
f g = R p n T + μ = R p n T ( 1 ln ϕ )
Equation (4) indicates that the free energy consumed during moisture evaporation is a state function dependent on absolute temperature and water activity. Equation (4) applies to any material and vaporization process, holding universal significance. All specific properties associated with the internal structure and processes of the substance can be confirmed as real interval characteristic parameters Rpn. Concerning the microstructure and process variables of the material, no assumptions were made [36]. Using internationally recognized steam tables and experimental data within the temperature range of 0 ≤ T −273.15 ≤ 100 °C at one-atmosphere pressure, the constant characteristic formula for the state change interval of water vapor in [273.15, T] is derived as Equation (5).
R p n = 2 × 10 9 t 4 4 × 10 7 t 3 + 4 × 10 5 t 2 0.0015 t + 2.3754 t = T 273.15
Expression (6) was derived by dividing Equation (4) by Equation (3), where the 1 ln ϕ term denoted by the symbol Cz represents the specific free energy consumption coefficient.
C z = f g f = 1 ln ϕ
Expression (6) demonstrates that water activity equals 1 in a free state, with the specific free energy consumption coefficient also equaling 1. In the presence of water-binding energy, water activity was less than 1, resulting in a free energy ratio greater than 1. When vaporized from the material, the same amount of water consumed 1 ln ϕ times heat compared to its free state.

3.2. The Free Energy Balance Equation for the Evaporation Process

Regardless of the circumstances, the vaporization of an equivalent amount of moisture from any material consistently required the same quantity of free energy. Assuming the material’s moisture content is M0 (% d.b.) at a water activity level of 1, the equilibrium moisture content is represented by Me (% d.b.). The maximum drying rate (% d.b.h−1) and drying rate (% d.b.h−1) within the drying interval M (MeMM0) can be represented by symbols ( d M d θ ) max and d M d θ , respectively. The reduction d M d θ relative to the ( d M d θ ) max was found to correspond to the additional free energy consumption required for the evaporation of water constrained within the material, leading to the derivation of equilibrium Equation (7) that describes the relationship between free energy consumption. Simplifying Equation (4) provided the inherent drying rate ratio state function expression (8) of the material.
f d M d θ + μ d M d θ = f ( d M d θ ) max
d M d θ ( d M d θ ) max = f f + μ
where f = R p n T is the free energy of 1 kg water vapor with 1 equivalent of water activity, kJ/kg; θ denotes drying time, and M represents the dry base moisture content of the material, % d.b.
By substituting Formulas (1) and (3) into Formula (8), an expression (9) was obtained for the drying rate ratio state function.
d M d θ ( d M d θ ) max = 1 1 ln φ ,   0 < max 1
Equations (7)–(9), d M d θ , and ( d M d θ ) max represent the theoretical drying rate and maximum drying rate inherent to the material itself, reflecting the drying rates that can theoretically be achieved by the material. The 1 1 ln φ in Equation (9) represents the theoretical drying rate ratio of the material in an ideal, perfectly dry medium under steady-state conditions, often referred to as the theoretical drying rate ratio.
Equation (9) reveals that the theoretical drying rate ratio is a single-valued function of the moisture content of the material. Under isothermal conditions, φ is a single-valued function of M, a mass ratio expressed in % d.b., and remains unaffected by the moisture distribution within the material. Therefore, Equation (9) can be used to express the drying rate ratio for any drying process and the corresponding moisture content state.

3.3. Characteristics of Theoretical Drying Rate Ratio

The theoretical drying rate ratio is a single-valued function of water activity. Therefore, the theoretical drying rate ratio can be determined by Equation (9), as illustrated in Figure 3, under steady-state drying conditions as a function of water activity within the range of [1, 0).
Water activity indicates the extent of moisture binding and resistance to moisture migration, which affects the consumption of free energy during the moisture migration process and is influenced by the dry matter components’ resistance.
As depicted in Figure 4, higher water activity values correspond to greater drying rate ratios, indicating a lower degree of binding and smaller binding energy between the moisture and dry matter components. Conversely, lower water activity values correspond to higher degrees of binding and larger binding energies between the moisture and dry matter components. As water activity approaches zero, the drying rate ratio similarly approaches zero. Throughout this process, even slight fluctuations in water activity can lead to notable variations in the drying rate ratio.

3.4. The Water Activity Model for Paddy Rice

Water activity is influenced by the moisture content and temperature of the material and serves as an inherent characteristic function of the material itself. Under steady-state conditions, the equilibrium moisture content represents the maximum moisture level that the material can attain during the drying process. It signifies the moisture content at which the material fully returns to its environmental state.
Various equations have been proposed by different research institutions and scholars across different countries to calculate the equilibrium moisture content for various types of materials. These equations are collectively referred to as equilibrium moisture content models [37,38,39].
Since equilibrium moisture content can be accurately determined under static experimental conditions, several precise formulas have been developed for various materials. Equation (10) is among the widely recognized equilibrium moisture content models used internationally, particularly for grains [40].
M e = ln ( 1 ϕ e ) A ( T a 273.15 + B ) 1 n
where T a represents the temperature of water vaporization, K. M e is the dry base equilibrium moisture content (%); A, B, and n represent the calculation parameters. For paddy, A = 1.9187 × 10−5, B = 51.161, n = 2.4451 [28].
Equilibrium moisture content is a state function and can theoretically be expressed as M e = f T , φ e , where subscripts indicate constant temperature and humidity drying conditions. Based on Equation (10), the expression for the equilibrium relative humidity φ e of the drying medium can be derived, as shown in Equation (11).
ϕ e = 1 exp M e n A ( T a 273.15 + B )
Under phase equilibrium conditions, the migration potential at each point within the material system was uniform, resulting in a zero mass migration gradient. The temperature and pressure inside and outside the material reach equilibrium, and all state parameters remain constant over time. Therefore, under phase equilibrium conditions, the water activity is numerically equal to the relative humidity of the medium, and it can be calculated based on the temperature and relative humidity of the medium. The drying phenomenon can be regarded as a transition from one equilibrium state to another when employing the general method of thermodynamic analysis. Therefore, Equation (12) can be used to represent the grain’s water activity during the drying process.
ϕ = 1 exp M n A ( T g 273.15 + B )
where T g is the material temperature, K; ϕ is the water activity; M represents the moisture content, % d.b.
The moisture content and temperature of the material were the two independent variables in Equations (11) and (12). These two parameters are both state parameters of the drying material system; however, the physical meanings of Equations (11) and (12) are distinct. The water activity of the rice grain is an intrinsic property of the grain and acts as a state function dependent on the grain’s temperature and moisture content. In contrast, the relative humidity of the surrounding medium is a state function determined by the temperature and moisture content of that medium.

4. Theoretical Analysis and Discussion

The free energy consumption ratio refers to the energy, expressed in kJ/kg, consumed by the material system to facilitate the movement of water vapor and exert external forces when 1 kg of water evaporates from the material. Energy consumption refers to the internal energy expended by the material system to perform this function.
Theoretically, it denotes the maximum technical work achievable after the evaporation of 1 kg of water from the material. During the drying process, the intrinsic energy of the rice grain is transformed into kinetic energy, which is necessary for the movement of water vapor [27].
Figure 4a illustrates the relationship between the free energy consumption ratio of rice grain during moisture evaporation and moisture content at various temperature conditions, as determined by Equations (4), (6), (9) and (12). Figure 4b illustrates the variation in the free energy coefficient with moisture content. Figure 5a illustrates the relationship between moisture content and water activity, while Figure 5b illustrates the theoretical drying rate ratio as determined from the analysis.
Based on Figure 4a, it was evident that the free energy of rice grain decreased as the moisture content increased during slow drying. The impact of moisture content was more pronounced than that of temperature, especially at lower moisture levels. However, when the moisture content exceeded 25% d.b., changes in free energy were more closely associated with temperature, leading to an increase in values as temperature increased, approaching a constant state.
At this stage, the change in free energy related to moisture in the material can be described as a single-valued function of temperature, and the moisture evaporation process is equivalent to free surface evaporation. Therefore, it can be inferred that the drying process of the material with high moisture content can be simplified as constant-rate drying. The curves in Figure 4a intersect at varying temperatures, suggesting that the moisture evaporation rate in the rice hull was higher than that of the rice grain during the initial drying stage. As the temperature increased, the moisture content in the rice hull decreased, leading to a corresponding increase in the consumption of free energy by the rice hull.
The consumption ratio of free energy is a universal property that reflects all changes in drying material states. According to Equation (4), when the water activity was equal to 1, it demonstrated the 1 equivalent free energy consumption. However, if there was water binding energy and the water activity dropped below 1, the free energy consumption increased beyond 1.
In any drying process system, there exists a maximum drying rate state point, which corresponds to the highest water activity achievable under the same environmental conditions. As per Equation (4), the free energy required to convert moisture into vapor depends on the temperature and water activity. The change in free energy consumption during the drying process is solely determined by the water activity, assuming the temperature remains constant. Therefore, the free energy coefficient can be determined at any drying temperature by using Equation (12), which establishes a correlation between moisture content, temperature, and water activity. Figure 4b shows the variation in the free energy coefficient of rice grain during moisture evaporation under different temperature conditions. It is evident from Figure 4b that the free energy coefficient of rice grain during moisture evaporation increased as the moisture content decreased, regardless of the temperature. Furthermore, the free energy coefficient decreases with increasing temperature at a given moisture content. This indicated that as the temperature rose, the drying driving force increased, thereby reducing the free energy consumption during moisture evaporation.
In this study, we examined the moisture content variations over time in fresh rice grains, initially set at 40% moisture, by utilizing a thin-layer dryer within various drying environments. The experimental setup included monitoring moisture levels, specifically in the central area of the drying chamber. As depicted in Figure 4c, our findings reveal that elevated temperatures significantly enhance the initial reduction in moisture content within the rice grains. However, during the later stages of the drying process, a marked deceleration in the rate of moisture loss was observed. This phenomenon can be attributed to the increased consumption of free energy when the rice grains reach lower moisture levels. The experimental outcomes align well with the predictions of our theoretical model, thereby validating the model’s accuracy in simulating the drying process of rice grains.
Water activity is also dimensionless, with a range of (0, 1). It is evident from Figure 5a that the moisture content change varied most significantly at 0 and progressively decreased toward the maximum value at the same water activity under different temperature conditions. Similarly, the change in water activity difference under varying temperature conditions at the same moisture content followed a pattern of initiating at 0, reaching a maximum, and then declining toward 0. When the moisture content of rice grain exceeded 30% d.b., the water activity of the material closely approached that of a free state. However, when the moisture content was below 5%, the water activity was very low and gradually decreased with decreasing moisture content. This indicated that the binding force of moisture to the dry matter components of rice grain rapidly increased with decreasing moisture content, following a nonlinear process.
In two geometrically similar drying systems, if they share the same material water activity or the difference between material water activity and the relative humidity of the medium, their drying rates relative to the maximum drying rate under those conditions should be identical. A lower drying rate ratio indicates greater binding energy between the material and moisture, reflecting an inherent property of the material itself. Therefore, the drying rate ratio expresses the intrinsic drying mechanism of the material itself, independent of dimensions, and comprehensively describes the principles of material drying, revealing its fundamental nature. Therefore, regardless of the distribution of water within the material or changes in other parameters of the drying mechanism, the drying phenomenon can be quantified using the state function of the drying rate ratio. Figure 5b illustrates the variation in the rice grain drying rate ratio with moisture content under different temperature conditions. It is apparent that at the same theoretical drying rate ratio level, the difference in moisture content change under various temperature conditions followed a pattern of starting at 0, reaching a peak, and then decreasing back towards 0.
Similarly, at identical moisture content, the variation in water activity difference under different temperature conditions followed a pattern of starting at 0, reaching a peak, and then decreasing towards 0. When the moisture content of rice grain surpassed 30% d.b., the theoretical drying rate ratio showed minimal variation. As the drying progressed into the falling rate stage with water activity below 1, the theoretical drying rate ratio decreased rapidly with decreasing moisture content. Moreover, at the same moisture content, the higher temperature led to a higher drying rate ratio. However, the influence of temperature on the drying rate ratio was less significant when the moisture content was reduced to about 5% d.b. The differences in both vertical and horizontal dimensions between the curves in Figure 5b distinctly illustrate that the drying process of rice grain was influenced by both temperature and moisture content concurrently. The influence of temperature and moisture content on the drying rate ratio exhibits a pattern of increase followed by decrease. These analytical findings offered a novel theoretical foundation and method for strategically adjusting drying condition parameters based on the material’s drying state, aiming to achieve efficient and energy-saving drying processes.
It is evident from the comparison of Figure 4a,b and Figure 5a,b that the free energy consumption of the material system was less affected by temperature and moisture content when the moisture content of rice grain exceeded 25% than by the free energy coefficient, water activity, and drying rate ratio. In a high-moisture state, the material has a water activity close to 1. This high water activity provides a strong driving force for water removal, making the drying process less dependent on temperature increases. During the falling rate drying stage, with water activity below 1 (for example, in the moisture content range of 25% d.b. to 10% d.b., as shown in Figure 5b), temperature and moisture content had a larger impact on free energy consumption compared to the free energy coefficient, water activity, and drying rate ratio. Moreover, all these factors displayed nonlinear characteristics. At this stage, the drying process can be strengthened by increasing the drying temperature, which can considerably increase the free energy, water activity, and drying rate ratio while simultaneously reducing the free energy coefficient.

5. Conclusions

This work presents a theoretical solution for the amount of energy required to bind moisture within rice grains. It also demonstrates a strong mathematical relationship between the rate at which rice grains dry and the amount of energy consumed during the evaporation of moisture at different levels of water activity. Consequently, a theoretical assessment of the thermal consumption associated with rice grain drying can be derived under a variety of environmental conditions. In practical applications, this study not only offers an analytical solution for the drying process of rice grains under specific drying conditions but also establishes a theoretical framework to understand how rice grains dry under various external conditions through analytical methods. This study provides a theoretical foundation for future research into objectively evaluating the drying efficiency of rice grains in practical drying operations and for the development of efficiency evaluation standards for energy utilization in scientifically engineered drying equipment systems.
This study revealed the mathematical expression for the free energy ratio consumption during moisture vaporization in the dry matter system. It obtained the specific free energy consumption coefficient and derived the intrinsic characteristic function of the drying rate ratio of paddy. The theory reveals the following facts:
  • The drying process can be quantitatively described by the transfer and conversion of free energy, regardless of the form of water within the material or changes in mechanical parameters.
  • The intrinsic theoretical drying rate ratio of the material itself is a single-valued function of the material’s water activity, with its mathematical expression as follows: d M d θ ( d M d θ ) max = 1 1 ln ϕ . Under isothermal conditions, water activity is a single-valued function of water content, where water content can be expressed as a mass ratio.
  • Reduced water activity results in an increase in the free energy consumption of water evaporation in paddy, which can be expressed mathematically as   f g = R p n T 1 ln φ .
  • The coefficient of specific free energy consumption was found to be equal to 1 when the water activity was 1 in a free state. Furthermore, it was observed that when binding energy between moisture and dry matter existed, the water activity was less than 1. The specific free energy consumption coefficient was found to be greater than 1, indicating that the same amount of water vaporized from the paddy consumed C z = 1 ln ϕ time heat compared to the free state. Its mathematical expression is C z = 1 ln ϕ .

Author Contributions

Conceptualization, C.L. and H.W.; methodology, C.L.; software, H.W.; validation, H.W., Y.Z. and K.A.M.A.; formal analysis, K.A.M.A.; investigation, H.W. and S.L.; resources, C.L.; data curation, H.W. and Y.Z.; writing—original draft preparation, C.L.; writing—review and editing, H.W.; visualization, K.A.M.A. and Y.Z.; supervision, K.A.M.A.; project administration, C.L. and W.W.; funding acquisition, C.L. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Natural Science Foundation of China, grant number 32171906.

Data Availability Statement

The original contributions presented in the study are included in the article, further inquiries can be directed to the corresponding author.

Conflicts of Interest

The authors declare no conflicts of interest.

Nomenclature

M moisture content on a dry basis, % d.b.
M e equilibrium moisture content, % d.b.
M c initial moisture content, % d.b.
M e l equilibrium moisture content in the first falling-rate drying stage, d.b.
M k limiting moisture content, % d.b.
φ water activity, decimal
P water vapor partial pressure, KPa
P s saturated vapor pressure, KPa
R p n characteristic constant of the water evaporation process, kJ/(kg·K)
C p specific heat capacity of water vapor at constant pressure, kJ/kg
f specific free energy consumption during the vaporization of free-state water, kJ/kg
C p i specific heat values at different pressures under isothermal conditions, kJ/kg
C z specific free energy consumption coefficient
T   moisture vaporization temperature, K
t Kelvin temperature, K
θ   drying time, hours

References

  1. Crank, J. The Mathematics of Diffusion; Clarendon Press: Oxford, UK, 1990. [Google Scholar]
  2. Le, K.H.; Tran, T.T.H.; Kharaghani, A.; Tsotsas, E. Experimental benchmarking of diffusion and reduced models for convective drying of single rice grains. Dry. Technol. 2020, 38, 200–210. [Google Scholar] [CrossRef]
  3. Wu, J.; Zhang, H.; Li, F. A study on drying models and internal stresses of the rice kernel during infrared drying. Dry. Technol. 2017, 35, 680–688. [Google Scholar] [CrossRef]
  4. Zhang, S.; Liu, N.; Pan, Y.; Wang, W.; Li, Y.; Zhu, Y. Three-dimensional modelling of two-phase flow and transport in a pilot centrifugal spray dryer. Chem. Phys. Lett. 2021, 765, 138309. [Google Scholar] [CrossRef]
  5. Whitaker, S. The Method of Volume Averaging; Springer: Berlin/Heidelberg, Germany, 1998. [Google Scholar]
  6. Whitaker, S. Simultaneous Heat, Mass, and Momentum Transfer in Porous Media: A Theory of Drying. In Advances in Heat Transfer; Hartnett, J.P., Irvine, T.F., Eds.; Elsevier: Amsterdam, The Netherlands, 1977; Volume 13, pp. 119–203. [Google Scholar]
  7. Arrieche, L.S.; Sartori, D.J.M. Fluid Flow Effect and Mechanical Interactions during Drying of a Deformable Food Model. Dry. Technol. 2007, 26, 54–63. [Google Scholar] [CrossRef]
  8. Roques, M.A. From Drying to Thermo-Hydro-Rheology—A Thirty-Year-Long Unfinished Maze. Dry. Technol. 2008, 26, 1172–1179. [Google Scholar] [CrossRef]
  9. Khalloufi, S.; Almeida-Rivera, C.; Janssen, J.; Van-Der-Vaart, M.; Bongers, P. Mathematical Model for Simulating the Springback Effect of Gel Matrixes During Drying Processes and Its Experimental Validation. Dry. Technol. 2011, 29, 1972–1980. [Google Scholar] [CrossRef]
  10. Changyou, L.; Yaojian, S.; Junichi, K. An Analytical Solution of the Granular Product in Deep-Bed Falling Rate Drying Process. Dry. Technol. 1998, 17, 1959–1969. [Google Scholar] [CrossRef]
  11. Ostanek, J.; Ileleji, K. Conjugate heat and mass transfer model for predicting thin-layer drying uniformity in a compact, crossflow dehydrator. Dry. Technol. 2020, 38, 775–792. [Google Scholar] [CrossRef]
  12. Xu, Z.; Pillai, K.M. A pore-network study on the factors influencing the isothermal drying of single- and dual-scale porous media. Dry. Technol. 2021, 39, 1294–1313. [Google Scholar] [CrossRef]
  13. Ozturk, O.K.; Takhar, P.S. Hybrid mixture theory-based modeling of moisture transport coupled with quality changes in strawberries and carrots. Dry. Technol. 2021, 39, 932–949. [Google Scholar] [CrossRef]
  14. Lígia Cavalcante de Almeida Farias Aires, K.; Pereira da Silva, W.; Everton de Farias Aires, J.; Freire da Silva Júnior, A.; Maria Diniz Pereira da Silva e Silva, C. Convective drying of osmotically dehydrated apples described by three-dimensional numerical solution of the diffusion equation with analysis of water effective diffusivity spatial distribution. Dry. Technol. 2019, 37, 2034–2046. [Google Scholar] [CrossRef]
  15. Aires, K.L.C.d.A.F.; da Silva, W.P.; de Farias Aires, J.E.; da Silva Júnior, A.F.; Silva, C.M.D.P.d.S.e. Apple osmotic dehydration described by three-dimensional numerical solution of the diffusion equation. Dry. Technol. 2018, 36, 1970–1981. [Google Scholar] [CrossRef]
  16. Nowicki, S.C.; Davis, H.T.; Scriven, L.E. Microscopic Determination of Transport Paraheters in Drying Porous Media. Dry. Technol. 1992, 10, 925–946. [Google Scholar] [CrossRef]
  17. Prakash, B.; Bingol, G.; Pan, Z. Moisture Diffusivity in Rice Components during Absorption and Desorption. Dry. Technol. 2011, 29, 939–945. [Google Scholar] [CrossRef]
  18. Beke, J.; Mujumdar, A.S.; Bosisio, R.G. Drying of Fresh and Rewetted Shelled Corn in Microwave Fields. Dry. Technol. 1995, 13, 463–475. [Google Scholar] [CrossRef]
  19. Liu, X.; Yang, B. Review and vista on drying theories of porous medium. J. China Agric. Univ. 2005, 10, 81–92. [Google Scholar] [CrossRef]
  20. Chen, G.; Wang, W.; Mujumdar, A.S. Theoretical study of microwave heating patterns on batch fluidized bed drying of porous material. Chem. Eng. Sci. 2001, 56, 6823–6835. [Google Scholar] [CrossRef]
  21. Budd, C.J.; Hill, A.D.C. A comparison of models and methods for simulating the microwave heating of moist foodstuffs. Int. J. Heat Mass Transf. 2011, 54, 807–817. [Google Scholar] [CrossRef]
  22. Mujumdar, A.S. Editorial: The Role of Drying Technology in Sustainable R&D and Innovation. Dry. Technol. 2012, 30, 1387. [Google Scholar] [CrossRef]
  23. Cohen, J.S.; Yang, T.C.S. Progress in food dehydration. Trends Food Sci. Technol. 1995, 6, 20–25. [Google Scholar] [CrossRef]
  24. Putranto, A.; Chen, X.D. An assessment on modeling drying processes: Equilibrium multiphase model and the spatial reaction engineering approach (S-REA). Chem. Eng. Res. Des. 2015, 94, 660–672. [Google Scholar] [CrossRef]
  25. Li, C.; Liu, J.; Chen, L. The Moisture Distribution of High Moisture Content Rough Rice during Harvesting, Storage, and Drying. Dry. Technol. 2003, 21, 1115–1125. [Google Scholar] [CrossRef]
  26. Li, C.; Ban, H.; Shen, W. Self-Adaptive Control System of Grain Drying Device. Dry. Technol. 2008, 26, 1351–1354. [Google Scholar] [CrossRef]
  27. Li, C.; Mai, Z.; Fang, Z. Analytical study of grain moisture binding energy and hot air drying dynamics. Trans. Chin. Soc. Agric. Eng. 2014, 30, 236–242. [Google Scholar]
  28. Li, C. Analytical theory study on latent heat coefficient of grain water vaporization. Dry. Technol. 2021, 39, 803–819. [Google Scholar] [CrossRef]
  29. Makeev, A.G.; Semendyaeva, N.L.; Slinko, M.M. Synergetic effect and oscillatory behavior of CO oxidation over a bimetallic composite catalyst. Chem. Eng. J. 2015, 282, 3–10. [Google Scholar] [CrossRef]
  30. Li, J.; Liu, Y.; Wang, X.; Zhao, G.; Ai, Y.; Han, B.; Wen, T.; Hayat, T.; Alsaedi, A.; Wang, X. Experimental and theoretical study on selenate uptake to zirconium metal–organic frameworks: Effect of defects and ligands. Chem. Eng. J. 2017, 330, 1012–1021. [Google Scholar] [CrossRef]
  31. Pan, B.; Chen, D.; Zhang, H.; Wu, J.; He, F.; Wang, J.; Chen, J. Stability of hydrous ferric oxide nanoparticles encapsulated inside porous matrices: Effect of solution and matrix phase. Chem. Eng. J. 2018, 347, 870–876. [Google Scholar] [CrossRef]
  32. Nguyen Tran, H. Comment on “simultaneous and efficient removal of Cr(VI) and methyl orange on LDHs decorated porous carbons”. Chem. Eng. J. 2019, 359, 810–812. [Google Scholar] [CrossRef]
  33. Li, C. Characteristic functions of drying material system and its theoretical solution. Trans. Chin. Soc. Agric. Eng. 2020, 36, 286–295. [Google Scholar] [CrossRef]
  34. Yang, W.; Yang, Y.; Yang, H.; Zhou, H. Regulating Water Activity for Rechargeable Zinc-Ion Batteries: Progress and Perspective. ACS Energy Lett. 2022, 7, 2515–2530. [Google Scholar] [CrossRef]
  35. Hosokawa, A.; Motohashi, K. Drying Characteristics of a Single Grain of Rough Rice. J. Jpn. Soc. Agric. Mach. 1971, 33, 53–59. [Google Scholar] [CrossRef]
  36. Li, C. Engineering Thermodynamics and Heat Transfer; Agricultural University Press: Beijing, China, 2012. (In Chinese) [Google Scholar]
  37. Yu, Y.; Wang, J. Modeling Equilibrium Moisture Content of γ-Ray Irradiated Rough Rice. Dry. Technol. 2006, 24, 671–676. [Google Scholar] [CrossRef]
  38. Reddy, B.S.; Chakraverty, A. Equilibrium Moisture Characteristics of Raw and Parboiled Paddy, Brown Rice, and Bran. Dry. Technol. 2004, 22, 837–851. [Google Scholar] [CrossRef]
  39. Chen, S.-F.; Danao, M.-G.C. Modeling the equilibrium moisture content (EMC) of Miscanthus sinensis, miscane, energy cane, and energy sorghum. Fuel 2015, 147, 18–26. [Google Scholar] [CrossRef]
  40. Li, C. Research progress in analytical theory of drying system. Trans. Chin. Soc. Agric. Eng. 2024, 40, 1–13. [Google Scholar] [CrossRef]
Figure 1. (a) The moisture content changed with RH. (b) The forms of moisture in Paddy Rice.
Figure 1. (a) The moisture content changed with RH. (b) The forms of moisture in Paddy Rice.
Agronomy 14 01728 g001
Figure 2. Two−stage characteristics of paddy drying Note: I and II denote the 1st and 2nd descending drying sections, respectively. Me represents the equilibrium moisture content (% d.b.), Me1 represents the equilibrium moisture content (% d.b.) of the 1st descending drying section, Mk is the cut−off moisture content (%), and Mc is the initial moisture content (% d.b.).
Figure 2. Two−stage characteristics of paddy drying Note: I and II denote the 1st and 2nd descending drying sections, respectively. Me represents the equilibrium moisture content (% d.b.), Me1 represents the equilibrium moisture content (% d.b.) of the 1st descending drying section, Mk is the cut−off moisture content (%), and Mc is the initial moisture content (% d.b.).
Agronomy 14 01728 g002
Figure 3. The changes in drying rate ratio with different water activities.
Figure 3. The changes in drying rate ratio with different water activities.
Agronomy 14 01728 g003
Figure 4. (a) The variation pattern of rice water activity. (b) The theoretical drying rate ratio of rice. (c) Change in moisture content in grain.
Figure 4. (a) The variation pattern of rice water activity. (b) The theoretical drying rate ratio of rice. (c) Change in moisture content in grain.
Agronomy 14 01728 g004
Figure 5. (a) The consumption of free energy during the drying process of rice. (b) The variation pattern of the ratio of free energy during the drying process of rice.
Figure 5. (a) The consumption of free energy during the drying process of rice. (b) The variation pattern of the ratio of free energy during the drying process of rice.
Agronomy 14 01728 g005
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Wang, H.; Li, S.; Zhang, Y.; Ali, K.A.M.; Wu, W.; Li, C. The Analytical Theory of Water Activity in the Paddy Rice Drying Process. Agronomy 2024, 14, 1728. https://doi.org/10.3390/agronomy14081728

AMA Style

Wang H, Li S, Zhang Y, Ali KAM, Wu W, Li C. The Analytical Theory of Water Activity in the Paddy Rice Drying Process. Agronomy. 2024; 14(8):1728. https://doi.org/10.3390/agronomy14081728

Chicago/Turabian Style

Wang, Han, Shengtao Li, Ye Zhang, Khaled Abdeen Mousa Ali, Weihui Wu, and Changyou Li. 2024. "The Analytical Theory of Water Activity in the Paddy Rice Drying Process" Agronomy 14, no. 8: 1728. https://doi.org/10.3390/agronomy14081728

APA Style

Wang, H., Li, S., Zhang, Y., Ali, K. A. M., Wu, W., & Li, C. (2024). The Analytical Theory of Water Activity in the Paddy Rice Drying Process. Agronomy, 14(8), 1728. https://doi.org/10.3390/agronomy14081728

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop