We present a new mesoscale model for ionic liquids based on a low Mach number fluctuating hydrodynamics formulation for multicomponent charged species. The low Mach number approach eliminates sound waves from the fully compressible equations leading to a computationally efficient incompressible formulation. The model uses a Gibbs free-energy functional that includes enthalpy of mixing, interfacial energy, and electrostatic contributions. These lead to a new fourth-order term in the mass equations and a reversible stress in the momentum equations. We calibrate our model using parameters for [DMPI+][F6P-], an extensively studied room temperature ionic liquid (RTIL), and numerically demonstrate the formation of mesoscopic structuring at equilibrium in two and three dimensions. In simulations with electrode boundaries the measured double-layer capacitance decreases with voltage, in agreement with theoretical predictions and experimental measurements for RTILs. Finally, we present a shear electroosmosis example to demonstrate that the methodology can be used to model electrokinetic flows.

Liquid-liquid phase separation is key to understanding aqueous two-phase systems (ATPS) arising throughout cell biology, medical science, and the pharmaceutical industry. Controlling the detailed morphology of phase-separating compound droplets leads to new technologies for efficient single-cell analysis, targeted drug delivery, and effective cell scaffolds for wound healing. We present a computational model of liquid-liquid phase separation relevant to recent laboratory experiments with gelatin-polyethylene glycol mixtures. We include buoyancy and surface-tension-driven finite viscosity fluid dynamics with thermally induced phase separation. We show that the fluid dynamics greatly alters the evolution and equilibria of the phase separation problem. Notably, buoyancy plays a critical role in driving the ATPS to energy-minimizing crescent-shaped morphologies and shear flows can generate a tenfold speedup in particle formation. Neglecting fluid dynamics produces incorrect minimum-energy droplet shapes. The model allows for optimization of current manufacturing procedures for structured microparticles and improves understanding of ATPS evolution in confined and flowing settings important in biology and biotechnology.