Dynamic Water Adsorption Characteristics of Distillers Dried Grains

0 downloads 0 Views 350KB Size Report
Distillers dried grains with solubles (DDGS) is one of the coproducts obtained from dry-grind ethanol manufacturing. As the ethanol industry is growing ...
Dynamic Water Adsorption Characteristics of Distillers Dried Grains with Solubles (DDGS) V. Ganesan,1 K. A. Rosentrater,2,3 and K. Muthukumarappan4 ABSTRACT

Cereal Chem. 84(6):548–555

Distillers dried grains with solubles (DDGS) is one of the coproducts obtained from dry-grind ethanol manufacturing. As the ethanol industry is growing exponentially, the production of DDGS has been significantly increasing as well. To optimize the use of DDGS, it has to be economically transported from one part of the country to other parts, and stored efficiently. But DDGS has some flow issues, which often makes storage and transportation very problematic. So the objective of this study was to investigate the dynamic water adsorption characteristics of DDGS with four soluble levels at four temperatures and four relative humidities. Three mathematical models were then used to fit the adsorption data

(Peleg, Pilosof, and Singh-Kulshrestha). As there was no model available for describing the water adsorption characteristics of DDGS with varying soluble levels at various temperature and relative humidity conditions, a new comprehensive model was developed. The new model, GanesanRosentrater-Muthu (GRM) model, encompassed soluble level, temperature, and relative humidity effects, along with time and moisture content. The GRM model (R2 = 0.94; F =16503.90) provided a good description of DDGS water adsorption behavior and can be used to predict the dynamic adsorption of water in DDGS for a broad range of storage conditions.

Distillers dried grains with solubles (DDGS) is one of the coproducts obtained from dry-grind ethanol milling. Distillers grains are energy dense because they often contain relatively high fat contents and are nearly devoid of starch. DDGS typically contains 86–93% (db) dry matter, 26–34% (db) crude protein, and 3–13% (db) fat (Rosentrater and Muthukumarappan 2006). It is an excellent source of energy and protein for animal feeds. As the U.S. ethanol industry has been growing tremendously during the past decade, the production of DDGS has been significantly increasing as well, and it is forecasted to be over 13 million metric tons in the fiscal year 2006-2007 (AAFC 2006). To effectively utilize DDGS in the domestic and international markets, it has to be transported greater distances. During transportation, DDGS can be subjected to varying temperature and humidity levels. DDGS is stored in bins and silos until use; sometimes this may be a relatively long time. DDGS storage and transportation can become problematic as it often hardens inside the storage structures and railcars. This can lead to severe damage when these vessels are unloaded. DDGS storage and flow behavior depend on physical and chemical properties of the material itself, as well as external conditions such as air temperature and relative humidity. At times, even small differences in these factors can cause caking and bridging (i.e., poor flowability) in DDGS. The transport and storage of DDGS, and the resulting flow behavior under the various storage conditions are important considerations, especially in the light of the growth of the industry. Most organic granular materials are hygroscopic in nature, so they will gain or lose moisture when they are exposed to various humidity conditions. This can lead to physical and chemical changes in the material itself. Each material absorbs moisture at different rates under different temperature and humidity conditions.

Hence, there have been many studies focused on characterizing the sorption process of biomaterials such as chickpeas (Turhan et al 2002; Wood and Harden 2006), dasheen leaves (Maharaj and Sankat 2000), yellow field peas (Wang et al 2003), corn kernels (Muthukumarappan and Gunasekaran 1990), canola (Thakor et al 1995), wheat and barley (Tagawa et al 2003), African breadfruit seeds (Shittu et al 2004), pearl millet and sorghum (Badau et al 2005). Moisture sorption is coupled with increased cohesiveness, chiefly because of interparticle liquid bridge formation. Moisture content is a significant variable that affects the cohesive strength and arching ability of bulk materials (Johanson 1978). As the moisture content of granular solids increases, the adhesion (Craik and Miller 1958) and cohesion (Moreyra and Peleg 1981) tend to increase. Even a small change in moisture content will significantly affect the frictional properties (e.g., wall friction angle, internal angle of friction) of a material (Marinelli and Carson 1992). A few empirical equations are available to model the water adsorption characteristics of food and biological materials. Of those models, the Peleg model (Peleg 1988) is the simplest and it has been applied to a variety of materials by many researchers (Sopade and Obekpa 1990; Sopade et al 1992, 1994; Sopade and Okonmah 1993; Maharaj and Sankat 2000; Turhan et al 2002; Shittu et al 2004; Badau et al 2005). Pilosof et al (1985) developed a model for dehydrated food powders and Maia and Cal-Vidal (1994) evaluated its applicability to water uptake by citrus juice powder under different relative humidity conditions. To date, reports evaluating its applicability to food grains are few. Also, the Singh-Kulshrestha model has received less attention but it has been used to model water absorption data for soybean and pigeon pea (Singh and Kulshrestha 1987). Understanding water adsorption in DDGS at various temperatures and relative humidities is of practical importance to the industry to optimize storage and transportation conditions. However, to date, there are no published reports on the water adsorption characteristics of DDGS. Hence, the objectives of this study were to 1) generate water adsorption data for DDGS with four soluble levels (10, 15, 20, and 25% db) at four practical temperatures (10, 20, 30, and 40°C) and four relative humidity levels (60, 70, 80, and 90%) using appropriate saturated salt solutions; 2) evaluate the applicability of three empirical models for the prediction of moisture content of DDGS in terms of time and initial moisture content; 3) develop a new adsorption model incorporating soluble level, temperature, and relative humidity effects in addition to time and initial moisture content.

1 Graduate

research assistant, Department of Agricultural and Biosystems Engineering, South Dakota State University, Brookings, SD. 2 Agricultural and Bioprocess Engineer, USDA-ARS, North Central Agricultural Research Laboratory, Brookings, SD. Mention of a trade name, propriety product or specific equipment does not constitute a guarantee or warranty by the United States Department of Agriculture and does not imply approval of a product to the exclusion of others that may be suitable. 3 Corresponding author. Phone: 604-693-3241. Fax: 605-693-5240. E-mail address: [email protected] 4 Professor, Department of Agricultural and Biosystems Engineering, South Dakota State University, Brookings, SD. doi:10.1094 / CCHEM-84-6-0548 © 2007 AACC International, Inc.

548

CEREAL CHEMISTRY

MATERIALS AND METHODS

Singh-Kulshrestha model (Singh and Kulshrestha 1987)

Sample Collection and Preparation Samples of condensed distillers solubles (CDS) and distillers dried grains (DDG) were obtained from a commercial ethanol plant (Dakota Ethanol, Wentworth, SD) and were stored in sealed plastic buckets (the DDG at room temperature and the CDS at 4 ± 1°C) until needed. The soluble content of both the DDG and the CDS were determined, and the DDGS with four soluble levels (10, 15, 20, and 25% db) were prepared using the methodology developed by Ganesan et al (2006a). Additionally, the moisture content of the DDG and CDS were determined using Approved Method 44-19 (AACC International 2000). Initially DDG had ≈10% (db) solubles and ≈6.1% (db) moisture content. To obtain the water adsorption data over a wide range of temperatures and relative humidities, it was essential that experimental samples had constant initial moisture content. So, before experimentation, the DDGS samples that had been prepared with the various soluble contents were dried at 100°C to a moisture content of 6.1% (db) using an air oven. The time for drying DDGS was 15–20 min. Water adsorption characteristics and the equilibrium moisture contents of the DDGS were determined at four constant temperatures (10, 20, 30, and 40°C) and four humidity levels (60, 70, 80, and 90%) using the static gravimetric method. Saturated salt solutions were used to maintain a constant relative humidity in each environment. The equilibrium relative humidity (ERH) values of the saturated salt solutions (Table I) were obtained from the reported data (Spencer 1926; Hodgman 1954, 1955; Greenspan 1977) and were verified using thermo hygrometers (model 800017, Sper Scientific). In-depth details about the experimental set up and procedures are found in Ganesan et al (2006b). The DDGS samples were exposed to different temperature and humidity conditions and weighed every 1 hr during an 8-hr period, and then at constant intervals (every 24 hr in the first week; every 48 hr in the second week; every 72 hr in the third week; and every 168 hr in rest of the weeks) for a period of five to six weeks. The experiment was stopped when the difference between successive weighings was ≤0.001 g. Triplicates were measured for each treatment. Water Adsorption Models In this study, the applicability of three common empirical models that describe moisture adsorption over time were evaluated for DDGS: Peleg model, Pilosof model, and Singh-Kulshrestha model. Peleg model (Peleg 1988) M = M0 +

t A + Bt

(1)

where M is the moisture content after time t (% db), M0 is the initial moisture content (% db), A is the water absorption rate constant (hr/%db), and B is the characteristic water absorption constant (1/% db). Pilosof model (Pilosof et al 1985) M =

At B+t

(2)

where M is the total moisture uptake at time t (g/100 g of solid), A is the water absorption capacity (g/g of dry solid), and B is the half life (hr).

M − M0 Bt = A − M 0 Bt + 1

(3)

where M is the moisture content at any time t (g/g of dry solid), M0 is the initial moisture content (g/g, db), A is the water absorption capacity (g/g of dry solid), and B is the water absorption rate constant (1/hr). Model constants were determined by iterative nonlinear regression on the collected data using PROC NLIN of SAS v.9.1 (SAS Institute, Cary, NC). The nonlinear regression procedure used the Gauss-Newton method and minimized the sum of squares of deviation between the observed and predicted moisture adsorption data to resolve the models. The goodness-of-fit between the experimentally observed and predicted amount of water absorbed was evaluated and compared using the standard error of the mean (SEM), the root mean square deviation (RMSD), the coefficient of determination (R2), and F-statistic values

where Mi and Mˆ i are the experimentally observed and predicted amount of absorbed water respectively; n is the number of data points; and DF is the degrees of freedom (number of data points minus the number of coefficients in each model). PROC NLIN directly provides the F-statistic and R2 values (model sum of squares/total sum of squares). Many researchers have used other criteria to test model performance. Some have used Chi square, R2, and RMSD (Sopade et al 1992; Sacchetti et al 2003; Shittu et al 2004; Badau et al 2005), and some have used mean relative error (MRE), R2, and P values (Maharaj and Sankat 2000; Turhan et al 2002; Wang et al 2003) for model performance analysis. In this study, a model was considered “good”, when it produced relatively small SEM, RMSD, and high F-statistic and R2 values. RESULTS AND DISCUSSION Water Adsorption Behavior The water adsorption behavior of DDGS with four soluble (10, 15, 20, and 25% db) levels at four temperatures (10, 20, 30, and 40°C) and four relative humidity (60, 70, 80, and 90%) levels are shown in Figs. 1 and 2. The water adsorption trends followed an asymptotic growth pattern for all treatment combinations. Similar patterns were reported for cereals, legumes, and seeds (Sopade et al 1992, 1994; Hung et al 1993; Turhan et al 2002; Shittu et al 2004; Badau et al 2005). The water adsorption of DDGS for all the temperature and humidity combinations increased rapidly in the first 24 hr, and then slowly approached steady-state, which was typically reached at 300–500 hr. It was observed that the moisture adsorption rate was very high initially and became slower in the later stages. Similar responses were found for cereals grains and legumes (Engels et al 1986, 1987; Sopade and Obekpa 1990; Sopade et al 1992, 1994; Sopade and Okonmah 1993). The water adsorption capacity of DDGS with a soluble level of 25% (db) at 40°C and 90% RH was greater compared with the other treatment

TABLE I Saturated Salt Solutions Used for Producing Control Environments (actual measured ERHa) RH (decimal) 0.60 0.70 0.80 0.90 a

10°C Sodium bromide (0.622) Lithium acetate dihydrate (0.720) Ammonium sulfate (0.821) Strontium nitrate (0.906)

20°C Sodium bromide (0.591) Lithium acetate dihydrate (0.700) Ammonium chloride (0.792) Barium chloride (0.910)

30°C Sodium bromide (0.560) Strontium chloride (0.691) Potassium bromide (0.803) Barium chloride (0.890)

40°C Sodium bromide (0.570) Strontium chloride (0.680) Potassium bromide (0.794) Potassium nitrate (0.890)

Equilibrium relative humidity values. Vol. 84, No. 6, 2007

549

combinations. Sopade et al (1992) observed that the water absorption of cereal grains increased with increasing temperature. Large differences were found at higher equilibrium relative humidity (ERH) values, often because of the inherent difficulty in determining equilibrium moisture content (EMC) at these elevated humidity levels. Similar issues were reported by Maroulis et al (1988) and Yu et al (1999).

Water Adsorption Modeling The Peleg model, Pilosof model, and Singh-Kulshrestha model were used to fit the experimentally observed water adsorption data. The Singh-Kulshrestha model did not fit the experimental data well, as the model coefficients did not converge with the iterative nonlinear regression procedure. So the model was not considered for further evaluation. Tables II and III show parameters obtained

Fig. 1. Water adsorption characteristics of DDGS with various soluble levels at 10 and 20°C, and 60, 70, 80, 90% RH conditions (± 1 SD). 550

CEREAL CHEMISTRY

from nonlinear regression analysis using the Peleg and Pilosof models. The R2 values varied from 0.71 to 0.99 and 0.95 to 0.99 for the Peleg and Pilosof models, respectively. Both models produced high F-statistic and R2 values and small error terms (SEM, RMSD). Overall, the Peleg model produced higher Fstatistic values and smaller error terms than did the Pilosof model, except for few treatment combinations. The DDGS did show some

instability at 40°C, and a similar effect was reported for African breadfruit seeds at higher temperatures (Shittu et al 2004). Both the Peleg and Pilosof models performed and fit the DDGS water adsorption data well. In commercial ethanol plants, CDS (often referred to as “syrup”) is added back to the wet-distillers grains, the combination of which is then dried to produce DDGS.

Fig. 2. Water adsorption characteristics of DDGS with various soluble levels at 30 and 40°C, and 60, 70, 80, 90% RH conditions (± 1 SD). Vol. 84, No. 6, 2007

551

The amount of CDS added will vary for each production cycle and is often not a systematic process. Considering this, the soluble level was subsequently incorporated into the study as a key parameter in the Peleg equation, and this modified equation was then evaluated for performance.

Modified Peleg Equation (MPE1) (6)

where M is moisture content after time t (% db), M0 is the initial

TABLE II Regression Parameters Obtained for DDGS with 10 to 25% (db) Solubles at Various Temperature and Relative Humidity Conditions Using Peleg’s Equation Soluble Level (% db) 10

Temp (°C) 10

20

30

40

15

10

20

30

40

20

10

20

30

40

25

10

20

30

40

552

CEREAL CHEMISTRY

Relative Humidity (%) 60 70 80 90 60 70 80 90 60 70 80 90 60 70 80 90 60 70 80 90 60 70 80 90 60 70 80 90 60 70 80 90 60 70 80 90 60 70 80 90 60 70 80 90 60 70 80 90 60 70 80 90 60 70 80 90 60 70 80 90 60 70 80 90

Constant A (hr/%db) 1.45 1.57 1.70 1.45 1.52 1.70 1.05 0.92 1.54 1.46 0.69 0.63 0.25 0.55 –1.98 0.36 1.16 1.37 1.15 1.63 0.64 0.82 0.78 0.74 0.50 0.40 0.33 0.38 0.15 –1.74 0.38 0.79 0.72 0.78 0.86 1.17 0.89 0.59 0.87 0.73 0.20 0.33 0.29 0.27 –4.16 0.32 0.29 0.45 1.05 1.23 1.29 1.17 0.55 0.51 0.66 0.53 0.27 0.34 0.28 0.34 0.13 0.28 –1.06 0.77

Constant B (1/%db) 0.08 0.07 0.05 0.03 0.25 0.11 0.06 0.04 0.43 0.09 0.07 0.04 0.25 0.07 0.12 0.03 0.12 0.08 0.05 0.03 0.06 0.08 0.03 0.03 0.07 0.04 0.04 0.03 0.16 0.11 0.05 0.01 0.08 0.06 0.04 0.02 0.08 0.04 0.02 0.02 0.11 0.05 0.04 0.02 0.23 0.04 0.04 0.01 0.08 0.05 0.01 0.02 0.04 0.03 0.02 0.02 0.11 0.05 0.04 0.02 0.09 0.04 0.07 0.01

F-Statistic × 104 1.89 5.11 2.24 0.89 1.12 1.66 3.45 1.73 7.72 0.76 0.60 3.43 0.96 0.48 0.01 0.77 6.38 3.94 3.62 1.56 5.47 10.60 1.58 9.54 0.51 0.67 0.32 0.80 0.69 0.01 0.39 0.20 2.63 3.29 2.90 1.21 3.69 8.57 1.55 2.04 2.50 1.51 2.60 3.91 0.01 0.65 0.78 0.34 4.14 3.90 0.67 1.05 0.10 9.64 1.65 4.53 2.02 14.84 7.48 3.78 0.25 1.65 0.01 0.07

R2 0.99 0.99 0.99 0.99 0.99 0.99 0.99 0.99 0.99 0.99 0.99 0.99 0.99 0.99 0.76 0.99 0.99 0.99 0.99 0.99 0.99 0.99 0.99 0.99 0.99 0.99 0.99 0.99 0.99 0.76 0.99 0.98 0.99 0.99 0.99 0.99 0.99 0.99 0.99 0.99 0.99 0.99 0.99 0.99 0.82 0.99 0.99 0.99 0.99 0.99 0.99 0.99 0.99 0.99 0.99 0.99 0.99 0.99 0.99 0.99 0.99 0.99 0.71 0.96

SEM 0.57 0.37 0.68 1.39 0.49 0.53 0.41 0.97 0.16 0.85 1.23 0.75 0.57 1.23 9.56 2.06 0.24 0.37 0.49 1.08 0.39 0.24 0.88 0.49 1.25 1.27 2.40 1.98 0.81 10.23 2.10 6.61 0.46 0.50 0.67 1.63 0.31 0.40 1.61 1.30 0.49 0.90 0.84 1.10 6.17 1.75 1.77 4.78 0.41 0.52 2.73 2.02 1.23 0.47 1.25 0.92 0.53 0.31 0.51 1.11 1.60 1.15 16.37 11.16

RMSD 0.07 0.05 0.08 0.17 0.06 0.06 0.05 0.11 0.02 0.10 0.14 0.08 0.07 0.15 1.14 0.24 0.03 0.04 0.06 0.12 0.05 0.03 0.09 0.06 0.16 0.14 0.30 0.25 0.10 1.20 0.21 0.83 0.06 0.07 0.09 0.21 0.04 0.05 0.20 0.17 0.06 0.12 0.10 0.14 0.75 0.21 0.21 0.61 0.05 0.06 0.34 0.25 0.16 0.06 0.17 0.12 0.07 0.04 0.06 0.14 0.21 0.14 1.94 1.45

moisture content (% db), S is soluble level (% db), A is time constant (hr), B is dimensionless constant, and C is water adsorption rate constant (hr/% db) The parameters obtained for this model are shown in Table IV. The model produced higher F-statistic and R2 values and smaller

error terms (SEM, RMSD). The R2 obtained varied from 0.96 to 0.99. The R2 and F-statistic values were smaller and error terms were higher when compared with the parameters obtained from the Peleg model (Eq. 1). But these values were not significantly different from the Peleg model. This was expected, however, as a

TABLE III Regression Parameters Obtained for DDGS with 10 to 25% (db) Solubles at Various Temperature and Relative Humidity Conditions Using Pilosof’s Equation Soluble Level (% db) 10

Temp (°C) 10

20

30

40

15

10

20

30

40

20

10

20

30

40

25

10

20

30

40

Relative Constant A (g/g Humidity (%) of dry solids) 60 70 80 90 60 70 80 90 60 70 80 90 60 70 80 90 60 70 80 90 60 70 80 90 60 70 80 90 60 70 80 90 60 70 80 90 60 70 80 90 60 70 80 90 60 70 80 90 60 70 80 90 60 70 80 90 60 70 80 90 60 70 80 90

17.88 19.67 25.24 35.40 9.73 14.24 20.11 30.82 8.30 15.66 20.68 32.92 10.14 19.66 23.54 42.28 13.67 17.82 24.64 41.92 21.30 17.62 38.62 42.00 19.23 28.35 31.43 43.53 12.50 25.09 27.53 83.34 17.40 22.75 31.03 55.83 16.75 31.06 63.83 55.34 15.16 27.50 30.40 53.79 15.38 29.64 33.17 83.99 18.19 24.13 75.23 65.94 31.04 40.99 58.78 60.04 15.49 27.73 33.18 58.17 17.16 32.35 37.34 126.70

Constant B (hr) 4.49 5.99 10.20 18.24 0.69 2.69 4.79 10.48 0.31 3.05 3.85 8.60 0.28 2.77 3.38 7.73 2.04 4.31 7.75 29.10 3.99 3.16 11.59 14.68 2.48 4.37 4.81 8.66 0.40 4.18 3.90 43.01 2.70 5.03 9.44 35.28 2.82 7.27 33.26 22.91 0.79 3.62 4.03 8.82 0.59 3.96 4.41 24.26 4.12 7.98 61.58 46.76 7.22 10.15 20.65 18.61 1.08 3.92 4.51 11.85 0.74 4.22 4.78 66.57

F-Statistic × 104

R2

0.14 0.15 0.12 0.12 0.12 0.14 0.13 0.26 0.11 0.14 0.20 0.40 0.15 0.14 0.15 0.31 0.12 0.12 0.15 0.11 0.25 0.20 0.13 0.29 0.18 0.17 0.21 0.34 0.19 0.14 0.20 0.11 0.16 0.17 0.18 0.16 0.09 0.29 0.24 0.30 0.28 0.27 0.47 0.70 0.22 0.29 0.33 0.16 0.22 0.21 0.16 0.18 0.32 0.37 0.18 0.37 0.25 0.47 0.64 0.55 0.14 0.48 0.39 0.05

0.98 0.98 0.98 0.97 0.97 0.98 0.98 0.99 0.97 0.98 0.98 0.99 0.98 0.98 0.98 0.99 0.98 0.98 0.98 0.97 0.99 0.99 0.98 0.99 0.98 0.99 0.99 0.99 0.98 0.98 0.98 0.97 0.98 0.98 0.98 0.98 0.98 0.99 0.99 0.99 0.99 0.99 0.99 1.00 0.99 0.99 0.99 0.98 0.99 0.99 0.98 0.98 0.99 0.99 0.99 0.99 0.99 0.99 1.00 0.99 0.98 0.99 0.99 0.95

SEM 2.07 2.13 2.83 3.77 1.49 1.81 2.09 2.47 1.36 2.00 2.12 2.21 1.44 2.25 2.80 3.23 1.71 2.07 2.41 3.95 1.81 1.73 3.03 2.78 2.11 2.52 2.96 3.05 1.55 3.23 2.93 8.75 1.88 2.16 2.68 4.49 1.93 2.19 4.11 3.36 1.45 2.12 1.96 2.60 1.73 2.63 2.72 6.93 1.77 2.23 5.61 4.87 2.17 2.43 3.79 3.20 1.48 1.72 1.74 2.90 2.14 2.12 2.81 13.52

RMSD 0.26 0.26 0.35 0.47 0.18 0.22 0.28 0.30 0.16 0.24 0.26 0.27 0.18 0.29 0.34 0.39 0.22 0.27 0.31 0.51 0.23 0.22 0.43 0.35 0.26 0.35 0.37 0.38 0.19 0.38 0.35 1.04 0.24 0.28 0.35 0.57 0.28 0.28 0.52 0.43 0.18 0.28 0.24 0.32 0.21 0.31 0.32 0.88 0.21 0.27 0.70 0.60 0.27 0.32 0.52 0.42 0.19 0.22 0.22 0.37 0.29 0.26 0.33 1.76

Vol. 84, No. 6, 2007

553

6.00, RMSD 0.10) than the original Peleg model (Eq. 1). Out of these three models (MPE 2, MPE 3, and MPE 4), the MPE 4 model performed the best and is called the Ganesan-RosentraterMuthu (GRM) model. The GRM model produced comparatively high R2 and F-statistic values (Table V). The F-statistic value produced by this model is nearly twice the F-statistic values produced by other models (MPE 2 and MPE 3). Also, the error terms are comparatively less than those of the MPE 2 and MPE 3 models. Moreover, the GRM model encompassed the effect of all soluble levels, temperatures, and relative humidity levels. The GRM model was superior to other models and is applicable over a wide range of environmental combinations. Thus the GRM model is recommended for characterizing dynamic DDGS water adsorption and should be further examined with DDGS from other ethanol plants.

greater number of varying data was included in the model equation. As discussed above, temperature and relative humidity can play a major role in the rate and amount of water absorbed by DDGS during transport and storage. Hence, the model (Eq. 6) was further modified to include the effects of temperature and relative humidity individually, as well as in combination. Thirty different combinations were iteratively tried; only the three models that performed well are discussed below. Modified Peleg Equation 2 (MPE 2) (7)

where M is the moisture content after time t (% db), M0 is the initial moisture content (% db), S is the soluble level (% db), T is the temperature level (°C), A is an empirical constant (% db), B is the time constant (hr), and C is an empirical constant (°C). Modified Peleg Equation 3 (MPE 3)

CONCLUSIONS Dynamic water adsorption characteristics of DDGS with four soluble levels (10, 15, 20, and 25% db) at four temperatures (10, 20, 30, and 40°C) and four relative humidities (60, 70, 80, and 90%) were studied. This was the first step taken toward quantifying the water adsorption characteristics of DDGS. DDGS water adsorption data followed an asymptotical growth pattern and showed higher experimental error at higher temperatures and relative humidity levels. The obtained water adsorption data were mathematically modeled using the three common empirical equations (Peleg, Pilosof, and Singh-Kulshrestha). The Peleg and Pilosof models gave better descriptions of DDGS water adsorption data. The Singh-Kulshrestha model did not perform well because the model parameters did not converge. The addition of CDS greatly affects the water adsorption characteristics of DDGS at various temperature and humidity conditions. So the Peleg equation was iteratively modified by including soluble, temperature, and relative humidity effects. The optimal model, the GRM model (R2 0.94; F 16503.90), encompassed soluble, temperature, and relative humidity effects, in addition to time and moisture content effects and is thus applicable over a variety of environmental conditions.

(8)

where M is the moisture content after time t (% db), M0 is the initial moisture content (% db), S is the soluble level (% db), R is the relative humidity level (%), A is an empirical constant (% db), B is a time constant (hr), and C is a dimensionless constant. Modified Peleg Equation 4 (MPE 4) (9)

where M is the moisture content after time t (% db), M0 is the initial moisture content (% db), S is the soluble level (% db), T is the temperature level (°C), R is the relative humidity level (%), A is a constant (% db), B is a time constant (hr), C is an empirical constant (°C), and D is a dimensionless constant. Table V shows the parameters obtained through modified Peleg equations (Eqs. 7–9). These models produced smaller R2 (0.94) and F-statistic (16503.90) values and higher error terms (SEM

TABLE IV Regression Parameters Obtained for DDGS at Various Temperature and Relative Humidity Conditions Using Modified Peleg Equation Temp. (°C) 10

Relative Humidity (%) 60 70 80 90 60 70 80 90 60 70 80 90 60 70 80 90

20

30

40

Constant A (-)

Constant B (1/hr)

–13.98 –10.65 64.65 –1.97 43.64 44.38 13.32 3.36 2.51 1.33 –0.32 0.35 2.24 1.22 –0.40 0.22

1.60 1.11 0.45 0.37 1.05 0.78 0.37 0.40 2.08 0.96 0.76 0.42 2.25 0.84 0.72 0.21

Constant C (hr/% db)

F-Statistic × 104

R2

SEM

RMSD

0.29 0.53 0.18 1.45 0.28 0.55 0.57 3.68 0.19 0.37 0.47 1.33 1.15 0.89 0.85 0.32

0.97 0.99 0.96 0.99 0.97 0.99 0.99 1.00 0.96 0.98 0.98 0.99 0.99 0.99 0.99 0.97

2.24 1.97 5.40 2.24 2.60 2.34 3.13 1.42 2.74 2.76 3.03 2.69 1.09 2.18 2.56 7.74

0.14 0.13 0.35 0.14 0.17 0.15 0.21 0.09 0.17 0.18 0.19 0.17 0.07 0.13 0.15 0.48

SEM

RMSD

1.93 1.81 –1.69 1.39 –1.24 –1.40 0.18 0.49 0.20 0.29 0.35 0.34 0.03 0.28 0.32 0.60

TABLE V Regression Parameters Obtained for DDGS Using Modified Peleg Equations Constants ID MPE 2 MPE 3 MPE 4a a

A (% db)

B (hr)

C (°C)

0.69 0.74 0.29

0.94 1.37 0.56

0.10 6.42 0.33

D (-)

F-Statistic × 104

R2

0.72

0.68 0.71 1.65

0.84 0.84 0.94

MPE 4 model has been termed the Ganesan-Rosentrater-Muthu or GRM model in this study.

554

CEREAL CHEMISTRY

10.20 10.02 6.00

0.16 0.16 0.10

The proposed GRM model can be used by ethanol producers to plausibly predict water adsorption of DDGS at any given time, specific temperature, and relative humidity conditions. Future work should investigate the sorption behavior of DDGS at lower relative humidity levels as well. The flowability of DDGS was not investigated in this study after exposing DDGS to various temperature and humidity levels. Thus, it is suggested to examine the flowability of DDGS at various levels of temperature and humidity. This would allow us to determine the optimum level of solubles, temperature, and RH for better flowability of DDGS during storage and transport. ACKNOWLEDGMENTS We would like to thank Dakota Ethanol (Wentworth, SD) which contributed DDG and CDS for this study. We would also like to thank the South Dakota Corn Utilization Council (SDCUC), the South Dakota Agricultural Experimental Station (AES), and the USDA-ARS for financial support. LITERATURE CITED AACC International. 2000. Approved Methods of the American Association of Cereal Chemists, 10th Ed. Method 44-19. The Association: St. Paul, MN. AAFC. 2006. Protein meal: Situation and outlook. Bi-weekly Bulletin of Agriculture and Agri-Food Canada 19(3):1-4. Badau, M. H., Nkama, I., and Jideani, I. A. 2005. Water-absorption characteristics of various pearl millet cultivars and sorghum grown in northern Nigeria. J. Food Process Eng. 28:282-298. Craik, D. J., and Miller, B. F. 1958. The flow properties of powders under humid conditions. J. Pharm. Pharmacol. 10:136-144. Engels, C., Hendrickx, M., De Saemblanx, S., De Gryze, I., and Tobback, P. 1986. Modelling water diffusion during long grain rice soaking. J. Food Eng. 5:55-73. Engels, C., Hendrickx, M., and Tobback, P. 1987. Limited multilayer absorption of brown, parboiled rice. Int. J. Food Sci. Technol. 22:219223. Ganesan, V., Rosentrater, K. A., and Muthukumarappan, K. 2006a. Methodology to determine soluble content in dry grind ethanol coproduct streams. Appl. Eng. Agric. 22:899-903. Ganesan, V., Muthukumarappan, K., and Rosentrater, K. A. 2006b. Sorption isotherm characteristics of distillers dried grains with solubles (DDGS). Paper No. 066165. ASABE: St. Joseph, MI. Greenspan, L. 1977. Humidity fixed points of binary saturated aqueous solutions. J. Res. Nat. Bur. Stand. A. Physics and Chem. 81A(1):89-102. Hodgman, C. D. 1955. Handbook of Chemistry and Physics. 2309-2310. CRC: Cleveland, OH. Hung, T. V., Liv, L. H., Black, R. G., and Trewhella, M. A. 1993. Water absorption in chickpea (C. arientum) and field pea (P. sativum) cultivars using Peleg’s model. J. Food Sci. 58:848-852. Johanson, J. R. 1978. Know your material—How to predict and use the properties of bulk solids. Chemical Engineering. Desk book issue: 9-17. Maharaj, V., and Sankat, C. K. 2000. The rehydration characteristics and quality of dehydrated dasheen leaves. Can. Agric. Eng. 42:81-85. Maia, M. C. A., and Cal-Vidal, J. 1994. Kinetics of water uptake by citrus juices in powder form. Int. J. Food Sci. Technol. 29:137-141.

Marinelli, J., and Carson, J. W. 1992. Solve solids flow problems in bins, hoppers and feeders. Chem. Eng. Progress. 88:22-28. Maroulis, Z. B., Tsami, M. E., and Marinos-Kouris, D. 1988. Application of the GAB model to the moisture sorption isotherm for dried fruits. J. Food Eng. 7:63-78. Moreyra, R., and Peleg, M. 1981. Effect of equilibrium water activity on bulk properties of selected food powders. J. Food Sci. 46:1918-1922. Muthukumarappan, K., and Gunasekaran, S. 1990. Vapor diffusivity and hygroscopic expansion of corn kernels during adsorption. Trans. ASAE 33:1637-1641. Peleg, M. 1988. An empirical model for description of moisture sorption curves. J. Food Sci. 53:1216-1219. Pilosof, A. M. R., Boquet, R., and Bartholomai, G. B. 1985. Kinetics of water uptake by food powders. J. Food Sci. 50:278-282. Rosentrater, K. A., and Muthukumarappan, K. 2006. Corn ethanol coproducts: Generation, properties, and future prospects. Int. Sugar J. 108:648-657. Sacchetti, G., Pittia, P., Biserni, M., Pinnavaia, G. G., and Rosa, M. D. 2003. Kinetic modelling of textural changes in ready-to-eat breakfast cereals during soaking in semi-skimmed milk. Int. J. Food Sci. Technol. 38:135-143. Shittu, T. A., Awonorin, S. O., and Raji, A. O. 2004. Evaluating some empirical models for predicting water absorption in African breadfruit (Treculia africana) seeds. Int. J. Food Prop. 7:585-602. Singh, B. P. N., and Kulshrestha, S. P. 1987. Kinetics of water sorption by soybean and pigeonpea grains. J. Food Sci. 52:1538-1541, 1544. Sopade, P. A., and Obekpa, J. A. 1990. Modeling water absorption in cowpeas, soyabean and peanut at three temperatures using Peleg’s equation. J. Food Sci. 55:1084-1087. Sopade, P. A., and Okonmah, G. N. 1993. The influence of physicochemical properties of foods on their water absorption behavior: A quantitative approach. ASEAN Food J. 8:32-36. Sopade, P. A., Ajisegiri, E. S., and Badau, M. H. 1992. The use of Peleg’s equation to model water absorption in some cereal grains during soaking. J. Food Eng. 1:269-283. Sopade, P. A., Ajisegiri, E. S., and Okonmah, G. N. 1994. Modeling water absorption characteristics of some Nigerian varieties of cowpea during soaking. Trop. Sci. 34:297-305. Spencer, H. M. 1926. Laboratory methods for maintaining constant humidity. Pages 67-68 in: International Critical Tables of Numerical Data, Physics, Chemistry and Technology. McGraw-Hill: New York. Tagawa, A., Muramatsu, Y., Nagasuna, T., Yano, A., Limoto, M., and Murata, S. 2003. Water absorption characteristics of wheat and barley during soaking. Trans. ASAE 46:361-366. Thakor, N. J., Sokhansanj, S., Patil, R. T., and Deshpande, S. D. 1995. Moisture sorption and volumetric changes of canola during hydrothermal processing. J. Food Process Eng. 18:233-242. Turhan, M., Seyar, S., and Gunasekaran, S. 2002. Application of Peleg model to study water absorption in chickpea during soaking. J. Food Eng. 53:153-159. Wang, N., Daun, J. K., and Malcolmson, L. J. 2003. Relationship between physicochemical and cooking properties, and effects of cooking on antinutrients, of yellow field peas (Pisum sativum). J. Sci. Food Agric. 83:1228-1237. Wood, J. A., and Harden, S. 2006. A method to estimate the hydration and swelling properties of chickpeas (Cicer arietinum L.). J. Food Sci. 71:E190-E195. Yu, L., Mazza, G., and Jayas, D. S. 1999. Moisture sorption characteristics of freeze dried, osmo-freeze-dried, and osmo-air-dried cherries and blueberries. Trans. ASAE 42:141-147.

[Received May 7, 2007. Accepted July 6, 2007.]

Vol. 84, No. 6, 2007

555