Splitting for Highly Dissipative Smoothed Particle Dynamics

The method of smoothed dissipative particle dynamics (SDPD) is a novel coarse grained method for simulation of complex fluids. It has some advantages over more traditional particles based methods (Espanol and Warren, Europhys. Lett. 30(4):191-196, 1995). But one of the problems common for particle based simulations of microfluid system takes place also for SDPD: it fails to realize Schmidt number of O(10(3)) typical of liquids. In present paper we apply the implicit numerical scheme that allows significantly increase time step in SDPD and perform simulation for larger Schmidt number. Simulations using this methods show close agreement with serial solutions for Couette and Poiseuille flows. The results of benchmarks based on temperature control are presented. The dependence of self-diffusion coefficient D on kinematic viscosity is examined and found to be in agreement with empirical observations (Li and Chang, J. Chem. Phys. 23(3):518-520, 1955).