A refined mechanistic model for swelling kinetics of starch granules

Abstract

This paper investigates the gelatinization of individual starch granules using numerical simulations, validated against experimental microscopy data from a ParCS apparatus. We show that the dynamics of starch-granule swelling can be captured by a diffusion equation for mass transfer of water into the granule, with the equilibrium water content captured by a Flory–Rehner theory of a cross-linked network in which the fraction of cross-linked chains is made to vary as an empirical function of temperature. Having the cross-link density vary with temperature is vital to capture the swelling behavior at large and small swelling extents (i.e., close to and far away from the gelatinization temperature). The theory produces excellent agreement with both equilibrium swelling data and dynamic swelling data for red bean starch. Further, we show that the model is able to reproduce a previous experimental finding that swelling data from different granules from red bean, chickpea, green lentil, and yellow pea starches can be collapsed onto a universal curve with only two empirical parameters. The simulations are then used to predict the relationship between the empirical parameters in the master curve and the true material properties. The modified theory presented here is a major step forward in the fundamental understanding of starch gelatinization and the ability to use predictive models for optimization of industrial manufacturing processes.

---

1 Introduction

Starch is the main source of carbohydrates for humans in their diet, and makes up a significant portion of their caloric intake.¹ Starch is obtained from agricultural products such as rice, corn, potato, and legumes,² and determines the texture, mouth-feel, shelf life and consistency of many foods.³ Starch naturally exists in the form of discrete, insoluble, semi-crystalline granules that differ in size and shape depending on the plant origin. These granules may range from around 1 μm to more than 100 μm,⁴ and the shape can be ovoid, disc-like, polyhedral or irregular.⁵ The granules themselves are mostly composed of unbranched amylose and highly branched amylopectin.¹ The relative amounts of these two polymers and their molecular structure lead to fundamental variability in the physical properties of different starch varieties,¹ and result in differences in how the different starches behave during processing, as well as the attributes of the final product.⁶

One of the major changes that occurs during starch processing is gelatinization, which is the sudden intake of water into their granules above a critical temperature that occurs between 60 °C and 80 °C for most starch varieties.^7,8 It is known that starch gelatinization is affected by many factors such as the heating rate, heating time, moisture, and presence of other ingredients.^9–11 Experimentally, gelatinization is characterized by many methods. Rheometry is widely used to obtain bulk measurements of gelatinization temperature, while differential scanning calorimetry is used to obtain this quantity as well as the gelatinization enthalpy.^12–16 To examine gelatinization at the individual granule level, optical microscopy can be used to monitor morphological changes during heating.^17–19 Recent studies have suggested that there is considerable variability in gelatinization between granules from the same starch sample, causing the study of starch gelatinization to be a complex undertaking.^19–21

Besides experimental analysis, theoretical analysis and quantitative prediction of starch gelatinization are in great demand and significant effort has been made in this area. For example, recently Desam et al. proposed a first-principles model for the swelling kinetics and gelatinization behavior of waxy maize starch.²² This model consists of three elements: (a) a Flory–Rehner theory^23–27 for the equilibrium swelling of starch granules (treated as a polymer network), (b) a heat and mass transport model for water diffusion into granules to describe swelling kinetics towards equilibrium, and (c) a population balance model to describe how the granule size distribution evolves over time. Subsequently, these ideas have been used to describe the swelling of many other starches (potato, rice),^17,22,28 as well as forecast the rheology of starch dispersions during the initial stages of swelling.²⁸ In other work, Evans et al. used the Flory theory to explain the effect of solutions on the gelatinization temperature of potato starch.²⁹

While these models are promising, more work is needed to validate them and improve them. For example, Flory theories of polymer swelling often assume reversible swelling, when in fact starch swelling is known to be irreversible. Additionally, Flory theories assume the granule is a loosely cross-linked polymer network, which is reasonable during the late stages of swelling but may not be the case during the early stages when the granule is tightly packed. Prior to the present work, it has not been clear if these considerations need to be taken into account to improve modeling.

With this motivation in mind, this study aims to make a systematic comparison between the theoretical Flory models of Desam et al.²² and recent experimental measurements from our group.¹⁹ The experimental study described in ref. 19 performed measurements of swelling kinetics for hundreds of individual granules during gelantinization from four different legume plants (red bean, chickpea, green lentil, and yellow pea). The data showed that all of the data from the hundreds of granules can be described by a single empirical function with only four free parameters, where two parameters correspond to the initial and final granule size. Because of this, the data obtained in that study present a highly convenient case for validation of theoretical models that can predict granule swelling. New experimental data were also gathered for the swelling of red bean, to be compared against existing models.

In this work, we will test the experimental data against the model from Desam et al.²² to see if the model can be validated. If necessary, the assumptions made in the model will be revisited and the model revised. We hypothesize that the theoretical model of Desam et al. can be non-dimensionalized in such a way that (with appropriate constraints) it will also result in a universal curve as found experimentally by Mo et al. in ref. 19. We further hypothesize that it will be possible to make a clear mathematical connection between the theoretical model parameters of Desam et al. (based on material properties) and the empirical model parameters of Mo et al.

2 Methods and materials

2.1 Materials

Starch from red bean (Vigna angularis) (President's Choice Canada) was isolated following a the method outlined in ref. 30 with some minor modifications. Briefly, 100 g of dried beans was soaked in 200 mL of 0.5 wt% sodium bisulfite solution for 20 h at 4 °C, then drained and blended with 500 mL of distilled water for 5 min. The mixture was filtered through a 125 μm sieve, then centrifuged at 1500g for 15 minutes. The supernatant was discarded and the starch was resuspended in distilled water for centrifuging again under the same conditions. The top non-white layer in the centrifuge tube was removed and discarded. Finally, the starch was dried at 40 °C for 24 h and gently broken into a powder using a mortar and pestle.

2.2 Equilibrium granule swelling measurements

Starch samples with 20 ± 5 mg were dispersed in 17 mL of deionized water in a 20 mL vial. After shaking by hand, a small amount of solution was injected into the ParCS chamber (custom instrument described in ref. 31) using a syringe. The temperature was monitored and controlled by a platinum series universal benchtop PID controller (OMEGA Engineering, Norwalk, CT, USA) to achieve a sequence of increasing temperatures.

The temperature was controlled according to the following steps. First, the chamber was equilibrated to 55 °C (328 K) and maintained for 4 h, after which we assumed that equilibrium was reached and the size of the granule was recorded. Then, the temperature was increased to 333 K, 338 K, 340.5 K, 343 K, 345.5 K, 353 K, 355.5 K, and 358 K in sequence. Each increment in the temperature was completed within 30 s and each temperature was maintained for 7 h before recording the size and then increasing the temperature.

The swelling process of the individual starch granules was captured by camera and image analysis was done using a custom Python script. Using the Python script, the boundaries of the granules were manually outlined and the projected area (A) was computed. The equivalent spherical diameter (d) of the starch granules was then calculated using the following formula:

It was noted during method development that there was a fair amount of variation in the equilibrium granule size for starch granules with different swelling ratios. Therefore, when preparing the data for comparison to theory, granules were only included in the data set if the swelling ratio (ratio of swollen and unswollen granule diameters) was 2.15 ± 0.24 (approximately 10% variation about the mean). Due to the painstaking nature of making this measurement, data from only 12 granules were included in this manuscript due to the need to restrict the range of swelling ratios, and this was found to be sufficient for the present purpose.

3 Theory and numerical methods

3.1 Equilibrium granule swelling (Flory–Rehner theory)

The Flory–Rehner theory states that the free energy of a polymer network is determined by the competition between water–polymer mixing and the elastic stretching of the polymer chains. In this work, we will assume that a starch granule can be modeled as a polymer network that can be swollen in water, throughout the gelatinization process. We will also assume that the granule is neutral and the internal granule architecture is uniform. A number of more refined assumptions related to the actual structure of starch granules could be made based on previous theoretical work.^32–45 However, incorporating these additional assumptions is not trivial in light of our specific objectives and hypotheses, and so will be reserved for future work.

At equilibrium, the chemical potential of water within the granule (μ₁) must be equal to its chemical potential outside the granule (μ_1,b). In this work we assume that the expression for the chemical potential difference is given by the following equation:

In eqn (2), ϕ is the volume fraction of starch in the granule, , T, χ(T), ν* denote the ideal gas constant, temperature, the Flory–Huggins parameter, and the fraction of cross-linked chains, respectively. The first two terms on the right hand side of eqn (2) are derived from mixing entropy, the third term from mixing enthalpy, and the fourth one from elastic entropy due to polymer stretching. In the elastic portion of the model, we have assumed that the reference state (unstretched state) corresponds to a fully dry granule (ϕ = 1), though other assumptions are possible and will be explored in future work.⁴⁰

Since starch swelling occurs mainly above the gelatinization temperature T_g, the molar enthalpy of interaction with water, χ, will be assumed to be constant below T_g, and only dependent on temperature. Others have shown that under certain circumstances, χ is effectively independent of temperature.³² However, in this work we have elected to use the temperature dependence from the lattice model used in Flory–Huggins solution theory:^46,47

where ΔH_mix is the enthalpy of mixing between starch segment and water. In this model, we will assume that the enthalpy of mixing ΔH_mix is equal to enthalpy of gelatinization ΔH, although this may overpredict the variation of χ with temperature.

Lastly, we also assume that χ is independent of composition, though we will revisit that assumption in future work. Thus, we can treat the Flory–Huggins parameter χ(T) as a piece-wise function as done in ref. 22:

where ΔH is the enthalpy of gelatinization (typically determined by DSC). An individual granule will reach an equilibrium state when the chemical potential inside and outside the granule are equal –, that is, the left-hand side of eqn (2) is zero. The equilibrium starch volume fraction is thus the solution of the following nonlinear equation:

Further, the conservation of starch inside the granule requires that:

where ϕ_eq is the equilibrium starch volume fraction, ϕ₀ is the initial starch volume fraction, and R₀ is the initial granule size. Solving the nonlinear eqn (4) for starch volume fraction and plugging into eqn (5) allows one to determine the equilibrium granule size. The inputs to this theory are the temperature T, initial starch composition ϕ₀, initial radius R₀, cross-linking density ν*, starch–water interaction parameter χ₀, and gelatinization temperature and enthalphy of gelatinization (T_g, ΔH).

3.2 Kinetics of granule swelling (diffusion equation)

When a starch granule is heated, gelatinization occurs. During this process, the difference in water's chemical potential between the outside and inside of the granule leads to mass transfer of water into the granule. This process occurs until the granule reaches chemical equilibrium. The temporal evolution of chemical potential of water inside a spherical granule μ₁ can be described by a diffusion equation:²²Here t and r denote the time and radial position, respectively; D is the empirical diffusion coefficient of water into the granule; and R(t) is the granule radius which increases in time.

The chemical potential μ₁ is related to the local starch volume fraction ϕ(r,t) through the Flory–Rehner equation (eqn (2)). The starch volume fraction is related to the granule size R(t) through conservation of mass:

where [small phi, Greek, macron] is the average starch volume fraction in the granule. Lastly, the effective diffusivity D is assumed to have a dependence on temperature that follows the Arrenhius equation, as well as a dependence on the porosity and tortuosity^45,48 as follows:Here, E = 1890 K⁻¹ and T_ref = 323 K are chosen to provide the best fit to diffusivity data for pure water (i.e., when [small phi, Greek, macron] = 0),⁴⁴ and can be regarded as constants in the model. The effective diffusivity is taken to be proportional to the porosity (1 − [small phi, Greek, macron]) and inversely proportional to tortuosity as is commonly done in a Darcy model of a porous medium, and the tortuosity is related to porosity according to an empirical model proposed in ref. 45, and here we set c = 0.77 as done in that work.

We assume that the temperature inside the starch granule is uniform and equal to the outside water temperature. This assumption is reasonable when the Lewis number (ratio of thermal diffusivity to mass diffusivity) is very large, as is expected here. In other words, the timescale of heat transfer is much smaller than that of molecular diffusion of water. Thus, the temperature equation describing the heat conduction in the granule is neglected in this study, and thermal equilibrium is assumed throughout the process. In turn, this means that the kinetics of swelling (as opposed to the equilibrium effects) can only be governed by mass transfer.

The following dimensionless variables are adopted to re-write eqn (6)–(8):

where μ_1,b denotes the chemical potential of pure solvent at the edge of granule, while μ₁₀ is at the initial state. Eqn (6) can then be recast as:which is subject to dimensionless initial and boundary conditions:In eqn (13), R*(t*) = R(t)/R₀ is the dimensionless granule radius at dimensionless time t* and is the result of assuming the granule surface is in equilibrium with water when the external resistance to mass transfer is negligible.

The chemical potential is related to the volume fraction through the Flory equation. In dimensionless form, we have:

while for conservation of mass we have:Eqn (8)–(15) must be solved simultaneously to obtain the granule size as a function of time as it approaches equilibrium.

3.3 Simulation details

All the simulations are performed using Python on a laptop with Intel(R) Core(TM) i7-9750H CPU@2.60 GHz and 16.0 GB RAM. In the calculation, dimensionless eqn (10) subject to the conditions in eqn (11)–(13) are solved firstly by adopting an implicit finite difference method:where t_n = nτ, n = 0, 1, …, r_j = jh, j = 0, 1, …, and τ and h are time step and space step sizes, respectively. At the beginning of the simulation, D*(T) was set to be the diffusion coefficient at the initial temperature T₀, while R*(t*) in the boundary condition, eqn (13), is the dimensionless granule radius at initial time t₀. The diffusion coefficient and granule radius were evaluated and updated at each t_n as discussed in details in the following paragraph.

Given a radius R*(t_n) and diffusion D*(t_n) at time t_n, we compute the chemical potential at time t_n+1 using the implicit finite difference method in eqn (16). With the chemical potential obtained from eqn (10), Flory–Rehner theory eqn (2) is utilized to obtain the starch volume fraction ϕ at time t_n+1 by using the “root” function from the scipy module in Python. The average starch volume fraction [small phi, Greek, macron](t*) at time t_n+1 is calculated as:

By substituting [small phi, Greek, macron](t*) into eqn (7), the starch granule size R*(t_n+1) at time t_n+1 is obtained. We compute D*(t_n+1) using eqn (8), and replace the diffusion coefficient D*(t_n) in eqn (10) and granule size R*(t_n) in eqn (13) with the new D*(t_n+1) and R*(t_n+1). Having computed all quantities at t_n+1 and updated all equations, we repeat the procedure to compute quantities at the next time step t_n+2, and so on. Thus, an iteration is constructed until the size of starch granule at a set time is acquired. The parameter values that are not changed in any of the simulations are given in Table 1 and reflect the values in ref. 22. The default values for all other parameters, when not being varied for analysis, are given in Table 2.

Table 1 Model parameters used in the simulations that are constant for all simulations

Parameter	Value
ϕ 0	0.586
R 0	3 × 10^−5 m
D ref	3.98 × 10^−9 m^2 s^−1
T ref	323 K

---

Table 2 Default values of the model parameters used in the simulations that are used when not being explicitly varied

Parameter	Value
ν ∞	0.012
a	0.5 K^−1
T 0	323 K
T g	339 K
ΔH	1.35 × 10^4 J mol^−1
χ 0	0.497
c	0.9

---

4 Results and discussions

4.1 Missing physics in standard Flory–Rehner theory

In the process of validating the theory outlined in Section 3, we observed a non-physical result that stems from the Flory–Rehner portion of the model. In this section, we present the results of some of the simulations of granule swelling (radius vs. time) that illustrate the deficiency of the model. In each of the simulations in this section we examined a starch granule with the physical parameters listed in Table 1 and a gelatinization temperature of T_g = 339 K. The granule was exposed to a temperature profile that increased linearly in time as T = T₀ + bt, with four different heating rates b, from an initial temperature T₀ = 323 K to a final temperature T_f = 363 K where the temperature was then held steady.

In Fig. 1(a), we see that the granule size undergoes two distinct stages. The first stage occurs when T < T_g and the granule swells rapidly to an intermediate size and then stops increasing. The second stage occurs when T > T_g and the granule swells again to its final size corresponding to the equilibrium size at T = T_f. In experiments, only the second stage is observed, and to our knowledge previous modeling efforts have also only shown the second stage (e.g., ref. 22). However, stage 1 shown here using the same physics as ref. 22, is inconsistent with experiments, as granule swelling should not occur appreciably below the gelatinization temperature.

Fig. 1 Demonstration of a problem with the standard Flory–Rehner theory when applied to starch granules. (a) Dimensionless granule radius R* versus dimensionless time t* showing the impact of changing the heating rate, where b is the heating rate (K s⁻¹) according to T = T₀ + bt. (b) Average starch volume fraction [small phi, Greek, macron] versus dimensionless time showing the impact of varying the heating rate. Both plots show that there is an initial swelling that occurs for T < T_g, when in actuality no appreciable swelling should occur.

To illustrate what is happening in stage 1 more clearly, we plot the average starch volume fraction [small phi, Greek, macron](t*) vs. time in Fig. 1(b) for two of the heating rates. Fig. 1(b) reveals that the volume fraction will drop rapidly from the given initial value ϕ₀ to a steady state value before the gelatinization temperature is reached. This suggests that the specified initial starch volume fraction of the granule does not correspond to the equilibrium value calculated from Flory–Rehner theory for temperatures below the gelatinization temperature.

To verify that the initial volume fraction does not correspond to an equilibrium condition, we plot the equilibrium starch volume fraction ϕ_eq as a function of temperature in Fig. 2(a). The values of ϕ_eq in this plot are calculated from the Flory–Rehner theory using eqn (4). The data show that below the gelatinization temperature T < T_g the equilibrium starch volume fraction is constant at ϕ_eq ≈ 0.17, which is much smaller than typical initial conditions of ϕ₀ > 0.5 (and ϕ₀ = 0.586 in this case) and hence drives strong swelling. In actuality, the equilibrium volume fraction for T < T_g should be very close to ϕ₀ so that no appreciable swelling will be observed.

Fig. 2 (a) Illustration of the narrow range of equilibrium starch volume fractions that are possible with the standard Flory–Rehner theory. (b) Illustration of the impact of cross-linking density parameter ν* and the solvent interaction parameter χ₀ on the equilibrium starch volume fraction at initial temperature T = T₀.

Before jumping to the conclusion that this points to a flaw in the theory, it is important to determine if this discrepancy is simply the result of poorly chosen material parameters (e.g., Table 1). To check this, we examined the effect of the two material parameters that influence ϕ_eq at T₀: the density of cross-links ν* and the starch–water interaction parameter χ₀. Fig. 2(b) shows that ϕ_eq at T₀ is rather insensitive to χ₀, but quite sensitive to ν*. However, we found that if ν* is adjusted so that ϕ_eq = ϕ₀ at T₀, swelling is severely restricted and it is not possible to reach realistic values of the final volume fraction. Thus, we can conclude that there is a deeper problem with the use of Flory–Rehner theory for this application.

We propose that the problem with the standard Flory–Rehner model is that it is likely a missing piece of the underlying physics. Generally, Flory swelling theories assume constant material properties (e.g., ν* and χ₀ are constant), while in reality the cross-link density in a starch granule, ν*, may decrease during swelling as some of the physical cross-links break, which break much more easily than chemical cross-links. In turn, the starch–water interaction parameter, χ₀, could also change since more starch–water sites would be exposed as the cross-links break. We assume that from a microstructural viewpoint, the granule can be considered to appear initially as a tightly-packed polymer glass that undergoes a phase transition to become a loosely-packed polymer network as the degree of gelatinization increases.^38,39

With this physical picture of a granule that undergoes a phase transition, it makes sense that a single cross-link density would be incapable of accurately describing the polymer network both before and after gelatinization as shown in Fig. 2. Specifically, if we fix χ₀ = 0.497, the cross-link density that results in the correct initial volume fraction for T < T_g is ν* = 0.23. Similarly, the cross-link density that results in the correct final volume fraction after gelatinization at 363 K is ν* = 0.0042. Therefore, we hypothesize that a simple modification to the Flory–Rehner theory that allows for a variable value of the cross-link density will allow us to qualitatively capture granule swelling from the unswollen size to the final granule size.

Here, we will test this hypothesis by examining four of the simplest possible dependencies of the cross-link density on temperature. These four functional relationships are depicted in Fig. 3(a) and given by:

Fig. 3 (a) Five different models that allow ν* to vary as a function of temperature. For the functions plotted here, the following parameters were used for convenience to make visualization easier: ν₀ = 0.23, ν_∞ = 0.0042, T_g = 339 K, a_linear = −0.04, a_exp = 0.23, and a = 0.82. (b) Plot of ϕ_eqvs. temperature calculated from each of the models in (a). For both plots, the following common parameters are used: χ₀ = 0.497, ΔH = 1.281 × 10⁴ J mol⁻¹.

(1) a step function

(2) a linear function

(3) a negative exponential function

and (4) a modified logistic function

In these four empirical models T_g is the gelatinization temperature, ν₀ is the cross-link density at T < T_g that sets the initial starch volume fraction ϕ₀, ν_∞ is the cross-linking density that sets the equilibrium swelling at the final temperature T = T_f, and a_linear, a_exp, and a are free parameters to be determined.

The advantage of these four models is that they have at most one free parameter that cannot be determined from material properties alone. The functions were chosen in such a way that they have increasing smoothness of the transition between different cross-link densities. Also, all of the models are designed to produce a relatively sharp transition in the cross-link density after T > T_g since prior work showed a reasonable agreement for the latter stages of granule swelling with a constant value of the cross-link density equal to ν_∞.²²Fig. 3(b) shows how each of these four models influences the equilibrium starch volume fraction ϕ_eq as a function of temperature.

4.2 Comparison of modified theory to experimental results

Next, we examined the ability of each of the four empirical models for cross-link density (eqn (18)–(21)) to improve the accuracy of the numerical simulations of granule swelling. In Fig. 4, we study how the four different models compare against granule swelling data from literature. In ref. 19, gelatinization of individual starch granules from four types of legume starches (chickpea, red bean, green lentil, and yellow pea) was observed, and it was shown that all of the data from hundreds of granules of each type could be described by a universal curve (after appropriate normalization). We will first attempt to validate our model against these dynamic swelling measurements.

Fig. 4 Comparison of the best fit to the experimental data using the four different empirical models for the two extreme values of the swelling ratio: (a) R_f/R₀ = 3.21 and (b) R_f/R₀ = 1.69. (c) and (d) show the error between the different models and the experimental data from (a) and (b), respectively.

In the experiments, the heating was done via a quadratic function of time as follows:

T = T₀ + T₁t − T₂t² (22)

where the initial and final temperatures are T₀ = 323 K and T_final = 363 K, respectively, T₁ = 0.1759 K s⁻¹, and T₂ = 1.92 × 10⁻⁴ K s⁻². Because the final granule size varies from granule to granule, we compare the simulations for the two extremes observed in the experiments of ref. 19, with granule swelling ratios (i.e., ratio of final granule size to initial granule size): R_f/R₀ = 3.21 and R_f/R₀ = 1.69. The free parameters used in each of the four empirical models, as well as all the other free parameters listed in Table 2 were optimized^49,50 to produce the best fit against the data for each model (i.e., the values were allowed to be different for each model). In the optimization scheme, the following cost function was minimized:where differences at n points between the numerical result R^num_i and experiments R^exp_i are accumulated. The values of S_global observed for the four models are shown in Table 3, and it can be clearly seen from both Fig. 4 and Table 3 that numerical simulations using the modified logistic function for varying the cross-link density give the best fit to the experimental data.

Table 3 Global differences between predictions from four empirical models and experimental data along with the best fit parameters for each model that resulted in these values. The units of each parameter are as follows: a_linear [=] K⁻¹; a_exp [=] K⁻¹; a [=] K⁻¹

Model (free parameters)	S global		Best fit parameter
Step	1.5394	1.9032	—	—
Linear (alinear)	1.2509	0.7216	−0.0625	−0.0230
Negative exponential (aexp)	1.2677	0.7540	0.2573	0.0857
Modified logistic (a)	0.0381	0.0445	0.6068	0.4884

---

Based on the comparison between simulations and data in Fig. 4, we can conclude that, indeed, allowing for the cross-link density to vary with temperature results in a significant improvement in the ability of the simulations to capture the dynamics of granule swelling. However, Flory–Rehner is inherently an equilibrium theory and therefore it is important to also compare this modification against equilibrium granule swelling. To do this we measured the equilibrium granule size as a function of temperature for 12 individual red-bean, starch granules that all had a similar swelling ratio (see Section 2.2). Because we do not have a way of directly measuring the initial starch volume fraction, we assumed that all granules start with an equilibrium starch volume fraction of 0.586. We then calculated the equilibrium starch volume fraction as a function of temperature and the results of these measurements are shown in Fig. 5.

Fig. 5 (a) Comparison between experimental values of the equilibrium volume fraction of red bean starch granules (labeled as “data” in the legend) and the prediction of three different models. The experimental data correspond to 12 different granules. (b) and (c) Comparison of ν*(T) and 1/ν*(T) from the modified logistic function and eqn (24) fit to the data. Note, the free parameters used in the modified logistic function shown here correspond to the best fit for a swelling ratio of 2.15 to match the data set.

As with the dynamic data from simulations, it is clear from Fig. 5(a) that the Flory–Rehner theory with a constant cross-link density is unable to adequately describe the data. However, it is also evident that a modified logistic function for the cross-link density is not a perfect fit either. Therefore, we also proposed the following fitting function for describing the data with a relatively simple analytical expression:

where the optimal parameter values are: A = −0.33, d = 0.78 K⁻¹, m = 0.0066 K⁻¹, b = 2.75, ΔT = 341.24 K.

By substituting this fitting function (eqn (24)) into eqn (4), we can also back-calculate ν*(T) to compare it to the modified logistic function for the cross-link density as shown in Fig. 5(b) and (c). As with a comparison of ϕ_eq(T), the data do not agree perfectly with the modified logistic function and the inherent differences in the functional form cause large differences in cross-link density both at low and high temperature. However, the data do support the idea that the current Flory–Rehner theory is inadequate and that using a variable cross-link density that varies sigmoidally with temperature is a significant improvement. Interestingly, the degree of gelatinization for rice starch was also found to vary sigmoidally with temperature through DSC experiments,³⁷ though in that work it is unclear if the sigmoidal character is related to the phase transition or due to intrasample variability (see ref. 19 for a discussion of intrasample variability).

Although, it would not be feasible as a general strategy, we also attempted to use eqn (24) directly for the simulations to see if this resulted in better predictions of the dynamics of swelling. To do this, we used the values of ν*(T) calculated from eqn (24) in the swelling simulations. Fig. 6 shows the results of this simulation and from this we learn the following. First, the swelling dynamics are extremely sensitive to minor variations in the value of ν* at later stages of swelling, as evidenced by the differences in the shape of the swelling curve for temperatures above 340 K where the differences in ν* are small between the data and the logistic function. We also learn that, in contrast to the effect at higher temperatures, very large differences in the value of ν* make only a small difference at lower temperatures when swelling is just beginning to occur.

Fig. 6 Comparison of the swelling process simulated using the modified logistic function for ν* with the swelling process simulated using eqn (24).

In summary, we have considered multiple options for correcting the missing physics in the model. It is clear that a variable cross-link density is necessary and that the modified logistic function captures this in a way that produces realistic swelling dynamics. Therefore, we conclude that it is reasonable to proceed to perform simulations where ν* varies sigmoidally with temperature and in the remainder of this paper will use the modified logistic function in all simulations.

4.3 Relationship between modified Flory model and experimental master curve

As discussed in the previous section, we determined that by using a modified logistic function for the cross-link density ν*(T) in the Flory–Rehner swelling model, one can reasonably predict dynamic experimental swelling data. In a previous paper, we also found that for a given heating profile (eqn (22)), the experimental swelling data for four types of starch granules follow a universal curve with only two empirical parameters:where the dimensionless parameter k_G is related to the characteristic time scale for gelatinization; while the dimensionless parameter t_G relates to the time – and by extension the temperature – at which swelling is rapid. Therefore, in this section we will attempt to determine if the empirical parameters (k_G and t_G) can be related to physical parameters in the theoretical model. To do this, we will vary the model parameters that remain after non-dimensionalizing the governing equations and attempt to collapse the resulting swelling curves. Unless specified otherwise, all model parameters used in the simulations for this section will be the same as those listed in Tables 1 and 2.

In the sections that follow, we will systematically vary the parameters, one at a time, and generate dynamic swelling curves. We will then attempt to collapse the resultant numerical data onto a master curve by shifting the time and/or rescaling the time according to eqn (25) in order to collapse the numerical data onto a master curve. Finally, we will examine the relationship between the empirical parameters and the theoretical parameters to determine what relationship may exist between them. Ideally, such relationships will be obtained as some specific, but as yet unknown, mathematical functions as follows:

t_G = f(a,ν_∞,ϕ₀,T_g,T₀,ΔH,χ₀,c) (26)

k_G = g(a,ν_∞,ϕ₀,T_g,T₀,ΔH,χ₀,c) (27)

4.3.1 Effect of empirical parameters used to vary cross-link density: a, ν_∞, and ϕ₀. In this section we vary the free parameter used in eqn (21), a, to determine what impact it has on the dynamics of swelling. While a is a fully empirical parameter, we may presume that it does represent material properties and therefore may vary from one type of starch to another. From Fig. 7(a) and (b) we see that, we are unable to successfully collapse all the data onto a master curve.

Fig. 7 (a) The swelling curves obtained when a is varied with all other parameters held constant and the swelling ratio is 1.69 (ν_∞ = 0.012), (b) the best possible rescaling and shifting of the time to collapse the curves in (a). (c) The relation between the optimal value of a and ν_∞ on log–log axes. (d) The swelling curves simulated when a is varied with ν_∞ to produce the best fit to the experimental data.

The inability to collapse the curves can be interpreted in at least two ways. First, it may suggest that the empirical function obtained in ref. 19 may not be applicable to all starches. Second, it may suggest that there is only one value of a that is physically realistic. However, recall that a is not based on a known material property and instead only used to smoothly vary the cross-link density as a function of temperature. Therefore, we elect to proceed in our analysis by adopting the latter interpretation and assume that a must be selected such that the swelling curve has the same qualitative shape as eqn (25).

Next we recall that in Section 4.2 we found that the optimal value of a was different for different values of ν_∞. From this we assume that a is a function of ν_∞ alone and determine this relationship by obtaining the best fit for a at four different values of ν_∞ as shown in Fig. 7(c). From the results in Fig. 7(c), we find that the relationship between the optimal value of a and ν_∞ can be approximated by the following power-law function:

a = 0.341ν_∞^−0.0750 (28)

We see that the resulting swelling curves shown in Fig. 7(d) are nearly identical with only a very small time shift, if any, required to collapse them. Thus for the remaining simulations, we will not consider a to be an independent parameter and control its value using eqn (28).

To illustrate the impact of varying ν_∞, we have replotted the data from Fig. 7(d), without rescaling the granule radius, in Fig. 8(a). As expected, we see that the final granule radius is directly controlled by the final cross-link density. Plotting the value of ν_∞ as a function of the swelling ratio R_f/R₀ in Fig. 8(b), we find that it is also possible to approximate their relationship using another power-law function:

Fig. 8 (a) Swelling curves from Fig. 7(d) replotted to show the impact of ν_∞ on the dimensionless radius. (b) Swelling ratio vs. ν_∞ on log–log axes.

The final parameter that is involved in the cross-link density is the initial starch volume fraction ϕ₀. Fig. 9(a) shows the swelling curves when ϕ₀ changes from 0.5 to 0.65 when ν_∞ = 0.012. The results show only slight differences, revealing that ϕ₀ has little impact on the swelling process. Applying a time shift, these curves collapse very well, as presented in Fig. 9(b). When the same time shifting procedure is applied at different values of ν_∞, Fig. 9(c) shows how the required time shift varies with ϕ₀ in each case.

Fig. 9 (a) Swelling curves simulated for various values of ϕ₀ for a swelling ratio of 1.69 (ν_∞ = 0.012). (b) Illustration of the collapse of the curves from (a) by changing t_G. (c) Illustration of the optimal values of t_G needed to collapse the swelling curves for different values of ϕ₀ at different values of ν_∞ on log–log axes. (d) Data in (c) collapsed onto a single curve given by eqn (31).

From the data in Fig. 9(c), we find that t_G appears to depend linearly on ϕ₀ and, after examining the best fits to a line at different values of ν_∞, we see that the dependence of both parameters is well approximated by the following type of function:

where the following values produce the best fit: n₁ = 539, n₂ = 0.0430, k₁ = 217, and k₂ = 0.0860. However, it is important to remember that with more parameters one can always obtain an artificially better fit. Therefore, seeing that n₂ is quite small, we instead propose that the functional dependence can be simplified to the following function:where the following values produce the best fit: C₁ = 682, k₁ = 105, and k₂ = 0.212.

4.3.2 Effect of gelatinization temperature T_g and initial temperature T₀. In this section we examine the influence of the model parameters related to the temperature. First, we vary the gelatinization temperature T_g while holding the initial temperature constant at T₀ = 323 K. It is known that the gelatinization temperature of the starch will impact swelling, and indeed Fig. 10(a), shows that T_g results in changes to the swelling curve. By applying a time shift to the different curves, we see in Fig. 10(b) that the data can still be collapsed onto a single curve. This result is qualitatively unsurprising since the difference in time corresponds directly to the difference in temperature at which swelling begins because of the way heating is performed.

Fig. 10 (a) Swelling curves simulated for various values of T_g for a swelling ratio of 3.21 (ν_∞ = 0.00038) and initial temperature of 323 K. (b) Illustration of the collapse of the curves from (a) by changing t_G. (c) Illustration of the optimal values of t_G needed to collapse the swelling curves for different values of T_g at different values of ν_∞ on linear axes.

As in Section 4.3.1, we see that the time shift required to collapse the curves depends both on the gelatinization temperature and the cross-link density as shown in Fig. 10(c). Quantitatively, we find that the relation between t_G and T_g can be well described by a second order polynomial where the prefactors each may depend on ν_∞ as follows:

where the following values produce the best fit: m₁ = 204, m₂ = 0.0905, m₃ = 35.1 K⁻¹, m₄ = 0.0319, m₅ = 0.284 K⁻², and m₆ = 0.0830. However, again we wish to avoid excessive use of free parameters to fit the data, and we find that we can simplify the function to the following form and still describe the data reasonably well:where the following values produce the best fit: C₂ = 182, C₃ = 0.0914, and k₃ = 37.1 K⁻¹.

Next, we vary the initial temperature T₀ while holding the gelatinization temperature at T_g = 339 K. Fig. 11(a) reveals the influence of initial temperature T₀ on the swelling process. Similar to T_g, changes in T₀ also only require time shifts to collapse onto a single curve in as shown in Fig. 11(b). Fig. 11(c) also reveals a very similar relationship between t_G and T_g (refer to Fig. 10(c)).

Fig. 11 (a) Swelling curves simulated for various values of T₀ for a swelling ratio of 1.69 (ν_∞ = 0.012) and gelatinization temperature of 339 K. (b) Illustration of the collapse of the curves from (a) by changing t_G. (c) Illustration of the optimal values of t_G needed to collapse the swelling curves for different values of T₀ at different values of ν_∞ on linear axes.

Therefore, we again find that the same type of second order polynomial in eqn (34) fits the data very well:

with m₁ = 189, m₂ = 0.0956, m₃ = 25.9 K⁻¹, m₄ = 0.00803, m₅ = 0.317 K⁻², and m₆ = 0.0868 providing the best fit. For the same reasons as above, we also simplify the function to the following form:with C₂ = 164, C₃ = 0.0969, and k₃ = 38.4 K⁻¹.

From the best fit parameters to eqn (33) and (35), we observe that the dependence of t_G is very similar for the same difference between T_g and T₀. For this reason, both Fig. 10(c) and 11(c) were plotted versus the temperature difference T_g − T₀, even though in each case one of the two values was held constant. Given the high degree of similarity when varying the two parameters independently, we combined all of the data from both cases and fit them to the following function:

When all of the data in both plots are combined, the following parameters provide the best fit: C₂ = 153, C₃ = 0.117, and k₃ = 36.5 K⁻¹.

4.3.3 Effect of other parameters: ΔH, χ₀, and c. Aside from the parameters directly related to gelatinization, there are three other material parameters that appear in the dimensionless equations: χ₀ and ΔH in the Flory–Huggins parameter expression (eqn (3)), as well as c in the diffusion coefficient (eqn (8)). From Fig. 12(a) we see that changing ΔH by several orders of magnitude reveals that a larger ΔH makes swelling faster, since more enthalpy change is involved in the interaction between starch and water, but this effect is very modest. Similarly, varying χ₀ and c within the range of physically meaningful values we see almost no change in the swelling curve (see Fig. 12(b) and (c)). Thus, we may safely say that these three parameters are unrelated to the empirical factors (t_G and k_G) used to collapse the data onto a master curve.

Fig. 12 Swelling curves simulated for various values of ΔH (a), χ₀ (b), and c (c) for a swelling ratio of 1.69 (ν_∞ = 0.012).

4.4 Summary of parametric analysis

In Section 4.3 we hypothesized that the empirical parameters that can be used to scale and shift the swelling data to produce a universal curve could be predicted from material properties and written as some function of the material parameters:

t_G = f(ν_∞,ϕ₀,a,T_g,T₀,ΔH,χ₀,c) (37)

and

k_G = g(ν_∞,ϕ₀,a,T_g,T₀,ΔH,χ₀,c). (38)

From our results we found that for a given value of ν_∞, only one value of a can produce a curve with the same qualitative behavior seen in experiments. Thus, we determined that a is not an independent parameter and we can eliminate a dependence on a from both functions. We also found that neither t_G nor k_G depend on ΔH, χ₀, or c, so those can also be eliminated from both functions.

Looking at the remaining parameters we see that none of them influence k_G, suggesting that the timescale of gelatinization is entirely determined by the diffusive time scale. This is interesting, because it suggests that by fitting experimental data, one can obtain an estimate of the diffusion coefficient that is relevant to gelatinization. In turn, this will help determine if assumptions in our model about the diffusivity are correct. This will be the subject of future work.

Further, t_G was found to be adequately described by a linear dependence on ϕ₀ and (T_g − T₀) with a slope and an intercept, respectively, that depend on ν_∞, which in turn depends on the swelling ratio. Thus, we can now revise eqn (37) and (38) to be:

and

k_G = 1. (40)

We can test the accuracy of eqn (39) by combining all of the simulation data and fitting as shown in Fig. 13. We find that the best fit results from the following values of the parameters: k₁ = 280, k₂ = 0.732, and k₃ = 37.8 K⁻¹, with a maximum error of 6.6% and an average error of 1.6%.

Fig. 13 The accuracy of eqn (39) by comparing its prediction for t_G with those from simulation.

Finally, we return to the functional dependence of the cross-link density on temperature. We found that there appears to be a power-law dependence of the parameter a on ν_∞ (eqn (28)), and hence on the swelling ratio. This means that we can rewrite the modified logistic function as follows:

From the previous fits to the data, we found the following values of the fitted parameters: k₄ = 0.201, k₅ = 5.38, k₆ = 0.409 K⁻¹, and k₇ = 0.900.

4.5 Influence of heating rate and heating function

As a final point of discussion, we wish to point out that the specific shape of the swelling curve must depend on the way in which the temperature of the system increases with time. This is a noteworthy point because our comparisons were done with experimental swelling data obtained where the temperature was increased quadratically with time, and other researchers might instead increase the temperature linearly with time. In Fig. 14, we have compared the swelling curves from simulations in which the temperature was increased linearly with time as T = 323 + bt at different rates, b in K s⁻¹, against the original parabolic heating function used in ref. 19. In each case, two different values of the swelling ratio were used.

Fig. 14 Comparison of the effect of the heating rate on the functional form of the swelling curve. Note that the parameter b corresponds to the linear heating ramp T = 323 + bt, while the parabolic profile refers to heating using eqn (22). The results from two different swelling ratios are shown: (a) R_f/R₀ = 3.21 and (b) R_f/R₀ = 1.69.

As expected, using different heating functions and rates has a large impact on the shape of the swelling curve. Even a slight difference in the heating function will result in different shapes of the swelling curve and the ability to collapse the data onto the function found in ref. 19 is clearly a special case. Nevertheless, the convenience of the functional dependence found in ref. 19 has allowed us to make substantial progress in understanding which material properties will influence the swelling process. Further, if other researchers utilize the same heating function, the present analysis will likely enable them to extract material properties directly from swelling data.

5 Conclusion

In this work we used numerical simulations based on a first-principles model of the swelling of starch granules to help understand a previously published experimental result¹⁹ in which a universal curve was seen for the granule size as a function of time. While validating the governing equations for the simulations, we found that the standard Flory–Rehner theory used in previous work²² was ill-suited to simulating the initial stages of starch swelling near the gelatinization temperature. However, we concluded that with a minor modification, in which we allow the fraction of cross-linked chains in the granule to decrease as temperature increases, that the Flory–Rehner theory can be used to successfully predict the swelling process of legume-starch granules.

Using the modified theory and systematically varying the model parameters, we found that the granule size as a function of time from numerical simulations can also be collapsed onto the same master curve found in experiments. We presented a simple relationship between the empirical scaling and shifting parameters used in experiments and the material properties used in the modified theory. Of particular note is that once the time is non-dimensionalized by the diffusion time scale, no further rescaling is needed to collapse the swelling data. Changes to all other model parameters only require a shift in the time, to collapse the swelling data.

In summary, this is the first study to demonstrate that the swelling of individual starch granules during gelatinization can be described by a universal function that is generated from first-principles. This will enable future experiments with granule swelling to extract material properties such as the diffusivity, gelatinization temperature, and initial starch volume fraction for individual starch granules. Further, the framework we developed will lend itself well to exploring the utility of using more sophisticated variants of the Flory theory in future work.^{32,34,36,40,41} Finally, this work will enable more accurate predictions of the behavior of starches during advanced processing operations, in which the gelatinization of starch is critical to the functional properties of the final product.

Author contributions

Botong Li: formal analysis; funding acquisition; software and validation; original draft. Lanxin Mo: experiment. Vivek Narsimhan: formal analysis; funding acquisition; methodology; review and editing. Ganesan Narsimhan: formal analysis; funding acquisition; methodology; review and editing. John M. Frostad: conceptualization; formal analysis; funding acquisition; methodology; project administration; original draft; review and editing.

Data availability

The datasets supporting this article have been uploaded as part of ESI.†

Conflicts of interest

The authors declare that they do not have any conflicts of interest.

Acknowledgements

The author B. Li is grateful for the China Scholarship Council (202106465015). Authors G. Narsimhan, V. Narsimhan, and J. M. Frostad would like to acknowledge support from USDA-NIFA (grant 2024-67017-42557). J. M. Frostad is also grateful for support through funding from the Canada Foundation of Innovation John R. Evans Leaders Fund, the British Columbia Knowledge Development Fund, and the Discovery Grant Program of the Natural Sciences and Engineering Research Council of Canada (NSERC) (Grant No. RGPIN-2017-04466).

References

1. S. Wang and L. Copeland, Molecular disassembly of starch granules during gelatinization and its effect on starch digestibility: a review, Food Funct., 2013, 4(11), 1564–1580,  10.1039/C3FO60258C.
2. L. Roy, Whistler and James N. BeMiller, Carbohydrate chemistry for food scientists, Eagan Press, 1997.
3. R. Hoover, The Impact of Heat-Moisture Treatment on Molecular Structures and Properties of Starches Isolated from Different Botanical Sources, Crit. Rev. Food Sci., 2010, 50(9), 835–847,  DOI:10.1080/10408390903001735.
4. S. Dhital, et al., Rice starch granule amylolysis – Differentiating effects of particle size, morphology, thermal properties and crystalline polymorph, Carbohydr. Polym., 2015, 115, 305–316,  DOI:10.1016/j.carbpol.2014.08.091.
5. P. Chen, et al., Morphology and Microstructure of Maize Starches with Different Amylose/Amylopectin Content, Starch-Stärke, 2006, 58(12), 611–615,  DOI:10.1590/1981-6723.0516.
6. N. Lindeboom, P. R. Chang and R. T. Tyler, Analytical, Biochemical and Physicochemical Aspects of Starch Granule Size, with Emphasis on Small Granule Starches: A Review, Starch-Stärke, 2004, 56(3–4), 89–99,  DOI:10.1002/star.200300218.
7. D. Lund and K. J. Lorenz, Influence of time, temperature, moisture, ingredients, and processing conditions on starch gelatinization, CRC Crit. Rev. Food Sci., 1984, 20(4), 249–273,  DOI:10.1080/10408398409527391.
8. Y. Liu, et al., Gelatinization behavior of starch: Reflecting beyond the endotherm measured by differential scanning calorimetry, Food Chem., 2019, 284, 53–59,  DOI:10.1016/j.foodchem.2019.01.095.
9. M. Beck, M. Jekle and T. Becker, Starch re-crystallization kinetics as a function of various cations, Starch-Stärke, 2011, 63(12), 792–800,  DOI:10.1002/star.201100071.
10. M. Schirmer, et al., Physicochemical and morphological characterization of different starches with variable amylose/amylopectin ratio, Food Hydrocolloids, 2013, 32(1), 52–63,  DOI:10.1016/j.foodhyd.2012.11.032.
11. M. Mariotti, et al., Effects of dispersing media and heating rates on pasting profiles of wheat and gluten-free samples in relation to their solvent retention capacities and mixing properties, LWT – Food Sci. Technol., 2016, 66, 201–210,  DOI:10.1016/j.lwt.2015.09.041.
12. M. Joshi, et al., Physicochemical and functional characteristics of lentil starch, Carbohydr. Polym., 2013, 92(2), 1484–1496,  DOI:10.1016/j.carbpol.2012.10.035.
13. M. S. Shang, et al., Effect of single and dual heat–moisture treatments on the gelatinization properties and crystalline structure of normal corn starch, Starch-Stärke, 2016, 68(11–12), 1196–1202,  DOI:10.1002/star.201500350.
14. T. J. Woodbury, et al., Mechanisms of the different effects of sucrose, glucose, fructose, and a glucose–fructose mixture on wheat starch gelatinization, pasting, and retrogradation, J. Food Sci., 2023, 88(1), 293–314,  DOI:10.1111/1750-3841.16414.
15. C. V. T. Rigueto, et al., Gelatin films from wastes: A review of production, characterization, and application trends in food preservation and agriculture, Food Res. Int., 2022, 162, 112114,  DOI:10.1016/j.foodres.2022.112114.
16. C. Li, Recent progress in understanding starch gelatinization – An important property determining food quality, Carbohydr. Polym., 2022, 293, 119735,  DOI:10.1016/j.carbpol.2022.119735.
17. G. P. Desam, et al., Prediction of swelling behavior of crosslinked maize starch suspensions, Carbohydr. Polym., 2018, 199, 331–340,  DOI:10.1016/j.carbpol.2018.07.020.
18. T. Zhang, et al., Gelatins as emulsifiers for oil-in-water emulsions: Extraction, chemical composition, molecular structure, and molecular modification, Trends Food Sci. Technol., 2020, 106, 113–131,  DOI:10.1016/j.tifs.2020.10.005.
19. L. Mo, J. Cheon and J. M. Frostad, Quantifying and modeling the gelatinization properties of individual pulse-starch granules by ParCS, Food Hydrocolloids, 2023, 135, 107896,  DOI:10.1016/j.foodhyd.2022.107896.
20. M. C. Allan, Q. D. Read and S. D. Johanningsmeier, Impact of sweetpotato starch structures, thermal properties, and granules sizes on sweetpotato fry textures, Food Hydrocolloids, 2023, 137, 108377,  DOI:10.1016/j.foodhyd.2022.108377.
21. A. Palanisamy, et al., Kinetic modelling of individual starch granules swelling, Food Struct., 2020, 26, 100150,  DOI:10.1016/j.foostr.2020.100150.
22. G. P. Desam, et al., A mechanistic model for swelling kinetics of waxy maize starch suspension, J. Food Eng., 2018, 222, 237–249,  DOI:10.1016/j.jfoodeng.2017.11.017.
23. P. J. Flory and J. Rehner, Statistical Mechanics of Cross-Linked Polymer Networks II. Swelling, J. Chem. Phys., 1943, 11(11), 521–526,  DOI:10.1063/1.1723792.
24. P. J. Flory, Statistical Mechanics of Swelling of Network Structures, J. Chem. Phys., 1950, 18(1), 108–111,  DOI:10.1063/1.1747424.
25. Y. Zhao and B. E. Eichinger, Study of solvent effects on the dilation modulus of poly(dimethylsiloxane), Macromolecules, 1992, 25(25), 6988–6995,  DOI:10.1021/ma00051a041.
26. Z.-Y. Lu and R. Hentschke, Swelling of model polymer networks with different crosslink densities: A computer simulation study, Phys. Rev. E:Stat., Nonlinear, Soft Matter Phys., 2002, 66(4), 041803,  DOI:10.1103/PhysRevE.66.041803.
27. C. Du and R. J. Hill, Swelling and Dissolution Transitions of DNA- and “Bis”-Cross-Linked Polyacrylamide, Macromolecules, 2020, 53(12), 4845–4854,  DOI:10.1021/acs.macromol.0c00471.
28. J. Li, et al., Methodology to predict the time-dependent storage modulus of starch suspensions during heating, Food Hydrocolloids, 2021, 113, 106463,  DOI:10.1016/j.foodhyd.2020.106463.
29. I. D. Evans and D. R. Haisman, The Effect of Solutes on the Gelatinization Temperature Range of Potato Starch, Starch-Stärke, 1982, 34(7), 224–231,  DOI:10.1002/star.19820340704.
30. M. Milene, et al., Physical and chemical properties of starch and flour from different common bean (Phaseolus vulgaris L.) cultivars, Braz. J. Food Technol., 2016, 19, e2016005.
31. D. T. Do, et al., A novel apparatus for time-lapse optical microscopy of gelatinisation and digestion of starch inside plant cells, Food Hydrocolloids, 2020, 104, 105551,  DOI:10.1016/j.foodhyd.2019.105551.
32. R. G. M. van der Sman and M. B. J. Meinders, Prediction of the state diagram of starch water mixtures using the Flory–Huggins free volume theory, Soft Matter, 2011, 7(2), 429–442,  10.1039/C0SM00280A.
33. R. G. M. van der Sman and M. B. J. Meinders, Moisture diffusivity in food materials, Food Chem., 2013, 138(2), 1265–1274,  DOI:10.1016/j.foodchem.2012.10.062.
34. R. G. M. van der Sman, Predicting the solubility of mixtures of sugars and their replacers using the Flory–Huggins theory, Food Funct., 2017, 8(1), 360–371,  10.1039/C6FO01497F.
35. R. G. M. van der Sman, Effects of viscoelasticity on moisture sorption of maltodextrins, Food Hydrocolloids, 2023, 139, 108481,  DOI:10.1016/j.foodhyd.2023.108481.
36. R. G. M. van der Sman, M. Curatolo and L. Teresi, Pore development in viscoelastic foods during drying, Soft Matter, 2024, 20(26), 5183–5194,  10.1039/D4SM00201F.
37. A. Briffaz, et al., Modelling starch phase transitions and water uptake of rice kernels during cooking, J. Cereal Sci., 2013, 58(3), 387–392,  DOI:10.1016/j.jcs.2013.08.001.
38. M. Alan and R. W. Hausslein, Elastic factors controlling sorption and transport properties of polyethylene, J. Polym. Sci., Part C, 1965, 10(1), 61–86.
39. M. Alan and R. W. Hausslein, Elastic constraints in solvent swollen polyethylene, J. Polym. Sci., Polym. Lett. Ed., 1965, 3(1), 61–62.
40. R. G. M. van der Sman, Biopolymer gel swelling analysed with scaling laws and Flory–Rehner theory, Food Hydrocolloids, 2015, 48, 94–101,  DOI:10.1016/j.foodhyd.2015.01.025.
41. W. Hong, Z. Liu and Z. Suo, Inhomogeneous swelling of a gel in equilibrium with a solvent and mechanical load, Int. J. Solids Struct., 2009, 46(17), 3282–3289,  DOI:10.1016/j.ijsolstr.2009.04.022.
42. C. Cai and C. Wei, In situ observation of crystallinity disruption patterns during starch gelatinization, Carbohydr. Polym., 2013, 92(1), 469–478,  DOI:10.1016/j.carbpol.2012.09.073.
43. W. David, et al., Characterizing internal cavity modulation of corn starch microcapsules, Heliyon, 2020, 6(10), e05294.
44. A. J. Easteal, W. E. Price and L. A. Woolf, Diaphragm cell for hightemperature diffusion measurements. Tracer Diffusion coefficients for water to 363 K, J. Chem. Soc., Faraday Trans. 1, 1989, 85(5), 1091–1097,  10.1039/F19898501091.
45. M. Matyka, A. Khalili and Z. Koza, Tortuosity-porosity relation in porous media flow, Phys. Rev. E:Stat., Nonlinear, Soft Matter Phys., 2008, 78(2), 026306,  DOI:10.1103/PhysRevE.78.026306.
46. P. J. Flory, Principles of polymer chemistry, Cornell University Press, 1953.
47. M. Rubinstein and R. H. Colby, Polymer Physics, Oxford University Press, USA, 2003.
48. J. J. M. Swinkels, Composition and Properties of Commercial Native Starches, Starch-Stärke, 1985, 37(1), 1–5,  DOI:10.1002/star.19850370102.
49. J. A. Nelder and R. Mead, A Simplex Method for Function Minimization, Comput. J., 1965, 7(4), 308–313,  DOI:10.1093/comjnl/7.4.308.
50. M. Wright, Direct search methods: Once scorned, now respectable. English (US), in Numerical analysis, ed. D. F. Griffiths and G. A. Watson, Addison-Wesley, 1996, pp. 191–208.

---

Footnote

† Electronic supplementary information (ESI) available. See DOI: https://doi.org/10.1039/d4sm00980k

---