Removal of Particle Bridges From a Porous Material by Ultrasonic Irradiation

Particle bridge formation during the flow of a liquid with particles through a porous material is a fouling mechanism that can block the pores and, hence, decrease the permeability of the material. Ultrasonic irradiation of the material is a cleaning method that can restore the permeability. We make a numerical study of this cleaning method using the lattice-Boltzmann method. We start from a pore blocked by two spherical particles attached to the pore wall by colloidal adhesion forces, thus forming a particle bridge. Next we calculate the hydrodynamic force exerted by a high-frequency acoustic wave on the two particles. By comparing the hydrodynamic force and the adhesion force we investigate, whether the particle bridge will be removed by the ultrasonic irradiation. A sensitivity study is carried out to investigate the influence of some relevant parameters, such as the acoustic wave amplitude, the acoustic frequency, the fluid flow velocity and the ratio of particle diameter and pore diameter. An upscaling procedure is applied to translate the microscopic results for the removal of the particles at the pore level to the permeability improvement of the material at the macroscopic level. A comparison is made between numerical results and experimental data. The agreement is reasonable.     

Structural relaxation and avalanche dynamics of particle piles under vertical vibration

Granular matter can exhibit solid or liquid behavior, which contains complex physical mechanisms. In this work, we experimentally investigated the structural relaxation and avalanche dynamics of particle piles under vertical vibration. The influence of vibration parameters on the avalanche process was studied. The morphological features of avalanches were recorded and classified using high-speed camera. The effects of vibration parameters and particle properties on the relaxation mode are obtained. It is found that the evolution of particle pile height with time can be described by an exponential decay function. The relaxation rate and avalanche characteristics of four types of particles with different sizes are discussed. At the same acceleration level, for two larger particles, a smaller amplitude (A = 0.025 mm) leads to a faster relaxation rate, while for two smaller particles, a large amplitude (A = 0.500 mm) leads to a faster relaxation rate. The analogy powder surface tension is introduced to address the cohesion and flowability evolution of particles under vibration.     

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.     

Assessment of force models on finite-sized particles at finite Reynolds numbers

Finite-sized inertial spherical particles are fully-resolved with the immersed boundary projection method(IBPM) in the turbulent open-channel flow by direct numerical simulation(DNS). The accuracy of the particle surface force models is investigated in comparison with the total force obtained via the fully-resolved method. The results show that the steady-state resistance only performs well in the streamwise direction, while the fluid acceleration force, the added-mass force, and the shear-induced Saffman lift can effectively compensate for the large-amplitude and high-frequency characteristics of the particle surface forces, especially for the wall-normal and spanwise components. The modified steady-state resistance with the correction effects of the acceleration and the fluid shear can better represent the overall forces imposed on the particles, and it is a preferable choice of the surface force model in the Lagrangian point-particle method.     

Particle image velocimetry investigation of a finite amplitude pressure wave

Particle image velocimetry is used to study the motion of gas within a duct subject to the passage of a finite amplitude pressure wave. The wave is representative of the pressure waves found in the exhaust systems of internal combustion engines. Gas particles are accelerated from stationary to 150 m/s and then back to stationary in 8 ms. It is demonstrated that gas particles at the head of the wave travel at the same velocity across the duct cross section at a given point in time. Towards the tail of the wave viscous effects are plainly evident causing the flow profile to tend towards parabolic. However, the instantaneous mean particle velocity across the section is shown to match well with the velocity calculated from a corresponding measured pressure history using 1D gas dynamic theory. The measured pressure history at a point in the duct was acquired using a high speed pressure transducer of the type typically used for engine research in intake and exhaust systems. It is demonstrated that these are unable to follow the rapid changes in pressure accurately and that they are prone to resonate under certain circumstances.     

Research on the effect of cylinder particles on the turbulent properties in particulate flows

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

Effect of poly(ethylene oxide) on the rheological behavior of silica suspensions

The steady-shear viscosity, dynamic viscoelasticity, and sedimentation behavior were measured for silica suspensions dispersed in aqueous solutions of poly(ethylene oxide) (PEO). For suspensions prepared with polymer solutions in which the transient network is developed by entanglements, the viscosity at a given shear rate decreases, shows a minimum, and then increases with increasing particle concentration. Because the suspensions are sterically stabilized under the conditions where the particle surfaces are fully covered with by a thick layer of adsorbed polymer, the viscosity decrease can be attributed to the reduction of network density in solution. But under the low coverage conditions, the particles are flocculated by bridging and this leads to a viscosity increase with shear-thinning profiles. The polymer chains with high molecular weights form flexible bridges between particles. The stress-dependent curve of storage modulus measured by a stress amplitude sweep shows an increase prior to a drastic drop due to structural breakdown. The increase in elastic responses may arise from the restoring forces of extended bridges with high deformability. The effect of PEO on the rheological behavior of silica suspensions can be explained by a combination of concentration reduction of polymer in solution and flocculation by bridging.     

Two-phase couette flow

A study is made of plane laminar Couette flow, in which foreign particles are injected through the upper boundary. The effect of the particles on friction and heat transfer is analyzed on the basis of the equations of two-fluid theory. A two-phase boundary layer on a plate has been considered in [1, 2] with the effect of the particles on the gas flow field neglected. A solution has been obtained in [3] for a laminar boundary layer on a plate with allowance for the dynamic and thermal effects of the particles on the gas parameters. There are also solutions for the case of the impulsive motion of a plate in a two-phase medium [4–6], and local rotation of the particles is taken into account in [5, 6]. The simplest model accounting for the effect of the particles on friction and heat transfer for the general case, when the particles are not in equilibrium with the gas at the outer edge of the boundary layer, is Couette flow. This type of flow with particle injection and a fixed surface has been considered in [7] under the assumptions of constant gas viscosity and the simplest drag and heat-transfer law. A solution for an accelerated Couette flow without particle injection and with a wall has been obtained in [6]. In the present paper fairly general assumptions are used to obtain a numerical solution of the problem of two-phase Couette flow with particle injection, and simple formulas useful for estimating the effect of the particles on friction and heat transfer are also obtained.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 3, pp. 42–46, May–June, 1976.     

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.     

水平摆振下环形巴西果分离效应模拟研究

为了研究水平摆振下混合颗粒体系的分离行为及其机理,利用离散单元法进行数值模拟,对二元混合颗粒振动下的运动行为进行了研究,观察到圆筒形容器内一种新的分离现象,即"环形巴西果"分离。通过对比不同振动参数、不同区域内混合颗粒体系的运动特性,分析了颗粒的整体流态、能量、离心力作用以及颗粒间的力链结构强度,对混合颗粒的中心偏析、边界堆积以及整体分离行为,揭示了"环形巴西果"分离构型的内在机理。结果表明:频率为1Hz、振幅为160°时,混合颗粒出现"环形巴西果"分离现象;频率为3Hz、振幅为40°时,混合颗粒出现"巴西果"分离现象,且通过调整振幅和频率可以控制中心区域大颗粒偏析程度,实现"环形巴西果"分离构型与"巴西果"分离构型之间的转换。     

柔性微粒介电泳分离过程的多尺度模拟

介电泳分离是一种高效的微细颗粒分离技术,利用非均匀电场极化并操纵分离微流道中的颗粒. 柔性微粒在介电泳分离过程中同时受多种物理场、多相流和微粒变形等复杂因素的影响,仅用单一的计算方法对其进行模拟存在一定的难度,本文采用有限单元——格子玻尔兹曼耦合计算的方法处理这一难题.介观尺度的格子玻尔兹曼方法将流体看成由大量微小粒子组成,在离散格子上求解玻尔兹曼输运方程,易于处理多相流及大变形问题,特别适合模拟柔性颗粒在介电泳分离过程中的变形情况.另一方面,介电泳分离过程的模拟需求解流体、电场和微粒运动方程,计算量相当庞大,通过有限单元法求解介电泳力,提高计算效率.利用这种多尺度耦合计算方法,对一款现有的介电泳芯片分离过程进行了模拟.分析了微粒在电场作用下产生的介电泳力,揭示了介电泳力与电场变化率等因素之间的关系.对微粒运动轨迹及其变形的情况进行了研究,发现微粒的变形主要与流体剪切作用有关.这种多尺度耦合计算方法,为复杂微流体的计算提供了一种有效的解决方案.      

A multiphase shock tube for shock wave interactions with dense particle fields

Currently there is a substantial lack of data for interactions of shock waves with particle fields having volume fractions residing between the dilute and granular regimes. To close this gap, a novel multiphase shock tube has been constructed to drive a planar shock wave into a dense gas–solid field of particles. A nearly spatially isotropic field of particles is generated in the test section by a gravity-fed method that results in a spanwise curtain of spherical 100-micron particles having a volume fraction of about 20%. Interactions with incident shock Mach numbers of 1.66, 1.92, and 2.02 are reported. High-speed schlieren imaging simultaneous with high-frequency wall pressure measurements are used to reveal the complex wave structure associated with the interaction. Following incident shock impingement, transmitted and reflected shocks are observed, which lead to differences in particle drag across the streamwise dimension of the curtain. Shortly thereafter, the particle field begins to propagate downstream and spread. For all three Mach numbers tested, the energy and momentum fluxes in the induced flow far downstream are reduced about 30–40% by the presence of the particle field.     

Experimental investigation of a buckled beam under high-frequency excitation

In this paper, we investigate theoretically and experimentally dynamics of a buckled beam under high-frequency excitation. It is theoretically predicted from linear analysis that the high-frequency excitation shifts the pitchfork bifurcation point and increases the buckling force. The shifting amount increases as the excitation amplitude or frequency increases. Namely, under the compressive force exceeding the buckling one, high-frequency excitation can stabilize the beam to the straight position. Some experiments are performed to investigate effects of the high-frequency excitation on the buckled beam. The dependency of the buckling force on the amounts of excitation amplitude and frequency is compared with theoretical results. The transient state is observed in which the beam is recovered from the buckled position to the straight position due to the excitation. Furthermore, the bifurcation diagrams are measured in the cases with and without high-frequency excitation. It is experimentally clarified that the high-frequency excitation changes the nonlinear property of the bifurcation from supercritical pitchfork bifurcation to subcritical pitchfork bifurcation and then the stable steady state of the beam exhibits hysteresis as the compressive force is reversed. This work was partially supported by the Japanese Ministry of Education, Culture, Sports, Science, and Technology, under Grants-in-Aid for Scientific Research 16560377.     

Approximate solutions for power-law fluid flow past a particle at low Reynolds numbers

A simple perturbation approximation is proposed for describing flow behaviour of particles immersed in a uniform flow and an extensional flow of power-law fluids. The present solution for particles in a uniform flow field is in good agreement with the numerical solution in the literature. Theoretical predictions indicate that the effect of pseudoplasticity on flow around particles in an extensional flow field is small compared with that for particles in a uniform flow field.From the viewpoint of perturbation techniques, existing analytical solutions based on linearization of the equations of motion for particle in a power-law fluid are re-examined. Mass transfer to a power-law fluid from a particle is also discussed.     

Determination of particle charge and size from the propagation of ultrasound in suspensions

The propagation of small perturbations in raulticomponent disperse media consisting of an uncharged dispersion fluid, positive and negative ions and charged particles or droplets of another fluid is investigated. When weak waves pass through emulsions and suspensions, because of the difference in the velocities of the ions and charged particles a non-uniform distribution of electric potential develops in the medium [1–3]. Expressions relating the amplitude of the electric potential and the amplitude of the fluid velocity in the wave, the particle charge and the parameters characterizing the medium are derived. Relations are obtained for the phase shift between the values of the electric potential and the fluid velocity. It is proposed to use the expressions obtained, which describe the propagation of ultrasound, for the experimental determination of the particle charge and other parameters of the disperse medium, in particular, the particle size.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 1, pp. 122–128, January–February, 1988.     

Numerical simulation of polygonal particles moving in incompressible viscous fluids

A two-dimensional coupled lattice Boltzmann immersed boundary discrete element method is introduced for the simulation of polygonal particles moving in incompressible viscous fluids. A collision model of polygonal particles is used in the discrete element method. Instead of a collision model of circular particles, the collision model used in our method can deal with particles of more complex shape and efficiently simulate the effects of shape on particle–particle and particle–wall interactions. For two particles falling under gravity, because of the edges and corners, different collision patterns for circular and polygonal particles are found in our simulations. The complex vortexes generated near the corners of polygonal particles affect the flow field and lead to a difference in particle motions between circular and polygonal particles. For multiple particles falling under gravity, the polygonal particles easily become stuck owing to their corners and edges, while circular particles slip along contact areas. The present method provides an efficient approach for understanding the effects of particle shape on the dynamics of non-circular particles in fluids.     

Development and validation of a reduced order history force model

The history force model accounts for temporal development in fluid gradients in the viscous region surrounding a particle in point particle methods. The calculation of the history force typically requires storing and using relative velocity information during the life time of the particle. For a large number of particles integrated over large times, history force calculation can become prohibitively expensive. The current work presents a new modeling approach to calculate the history force in which a decay function is applied to a stored cumulative value of the history force. The proposed formulation is equivalent to applying the same function obtained from a constant acceleration assumption to a running average of the acceleration within the memory time of the particle. The new force model is validated with experimental measurements of settling spheres at Reynolds numbers ranging from around one to a few hundreds and at density ratios from 1.2 to about 9.32. More validation work was carried-out with experimental measurements of oscillating spheres at different frequencies and amplitudes, as well as bouncing spheres at different Reynolds numbers and density ratios. The model shows very good agreement with the experiments of settling spheres and reasonable/good agreement with oscillating and bouncing sphere experiments. The proposed model significantly reduces the computational resources required to calculate the history force especially when large number of particles need to be integrated over long times.     

DEM simulation of influence of parameters on energy and degree of mixing

Powder mixing is being modeled using a simulation based on Newtonian mechanics. Variables under consideration include particle friction and the amplitude, frequency, and direction of shaking. Trajectories for each particle were recorded, and a mixing degree was calculated for each simulation, for which the average energy transferred into the powder system was recorded and compared to the resulting mixing degree. Mixing of particles originally located near the bottom was studied separately, as was the mixing of particles near the surface. This study shows that choosing the proper mixing parameters not only enhances the final result of mixing, but also yields good results with less strain on the material mixed and on the mixing device.     

The elastic response of a cohesive aggregate—a discrete element model with coupled particle interaction

A model is presented for the deformation of a cohesive aggregate of elastic particles that incorporates two important effects of large-sized inter-particle junctions. A finite element model is used to derive a particle response rule, for both normal and tangential relative deformations between pairs of particles. This model agrees with the Hertzian contact theory for small junctions, and is valid for junctions as large as half the nominal particle size. Further, the aggregate model uses elastic superposition to account for the coupled force–displacement response due to the simultaneous displacement of all of the neighbors of each particle in the aggregate. A particle stiffness matrix is developed, relating the forces at each junction to the three displacement degrees of freedom at all of the neighboring-particle junctions. The particle response satisfies force and moment equilibrium, so that the model is properly posed to allow for rigid rotation of the particle without introducing rotational degrees of freedom. A computer-simulated sintering algorithm is used to generate a random particle packing, and the stiffness matrix is derived for each particle. The effective elastic response is then estimated using a mean field or affine displacement calculation, and is also found exactly by a discrete element model, solving for the equilibrium response of the aggregate to uniform-strain boundary conditions. Both the estimate and the exact solution compare favorably with experimental data for the bulk modulus of sintered alumina, whereas Hertzian contact-based models underestimate the modulus significantly. Poisson's ratio is, however, accurately determined only by the full equilibrium discrete element solution, and shown to depend significantly on whether or not rigid particle rotation is permitted in the model. Moreover, this discrete element model is sufficiently robust, so it can be applied to problems involving non-homogeneous deformations in such cohesive aggregates.     

孤立波与半无限复合材料体耦合作用研究

基于一维颗粒链中产生的高度非线性孤立波,研究孤立波与半无限复合材料体的耦合作用。根据赫兹定律推导了一维颗粒链中颗粒间相互作用的运动微分方程,建立了颗粒链与半无限复合材料体的接触模型。对于颗粒与复合材料的接触,采用已有文献中修正后的赫兹定律,研究了高度非线性孤立波与半无限复合材料体的耦合力学作用机理,推导了颗粒链与半无限复合材料体的相互耦合运动微分方程组,通过数值计算,得到了各颗粒的内力、速度、位移曲线。分析了材料属性对回弹孤立波出现的时间、幅值的影响。结果表明:随着纤维方向弹性模量的增大,次级回弹波出现的时间和波幅都逐渐增大,随着垂直纤维方向弹性模量的增大,次级回弹波出现的时间先减小后增大,次级回弹波的幅值逐渐减小直至消失。