No Arabic abstract
We study the dynamics of the solid to liquid transition for a model material made of elastic particles immersed in a viscous fluid. The interaction between particle surfaces includes their viscous lubrication, a sharp repulsion when they get closer than a tuned steric length and their elastic deflection induced by those two forces. We use Soft Dynamics to simulate the dynamics of this material when it experiences a step increase in the shear stress and a constant normal stress. We observe a long creep phase before a substantial flow eventually establishes. We find that the typical creep time relies on an internal relaxation process, namely the separation of two particles driven by the applied stress and resisted by the viscous friction. This mechanism should be relevant for granular pastes, living cells, emulsions and wet foams.
The critical dimension necessary for a flame to propagate in suspensions of fuel particles in oxidizer is studied analytically and numerically. Two types of models are considered: First, a continuum model, wherein the individual particulate sources are not resolved and the heat release is assumed spatially uniform, is solved via conventional finite difference techniques. Second, a discrete source model, wherein the heat diffusion from individual sources is modeled via superposition of the Greens function of each source, is employed to examine the influence of the random, discrete nature of the media. Heat transfer to cold, isothermal walls and to a layer of inert gas surrounding the reactive medium are considered as the loss mechanisms. Both cylindrical and rectangular (slab) geometries of the reactive medium are considered, and the flame speed is measured as a function of the diameter and thickness of the domains, respectively. In the continuum model with inert gas confinement, a universal scaling of critical diameter to critical thickness near 2:1 is found. In the discrete source model, as the time scale of heat release of the sources is made small compared to the interparticle diffusion time, the geometric scaling between cylinders and slabs exhibits values greater than 2:1. The ability of the flame in the discrete regime to propagate in thinner slabs than predicted by continuum scaling is attributed to the flame being able to exploit local fluctuations in concentration across the slab to sustain propagation. As the heat release time of the sources is increased, the discrete source model reverts back to results consistent with the continuum model. Implications of these results for experiments are discussed.
We investigated the transient relaxation of a Discontinuous Shear Thickening suspension of cornstarch in water. We performed 2 types of relaxation experiments starting from a steady shear in a parallel plate rheometer, followed by either stopping the top plate rotation and measuring the transient torque relaxation, or removing the torque on the plate and measuring the transient tool rotation. We found that at low weight fraction $phi_{eff}<58.8pm0.4%$, the suspensions exhibited a relaxation behavior consistent with a generalized Newtonian fluid. However, for larger weight fraction $58.8% < phi_{eff} < 61.0%$, near the liquid-solid transition $phi_c=61.0pm0.7%$, we found relaxation behaviors different from the generalized Newtonian model. The relaxation time in this range scales with the inverse of the critical shear rate at the onset of shear thickening. In this range the relaxation time was the same in both stress and rate controlled experiments, rather than the viscosity calculated from the relaxation time which is expected to be intrinsic material parameter in the generalized Newtonian model. The discrepancy between the measured relaxation times and the generalized Newtonian prediction was found to be up to $10^4$, and extrapolations diverge in the limit of $phi_c$ as the generalized Newtonian prediction approaches 0. At the highest weight fractions, the relaxation time scales were measured to be on the order of $sim 1$ s. The fact that this timescale is resolvable by the naked eye may be important to understanding some of the dynamic phenomenon commonly observed in these systems. We also showed that using the critical shear rate $dotgamma_c$ at the onset of shear thickening to characterize the effective weight fraction can more precisely characterize material properties near the critical point $phi_c$, allowing us to resolve this transition so close to $phi_c$.
An extremely broad and important class of phenomena in nature involves the settling and aggregation of matter under gravitation in fluid systems. Some examples include: sedimenting marine snow particles in lakes and oceans (central to carbon sequestration), dense microplastics in the oceans (which impact ocean ecology and the food chain), and even iron snow on Mercury (conjectured as its magnetic field source). These fluid systems all have stable density stratification, which is known to trap particulates through upper lightweight fluid coating the sinking particles, thus providing transient buoyancy. The current understanding of aggregation of such trapped matter involves collisions (due to Brownian motion, shear, and differential settling) and adhesion. Here, we observe and rationalize a new fundamental effective attractive mechanism by which particles suspended within stratification may self-assemble and form large aggregates without need for short range binding effects. This phenomenon arises through a complex interplay involving solute diffusion, impermeable boundaries, and aggregate geometry, which produces toroidal flows. We show that these toroidal flows yield attractive horizontal forces between particles. We observe that many particles demonstrate a collective motion revealing a system which self-assembles, appearing to solve jigsaw-like puzzles on its way to organizing into a disc-like shape, with the effective force increasing as the collective disc radius grows. Control experiments with two objects isolate the individual dynamics, which are quantitatively predicted through numerical integration of the underlying equations of motion. This new mechanism may be an important process in formation of marine snow aggregates and distribution of phytoplankton in lakes and oceans. Further, it potentially provides a new mechanism for general sorting and packing of layered material.
A multiple-relaxation-time discrete Boltzmann model (DBM) is proposed for multicomponent mixtures, where compressible, hydrodynamic, and thermodynamic nonequilibrium effects are taken into account. It allows the specific heat ratio and the Prandtl number to be adjustable, and is suitable for both low and high speed fluid flows. From the physical side, besides being consistent with the multicomponent Navier-Stokes equations, Ficks law and Stefan-Maxwell diffusion equation in the hydrodynamic limit, the DBM provides more kinetic information about the nonequilibrium effects. The physical capability of DBM to describe the nonequilibrium flows, beyond the Navier-Stokes representation, enables the study of the entropy production mechanism in complex flows, especially in multicomponent mixtures. Moreover, the current kinetic model is employed to investigate nonequilibrium behaviors of the compressible Kelvin-Helmholtz instability (KHI). It is found that, in the dynamic KHI process, the mixing degree and fluid flow are similar for cases with various thermal conductivity and initial temperature configurations. Physically, both heat conduction and temperature exert slight influences on the formation and evolution of the KHI.
For over 70 years it has been assumed that scalar wave propagation in (ensemble-averaged) random particulate materials can be characterised by a single effective wavenumber. Here, however, we show that there exist many effective wavenumbers, each contributing to the effective transmitted wave field. Most of these contributions rapidly attenuate away from boundaries, but they make a significant contribution to the reflected and total transmitted field beyond the low-frequency regime. In some cases at least two effective wavenumbers have the same order of attenuation. In these cases a single effective wavenumber does not accurately describe wave propagation even far away from boundaries. We develop an efficient method to calculate all of the contributions to the wave field for the scalar wave equation in two spatial dimensions, and then compare results with numerical finite-difference calculations. This new method is, to the authors knowledge, the first of its kind to give such accurate predictions across a broad frequency range and for general particle volume fractions.