Abstract

Biltong is a dried meat product that is widely consumed in South Africa. The marinated meat is traditionally dried under ambient winter conditions while commercial biltong producers use hot air driers. Hot air drying is time-consuming and energy-intensive. A combined infrared and hot air drying (IRHAD) is an alternative method of drying meat during biltong processing. The aim of this study was to establish the effect of the infrared (IR) power, the temperature, and velocity of the drying air on the drying kinetics of marinated beef and subsequently select the best thin-layer drying model for IRHAD during biltong processing. Marinated beef samples were dried at IR power levels of 500, 750, and 1000 W; drying air temperatures of 30, 35, and 40°C; and air velocity of 1.5 and 2.5 m∙s-1. Results indicate that increasing the IR power and the drying air temperature increased the IR emitter temperature and the core temperature of the marinated beef sample. Consequently, increasing the drying rate thus reduced drying time. The air velocity had an inverse relationship with the IR emitter temperature, the core temperature of the marinated beef sample, and the drying rate. The drying process was characterised by a rising rate period in the first half an hour, followed by a falling rate period which implies that moisture transport occurred partly by surface evaporation and predominantly by diffusion. The effective moisture diffusivity ranged from to , while, the activation energy ranged between 40.97 and 59.16 kJ∙mol-1. The IRHAD of marinated beef during its processing to biltong was best described by the two-term model since it had the highest (0.9982-0.9993) and the lowest RMSE (0.0062-0.0099). The power level of the IR emitter of 1000 W combined with a drying air temperature and velocity of 40°C and 1.5 m∙s-1, respectively, showed the highest improvement in the drying kinetics and the lowest drying time of hours; hence, it is recommended as a possible drying alternative for the processing of biltong.

1. Introduction

Biltong is a dried meat product that is predominantly consumed in South Africa. Biltong is made by drying thin slices of meat that is marinated with a mixture of salt, vinegar and dry spices (coriander, black pepper, and brown sugar) [1]. The marinated slices of meat are traditionally dried under ambient conditions during the winter season (), whereas commercial biltong producers use temperature-controlled hot air dryers with the drying air temperature ranging between 20 and 40°C ([2, 3]; Burfoot et al., 2010; [4]).

The market value of biltong in South Africa was estimated at US$ 170 million in 2015 [5]. The market for biltong has expanded both locally and internationally prompting its increased production [6, 7]. The rising popularity of biltong has spurred the development of a variety of hot air biltong dryers. The hot air biltong dryers range from simple domestic types to high capacity commercial dryers.

The conventional hot air drying (HAD) method used in the commercial production of biltong is an energy-intensive drying method [8]. The low thermal conductivity of agricultural products combined with the case hardening of these products during HAD decelerates the moisture migration which results in longer drying time and increased energy consumption [9, 10]. HAD degrades the quality of agricultural products through the loss of colour, loss of heat-sensitive nutrients, deformation, and internal structure damage [1113]. Moreover, the temperature range of 20–40°C used in conventional HAD of meat during biltong production is not sufficient to achieve the recommended microbial reduction in the resultant biltong [3]. These concerns underscore the need for alternative drying methods for biltong production.

Alternative heating technologies like microwave, inductive heating, radiofrequency, and infrared (IR) provide volumetric heating that positively impacts on the energetic, exergetic, and heating efficiency [14]. According to Li et al. [15], IR radiation improves the dehydration efficiency of beef jerky by promoting protein denaturation which transforms immobilised water to free water, consequently reducing the activation energy. Cherono et al. [16] noted that IR drying of beef reduced the microbial count on the resulting biltong. However, moisture condensation observed on the surface of biltong during IR drying highlights the inability of natural convection to cope with the improved drying rates [17]. A combined infrared and hot air drying (IRHAD) could accelerate the removal of moisture from the meat surface to sustain the high drying rates. Drying agricultural products using IRHAD requires less energy and produces dried products of higher quality compared to using IR drying or HAD independently [18, 19].

The application of IRHAD as a possible alternative to HAD in the making of biltong requires a quantitative understanding of the heat and mass transfer in meat subjected to IRHAD. According to Trujillo et al. [20], Fickian diffusion is the predominant mode of moisture transport in meat during drying. The Page model and thin layer models derived from Fick’s Second Law of diffusion suitably describe the drying behaviour of agricultural products. Abe and Afzal [21] and Das et al. [22] identified the Page model as the best model for describing the thin-layer IR drying of rough rice. According to Toğrul [23], Puente-Díaz et al. [24], Sui et al. [25], and Sadin et al. [26], the Midilli model is the most suitable model for predicting the drying kinetics of apple, murta berries, wine grape pomace, and tomato slices, respectively, when subjected to IRHAD. The Logarithmic model was identified by Doymaz [27] as the best model for the IR drying of sweet potatoes. The study by Cherono [17] identified the approximation of diffusion model as the best model for the IR drying of marinated beef. Nonetheless, the two-term model and the Midilli model also showed good accuracy in predicting the IR drying of marinated beef.

A complete characterisation of the drying kinetics is critical in modelling the drying process [28]. The aforementioned studies highlight the research done in modelling the drying of a number of agricultural products subjected to IR drying and IRHAD. However, there is no literature on the drying characteristics of marinated beef subjected to IRHAD. A good understanding of the drying kinetics of marinated beef subjected to IRHAD would guide the application of IRHAD as a possible alternative to HAD in the processing of biltong. Therefore, the aim of this study was to establish the effect of the IR power, the temperature, and velocity of the drying air on the drying kinetics of marinated beef and subsequently select the best thin-layer drying model.

2. Materials and Methods

2.1. Sample Preparation

The samples were prepared as outlined by Muga et al. [1]. Beef portions procured from a local butchery (Pick n Pay, Pietermaritzburg, South Africa) were sliced along the muscle fibres to dimensions of [2]. The moisture content of each sample was determined by drying a 10 g portion of the sample at 105°C for 24 hours in a hot air oven (AX 60, Carbolite Gero Ltd., Hope Valley, UK) [29]. The experimental samples were put in sealed polythene bags and stored in a refrigerator (Defy C250, Defy Appliances (Pty) Ltd., Durban, South Africa) set at 4°C, awaiting marination.

The biltong marinade was made by mixing sodium chloride and vinegar in a ratio of 1 : 2, respectively, based on their mass. The experimental beef samples were marinated by adding 0.075 kg of marinade for every kg of beef [30, 31]. The marinade was uniformly spread over the beef samples to ensure that all sides of the beef samples were coated with the marinade. The marinated samples were refrigerated at 4°C, for 24 hours. The samples were turned every six hours during the refrigerated storage period to ensure uniform distribution of marinade [32].

2.2. The Drying Unit

The infrared-assisted hot air dryer was made by retrofitting an existing hot air cabinet dryer with an IR emitter. A schematic of the IRHAD experimental set-up is shown in Figure 1. The drying chamber is made of stainless steel and has a 100 mm thick polyurethane insulation. A drying platform () made of stainless steel mesh is suspended on a weighing balance (Shimadzu UW6200H, Shimadzu Corporation, Kyoto, Japan) that is placed on top of the drying cabinet. The weighing balance is connected to a computer via an RS-232C connector to log the mass of the sample being dried. A complete description of the hot air cabinet drier is provided in [1].

A ceramic IR emitter (T-FSR, Elstein-Werk M. Steinmetz GmbH & Co. KG, Northeim, Germany) is attached to the roof of the drying chamber, directly above the drying platform (Figure 1). The distance between the IR emitter and the drying platform is 280 mm. The IR emitter has a power rating of 1000 W and is of dimensions 250 by 62.5 mm. The wavelength of the IR radiation from the IR emitter ranges from 2 to 10 μm [33]. The power to the IR emitter was varied between 500 and 1000 W by regulating the current using an AC light dimmer circuit. The IR emitter has an inbuilt K-type thermocouple that was connected to a data logger (OM-DAQ-USB-2401, Omega, UK) to record its temperature.

The temperature and velocity of the drying air are monitored using a supervisory control and data acquisition (SCADA) system that is installed in the dryer. The core temperature of the sample placed on the drying platform is measured using a K-type thermocouple (Temperature Controls (Pty), Pinetown, South Africa). The K-type thermocouple is connected to a data logger (OM-DAQ-USB-2401, Omega, UK) to record the core temperature data.

2.3. Drying Experiments

A marinated beef sample was retrieved from refrigerated storage and allowed to equilibrate to room temperature for 60 minutes. Synchronously, the dryer was preheated for 60 minutes before the commencement of the drying experiments. A five-gram portion of meat was cut from the beef sample for the determination of the moisture content prior to drying [29]. Thereafter, a K-type thermocouple was inserted into the centre of the beef sample, and the sample was placed on the drying platform as shown in Figure 1. The IR emitter was instantaneously switched on, and the beef sample allowed to dry until 50% of its mass was lost. The core temperature and mass of the beef sample, as well as the temperature of the IR emitter, were recorded continuously throughout the drying experiment.

The drying experiment employed a three-factor, full-factorial experimental design. The first two factors had three levels, and the third factor had two levels. The factors studied were the IR power (500, 750, and 1000 W), the temperature of the drying air (30, 35, and 40°C), and the velocity of the drying air (1.5 and 2.5 m.s-1).

The data obtained from the experiments was subjected to the analysis of variance (ANOVA) at 5% level of significance. Where a significant ANOVA result was found, the mean comparison was done using Fisher’s Unprotected LSD method. All data analysis was done using GenStat® 18th Edition (VSN International Ltd., Hemel Hempstead, United Kingdom).

2.4. Evaluation of the Drying Characteristics

The moisture loss data obtained from the drying experiments was used to determine the instantaneous moisture content of the beef sample during drying. Thereafter, the moisture ratio (MR) was calculated using equation (1), and the drying rate () was obtained from the derivative of the MR with respect to drying time (equation (2)). where:

 = drying time (s)

 = instantaneous moisture content (kg of water/kg of dry matter)

 = the initial moisture content (kg of water/kg of dry matter)

The effective moisture diffusivity () was obtained from the linear form of the solution to Fick’s Second Law of diffusion as shown in equation (3) [34]. where:

 = the half-thickness of the beef sample (m)

The activation energy () was determined from the Arrhenius-type relationship between the drying temperature and the effective moisture diffusivity (equation (4)) [35]. where:

 = the preexponential factor of Arrhenius equation equivalent to diffusivity at the maximum temperature (m2∙s-1)

 = the activation energy (kJmol-1)

 = the universal gas constant (kJmol-1∙K-1)

 = the average beef sample temperature in (°K)

The activation energy () was obtained from the linear form of equation (4) (equation (5)) by determining the slope of the graph of against .

2.5. Selection of the Best Thin-Layer Drying Model

The thin-layer drying models considered in this study are outlined in Table 1. The models were selected based on previous research as highlighted in the introduction.

The experimental data was fitted to the selected thin layer dying models using the nonlinear least square analysis in MATLAB (MATLAB R18.2b, Mathworks, Inc., Natick, MA, USA). The best model was chosen based on a combination of the highest and the lowest RMSE.

3. Results and Discussion

3.1. IR Emitter Temperature

The temperature of the IR emitter recorded for the tested drying conditions is presented in Table 2. The temperature of the IR emitter was significantly () affected by the power level of the IR emitter. The highest and lowest IR emitter temperatures of and were observed at an IR power level of 1000 W and 500 W, respectively. Increasing the IR emitter power consumption increased the temperature of the IR emitter. The increase in the emitter temperatures is caused by the increase in ohmic losses as the power level of the IR emitter increases [9]. These results are consistent with previous findings by Ali et al. [39] and Ott et al. [40].

A change in the velocity of the drying air resulted in significant () variation in the temperature of the IR emitter. Higher IR emitter temperatures were observed at a low drying air velocity of 1.5 m∙s-1, while the drying air velocity of 2.5 m∙s-1 resulted in low IR emitter temperatures. The reduction in the temperature of the IR emitter at higher drying air velocity is attributed to the cooling effect induced by the drying air on the IR emitter due to convective heat losses. According to Das and Das [41], the radiation efficiency of mid to far IR emitters ranges between 40 and 60% with some heat lost via convection. Increasing the drying air velocity increases the mass flow rate of air which increases the convective heat losses, thus reducing the IR emitter temperature.

The temperature of the drying air also had a significant () effect on the temperature of the IR emitter. The temperature of the IR emitter increased with an increase in the drying air temperature. An increase in the drying air temperature reduced the temperature gradient between the IR emitter and the drying air, thus reducing the convective heat losses and vice versa. Consequently, the IR emitter temperatures increased with an increase in drying air temperature.

The power level of the IR emitter had a synergistic interaction with the temperature and velocity of the drying air that significantly () affected the temperature of the IR emitter. However, the two-way interaction between the drying air temperature and velocity, and the three-way interaction of all the experimental factors, had no significant effect on the IR emitter temperature.

3.2. Core Temperature of the Beef Sample during Drying

The average core temperature of the beef sample over the entire drying period ranged from to (Table 3). The core temperature of the beef sample was significantly () affected by the power level of the IR emitter. Higher core temperatures of the sample were observed at higher power levels of IR emitter. The increase in the core temperature of the beef sample with an increasing power level of the IR emitter is caused by the high temperatures of the IR emitter achieved at high power levels of the IR emitter.

The temperature of the IR emitter determines the quality and intensity of IR radiation [42]. According to Planck’s law, higher temperatures of the IR emitter shift the peak IR energy wavelengths towards the shorter wavelength [43]. The short IR wavelengths are associated with higher IR energy flux and vice versa [44]. The peak IR wavelengths calculated using the Wien’s displacement law and the average temperatures of the IR emitter were , , and , at 1000, 750, and 500 W, respectively. The shorter wavelength of observed at and an IR power level of 1000 W indicates an increase in the IR energy flux into the beef sample. The increased IR energy flux into the beef sample increases the core temperature of the beef sample at the IR power level of 1000 W compared to the IR power level of 750 and 500 W. These results concur with previous research that reported an increase in the core temperature of the sample with an increase in IR power [17, 45].

The velocity of the drying air caused significant () variation in the average core temperature of the beef sample. A higher core temperature of the beef samples was observed at a drying air velocity of 1.5 m∙s-1, whereas lower core temperatures were observed at a drying air velocity of 2.5 m∙s-1. An increase in the air velocity is associated with low IR emitter temperature which reduces the radiant energy flux into the sample, resulting in low core temperature of the samples. The inverse relationship between the drying air velocity and the core temperature of the beef sample can also be attributed to the evaporative cooling due to the increased mass flow rate of the drying air as observed by Kocabiyik and Tezer [46] during the IR drying of carrots.

The core temperature of the beef sample also varied significantly () with the temperature of the drying air. Increasing the temperature of the drying air increased the core temperature of the sample. The variation in core temperature of the sample became pronounced over time as shown by the divergence of the sample temperature curves in the later stages of drying (Figure 2). This variation can be attributed to the synergistic transfer of energy from both the IR and the hot drying air to the beef sample from the onset of drying until the sample temperature reaches the drying air temperature [19]. Thereafter, there is simultaneous energy transfer into and out of the beef sample. The energy influx is primarily due to the radiant flux from the IR emitter while the energy outflux is driven by the temperature gradient between the sample and the drying air. The temperature gradient between the beef sample and the drying air increases with a decrease in the drying air temperature, thus promoting a high energy outflux which results in lower core temperature of the beef sample.

3.3. Drying Time

The drying time is the duration spent in drying the beef sample until 50% of the original mass was lost. The drying time ranged from to hours (Table 4). This is significantly lower than the 24–230 hours required to dry biltong using the conventional hot air drying [3, 4, 17, 47]. The drying time achieved in this study is also lower than the range of drying time of 10.25–36 hours, achieved by Cherono [17] that solely used IR drying to make biltong. A similar reduction in drying times has been reported by Sharma et al. [48] and Li et al. [15] on IRHAD of onions and beef jerky, respectively. According to Hebbar et al. [19], the synergistic interaction between IR and hot air promotes rapid heating of the product being dried which accelerates the rate of mass transfer, thus improving the dehydration efficiency and subsequently lowering the drying time.

The variation in MR over the drying period is presented in Figure 3. The IR power significantly () affected the drying time. Shorter drying times were observed at 1000 W while the longest drying time was observed at an IR power level of 500 W. Similarly, the drying time decreased significantly () with an increase in temperature of the drying air. These observations are in agreement with the findings of Sharma et al. [48], Kumar et al. [49], and Nasiroglu and Kocabiyik [50]. An increase in both the IR power and the drying air temperature increases the core temperature of the beef sample. The increase in the core temperature of the beef sample indicates the increased heat influx which accelerates the vapourisation of water in the beef sample, thus increasing the vapour pressure and expediting the moisture transport out of the sample [51]. According to Jaturonglumlert and Kiatsiriroat [52], increasing the power of the IR emitter and the temperature of the drying air increases the ratio of heat and mass transfer coefficient, thus leading to the shortening of the drying time.

Contrary to the effect of the power level of the IR emitter and the temperature of the drying air, increasing the velocity of the drying air increased the drying time. Similar findings were reported by Afzal et al. [18] and Kocabiyik and Tezer [46]. This observation can be attributed to the negative effect of the drying air velocity on the core temperature of the beef sample. Increasing the drying air velocity decreases the ratio of heat and mass transfer coefficient, thus lengthening the drying time [52].

3.4. Drying Rate

The drying rate refers to the rate of moisture removal from the beef sample. The variation in the drying rate during the drying process is shown in Figure 4. The drying rate curves at the IR power level of 500 and 750 W indicated a short rising rate drying period followed by the first and second falling rate drying periods which is typical of drying rate curves obtained from IRHAD [53]. The drying rate curve obtained at the IR power level of 1000 W showed a short constant rate drying period immediately after the rising rate period and just before the first falling rate period. This observation is consistent with the results obtained by Mongpraneet et al. [54] when drying welsh onions under IR radiation.

IR radiation has an efficient heat transfer mechanism that promotes a rapid increase in the core temperature of the beef sample at the onset of drying, thus increasing the drying rate [43]. The rising rate drying period occurred within the first 30 minutes of drying and coincided with an increase in the rate of change in the core temperature of the beef sample (Figure 4). The moisture transport during this period can be attributed to surface evaporation [53].

The rising rate drying period is superseded by the first falling rate drying period. The first falling rate drying period coincides with a decrease in the rate of change in the core temperature of the beef sample (Figure 4). This drying period occurred between 30 minutes and 6 hours of drying at IR power levels of 500 and 750 W, while the same occurred at between 30 minutes and 2 hours of drying at IR power level of 1000 W. According to Trujillo et al. [20], diffusion is the predominant mode of moisture transport during the first falling rate drying period. The decrease in drying rate can be attributed to the reduced heat flux into the sample which lowers the vapour pressure gradient and subsequently reduces the moisture transport out of the sample.

The second falling rate period occurred after the first 6 hours of drying at IR power levels of 500 and 750 W, while the same occurred after the first 2 hours of drying at an IR power level of 1000 W. The second falling rate period is characterised by a low and fluctuating drying rate. According to Kocabiyik [53], most of the free water is lost during the first falling rate period leaving the bound water to play an active role in the second falling rate period. Bound water requires more energy to extract from the food matrix. However, there is minimal heat flux into the beef sample during this drying period as indicated by the relatively constant product temperature. The low drying rate may also be attributed to the additional internal resistance to moisture movement as the sample moisture content decreases [55].

The drying rate was significantly () affected by the power level of the IR emitter, the temperature, and the velocity of the drying air. This is attributed to the effect that these factors had on the core temperature of the beef sample, consequently affecting the amount of heat flux into the sample. Feyissa et al. [56] reported that heat transfer is the key driver of mass transfer during the roasting of chicken meat. The coupled relationship between the core temperature of the beef sample and the moisture transport is illustrated by the synchrony between the drying rate curve and the curve depicting the rate of increase in sample temperature (Figure 4).

3.5. Effective Moisture Diffusivity

The effective moisture diffusivity () of marinated beef subjected to infrared assisted hot air drying ranged between and (Table 5). These values are within the range of values obtained by Li et al. [15] for beef jerky dried under mid and far IR radiation. The values obtained in this study are greater than those reported in previous research on hot air drying of meat product such kaddid [57], eland jerky [58], chicken breast meat [59], and beef [60]. According to Li et al. [15], the selective heating of water molecules by the IR radiation causes the immobilized water in the myofibrillar network to migrate out of the network. The migration of the immobilized water out of the myofibrillar network increases the internal free water content, consequently expediting the transport of the free water to the surface by diffusion.

The IR power level caused significant () variations in under all drying conditions. The temperature of the drying air only caused significant variations between 30 and 40 °C, whereas the drying air velocity had a significant effect on the at IR emitter power levels of 1000 and 750 W but not at 500 W. Khir et al. [61] reported that increases with an increase in sample temperature during IR drying. The effect of the IR power, the temperature, and the velocity of the drying air on the in this study is attributed to their influence on the core temperature of the beef sample.

3.6. Activation Energy

The activation energy (Ea) ranged between 40.97 and 59.16 kJ∙mol-1 (Figure 5). The Ea obtained in this study is comparable to the Ea of 32.8 kJ∙mol-1 obtained by Li et al. [15] during mid- and far-IR drying of beef. The difference between the Ea obtained in this study and that reported by Li et al. [15] may be due to the difference in the range of IR wavelength used in both studies. The wavelength of the IR radiation used by Li et al. [15] ranged between 2.9–3.1 and 5.8–6.2 μm, whereas the wavelength of the IR radiation used in study ranged between and . The IR radiation wavelengths used by Li et al. [15] can effectively heat up the water molecules due to their proximity to the maximum absorption wavelengths of water of 3 μm and 6 μm [9]. Consequently, the availability of free water in the food sample increases. Increased availability of free water in the food matrix is associated with reduced Ea [62].

The inverse relationship between the Ea and the power level of the IR emitter observed in this study is attributed to the change in the wavelength of the IR radiation. Increasing the power of the IR emitter shortens the wavelength of the IR radiation from to . The absorption of the IR radiation by the water molecules in the beef sample increases as the radiation wavelength decreases towards the maximum absorption wavelength of water of 6 μm, thus lowering the activation energy.

3.7. Selection of the Best Thin-Layer Drying Model

The average coefficient of determination and the average root mean square error for the five thin-layer drying models considered in this study are shown in Table 6. All the five models had an and a . The high and the low RMSE indicate that all the five models can predict the changes in MR of marinated beef subjected to IRHAD with acceptable accuracy [63]. Overall, the two-term model had the highest (0.9982-0.9993) and the lowest RMSE (0.0062-0.0099) across all the IRHAD conditions. Therefore, the two-term model is the most suitable model in predicting the drying behaviour of marinated beef subjected to IRHAD.

A summary of the two-term model coefficients and the statistical parameters for all the experimental drying conditions is presented in Table 7.

The plot of the experimental MR and the MR predicted using the two-term model is shown in Figure 6. The data points lie along the 45° line which confirms the high degree of agreement between the experimental and predicted MR values.

4. Conclusion

The IRHAD of marinated beef is significantly affected by the power level of the IR emitter, the temperature, and the velocity of the drying air. The power level of the IR emitter is the predominant factor that influences the drying process of marinated beef subjected to IRHAD. The power level of the IR emitter determines the temperature of the emitter, thus influencing the core temperature of the marinated beef sample. The drying air temperature and drying air velocity have a direct and inverse relationship, respectively, with IR emitter temperature and the core temperature of the beef samples. The core temperature of the beef sample is the most critical parameter that influences the drying kinetics of marinated beef subjected to IRHAD. The combined effect of the power level of the IR emitter, temperature, and velocity of the drying air on the drying time, drying rate, effective moisture diffusivity, and the activation energy is directly related to how these factors influence the core temperature of the beef sample. The drying of marinated beef at IR power levels of 500 W and 750 W has a short rising rate drying period followed by the falling rate drying period. The drying process at 1000 W has a short constant rate period sandwiched between the rising rate and the first falling rate period. The rising and falling rate drying periods imply that the moisture transport during IRHAD of marinated occurs partly by surface evaporation and predominantly by diffusion. The two-term thin-layer drying model best describes the drying of marinated beef under IRHAD. The effective moisture diffusivity ranges from to , while the activation energy ranges between 40.97 and 59.16 kJ∙mol-1. The results from this study can inform the application of IRHAD in biltong processing. A power level of 1000 W with an IR radiation wavelength of resulted in the shortest drying time, the highest effective moisture diffusivity, and the lowest activation energy. The power level of the IR emitter of 1000 W combined with a temperature and velocity of the drying air of 40°C and 1.5 m∙s-1, respectively, highly improved the drying kinetics of biltong, hence is recommended as a possible drying alternative for biltong processing.

Data Availability

The data is available on request through the corresponding author.

Conflicts of Interest

The authors declare that there is no conflict of interest regarding the publication of this paper.

Acknowledgments

The authors appreciate the financial support from the University of KwaZulu-Natal towards this research.