A microscale thermophoretic turbine driven by external diffusive heat flux

Thermal conductivity effect on thermophoresis of charged spheroidal colloids in aqueous media

Thermophoresis of spheroidal colloids in aqueous media under the thermal conductivity effect is analyzed. The thermophoretic velocity and the thermodiffusion coefficient of spheroidal colloids have been formulated for extremely thin electric double layer (EDL) cases. Furthermore, a numerical thermophoretic model is built for arbitrary EDL thickness cases. The parametric studies show that the thermal conductivity mismatch of particle and liquid gives rise to a nonlinear temperature region around the spheroid, with the thickness close to the minor semiaxis. When the EDL region is thin relative to such nonlinear temperature region, the thermal conductivity effect on the thermophoresis of spheroidal colloids is significant, which strongly depends on the ratio of the minor semiaxis to the EDL thickness, the thermal conductivity ratio of particle to liquid, and the particle aspect ratio. Finally, to estimate the thermodiffusion coefficient of spheroidal colloids with arbitrary thermal conductivity, electrolyte concentration, and particle shape, the average dimensionless axial temperature gradient on the spheroidal equator plane in the EDL region is proposed.

Analytical analysis of anisotropic thermophoresis of a charged spheroidal colloid in aqueous media for extremely thin EDL cases

Thermophoresis of charged spheroids has been widely applied in biology and medical science. In this work, we report an analysis of the anisotropic thermophoresis of diluted spheroidal colloids in aqueous media for extremely thin EDL cases. Under the boundary layer approximation, we formulate the thermophoretic velocity, the thermophoretic force, and the thermodiffusion coefficient of a randomly dispersed spheroid. The parametric studies show that under the aforementioned conditions, the thermophoresis is anisotropic and its thermodiffusion coefficient should be considered as a vector, D_T . The thermodiffusion coefficient values and directions of D_T are strongly related to the aspect ratio and the angle θ between the externally applied temperature gradient and the particle's axis of revolution: The increasing aspect ratio enlarges the thermodiffusion coefficient value D_T of prolate (oblate) spheroids to a constant value when θ < 60° (θ > 45°), and it reduces D_T of prolate (oblate) spheroids to a constant value when θ > 60° (θ < 45°). The thermodiffusion coefficient direction of both prolate and oblate spheroids deviates slightly from -∇T_∞ for a small aspect ratio, and such deviation becomes serious for a large aspect ratio.

Dissipative particle dynamics with energy conservation: Isoenergetic integration and transport properties

Simulations of nano- to micro-meter scale fluidic systems under thermal gradients require consistent mesoscopic methods accounting for both hydrodynamic interactions and proper transport of energy. One such method is dissipative particle dynamics with energy conservation (DPDE), which has been used for various fluid systems with non-uniform temperature distributions. We propose an easily parallelizable modification of the velocity-Verlet algorithm based on local energy redistribution for each DPDE particle such that the total energy in a simulated system is conserved up to machine precision. Furthermore, transport properties of a DPDE fluid are analyzed in detail. In particular, an analytical approximation for the thermal conductivity coefficient is derived, which allows its a priori estimation for a given parameter set. Finally, we provide approximate expressions for the dimensionless Prandtl and Schmidt numbers, which characterize fluid transport properties and can be adjusted independently by a proper selection of model parameters. In conclusion, our results strengthen the DPDE method as a very robust approach for the investigation of mesoscopic systems with temperature inhomogeneities.

Thermophoresis of biological and biocompatible compounds in aqueous solution

With rising popularity of microscale thermophoresis for the characterisation of protein-ligand binding reactions and possible applications in microfluidic devices, there is a growing interest in considering thermodiffusion in the context of life sciences. But although the understanding of thermodiffusion in non-polar mixtures has grown rapidly in recent years, predictions for associated mixtures like aqueous solutions remain challenging. This review aims to give an overview of the literature on thermodiffusion in aqueous systems, show the difficulties in theoretical description that arise from the non-ideal behaviour of water-mixtures, and highlight the relevance of thermodiffusion in a biological context. We find that the thermodiffusion in aqueous systems is dominated by contributions from heat of transfer, hydrogen bond interactions and charge effects. However, the separation of these effects is often difficult, especially in case of biological systems where a systematic exclusion of contributions may not be feasible.

Thermoosomosis in microfluidic devices containing a temperature gradient normal to the channel walls

We analyze a microfluidic pump from the literature that utilizes a flat channel with boundary walls at different temperatures and tilted elongated pillars within in order to construct an adequate theory for designing devices in which the temperature gradient between channel walls is transformed into a longitudinal temperature gradient along the channel length. The action of the device is based on thermoosmosis in the secondary longitudinal temperature gradient associated with the specific geometry of the device, which can be described using physicochemical hydrodynamics without invoking the concept of thermophoretic force. We also describe a rotating drive device based on the same principle and design.

Thermo-orientation in fluids of arbitrarily shaped particles

Recent nonequilibrium Molecular Dynamics (NEMD) simulations revealed preferential orientation, induced by a temperature gradient, in fluids of uncharged dumbbell-like particles. The magnitude of this phenomenon, called thermo-orientation, was found to be linear in the applied temperature gradient and to increase with the difference in shape or mass between the two beads of the particles. The underlying mechanism and the microscopic determinants of the phenomenon are not obvious. Here, after examination of the general symmetry requirements for thermo-orientation, we have extended the NEMD simulations to uncharged particles of various shapes and mass distribution, including chiral cases. The numerical results are rationalized by a microscopic model, based on the assumption of local equilibrium. This allows us to correlate the thermo-orientation response of arbitrarily shaped particles to quantities that characterize their shape and mass distribution.

Dynamics of a colloidal particle near a thermoosmotic wall under illumination

The effects of light on colloidal particles in solution are multiple, including transfer of photonic linear/angular momentum and heating of the particles or their surroundings. The temperature increase around colloidal particles due to light heating can drive a thermoosmotic flow along a nearby boundary wall, which significantly influences the motion of the particles. Here we perform mesoscopic dynamics simulations to study two typical scenarios, where thermoosmosis is critical. In the first scenario, we consider a light-absorbing colloidal particle heated by uniform illumination. Depending on the fluid-wall interactions, the thermoosmotic flow generated by the wall can exert a long-ranged hydrodynamic attraction or repulsion on the "hot" Brownian particle, which leads to a strong accumulation/depletion of the particle to/from the boundary. In the second scenario, we investigate the motion of a colloidal particle confined by an optical tweezer in a light-absorbing solution. In this case, thermoosmosis can induce a unidirectional rotation of the trapped particle, whose direction is determined by thermoosmotic properties of the wall. We show that colloidal particles near a thermoosmotic wall exhibit rich dynamics. Our results can be applied for the manipulation of colloidal particles in microfluidics.

Mesoscopic modelling and simulation of soft matter

The deformability of soft condensed matter often requires modelling of hydrodynamical aspects to gain quantitative understanding. This, however, requires specialised methods that can resolve the multiscale nature of soft matter systems. We review a number of the most popular simulation methods that have emerged, such as Langevin dynamics, dissipative particle dynamics, multi-particle collision dynamics, sometimes also referred to as stochastic rotation dynamics, and the lattice-Boltzmann method. We conclude this review with a short glance at current compute architectures for high-performance computing and community codes for soft matter simulation.

Anisotropic thermophoresis

Colloidal migration in a temperature gradient is referred to as thermophoresis. In contrast to particles with a spherical shape, we show that elongated colloids may have a thermophoretic response that varies with the colloid orientation. Remarkably, this can translate into a non-vanishing thermophoretic force in the direction perpendicular to the temperature gradient. Opposite to the friction force, the thermophoretic force of a rod oriented with the temperature gradient can be larger or smaller than when oriented perpendicular to it. The precise anisotropic thermophoretic behavior clearly depends on the colloidal rod aspect ratio, and also on its surface details, which provides an interesting tunability to the devices constructed based on this principle. By means of mesoscale hydrodynamic simulations, we characterize this effect for different types of rod-like colloids.

An Optically Driven Bistable Janus Rotor with Patterned Metal Coatings

Bistable rotation is realized for a gold-coated Janus colloidal particle in an infrared optical trap. The metal coating on the Janus particles are patterned by sputtering gold on a monolayer of closely packed polystyrene particles. The Janus particle is observed to stably rotate in an optical trap. Both the direction and the rate of rotation can be experimentally controlled. Numerical calculations reveal that the bistable rotation is the result of spontaneous symmetry breaking induced by the uneven curvature of the coating patterns on the Janus sphere. Our results thus provide a simple method to construct large quantities of fully functional rotary motors for nano- or microdevices.

A microscale turbine driven by diffusive mass flux

An external diffusive mass flux is shown to be able to generate a mechanical torque on a microscale object based on anisotropic diffusiophoresis. In light of this finding, we propose a theoretical prototype micro-turbine driven purely by diffusive mass flux, which is in strong contrast to conventional turbines driven by convective mass flows. The rotational velocity of the proposed turbine is determined by the external concentration gradient, the geometry and the diffusiophoretic properties of the turbine. This scenario is validated by performing computer simulations. Our finding thus provides a new type of chemo-mechanical response which could be used to exploit existing chemical energies at small scales.

Effect of angular momentum conservation on hydrodynamic simulations of colloids

In contrast to most real fluids, angular momentum is not a locally conserved quantity in some mesoscopic simulation methods. Here we quantify the importance of this conservation in the flow fields associated with different colloidal systems. The flow field is analytically calculated with and without angular momentum conservation for the multiparticle collision dynamics (MPC) method, and simulations are performed to verify the predictions. The flow field generated around a colloidal particle moving under an external force with slip boundary conditions depends on the conservation of angular momentum, and the amplitude of the friction force is substantially affected. Interestingly, no dependence on the angular momentum conservation is found for the flow fields generated around colloids under the influence of phoretic forces. Moreover, circular Couette flow between a no-slip and a slip cylinder is investigated, which allows us to validate one of the two existing expressions for the MPC stress tensor.