1i] acting on particles at time t, the fundamental QDPD cycle, repeated for as many timesteps as are in a simulation, is
3] coordinates make sure that particles are assigned to one of the cells in the QDPD simulation box (i.
In treating edge cells, a "ghost" layer of cells is added to the QDPD simulation box.
Because QDPD in its present form is being used to study the steady-shear viscosity of a suspension of solid inclusions in a Newtonian fluid, there is a shear boundary condition at the topmost layer of the QDPD simulation box, implemented with the Lees-Edwards boundary conditions  (pp.
QDPD was originally written in Fortran 77 as a serial program.
This extra communication step is necessary in the QDPD method since the interparticle force calculation involves the use of a random number for thermal effects and momentum conservation requires that the same random number be used in the equal and opposite force calculation.
While the replicated data approach of the previous section has been the workhorse of QDPD work for some time now, it has had its disadvantages.
Second, the QDPD technique is being applied to suspensions, so there are two types of particles: "free" particles and particles belonging to solid inclusions such as ellipsoids.
6, Kirk Kern of SGI for the 64 processor SGI runs cited in the QDPD section (Sec.