search

erizmr / SPH_Taichi

A high-performance implementation of SPH in Taichi.
195 stars 39 forks source link

WCPSH dynamic rigid body may be error #29

Open jiajun-c opened 6 months ago

jiajun-c commented 6 months ago 
        elif self.ps.material[p_j] == self.ps.material_solid:
            boundary_viscosity = 0.0
            # Boundary neighbors
            ## Akinci2012
            f_v = d * boundary_viscosity * (self.density_0 * self.ps.m_V[p_j] / (self.ps.density[p_i])) * v_xy / (
                r.norm()**2 + 0.01 * self.ps.support_radius**2) * self.cubic_kernel_derivative(r)
            ret += f_v

If boundary_viscosity is zero, the f_v will be zero too.....