Contact Changes of Sheared Systems: Scaling, Correlations, and Mechanisms

We probe the onset and effect of contact changes in 2D soft harmonic particle packings which are sheared quasistatically under controlled strain. First, we show that in the majority of cases, the first contact changes correspond to the creation or breaking of contacts on a single particle, with contact breaking overwhelmingly likely for low pressures and/or small systems, and contact making and breaking equally likely for large pressures and in the thermodynamic limit. The statistics of the corresponding strains are near-Poissonian. The mean characteristic strains exhibit scaling with the number of particles N and pressure P, and reveal the existence of finite size effects akin to those seen for linear response quantities. Second, we show that linear response accurately predicts the strains of the first contact changes, which allows us to study the scaling of the characteristic strains of making and breaking contacts separately. Both of these show finite size scaling, and we formulate scaling arguments that are consistent with the observed behavior. Third, we probe the effect of the first contact change on the shear modulus G, and show in detail how the variation of G remains smooth and bounded in the large system size limit: even though contact changes occur then at vanishingly small strains, their cumulative effect, even at a fixed value of the strain, are limited, so that effectively, linear response remains well-defined. Fourth, we explore multiple contact changes under shear, and find strong and surprising correlations between alternating making and breaking events. Fifth, we show that by making a link with extremal statistics, our data is consistent with a very slow crossover to self averaging with system size, so that the thermodynamic limit is reached much more slowly than expected based on finite size scaling of elastic quantities or contact breaking strains

Direct Numerical Simulation of Complex Gas-Liquid-Solid Flows using a Combined\ud Immersed Boundary (IB) and Volume Of Fluid (VOF) Approach

In this paper a simulation model is presented for the Direct Numerical Simulation (DNS) of complex multi-fluid flows in which simultaneously (moving) deformable (drops or bubbles) and non-deformable (moving) elements (particles) are present, possibly with the additional presence of free surfaces. Our model combines the VOF model developed by van Sint Annaland et al. (2005) and the Immersed Boundary (IB) model The Volume of Fluid (VOF) part features i) an interface reconstruction technique based on piecewise linear interface representation ii) a three-dimensional version of the CSF model of Brackbill et al. (1992). The Immersed Boundary (IB) part incorporates both particle-fluid and particle-particle interaction via a Direct Forcing Method (DFM) and a hard sphere Discrete Particle (DP) approach. In our model a fixed (Eulerian) grid is utilized to solve the Navier-Stokes equations for the entire computational domain. The no-slip condition at the surface of the moving particles is enforced via a momentum source term which only acts in the vicinity of the particle surface. Specifically Lagrangian force points are used which are distributed evenly over the surface of the particle. Dissipative particle-particle and/or particle-wall collisions are accounted via a hard sphere DP approach using a three-parameter particle-particle interaction model accounting for normal and tangential restitution and tangential friction. The capabilities of the hybrid VOF-IB model are demonstrated with a number of examples in which complex topological changes in the interface are encountered

Numerical Simulation of Heat Transport in Dispersed Gas-Liquid Two-Phase Flow using a Front Tracking Approach

In this paper a simulation model is presented for the Direct Numerical Simulation (DNS) of heat transport in dispersed gas-liquid two-phase flow using the Front Tracking (FT) approach. Our model extends the FT model developed by van Sint Annaland et al. (2006) to non-isothermal conditions. In FT an unstructured dynamic mesh is used to represent and track the interface explicitly by a number of interconnected marker points. The Lagrangian representation of the interface avoids the necessity to reconstruct the interface from the local distribution of the fractions of the phases and, moreover, allows a direct and accurate calculation of the surface tension force circumventing the (problematic) computation of the interface curvature. The extended model is applied to predict the heat exchange rate between the liquid and a hot wall kept at a fixed temperature. It is found that the wall-to-liquid heat transfer coefficient exhibits a maximum in the vicinity of the bubble that can be attributed to the locally decreased thickness of the thermal boundary layer

Direct Numerical Simulation of Complex Multi-Fluid Flows Using a Combined Volume of Fluid and Immersed Boundary Method

In this paper a simulation model is presented for the Direct Numerical Simulation (DNS) of complex multi-fluid flows in which simultaneously (moving) deformable (drops or bubbles) and non-deformable (moving) elements (particles) are present, possibly with the additional presence of free surfaces. Our model combines the VOF model developed by van Sint Annaland et al. (2005) and the Immersed Boundary (IB) model developed by van der Hoef et al. (2006). The Volume of Fluid (VOF) part features i) an interface reconstruction technique based on piecewise linear interface representation ii) a three-dimensional version of the CSF model of Brackbill et al. (1992). The Immersed Boundary (IB) part incorporates both particle-fluid and particle-particle interaction via a Direct Forcing Method (DFM) and a hard sphere Discrete Particle (DP) approach. In our model a fixed (Eulerian) grid is utilized to solve the Navier-Stokes equations for the entire computational domain. The no-slip condition at the surface of the moving particles is enforced via a momentum source term which only acts in the vicinity of the particle surface. For the enforcement of the no-slip condition Lagrangian force points are used which are distributed evenly over the surface of the particle. Dissipative particle-particle and/or particle-wall collisions are accounted via a hard sphere DP approach (Hoomans et al., 1996) using a three-parameter particle-particle interaction model accounting for normal and tangential restitution and tangential friction. The capabilities of the hybrid VOF-IB model are demonstrated with a number of examples in which complex topological changes in the interface are encountered

Gas-solid two-phase turbulent flow in a circulating fluidized bed riser: an\ud experimental and numerical study

Hydrodynamics of gas-particle two-phase turbulent flow in a circulating fluidized bed riser is studied experimentally by Particle Image Velocimetry (PIV) and numerically with the use of a 3D discrete hard sphere particle model (DPM). Mean particle velocities and RMS velocities are obtained and the influence of turbulence on the flow is investigated. The experimental data are analyzed and compared with the numerical results showing a reasonable agreement

The essence of fertilization: oocyte meets sperm

The problem of reduced fertility in high yielding dairy cattle is a very complicated one, and the relationship between various measures of fertility and level of milk production remains controversial. In this brief review the essence of the problem is considered: what is the oocyte's and the sperm's contribution, and what is the importance of the resulting embryo in the declining fertility of the Holstein Friesian cow

A discrete particle simulation study on the influence of restitution coefficient on spout fluidized bed dynamics

In this paper the influence of the normal restitution coefficient on the bed dynamics\ud during different flow regimes was studied using the discrete element model. The three flow\ud regimes comprise the intermediate/spout-fluidization, spouting-with-aeration regime and the jetin-\ud fluidized-bed. It is shown that with increase of the restitution coefficient the average bed\ud height decreases for all flow regime cases. To study the influence of the granule impact velocity\ud and the liquid film on the wall surface the free-fall experiments in the velocity range of 0.5-\ud 4.5 m/s were performed with spherical γ-Al2O3 granules. During the free-fall tests, the impact\ud angle from 0° to 80° and the viscosity of thin liquid film in the range of 1-300 mPa·s were\ud varied

Direct numerical simulation of heat transport in dispersed gas-liquid two-phase flow using a front tracking approach

In this paper a simulation model is presented for the Direct Numerical Simulation (DNS) of heat transport in dispersed gas-liquid two-phase flow using the Front Tracking (FT) approach. Our model extends the FT model developed by van Sint Annaland et al. (2006) to non-isothermal conditions. In FT an unstructured dynamic mesh is used to represent and track the interface explicitly by a number of interconnected marker points. The Lagrangian representation of the interface avoids the necessity to reconstruct the interface from the local distribution of the fractions of the phases and, moreover, allows a direct and accurate calculation of the surface tension force circumventing the (problematic) computation of the interface curvature. The extended model is applied to predict the heat exchange rate between the liquid and a hot wall kept at a fixed temperature. It is found that the wall-to-liquid heat transfer coefficient exhibits a maximum in the vicinity of the bubble that can be attributed to the locally decreased thickness of the thermal boundary layer