Smoothed particle hydrodynamics for coarse-grained modeling of rapid granular flow

We simulated rapid flow in transient plane Couette flows of granular particles using the smoothed particle hydrodynamics(SPH) solutions of a set of continuum equations.This simulation was performed to test the viability of SPH in solving the equations for the solid phase of the two-fluid model associated with fluidization.We found that SPH requires the handling of fewer particles in simulating the collective behavior of rapid granular flow,thereby bolstering expectations of solving the equations for the solid phase in the two-fluid modeling of fluidization.Further work is needed to investigate the effect of terms describing pressure and viscous stress of solids on stability in simulations.     

SPH算法在长杆弹侵彻多层间隔靶中的应用

基于Autodyn软件,应用SPH算法对长杆弹高速侵彻多层间隔金属靶板的过程进行数值模拟,分析了弹体侵彻靶板后的剩余速度、剩余动能、以及形成的破片及其飞散情况,同时模拟初始速度及弹体质量对侵彻过程的影响.通过与实验结果的比较,表明数值模拟的结果是合理的,能够准确模拟长杆弹对多层间隔靶的高速侵彻过程,对预测弹体侵彻件能以...     

Smoothed particle hydrodynamics formulation for penetrating impacts on ballistic gelatine

A particle method is applied to the investigation of impact biomechanics in the case of penetrating ballistic. A three dimensional model is proposed using the Smoothed Particle Hydrodynamics (SPH) method combined with Finite Elements (FE) method. The problem consists in the violent impact of a steel sphere on soft tissues, simulated by 20% ballistic gelatine (BG) material which is considered as a very interesting human tissue surrogate. Comparisons with experimental data are established to validate the proposed model. The results, in terms of penetrating curves, show very promising results. The use of particle methods appears to be an interesting way to model high speed loading, especially penetrating ballistic impact whose classical FE modelling can bring some important limitations in terms of mesh and element distortions.     

三维溃坝流的光滑粒子动力学方法模拟

采用光滑粒子动力学SPH(Smoothed Particle Hydrodynamics)方法对三维溃坝流问题进行了数值模拟。为了逼真地模拟出坝内水体与壁面间相互作用而产生的水花飞溅、融合以及近壁面流动等现象,加入了混合长度形式的湍流模型。为了有效地防止粒子穿透固壁,提出了一种新型的适合三维数值模拟的固壁边界处理方法。应用SPH方法对三维溃坝流进行了数值模拟,并分别考虑了未添加障碍物和添加圆柱障碍物两种情形。计算结果表明,改进SPH方法能够精细地捕捉溃坝流在不同时刻的自由液面,并获得稳定而精确的数值结果。     

Incompressible smoothed particle hydrodynamics simulation of multifluid flows

This paper presents an incompressible SPH (ISPH) technique to simulate multifluid flows. The SPH method is a mesh‐free particle modeling approach that can treat free surfaces and multi‐interfaces in a simple and efficient manner. The ISPH method employs an incompressible hydrodynamic formulation to solve the fluid pressure that ensures a stable pressure field. Two multifluid ISPH models are proposed following different interface treatments: the coupled ISPH model does not distinguish the different fluid phases and applies the standard ISPH technique across the interface, whereas the decoupled ISPH model first treats each fluid phase separately and then couples the different phases by applying pressure and shear stress continuities across the interface. The two proposed models were used to investigate a gravity underflow with a low density ratio in a Generalized Reservoir Hydrodynamics (GRH) flume and a horizontal lock exchange flow with a high density ratio. Comparisons with data and relevant numerical error analysis indicated that the decoupled model performed well in cases of both low and high density ratios because of the accurate treatment of interface boundaries. Copyright © 2011 John Wiley & Sons, Ltd.     

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.     

Transient non-Newtonian flow of a suspension with a compressible particle phase

Equations governing transient two-phase fluid-particle laminar flow over an infinite porous flat plate are developed. Both phases are assumed to behave as non-Newtonian power-law fluids. The mathematical model accounts for particle-phase viscous and diffusive effects. The particles are assumed spherical in shape and having a non-uniform density distribution. The resulting governing equations are nondimensionalized and solved numerically subject to appropriate initial and boundary conditions using an iterative, implicit, tri-diagonal finite-difference method. Graphical results for the displacement thicknesses and the skin-friction coefficients for both the fluid and particle phases are presented and discussed to illustrate special trends of the solutions.     

Numerical simulation of landslide impulsive waves by incompressible smoothed particle hydrodynamics

An incompressible‐smoothed particle hydrodynamics (I‐SPH) formulation is presented to simulate impulsive waves generated by landslides. The governing equations, Navier–Stokes equations, are solved in a Lagrangian form using a two‐step fractional method. Landslides in this paper are simulated by a submerged mass sliding along an inclined plane. During sliding, both rigid and deformable landslides mass are considered. The present numerical method is examined for a rigid wedge sliding into water along an inclined plane. In addition solitary wave generated by a heavy box falling inside water, known as Scott Russell wave generator, which is an example for simulating falling rock avalanche into artificial and natural reservoirs, is simulated and compared with experimental results. The numerical model is also validated for gravel mass sliding along an inclined plane. The sliding mass approximately behaves like a non‐Newtonian fluid. A rheological model, implemented as a combination of the Bingham and the general Cross models, is utilized for simulation of the landslide behaviour. In order to match the experimental data with the computed wave profiles generated by deformable landslides, parameters of the rheological model are adjusted and the numerical model results effectively match the experimental results. The results prove the efficiency and applicability of the I‐SPH method for simulation of these kinds of complex free surface problems. Copyright © 2007 John Wiley & Sons, Ltd.     

Variational principles in hydrodynamics of a non-Newtonian fluid

I.IntroductionTileresearchworkonvariationalprincipleinhydrodynamicsstartedasearlyas1940's.Referencesl:~6]aretheoneswhichmakemoreimportantcolltributions.l-lowevcr,mostoftheseworksstudied11on-viscousflow,putemphasisontheproblemsofexternalllcld,andworkedontilebasisofBernoulli'sequation.Prot'essorChienWeizang11I.slesttlblishcdvariationalprincipleforhydrodynalllicpl-oblemofviscousfluid,thatistheprillciplcofmaxilnLllllpowerlosses,alldgeneralizedvariational(stationary)pl.inciplesoiltilebiLsisofNav…     

Smoothed particle hydrodynamics techniques for the solution of kinetic theory problems: Part 2. The effect of flow perturbations on the simple shear behavior of LCPs

The effect of slight perturbations to simple shear flow of liquid-crystalline polymers (LCPs) is explored by using the SPH technique to solve the unapproximated orientation distribution function equation arising from the Doi Theory. First, the case of simple shear flow is outlined, and it is shown that skewed distributions play an important role in the transition from periodic to steady behavior as the shear rate is increased. Next, we consider perturbations to flows that are slightly more extensional than simple shear, parametrized by the flow type parameter α. They are shown to eliminate all periodic director behavior (tumbling and wagging), even when the relative increment in flow type is small. At lower shear rates (or more properly, lower Peclet number Pe based upon the rotational diffusivity), the elimination occurs through a homoclinic bifurcation, the transition being rather abrupt as the flow type is changed. At higher Pe, periodic behavior is suppressed more gradually through a Hopf bifurcation, with tumbling being replaced by wagging and negative θ flow-aligning, where θ is the angle of the director in the shear plane. The effect of these perturbations on rheological behavior is also explored. As the flow is made slightly more extensional, the zero-shear rate limiting value of the generalized viscosity η decreases dramatically, due to the slowing down of tumbling as the system approaches a homoclinic orbit; as Pe is increased, the viscosity rises again before falling, due to the induction of wagging behavior where tumbling would normally prevail in simple shear. Finally, it is found when the flow type is changed sufficiently, the interesting, non-monotonic behavior of rheological functions seen in simple shear of LCPs is replaced by monotonic behavior, even though the flow is still relatively close to simple shear.     

Numerical simulation of pulsatile non-Newtonian flow in the carotid artery bifurcation

Both clinical and post mortem studies indicate that, in humans, the carotid sinus of the carotid artery bifurcation is one of the favored sites for the genesis and development of atherosclerotic lesions. Hemodynamic factors have been suggested to be important in atherogenesis. To understand the correlation between atherogenesis and fluid dynamics in the carotid sinus, the blood flow in artery was simulated numerically. In those studies, the property of blood was treated as an incompressible, Newtonian fluid. In fact, however, the blood is a complicated non-Newtonian fluid with shear thinning and viscoelastic properties, especially when the shear rate is low. A variety of non-Newtonian models have been applied in the numerical studies. Among them, the Casson equation was widely used. However, the Casson equation agrees well only when the shear rate is less than 10 s-1. The flow field of the carotid bifurcation usually covers a wide range of shear rate. We therefore believe that it may not be sufficient to describe the property of blood only using the Casson equation in the whole flow field of the carotid bifurcation. In the present study, three different blood constitutive models, namely, the Newtonian, the Casson and the hybrid fluid constitutive models were used in the flow simulation of the human carotid bifurcation. The results were compared among the three models. The results showed that the Newtonian model and the hybrid model had verysimilar distributions of the axial velocity, secondary flow and wall shear stress, but the Casson model resulted in significant differences in these distributions from the other two models. This study suggests that it is not appropriate to only use the Casson equation to simulate the whole flow field of the carotid bifurcation, and on the other hand, Newtonian fluid is a good approximation to blood for flow simulations in the carotid artery bifurcation.     

Numerical simulation of fluid–particle hydrodynamics in a rectangular spouted vessel

A three-dimensional, Eulerian simulation was developed to describe isothermal, two-phase flow of the continuous (water) and dispersed (solid particles) phases in a rectangular spouted vessel. The mass and momentum conservation equations for each phase were solved using the finite volume technique, which treats each phase separately, while coupling them through drag, turbulence, and energy dissipation due to particle fluctuations. Particle–particle interactions via friction were also included.     

Curvilinear smoothed particle hydrodynamics

We suggest a new set of equations to employ smoothed particle hydrodynamics (SPH) in a curvilinear space, and we refer to it as curvSPH. In classical SPH, the horizontal and vertical resolution of discretization is supposed to be equal for fluid particles. However, curvSPH makes the horizontal and vertical resolutions independent from each other. This is performed by transformation of physical space into an appropriate computational space with a different scale in horizontal and vertical directions. Solving a problem using SPH in a curvilinear space also provides capability to model curved boundaries as straight lines. In classical SPH, special care is needed to reach a uniform mass distribution along curved boundaries; however, producing uniform mass distribution along a line using curvSPH is straight forward. Different simulations, including simulation of a flip bucket are performed to demonstrate the applicability of the proposed method. Good agreement of results with experimental data and classical SPH confirms the capabilities of curvSPH. Copyright © 2016 John Wiley & Sons, Ltd.     

Numerical simulation of the non-Newtonian blood flow through a mechanical aortic valve

This work focuses on the comparison between Newtonian and non-Newtonian blood flows through a bileaflet mechanical heart valve in the aortic root. The blood, in fact, is a concentrated suspension of cells, mainly red blood cells, in a Newtonian matrix, the plasma, and consequently its overall behavior is that of a non-Newtonian fluid owing to the action of the cells’ membrane on the fluid part. The common practice, however, assumes the blood in large vessels as a Newtonian fluid since the shear rate is generally high and the effective viscosity becomes independent of the former. In this paper, we show that this is not always the case even in the aorta, the largest artery of the systemic circulation, owing to the pulsatile and transitional nature of the flow. Unexpectedly, for most of the pulsating cycle and in a large part of the fluid volume, the shear rate is smaller than the threshold level for the blood to display a constant effective viscosity and its shear thinning character might affect the system dynamics. A direct inspection of the various flow features has shown that the valve dynamics, the transvalvular pressure drop and the large-scale features of the flow are very similar for the Newtonian and non-Newtonian fluid models. On the other hand, the mechanical damage of the red blood cells (hemolysis), induced by the altered stress values in the flow, is larger for the non-Newtonian fluid model than for the Newtonian one.     

钢箱内部爆炸破坏的SPH数值模拟

随着恐怖袭击的不断演化,船舶、桥梁等以钢箱为主要支撑的战略性结构逐渐成为恐怖袭击和敌方军事打击的重要目标。本文中采用光滑粒子流体动力学(smoothed particle hydrodynamics, SPH)方法对钢箱内部爆炸过程进行了数值模拟,对比实验,分析了钢箱内部爆炸时钢箱表面的变形过程,得到了钢箱表面挠度值的变化趋势、不同时刻钢箱表面压强和von Mises应力的分布情况、钢箱表面中心点处速度和压强的变化趋势,验证了SPH方法在模拟钢箱内部爆炸问题上的有效性。通过进一步数值模拟,探讨了钢箱内部不同位置爆炸时钢箱的破坏形式和损伤程度,结果表明:炸药在钢箱内部角隅处爆炸时,钢箱的损伤程度最严重;炸药在钢箱内部正中心爆炸时,钢箱的损伤程度最轻。     

Three‐dimensional remeshed smoothed particle hydrodynamics for the simulation of isotropic turbulence

We present a remeshed particle‐mesh method for the simulation of three‐dimensional compressible turbulent flow. The method is related to the meshfree smoothed particle hydrodynamics method, but the present method introduces a mesh for efficient calculation of the pressure gradient, and laminar and turbulent diffusion. In addition, the mesh is used to remesh (reorganise uniformly) the particles to ensure a regular particle distribution and convergence of the method. The accuracy of the presented methodology is tested for a number of benchmark problems involving two‐ and three‐dimensional Taylor‐Green flow, thin double shear layer, and three‐dimensional isotropic turbulence. Two models were implemented, direct numerical simulations, and Smagorinsky model. Taking advantage of the Lagrangian advection, and the finite difference efficiency, the method is capable of providing quality simulations while maintaining its robustness and versatility.     

An incompressible smoothed particle hydrodynamics scheme for Newtonian/non-Newtonian multiphase flows including semi-analytical solutions for two-phase inelastic Poiseuille flows

An incompressible smoothed particle hydrodynamics (ISPH) method is developed for the modeling of multiphase Newtonian and inelastic non-Newtonian flows at low density ratios. This new method is the multiphase extension of Xenakis et al, J. Non-Newtonian Fluid Mech., 218, 1-15, which has been shown to be stable and accurate, with a virtually noise-free pressure field for single-phase non-Newtonian flows. For the validation of the method a semi-analytical solution of a two-phase Newtonian/non-Newtonian (inelastic) Poiseuille flow is derived. The developed method is also compared with the benchmark multiphase case of the Rayleigh Taylor instability and a submarine landslide, thereby demonstrating capability in both Newtonian/Newtonian and Newtonian/non-Newtonian two-phase applications. Comparisons with analytical solutions, experimental and previously published results are conducted and show that the proposed methodology can accurately predict the free-surface and interface profiles of complex incompressible multi-phase flows at low-density ratios relevant, for example, to geophysical environmental applications.     

Turbulent flow of non-Newtonian substances

Summary A unique shear stress-shear rate relationship exists for laminar flow of any time independent substance in a tube, whereas this is not the case for turbulent flow. In order to obtain a unique relationship for turbulent flow, a new approach based on the elementary theoretical interpretation of experimental data is adopted in the present paper. In particular, wall shear stress is found to be a unique function of a new turbulent pseudo shear rate term. In this relationship therè are two parameters which characterize a given substance — the limiting viscosity at high shear rateµ _m and a factor _m which takes into account modification of turbulent structure by the non-Newtonian properties. Both of these parameters must be determined experimentally. Methods of predicting pressure gradients and of scaling up are outlined. In applying the approach to suspensions in which the solid phase has a density greater than that of the liquid medium, it may be important to determine the increment in shear stress equivalent to the energy required to maintain the solid particles in suspension.The validity of this approach is confirmed by data for the flow of a variety of substances including kaolin suspensions and Carbopol solutions in tubes ranging in diameter from 1.5 to 20 mm. Nomenclature C volume fraction solid in suspension - D tube diameter - f Darcy-Weisbach friction factor - g gravitational acceleration - K _s proportionality constant defined by eq. [10] - L length of tube - P pressure - Re Reynolds number - t exponent defined by eq. [1] - V mean velocity - V ^* volume of particles in pipe lengthL - W settling velocity of particles - _m factor defined by eq. [1] - shear rate - turbulent pseudo shear rate defined by eqs. [8] and [9] - _w wall shear stress - ( _w) _s increment in wall shear stress due to presence of settling particles - µ _m limiting viscosity at high rate of shear - ₁ density of carrier liquid - _m density of mixture - _s density of solid Professor of Chemical Engineering, University of Toronto and scientific advisor to Worthington (Canada) Ltd.With 8 figures     

Three-dimensional simulation of two viscoelastic droplets impacting onto a rigid plate using smoothed particle hydrodynamics

In this article, a smoothed particle hydrodynamics method is developed to simulate the dynamic process of the impact of two viscoelastic droplets onto a rigid plate. The Oldroyd-B fluid is considered as the rheological model to describe the viscoelastic characteristics. An artificial stress is added into the momentum equation to remove the tensile instability. The solution of the problem of two successive impacts of droplets are demonstrated to be in good agreement with the literature data. The problem of two droplets impacting simultaneously onto a rigid plate is investigated.