Pulse propagation in particulate suspensions

The propagation of finite pulses of high-frequency in a suspension composed of a distribution of incompressible particles in an incompressible fluid is examined by using the technique of modulated simple-wave theory. The differential equation which governs the propagation of high-frequency pulses is derived and its consequences are examined. The properties of small amplitude high-frequency pulses are examined in detail.     

Some unsteady parallel flows of particulate suspensions

Analytical and numerical solutions are obtained for five unsteady parallel flows of particulate suspensions which exhibit boundary-layer characteristics. The governing equations are based on a continuum theory of particulate suspension behavior. An attempt is made to infer some of the pertinent features of two-phase and multiphase boundary-layer flows from the solutions.     

Development and validation of SuperDEM for non-spherical particulate systems using a superquadric particle method

This article presents the development and validation of the Superquadric Discrete Element Method(SuperDEM) for non-spherical particle simulation using a superqu...     

Development and validation of SuperDEM for non-spherical particulate systems using a superquadric particle method

This article presents the development and validation of the Superquadric Discrete Element Method (SuperDEM) for non-spherical particle simulation using a superquadric particle method in open-source CFD suite MFiX. A superquadric particle–particle contact algorithm with accelerating and stabilizing strategy was developed. A superquadric particle–arbitrary wall contact algorithm was developed, which enables the simulation in complex geometry. The solver was validated by comparing with experimental data generated in this study or available in the literature. Tests include cylinder contacting with a wall, static packing of M&M chocolate candies in a cylindrical container, static packing of cylinders in a cylindrical container, dynamic angle of repose of cylinders in a rotating drum, and discharging of chocolate candies from a hopper. Besides, MPI parallelization of the solver was implemented and the parallel performance of the solver using MPI was assessed through large-scale simulations of 1 million, 10 million, and 100 million particles on up to 6800 cores, which demonstrates that the SuperDEM solver has great potential for industrial-scale systems simulation.     

非球弹丸超高速撞击航天器防护结构数值模拟

采用AUTODYN软件对非球弹丸超高速正撞击航天器单防护屏防护结构进行了数值模拟，给出了2维及3维模拟的结果。研究了在相同质量和速度的条件下，不同形状弹丸长径比、撞击方向等对超高速撞击防护结构所产生碎片云特性及舱壁损伤尺寸的影响，并与球形弹丸撞击所产生的碎片云及舱壁损伤进行了比较。结果表明:弹丸的长径比越大，弹丸的穿孔能力越强；非球弹丸的撞击方向不同，所产生的碎片云形状、质量分布、破碎的程度和穿孔的能力是不同的。     

Self-consistent particle simulation of model-stabilized colloidal suspensions

Dynamics of model-stabilized colloidal suspensions were investigated by the self-consistent particle simulation method (SC), a new simulation algorithm that takes into account the interaction between the particles and suspending fluid. In this method, the fluid-particle interaction is introduced self-consistently by combining the finite element method (FEM) for fluid motion with Brownian dynamics (BD) for particle dynamics. To validate the reliability of the proposed algorithm, the shear dynamics of the stable particle suspensions were investigated. Relative viscosity and microstructure as a function of dimensionless shear rate at different volume fractions were in good agreement with previous observations. The robustness of the method was also verified through numerical convergence test. The effect of the fluid-particle interaction was well represented in simulations of two model problems, pressure-driven channel flow and rotating Couette flow. Plug-shaped velocity profile was observed in pressure-driven channel flow, which arised from shear thinning behavior of the stable suspension. In rotating Couette flow, shear banded nonlinear flow profile was observed. Although full hydrodynamic interaction (HI) was not rigorously taken into account, it successfully captured the macroscopic structure-induced flow field. It also takes advantage of the geometrical adaptability of FEM and computational efficiency of BD. We expect this newly developed simulation platform to be useful and efficient for probing the complex flow dynamics of particle systems as well as for practical applications in the complex flow of complex fluids.     

Structural instability in the sedimentation of particulate suspensions through viscoelastic fluids

A theory is developed to describe a structural instability that has been observed during the sedimentation of particulate suspensions through viscoelastic fluids. The theory is based on the assumption that the influence of hydrodynamic interactions in viscoelastic fluids, which tend to cause particles to aggregate, is in competition with hydrodynamic dispersion, which acts to maintain a homogeneous microstructure. In keeping with the experimental observations, it predicts that the suspension structure will stratify into vertical columns when a dimensionless stability parameter exceeds a critical value. The column-to-column separation, measured in particle radii, is predicted to be proportional to the square root of the ratio of the dimensionless dispersion coefficient to the product of the particle volume fraction and the Deborah number. The time for the formation of the columns is predicted to scale with the inverse of the average volume fraction. These predictions are in agreement with experimental data reported in the literature.     

Dynamic properties of shear thickening colloidal suspensions

The transient shear rheology (i.e., frequency and strain dependence) is compared to the steady rheology for a model colloidal dispersion through the shear thickening transition. Reversible shear thickening is observed and the transition stress compares well to theoretical predictions. Steady and transient shear thickening are observed to occur at the same value of the average stress. The critical strain for shear thickening is found to depend inversely on the frequency at fixed applied stress for low frequencies (high strains), but is limited to an apparent minimum critical strain at higher frequencies. This minimum critical strain is shown to be an artifact of slip. Lissajous plots illustrate the transition in material properties through the shear thickening transition, and the energy dissipated by a shear thickening suspension is analyzed as a function of strain amplitude.     

A stochastic model for filtration of particulate suspensions with incomplete pore plugging

A population balance model for particulate suspension transport with capture of particles by porous medium accounting for complete and incomplete plugging of pores by retained particles is derived. The model accounts for pore space accessibility, due to restriction on finite size particle movement through the overall pore space, and for particle flux reduction, due to transport of particles by the fraction of the overall flux. The novel feature of the model is the residual pore conductivity after the particle retention in the pore and the possibility of one pore to capture several particles. A closed system of governing stochastic equations determines the evolution of size distributions for suspended particles and pores. Its averaging results in the closed system of hydrodynamic equations accounting for permeability and porosity reduction due to plugging. The problem of deep bed filtration of a single particle size suspension through a single pore size medium where a pore can be completely plugged by two particles allows for an exact analytical solution. The phenomenological deep bed filtration model follows from the analytical solution.     

An extended unresolved CFD-DEM coupling method for simulation of fluid and non-spherical particles

This study develops an extended unresolved CFD-DEM coupling method for simulation of the fluid–solid flow with non-spherical particles. The limitation of fluid grid size is discussed, by simulating the settling of a cylinder in a Newtonian fluid based on the resolved and unresolved CFD-DEM coupling method. Then, the calculation of porosity and the fluid–particle relative velocity based on the particle shape enlargement method for simulation of non-spherical particles is proposed. The availability of the particle shape enlargement method for the simulation of non-spherical particles with different sphericity is discussed in this work, by comparing it with the results from the equivalent diameter enlargement method. The limitation of the equivalent diameter enlargement method for non-spherical particles is revealed from the simulation results. Several typical cases are employed to elaborate and verify the extended unresolved CFD-DEM method based on particle shape enlargement method, by presenting a good consistency with the experimental results. It proves that the extended unresolved CFD-DEM method is suitable for different CFD grid size ratios, and consolidates that it is a universal calculation method for CFD-DEM coupling simulation.     

Dynamic and rheological properties of soft biological cell suspensions

Quantifying dynamic and rheological properties of suspensions of soft biological particles such as vesicles, capsules, and red blood cells (RBCs) is fundamentally important in computational biology and biomedical engineering. In this review, recent studies on dynamic and rheological behavior of soft biological cell suspensions by computer simulations are presented, considering both unbounded and confined shear flow. Furthermore, the hemodynamic and hemorheological characteristics of RBCs in diseases such as malaria and sickle cell anemia are highlighted.     

Sphericities of non-spherical objects

Sphericity, as one of the most important shape parameter for non-spherical objects, is extensively applied in evaluating the porosity or packing density of particles. In this paper, the sphericities of common non-spherical objects are deduced and investigated. Maximum sphericities and optimum shapes of these objects are presented as well. A decreasing order of sphericity from sphere (1.0) to regular tetrahedron (0.671) for objects with constant sphericity is given. Similar trends are found in most sphericity–aspect ratio relationships, which exhibit single peak and the sphericity increases with the growth of aspect ratio before the peak point and decreases afterward. The peak loci of aspect ratio are all around 1.0 which makes the shape approaching to a sphere. The information in the paper could be useful as literature for general application.     

Computer simulation of the rheology of concentrated star polymer suspensions

We use particle-based computer simulations to study the rheology of suspensions of high-functionality star polymers with long entangled arms. Such particles have properties which are intermediate between those of soft colloidal particles and entangled polymer chains. In the simulations, each star polymer is coarse-grained to a single particle. In order to faithfully reproduce dynamical properties, it is very important to not only include time-averaged interactions (potentials of mean force) but to also account for transient interactions induced by entanglements between the arms of different star polymers. Using a model which has all these features, it is found that, for sufficiently high shear rates, the start-up shear stress displays an overshoot. With increasing concentration, the core interactions increasingly dominate the initial stress response, leading to a maximum in the stress overshoot at relatively low strain values (0.1 to 0.5). Transient forces start to dominate after this initial stage. In a simulated experiment in which the shear rate is suddenly stepped-down from a high to a lower value, the stress shows a clear undershoot, with the minimum stress again at a relatively low strain value (based on the new shear rate). Finally, it is shown that a stress plateau develops in the flow curve. This plateau is absent when the transient forces between the polymer stars are not taken into account.     

NUMERICAL SIMULATION OF ORIENTATION DISTRIBUTION FUNCTION OF CYLINDRICAL PARTICLE SUSPENSIONS

ＩｎｔｒｏｄｕｃｔｉｏｎＴｈｅｃｙｌｉｎｄｒｉｃａｌｐａｒｔｉｃｌｅｓｕｓｐｅｎｓｉｏｎｓｉｎｖｏｌｖｅｓｏｍｅｄｉｆｆｉｃｕｌｔｐｒｏｂｌｅｍｓｉｎｔｈｅｔｈｅｏｒｅｔｉｃａｌｒｅｓｅａｒｃｈｏｆｍｕｌｔｉｐｈａｓｅｆｌｏｗａｎｄｎｏｎ＿Ｎｅｗｔｏｎｉａｎｆｌｏｗ ,ａｎｄｈａｖｅｗｉｄｅａｐｐｌｉｃａｔｉｏｎｉｎｔｈｅｉｎｄｕｓｔｒｉａｌｅｎｇｉｎｅｅｒｉｎｇ .Ｆｏｒｅｘａｍｐｌｅ,ｄｕｒｉｎｇｔｈｅｐｒｏｃｅｓｓｉｎｇｏｆｓｏｍｅｍａｔｅｒｉａｌｓ,ｐｒｉｍｉｔｉｖｅｍａｔｅｒｉａｌｐｒｏｐｅ…     

Brownian dynamics simulation for suspensions of oblong-particles under shear flow

Computer simulations of the shear flow for the suspensions of oblong-particles were performed using nonequilibrium Brownian Dynamics (BD). The model particle is a rigid body made up of linearly connected spheres with the interparticle potential of a repulsive Lennard-Jones potential. The length-over-width ratios of the oblong-particles used in the present calculations are 5/3 and 3. In the concentrated suspensions high orientation is easily induced by shear at low shear rates. The systems of the oblong-particles exhibit the structural transition that causes the significant change in the rheological properties at high shear rates. Furthermore, the dependence of the length-over-width ratio of the particle is examined. Received: 16 June 1997 Accepted: 3 February 1998     

Direct-forcing immersed boundary lattice Boltzmann simulation of particle/fluid interactions for spherical and non-spherical particles

The lattice Boltzmann method (LBM) is a useful technique for simulating multiphase flows and modeling complex physics. Specifically, we use LBM combined with a direct-forcing (DF) immersed boundary (IB) method to simulate fluid–particle interactions in two-phase particulate flows. Two grids are used in the simulation: a fixed uniform Eulerian grid for the fluid phase and a Lagrangian grid that is attached to and moves with the immersed particles. Forces are calculated at each Lagrangian point. To exchange numerical information between the two grids, discrete delta functions are used. The resulting DF IB-LBM approach is then successfully applied to a variety of reference flows, namely the sedimentation of one and two circular particles in a vertical channel, the sedimentation of one or two spheres in an enclosure, and a neutrally buoyant prolate spheroid in a Couette flow. This last application proves that the developed approach can be used also for non-spherical particles. The three forcing schemes and the different factors affecting the simulation (added mass effect, corrected radius) are also discussed.     

Numerical simulation of thermo-mechanical fatigue properties for particulate reinforced composites

In this paper, a two dimensional Voronoi cell element, formulated with creep, thermal and plastic strain, is applied for the numerical simulation of thermo-mechanical fatigue behavior for particulate reinforced composites. The relation between mechanical fatigue phases and thermal fatigue phases influences the thermo-mechanical fatigue behavior and cyclic creep damage. The topological features of micro-structure in particulate reinforced composites, such as the orientation, depth-width ratio, distribution and volume fraction of inclusions, have a great influence on thermo-mechanical behavior. Some related conclusions are obtained by examples of numerical simulation.The project supported by the Special Funds for the National Major Fundamental Research Projects (2004CB619304), the National Natural Science Foundation of China (10276020 and 50371042), the Key Grant Project of Chinese Ministry of Education (0306)     

A boundary element method for the simulation of non-spherical bubbles and their interactions near a free surface

The basic principle and numerical technique for simulating two three-dimensional bubbles near a free surface are studied in detail by using boundary element method. The singularities of influence coefficient matrix are eliminated using coordinate transformation and so-called 4 π rule. The solid angle for the open surface is treated in direct method based on its definition. Several kinds of configurations for the bubbles and free surface have been investigated. The pressure contours during the evolution of bubbles are obtained in our model and can better illuminate the mechanism underlying the motions of bubbles and free surface. The bubble dynamics and their interactions have close relation with the standoff distances, buoyancy parameters and initial sizes of bubbles. Completely different bubble shapes, free surface motions, jetting patterns and pressure distributions under different parameters can be observed in our model, as demonstrated in our calculation results.     

非球形颗粒两相流的数值模拟研究进展

非球形颗粒两相流是多相流的重要研究方向之一, 常见于自然界及工业生产过程中. 不同于球形颗粒, 由于非球形颗粒形状的各向异性, 除了颗粒平动行为, 还需要考虑颗粒的转动与取向行为, 颗粒的取向与转动行为会影响颗粒所受的力和力矩. 为了准确模拟非球形颗粒的运动行为, 目前非球形颗粒两相流的数值模拟研究主要基于欧拉?拉格朗日的求解框架展开, 常见的非球形颗粒两相流数值模拟方法主要包括点颗粒法与全分辨颗粒法. 本文将对这两类方法进行介绍, 同时会全面介绍非球形颗粒两相流研究的基础理论模型, 并系统总结非球形颗粒在简单基本流和复杂湍流中的研究进展, 包括对于非球形颗粒在湍流中的取向与转动行为机理, 以及颗粒对湍流减阻调制作用的研究. 最后, 本文提出了非球形颗粒两相流研究存在的问题及未来研究方向.