TWO-DIMENSIONAL NONISOTHERMAL NUMERICAL ANALYSIS OF SPEED RATIO EFFECTS ON DISPERSION AND DISTRIBUTION IN HIGH-VISCOSITY PARTIALLY FILLED RUBBER MIXING

2019 ◽  
Vol 92 (1) ◽  
pp. 168-185
Author(s):  
Istiaque Ahmed ◽  
Hari Poudyal ◽  
Abhilash J. Chandy

ABSTRACT Two-dimensional, transient, and nonisothermal computational fluid dynamics simulations are conducted for high-viscosity rubber mixing in a two-wing rotor-equipped partially filled chamber of fill factor 75%. Calculations presented assess the effect of three differential speeds or speed ratios of the two rotors for the rubber mixing process: 1.0 (also called even speed), 1.125, and 1.5. A Eulerian multiphase model, the volume of fluid technique, is employed to simulate two different phases, rubber and air, by calculating the free surface between the two phases, in addition to the main governing equations such as the continuity, momentum, and energy equations. To characterize the non-Newtonian, highly viscous rubber under nonisothermal conditions, the shear rate–dependent Carreau-Yasuda model along with an Arrhenius function are employed. A set of massless particles is introduced into the chamber to calculate several parameters related to dispersive and distributive mixing characteristics. Specifically, the mixing index and maximum shear stress are analyzed for the dispersive nature, whereas cluster distribution index and length of stretch are calculated for investigating the distributive nature of the mixing process. Also, the temporal viscous heat generation rate, a good indication of the temperature rise throughout the domain, which is critical in the process and equipment design, is analyzed here. Results showed that the 1.125 speed ratio was the most efficient in terms of distributive mixing and heat generation.

2019 ◽  
Vol 92 (1) ◽  
pp. 152-167 ◽  
Author(s):  
Hari Poudyal ◽  
Suma R. Das ◽  
Abhilash J. Chandy

ABSTRACT A finite volume technique is used to analyze the isothermal and non-isothermal flow behavior for the rubber mixing process in a two-dimensional, partially filled (75%) internal mixer, which consists of two counterrotating rotors rotating at 20 rpm. In order to capture the interface between air and rubber, an Eulerian multiphase model called volume of fluid (VOF) has been employed here. The transient flow behavior was accomplished by a sliding mesh technique, and the highly viscous, non-Newtonian properties of the rubber have been characterized using the Bird–Carreau model. Most of the previous computational fluid dynamic (CFD)-based investigations of rubber mixing assumed isothermal flow, and consequently negligible viscous heat generation, temperature rise, and viscosity drop associated with heat generation. Hence, a non-isothermal simulation is carried out, and results are compared with those of an equivalent isothermal simulation. In addition, dispersive and distributive mixing characteristics are assessed using statistics calculated from particle tracks generated by a set of massless and neutral particles that have been injected in the simulation. For this purpose, quantities such as the cumulative distribution of maximum shear stress, length of stretch, and cluster distribution index are calculated and compared between isothermal and non-isothermal conditions. Results showed a significant difference between the isothermal and non-isothermal simulations, thus making the isothermal assumption critical. Also, the non-isothermal simulation predicted better mixing during the entire mixing cycle.


2016 ◽  
Vol 89 (3) ◽  
pp. 371-391 ◽  
Author(s):  
Suma R. Das ◽  
Pashupati Dhakal ◽  
Hari Poudyal ◽  
Abhilash J. Chandy

ABSTRACT Three-dimensional, transient, isothermal, and incompressible computational fluid dynamics (CFD) simulations are carried out for rubber mixing with two counter-rotating rotors in a partially filled chamber in order to assess the effect of different speed ratios. The three different speed ratios that are investigated include 1.0, 1.125, and 1.5. In addition to the solution of the incompressible continuity and momentum equations, a Eulerian multiphase model is employed to simulate two phases, rubber and air, and the volume of fluid (VOF) technique is used to calculate the free surface flow between the phases. The Bird–Carreau model is used to characterize the non-Newtonian highly viscous rubber. Massless particles are injected in the simulations to obtain data required for statistical calculations related to dispersive and distributive mixing characteristics. Specifically, joint probability density functions of mixing index and shear rate, and cumulative distribution functions of maximum shear stress are calculated to assess dispersive mixing, while distributive mixing capabilities are evaluated using various quantities such as cluster distribution index, axial distribution, interchamber particle transfer, and segregation scale. Results showed the speed ratio 1.125 to be consistently superior to 1.5 and 1.0, in terms of both dispersive and distributive mixing performance. The large speed difference between the rotors in the case of 1.5 caused it to perform the worst.


1992 ◽  
Vol 296 ◽  
Author(s):  
Robert S. Sinkovits ◽  
Lee Phillips ◽  
Elaine S. Oran ◽  
Jay P. Boris

AbstractThe interactions of shocks with defects in two-dimensional square and hexagonal lattices of particles interacting through Lennard-Jones potentials are studied using molecular dynamics. In perfect lattices at zero temperature, shocks directed along one of the principal axes propagate through the crystal causing no permanent disruption. Vacancies, interstitials, and to a lesser degree, massive defects are all effective at converting directed shock motion into thermalized two-dimensional motion. Measures of lattice disruption quantitatively describe the effects of the different defects. The square lattice is unstable at nonzero temperatures, as shown by its tendency upon impact to reorganize into the lower-energy hexagonal state. This transition also occurs in the disordered region associated with the shock-defect interaction. The hexagonal lattice can be made arbitrarily stable even for shock-vacancy interactions through appropriate choice of potential parameters. In reactive crystals, these defect sites may be responsible for the onset of detonation. All calculations are performed using a program optimized for the massively parallel Connection Machine.


2015 ◽  
Vol 82 (7) ◽  
Author(s):  
H. Song ◽  
R. J. Dikken ◽  
L. Nicola ◽  
E. Van der Giessen

Part of the friction between two rough surfaces is due to the interlocking between asperities on opposite surfaces. In order for the surfaces to slide relative to each other, these interlocking asperities have to deform plastically. Here, we study the unit process of plastic ploughing of a single micrometer-scale asperity by means of two-dimensional dislocation dynamics simulations. Plastic deformation is described through the generation, motion, and annihilation of edge dislocations inside the asperity as well as in the subsurface. We find that the force required to plough an asperity at different ploughing depths follows a Gaussian distribution. For self-similar asperities, the friction stress is found to increase with the inverse of size. Comparison of the friction stress is made with other two contact models to show that interlocking asperities that are larger than ∼2 μm are easier to shear off plastically than asperities with a flat contact.


Author(s):  
Bharti bharti ◽  
Debabrata Deb

We use molecular dynamics simulations to investigate the ordering phenomena in two-dimensional (2D) liquid crystals over the one-dimensional periodic substrate (1DPS). We have used Gay-Berne (GB) potential to model the...


Polymers ◽  
2018 ◽  
Vol 10 (8) ◽  
pp. 926 ◽  
Author(s):  
Richard Schwarzl ◽  
Roland Netz

We study collapsed homo-polymeric molecules under linear shear flow conditions using hydrodynamic Brownian dynamics simulations. Tensile force profiles and the shear-rate-dependent globular-coil transition for grafted and non-grafted chains are investigated to shine light on the different unfolding mechanisms. The scaling of the critical shear rate, at which the globular-coil transition takes place, with the monomer number is inverse for the grafted and non-grafted scenarios. This implicates that for the grafted scenario, larger chains have a decreased critical shear rate, while for the non-grafted scenario higher shear rates are needed in order to unfold larger chains. Protrusions govern the unfolding transition of non-grafted polymers, while for grafted polymers, the maximal tension appears at the grafted end.


1993 ◽  
Vol 321 ◽  
Author(s):  
M. Li ◽  
W. L. Johnson ◽  
W. A. Goddard

ABSTRACTThermodynamic properties, structures, defects and their configurations of a two-dimensional Lennard-Jones (LJ) system are investigated close to crystal to glass transition (CGT) via molecular dynamics simulations. The CGT is achieved by saturating the LJ binary arrays below glass transition temperature with one type of the atoms which has different atomic size from that of the host atoms. It was found that for a given atomic size difference larger than a critical value, the CGT proceeds with increasing solute concentrations in three stages, each of which is characterized by distinct behaviors of translational and bond-orientational order correlation functions. An intermediate phase which has a quasi-long range orientational order but short range translational order has been found to exist prior to the formation of the amorphous phase. The destabilization of crystallinity is observed to be directly related to defects. We examine these results in the context of two dimensional (2D) melting theory. Finite size effects on these results, in particular on the intermediate phase formation, are discussed.


Sign in / Sign up

Export Citation Format

Share Document