A novel hybrid solid-like fluid-like (SLFL) method for the simulation of dry granular flows

This paper proposes a novel hybrid method to simulate the dry granular flow of materials over a wide range of inertial numbers that simultaneously covers the quasi-static and dense granular flow regimes. To overcome the lack of incremental objectivity whenever large deformations occur in solid-like regimes and to remove computational singularities in fluid-like regimes close to rest, the elastic–perfectly plastic theory based on the Drucker–Prager yield criterion is combined with the theory of dense granular flows. By implementing some new modifications at the boundaries and removing all ghost particles, smoothed particle hydrodynamics (SPH) is used as the framework for the method. A number of benchmark problems have been solved to show the capabilities of the new modified SPH method. Precise prediction of both location and pressure makes the modifications comparable with the previous works on SPH. Finally, the method is used to solve the classic 2D dry granular cliff collapse problem and to model dry granular material flow inside a rotary drum. The outcomes of the numerical simulation show good agreement with tabletop experiments and published results.     

Numerical simulations of granular free-surface flows using smoothed particle hydrodynamics

This paper presents a two-dimensional SPH model designed to simulate free-surface flows of dense granular materials. Smoothed particle hydrodynamics (SPH) is a mesh-free numerical method based on a Lagrangian discretization of the continuum mass and momentum conservation equations. The rheology of dense granular materials is modelled using a new local constitutive law recently proposed by Jop et al. (Nature, 2006). Of the viscoplastic class, this law is characterized by an apparent viscosity depending both on the local strain rate and local pressure. Validation test cases of the model in steady and unsteady configurations are presented. For steady cases (vertical chute flow and uniform free-surface layers on inclines), excellent agreement with analytical predictions is obtained. In the unsteady case, the simulations satisfactorily capture the dynamics of gravity-driven surges observed in experiments, including behaviours that are very specific to granular materials. Among the various parameters involved in the computations, the influence of SPH particle configuration within the flow and of the threshold viscosity used in the regularization of the constitutive yield criterion are particularly discussed.     

A three-invariant cap model with isotropic–kinematic hardening rule and associated plasticity for granular materials

In this paper, a three-invariant cap model is developed for the isotropic–kinematic hardening and associated plasticity of granular materials. The model is based on the concepts of elasticity and plasticity theories together with an associated flow rule and a work hardening law for plastic deformations of granulars. The hardening rule is defined by its decomposition into the isotropic and kinematic material functions. The constitutive elasto-plastic matrix and its components are derived by using the definition of yield surface, material functions and non-linear elastic behavior, as function of hardening parameters. The model assessment and procedure for determination of material parameters are described. Finally, the applicability of proposed plasticity model is demonstrated in numerical simulation of several triaxial and confining pressure tests on different granular materials, including: wheat, rape, synthetic granulate and sand.     

A CORRELATION EQUATION FOR CALCULATING INCLINED JET PENETRATION LENGTH IN A GAS-SOLID FLUIDIZED BED

Numerical simulation of gas-solid flow in a two-dimensional fluidized bed with an inclined jet was performed. The numerical model is based on the two-fluid model of gas and solids phase in which the solids constitutive equations are based on the kinetic theory of granular flow. The improved ICE algorithm, which can be used for both low and high-velocity fluid flow, were used to solve the model equations. The mechanism of jet formation was analyzed using both numerical simulations and experiments. The emergence and movement of gas bubbles were captured numerically and experimentally. The influences of jet velocity, nozzle diameter, nozzle inclination and jet position on jet penetration length were obtained. A semi-empirical expression was derived and the parameters were correlated from experimental data. The correlation equation, which can be easily used to obtain the inclined jet penetration length, was compared with our experimental data and published correlation equations.     

Evaluation of the effect of wall boundary conditions on numerical simulations of circulating fluidized beds

A computational fluid dynamics (CFD) modeling of the gas–solids two-phase flow in a circulating fluidized bed (CFB) riser is carried out. The Eularian–Eularian method with the kinetic theory of granular flow is used to solve the gas–solids two-phase flow in the CFB riser. The wall boundary condition of the riser is defined based on the Johnson and Jackson wall boundary theory (Johnson & Jackson, 1987) with specularity coefficient and particle–wall restitution coefficient. The numerical results show that these two coefficients in the wall boundary condition play a major role in the predicted solids lateral velocity, which affects the solid particle distribution in the CFB riser. And the effect of each of the two coefficients on the solids distribution also depends on the other one. The generality of the CFD model is further validated under different operating conditions of the CFB riser.     

GPU-accelerated smoothed particle hydrodynamics modeling of jet formation and penetration capability of shaped charges

The prediction of the penetration of three-dimensional (3D) shaped charge into steel plates is a challenging task. In this paper, the smoothed particle hydrodynamics (SPH) method is applied to simulate the jet formation generated by the shaped charge detonation and its damage to steel plates. The Jones–Wilkins–Lee (JWL) equation of state (EOS), Tillotson EOS, and elastic–perfectly plastic constitutive model were incorporated into SPH for the modeling of explosive detonation and dynamic behavior of metal material. The compute unified device architecture (CUDA) parallel programming interface has been employed in SPH to improve the computational efficiency of SPH. Firstly, the constitutive models and EOSs are validated by 3D TNT slab detonation and aluminum–aluminum (Al–Al) high velocity impact. Then the jet formation of the shaped charge detonation and its penetration into the steel plates are investigated using the graphics processing unit (GPU)-accelerated SPH methodology. The numerical results of these test cases are compared against the published experimental data or analytical result, which shows that the GPU-accelerated SPH methodology is capable of tackling the 3D shaped charge detonation and penetration involving millions of particles with high computational efficiency.     

A microstructural elastoplastic model for unsaturated granular materials

The homogenization technique is used to obtain an elastoplastic stress–strain relationship for dry, saturated and unsaturated granular materials. Deformation of a representative volume of material is generated by mobilizing particle contacts in all orientations. In this way, the stress–strain relationship can be derived as an average of the mobilization behavior of these local contact planes. The local behavior is assumed to follow a Hertz–Mindlin’s elastic law and a Mohr–Coulomb’s plastic law. For the non-saturated state, capillary forces at the grain contacts are added to the contact forces created by an external load. They are calculated as a function of the degree of saturation, depending on the grain size distribution and on the void ratio of the granular assembly. Numerical simulations show that the model is capable of reproducing the major trends of a partially saturated granular assembly under various stress and water content conditions. The model predictions are compared to experimental results on saturated and unsaturated samples of silty sands under undrained triaxial loading condition. This comparison shows that the model is able to account for the influence of capillary forces on the stress–strain response of the granular materials and therefore, to reproduce the overall mechanical behavior of unsaturated granular materials.     

An explicit formulation of a three-dimensional material damping model with transverse isotropy

A constitutive three-dimensional (3D) damping model is derived for transversely isotropic material symmetry, using the augmented Hooke's law [Intl. J. Solids Struct. 32 (1995) 2835] as a starting point. The proposed material model is tested numerically, via finite-element techniques, on a laminate structure built from stacked aluminium and Plexiglas plates. Effective 3D transversely isotropic material properties are given in terms of homogeneous material damping functions in connection with homogenised elastic laminate properties. Comparisons made between the results from the elastic (undamped) eigenvalue problem of the detailed (layerwise) model of the laminate and the effective 3D elastic model show that the homogenised model is reasonably accurate, in terms of predicted elastic eigenfrequencies for the first 20 modes. The dynamic homogenisation process, with damping included, is evaluated in terms of forced vibration response for the laminate structure, using effective transversely isotropic frequency dependent material properties. The dynamic 3D effective homogeneous material model is found to simulate very closely the detailed model in the studied frequency interval for the numerical test case.     

基于表面力模型的液滴冲击弹性基底SPH模拟

采用光滑粒子动力学SPH方法建立液滴冲击弹性基底的流固耦合数值模型,给出描述粘性流体和弹性固体运动的SPH离散方程和数值处理格式,引入人工耗散项来抑制标准SPH方法的数值震荡。为模拟液滴的表面张力效应,通过精确检测边界粒子,采用拉格朗日插值方法计算表面法向量和曲率,结合界面理论中的连续表面力CSF方法,建立了适用于自由表面液滴的表面力模型,方形液滴变形的模拟结果与拉普拉斯理论解吻合较好。随后,采用SPH流固耦合模型模拟1.0 mm直径水滴以不同速度(0.2 m/s～3.0 m/s)冲击两种薄板型基底,分析了基底弹性变形对液滴铺展、收缩以及回弹行为的影响。     

Incompressible SPH simulation of water entry of a free‐falling object

In this paper an incompressible smoothed particle hydrodynamics (Incom‐SPH) model is used to simulate the interactions between the free surface flow and a moving object. Incom‐SPH method is a two‐step semi‐implicit hydrodynamic formulation of the SPH algorithm and is capable of accurately treating the free surface deformations and impact forces during the solid–fluid interactions. For a free‐falling object, its motion is tracked by an additional Lagrangian algorithm based on Newton's law to couple with the Incom‐SPH program. The developed model is employed to investigate the water entry of a free‐falling wedge. The accuracy of the computations is validated by the good agreement in predicting the relevant hydrokinematic and hydrodynamic parameters. Finally, a numerical test is performed to study the influence of spatial resolution on the water entry features. The Incom‐SPH modeling coupled with the solid–fluid interaction algorithm could provide a promising computational tool to predict the slamming problems in coastal and offshore engineering. Copyright © 2008 John Wiley & Sons, Ltd.     

Investigation on effect of drag models on flow behavior of power-law fluid–solid two-phase flow in fluidized bed

In this study, a Eulerian-Eulerian two-fluid model combined with the kinetic theory of granular flow is adopted to simulate power-law fluid–solid two-phase flow in the fluidized bed. Two new power-law liquid–solid drag models are proposed based on the rheological equation of power-law fluid and pressure drop. One called model A is a modified drag model considering tortuosity of flow channel and ratio of the throat to pore, and the other called model B is a blending drag model combining drag coefficients of high and low particle concentrations. Predictions are compared with experimental data measured by Lali et al., where the computed porosities from model B are closer to the measured data than other models. Furthermore, the predicted pressure drop rises as liquid velocity increases, while it decreases with the increase of particle size. Simulation results indicate that the increases of consistency coefficient and flow behavior index lead to the decrease of drag coefficient, and particle concentration, granular temperature, granular pressure, and granular viscosity go down accordingly.     

Determination of material properties using nanoindentation and multiple indenter tips

Nanoindentation is a useful method to probe the material properties of a solid. Its effective use lies in interpreting the data collected from a nanoindentation experiment with an associated analytical/numerical solution of the corresponding problem configuration. In this paper, a parametric finite element study has been performed to develop a new procedure for extracting elastic–plastic properties of a material through nanoindentation experiments with a substantially improved accuracy for the elastic properties of a elastic–plastic solid. The procedure involves data collected through the use of two, different, nanoindenter tips. Non-dimensional functions were constructed for two different indenter geometries to show that test results from multiple indenters, when appropriately manipulated, deliver superior results, compared to using one indenter. The material was assumed to be an isotropic elastic–plastic solid with power law hardening. Friction between the indenter and the material was included in the cases studied. The ratio of yield strength to elastic modulus was assumed to be in the range 0.0005–0.02 and the hardening coefficient was assumed to be between 0 and 0.4. Poisson’s ratio was fixed at 0.3.     

Investigation of proper modeling of very dense granular flows in the recirculation system of CFBs

The aim of this paper is the development of new models and/or the improvement of existing numerical models, used for simulating granular flow in CFB (circulating fluidized bed) recirculation systems. Most recent models follow the TFM (two-fluid model) methodology, but they cannot effectively simulate the inter-particle friction forces in the recirculation system, because the respective stress tensor does not incorporate compressibility of flow due to change of effective particle density. As a consequence, the induced normal and shear stresses are not modeled appropriately during the flow of the granular phase in the CFB recirculation system. The failure of conventional models, such as that of von Mises/Coulomb, is mainly caused by false approximation of the yield criterion which is not applicable to the CFB recirculation system. The present work adopts an alternative yield function, used for the first time in TFM Eulerian modeling. The proposed model is based on the Pitman–Schaeffer–Gray–Stiles yield criterion. Both the temporal deformation of the solid granular phase and the repose angle that the granular phase forms are more accurately simulated by this model. The numerical results of the proposed model agree well with experimental data, implying that frictional forces are efficiently simulated by the new model.     

An elasto-plastic model for granular materials with microstructural consideration

In this paper, we have extended the granular mechanics approach to derive an elasto-plastic stress–strain relationship. The deformation of a representative volume of the material is generated by mobilizing particle contacts in all orientations. Thus, the stress–strain relationship can be derived as an average of the mobilization behavior of these local contact planes. The local behavior is assumed to follow a Hertz–Mindlin’s elastic law and a Mohr–Coulomb’s plastic law. Essential features such as continuous displacement field, inter-particle stiffness, and fabric tensor are discussed. The predictions of the derived stress–strain model are compared to experimental results for sand under both drained and undrained triaxial loading conditions. The comparisons demonstrate the ability of this model to reproduce accurately the overall mechanical behavior of granular media and to account for the influence of key parameters such as void ratio and mean stress. A part of this paper is devoted to the study of anisotropic specimens loaded in different directions, which shows the model capability of considering the influence of inherent anisotropy on the stress–strain response under a drained triaxial loading condition.     

A generalized constitutive elasticity law for GLPD micromorphic materials,with application to the problem of a spherical shell subjected to axisymmetric loading conditions

In this work we propose to replace the GLPD hypo-elasticity law by a more rigorous generalized Hooke's law based on classical material symmetry characterization assumptions. This law introduces in addition to the two well-known Lame's moduli, five constitutive constants. An analytical solution is derived for the problem of a spherical shell subjected to axisymmetric loading conditions to illustrate the potential of the proposed generalized Hooke's law.     

Two-fluid smoothed particle hydrodynamics simulation of submerged granular column collapse

In this paper, a two-fluid smoothed particle hydrodynamics (SPH) model, based on the mixture theory, is employed to investigate the complex interactions between the solid particles and the ambient water during the process of submerged granular column collapse. From the simulation, two regimes of the collapse, one being quick and the other being slow, are identified and the reasons of formation are analyzed. It is found that, a large internal friction angle of the granular phase, representing large drag force between solid particles, helps form the slow regime. Small hydraulic conductivity, representing large inter-phase drag force, also retards the collapse dramatically. Good agreements between our numerical results and other researchers’ numerical and experimental results are observed, which demonstrates the capability of the proposed two-fluid SPH approach in dealing with saturated water–soil mixture flows.     

Analytical modelling of thermal stresses in anisotropic composites

This paper represents a continuation of the author's previous work which deals with an analytical model of thermal stresses which originate during a cooling process of an anisotropic solid elastic continuum. This continuum consists of anisotropic spherical particles which are periodically distributed in an anisotropic infinite matrix. The infinite matrix is imaginarily divided into identical cubic cells with central particles. This multi-particle–matrix system represents a model system which is applicable to two-component materials of the precipitate–matrix type. The thermal stresses, which originate due to different thermal expansion coefficients of components of the model system, are determined within the cubic cell. The analytical modelling results from fundamental equations of continuum mechanics for solid elastic continuum (Cauchy's, compatibility and equilibrium equations, Hooke's law). This paper presents suitable mathematical procedures which are applied to the fundamental equations. These mathematical procedures lead to such final formulae for the thermal stresses which are relatively simple in comparison with the final formulae presented in the author's previous work which are extremely extensive. Using these new final formulae, the numerical determination of the thermal stresses in real two-component materials with anisotropic components is not time-consuming.     

VISCOPLASTIC SOLUTION TO FIELD AT STEADILY PROPAGATING CRACK TIP IN LINEAR- HARDENING MATERIALS

An elastic-viscoplastic constitutive model was adopted to analyze asymptotically the tip-field of moving crack in linear-hardening materials under plane strain condition. Under the assumption that the artificial viscosity coefficient was in inverse proportion to power law of the rate of effective plastic strain, it is obtained that stress and strain both possess power law singularity and the singularity exponent is uniquely determined by the power law exponent of the rate of effective plastic strain. Variations of zoning structure according to each material parameter were discussed by means of numerical computation for the tip-field of mode Ⅱ dynamic propagating crack, which show that the structure of crack tip field is dominated by hardening coefficient rather than viscosity coefficient. The secondary plastic zone can be ignored for weak hardening materials while the secondary plastic zone and the secondary elastic zone both have important influence on crack tip field for strong hardening materials. The dynamic solution approaches to the corresponding quasi-static solution when the crack moving speed goes to zero, and further approaches to the HR (Hui-Riedel) solution when the hardening coefficient is equal to zero.     

基于两相紊流模型的非平衡输沙研究

基于雷诺平均的水沙两相流方程, 建立了一个非平衡全沙输移二维数学模型. 模型考虑相间相对运动以及多颗粒之间的相互影响, 通过相间作用力进行两相耦合.和传统的单相流模型以及低浓度两相流模型相比, 该模型摆脱了依赖经验公式给定床面边界条件的局限性. 针对明渠净冲刷问题, 在合理给定水相和泥沙相边界条件的前提下计算了泥沙浓度分布的沿程变化, 并利用物理模型实验的结果和理论解验证了数学模型的正确性,同时也分析了明渠净冲刷问题中紊动扩散和重力沉降现象的特征.      

含液颗粒材料液固耦合分析的离散颗粒模型及特征线SPH法

对含液颗粒材料流固耦合分析建议了一个基于离散颗粒模型与特征线SPH法的显式拉格朗日-欧拉无网格方案。在已有的用以模拟固体颗粒集合体的离散颗粒模型[1]基础上,将颗粒间间隙内的流体模型化为连续介质,对其提出并推导了基于特征线的SPH法。数值例题显示了所建议方案在模拟颗粒材料与间隙流相互作用的能力和性能以及间隙流体对颗粒结构承载能力及变形的影响。