The processing of soybean meal, a primary by-product of soybean oil extraction and a crucial protein source in animal feed, presents significant challenges in achieving optimal nutritional quality. The inherent anti-nutritional factors and the need for improved protein digestibility necessitate advanced processing techniques. Extrusion technology, particularly using multi-screw extruders, has emerged as a key method for enhancing the functional and nutritional properties of soybean meal. While triple-screw extruders offer advantages over twin-screw counterparts in terms of processing capacity, mixing area, and handling of high-viscosity materials, traditional designs often fall short in delivering the intensive shear and distributive mixing required for homogeneous product quality.
The primary defect of conventional co-rotating triple-screw extruders lies in their relatively uniform and laminar flow patterns, which limit the residence time distribution and the intensity of shear deformation imposed on the material. This can lead to insufficient breakdown of components and uneven mixing. To overcome this limitation, this analysis focuses on the design and evaluation of a novel variable-speed triple-screw extruder. The core innovation involves the integration of a coaxial speed-changing mechanism based on bevel gears and the strategic placement of specialized mixing elements. The objective is to deliberately manipulate the flow field—introducing flow stagnation, reversal, and reorientation—to significantly enhance shear and dispersive mixing performance for soybean meal processing.
The design and evaluation process follows a structured computational engineering approach. First, critical screw geometry parameters are optimized. Then, a detailed three-dimensional model of the flow domain is created. Using Computational Fluid Dynamics (CFD) based on finite element methods, the non-Newtonian flow of soybean meal is simulated. The results from the novel variable-speed design are systematically compared against a baseline conventional triple-screw extruder through analysis of pressure fields, velocity vector fields, velocity streamlines, and axial velocity profiles. Finally, experimental validation is discussed to corroborate the simulation findings.
1. Geometric Model Optimization and Design
The performance of an extruder screw is fundamentally governed by its geometric parameters, which directly influence conveying capacity, shear rate, and residence time. For a given material and screw speed, the volumetric conveying capacity in the feed section can be expressed as a function of the screw’s cross-sectional area and its axial feed velocity.
The volumetric flow rate \( Q \) is given by:
$$ Q = A \cdot v $$
where \( \rho \) is material density, \( A \) is the effective cross-sectional area of the screw channel, and \( v \) is the axial feed velocity.
The cross-sectional area \( A \) for a screw with a rectangular channel approximation is:
$$ A = \frac{\pi}{4} (D^2 – d^2) – \frac{\delta H}{\sin \alpha} $$
where \( D \) is the screw outer diameter, \( d \) is the screw root diameter, \( \delta \) is the flight width, \( H \) is the channel depth, and \( \alpha \) is the helix angle.
The axial feed velocity \( v \) is related to the screw rotation speed \( n \) and the pitch \( S \) by:
$$ v = \frac{S \cdot n}{60} $$
The helix angle \( \alpha \) is related to the pitch and outer diameter:
$$ \alpha = \arctan\left(\frac{S}{\pi D}\right) $$
Combining these equations, the conveying capacity can be modeled as:
$$ Q = k \cdot \left[ \frac{\pi}{4}\left(D^2 – (D – 2H)^2\right) – \frac{\delta H}{\sin\left(\arctan\frac{S}{\pi D}\right)} \right] \cdot S $$
where \( k \) is a constant incorporating rotational speed and material factors.
For optimization, the channel depth \( H \) and pitch \( S \) were selected as design variables. The Particle Swarm Optimization (PSO) algorithm was employed due to its rapid convergence, simple implementation, and effectiveness in navigating parameter spaces. The algorithm was set to find the combination of \( H \) (range: 1–8 mm) and \( S \) (range: 5–25 mm) that maximizes \( Q \) within the mechanical constraints of a triple-screw configuration.
The PSO process iteratively updates a swarm of candidate solutions. The fitness of each particle (parameter set) is evaluated using the derived \( Q \) function. Individual and global best positions guide the swarm’s movement through the design space until convergence. The optimization yielded an optimal channel depth of 7.4 mm and a pitch of 19.6 mm. These values were used for the screw design in both the conventional and the novel extruder for a consistent baseline comparison of conveying capacity.
| Parameter | Symbol | Value | Unit |
|---|---|---|---|
| Total Screw Length | L | 300 | mm |
| Outer Diameter | D | 80 | mm |
| Root Diameter | d | 65.2 | mm |
| Optimal Channel Depth | H | 7.4 | mm |
| Optimal Pitch | S | 19.6 | mm |
| Flight Width | δ | 8 | mm |
2. Novel Extruder Design with Integrated Bevel Gear Transmission
The novel extruder design features a segmented screw configuration assembled on a single shaft. The total length is divided into three functional zones: a conventional conveying zone, a mixing zone with kneading blocks, and a variable-speed conveying zone. The key innovation resides within the variable-speed zone.
This zone houses an internally mounted coaxial transmission system utilizing bevel gears. The mechanism consists of two driving bevel gears mounted on the central shaft, each meshing with a larger driven bevel gear connected to the screw profile in that zone. The gear train is completely enclosed within the barrel. With a designed gear ratio of 2:1 (driven gear teeth to driving gear teeth), this system achieves two critical functions simultaneously: speed reduction and direction reversal of the screw rotation within that specific segment.
| Component | Number of Teeth | Module | Function |
|---|---|---|---|
| Driving Bevel Gear (on shaft) | 17 | 1.25 | Input from main drive |
| Driven Bevel Gear (on screw segment) | 34 | Output to variable-speed segment | |
| Gear Ratio | 2:1 (34/17) | Halves speed, reverses direction | |
The mixing zone, located prior to the variable-speed zone, is equipped with neutral kneading disc blocks. These blocks are arranged with a staggered angle (e.g., 60°) and have a profile matching the triple-screw cross-section. They provide intense distributive mixing through repeated splitting and recombination of the material stream but contribute negligible axial pumping force.
Thus, the material progression encounters: 1) Forward conveying and mild shear, 2) Intensive distributive mixing with flow retardation in the kneading block zone, 3) A sudden reduction in forward conveying velocity and the onset of flow hindrance due to the reversed rotation in the variable-speed zone, effectively acting as a dynamic restrictive element, and 4) Resumed forward conveying in the final section. This sequenced disruption is designed to maximize total shear strain and residence time.
3. Mathematical Modeling and Governing Equations
To simulate the flow, soybean meal in the molten state within the extruder barrel is modeled as an incompressible, non-Newtonian fluid under isothermal, laminar flow conditions. The following assumptions are made: the flow is fully developed and fills the channel; inertia and body forces are negligible; and no-slip conditions exist at all solid boundaries (screw and barrel surfaces).
The flow is governed by the conservation equations of mass and momentum. The constitutive equation for a power-law fluid model is used to represent the shear-thinning behavior typical of biopolymer melts like soybean meal.
Continuity Equation (Conservation of Mass):
$$ \nabla \cdot \vec{v} = 0 $$
or, in Cartesian components:
$$ \frac{\partial v_x}{\partial x} + \frac{\partial v_y}{\partial y} + \frac{\partial v_z}{\partial z} = 0 $$
where \( \vec{v} \) is the velocity vector.
Cauchy Momentum Equation (Conservation of Momentum):
$$ -\nabla p + \nabla \cdot \vec{\tau} = 0 $$
where \( p \) is the isotropic pressure and \( \vec{\tau} \) is the deviatoric stress tensor. In component form:
$$
\begin{aligned}
-\frac{\partial p}{\partial x} + \frac{\partial \tau_{xx}}{\partial x} + \frac{\partial \tau_{yx}}{\partial y} + \frac{\partial \tau_{zx}}{\partial z} &= 0 \\
-\frac{\partial p}{\partial y} + \frac{\partial \tau_{xy}}{\partial x} + \frac{\partial \tau_{yy}}{\partial y} + \frac{\partial \tau_{zy}}{\partial z} &= 0 \\
-\frac{\partial p}{\partial z} + \frac{\partial \tau_{xz}}{\partial x} + \frac{\partial \tau_{yz}}{\partial y} + \frac{\partial \tau_{zz}}{\partial z} &= 0
\end{aligned}
$$
Power-Law Constitutive Equation:
$$ \vec{\tau} = \mu(\dot{\gamma}) \cdot \dot{\vec{\gamma}} $$
$$ \mu(\dot{\gamma}) = m \cdot \dot{\gamma}^{\,n-1} $$
where \( \dot{\vec{\gamma}} \) is the shear rate tensor, \( \mu \) is the apparent viscosity, \( m \) is the consistency coefficient, and \( n \) is the power-law index. The magnitude of the shear rate is \( \dot{\gamma} = \sqrt{\frac{1}{2} (\dot{\vec{\gamma}}:\dot{\vec{\gamma}})} \).
| Parameter | Value | Unit |
|---|---|---|
| Material | Molten Soybean Meal | – |
| Density (ρ) | 665 | kg/m³ |
| Consistency Coefficient (m) | 1930 | Pa·sⁿ |
| Power-Law Index (n) | 0.5 (assumed) | – |
| Barrel Temperature | 140 | °C |
| Conventional Zone Screw Speed | 120 | rpm |
| Variable-Speed Zone Screw Speed (via bevel gears) | 60 (reverse direction) | rpm |
| Inlet Velocity | 0.05 | m/s |
| Outlet Pressure | 3 | MPa |
| Barrel Wall Condition | No-Slip (0 m/s) | – |
| Screw Surface Condition | No-Slip (Rotating speed defined per zone) | – |
4. CFD Simulation Results and Analysis
The three-dimensional flow domain for both the conventional and the novel variable-speed triple-screw extruder was meshed using tetrahedral elements. A mesh independence study was conducted to ensure solution accuracy. The ANSYS CFX solver was used to solve the system of governing equations with the specified boundary conditions.
4.1 Pressure Field Analysis
The pressure field is a direct indicator of the screw’s pumping or build-up capability. A steady, progressive pressure rise along the extrusion direction is characteristic of efficient conveying.
In the conventional triple-screw extruder, the pressure field shows a smooth, monotonic increase from inlet to outlet. The isobars are relatively uniform and parallel to the screw flights, indicating a steady, predictable pumping action with minimal flow disruption. This is efficient for conveying but suggests limited mixing enhancement.
In the novel variable-speed extruder, the pressure field exhibits distinct zones corresponding to its functional segments. The pressure increases in the initial conveying section. Upon entering the neutral kneading block zone, the pressure gradient flattens significantly, as these elements provide mixing but not positive displacement. The most critical change occurs in the variable-speed zone. Here, due to the action of the bevel gear transmission which reverses screw rotation, this segment functions as a dynamic restriction or a “reverse screw.” Consequently, the pressure decreases across this zone, as the reversed rotation actively opposes the forward flow, requiring the upstream sections to build sufficient pressure to overcome this barrier. Finally, in the last conveying segment, the pressure increases again towards the die. This non-monotonic pressure profile—featuring a plateau and a local drop—directly demonstrates the flow manipulation introduced by the design, promoting longer residence time and increased shear work.
| Metric | Conventional Triple-Screw | Novel Variable-Speed Triple-Screw | Implication |
|---|---|---|---|
| Pressure Profile | Smooth, monotonic increase | Non-monotonic: Rise, plateau, local drop, final rise | Variable-speed zone creates a dynamic flow restriction. |
| Axial Velocity in Conveying Zones | ~0.12 m/s (steady) | ~0.12 m/s (forward zones), ~ -0.06 m/s (variable-speed zone) | Bevel gears cause speed halving and reversal, inducing backflow. |
| Flow Streamlines | Continuous, laminar, aligned with flights | Disrupted: Stagnation in kneading zone, folding/breaking in variable-speed zone | Greatly enhanced distributive mixing and interfacial renewal. |
| Primary Mixing Mechanism | Laminar shear, limited reorientation | Distributive mixing (splitting/recombination) + Dispersive mixing (reversed flow high shear) | Synergistic effect of kneading blocks and bevel gear-driven speed differential. |
4.2 Velocity Field and Streamline Analysis
The velocity field provides insight into local shear rates and material trajectory. The axial velocity component is particularly telling for understanding conveying efficiency and backflow.
A plot of axial velocity along a line through the center of the flow domain reveals stark differences. For the conventional screw, the axial velocity remains positive and fairly constant in conveying sections, with minor dips at intermeshing regions. For the novel screw, the velocity drops to near zero in the kneading block zone (no conveying), plummets to a significant negative value (approximately -0.06 m/s) in the bevel gear-driven variable-speed zone, and recovers to a positive value in the final section. This negative axial velocity is unequivocal evidence of backflow, where material is pushed backward relative to the extrusion direction, creating a powerful shearing and remixing region.
Velocity vector plots and streamlines offer a visual representation of mixing quality. In the conventional extruder, streamlines are smooth, continuous, and largely follow the helical path of the screws, indicating a predominance of laminar shear with little chaotic motion. In the novel extruder, the streamlines become highly complex. In the kneading block zone, they show multiple splitting and recombining paths, indicative of good distributive mixing. Most dramatically, at the entrance to the variable-speed zone, the sudden reversal of screw motion causes the streamlines to fold, break, and form localized vortices. This flow separation and re-circulation is a hallmark of efficient dispersive mixing, as it repeatedly stretches and folds fluid elements, exposing more material to high shear stress and improving homogeneity.
The synergy between the neutral kneading blocks and the bevel gear-activated variable-speed segment is clear. The kneading blocks first slow down and mix the material distributively. The subsequent variable-speed segment, with its reversed rotation, then acts on this pre-mixed material, subjecting it to intense shear and elongational flows due to the conflicting directions of motion, further homogenizing the melt.
5. Experimental Validation
To validate the computational findings, a prototype variable-speed screw segment incorporating the bevel gear mechanism was manufactured and tested against a conventional screw in a triple-screw extruder setup under controlled conditions. Soybean meal was used as the feedstock. The primary measured response was the mean residence time, a critical parameter directly related to mixing effectiveness.
The experiment was conducted at various screw speeds while maintaining constant feed rate (5 kg/h) and barrel temperature (140°C). A tracer method was used to determine the time from feeding to extrusion. The results consistently showed that the novel extruder configuration with the variable-speed segment yielded significantly longer residence times compared to the conventional configuration at the same nominal screw speed.
| Main Screw Speed (rpm) | Residence Time – Conventional (s) | Residence Time – Novel Variable-Speed (s) | Increase (%) |
|---|---|---|---|
| 80 | 45.5 | 54.5 | 19.8 |
| 100 | 35.7 | 44.6 | 24.9 |
| 120 | 30.9 | 39.7 | 28.5 |
| 140 | 22.4 | 29.1 | 29.9 |
| 160 | 17.7 | 21.2 | 19.8 |
The data confirms the core simulation prediction: the integrated design, centered on the bevel gear transmission, successfully extends material residence time by 20% to 30%. This increased dwell time within the shear-intensive zones (kneading blocks and reversed-flow segment) directly correlates with improved mixing and shear treatment, leading to a more uniform and potentially higher-quality extruded soybean meal product.
6. Conclusion
This comprehensive analysis demonstrates the significant enhancement in mixing and shear performance achievable in a triple-screw extruder through the innovative integration of a coaxial variable-speed mechanism based on bevel gears. The systematic approach—encompassing geometric parameter optimization via PSO, detailed 3D CFD simulation, and experimental validation—provides a robust evaluation of the design.
The key findings are:
- The optimized screw geometry (H=7.4 mm, S=19.6 mm) provides a baseline for efficient conveying.
- The incorporation of neutral kneading blocks effectively slows the material flow and promotes distributive mixing through repeated flow division and recombination.
- The central innovation, the bevel gear transmission system, successfully creates a coaxial speed differential and directional reversal in a downstream screw segment. This segment acts as a dynamic restriction, generating measurable backflow (negative axial velocity ~ -0.06 m/s) and complex, turbulent-like flow patterns including vortex formation and streamline breaking.
- CFD results clearly show the novel design transforms the pressure field from monotonic to non-monotonic and the flow field from laminar-dominated to chaos-enhanced, indicating superior mixing.
- Experimental tests confirm a 20-30% increase in mean residence time, physically validating the simulated flow retardation and enhanced mixing efficacy.
Therefore, the variable-speed triple-screw extruder utilizing bevel gears presents a mechanically elegant and functionally superior alternative to conventional designs for applications like soybean meal processing, where intensive shear and homogeneous mixing are critical for product quality. The principles demonstrated here—using internal gearing like bevel gears to locally manipulate screw kinematics—offer a powerful strategy for advanced polymer and food processing extruder design.