In this work, an enhanced treatment of the solid boundaries is proposed for smoothed particle hydrodynamics with implicit time integration scheme （Implicit SPH）. Three types of virtual particles, i.e., boundary part...In this work, an enhanced treatment of the solid boundaries is proposed for smoothed particle hydrodynamics with implicit time integration scheme （Implicit SPH）. Three types of virtual particles, i.e., boundary particles, image particles and mirror particles, are used to impose boundary conditions. Boundary particles are fixed on the solid boundary, and each boundary particle is associated with two fixed image particles inside the fluid domain and two fixed mirror particles outside the fluid domain. The image particles take the flow properties through fluid particles with moving least squares （MLS） interpolation and the properties of mirror particles can be obtained by the corresponding image particles. A repulsive force is also applied for boundary particles to prevent fluid particles from unphysical penetra- tion through solid boundaries. The new boundary treatment method has been validated with five numerical examples. All the numerical results show that Implicit SPH with this new boundary-treatment method can obtain accurate results for non-Newtonian fluids as well as Newtonian fluids, and this method is suitable for complex solid boundaries and can be easily extended to 3D problems.展开更多
In this paper, first we calculate finite-difference coefficients of implicit finite- difference methods （IFDM） for the first and second-order derivatives on normal grids and first- order derivatives on staggered gri...In this paper, first we calculate finite-difference coefficients of implicit finite- difference methods （IFDM） for the first and second-order derivatives on normal grids and first- order derivatives on staggered grids and find that small coefficients of high-order IFDMs exist. Dispersion analysis demonstrates that omitting these small coefficients can retain approximately the same order accuracy but greatly reduce computational costs. Then, we introduce a mirrorimage symmetric boundary condition to improve IFDMs accuracy and stability and adopt the hybrid absorbing boundary condition （ABC） to reduce unwanted reflections from the model boundary. Last, we give elastic wave modeling examples for homogeneous and heterogeneous models to demonstrate the advantages of the proposed scheme.展开更多
Implicit smoothed particle hydrodynamics method has been proposed to overcome the problem that only very small time steps can be used for high viscosity fluids(such as power law fluids)in order to obtain a stable simu...Implicit smoothed particle hydrodynamics method has been proposed to overcome the problem that only very small time steps can be used for high viscosity fluids(such as power law fluids)in order to obtain a stable simulation.However,the pressure field is difficult to simulate correctly with this method because the numerical high-frequency noise on the pressure field cannot be removed.In this study,several improvements,which are the diffusive term in the continuity equation,the artificial viscosity and a simplified physical viscosity term in the momentum equation,are introduced,and a new boundary treatment is also proposed.The linear equations derived from the momentum equation are large-scale,sparse and positive definite but unsymmetrical,therefore,Conjugate Gradient Squared(CGS)is used to solve them.For the purpose of verifying the validity of the proposed method,Poiseuille flows with Newtonian and power law fluids are solved and compared with exact solution and traditional SPH.Drops of different fluid properties impacting a rigid wall are also simulated and compared with VOF solution.All the numerical results obtained by the proposed method agree well with available data.The proposed method shows the higher efficiency than traditional SPH and the less numerical noise on the pressure field and better stability than implicit SPH.展开更多
基金supported by the National Natural Science Foundation of China(51276192)
文摘In this work, an enhanced treatment of the solid boundaries is proposed for smoothed particle hydrodynamics with implicit time integration scheme （Implicit SPH）. Three types of virtual particles, i.e., boundary particles, image particles and mirror particles, are used to impose boundary conditions. Boundary particles are fixed on the solid boundary, and each boundary particle is associated with two fixed image particles inside the fluid domain and two fixed mirror particles outside the fluid domain. The image particles take the flow properties through fluid particles with moving least squares （MLS） interpolation and the properties of mirror particles can be obtained by the corresponding image particles. A repulsive force is also applied for boundary particles to prevent fluid particles from unphysical penetra- tion through solid boundaries. The new boundary treatment method has been validated with five numerical examples. All the numerical results show that Implicit SPH with this new boundary-treatment method can obtain accurate results for non-Newtonian fluids as well as Newtonian fluids, and this method is suitable for complex solid boundaries and can be easily extended to 3D problems.
基金supported by the National Natural Science Foundation of China(NSFC)(Grant No. 41074100)the Program for New Century Excellent Talents in University of Ministry of Education of China(Grant No. NCET-10-0812)
文摘In this paper, first we calculate finite-difference coefficients of implicit finite- difference methods （IFDM） for the first and second-order derivatives on normal grids and first- order derivatives on staggered grids and find that small coefficients of high-order IFDMs exist. Dispersion analysis demonstrates that omitting these small coefficients can retain approximately the same order accuracy but greatly reduce computational costs. Then, we introduce a mirrorimage symmetric boundary condition to improve IFDMs accuracy and stability and adopt the hybrid absorbing boundary condition （ABC） to reduce unwanted reflections from the model boundary. Last, we give elastic wave modeling examples for homogeneous and heterogeneous models to demonstrate the advantages of the proposed scheme.
基金supported by the National Natural Science Foundation of China(Grant No.51276192)
文摘Implicit smoothed particle hydrodynamics method has been proposed to overcome the problem that only very small time steps can be used for high viscosity fluids(such as power law fluids)in order to obtain a stable simulation.However,the pressure field is difficult to simulate correctly with this method because the numerical high-frequency noise on the pressure field cannot be removed.In this study,several improvements,which are the diffusive term in the continuity equation,the artificial viscosity and a simplified physical viscosity term in the momentum equation,are introduced,and a new boundary treatment is also proposed.The linear equations derived from the momentum equation are large-scale,sparse and positive definite but unsymmetrical,therefore,Conjugate Gradient Squared(CGS)is used to solve them.For the purpose of verifying the validity of the proposed method,Poiseuille flows with Newtonian and power law fluids are solved and compared with exact solution and traditional SPH.Drops of different fluid properties impacting a rigid wall are also simulated and compared with VOF solution.All the numerical results obtained by the proposed method agree well with available data.The proposed method shows the higher efficiency than traditional SPH and the less numerical noise on the pressure field and better stability than implicit SPH.
