Fill Factor Effects in Highly-Viscous Non-Isothermal Rubber Mixing Simulations

Abstract

A finite volume technique in a commercial computational fluid dynamics (CFD) code is employed in this study to simulate transient, incompressible, non-Newtonian and non-isothermal rubber mixing. The simulation processes are conducted in a two-dimensional(2D) domain, where a mixing chamber partially-filled with rubber is equipped with a pair of two-wing non-intermeshing counter-rotating rotors. The main objective is to assess the effect of different fill factors of rubber on dispersive and distributive mixing characteristics by simulating 15 revolutions of the rotors rotating at 20 min−1. 50%, 60%, 70%, 75%, 80% and 90% are the six different fill factors chosen for the study. An Eulerian multiphase method has been applied to solve for the two different phases, rubber and air. The non-Newtonian property of rubber is handled using the shear rate dependent Carreau-Yasuda model, along with an Arrhenius function to include the temperature dependency. In addition to the governing equations related to the conservation of mass, momentum and energy, the volume of fluid (VOF) method is chosen to track the interface between air and rubber. With regard to the results, flow patterns, thermal distributions, viscosity behavior and volume fraction are analyzed for the different fill factors. In addition, dispersive and distributive mixing behavior is also assessed in detail using Lagrangian statistics, such as mixing index, cumulative distribution of maximum shear stress, cluster distribution index (CDI), scale of segregation (SOS) and length of stretch (LOS), calculated from massless particles. Both the Eulerian and Lagrangian results showed that fill factors between 70% and 80% presented the most reasonable and efficient mixing scenario, and also exhibited the best dispersive and distributive mixing characteristics combined.

Non-isothermal non-Newtonian fluids: the stationary case

“…Remark 2.1. Notice that the definition is well posed, since, when d = 3, φ ∈ V 0 ∩W 1,3 (Ω) ensures that the convective term in (2.6) is finite, by the embeddings V p div → L q (Ω), for some q > 2, and V 0 → L 6 (Ω). On the contrary, in the case d = 2 it is enough to consider φ ∈ V 0 , by the embeddings V p div → L q (Ω), for some q > 2, and V 0 → L r (Ω) for any r ∈ [2, ∞).…”

“…Moreover, most of them are carried out under non-isothermal conditions. Motivated by these important applications, in the present paper we focus on a class of nonisothermal non-Newtonian fluid models (see, for instance, [3,4,25,32,37,35,13] and references therein). We recall that the description of a non-Newtonian fluid behavior is based on a power-law ansatz.…”

Numerical Optimization Simulation of Synchronous Four-Wing Rotor

“…For polymer melt flow, the basic governing equation can be expressed as follows [20]: The conservation of continuity equation:…”

Numerical Modeling of the Mixing of Highly Viscous Polymer Suspensions in Partially Filled Sigma Blade Mixers

“…One of the pioneers within the field of mixing is Ica Manas-Zlaczower [6][7][8][9][10][11]. Her famous dispersive mixing index from 1992 is still used today [12][13][14][15][16][17], along with her distributive Cluster Distribute Index [12,17,18]. There have not been many alternatives to Zlaczower's dispersive mixing index.…”

“…Heating is an important aspect of mixing, since it can reduce viscosity drastically [12,17,24,25]. The advantage of reducing the viscosity is that it leads to an increased Reynolds Number, which improves mass transfer (i.e., mixing) [25].…”