My PhD work consists of two simulation projects. One is to study the effect of polydispersity on the dense 2D granular system under steady shear. We use the Durian bubble model together with Lees-Edwards boundary conditions to generate the shear on our systems with exponential size distributions with various size spans. Then we compare the results with conventional bidisperse system. Shear produces a mean affine flow, and nonaffine plastic deformations resulting from local rearrangements. We calculate the deviation from the affine flow to quantify the nonaffinity for individual particles. We also calculate the deviation from the local group affine motions to quantify the local plastic deformation. We find that both of them significantly depend on the particle size as well as the positions to other particles within the system. In contrast to bidisperse systems, the large particles in our simulations cause a new flow pattern for the relatively smaller particles. This flow pattern leads to more complicated ways of rearrangements that are the origins of the new behaviors we find. We further demonstrate a quantitative way to distinguish between “large” and “small” particles. Finally, we show how these results become increasingly important as the particle size distribution broadens. These findings are qualitatively different than previously found in bidisperse systems.
The other project is to apply the isomorph theory on the glassy simple system under simple shear. After cooling the Kob-Andersen binary Lennard-Jones system below the glass transition, we generate a so-called isomorph from the fluctuations of potential energy and virial in the NVT ensemble: a set of density, temperature pairs for which structure and dynamics are identical when expressed in appropriate reduced units. To access dynamical features, we shear the system using the SLLOD algorithm coupled with Lees-Edwards boundary conditions and study the statistics of stress fluctuations and the particle displacements transverse to the shearing direction in steady state. We find good collapse of the statistical data, showing that isomorph theory works well in this regime. The analysis of the distribution of stress fluctuations allows us to identify a clear signature of avalanche behavior in the form of an exponential tail on the negative side. This feature is also isomorph invariant. We then investigate further and turn our focus on the transient part of the stress and strain curve when the system yields. For the study here, we investigate a much larger density span over which the performance of various isomorph generating methods needs to be examined. Comparisons and comments on these methods are provided. We then shear the system along the identified isomorph. Here since the transient part depends on the thermal history, we shear the glassy samples generated by different cooing rates with different strain rates. Excellent collapsing quality again for steady state stress is verified. We notice however, that the peak stress at the transient part on the stress strain curve is not invariant, but decreases by a few percent for each ten percent increase in density, although the differences decrease with increasing density.
Table of Contents
This table of contents is under embargo until 26 May 2024
About this Dissertation
|Committee Chair / Thesis Advisor|
|File download under embargo until 26 May 2024||2022-03-24 14:17:42 -0400||File download under embargo until 26 May 2024|