Generalization of a lumped parameters model ... - Wiley Online Library

7 downloads 0 Views 484KB Size Report
Sep 25, 2017 - sis of the residues was made through the Shapiro-Wilk normality test, with 95% confidence. If accepted, the null hypothesis of this test sug-.
Received: 22 June 2017

|

Revised: 25 September 2017

|

Accepted: 11 October 2017

DOI: 10.1111/jfpe.12641

ORIGINAL ARTICLE

Generalization of a lumped parameters model using fractional derivatives applied to rice hydration Douglas Junior Nicolin1

|

Thaisa Carvalho Volpe Balbinoti2 |

Regina Maria Matos Jorge2,3 | Luiz Mario de Matos Jorge3 1 Chemical Engineering Department, Federal Technological University of Parana (UTFPR), Linha Santa Barbara, PO Box 135, Francisco Beltr~ao, 85601-970, Parana, Brazil 2

Abstract This study presents a generalization of a lumped parameters model through the use of fractional derivatives of arbitrary order. The conventional model has an analytical solution in terms of expo-

Chemical Engineering Department, Food Engineering Graduate Program, Federal University of Parana (UFPR), Cel. Francisco H. dos Santos Avenue, 210, Curitiba, 81531-970, Parana, Brazil

nential function. However, when applied to kinetic data on moisture absorption in rice grains, its

3

tion of moisture values in relation to time and process temperature.

Chemical Engineering Graduate Program (PEQ), State University of Maring a (UEM), Colombo Avenue, 5790, Bl. E-46, Maring a, 87020-900, Brazil Correspondence Douglas Junior Nicolin, Chemical Engineering Department, Federal Technological University of Parana (UTFPR), Linha Santa Barbara, PO Box 135, Francisco Beltr~ao, 85601-970, PR, Brazil. Email: [email protected] and Luiz Mario M. Jorge, Chemical Engineering Graduate Program (PEQ), State University of Maringa (UEM), Colombo Avenue, 5790, Bl. E-46, Maringa, 87020-900, Brazil. Email: [email protected]

performance was worse when compared to the proposed generalized model. The results were analyzed for four hydration temperatures. Despite the similarity of the mass transfer coefficients for both models, the proposed generalization provided a statistically significant improvement in predic-

Practical applications The model proposed in this paper is applicable in describing the hydration kinetics of rice which is an important step in parboiling process, for instance. The simplicity of the analytical solution of the model and the improvement of its prediction capability makes the presented approach useful to project of equipment intended to process rice.

1 | INTRODUCTION

2005), minimization of grain breakage (Igathinathane, Chattopadhyay, & Pordesimo, 2005), as well as an increase in shelf life (Dutta & Mahanta,

Rice is considered one of the principal cereals being produced and con-

2014).

sumed around the world. It is estimated that half of the world’s popula-

The parboiling process usually consists of three stages: soaking the

tion consumes rice on a daily basis (Mir, Bosco, Shah, & Mir, 2016; Oli,

grain until reaching saturation, gelatinization of the starch of the grain

Ward, Adhikari, & Torley, 2014). Given the great economic and nutri-

by heating of the hydrated grains, and drying of the product until an

tional importance of rice, studies dedicated to its various processing

appropriate moisture level for grinding or storage is reached. According

methods are relevant to the development of the final product quality.

to Pascual et al. (2013) and Thammapat, Meeso, and Siriamornpun

Parboiling is one of the processes in which the main stage is the

(2016), hydration has the aim of providing the rice with the necessary

absorption of moisture by rice grains. Numerous studies have confirmed

moisture to induce starch gelatinization. Gelatinization provides mor-

the beneficial changes that the grain undergoes when parboiled, such as

phological transformations in the physical structure of the food (Ali &

alteration of texture (Buggenhout, Brijs, Celus, & Delcour, 2013), altera-

Pandya, 1974; Sittipod & Shi, 2016) involving the elimination of preex-

tion of color (Lamberts, Brijs, Mohamed, Verhelst, & Delcour, 2006),

isting fissures in the grain (Chen & Chen, 1986).

alteration in nutritional content (Paiva et al., 2016; Setyaningsih, Saputro,

Given the importance of the soaking stage in the processing of

Barbero, Palma, & García Barroso, 2015; Storck, Silva, & Alves Fagundes,

parboiled rice, study of the water absorption kinetics over time is of

J Food Process Eng. 2018;41:e12641. https://doi.org/10.1111/jfpe.12641

wileyonlinelibrary.com/journal/jfpe

C 2017 Wiley Periodicals, Inc. V

|

1 of 7

2 of 7

|

NICOLIN

ET AL.

particular interest, so as to enable modeling and optimization of such

within 3 hr and subsequently at intervals of 1 hr to complete the tests.

process.

The moisture gain of the samples were evaluated according to AOAC

One way to mathematically model hydration kinetics is the proposi-

International (2000).

tion of lumped parameters models. In this methodology, moisture is considered to be distributed homogeneously inside the food, without the presence of moisture gradients inside the grains (Hangos & Cameron,

2.2 | Theory

2001). Such approach was efficiently used in hydration modeling of peas (Omoto et al., 2009) and soybean (Coutinho, Conceiç~ao, Omoto, Andrade, & Jorge, 2007; Coutinho, Omoto, Andrade, & Jorge, 2005; Nicolin, Neto, Paraíso, Jorge, and Jorge, 2015). On account of the simple nature of these models, the equation proposed by Omoto et al. (2009) enables easy practical application. However, given the fact that the analytical solution of these models has an exponential nature, it can fail in the description of kinetic curves that are nonexponential. An alternative for solving such situations would be the use of a fractional calculus, which generalizes the model and introduces improved performance. The fractional calculus has the characteristic of capturing behavior that traditional calculus often cannot, by considering operators such as

The generalized model of the present study was based on the model originally proposed by Omoto et al. (2009). In their study, the authors proposed a model in which the rate of moisture accumulation, inside spherical pea grains, was caused by the convective flow of moisture around the grains, which were submersed in water. The convective flow of water was given by NA 5KS Aðqeq 2qA Þ, where NA is the flow of water (g/cm2hr), A the surface area of the grain (cm2), qA water concentration in the grain (g/cm3), and qeq water concentration in the grain at the point of equilibrium (g/cm3). Therefore, the obtained mass balance generated a lumped parameters model (Equation 1), in which V designates the volume (cm3) and qA0 the initial moisture of the grain (g/cm3). dqA AKS 5 ðqeq 2qA Þ dt V

integrals and derivatives of arbitrary order (Oldham & Spanier, 1974; Podlubny, 1999). Depending on the fractional order, the parameters of

(1)

this kind of model may lose their physical reality, although generally

The initial condition used for the solution of the model considered

they supply a better representation of the hydration process. For

that the concentration of water in the grains was known at the begin-

kinetic models, for example, this mathematical technique can be carried

ning of the hydration process, such that qA ð0Þ5qA0 .

out considering a derivative of fractional order to represent the varia-

The original work considered spherical geometry for pea grains as

tion of the property to be described. In recent studies on the adsorp-

a hypothesis for the model. However, the grains of rice investigated in

tion kinetics of heavy metals, Friesen, Leitoles, Gonçalves, Lenzi, and

this study were considered cylindrical. Thus, the model for rice grains

jo, Lenzi, Silva, and Lenzi (2013) proLenzi (2015) and Gomes, Arau

becomes:

posed lumped parameters models of fractional order and obtained satisfactory results in the description of the experimental data analyzed. In this context, the present study has the objective of using fractional calculation to generalize the exponential model proposed by

dqA 5kðqeq 2qA Þ dt

(2)

qA ð0Þ5qA0

(3)

Omoto et al. (2009), and of evaluating the benefits of this generalization

ðh1RÞ , h being the length of the cylinder (cm) and R Where k5 2KS hR

on the performance of the original model, when applied to water

the radius of the cylinder (cm). The constant was calculated using the

absorption data of rice grains. Furthermore, it is the intention of this

equations for surface area and for cylinder volume, A52phR12pR2

study to analyze the behavior of the main parameters present in the

and V5pR2 h. A dimensional analysis of the constant k allowed to con-

conventional model and fractional model as a function of temperature at

clude that it has units of hour21 (hr21).

which the hydration process occurs. Omoto’s original equation and the proposed fractional order model were compared statistically in order to identify the best performance in describing the experimental data.

2 | MATERIALS AND METHODS

The model presented by Equations 2 and 3 has an analytical solution as follows: qA ðtÞ5qeq 1ðqA0 2qeq Þe2kt

(4)

The generalization proposed in the present study replaces the first derivative present in Equation 2 with a derivative of arbitrary order (or

2.1 | Obtaining experimental data In order to obtain moisture of rice as a function of time tests were performed with paddy rice variety SCS117CL (0.177 6 0.001 kg/kg dry basis; 1.207 6 0.063 mm radius; and 7.162 6 0.327 mm long) grown in

of fractional order) as presented below: d a qA 5kðqeq 2qA Þ dta

(5)

qA ð0Þ5qA0

(6)

Navegantes, SC, Brazil, belonging to the 2014/2015 harvest. For grain

Where a is the fractional order of the derivative (dimensionless).

moisturizing 1 kg of paddy rice was submitted in triplicate to the hydra-

For practical applications, the best definition for the fractional

tion process done with thermostatic bath (SOLAB/SL-143155/22/Bra-

derivative of order a is the Caputo derivative (Caputo, 1967). This defi-

zil) containing 4 L of distilled water. The test conditions were 35.5,

nition involves the initial values of the differential equation of fractional

45.5, 55.5, and 65.5 8C for 10 hr, in which samples of approximately

order, which makes it adequate for application to practical problems.

20 g (Shimadzu/Scale AY220/Japan) were collected every 30 min

Furthermore, the Caputo derivative has a known Laplace transform

NICOLIN

|

ET AL.

(Oldham & Spanier, 1974; Podlubny, 1999). The Caputo derivative has the following definition:

3 of 7

Where / is the objective function (g2/cm6) and Nexp the number of experimental data (dimensionless).

ðt ðmÞ qA ðsÞ da qA 1 5 ds Cða2mÞ ðt2sÞa1m11 dta

The adjustment of the models to the experimental data was qualifor m21