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.展开更多
Based on the improved interpolating moving least-squares (ⅡMLS) method and the Galerkin weak form, an improved interpolating element-free Galerkin (ⅡEFG) method is presented for two-dimensional elasticity proble...Based on the improved interpolating moving least-squares (ⅡMLS) method and the Galerkin weak form, an improved interpolating element-free Galerkin (ⅡEFG) method is presented for two-dimensional elasticity problems in this paper. Compared with the interpolating moving least-squares (IMLS) method presented by Lancaster, the ⅡMLS method uses the nonsingular weight function. The number of unknown coefficients in the trial function of the ⅡMLS method is less than that of the MLS approximation and the shape function of the ⅡMLS method satisfies the property of Kronecker δ function. Thus in the ⅡEFG method, the essential boundary conditions can be applied directly and easily, then the numerical solutions can be obtained with higher precision than those obtained by the interpolating element-free Galerkin (IEFG) method. For the purposes of demonstration, four numerical examples are solved using the ⅡEFG method.展开更多
基金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.
基金Project supported by the National Natural Science Foundation of China(Grant No.11171208)the Shanghai Leading Academic Discipline Project,China(Grant No.S30106)
文摘Based on the improved interpolating moving least-squares (ⅡMLS) method and the Galerkin weak form, an improved interpolating element-free Galerkin (ⅡEFG) method is presented for two-dimensional elasticity problems in this paper. Compared with the interpolating moving least-squares (IMLS) method presented by Lancaster, the ⅡMLS method uses the nonsingular weight function. The number of unknown coefficients in the trial function of the ⅡMLS method is less than that of the MLS approximation and the shape function of the ⅡMLS method satisfies the property of Kronecker δ function. Thus in the ⅡEFG method, the essential boundary conditions can be applied directly and easily, then the numerical solutions can be obtained with higher precision than those obtained by the interpolating element-free Galerkin (IEFG) method. For the purposes of demonstration, four numerical examples are solved using the ⅡEFG method.