Local electroneutrality breakdown for electrolytes within varying-section nanopores.

We determine the local charge dynamics of a [Formula: see text] electrolyte embedded in a varying-section channel. By means of an expansion based on the length scale separation between the axial and transverse direction of the channel, we derive closed formulas for the local excess charge for both, dielectric and conducting walls, in 2D (planar geometry) as well as in 3D (cylindrical geometry). Our results show that, even at equilibrium, the local charge electroneutrality is broken whenever the section of the channel is not homogeneous for both dielectric and conducting walls as well as for 2D and 3D channels. Interestingly, even within our expansion, the local excess charge in the fluid can be comparable to the net charge on the walls. We critically discuss the onset of such local electroneutrality breakdown in particular with respect to the correction that it induces on the effective free energy profile experienced by tracer ions.

Enhancement of bubble transport in porous electrodes and catalysts.

We investigate the formation and transport of gas bubbles across a model porous electrode/catalyst using lattice Boltzmann simulations. This approach enables us to systematically examine the influence of a wide range of morphologies, flow velocities, and reaction rates on the efficiency of gas production. By exploring these parameters, we identify critical parameter combinations that significantly contribute to an enhanced yield of gas output. Our simulations reveal the existence of an optimal pore geometry for which the product output is maximized. Intriguingly, we also observe that lower flow velocities improve gas production by leveraging coalescence-induced bubble detachment from the electrode/catalyst.

Cell-free layer development and spatial organization of healthy and rigid red blood cells in a microfluidic bifurcation.

Bifurcations and branches in the microcirculation dramatically affect blood flow as they determine the spatiotemporal organization of red blood cells (RBCs). Such changes in vessel geometries can further influence the formation of a cell-free layer (CFL) close to the vessel walls. Biophysical cell properties, such as their deformability, which is impaired in various diseases, are often thought to impact blood flow and affect the distribution of flowing RBCs. This study investigates the flow behavior of healthy and artificially hardened RBCs in a bifurcating microfluidic T-junction. We determine the RBC distribution across the channel width at multiple positions before and after the bifurcation. Thus, we reveal distinct focusing profiles in the feeding mother channel for rigid and healthy RBCs that dramatically impact the cell organization in the successive daughter channels. Moreover, we experimentally show how the characteristic asymmetric CFLs in the daughter vessels develop along their flow direction. Complimentary numerical simulations indicate that the buildup of the CFL is faster for healthy than for rigid RBCs. Our results provide fundamental knowledge to understand the partitioning of rigid RBC as a model of cells with pathologically impaired deformability in complex in vitro networks.

Closed Formula for Transport across Constrictions.

In the last decade, the Fick-Jacobs approximation has been exploited to capture transport across constrictions. Here, we review the derivation of the Fick-Jacobs equation with particular emphasis on its linear response regime. We show that, for fore-aft symmetric channels, the flux of noninteracting systems is fully captured by its linear response regime. For this case, we derive a very simple formula that captures the correct trends and can be exploited as a simple tool to design experiments or simulations. Lastly, we show that higher-order corrections in the flux may appear for nonsymmetric channels.

Technological Pathways to Produce Compressed and Highly Pure Hydrogen from Solar Power.

Hydrogen (H2 ) produced from renewables will have a growing impact on the global energy dynamics towards sustainable and carbon-neutral standards. The share of green H2 is still too low to meet the net-zero target, while the demand for high-quality hydrogen continues to rise. These factors amplify the need for economically viable H2 generation technologies. The present article aims at evaluating the existing technologies for high-quality H2 production based on solar energy. Technologies such as water electrolysis, photoelectrochemical and solar thermochemical water splitting, liquid metal reactors and plasma conversion utilize solar power directly or indirectly (as carbon-neutral electrons) and are reviewed from the perspective of their current development level, technical limitations and future potential.

Hydrodynamic simulations of sedimenting dilute particle suspensions under repulsive DLVO interactions.

We present guidelines to estimate the effect of electrostatic repulsion in sedimenting dilute particle suspensions. Our results are based on combined Langevin dynamics and lattice Boltzmann simulations for a range of particle radii, Debye lengths and particle concentrations. They show a simple relationship between the slope K of the concentration-dependent sedimentation velocity and the range χ of the electrostatic repulsion normalized by the average particle-particle distance. When χ → 0, the particles are too far away from each other to interact electrostatically and K = 6.55 as predicted by the theory of Batchelor. As χ increases, K likewise increases as if the particle radius increased in proportion to χ up to a maximum around χ = 0.4. Over the range χ = 0.4-1, K relaxes exponentially to a concentration-dependent constant consistent with known results for ordered particle distributions. Meanwhile the radial distribution function transitions from a disordered gas-like to a liquid-like form. Power law fits to the concentration-dependent sedimentation velocity similarly yield a simple master curve for the exponent as a function of χ, with a step-like transition from 1 to 1/3 centered around χ = 0.6.

Transport of neutral and charged nanorods across varying-section channels.

We study the dynamics of neutral and charged rods embedded in varying-section channels. By means of systematic approximations, we derive the dependence of the local diffusion coefficient on both the geometry and charge of the rods. This microscopic insight allows us to provide predictions for the permeability of varying-section channels to rods with diverse lengths, aspect ratios and charge. Our analysis shows that the dynamics of charged rods is sensitive to the geometry of the channel and that their transport can be controlled by tuning both the shape of the confining walls and the charge of the rod. Interestingly, we find that the channel permeability does not depend monotonically on the charge of the rod. This opens the possibility of a novel mechanism to separate charged rods.

Regimes of motion of magnetocapillary swimmers.

The dynamics of a triangular magnetocapillary swimmer is studied using the lattice Boltzmann method. We extend on our previous work, which deals with the self-assembly and a specific type of the swimmer motion characterized by the swimmer's maximum velocity centred around the particle's inverse viscous time. Here, we identify additional regimes of motion. First, modifying the ratio of surface tension and magnetic forces allows to study the swimmer propagation in the regime of significantly lower frequencies mainly defined by the strength of the magnetocapillary potential. Second, introducing a constant magnetic contribution in each of the particles in addition to their magnetic moment induced by external fields leads to another regime characterized by strong in-plane swimmer reorientations that resemble experimental observations.

Two-dimensional Cahn-Hilliard simulations for coarsening kinetics of spinodal decomposition in binary mixtures.

The evolution of the microstructure due to spinodal decomposition in phase separated mixtures has a strong impact on the final material properties. In the late stage of coarsening, the system is characterized by the growth of a single characteristic length scale L ∼ Ctα. To understand the structure-property relationship, the knowledge of the coarsening exponent α and the coarsening rate constant C is mandatory. Since the existing literature is not entirely consistent, we perform phase field simulations based on the Cahn-Hilliard equation. We restrict ourselves to binary mixtures using a symmetric Flory-Huggins free energy and a constant composition-independent mobility term and show that the coarsening for off-critical mixtures is slower than the expected t1/3-growth. Instead, we find α to be dependent on the mixture composition and associate this with the observed morphologies. Finally, we propose a model to describe the complete coarsening kinetics including the rate constant C.

Self-Similar Liquid Lens Coalescence.

A basic feature of liquid drops is that they can merge upon contact to form a larger drop. In spite of its importance to various applications, drop coalescence on prewetted substrates has received little attention. Here, we experimentally and theoretically reveal the dynamics of drop coalescence on a thick layer of a low viscosity liquid. It is shown that these so-called "liquid lenses" merge by the self-similar vertical growth of a bridge connecting the two lenses. Using a slender analysis, we derive similarity solutions corresponding to the viscous and inertial limits. Excellent agreement is found with the experiments without any adjustable parameters, capturing both the spatial and temporal structures of the flow during coalescence. Finally, we consider the crossover between the two regimes and show that all data of different lens viscosities collapse on a single curve capturing the full range of the coalescence dynamics.

Desorption energy of soft particles from a fluid interface.

The efficiency of soft particles to stabilize emulsions is examined by measuring their desorption free energy, i.e., the mechanical work required to detach the particle from a fluid interface. Here, we consider rubber-like elastic as well as microgel particles, using coarse-grained molecular dynamics simulations. The energy of desorption is computed for two and three-dimensional configurations by means of the mean thermodynamic integration method. It is shown that the softness affects the particle-interface binding in two opposing directions as compared to rigid particles. On the one hand, a soft particle spreads at the interface and thereby removes a larger unfavorable liquid-liquid contact area compared to rigid particles. On the other hand, softness provides the particle with an additional degree of freedom to get reshaped instead of deforming the interface, resulting in a smaller restoring force during the detachment. It is shown that the first effect prevails so that a soft spherical particle attaches to the fluid interface more strongly than rigid spheres. Finally, we consider microgel particles both in the swollen and in the collapsed state. Surprisingly, we find that the latter has a larger binding energy. All results are rationalised using thermodynamic arguments and thereby offer detailed insights into the desorption energy of soft particles from fluid interfaces.

Capillary interactions between soft capsules protruding through thin fluid films.

When a suspension dries, the suspending fluid evaporates, leaving behind a dry film composed of the suspended particles. During the final stages of drying, the height of the fluid film on the substrate drops below the particle size, inducing local interface deformations that lead to strong capillary interactions among the particles. Although capillary interactions between rigid particles are well studied, much is still to be understood about the behaviour of soft particles and the role of their softness during the final stages of film drying. Here, we use our recently-introduced numerical method that couples a fluid described using the lattice Boltzmann approach to a finite element description of deformable objects to investigate the drying process of a film with suspended soft particles. Our measured menisci deformations and lateral capillary forces, which agree well with previous theoretical and experimental works in case of rigid particles, show that the deformations become smaller with increasing particle softness, resulting in weaker lateral interaction forces. At large interparticle distances, the force approaches that of rigid particles. Finally, we investigate the time dependent formation of particle clusters at the late stages of the film drying.

A phase-field model for the evaporation of thin film mixtures.

The performance of solution-processed solar cells strongly depends on the geometrical structure and roughness of the photovoltaic layers formed during film drying. During the drying process, the interplay of crystallization and liquid-liquid demixing leads to structure formation on the nano- and microscale and to the final rough film. In order to better understand how the film structure can be improved by process engineering, we aim at theoretically investigating these systems by means of phase-field simulations. We introduce an evaporation model based on the Cahn-Hilliard equation for the evolution of the fluid concentrations coupled to the Allen-Cahn equation for the liquid-vapour phase transformation. We demonstrate its ability to match the experimentally measured drying kinetics and study the impact of the parameters of our model. Furthermore, the evaporation of solvent blends and solvent-vapour annealing are investigated. The dry film roughness emerges naturally from our set of equations, as illustrated through preliminary simulations of spinodal decomposition and film drying on structured substrates.

The effect of the liquid layer thickness on the dissolution of immersed surface droplets.

Droplets on a liquid-immersed solid surface are key elements in many applications, such as high-throughput chemical analysis and droplet-templated porous materials. Such surface droplets dissolve when the surrounding liquid is undersaturated and the dissolution process is usually treated analogous to a sessile droplet evaporating in air. Typically, theoretical models predict the mass loss rate of dissolving droplets as a function of droplet geometrical factors (radius, constant angle), and droplet material properties (diffusion constant and densities), where the thickness of the surrounding liquid layer is neglected. Here, we investigate, both numerically and theoretically, the effect of the liquid layer thickness on the dissolution of surface droplets. We perform 3D lattice Boltzmann simulations and obtain the density distribution and time evolution of droplet height during dissolution. Moreover, we find that the dissolution slows down and the lifetime linearly increases with increasing the liquid layer thickness. We propose a theoretical model based on a quasistatic diffusion equation which agrees quantitatively with simulation results for thick liquid layers. Our results offer insight to the fundamental understanding of dissolving surface droplets and can provide valuable guidelines for the design of devices where the droplet lifetime is of importance.

Correction: Capillary assemblies in a rotating magnetic field.

Correction for 'Capillary assemblies in a rotating magnetic field' by Galien Grosjean et al., Soft Matter, 2019, DOI: .

Capillary assemblies in a rotating magnetic field.

Small objects floating on a fluid have a tendency to aggregate due to capillary forces. This effect has been used, with the help of a magnetic induction field, to assemble submillimeter metallic spheres into a variety of structures, whose shape and size can be tuned. Under time-varying fields, these assemblies can propel themselves due to a breaking of time reversal symmetry in their adopted shapes. In this article, we study the influence of an in-plane rotation of the magnetic field on these structures. Various rotational modes have been observed with different underlying mechanisms. The magnetic properties of the particles cause them to rotate individually. Dipole-dipole interactions in the assembly can cause the whole structure to align with the field. Finally, non-reciprocal deformations can power the rotation of the assembly. Symmetry plays an important role in the dynamics, as well as the frequency and amplitude of the applied field. Understanding the interplay of these effects is essential, both to explain previous observations and to develop new functions for these assemblies.

How antagonistic salts cause nematic ordering and behave like diblock copolymers.

We present simulation results and an explanatory theory on how antagonistic salts affect the spinodal decomposition of binary fluid mixtures. We find that spinodal decomposition is arrested and complex structures form only when electrostatic ion-ion interactions are small. In this case, the fluid and ion concentrations couple and the charge field can be approximated as a polynomial function of the relative fluid concentrations alone. When the solvation energy associated with transferring an ion from one fluid phase to the other is of the order of a few kBT, the coupled fluid and charge fields evolve according to the Ohta-Kawasaki free energy functional. This allows us to accurately predict structure sizes and reduce the parameter space to two dimensionless numbers. The lamellar structures induced by the presence of the antagonistic salt in our simulations exhibit a high degree of nematic ordering and the growth of ordered domains over time follows a power law. This power law carries a time exponent proportional to the salt concentration. We qualitatively reproduce and interpret neutron scattering data from previous experiments of similar systems. The dissolution of structures at high salt concentrations observed in these experiments agrees with our simulations, and we explain it as the result of a vanishing surface tension due to electrostatic contributions. We conclude by presenting 3D results showing the same morphologies as predicted by the Ohta-Kawasaki model as a function of volume fraction and suggesting that our findings from 2D systems remain valid in 3D.

Optimal motion of triangular magnetocapillary swimmers.

A system of ferromagnetic particles trapped at a liquid-liquid interface and subjected to a set of magnetic fields (magnetocapillary swimmers) is studied numerically using a hybrid method combining the pseudopotential lattice Boltzmann method and the discrete element method. After investigating the equilibrium properties of a single, two, and three particles at the interface, we demonstrate a controlled motion of the swimmer formed by three particles. It shows a sharp dependence of the average center-of-mass speed on the frequency of the time-dependent external magnetic field. Inspired by experiments on magnetocapillary microswimmers, we interpret the obtained maxima of the swimmer speed by the optimal frequency centered around the characteristic relaxation time of a spherical particle. It is also shown that the frequency corresponding to the maximum speed grows and the maximum average speed decreases with increasing interparticle distances at moderate swimmer sizes. The findings of our lattice Boltzmann simulations are supported by bead-spring model calculations.

Blood Crystal: Emergent Order of Red Blood Cells Under Wall-Confined Shear Flow.

Driven or active suspensions can display fascinating collective behavior, where coherent motions or structures arise on a scale much larger than that of the constituent particles. Here, we report numerical simulations and an analytical model revealing that deformable particles and, in particular, red blood cells (RBCs) assemble into regular patterns in a confined shear flow. The pattern wavelength concurs well with our experimental observations. The order is of a pure hydrodynamic and inertialess origin, and it emerges from a subtle interplay between (i) hydrodynamic repulsion by the bounding walls that drives deformable cells towards the channel midplane and (ii) intercellular hydrodynamic interactions that can be attractive or repulsive depending on cell-cell separation. Various crystal-like structures arise depending on the RBC concentration and confinement. Hardened RBCs in experiments and rigid particles in simulations remain disordered under the same conditions where deformable RBCs form regular patterns, highlighting the intimate link between particle deformability and the emergence of order.

From Dot to Ring: The Role of Friction in the Deposition Pattern of a Drying Colloidal Suspension Droplet.

The deposition of particles on a substrate by drying a colloidal suspension droplet is at the core of applications ranging from traditional printing on paper to printable electronics or photovoltaic devices. The self-pinning induced by the accumulation of particles at the contact line plays an important role in the formation of a deposit. In this article, we investigate, both numerically and theoretically, the effect of friction between the particles and the substrate on the deposition pattern. Without friction, the contact line shows a stick-slip behavior and a dotlike deposit is left after the droplet is evaporated. By increasing the friction force, we observe a transition from a dotlike to a ringlike deposit. We propose a theoretical model to predict the effective radius of the particle deposit as a function of the friction force. Our theoretical model predicts a critical friction force when self-pinning happens and the effective radius of deposit increases with increasing friction force, confirmed by our simulation results. Our results can find implications for developing active control strategies for the deposition of drying droplets.