Next Article in Journal
Spoilage Potential of Contaminating Yeast Species Kluyveromyces marxianus, Pichia kudriavzevii and Torulaspora delbrueckii during Cold Storage of Skyr
Next Article in Special Issue
Rheological Properties, Particle Size Distribution and Physical Stability of Novel Refined Pumpkin Seed Oil Creams with Oleogel and Lucuma Powder
Previous Article in Journal
Microwave Irradiation: Effects on the Change of Colour Characteristics and Main Phenolic Compounds of Cabernet Gernischt Dry Red Wine during Storage
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Extrusion Simulation for the Design of Cereal and Legume Foods

by
Magdalena Kristiawan
1,*,
Guy Della Valle
1 and
Françoise Berzin
2
1
INRAE, UR 1268 Biopolymers Interactions and Assemblies (BIA), 44316 Nantes, France
2
Université de Reims Champagne Ardenne, INRAE, FARE, UMR A 614, 51097 Reims, France
*
Author to whom correspondence should be addressed.
Foods 2022, 11(12), 1780; https://doi.org/10.3390/foods11121780
Submission received: 11 May 2022 / Revised: 8 June 2022 / Accepted: 12 June 2022 / Published: 16 June 2022
(This article belongs to the Special Issue Food Rheology and Applications in Food Product Design)

Abstract

:
A 1D global twin-screw extrusion model, implemented in numerical software, Ludovic®, was applied to predict extrusion variables and, therefore, to design various starchy products with targeted structure and properties. An experimental database was built with seven starchy food formulations for manufacturing dense and expanded foods made from starches, starch blends, breakfast cereals, pulse crop ingredients such as pea flour, fava bean flour, and fava bean starch concentrated, and wheat flour enriched with wheat bran. This database includes the thermal and physical properties of the formulations at solid and molten states, melt viscosity model, extruder configurations and operating parameters, and extruded foods properties. Using extrusion and viscosity models, melt temperature (T) and specific mechanical energy (SME) were satisfactorily predicted. A sensitivity analysis of variables at die exit was performed on formulation, extruder configuration, and operating parameters, generating the extruder operating charts. Results allowed the establishment of relationships between predicted variables (T, SME, melt viscosity) and product features such as starch and protein structural change, density and cellular structure, and functional properties. The extrusion operating conditions leading to targeted food properties can be assessed from these relationships and also the relationship between extrusion operating parameters and variables provided by simulation.

1. Introduction

Many starch-based foods, such as texturized ingredients, breakfast cereals, snacks, crackers, pasta, noodles, pet food, and many others, are produced by extrusion processes [1]. Starches and their blends are extruded to gain desirable end-usage properties. Ready-to-eat breakfast cereals are still the market leader in extruded foods. Society is leaning increasingly towards more sustainable and healthy diets by transitioning from animal to plant protein sources and enriching formulations with legume proteins and dietary fibers. Despite their nutritional advantages (high protein content, low glycaemic index) compared to cereal-based extruded foods, legume-based extruded foods are still rare. The growing demand for extruded snacks (4.4% per year in the market) makes the development of healthy extruded snacks even more desirable [2].
Under mechanical stresses and heat developed between screw and barrel during the extrusion process, the biopolymers undergo structural modifications: disruption of starch granules, melting of starch crystals, depolymerization of starch macromolecules, denaturation and aggregation of proteins, depolymerization of dietary fibers, and product browning due to Maillard reactions [3,4,5]. These changes are very important because they provide extruded starchy foods with an attractive appearance, structure, and relevant functional properties.
Due to an abrupt pressure drop at the die exit, the flash evaporation of superheated moisture forms a porous and expanded melt structure that is set after crossing its glass transition temperature. The structural features of the expanded foods govern their texture, functional and nutritional properties, and consumer acceptability [6,7,8,9,10,11,12]. The structural features consist of density, cellular structure, and properties of constitutive material, i.e., cell walls [13]. A phenomenological model of expansion allowed us to link extruded food features, such as expansion indices and cellular fineness to melt shear viscosity [14,15]. The melt viscosity takes into account the effect of formulation, such as fiber, protein, and starch contents, amylose to amylopectin ratio, and biopolymer structural changes during processing [16].
The development of mechanistic models of extrusion, based on the continuum mechanics approach, provides many benefits: a better understanding of the extrusion process and the development of computer-aided tools in simulation, optimization, and scale-up. In twin-screw extrusion modeling, two main approaches can be considered. In the first approach, local flow can be described with accuracy in a limited portion of the extruder by sophisticated numerical 2D or 3D models based on the finite element method (FEM) [17,18,19,20,21,22]. It can provide a very accurate description of the flow field only in a specific part of the extruder, and it is expensive in terms of computational resources and time.
The second approach is a global one; it is based on a simplified flow assumption and provides average values of the local flow variables (such as temperature, total dissipated energy, pressure, viscosity, filling ratio, and mean residence time) of the whole process, along the screws, from the hopper to the die exit [23,24]. It is based on various assumptions: (1) the flow is locally 1D or 2D; (2) the melting occurs instantaneously just before the first restrictive screw element; (3) the molten material behaves locally as a Newtonian fluid [24]. Computation of the flow parameters is done separately for each type of screw element and for the die components. The elementary models are linked together to obtain a global description of the flow along the extruder. This 1D global model has been implemented into a software called Ludovic® v7.0.0 Classic Edition (Sciences Computers Consultants, Saint Etienne, France). It is quite straightforward to provide predictions of flow variables with simple computing resources and time: about a few seconds per simulation using a laptop configured with Intel® Core™ i7-7600U processor (dual core). A fairly good agreement was obtained between simulation results from this 1D global model and a full 3D FEM model [18].
Despite significant progress in modeling and simulation, the design of extruded products at the industrial level is still based on a trial-and-error approach. One of the main challenges is determining the viscous behavior of melts under extrusion-like conditions that require specific rheometers, e.g., an in-line slit die rheometer or pre-shearing capillary rheometer (Rheoplast®) (Société Courbon, Saint-Etienne, France) [25,26]. This difficulty can be circumvented, in the first approach, by assuming that starch viscosity is a function of amylose content, and it can be interpolated from the published data on maize starches having several amylose contents [27].
The 1D global model has been widely used for the design of experiments of processing, but it is still limited to simple formulation (maize and wheat starches and wheat flour) [23,28,29,30]. In process design, the model is exploited as a numerical screening tool of operating conditions to define and optimize the experimental domain in terms of flow variables (melt temperature, specific mechanical energy (SME), pressure). Its implementation in product design, for a wide range of complex formulations encompassing various structural and functional properties for one product, has never been tackled before. Until now, a simple link between a single property of extruded starchy products (for example, hydrosolubility) and a single process variable (for example, SME) was considered [29]. In this context, the aim of this paper is to test whether the 1D global extrusion model, implemented with an appropriate viscosity model, can be deployed as a prototype of a computer-aided tool for predicting properties of extruded foods and structural changes of biopolymers based on realistic formulations, from operating conditions. In other words, it will serve for product design applications.
For this purpose, we selected seven starchy food formulations to manufacture dense and expanded foods: cereal and tuberous-based pure starches and their blends with a wide range of amylose–amylopectin ratio, legumes-based ingredients (pea flour, fava bean flour, and fava bean starch concentrate), breakfast cereals, and blends of wheat flour–wheat bran. A detailed experimental database was built, including input and output data for extrusion simulation and extruded product structural and functional features. After describing the 1D extrusion model, virtual extrusion trials were performed, and validation steps were systematically detailed. The extruder operating curves representing a variety of predicted extrusion variables in the function of operating parameters were built. Then we strived to show how products with targeted structural and functional properties can be manufactured using the relations between extruder operating charts, predicted extrusion variables, and extruded product features, including biopolymer structural changes.

2. Materials and Methods

2.1. Building Database

Database for twin-screw extrusion simulation has been collected from literature and our own experimental results [31,32,33,34,35,36]. Due to the large number of experimental works, case studies were selected by applying several criteria: (1) presence of melt viscosity model and melting transition data as input, (2) existence of experimental data of melt temperature and specific mechanical energy (SME), and (3) properties of final products exhibiting clear relations with extrusion variables.

2.1.1. Raw Materials

Seven food recipes have been selected, reflecting a wide range of complexity of food formulations. The recipes were extruded to manufacture foods from the simplest ones, e.g., dense potato starch [32] and expanded starches (high amylose maize, waxy maize, and their blends revealing a broad variation of amylose content [27,33], and potato [31], to complex ones, e.g., breakfast cereals, and expanded snacks made from (1) legumes-based ingredients (pea flour [34], fava bean flour, and fava bean starch concentrate), and (2) blends of wheat flour and wheat bran [35,36,37] (Table 1). The extrusion data set of fava bean ingredients and breakfast cereals were produced from our original works. For these cases, raw material chemical composition was determined according to standard methods [34].
Two breakfast cereal foods were obtained by extrusion of a cacao-based formulation (called CB/Choco Ball) composed of a blend of rice flour, wheat flour, cacao powder, sugar, and salt and of a honey-based formulation (HB/Honey Ball) composed of a blend of corn flour, malt extract, sugar, and salt.
An overview of the chemical composition of all extrusion feeds is given in Table 1.

2.1.2. Extrusion Trials

The extrusion variables cover a broad spectrum of moisture content (15–36% wet basis), melt temperature (85–190 °C), and specific mechanical energy (100–5500 kJ/kg). An overview of extrusion systems for all feeds is given in Table 2, including equipment specificities.
Maize and potato-based expanded starches (Maize A, Maize B, Maize C, Maize D, and Potato 1; case n°1) and wheat-based snacks (WF, WF-LB, WF-HB; case n°7) were manufactured using pre-pilot extruders having a screw diameter between 25 and 55.5 mm [27,31,33,35,36]. The screw profile included conveying screw elements followed by a left-handed screw element. While maize starch blends were extruded through a long-slit die, the extrusion of potato starch and wheat flour-wheat bran used a short circular die. Pea flour extrusion (PF; case n°5) [34] used a more complex screw profile: conveying elements followed by three kneading discs staggered at 45°, and then three left-handed elements. Several die geometries were used, circular (D = 3 mm) of various lengths (5, 10, and 15 mm), attached or not with a multi-step slit die rheometer, as described in detail in Appendix A (Figure A1).
Dense pea and potato starch (PS, Potato 2; case n°3) [32] and fava bean-based snacks (FF, FSC; case n°6) were extruded using laboratory-scale extruders with a screw diameter of 11 and 16 mm and equipped with a cylindrical die. The screw profile included conveying elements, followed by kneading discs staggered at 90° and one left-handed screw element. For manufacturing dense (non-expanded) products (PS, Potato 2), the temperatures of the last barrel (Tb) and die (Td) were regulated at 80–95 °C, otherwise, the temperatures were set at melting temperature level.
Extrusion of breakfast cereal formulations (HB and CB; case n°4) was carried out on a pre-pilot co-rotating twin-screw extruder (Coperion Werner & Pfleiderer ZSK 26Mc) equipped with a circular die. The same extruder and screw configuration, and temperature profile of barrel and die were used, as for pea snacks manufacturing [34]. The die had a diameter of 3 mm and a length of 15 mm.

2.2. Simulation of Extrusion Trials Using Ludovic® Software

2.2.1. Extrusion Modelling

The extrusion simulation is based on a global model developed by Vergnes et al. [24]. This model uses a 1-D approach under the assumption that molten material behaves locally as a Newtonian fluid.
The flow is described by the Stokes equations in cylindrical coordinates (r, θ , z). Each conveying element is considered as a C-chamber associated with an angular sector representing the intermeshing area. Assuming that temperature and viscosity are locally constant, the velocity field can be determined in each element, and the volumetric flow rate Q v through the screw channel is expressed as:
Q v = A Ω + B 1 η Δ P Δ θ
where A and B are coefficients depending on screw geometry, Ω is the screw rotation speed (rad·s−1), Δ P is the pressure variation (in Pa) along angular screw section Δ θ . The first term on the right side stands for drag flow, and the second one for pressure-driven flow. The non-Newtonian behavior of the molten material is taken into account by applying an appropriate model of shear viscosity η ( γ ˙ ) into Equation (1), where γ ˙ is the local shear rate. Specific models have also been developed for reverse screw elements [38] and blocks of kneading discs [39].
Flow equations are solved for each section of the screw and the die. These solutions are linked together to obtain a global description of the molten material flow along the extruder. The temperature changes Δ T in a section are obtained from the energy balance of accumulation, viscous dissipation E v , and heat transfer by conduction to the barrel E cd :
Cp   Δ T + E cd = E v = W ˙ ρ m Q v
with   W ˙ = V η γ ˙ 2 dV
where W ˙ is the power dissipated in the volume V of a screw element, ρ is the material melt density, and C p is its specific heat.
The melting zone is determined according to experimental observations: it is located at the first restrictive element (left-handed element, kneading disc), and melting is assumed to be instantaneous [40]. It is supposed that, at the point where the pressure starts to increase, upstream of the first restrictive element, the material is in the molten state and its temperature is equal to the melting temperature (Tm).
The computation starts from the die exit and proceeds backwards. As the final product temperature is unknown, an iterative procedure is used. A final temperature of the product at the die exit is chosen, and the pressure drop inside the die is calculated to go back to the upstream pressure value (to the upstream screw elements). The evolutions of pressure and temperature are thus determined element by element (section by section), from downstream to upstream until the melting point. The temperature calculated at this point is then compared with the melting temperature: if they are equal, the calculation is completed; otherwise, the final temperature is modified, and the calculation is repeated until convergence.

2.2.2. Input Variables

The input data are divided into three groups. The first group refers to extruder configuration. It covers screws, barrels, die, and also the position of specific zones (feeding, degassing). Examples of the screw profiles, for all feeds composed of right-handed and left-handed screw elements, and blocks of kneading discs are given in Appendix A (Figure A2). Dies are characterized by elementary elements (e.g., slit, pipe, associated in series or parallel).
The second group refers to food properties:
(1)
Physical and thermal properties
They cover density ρ , heat capacity C p , and thermal conductivity k in the solid and molten state, and melting temperature and enthalpy. Their dependence on moisture content (wet basis) can be calculated according to additive rules from the values for water and dry molten starches (Appendix A, Table A1) [41,42,43]. Most results on the melting temperature (Tm) as a function of moisture content of studied materials can be retrieved from the literature [32,34,42,44,45,46,47] (Appendix A, Figure A3). The melting point of maize starches, having different amylose content (5%, 23%, 47%, and 70%), were obtained from the data set of waxy maize (amylose content of 5%); wheat (25%), smooth pea (35%), and wrinkled pea (75%) starches. The melting temperature dataset of breakfast cereal formulations (HB, CB), and fava bean flour and starch concentrate were determined experimentally as the peak of the DSC melting endotherm using standard procedure [32,48] (Appendix A, Figure A3).
(2)
Rheological properties of molten starchy foods
The variation of starchy melt viscosity ( η ) with shear rate ( γ ˙ ) is generally described by the Ostwald–de Waele power law [25]:
η = K γ ˙ n 1
where n is the flow index (0 ≤ n ≤ 1) and K (Pa·sn) the consistency.
The parameters K and n were expressed as a function of temperature (T, in °C), moisture content (MC), and SME according to [27]:
K = K 0 exp ( E R ( 1 T a 1 T 0 ) α ( MC MC 0 ) β ( SME SME 0 ) )
n = n 0 + α 1   T + α 2   MC + α 3   SME + α 4   T   MC + α 5   T   SME + α 6   MC   SME
where K 0 is the consistency at reference conditions, E is the activation energy, R is the gas constant, T a is the absolute temperature (in K), α is the water plasticization coefficient, β reflects the thermomechanical history coefficient, and n 0 and α i are constants accounting for the dependences of the flow index on T, MC, and SME and their interactions. The other constants, MC 0 , T 0 , and SME 0 , stand for reference moisture content (10%), temperature (353 K), and SME (350 kJ/kg).
The viscosity model parameters of starchy melts, e.g., maize starches having different amylose content, potato starch, breakfast cereals, and blend of wheat flour–wheat bran were retrieved from the literature and are presented in Table 3. Under the hypothesis that viscous behavior is governed only by the amylose content of the starch, the viscosity of fava bean starch concentrate was approached by maize starch having an amylose content of 47%. The viscosity behavior of pea starch and fava bean flour was approached by the pea flour. Indeed, a previous study revealed that the master curve of shear viscosity of these feeds, determined by a pre-shearing capillary rheometer (Rheoplast®) (Société Courbon, Saint-Etienne, France), were very close to each other [49].
Finally, the third group of data concerns the operating conditions of the extruder. This data covers the barrel and die temperature profile, the moisture content, the total flow rate, the screw speed, and the heat transfer coefficients between molten material and barrel, screws and die, expressed through a Nusselt number. The data was overviewed in Table 2 for all recipes. For all simulations, regardless of the recipes, the Nusselt number was set to 20.

2.2.3. Output Variables

The main variables predicted by the 1D global extrusion model (Ludovic® v7.0.0 Classic Edition) (Sciences Computers Consultants, Saint Etienne, France) are melt temperature T, SME, pressure, and viscosity. These variables can be linked to the biopolymers transformation (e.g., starch melting and depolymerization, and protein aggregation), structural (e.g., expansion) and functional properties (e.g., water solubility index of starch) of extruded foods, considered as indirect output variables. Starch structural changes were represented by some features, such as the number of remnant starch granules/unit area, residual molar mass of starch, and intrinsic viscosity. Protein cross-linking by covalent bonds other than disulphide bonds was quantified by the amount of in-extractable proteins in sodium dodecyl sulphate (SDS) and dithioerythritol (DTE). Expansion behavior was characterized by density, sectional expansion index, and cell density (number of cells per volume unit). These properties have been listed in Table 4, hence completing the database, with their origin (literature, this work) and measurement methods. Actually, when these features were not available in the literature, they were characterized especially in this present work. In this case, methods have been detailed in Appendix A.2.

3. Results and Discussion

3.1. Validation of Extrusion Simulation

The reliability of the extrusion simulation was confirmed: (1) by comparing computed extrusion variables with the measured ones (T and SME) (Figure 1), (2) by studying the realistic feature of the axial profile of computed flow variables along the screw (Figure 2 and Figure 3), and (3) by examining in detail the responses of computed flow variables, upon the change of operating conditions (Figure 3 and Figure 4).

3.1.1. Comparison between Experimental and Simulation Results

Like in the experimental trials, the predicted melt temperature T_com was retrieved at the die exit. For most simulations, the computed temperature is overestimated by 5–10% compared to the measured temperature. The uncertainty of measurement can explain this difference. Temperature sensors are flush-mounted on the die wall with little penetration into the melt; consequently, the measured temperature differs from the actual one at the core of the product. Nevertheless, the computed temperature correlates well with the measured temperature, regardless of the formulation and extrusion system (R2 = 0.88) (Figure 1a). Good correlations were also obtained between predicted and measured SME values (R2 = 0.71–0.96), depending on the formulation (Figure 1b). The predicted SME is underestimated by 35–65% because the extrusion model only considers the melting energy and the viscous dissipation energy after the melting section, and neglects the solid transport and solid friction energy (inter-particle and particle/metal, metal/metal), which depend on raw material. All these trends in predictions: underestimation of SME and overestimation of temperature are in agreement with published results on extrusion simulation of starches and wheat flour [28,30]. The fair predictions of T and SME confirm the pertinence of the extrusion model for the design of extruded starchy products.

3.1.2. Axial Profile of Predicted Flow Variables

Simulation results were illustrated for two recipes of breakfast cereals taken as examples by presenting the axial profiles of the predicted variables (melt temperature, specific mechanical energy (SME_com), pressure, and viscosity) along the screws (Figure 2a) for typical moisture, feed rate and screw speed conditions (Figure 2 and Figure 3). These profiles give us insight into thermomechanical conditions along the screws. The results obtained for other recipes exhibit similar trends and are not presented here.
The temperature increased progressively to values superior to the melting temperature Tm (Tm = 133 °C for HB and Tm = 165 °C for CB) along the kneading discs and rose sharply along the left-handed elements (Figure 2b). The melt flow across the kneading discs caused a temperature increase of about 5 to 10 °C, depending on the material and flow conditions. A more severe temperature step was observed for left-handed elements: 35 to 60 °C for the three elements. The increase in temperature was essentially due to intense shear stresses across the restrictive elements causing high viscous dissipation. The temperature reached a maximum at the outlet of the last left-handed element. Then it decreased slightly in the melt conveying section before a final slight increase through the die.
The cumulative SME_com was calculated by summing values in each screw element. The melt flow across the restrictive elements (kneading discs and left-handed screw elements) increased the SME_com (Figure 2c), particularly in the left-handed elements. The sudden increase in SME_com was concomitant with the step increase in temperature, underlining the importance of viscous dissipation, which was mainly governed by the restrictive screw elements.
The axial pressure profile shows that the main part of the extruder was only partially filled since the relative pressure was zero for screw length between 50 and 200 mm (measured from the die) and beyond 300 mm. The left-handed screw elements and the right-handed screw element in front of the die were the only zones where the melt was under pressure (Figure 2d).
Since the viscosity depends on temperature and SME, the melt viscosity profile along the screws is also related to temperature and SME profiles (Figure 2b,c,e). At the melting section (in the first kneading bloc), the product temperature and SME were relatively low, and consequently, the local viscosity was high. In the following restrictive elements, and the conveying section nearby die inlet, there was an increase in temperature and SME due to the increase in viscous dissipation, leading to a decrease in viscosity. The knowledge of melt viscosity, mainly at the die exit, is very important, because it can be hardly measured experimentally and it partly controls the melt expansion at die exit [14].
The singularities of variable profiles (pressure peak, changes of temperature slope and energy increase and viscosity decrease) meet exactly the position of restrictive elements. This fact underlines the importance of viscous dissipation, mainly governed by the location of restrictive screw elements. Similar trends were observed for extrusion simulation of starches [28,29]. The CB recipe, once fully molten, exhibits higher temperature, SME, viscosity, and pressure along the whole screw length. Brugger et al. [50] observed that at close T and MC, the CB exhibited higher viscosity than the HB.
According to Figure 2b, the maximum temperature is predicted at the last restrictive element, i.e., at the position where it is difficult to measure it. For CB, at the selected operating condition (MC = 15%, N = 500 rpm, Q = 20 kg/h), the maximum temperature was high (230 °C) which suggests that food degradation may occur. For HB, under the same conditions, the maximum temperature was lower, but still very high for food stuffs (185 °C).
The importance of viscous dissipation on thermo-mechanical history was also proven by the negative effect of moisture content on the axial profile of all major flow variables (Figure 3a–c). When the moisture content is increased from 19% to 23%, the magnitude of temperature and SME_com is reduced all along the screws by a factor of 1.5 and 3, and pressure underwent a larger drop at the die entrance, by a factor of 4.5. The viscosity dependence on moisture content can explain this. Higher moisture content leads to a lower melt viscosity, which decreases energy dissipation and then lowers temperature and pressure.
Another important variable for a twin-screw extruder is the filling ratio, which is proportional to Q/N (Figure 3d). For each moisture content, the predicted SME_com of the CB recipe is negatively correlated with Q/N, following a power function. The same trend describing the effects of Q/N on SME was observed for experimental results of twin-screw extrusion of starchy products [28,32,34]. In fact, SME is usually proportional to N/Q: an increase of screw speed (N) or a decrease in feed rate (Q) will induce torque reduction, due to the shear-thinning behavior of melt, as well as the lower filling ratio of the extruder and more important viscous heating. According to Figure 3d, moisture content also affects SME-filling ratio curves negatively, which is explained by the influence of the water content on the viscosity.

3.1.3. Sensitivity Analysis

The impact of changing one operating conditions of the three ones (screw speed N, total flow rate Q, last barrel temperature Tb) on extrusion variables was examined systematically, while keeping other conditions (i.e., extruder configuration, screw profile) constant. Results are presented as 2D-operating charts and illustrated in Figure 4a–d for the manufacturing of breakfast cereal (HB recipe) at a moisture content of 19%.
As indicated previously, SME_com increases with screw speed and decreases with flow rate (Figure 4a). For each Q, SME increases by 20–30% when N increases from 150 to 700 rpm. At any N, SME drops by 35–40% when Q increases from 10 to 50 kg/h, with the highest decline occurring at the beginning of Q change (from 10 to 20 kg/h). The variations of T at die exit (T_com), when N and Q are changed, follow the same trend as SME_com, but it is quite insensitive to the change of Q (Figure 4b). The inverse trend with N and Q was observed for computed die entrance pressure (P_com): it decreases with N and increases with Q by a factor of 1.75–2.5 (Figure 4c). The negative impact of Q was also observed on predicted melt viscosity at the die exit: η_com drops by a factor of 2–2.25 (Figure 4e). However, an increase of Tb from 90 to 140 °C leads to a decline in η_com by a factor of 1.2–1.35. The increase of screw speed and the decrease of flow rate positively affects the viscous dissipation and, therefore, the SME and melt temperature. Indeed, the melt temperature increase leads to a decrease in the viscosity, and consequently, pressure.
Finally, the effect of conductive heating was studied on melt temperature. Regardless of N, an increase in barrel temperature (Tb), from 90 to 140 °C, leads to melt temperature (T_com) rising, by 5–10%, following linear trends (Figure 4e). However, the melt temperature is always higher than that of the barrel. It proves that viscous dissipation has more impact than the heat transfer between barrel and melt.
The set of Figure 4b–d can be used to evaluate the upper limit of operating pressure and melt temperature. For example, let us fix the upper limit of maximum melt temperature and pressure to be 165 °C and 10 MPa. The feasible operating domain for barrel temperature (Tb) of 90 °C at the highest flow rate (Q = 50 kg/h) is then screw speed in the interval of 300–550 rpm.
The order of magnitude of predicted extrusion variables and the trend of their responses to changes of extruder operating conditions and screw elements are similar to those observed for the extrusion of standard cereal and starches [28,29,56]. The similarity in the trend of responses was found for all studied starch-based recipes, likely because the responses are mainly governed by the viscous (shear thinning) behavior of the molten starch. These facts validate the 1D global extrusion model for the design of starchy foods. In addition, the 1D global model of twin-screw extrusion allows us to better understand the machine (extruder) and extrusion process itself. The knowledge of temperature profile and the possibility of adjustment of processing conditions using simulation appears thus advantageous, particularly for thermo-sensitive foods.

3.2. Application in Food Design

In this section, extrusion modeling is applied to determine the operating conditions that lead to a product with desired structure and properties. This approach is based on: (1) relationships between computed extrusion variables and product features (Figure 5a, Figure 6a,b, Figure 7a,b and Figure 8a–c, and Appendix A, Table A2 for equations) and (2) using these relations and the global extrusion model to select an interval of operating conditions that permit to obtain the desired properties

3.2.1. Intrinsic Viscosity of Expanded Starches

Starch transformation was assessed by the intrinsic viscosity, reflecting the extruded product’s molar mass. It was normalized to the property of native starch, resulting in a relative drop of intrinsic viscosity ΔVis (%). The higher ΔVis, the larger the extent of degradation, the more depolymerized extruded starch, and the lower the molar mass is. Extrusion conditions for expanded maize and potato starches (Maize A, Maize B, Maize C, Maize D, and Potato1) are indicated in Table 2.
For all recipes, an increase in computed SME (SME_com) led to a drop in intrinsic viscosity ΔVis between 20 to 80%, which can be interpreted by starch depolymerization (Figure 5a). Conversely, no clear trend was observed for the impact of melt temperature (not shown). The trends are in accordance with experimental results (Appendix A, Figure A4a). The effect of viscous dissipation favors the disruption of starch granules and macromolecular chain splitting has been observed many times before, e.g., by Li et al. [3]. The relationships were manifested fairly well by power trends (R2 ≈ 0.7–0.85), the parameters of which depend on starch origin (Figure 5a). Their numerical values are given in Appendix A, Table A2. It can also be seen that, in the SME_com interval [250, 500 kJ/kg], the values of ΔVis rank inversely to their amylose content for maize starches. This is due to the larger molar mass of amylopectin, a highly short-branched macromolecule, which makes it more sensitive to shear, therefore to SME, than amylose, a linear macromolecule (108 vs. 2 × 106 g/mol). The extent of degradation ΔVis is the most important in the case of potato starch, from 60% to 80% in the interval SME_com [250, 350 kJ/kg] (Figure 5a). This trend can be explained by the higher value of intrinsic viscosity for its native starch (450 mL/g) compared to the values given for normal maize starch (amylose content of 23%) in the native state (22.0 mL/g) [31]. Therefore, the larger size (molar mass) of potato starch macromolecules, amplified the sensitivity of this starch to shear.
Since it is possible to predict macromolecular degradation from a computed extrusion variable, it is interesting to predict the conditions for processing maize starch, e.g., having an amylose content of 47%, with a targeted intrinsic viscosity drop ΔVis = 50 ± 5%. According to Figure 5a, SME_com should be in the interval [625, 875 kJ/kg]. Besides operating conditions (Q, N), viscous dissipation, reflected by SME, is sensitive to the profile and geometry of screw elements, mainly restrictive elements. Therefore, simulations were conducted to study the effect of the length of the left-handed screw element: (a) 0 mm, (b) 25 mm, and (c) 50 mm. According to the response surfaces, or operating chart (Figure 5b), at feed MC = 20%, when the restrictive element was absent (0 mm, case a), even by using most severe conditions (highest N and lowest Q), SME_com was lower than 310 kJ/kg, hence, lower than the target. In case b, when its length is 25 mm, the extruder can only operate in a very restrictive interval of Q and N values, represented by the hatched area: N ≥ 650 rpm and Q < 11 kg/h, which leads to SME_com lower than 665 kJ/kg, hence in the target interval (Figure 5c). Finally, for the longest restrictive element (50 mm, case c), the operating window was larger, with feed rate Q possibly reaching 20 kg/h for larger N values (750 rpm) (Figure 5d). This example underlines how the simulation can be used to select the relevant screw profile, which is helpful from a practical viewpoint, e.g., before implementing experimental trials.

3.2.2. Partially Melted Non-Expanded Starches: Study of Starch Structural Modification

Here, we deal with specific products, namely partially melted (partially gelatinized) dense extruded starches, i.e., Potato 2 and pea (PS). This kind of product was obtained through extrusion at the moderate moisture content (30%) and at the low barrel and die temperatures, i.e., 75–95 °C (Table 2, Appendix A, Figure A3). The starch transformation indicator was accessed at two levels: (1) granular one, expressed by the number of granular starch remnants per mm2 (Ng), and (2) molecular one, revealed by the residual molar mass (RM, %) that is obtained by normalization to native starch property.
For both starches, the Ng, reflecting the extent of starch granule breakdown or starch crystallite melting, diminished with an increase in SME_com following exponential functions (Figure 6a). This phenomenon is common for starches extrusion [3]. In contrast with experimental results (Appendix A, Figure A4b), the Ng prediction of extruded pea starch was poor (R2 = 0.56) and the fitting parameters depended on the raw material; this is likely due to the underestimation of SME prediction of pea starch extrusion. This inaccurate prediction comes from approaching the viscous behavior of pea starch melt with pea flour as an input of the 1D global extrusion model. Indeed, a previous study revealed that the shear viscosity of these feeds, determined by a pre-shearing capillary rheometer (Rheoplast®) (Société Courbon, Saint-Etienne, France), operated at SME < 250 kJ/kg, was very close [49]. Prediction improvement could be achieved by determining the viscous behavior of pea starch on a larger SME interval. For potato starch, the Ng dropped by 95%, while SME_com increased by a factor of 4 (R2 = 0.73).
When the SME_com increases, a drop in RM occurs concomitantly with a decrease in Ng (Figure 6a,b). It means that the starch macromolecular chain splitting phenomenon is also involved in the starch crystal melting process; both phenomena being intensified with the more severity of shear stress [32]. The trends, specific for each product, can be fitted fairly well with power functions (R2 = 0.7–0.89). The general tendency is in accordance with experimental results; however it is obviously seen experimentally that potato starch is more sensitive to molecular degradation than pea starch (Appendix A, Figure A4c). The RM of potato product dropped by a factor of 4 when measured SME was increased by a factor 1.5. On the contrary, the RM of pea product diminished by a factor of only 1.5 when the SME_mes was intensified by a factor 4. This trend was not visible for the relationship between RM and computed SME due to underestimation of SME prediction. The higher sensitivity of potato starch macromolecules to shear can be explained by (1) its higher amylopectin content, which leads to a higher molar mass and lesser tendency to entangle and (2) higher extrusion temperature: [Tm, Tm + 10 °C] vs. [Tm−26, Tm−16 °C].
One example of a commercial application of this kind of product is dried noodles. The functional properties of noodles, such as optimal cooking time, water uptake, and cooking loss, depend on the state of starch transformation. During cooking (boiling), the remnant starch granules of noodles uptake water and swell, and cooking loss occurs through the diffusion (leaching) of depolymerized amylose and amylopectin outward the ghost starch granules, from the core to the surface of noodles. Cooking loss can be minimized by allowing starch swelling without excessive breakdown/leaching of its macromolecules. In other words, a compromise between two properties, Ng and RM, shall be sought, based on the assumption that the higher the RM, the more enhanced the leaching.
As an example, we want to manufacture partially gelatinized potato starch with the following features: (1) highly amount of undisrupted starch granule: 150 < Ng < 300; (2) moderate level of RM (RM = 50 ± 10%). According to the relationship described in Figure 6a,b, extrusion should be conducted at SME_com = 525 ± 25 kJ/kg. The operating chart of the extruder is given in Figure 6c. The feasible domain of screw speed and flow rate conditions is presented as a hatched area. The targeted product can be obtained at a low flow rate (0.6 kg/h) by tuning the screw speed N at 40–50 rpm. When a high throughput is required (Q = 1.8 kg/h), the screw speed should be set at a higher velocity: 140–165 rpm.

3.2.3. Legumes Based Snacks: Compromise between Protein Transformation and Expansion

Here, the features of protein-rich snacks, extruded from pea and fava bean ingredients, were represented by protein solubility, in terms of insoluble proteins and density. While the latter reveals the macrostructure, the first can reflect the structure at supramolecular level. The gain in insoluble proteins indicates a higher amount of proteins crosslinked by covalent bonds other than S-S ones, i.e., isopeptide bonds. The higher is the density and the less expanded is the snack.
The rising of melt temperature, at die exit, observed in the pea flour extrusion simulation promoted the protein cross-linking by non S-S bond (R2 ≈ 80, Figure 7a), in accord with experimental results (Appendix A, Figure A4d). The onset of protein cross-linking occurred at a computed melt temperature (T_com) of 155 °C. The snack density decreased from 800 to 100 kg.m3, in other words, expansion was magnified by a factor of eight, when T_com increased from 125 to 175 °C, regardless of formulation (R2 = 0.75, Figure 6b). The positive impact of temperature on legumes expansion, in agreement with experimental results (Appendix A, Figure A4e), is not common for extrusion of starches, where the inverse trend is usually observed [14]. The expansion mechanisms are likely governed by melt rheological changes linked to a certain morphology of starch-protein composite that is probably tuned by structural protein modification [34]. Numerical values of all fitting parameters are summarized in Appendix A, Table A2. The impact of SME_com on snack density and protein transformation is not clear (R2 = 0.45 and 0.55, result not shown).
A low density (lower than 200 kg/m3), resulting from high-temperature extrusion (155 < T_com °C < 175 °C) at low MC (18% and 21%), is indispensable for obtaining crispy snacks. However, extrusion at these temperatures will be accompanied by a higher amount of insoluble proteins (from 5 to 25%). In addition, according to some authors, the protein cross-linking can deteriorate nutritional properties, such as, i.e., protein digestibility [10], and cause a decline in customer acceptance because protein cross-linking increases viscosity of chewed food [9], leading to more difficulty in food swallowing. Therefore, the design of protein-based snacks like pea snacks requires a compromise between density and protein cross-linking.
An example of compromise is set as follows: snacks having a density lower than 200 kg/m3 and insoluble proteins lower than 10% can be extruded at computed temperature (T_com): 155 < T_com °C < 165 (see correlation charts, Figure 7a,b). The operating charts of extruder for extrusion at MC of 21% are given in Figure 7c,d. The hatched area indicates possible combinations of Q and N, for Tb = 90 °C (Figure 7c). and sets of Tb and N, for Q = 20 kg/h (Figure 7d), to achieve targeted T_com. When Tb is fixed at 90 °C, the domain of N for a given Q are 600–650 rpm at the lowest Q (10 kg/h) and 625–700 rpm at the highest Q (50 kg/h). When Q is fixed at 20 kg/h, N can be tuned as: N = 625–700 rpm for low Tb (90 °C) and N = 365–500 rpm for high Tb (140 °C).

3.2.4. Expansion of Wheat Snacks Enriched with Wheat Bran (Fiber)

The features of wheat-based snacks are represented by hydro-solubility of starch (WSIstarch), macrostructure described by sectional (radial) expansion index (SEI), and cellular structure expressed as cell density per cm3 (NC). The higher the WSIstarch, the more soluble are the starch molecules. This variable can thus indicate the state of starch depolymerization induced by extrusion. The higher the NC, the more porous the expanded structure.
Regardless of fiber (bran) content, all wheat flour-based snacks exhibited a unique linear and positive relationship between WSI of starch and SME_com (Figure 8a), in agreement with experimental results (Appendix A, Figure A4f). The higher the shear energy, the higher the extent of starch depolymerization, which is a common trend in starch extrusion [3]. More soluble and smaller starch molecules are more accessible to digestive enzymes during food digestion; however, a highly digested starch can lead to a high value of the glycaemic index.
The prediction of WSIstarch from both computed and measured SMEs is poor (Figure 8a and Figure A4f). If the fiber (wheat bran) is water-soluble, the data dispersion can be likely due to competition in water absorption between starch and fiber; otherwise, the dispersion origin can be explained by filler content (insoluble fiber). The increase of filler content in continuous molten starch can increase melt viscosity, affect viscous dissipation energy, and finally, the degree of starch transformation in the extruder [16]. The prediction could be improved by several means: First, presenting the WSI of total solids, considering fiber and starch, instead of the WSI of starch only. Second, determining the quantity of soluble starch experimentally using a more sophisticated method, such as orcinol titration [34] rather than the gravimetric-based method used by Robin et al. [36]. Third, improving the viscosity model, the effect of filler volume fraction could be taken into account.
The melt temperature and viscosity at the die exit was predicted by simulation since experimental data was not available. Acceptable correlations following exponential functions (R2 = 0.7–0.8) were found between T_com and structure, such as SEI (indicator of expansion in the radial direction) (Figure 8b) and cell density (porosity indicator) (Figure 8c). However, the bran inclusion led to a similar trend but different correlations compared to the flour without bran. We observed that melt temperature and bran enrichment negatively affected SEI. In contrast to SEI, bran inclusion and melt temperature positively affected cell density (NC). For bran-enriched snacks, NC escalated by a factor of 10 when T_com was raised from 125 to 195 °C. The onset of porous structure formation took place at close T_com (155 and 160 °C) for snacks made without and with bran inclusion, respectively. At this onset temperature, bran inclusion surged NC by a factor of 6. Wheat bran can be considered as a filler that can enhance nucleation and hence expansion [16]. An increase in melt temperature intensifies bubble nucleation and bubble growth and thus promotes expansion and porous structure (high NC) [14]. The cell density also negatively correlated to predicted melt viscosity (η_com) (R2 > 0.95, Figure 8d), with distinct power function expressed to each bran content. This trend can mainly be attributed to the negative effect of shear viscosity on the expansion mechanisms such as nucleation, bubble growth, and setting [14]. Higher porosity can make a snack texture more fragile/crispy. Therefore, the design of wheat-based snacks requires a compromise between hydrosolubility of starch and porosity of expanded structure.
As an example, we want to manufacture wheat snacks having the following characteristics: (1) moderately expanded (SEI > 4.5), (2) moderately porous, i.e., medium level of cell density (1000 < NC < 2000), and (3) moderate hydrosolubility of starch (50% < WSIstarch < 60%). From the correlations between extrusion variables (T_com, η_com) and structure (Figure 8b–d) and the relationship between SME_com and WSIstarch (Figure 8a), this kind of snack can be obtained by extruding LB-based recipe, containing 12.6% bran, at T_com of 165–170 °C, η_com < 100 Pa.s, and SME_com > 375 kJ/kg. By examining the feasible region in extruder operating charts (Figure 8e–g), for flow rate Q = 20 kg/h, the targeted T_com, η_com, and SME_com can be achieved by tuning the parameters as follow: a single point of N (700 rpm) for low Tb (125 °C) and small windows of N (630–700 rpm) for higher Tb (140 °C).
In general observation, the 1D global extrusion model is more accurate in the prediction of extruded legumes features than extruded starch ones. The protein structural changes of extruded legumes were predicted from the computed temperature at die exit (Figure 7a,b), while the starch structural changes from computed SME (Figure 5a, Figure 6a,b and Figure 8a). The overestimation of temperature prediction by the extrusion model is much lower than the underestimation of SME prediction (5–10% vs. 35–65%) (Figure 1) and, therefore, a better prediction of extruded legumes features can be expected. The extrusion model neglects the solid transport and solid friction energy and material sticking, leading to SME underestimation. The friction and sticking phenomena depend on the raw materials. It explains why the prediction of intrinsic viscosity of extruded potato starch, through computed SME, is poorer than that of other starches (Figure 5a). Better SME prediction, and thus better product features prediction, could be achieved by improving the extrusion model by taking into account the friction and sticking coefficient.

4. Conclusions

Twin-screw extrusion processing of starchy products was simulated using the 1D global extrusion model (Ludovic® software v7.0.0 Classic Edition (Sciences Computers Consultants, Saint Etienne, France)), and melt rheological model as input. A large data set has been built, including results for extruded products’ properties and extruder variables (temperature and specific mechanical energy), taken as output variables. Satisfactory correlations between process parameters predicted extrusion variables and features of extruded foods, such as structural and functional ones, justify the use of the extrusion model as a computer-aided tool for designing starchy food by extrusion. Moreover, the capacity for predicting these features is extremely vital to reducing time and labor costs in industrial foods R&D.
Flow variables at the die exit, such as melt temperature, specific mechanical energy and melt viscosity, can be computed and used, in turn, as the input variables for a phenomenological model of expansion, allowing predicting expansion indices and cellular structure of expanded starchy foods. The possibility to compute melt viscosity is particularly interesting since it is difficult to measure, either in-line or off-line and because it governs the expansion phenomenon, at least for low hydrated starchy foods. In addition, the knowledge of temperature profile and processing conditions resulting from the simulation is useful in setting extruder configuration and operating parameters for extrusion of thermo-sensitive foods.
As a future prospect, extrusion simulation will be applied to the design of a wide range of starch-protein blends and expanded snacks from pulse crops. Further possible extension includes the processing of plant proteins for meat analogues. Besides classical structural and functional properties, the textural and nutritional features, such as starch and protein digestibility, inactivation of anti-nutritional factors, and destruction of pathogenic microbes and vitamins, will be tackled. The rheology database will be updated with melt viscosity models encompassing a wide range of moisture content, temperatures, and SME, taking into account protein and fiber content, as well as protein structural changes and starch-protein morphology.

Author Contributions

Conceptualization, M.K., F.B. and G.D.V.; methodology, M.K., G.D.V. and F.B.; software, G.D.V.; validation, formal analysis, investigation, resources, data curation, and writing—original draft preparation, M.K.; writing—review and editing, G.D.V. and F.B.; visualization, M.K. and F.B.; supervision, G.D.V.; project administration, G.D.V. All authors agree to be personally accountable for the author’s own contributions and work integrity. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available in Appendix A, and the dataset can also be extracted from the Results Section (Figure 1, Figure 2, Figure 3, Figure 4, Figure 5, Figure 6, Figure 7 and Figure 8, Table 1, Table 2, Table 3 and Table 4).

Acknowledgments

We are thankful to Imen Jebalia (INRAE) for the density measurement of fava bean based-snacks.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations and Nomenclature

n 0 , α i parameters in viscosity model
A , B coefficients depending on screw geometry for an equation expressing volumetric flow rate through a screw channel
AF4asymmetrical flow field flow fractionation
CBcacao-based extruded breakfast cereals
Cp specific heat of hydrated molten starch (J kg−1·°C−1)
Cp s specific heat of anhydrous starch (J·kg−1·°C−1)
Cp w specific heat of water (J·kg−1·°C−1)
Ddie or screw diameter (mm)
DSCdifferential scanning calorimetry
DTE dithioerythritol
E activation energy (J·mol−1)
E v specific energy dissipated by molten starch (kJ·kg−1)
E cd specific energy supplied by conduction from the surfaces of the barrels and screws (kJ·kg−1)
FFfava bean flour
FSCfava bean starch concentrate
HBmalt-based extruded breakfast cereals
Insoltotal amount of in-extractable proteins in SDS and DTE (%)
K flow consistency index (Pa·sn)
K 0 consistency at reference temperature T 0 and reference moisture content MC 0 (Pa sn)
k thermal conductivity of hydrated starch (W·m−1 K−1)
k s thermal conductivity of anhydrous starch (W·m−1 K−1)
k w thermal conductivity of water (W·m−1 K−1)
Ldie or screw length (mm)
L/D length to diameter ratio (dimensionless)
Maize Ahigh amylose maize starch (Eurylon 7)
Maize Bblend of Euroylon 7 and Waxilys 100 with a ratio 2 to 1
Maize Cblend of Euroylon 7 and Waxilys 100 with a ratio 1 to 1
Maize Dwaxy maize starch (Waxilys 100)
MALLSmulti-angle laser light scattering
MCmoisture content (% in wet basis)
MC 0 reference moisture content (%)
n flow behavior index (dimensionless)
N screw speed (rpm)
NCcell density (cm−3)
Ngstarch granule number per unit area (mm−2)
NuNusselt number (dimensionless)
P_mes pressure, measured at die entrance of extruder (Pa)
P_comcomputed pressure at die entrance (Pa)
PFpea flour
Potato Apotato starch studied in Della Valle et al. [31]
Potato Bpotato starch studied in Logié et al. [32]
PSpea starch
R gas constant (J·mol−1·K−1)
SDSsodium dodecyl sulphate
SE-HPLCsize-exclusion high-performance liquid chromatography
RMstarch residual molecular mass (%)
Qtotal feed flow rate (kg·h−1)
Q v volumetric flow rate in a screw channel (m3·s−1)
Q/Nfilling ratio of extruder
SEIsectional expansion index
SMEspecific mechanical energy (kJ·kg−1)
SME_mespredicted SME (kJ·kg−1)
SME_comspecific mechanical energy, retrieved experimentally (kJ·kg−1)
Tmelt temperature (°C)
T a absolute temperature (K)
Tbimposed barrel temperature (°C)
Tdimposed die temperature (°C)
Tmmelting temperature (°C)
T 0 reference temperature (°C)
T_mesmelt temperature, measured at die exit (°C)
T_compredicted melt temperature at die exit (°C)
V free volume of a screw element (m3)
w.b.wet basis
WFrefined wheat flour
WF-LBrefined wheat flour containing low content of wheat bran (12.6%)
WF-HBrefined wheat flour containing high content of wheat bran (24.4%)
W ˙ power dissipated in the volume V of a screw element (Watt)
WSIstarchwater solubility index of starch (in kg/kg dry starch in flour or %)
X w mass fraction of water (kg·kg−1)
X s mass fraction of starch (kg·kg−1)
α water plasticization coefficient (dimensionless)
β thermomechanical history coefficient (kg·kJ−1)
η melt viscosity (Pa·s)
η_compredicted melt viscosity at die exit (Pa·s)
γ ˙ local or true shear rate (s−1)
Δ P / Δ θ pressure variation along angular screw section (Pa)
Δ T temperature change in a screw section (°C)
ΔVisdrop of intrinsic viscosity (%)
Ω screw speed (rad·s−1)
ρ density of snack or molten feed (kg·m−3)
ρ s density of anhydrous starch or flour (kg·m−3)
ρ w density of water (kg·m−3)
ν water volume fraction

Appendix A

Appendix A.1. Extruder Configuration

Figure A1. Die profile used in pea flour extrusion. It consists of a combination of a rheometric die and a cylindrical die. The schema is drawn with Ludovic® software according to description given by Horvat et al. [26]. The rheometric die is a slit-die built up by a triple-step insertion module in a slit channel having a width of 30 mm. The geometry (length and height) of each step is given as follows: step 1—170 × 7 mm2; step 2—120 × 5 mm2, and step 3—80 × 3 mm2. The cylindrical die is built up by a plate with a hollow/cylindrical channel having a diameter of 3 mm and three options of length (5, 10, and 15 mm).
Figure A1. Die profile used in pea flour extrusion. It consists of a combination of a rheometric die and a cylindrical die. The schema is drawn with Ludovic® software according to description given by Horvat et al. [26]. The rheometric die is a slit-die built up by a triple-step insertion module in a slit channel having a width of 30 mm. The geometry (length and height) of each step is given as follows: step 1—170 × 7 mm2; step 2—120 × 5 mm2, and step 3—80 × 3 mm2. The cylindrical die is built up by a plate with a hollow/cylindrical channel having a diameter of 3 mm and three options of length (5, 10, and 15 mm).
Foods 11 01780 g0a1
Figure A2. Illustration of extruder screw profiles for manufacturing of (a) Expanded sample from maize starches having different amylose content [27,33] and potato starch [31]. For potato starch, two cylindrical dies in parallel were used, (b) Dense sample from potato and pea starches [32] and expanded samples from fava bean flour and starch concentrate, and from blends of tapioca, Empure and Merizet starches, (c) expanded snacks from pea flour [34] and breakfast cereals products. For pea flour, the presented extruder was equipped with a cylindrical die only. (d) Expanded snacks from bran-enriched wheat flour [35,36]. All screw representations were build-up using Ludovic® software according to description given in the literatures.
Figure A2. Illustration of extruder screw profiles for manufacturing of (a) Expanded sample from maize starches having different amylose content [27,33] and potato starch [31]. For potato starch, two cylindrical dies in parallel were used, (b) Dense sample from potato and pea starches [32] and expanded samples from fava bean flour and starch concentrate, and from blends of tapioca, Empure and Merizet starches, (c) expanded snacks from pea flour [34] and breakfast cereals products. For pea flour, the presented extruder was equipped with a cylindrical die only. (d) Expanded snacks from bran-enriched wheat flour [35,36]. All screw representations were build-up using Ludovic® software according to description given in the literatures.
Foods 11 01780 g0a2

Appendix A.2. Supplementary Experiment

Melting transition in function of moisture content
The raw materials were hydrated to acquire moisture content between 15 and 30% wet basis, by conditioning them under a controlled relative humidity of 11 to 97% over saturated salt solutions at 20 °C for 2 weeks. The melting point and enthalpy were determined using DSC Q10 (TA Instrument, France) apparatus following the method of Nunez et al. [48] and Logié et al. [32]. The hydrated raw materials were heated from 5 to 200 °C at a rate of 3 °C/min, for breakfast cereals formulations and at 5 °C/min for fava bean ingredients. Melting temperature (Tm) was defined as the maximum temperature of the endothermic peak. The measurements were done in triplicate, with a relative error of 10%.
Density of fava bean-based snacks
The bulk density was obtained using the glass beads displacement method. The extrudate density was calculated from the mass of extrudate and its volume measured through the displacement of glass beads (diameter of 1 mm). The measurements were done in triplicate, with a relative error of 10%.

Appendix A.3. Melting Transition

Figure A3. Inverse of melting temperature (Tm) as a function of water volume fraction v. The points refer to melting temperature obtained as the maximum temperature of an endothermic peak of DSC. The curves represent the best fit using the Flory–Huggins equation [57].
Figure A3. Inverse of melting temperature (Tm) as a function of water volume fraction v. The points refer to melting temperature obtained as the maximum temperature of an endothermic peak of DSC. The curves represent the best fit using the Flory–Huggins equation [57].
Foods 11 01780 g0a3

Appendix A.4. Data Set of Product Properties in Function of Measured Processing Variables

Experimental dataset, with the exception of Figure A4f (fava bean ingredients), was retrieved from the literature (Table 2 and Table 4).
Figure A4. Impact of measured extrusion variables on product features. (a) Variation of the drop of intrinsic viscosity of expanded maize and potato starches with measured SME, (b) variation of remnant granule number per unit area of non-expanded starches with measured SME, (c) variation of residual molecular mass of non-expanded starches with measured SME, (d) variation of insoluble proteins and density (e) of legumes-based snack with measured melt temperature, (f) change of wheat snacks’ water solubility index of starch (WSIstarch) with measured SME.
Figure A4. Impact of measured extrusion variables on product features. (a) Variation of the drop of intrinsic viscosity of expanded maize and potato starches with measured SME, (b) variation of remnant granule number per unit area of non-expanded starches with measured SME, (c) variation of residual molecular mass of non-expanded starches with measured SME, (d) variation of insoluble proteins and density (e) of legumes-based snack with measured melt temperature, (f) change of wheat snacks’ water solubility index of starch (WSIstarch) with measured SME.
Foods 11 01780 g0a4

Appendix A.5. Supplementary Data Set

Table A1. Physical and thermal properties used in the extrusion simulation.
Table A1. Physical and thermal properties used in the extrusion simulation.
PropertiesGranular StateMolten StateReference
Density (kg/m3)700 1 ρ = 1 MC ρ s + MC ρ w [42]
ρ s = 1450 kg/m3
ρ w = 1000 kg/m3
Specific heat (J/(kg °C))1700 Cp = Cp s ( 1 MC ) + Cp w ) MC [42]
Cp w = 4188 kg/m3
Cp s = 1328 + 5.737 T
T in °C
Thermal conductivity (W/(m K)) 1 k = 1 v k s + v k w 1 k = 1 v k s + v k w [41,42,43]
v = X w ρ w X w ρ w + X s ρ s v = X w ρ w X w ρ w + X s ρ s
X i = mass fraction of component i X i = mass fraction of component i
k s = 2.01 × 10 1 + 1.39 × 10 3 T 4.33 × 10 6 T 2 k s = 0.21 + 0.41 × 10 3 ( T + 273 )
k w = 0.571 + 1.76 × 10 3 T 6.70 × 10 6 T 2 k w = 0.571 + 1.76 × 10 3 T 6.70 × 10 6 T 2
Melting latent heat (J/g) 5–20[42,58,59]
Table A2. Relationship between computed and predicted extrusion variables (Figure 1).
Table A2. Relationship between computed and predicted extrusion variables (Figure 1).
FeedsEquationR2
Melt temperature
AllT_com = 1.02 × T_mes0.88
SME
Extruded starches
(Expanded and Dense)
Maizes A, B, C, DSME_com = 0.34 × SME_mes + 156.140.71
Potato 1, Potato 2SME_com = 0.28 × SME_mes + 108.140.96
Expanded snacks
Starch-protein blendsSME_com = 0.61 × SME_mes + 49.440.86
(PF, FF, FSC)
Bran enriched SME_com = 0.31 × SME_mes + 173.450.53
wheat flour
(WF, WF-LB, WF-HB)

References

  1. Ek, P.; Baner, J.M.; Ganjyal, G.M. Chapter 8—Extrusion Processing of Cereal Grains, Tubers, and Seeds. In Extrusion Cooking: Cereal Grains Processing; Ganjyal, G.M., Ed.; Woodhead Publishing: Cambridge, MA, USA, 2020; pp. 225–263. ISBN 978-0-12-815360-4. [Google Scholar]
  2. Markets and Markets. Extruded Snacks Market by Type (Simply Extruded, Expanded, Co-Extruded), Raw Material (Wheat, Potato, Corn, Oats, Rice, Multigrain), Manufacturing Method (Single-Screw, Twin-Screw), Distribution Channel, and Region–Global Forecast to 2026; Markets and Markets: Pune, India, 2020. [Google Scholar]
  3. Li, M.; Hasjim, J.; Xie, F.; Halley, P.J.; Gilbert, R.G. Shear Degradation of Molecular, Crystalline, and Granular Structures of Starch during Extrusion. Starch Stärke 2014, 66, 595–605. [Google Scholar] [CrossRef]
  4. Ralet, M.C.; Della Valle, G.; Thibault, J.F. Raw and Extruded Fibre from Pea Hulls. Part I: Composition and Physico-Chemical Properties. Carbohydr. Polym. 1993, 20, 17–23. [Google Scholar] [CrossRef]
  5. Stanley, W.D. Protein Reactions during Extrusion Processing. In Extrusion Cooking; Mercier, C., Linko, P., Harper, J.M., Eds.; American Association of Cereal Chemists, Inc.: St. Paul, MN, USA, 1989; pp. 321–341. [Google Scholar]
  6. Alam, S.A.; Pentikäinen, S.; Närväinen, J.; Holopainen-Mantila, U.; Poutanen, K.; Sozer, N. Effects of Structural and Textural Properties of Brittle Cereal Foams on Mechanisms of Oral Breakdown and in Vitro Starch Digestibility. Food Res. Int. 2017, 96, 1–11. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  7. Jebalia, I.; Maigret, J.-E.; Réguerre, A.-L.; Novales, B.; Guessasma, S.; Lourdin, D.; Della Valle, G.; Kristiawan, M. Morphology and Mechanical Behaviour of Pea-Based Starch-Protein Composites Obtained by Extrusion. Carbohydr. Polym. 2019, 223, 115086. [Google Scholar] [CrossRef]
  8. Karkle, E.L.; Alavi, S.; Dogan, H. Cellular Architecture and Its Relationship with Mechanical Properties in Expanded Extrudates Containing Apple Pomace. Food Res. Int. 2012, 46, 10–21. [Google Scholar] [CrossRef]
  9. Kristiawan, M.; Della Valle, G.; Réguerre, A.L.; Micard, V.; Salles, C. Artificial Oral Processing of Extruded Pea Flour Snacks. Food Eng. Rev. 2021, 13, 247–261. [Google Scholar] [CrossRef]
  10. Laleg, K.; Barron, C.; Santé-Lhoutellier, V.; Walrand, S.; Micard, V. Protein Enriched Pasta: Structure and Digestibility of Its Protein Network. Food Funct. 2016, 7, 1196–1207. [Google Scholar] [CrossRef]
  11. Li, C.; Kowalski, R.J.; Li, L.; Ganjyal, G.M. Extrusion Expansion Characteristics of Samples of Select Varieties of Whole Yellow and Green Dry Pea Flours. Cereal Chem. J. 2016, 94, 385–391. [Google Scholar] [CrossRef]
  12. Robin, F.; Schuchmann, H.P.; Palzer, S. Dietary Fiber in Extruded Cereals: Limitations and Opportunities. Trends Food Sci. Technol. 2012, 28, 23–32. [Google Scholar] [CrossRef]
  13. Guessasma, S.; Chaunier, L.; Della Valle, G.; Lourdin, D. Mechanical Modelling of Cereal Solid Foods. Trends Food Sci. Technol. 2011, 22, 142–153. [Google Scholar] [CrossRef]
  14. Kristiawan, M.; Chaunier, L.; Della Valle, G.; Ndiaye, A.; Vergnes, B. Modeling of Starchy Melts Expansion by Extrusion. Trends Food Sci. Technol. 2016, 48, 13–26. [Google Scholar] [CrossRef]
  15. Kristiawan, M.; Della Valle, G.; Kansou, K.; Ndiaye, A.; Vergnes, B. Validation and Use for Product Optimization of a Phenomenological Model of Starch Foods Expansion by Extrusion. J. Food Eng. 2019, 246, 160–178. [Google Scholar] [CrossRef]
  16. Kristiawan, M.; Della Valle, G. Chapter 6–Transport Phenomena and Material Changes during Extrusion. In Extrusion Cooking: Cereal Grains Processing; Ganjyal, G.M., Ed.; Woodhead Publishing: Cambridge, MA, USA, 2020; pp. 179–204. ISBN 978-0-12-815360-4. [Google Scholar]
  17. Connelly, R.; Kokini, J. Examination of the Mixing Ability of Single and Twin Screw Mixers Using 2D Finite Element Method Simulation with Particle Tracking. J. Food Eng. 2007, 79, 956–969. [Google Scholar] [CrossRef]
  18. Durin, A.; De Micheli, P.; Nguyen, H.C.; David, C.; Valette, R.; Vergnes, B. Comparison between 1D and 3D Approaches for Twin-Screw Extrusion Simulation. Int. Polym. Process. 2014, 29, 641–648. [Google Scholar] [CrossRef]
  19. Emin, M.A.; Schuchmann, H.P. Analysis of the Dispersive Mixing Efficiency in a Twin-Screw Extrusion Processing of Starch Based Matrix. J. Food Eng. 2013, 115, 132–143. [Google Scholar] [CrossRef]
  20. Emin, M.A.; Karbstein, H. Droplet Breakup and Coalescence in a Twin-Screw Extrusion Processing of Starch Based Matrix. J. Food Eng. 2013, 116, 118–129. [Google Scholar] [CrossRef]
  21. Ficarella, A.; Milanese, M.; Laforgia, D. Numerical Study of the Extrusion Process in Cereals Production: Part I. Fluid-Dynamic Analysis of the Extrusion System. J. Food Eng. 2006, 73, 103–111. [Google Scholar] [CrossRef]
  22. Valette, R.; Coupez, T.; David, C.; Vergnes, B. A Direct 3D Numerical Simulation Code for Extrusion and Mixing Processes. Int. Polym. Process. 2009, 24, 141–147. [Google Scholar] [CrossRef]
  23. Tayeb, J.; Vergnes, B.; Della Valle, G. A Basic Model for a Twin-Screw Extruder. J. Food Sci. 1989, 54, 1047–1056. [Google Scholar] [CrossRef]
  24. Vergnes, B.; Della Valle, G.; Delamare, L. A Global Computer Software for Polymer Flows in Corotating Twin Screw Extruders. Polym. Eng. Sci. 1998, 38, 1781–1792. [Google Scholar] [CrossRef]
  25. Xie, F.; Halley, P.J.; Avérous, L. Rheology to Understand and Optimize Processibility, Structures and Properties of Starch Polymeric Materials. Prog. Polym. Sci. 2012, 37, 595–623. [Google Scholar] [CrossRef]
  26. Horvat, M.; Emin, M.A.; Hochstein, B.; Willenbacher, N.; Schuchmann, H.P. A Multiple-Step Slit Die Rheometer for Rheological Characterization of Extruded Starch Melts. J. Food Eng. 2013, 116, 398–403. [Google Scholar] [CrossRef]
  27. Della Valle, G.; Colonna, P.; Patria, A.; Vergnes, B. Influence of Amylose Content on the Viscous Behavior of Low Hydrated Molten Starches. J. Rheol. 1996, 40, 347–362. [Google Scholar] [CrossRef]
  28. Berzin, F.; Tara, A.; Tighzert, L.; Vergnes, B. Importance of Coupling between Specific Energy and Viscosity in the Modeling of Twin-Screw Extrusion of Starchy Products. Polym. Eng. Sci. 2010, 50, 1758–1766. [Google Scholar] [CrossRef]
  29. Della Valle, G.; Barrès, C.; Plewa, J.; Tayeb, J.; Vergnes, B. Computer Simulation of Starchy Products’ Transformation by Twin-Screw Extrusion. J. Food Eng. 1993, 19, 1–31. [Google Scholar] [CrossRef]
  30. Robin, F.; Engmann, J.; Pineau, N.; Chanvrier, H.; Bovet, N.; Valle, G. Della Extrusion, Structure and Mechanical Properties of Complex Starchy Foams. J. Food Eng. 2010, 98, 19–27. [Google Scholar] [CrossRef]
  31. Della Valle, G.; Boché, Y.; Colonna, P.; Vergnes, B. The Extrusion Behaviour of Potato Starch. Carbohydr. Polym. 1995, 28, 255–264. [Google Scholar] [CrossRef]
  32. Logié, N.; Della Valle, G.; Rolland-Sabaté, A.; Descamps, N.; Soulestin, J. How Does Temperature Govern Mechanisms of Starch Changes during Extrusion? Carbohydr. Polym. 2018, 184, 57–65. [Google Scholar] [CrossRef]
  33. Della Valle, G.; Vergnes, B.; Colonna, P.; Patria, A. Relations between Rheological Properties of Molten Starches and Their Expansion Behaviour in Extrusion. J. Food Eng. 1997, 31, 277–295. [Google Scholar] [CrossRef]
  34. Kristiawan, M.; Micard, V.; Maladira, P.; Alchamieh, C.; Maigret, J.-E.; Réguerre, A.-L.; Emin, M.A.; Della Valle, G. Multi-Scale Structural Changes of Starch and Proteins during Pea Flour Extrusion. Food Res. Int. 2018, 108, 203–215. [Google Scholar] [CrossRef]
  35. Robin, F.; Dubois, C.; Pineau, N.; Schuchmann, H.P.; Palzer, S. Expansion Mechanism of Extruded Foams Supplemented with Wheat Bran. J. Food Eng. 2011, 107, 80–89. [Google Scholar] [CrossRef]
  36. Robin, F.; Théoduloz, C.; Gianfrancesco, A.; Pineau, N.; Schuchmann, H.P.; Palzer, S. Starch Transformation in Bran-Enriched Extruded Wheat Flour. Carbohydr. Polym. 2011, 85, 65–74. [Google Scholar] [CrossRef]
  37. Robin, F.; Bovet, N.; Pineau, N.; Schuchmann, H.P.; Palzer, S. Online Shear Viscosity Measurement of Starchy Melts Enriched in Wheat Bran. J. Food Sci. 2011, 76, E405–E412. [Google Scholar] [CrossRef] [PubMed]
  38. Tayeb, J.; Vergnes, B.; Della Valle, G. Theoretical Computation of the Isothermal Flow through the Reverse Screw Element of a Twin Screw Extrusion Cooker. J. Food Sci. 1988, 53, 616–625. [Google Scholar] [CrossRef]
  39. Noe, I. Etude Des Écoulements de Polymères Dans Une Extrudeuse Bivis Corotative. Ph.D. Thesis, École Nationale Supérieure des Mines de Paris, Paris, France, 1992. [Google Scholar]
  40. Tayeb, J.; Della Valle, G.; Vergnes, B. Modélisation Thermomécanique Du Procédé d’extrusion Bi-Vis. Exemples d’applications. In La Cuisson-Extrusion; Colonna, P., Della Valle, G., Eds.; Lavoisier Tec & Doc: Paris, France, 1994; pp. 489–536. [Google Scholar]
  41. Choi, Y.; Okos, M.R. Effects of Temperature and Composition on the Thermal Properties of Foods. In Food Engineering and Process Applications, Vol. 1. Transport Phenomena; Maguer, M.L., Jelen, P., Eds.; Elsevier Applied Science Publishers: New York, NY, USA, 1986; pp. 93–101. [Google Scholar]
  42. Della Valle, G.; Vergnes, B. Propriétés Thermophysiques et Rhéologiques Des Substrats Utilisés En Cuisson-Extrusion. In La Cuisson-Extrusion; Colonna, P., Della Valle, G., Eds.; Lavoisier Tec & Doc: Paris, France, 1994; pp. 440–467. [Google Scholar]
  43. Muramatsu, Y.; Tagawa, A.; Kasai, T. Effective Thermal Vonductivity of Rice Flour and Whole and Skim Milk Powder. J. Food Sci. 2005, 70, E279–E287. [Google Scholar] [CrossRef]
  44. Barron, C. Déstructuration d’amidons Peu Hydratés Sous Cisaillement. Ph.D. Thesis, Université de Nantes, Nantes, France, 1999. [Google Scholar]
  45. Bouvier, J.M.; Campanella, O.H. The Generic Extrusion Process. III. Thermomechanical Cooking and Food Product Texturation. In Extrusion Processing Technology: Food and Non-Food Biomaterials; John Wiley and Sons.: Hoboken, NJ, USA, 2014; pp. 243–310. [Google Scholar]
  46. Colonna, P.; Buléon, A. Transformations Structurales d’amidons. In La Cuisson-Extrusion; Colonna, P., Della Valle, G., Eds.; Lavoisier Tec & Doc: Paris, France, 1994; pp. 18–43. [Google Scholar]
  47. Logié, N. Modifications Structurales et Comportement Rhéologique d’amidons Faiblement Hydratés Sous Traitement Thermomécanique. Ph.D. Thesis, Université de Lille 1, Lille, France, 2017. [Google Scholar]
  48. Núñez, M.; Sandoval, A.J.; Müller, A.J.; Della Valle, G.; Lourdin, D. Thermal Characterization and Phase Behavior of a Ready-to-Eat Breakfast Cereal Formulation and Its Starchy Components. Food Biophys. 2009, 4, 291–303. [Google Scholar] [CrossRef]
  49. Jebalia, I. Elaboration et Comportement Mécanique de Matériaux Composites Amylo-Protéiques. Ph.D. Thesis, Université de Nantes, Nantes, France, 2020. [Google Scholar]
  50. Brugger, A.; Chaunier, L.; Della Valle, G.; Dellemann, J.-P.; Kristiawan, M. Détermination de Comportements Visqueux Pour l’extrusion de Formulations de Céréales Petit-Déjeuner. Ind. Des Céréales 2013, Aout 2013, 18–22. [Google Scholar]
  51. Vergnes, B.; Della Valle, G.; Tayeb, J. A Specific Slit Die Rheometer for Extruded Starchy Products. Design, Validation and Application to Maize Starch. Rheol. Acta 1993, 32, 465–476. [Google Scholar] [CrossRef]
  52. Barron, C.; Buleon, A.; Colonna, P.; Della Valle, G. Structural Modifications of Low Hydrated Pea Starch Subjected to High Thermomechanical Processing. Carbohydr. Polym. 2000, 43, 171–181. [Google Scholar] [CrossRef]
  53. Rolland-Sabaté, A.; Colonna, P.; Mendez-Montealvo, M.G.; Planchot, V. Branching Features of Amylopectins and Glycogen Determined by Asymmetrical Flow Field Flow Fractionation Coupled with Multiangle Laser Light Scattering. Biomacromolecules 2007, 8, 2520–2532. [Google Scholar] [CrossRef]
  54. Vergnes, B.; Villemaire, J.P. Rheological Behaviour of Low Moisture Molten Maize Starch. Rheol. Acta 1987, 26, 570–576. [Google Scholar] [CrossRef]
  55. Morel, M.H.; Dehlon, P.; Autran, J.C.; Leygue, J.P.; Bar-L’Helgouac’h, C. Effects of Temperature, Sonication Time, and Power Settings on Size Distribution and Extractability of Total Wheat Flour Proteins as Determined by Size-Exclusion High-Performance Liquid Chromatography. Cereal Chem. 2000, 77, 685–691. [Google Scholar] [CrossRef]
  56. Berzin, F.; Tara, A.; Tighzert, L. In-Line Measurement of the Viscous Behaviour of Wheat Starch during Extrusion. Application to Starch Cationisation. Appl. Rheol. 2007, 17, 21221–21222. [Google Scholar] [CrossRef]
  57. Donovan, J.W. Phase Transitions of the Starch–Water System. Biopolymers 1979, 18, 263–275. [Google Scholar] [CrossRef]
  58. Maache-Rezzoug, Z.; Zarguili, I.; Loisel, C.; Queveau, D.; Buléon, A. Structural Modifications and Thermal Transitions of Standard Maize Starch after DIC Hydrothermal Treatment. Carbohydr. Polym. 2008, 74, 802–812. [Google Scholar] [CrossRef] [Green Version]
  59. Liu, H.; Yu, L.; Xie, F.; Chen, L. Gelatinization of Cornstarch with Different Amylose/Amylopectin Content. Carbohydr. Polym. 2006, 65, 357–363. [Google Scholar] [CrossRef]
Figure 1. Comparison between predicted and measured extrusion variables: temperature (a) and specific mechanical energy (SME) (b). The solid lines following linear equations stand for correlation between measured and predicted variables, while the dotted line stands for the equity between them.
Figure 1. Comparison between predicted and measured extrusion variables: temperature (a) and specific mechanical energy (SME) (b). The solid lines following linear equations stand for correlation between measured and predicted variables, while the dotted line stands for the equity between them.
Foods 11 01780 g001
Figure 2. Profiles of screw (a) and computed flow variables along the screws for breakfast cereals (CB -) and (HB,--): (b) product temperature, (c) total dissipated energy, (d) pressure, and (e) viscosity. Processing conditions: flow rate Q = 20 kg/h, screw speed N = 500 rpm, moisture content MC = 15% and barrel temperature profile: 20, 40, 60, 90, 90, 90, 90 °C.
Figure 2. Profiles of screw (a) and computed flow variables along the screws for breakfast cereals (CB -) and (HB,--): (b) product temperature, (c) total dissipated energy, (d) pressure, and (e) viscosity. Processing conditions: flow rate Q = 20 kg/h, screw speed N = 500 rpm, moisture content MC = 15% and barrel temperature profile: 20, 40, 60, 90, 90, 90, 90 °C.
Foods 11 01780 g002
Figure 3. Effect of moisture content on computed temperature (a), total dissipated energy (b), viscosity (c), and pressure (d) along the screw and variation of SME with filling ratio (Q/N) (e). Study case: breakfast cereal: CB formulation. For (ad), processing conditions: Q = 20 kg/h, N = 500 rpm, Tb = 90 °C.
Figure 3. Effect of moisture content on computed temperature (a), total dissipated energy (b), viscosity (c), and pressure (d) along the screw and variation of SME with filling ratio (Q/N) (e). Study case: breakfast cereal: CB formulation. For (ad), processing conditions: Q = 20 kg/h, N = 500 rpm, Tb = 90 °C.
Foods 11 01780 g003
Figure 4. Sensitivity study of the effect of extrusion operating conditions on predicted flow variables for breakfast cereals. HB recipe, MC = 19%. (a) variation of SME_com, (b) melt temperature at die exit T_com and (c) die entrance pressure P_com with screw speed N for various flow rates Q. Last barrel temperature Tb = 90 °C, (d) variation of T_com with Tb for various screw speeds N. Q = 20 kg/h, (e) variation of viscosity η_com with Q for various Tb. N = 500 rpm.
Figure 4. Sensitivity study of the effect of extrusion operating conditions on predicted flow variables for breakfast cereals. HB recipe, MC = 19%. (a) variation of SME_com, (b) melt temperature at die exit T_com and (c) die entrance pressure P_com with screw speed N for various flow rates Q. Last barrel temperature Tb = 90 °C, (d) variation of T_com with Tb for various screw speeds N. Q = 20 kg/h, (e) variation of viscosity η_com with Q for various Tb. N = 500 rpm.
Foods 11 01780 g004
Figure 5. Application of extrusion simulation in the design of expanded maize starches having targeted starch transformation. Step 1. Relationship between the drop of intrinsic viscosity and computed SME (SME_com) (a). The arrow indicates the lower and upper limits of targeted SME_com and ΔVis. Step 2. Creation of extruder operating chart by 1D extrusion model. The effect of length of the left-handed element was studied on the SME_com level: (b) 0 mm, (c) 25 mm, and (d) 50 mm. Step 3. Determination of range of extrusion conditions (Q, N), delimited by hatched area, leading to SME_com between 625 and 875 kJ/kg according to the goal: drop of intrinsic viscosity = 50 ± 5% for feed MC = 20% and amylose content = 47%.
Figure 5. Application of extrusion simulation in the design of expanded maize starches having targeted starch transformation. Step 1. Relationship between the drop of intrinsic viscosity and computed SME (SME_com) (a). The arrow indicates the lower and upper limits of targeted SME_com and ΔVis. Step 2. Creation of extruder operating chart by 1D extrusion model. The effect of length of the left-handed element was studied on the SME_com level: (b) 0 mm, (c) 25 mm, and (d) 50 mm. Step 3. Determination of range of extrusion conditions (Q, N), delimited by hatched area, leading to SME_com between 625 and 875 kJ/kg according to the goal: drop of intrinsic viscosity = 50 ± 5% for feed MC = 20% and amylose content = 47%.
Foods 11 01780 g005
Figure 6. Application of extrusion simulation in the design of dense extruded starches having targeted starch transformation. Step 1. Variation of the remnant granule number per unit area (Ng) (a) and residual molecular mass (RM) (b) with SME_com. The arrows indicate lower and upper limits of targeted SME_com, Ng, and RM. Step 2. Creation of extruder operating chart by the 1D extrusion model and determination of the range of extrusion conditions (Q, N), delimited with a hatched area, leading to SME_com between 500 and 550 kJ/kg according to the goal: RM = 50 ± 10% and 150 ≤ Ng ≤ 300 (c). Fixed conditions: MC = 30%, last barrel temperature = die temperature = 95 °C.
Figure 6. Application of extrusion simulation in the design of dense extruded starches having targeted starch transformation. Step 1. Variation of the remnant granule number per unit area (Ng) (a) and residual molecular mass (RM) (b) with SME_com. The arrows indicate lower and upper limits of targeted SME_com, Ng, and RM. Step 2. Creation of extruder operating chart by the 1D extrusion model and determination of the range of extrusion conditions (Q, N), delimited with a hatched area, leading to SME_com between 500 and 550 kJ/kg according to the goal: RM = 50 ± 10% and 150 ≤ Ng ≤ 300 (c). Fixed conditions: MC = 30%, last barrel temperature = die temperature = 95 °C.
Foods 11 01780 g006
Figure 7. Application of extrusion simulation in the design of pea snacks. Step 1. Variation of in-extractible proteins (a) and density (b) of snacks with T_com. The arrows indicate the upper limits of targeted T_com and snack properties. Step 2. Creation of extruder operating chart by the 1D extrusion model and determination of the range of extrusion conditions (Q, N, Tb), delimited with a hatched area, leading to predicted melt temperature between 155 and 165 °C according to the goals: snack density <200 kg/m3, and protein structural change: insoluble proteins <10% for feed moisture contents of 21% (c,d).
Figure 7. Application of extrusion simulation in the design of pea snacks. Step 1. Variation of in-extractible proteins (a) and density (b) of snacks with T_com. The arrows indicate the upper limits of targeted T_com and snack properties. Step 2. Creation of extruder operating chart by the 1D extrusion model and determination of the range of extrusion conditions (Q, N, Tb), delimited with a hatched area, leading to predicted melt temperature between 155 and 165 °C according to the goals: snack density <200 kg/m3, and protein structural change: insoluble proteins <10% for feed moisture contents of 21% (c,d).
Foods 11 01780 g007
Figure 8. Application of extrusion simulation in the design of wheat snacks. Feed: Wheat flour enriched with low bran content. Step 1. Variation of water solubility index of starch (WSIstarch) with SME_com (a), and variation of sectional expansion index (SEI) (b) and cell density (NC) (c,d) with T_com and η_com. In (a), the arrows indicate the lower limit of targets, while in (c), they indicate lower and upper bounds. Step 2. Creation of extruder operating chart and determination of extrusion parameters range. Fixed operating conditions: MC = 22% and Q = 20 kg/h. The hatched area in (eg) indicate the range of extrusion conditions (N, Tb) leading to T_com and SME_com according to the goals: a porous and expanded structure (1000 ≤ NC ≤ 2000 and SEI ≥ 4.5); and 50% ≤ WSIstarch ≤ 60%.
Figure 8. Application of extrusion simulation in the design of wheat snacks. Feed: Wheat flour enriched with low bran content. Step 1. Variation of water solubility index of starch (WSIstarch) with SME_com (a), and variation of sectional expansion index (SEI) (b) and cell density (NC) (c,d) with T_com and η_com. In (a), the arrows indicate the lower limit of targets, while in (c), they indicate lower and upper bounds. Step 2. Creation of extruder operating chart and determination of extrusion parameters range. Fixed operating conditions: MC = 22% and Q = 20 kg/h. The hatched area in (eg) indicate the range of extrusion conditions (N, Tb) leading to T_com and SME_com according to the goals: a porous and expanded structure (1000 ≤ NC ≤ 2000 and SEI ≥ 4.5); and 50% ≤ WSIstarch ≤ 60%.
Foods 11 01780 g008aFoods 11 01780 g008b
Table 1. Chemical composition of raw material, mostly provided by the supplier.
Table 1. Chemical composition of raw material, mostly provided by the supplier.
NoFormulationSupplierChemical Composition Starch Components
StarchProteinFiberLipidAshSugarSaltAmylose Amylopectin
Expanded starches
1Blends of high amylose Eurylon 7 (Maize A) and waxy maize Waxilys 100 (Maize D) starches [27,33]
Maize ARoquette (France)93.4%0.15%<1%1.3%0.15%n.g.<0.01%70%30%
Maize B: A/D = 2/1 93.0%0.18%<1%1.1%0.13%n.g.<0.01%47%53%
Maize C: A/D = 1/1 92.9%0.20%<1%0.67%0.13%n.g.<0.01%23.5%76.5%
Maize DRoquette (France)92.3%0.25%<1%0.20%0.10%n.g.<0.01%1%99%
2Potato starch: Potato 1 [31]Roquette (France)95.7%0.08%<1%0.06%0.60%n.g.<0.01%21%79%
Dense starch
3Potato starch (Potato 2) [32]Roquette (France)95.7%0.08%<1%0.06%0.6%n.g.<0.01%21%79%
Pea starch [32]Roquette (France)97.7%0.49%<1%0.79%0.5%n.g<0.01%35%65%
4Breakfast cereals (this work; confidential)
Malt extract-based (HB) 82.3%4.6%0.52%1.4%0.17%10.0%1.0%20–28%72–80%
Cacao based (CB) 72.7%6.3%1.3%2.8%0.8%10.0%1.0%20–28%72–80%
Expanded snacks
a) Leguminous starch-protein blends
5Pea flour (PF) [34]Sotexpro (France)46.3%23.9%25.7%2.0%2.1%n.g.n.g.35.0%65.0%
6Fava bean ingredients (this work)
Fava bean flour (FF)Valorex (France) 45.3%32.8%6.3%4.0%3.7%3.0%n.g.33.5%66.5%
Fava bean starch concentrate (FSC)Valorex (France) 59.2%20.6%7.3%1.7%2.7%2.8%n.g.39.9%60.1%
7Cereal based: Wheat flour type 550 enriched with wheat bran [35,36]
WF: Wheat flourProvimi Kliba S.A. (Switzerland)78.5%13.0%2.8%1.1%0.8%n.g.n.g.n.g.n.g.
WF-LB: WF + low bran contentProvimi Kliba S.A. (Switzerland)69.7%13.9%12.6%1.6%1.8%n.g.n.g.n.g.n.g.
WF-HB: WF + high bran contentProvimi Kliba S.A. (Switzerland)55.5%15.0%24.4%2.2%3.2%n.g.n.g.n.g.n.g.
Table 2. Database of extrusion operating conditions and variables.
Table 2. Database of extrusion operating conditions and variables.
Formulation Number and Products Code
Specifications1234
Expanded starchDense extruded starchBreakfast cereal
Maize APotato 1Potato 2CB
Maize B PSHB
Maize C
Maize D
References[27,33][31][32]This work
Extruder brandClextralClextralThermo-scientificCoperion Werner & Pfleiderer
BC 45BC 45HAAKE™ Polysoft OSZSK 26Mc
Screw diameter D (mm)55.555.51625.5
Screw length (L)10001000632740
Screw profileconveying elementsconveying elementsconveying elementsconveying elements
from hopper to diereverse screw elementreverse screw element90° bilobal kneading discsthree blocks of 45° kneading discs
reverse screw elementthree reverse screw elements
conveying elementconveying elements
Die configurationa twin-slittwo circular diesconical diecircular die
rheometric die
Moisture content MC (w.b.)20 < MC < 36%22 < MC < 35%25 and 30%15; 19; 23%
Melt temperature °C135 < T < 190120 < T < 19085 < T < 105confidential
SME (kJ/kg)250 < SME < 2000350 < SME < 1250850 < SME < 5500confidential
Formulation Number and Products Code
Specifications567
Expanded snacks
PFFFWF
FSCWF-LB
WF-HB
References[34]This work[35,36]
Extruder brandCoperion Werner & PfleidererThermo Scientific™Clextral
ZSK 26McProcess 11Evolum 25
Screw diameter D (mm)25.51125
Screw length (L)740269.5400
Screw profileconveying elementsconveying elementsconveying elements
from hopper to diethree blocks of 45° kneading discs90° bilobal kneading discsmixing elements
three reverse screw elementsreverse screw elementreverse screw element
conveying elementsconveying element
Die configuration(a) multiple slit die + circular diecircular diecircular die
(b) circular die
Moisture content MC (w.b.)18% ≤ MC ≤ 26%15; 25; 35%18 and 22%
Melt temperature °C120 < T < 165120 < T < 165n.g.
SME (kJ/kg)450 < SME < 1200135 < SME < 3500n.g.
Table 3. Parameters of viscosity model for various molten starchy products.
Table 3. Parameters of viscosity model for various molten starchy products.
CoefficientsMaize Starch with Amylose Content ofPotatoBreakfast CerealsPea FlourWheat Flour + Wheat Bran *
70%47%23.50%1%Starch
Product CodeMaize AMaize BMaize CMaize DPotato 1CBHBPFWFLBHB
References[27] [31] [50][49][37]
K 0 (Pa·sn)1.21 × 1071.91 × 1064.57 × 1074.72 × 1063.67 × 1061.71 × 1061.22 × 1078.1 × 1058.10 × 10−21.60 × 10−12.00 × 10−3
E R (K)11,440798310,85092356421676511,15342107343755311,201
α 14.617.929.526.114.913.613.711.423.527.434.5
β (kg kJ−1)-7.44 × 10−51.06 × 10−31.65 × 10−33.08 × 10−3---1.67 × 10−31.94 × 10−37.50 × 10−3
n 0 −1.161.813.54−1.021 × 10−2--0.290.1660.1650.129
α 1 (°C)−17.93 × 10−3−7.36 × 10−3−1.54 × 10−27.20 × 10−31.55 × 10−33.28 × 10−33.70 × 10−3----
α 2 1.31−1.93−3.192.549.96 × 10−2−2.18−2.55----
α 3 (kg·kJ−1)-−4.89 × 10−3−8.72 × 10−31.56 × 10−42.75 × 10−4------
α 4 (°C)−1----1.53 × 10−18.05 × 10−39.5 × 10−3----
α 5 (kg·°C−1·kJ−1)-1.98 × 10−54.10 × 10−5-−3.47 × 10−5------
α 6 (kg·kJ−1)-9.19 × 10−31.12 × 10−2-−1.27 × 10−7------
* The authors computed model parameters without taking into account reference variables ( T 0 , MC 0 ). Therefore, the viscosity model took the form as; K = K o exp [ E R ( 1 T ) α   ( MC ) β   ( SME ) ] .
Table 4. Data base of measured structure and properties of extruded products.
Table 4. Data base of measured structure and properties of extruded products.
NoProduct CodeProduct FeaturesMeasurement Method and ReferenceDatabase Source
Expanded starchesStarch destructurationUbbelhode viscometer [27][27]
1Maize AIntrinsic viscosity
Maize BMelt rheology
Maize CShear viscosityOn-line twin-slit rheometer[51][27]
Maize D (Rheopac)
Melting temperatureDSC[52][45,46] *
2Potato 1Intrinsic viscosityUbbelhode viscometer [31][31]
Dense starchStarch depolymerization
3Potato 2Residual molarAsymmetrical flow field [32,53][32]
Pea (PS) massflow fractionation (AF4)
coupled to MALLS
Starch melting
Granule number Light microscopy and image[32]
per unit areaanalysis
Melting temperatureDSC[32]
4Breakfast cerealsMelt rheology
Malt-based (HB)Shear viscosityOff-line preshearing[54][50]
Cacao-based (CB) capillary rheometer
(Rheoplast®)
Melting temperatureDSC[48]This work
Expanded snacks
Legume-based
5Pea flour (PF)Protein aggregation
In-extractible proteinsProtein solubility in [55][34]
SDS + DTE reagents
Structure
DensityBead displacement [34][34]
Melt rheology
Shear viscosityRheoplast® rheometer[54][49]
Melting temperatureDSC[52][34]
6Fava bean ingredientsStructure
FF, FSCDensityBead displacement [34]This work
Melt rheology
Shear viscosityRheoplast® rheometer[54][49]
Melting temperatureDSC[32]This work
Cereal-based
7Wheat flour enriched Structure
with wheat bran Radial expansionMeasurement of diameter[35][35]
WF, WF-LB, WF-HB index (SEI)of extrudate by Caliper
Cellular structureX-ray microtomography and[35][35]
(Cell number 3D image analyze
per cm3)(Component labeling operation)
Functional property
Water solubilitySolubilisation and centrifugation[36][36]
index of starch
Melt rheology
Shear viscosityOn-line twin-slit [37][37]
rheometer (TSAR)
Melting temperatureDSC[52][42] **
* Melting temperature for maize starches having various amylose content was approached by that of other starches found in the literature: waxy maize starch (amylose content of 1%), wheat starch (25%), smooth pea starch (35%), and wrinkled pea starch (75%). ** Melting temperature of bran enriched wheat flour was approached by that of wheat flour found in the literature.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Kristiawan, M.; Della Valle, G.; Berzin, F. Extrusion Simulation for the Design of Cereal and Legume Foods. Foods 2022, 11, 1780. https://doi.org/10.3390/foods11121780

AMA Style

Kristiawan M, Della Valle G, Berzin F. Extrusion Simulation for the Design of Cereal and Legume Foods. Foods. 2022; 11(12):1780. https://doi.org/10.3390/foods11121780

Chicago/Turabian Style

Kristiawan, Magdalena, Guy Della Valle, and Françoise Berzin. 2022. "Extrusion Simulation for the Design of Cereal and Legume Foods" Foods 11, no. 12: 1780. https://doi.org/10.3390/foods11121780

APA Style

Kristiawan, M., Della Valle, G., & Berzin, F. (2022). Extrusion Simulation for the Design of Cereal and Legume Foods. Foods, 11(12), 1780. https://doi.org/10.3390/foods11121780

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

Article Metrics

Back to TopTop