When ordinary Smoothed Particle Hydrodynamics （SPH） method is used to simulate wave propagation in a wave tank, it is usually observed that the wave height decays and the wave length elongates along the direction of...When ordinary Smoothed Particle Hydrodynamics （SPH） method is used to simulate wave propagation in a wave tank, it is usually observed that the wave height decays and the wave length elongates along the direction of wave propagation. Accompanied with this phenomenon, the pressure under water decays either and shows a big oscillation simultaneously. The reason is the natural potential tensile instability of modeling water motion with ordinary SPH which is caused by particle negative stress in the computation. I＇o deal with the problems, a new sextic kernel function is proposed to reduce this instability. An appropriate smooth length is given and its computation criterion is also suggested. At the same time, a new kind dynamic boundary condition is introduced. Based on these improvements, the new SPH method named stability improved SPH （SISPH） can simulate the wave propagation well. Both the water surface and pressure can be well expressed and the oscillation of pressure is nearly eliminated. Compared with other improved methods, SISPH can truly reveal the physical reality without bringing some new problems in a simple way.展开更多
Two kinds of fractures can be observed in the SPH （smoothed particle hydrodynamics） simulations, which are the physical fracture and the numerical fracture. The physical one exists in reality, while the numerical on...Two kinds of fractures can be observed in the SPH （smoothed particle hydrodynamics） simulations, which are the physical fracture and the numerical fracture. The physical one exists in reality, while the numerical one is fictitious. This paper presents the effects of both fractures and proposes a simple adding particle technique to avoid the numerical fracture. The real physical fracture is then figured out by using an applicable fracture criterion. Firstly, the effect of the numerical fracture on the computational accuracy is investigated by introducing the artificial fracture in a model of wave propagation. Secondly, a simple adding particle technique is proposed and validated by a three dimensional bending test. Finally, the experiments of penetration on the skin of aircrafts are simulated by both the initial SPH method and the improved method with the adding particle technique. The results show that the improved SPH method can describe the physical fracture very well with better accuracy.展开更多
In this paper,the smoothed particle hydrodynamics(SPH)method is employed in modeling and numerical simulation of droplet coalescence.Considering the effect of tangential force on boundary material,besides normal force...In this paper,the smoothed particle hydrodynamics(SPH)method is employed in modeling and numerical simulation of droplet coalescence.Considering the effect of tangential force on boundary material,besides normal force,tangential force is also introduced in the continuum surface force(CSF)model.The formation of droplet,the coalescence processes of two droplets and three droplets are simulated by the modified CSF model.The validity of the modified model is verified from the aspects of the morphological change of the droplet,the smoothness of free surface and the conservation of the centroid of the system.Compared with finite element method,the results of the modified CSF model show that tangential force plays a crucial role in the CSF model when dealing with model boundary with curves and sharp angles.展开更多
This paper investigates the effects of charge parameters of the underwater contact explosion based on the axisymmetric smoothed particle hydrodynamics （SPH） method. The dynamic boundary particle is proposed to impro...This paper investigates the effects of charge parameters of the underwater contact explosion based on the axisymmetric smoothed particle hydrodynamics （SPH） method. The dynamic boundary particle is proposed to improve the pressure fluctuation and numerical accuracy near the symmetric axis. An in-depth study is carried out over the influence of charge shapes and detonation modes on the near-field loads in terms of the peak pressure and impulse of shock waves. For different charge shapes, the cylindrical charge with different length-diameter ratios may cause strong directivity of peak pressure and impulse in the near field. Compared with spherical charge, the peak pressure of cylindrical charge may be either weakened or enhanced in different directions. Within a certain range, the greater the length-diameter ratio is, the more obvious the effect will be. The weakened ratio near the detonation end may reach 25% approximately, while the enhanced ratio may reach around 20% in the opposite direction. However, the impulse in different directions seems to be uniform. For different detonation modes, compared with point-source explosion, the peak pressure of plane-source explosion is enhanced by about 5%. Besides, the impulse of plane-source explosion is enhanced by around 5% near the detonation end, but close to those of the point-source explosion in other directions. Based on the material constitutive relation in the axisymmetric coordinates, a simple case of underwater contact explosion is simulated to verify the above conclusions, showing that the charge parameters of underwater contact explosion should not be ignored.展开更多
The water mitigation effect on the propagation of shock wave was investigated numerically. The traditional smoothed particle hydrodynamics (SPH) method was modified based on Riemann solution. The comparison of numeric...The water mitigation effect on the propagation of shock wave was investigated numerically. The traditional smoothed particle hydrodynamics (SPH) method was modified based on Riemann solution. The comparison of numerical results with the analytical solution indicated that the modified SPH method has more advantages than the traditional SPH method. Using the modified SPH algorithm, a series of one-dimensional planar wave propagation problems were investigated, focusing on the influence of the air-gap between the high-pressure air and water and the thickness of water. The numerical results showed that water mitigation effect is significant. Up to 60% shock wave pressure reduction could be achieved with the existence of water, and the shape of shock wave was also changed greatly. It is seemly that the small air-gap between the high-pressure air and water has more influence on water mitigation effect.展开更多
In this paper,we study a class of contact smoothed particle hydrodynamics(SPH)by introducing Riemann solvers and using high-order limiters.In particular,a promising concept ofWENO interpolation as limiter is presented...In this paper,we study a class of contact smoothed particle hydrodynamics(SPH)by introducing Riemann solvers and using high-order limiters.In particular,a promising concept ofWENO interpolation as limiter is presented in the reconstruction process.The physical values relating interactional particles used as the initial values of the Riemann problemcan be reconstructed by the Taylor series expansion.The contact solvers of the Riemann problem at contact points are incorporated in SPH approximations.In order to keep the fluid density at the wall rows to be consistent with that of the inner fluid wall boundaries,several lines of dummy particles are placed outside of the solid walls,which are assigned according to the initial configuration.At last,the method is applied to compressible flows with sharp discontinuities such as the collision of two strong shocks and the interaction of two blast waves and so on.The numerical results indicate that the method is capable of handling sharp discontinuity and efficiently reducing unphysical oscillations.展开更多
Based on the SPH method,this paper investigates a number of key issues in shock dynamics,including the nonlinear phenomenon and its influencing factors in the propagation process of shock waves,and the process and mec...Based on the SPH method,this paper investigates a number of key issues in shock dynamics,including the nonlinear phenomenon and its influencing factors in the propagation process of shock waves,and the process and mechanism of the interaction between shock waves and the surface of an object.Through the design of the experiments and the analysis of the results,the main conclusions and contributions of this study are drawn.At the same time,deficiencies and directions for further improvement are also pointed out,and prospects and suggestions for future research on impact dynamics are presented.This study is of great significance in promoting the research of impact dynamics.展开更多
This paper applies the meshfree Smoothed Particle Hydrodynamics （SPH） method with Graphical Processing Unit （GPU） parallel computing technique to investigate the highly complex 3-D dam-break flow in urban areas in...This paper applies the meshfree Smoothed Particle Hydrodynamics （SPH） method with Graphical Processing Unit （GPU） parallel computing technique to investigate the highly complex 3-D dam-break flow in urban areas including underground spaces. Taking the advantage of GPUs parallel computing techniques, simulations involving more than 107 particles can be achieved. We use a virtual geometric plane boundary to handle the outermost solid wall in order to save considerable video card memory for the GPU computing. To evaluate the accuracy of the new GPU-based SPH model, qualitative and quantitative comparison to a real flooding experiment is performed and the results of a numerical model based on Shallow Water Equations （SWEs） is given with good accuracy. With the new GPU-based SPH model, the effects of the building layouts and underground spaces on the propagation of dambreak flood through an intricate city layout are examined.展开更多
In this paper,a novel unconditionally energy stable Smoothed Particle Hydrodynamics(SPH)method is proposed and implemented for incompressible fluid flows.In this method,we apply operator splitting to break the momentu...In this paper,a novel unconditionally energy stable Smoothed Particle Hydrodynamics(SPH)method is proposed and implemented for incompressible fluid flows.In this method,we apply operator splitting to break the momentum equation into equations involving the non-pressure term and pressure term separately.The idea behind the splitting is to simplify the calculation while still maintaining energy stability,and the resulted algorithm,a type of improved pressure correction scheme,is both efficient and energy stable.We show in detail that energy stability is preserved at each full-time step,ensuring unconditionally numerical stability.Numerical examples are presented and compared to the analytical solutions,suggesting that the proposed method has better accuracy and stability.Moreover,we observe that if we are interested in steady-state solutions only,our method has good performance as it can achieve the steady-state solutions rapidly and accurately.展开更多
基金financially supported by the National Natural Science Foundation of China(Grant Nos.51579038 and 51490672)the National Basic Research Program of China(Grant No.2013CB036101)
文摘When ordinary Smoothed Particle Hydrodynamics （SPH） method is used to simulate wave propagation in a wave tank, it is usually observed that the wave height decays and the wave length elongates along the direction of wave propagation. Accompanied with this phenomenon, the pressure under water decays either and shows a big oscillation simultaneously. The reason is the natural potential tensile instability of modeling water motion with ordinary SPH which is caused by particle negative stress in the computation. I＇o deal with the problems, a new sextic kernel function is proposed to reduce this instability. An appropriate smooth length is given and its computation criterion is also suggested. At the same time, a new kind dynamic boundary condition is introduced. Based on these improvements, the new SPH method named stability improved SPH （SISPH） can simulate the wave propagation well. Both the water surface and pressure can be well expressed and the oscillation of pressure is nearly eliminated. Compared with other improved methods, SISPH can truly reveal the physical reality without bringing some new problems in a simple way.
基金supported by the National Natural Science Foundation of China (No.10577016)the 111 Project (No. B07050)the program for 2008 New Century Excellent Talents in University (No.NCET080454)
文摘Two kinds of fractures can be observed in the SPH （smoothed particle hydrodynamics） simulations, which are the physical fracture and the numerical fracture. The physical one exists in reality, while the numerical one is fictitious. This paper presents the effects of both fractures and proposes a simple adding particle technique to avoid the numerical fracture. The real physical fracture is then figured out by using an applicable fracture criterion. Firstly, the effect of the numerical fracture on the computational accuracy is investigated by introducing the artificial fracture in a model of wave propagation. Secondly, a simple adding particle technique is proposed and validated by a three dimensional bending test. Finally, the experiments of penetration on the skin of aircrafts are simulated by both the initial SPH method and the improved method with the adding particle technique. The results show that the improved SPH method can describe the physical fracture very well with better accuracy.
文摘In this paper,the smoothed particle hydrodynamics(SPH)method is employed in modeling and numerical simulation of droplet coalescence.Considering the effect of tangential force on boundary material,besides normal force,tangential force is also introduced in the continuum surface force(CSF)model.The formation of droplet,the coalescence processes of two droplets and three droplets are simulated by the modified CSF model.The validity of the modified model is verified from the aspects of the morphological change of the droplet,the smoothness of free surface and the conservation of the centroid of the system.Compared with finite element method,the results of the modified CSF model show that tangential force plays a crucial role in the CSF model when dealing with model boundary with curves and sharp angles.
基金supported by the National Natural Science Foundation of China(No.51379039)the Excellent Young Scientists Fund(No.51222904)
文摘This paper investigates the effects of charge parameters of the underwater contact explosion based on the axisymmetric smoothed particle hydrodynamics （SPH） method. The dynamic boundary particle is proposed to improve the pressure fluctuation and numerical accuracy near the symmetric axis. An in-depth study is carried out over the influence of charge shapes and detonation modes on the near-field loads in terms of the peak pressure and impulse of shock waves. For different charge shapes, the cylindrical charge with different length-diameter ratios may cause strong directivity of peak pressure and impulse in the near field. Compared with spherical charge, the peak pressure of cylindrical charge may be either weakened or enhanced in different directions. Within a certain range, the greater the length-diameter ratio is, the more obvious the effect will be. The weakened ratio near the detonation end may reach 25% approximately, while the enhanced ratio may reach around 20% in the opposite direction. However, the impulse in different directions seems to be uniform. For different detonation modes, compared with point-source explosion, the peak pressure of plane-source explosion is enhanced by about 5%. Besides, the impulse of plane-source explosion is enhanced by around 5% near the detonation end, but close to those of the point-source explosion in other directions. Based on the material constitutive relation in the axisymmetric coordinates, a simple case of underwater contact explosion is simulated to verify the above conclusions, showing that the charge parameters of underwater contact explosion should not be ignored.
基金Supported by National Natural Science Foundation of China(No.50638030 and 50525825)National Science and Technology Support Program(No.2006BAJ13B02)
文摘The water mitigation effect on the propagation of shock wave was investigated numerically. The traditional smoothed particle hydrodynamics (SPH) method was modified based on Riemann solution. The comparison of numerical results with the analytical solution indicated that the modified SPH method has more advantages than the traditional SPH method. Using the modified SPH algorithm, a series of one-dimensional planar wave propagation problems were investigated, focusing on the influence of the air-gap between the high-pressure air and water and the thickness of water. The numerical results showed that water mitigation effect is significant. Up to 60% shock wave pressure reduction could be achieved with the existence of water, and the shape of shock wave was also changed greatly. It is seemly that the small air-gap between the high-pressure air and water has more influence on water mitigation effect.
基金supported by the Fundamental Research Funds of Central University(B103093)Key Laboratory Foundation of Coastal Disasters.
文摘In this paper,we study a class of contact smoothed particle hydrodynamics(SPH)by introducing Riemann solvers and using high-order limiters.In particular,a promising concept ofWENO interpolation as limiter is presented in the reconstruction process.The physical values relating interactional particles used as the initial values of the Riemann problemcan be reconstructed by the Taylor series expansion.The contact solvers of the Riemann problem at contact points are incorporated in SPH approximations.In order to keep the fluid density at the wall rows to be consistent with that of the inner fluid wall boundaries,several lines of dummy particles are placed outside of the solid walls,which are assigned according to the initial configuration.At last,the method is applied to compressible flows with sharp discontinuities such as the collision of two strong shocks and the interaction of two blast waves and so on.The numerical results indicate that the method is capable of handling sharp discontinuity and efficiently reducing unphysical oscillations.
文摘Based on the SPH method,this paper investigates a number of key issues in shock dynamics,including the nonlinear phenomenon and its influencing factors in the propagation process of shock waves,and the process and mechanism of the interaction between shock waves and the surface of an object.Through the design of the experiments and the analysis of the results,the main conclusions and contributions of this study are drawn.At the same time,deficiencies and directions for further improvement are also pointed out,and prospects and suggestions for future research on impact dynamics are presented.This study is of great significance in promoting the research of impact dynamics.
基金Project supported by the National Basic Research Development Program of China(973 Program,No.2012CB719705)the National Natural Science Foundation of China(Grant Nos.91024032,70833003)
文摘This paper applies the meshfree Smoothed Particle Hydrodynamics （SPH） method with Graphical Processing Unit （GPU） parallel computing technique to investigate the highly complex 3-D dam-break flow in urban areas including underground spaces. Taking the advantage of GPUs parallel computing techniques, simulations involving more than 107 particles can be achieved. We use a virtual geometric plane boundary to handle the outermost solid wall in order to save considerable video card memory for the GPU computing. To evaluate the accuracy of the new GPU-based SPH model, qualitative and quantitative comparison to a real flooding experiment is performed and the results of a numerical model based on Shallow Water Equations （SWEs） is given with good accuracy. With the new GPU-based SPH model, the effects of the building layouts and underground spaces on the propagation of dambreak flood through an intricate city layout are examined.
基金This work is partially supported by King Abdullah University of Science and Technology(KAUST)through the grants BAS/1/1351-01,URF/1/4074-01,and URF/1/3769-01.
文摘In this paper,a novel unconditionally energy stable Smoothed Particle Hydrodynamics(SPH)method is proposed and implemented for incompressible fluid flows.In this method,we apply operator splitting to break the momentum equation into equations involving the non-pressure term and pressure term separately.The idea behind the splitting is to simplify the calculation while still maintaining energy stability,and the resulted algorithm,a type of improved pressure correction scheme,is both efficient and energy stable.We show in detail that energy stability is preserved at each full-time step,ensuring unconditionally numerical stability.Numerical examples are presented and compared to the analytical solutions,suggesting that the proposed method has better accuracy and stability.Moreover,we observe that if we are interested in steady-state solutions only,our method has good performance as it can achieve the steady-state solutions rapidly and accurately.
