We apply smoothed dissipative particle dynamics (SDPD) [Español and Revenga, Phys. Rev. E 67, 026705 (2003)] to model solid particles in suspension. SDPD is a thermodynamically consistent version of smoothed particle hydrodynamics (SPH) and can be interpreted as a multiscale particle framework linking the macroscopic SPH to the mesoscopic dissipative particle dynamics (DPD) method. Rigid structures of arbitrary shape embedded in the fluid are modeled by frozen particles on which artificial velocities are assigned in order to satisfy exactly the no-slip boundary condition on the solid-liquid interface. The dynamics of the rigid structures is decoupled from the solvent by solving extra equations for the rigid body translational/angular velocities derived from the total drag/torque exerted by the surrounding liquid. The correct scaling of the SDPD thermal fluctuations with the fluid-particle size allows us to describe the behavior of the particle suspension on spatial scales ranging continuously from the diffusion-dominated regime typical of sub-micron-sized objects towards the non-Brownian regime characterizing macro-continuum flow conditions. Extensive tests of the method are performed for the case of two/three dimensional bulk particle-system both in Brownian/non-Brownian environment showing numerical convergence and excellent agreement with analytical theories. Finally, to illustrate the ability of the model to couple with external boundary geometries, the effect of confinement on the diffusional properties of a single sphere within a micro-channel is considered, and the dependence of the diffusion coefficient on the wall-separation distance is evaluated and compared with available analytical results.

Separation of electrolysis products using fluidic inertial forces in a 3D printed flow cell.

We preset a computational study of bending models for the curvature elasticity of lipid bilayer membranes that are relevant for simulations of vesicles and red blood cells. We compute bending energy and forces on triangulated meshes and evaluate and extend four well established schemes for their approximation: Kantor and Nelson [1], Jülicher [2], Gompper and Kroll [3] and Meyer et. al. [4], termed A, B, C, D. We present a comparative study of these four schemes on the minimal bending model and propose extensions for schemes B, C and D. These extensions incorporate the reference state and non-local energy to account for the spontaneous curvature, bilayer coupling, and area-difference elasticity models. Our results indicate that the proposed extensions enhance the models to account for shape transformation including budding/vesiculation as well as for non-axisymmetric shapes. We find that the extended scheme B is superior to the rest in terms of accuracy, and robustness as well as simplicity of implementation. We demonstrate the capabilities of this scheme on several benchmark problems including the budding-vesiculating process and the reproduction of the phase diagram of vesicles.

We present a model for a polymer molecule in solution based on smoothed dissipative particle dynamics (SDPD) [Español and Revenga, Phys. Rev. E 67, 026705 (2003)]. This method is a thermodynamically consistent version of smoothed particle hydrodynamics able to discretize the Navier-Stokes equations and, at the same time, to incorporate thermal fluctuations according to the fluctuation-dissipation theorem. Within the framework of the method developed for mesoscopic multiphase flows by Hu and Adams [J. Comput. Phys. 213, 844 (2006)], we introduce additional finitely extendable nonlinear elastic interactions between particles that represent the beads of a polymer chain. In order to assess the accuracy of the technique, we analyze the static and dynamic conformational properties of the modeled polymer molecule in solution. Extensive tests of the method for the two-dimensional (2D) case are performed, showing good agreement with the analytical theory. Finally, the effect of confinement on the conformational properties of the polymer molecule is investigated by considering a 2D microchannel with gap H varying between 1 and 10 microm , of the same order as the polymer gyration radius. Several SDPD simulations are performed for different chain lengths corresponding to N=20-100 beads, giving a universal behavior of the gyration radius R_{G} and polymer stretch X as functions of the channel gap when normalized properly.

Adsorption on heterogeneous surfaces with three basic energy distribution models (uniform model, exponential model, and normal-like model) is studied. Exact analytical solutions of the adsorption isotherms and the heats of adsorption are derived for the uniform and exponential models, and, with these solutions including a numerical solution for the normal-like model, the behavior of the differential heat of adsorption and the "apparent" standard adsorption entropy concerning the overall surface is described as a function of coverage and temperature. The approximations underlying the isotherms and heats of adsorption in the Temkin, Freundlich, and Langmuir-Freundlich types of adsorption are rationalized. By comparing these empirical formulas to the exact solutions, the level of these approximations is found to be identical, which is similar to the "condensation approximation". Their preconditions are that either the temperature is low enough, or the surface is strongly heterogeneous. Generally, they are suitable for the middle coverage range. The exact solutions provide a method to obtain more information on the heats, entropy, and heterogeneity of the catalyst surface from the calorimetric measurement of the heat of adsorption.