Source author record

Halim Kusumaatmaja

Halim Kusumaatmaja appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

26works
10topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

26 published item(s)

preprint2026arXiv

Stretching and Compressing Capillary Bridges on Hydrophilic, Hydrophobic, and Liquid-infused Surfaces

Aqueous capillary liquid bridges are ubiquitous in nature and in technological processes. Here, we comparatively investigate capillary bridges formed between three distinct types of surfaces: (i) hydrophilic glass, (ii) hydrophobic dichlorodimethylsilane (DMS)-functionalized glass, and (iii) silicone-oil-infused LIS. We combine experimental measurements and computer simulations of the capillary bridge evolution upon changes in the gap size between the surfaces, deriving in each case the bridge geometry and the resulting capillary force. The results, also compared with predictions from the existing theory, follow expected trends on glass and DMS-functionalized surfaces: contact line pinning dominates the bridge behavior on glass with a characteristic stick-slip motion, whereas a pronounced advancing and receding hysteresis is observed on DMS surfaces. On LIS, the absence of pinning leads to minimal force variation, gravity-driven breaking of the bridge symmetry, and possible liquid exchange between LIS through bridge cloaking. These effects become particularly significant in asymmetric bridge configurations combining LIS and DMS surfaces, where the transfer of lubricant from LIS to DMS modifies the effective surface tension and alters bridge-surface interactions. Our systematic comparison of the capillary bridge behavior across solid and liquid interfaces with varying wettability provides a foundation for designing functional surface applications with controlled bridge-surface interactions.

preprint2022arXiv

Spontaneous Phase Separation of Ternary Fluid Mixtures

We computationally study the spontaneous phase separation of ternary fluid mixtures using the lattice Boltzmann method, both when all the surface tensions are equal and when they have different values. Previous theoretical works typically rely on analysing the sign of the eigenvalues resulting from a simple linear stability analysis, but we find this does not explain the fluid morphologies observed. Here, by combining systematic computer simulations over the full range of the composition space and theoretical analysis on the eigenvalues and eigenvectors of the unstable modes, we identify four fundamental phase separation pathways. In particular, we highlight a dominant but so-far overlooked mechanism involving enrichment and instability of the minor component at the fluid-fluid interface

preprint2020arXiv

Axisymmetric flows on the torus geometry

We present a series of analytically solvable axisymmetric flows on the torus geometry. For the single-component flows, we describe the propagation of sound waves for perfect fluids, as well as the viscous damping of shear and longitudinal waves for isothermal and thermal fluids. Unlike the case of planar geometry, the non-uniform curvature on a torus necessitates a distinct spectrum of eigenfrequencies and their corresponding basis functions. This has several interesting consequences, including breaking the degeneracy between even and odd modes, a lack of periodicity even in the flows of perfect fluids and the loss of Galilean invariance for flows with velocity components in the poloidal direction. For the multi-component flows, we study the equilibrium configurations and relaxation dynamics of axisymmetric fluid stripes, described using the Cahn-Hilliard equation. We find a second-order phase transition in the equilibrium location of the stripe as a function of its area $ΔA$. This phase transition leads to a complex dependence of the Laplace pressure on $ΔA$. We also derive the underdamped oscillatory dynamics as the stripes approach equilibrium. Furthermore, relaxing the assumption of axial symmetry, we derive the conditions under which the stripes become unstable. In all cases, the analytical results are confirmed numerically using a finite-difference Navier-Stokes solver.

preprint2020arXiv

Bidirectional Motion of Droplets on Gradient Liquid Infused Surfaces

We demonstrate spontaneous bidirectional motion of droplets on liquid infused surfaces in the presence of a topographical gradient, in which the droplets can move either toward the denser or the sparser solid fraction area. Our analytical theory explains the origin of this bidirectional motion. Furthermore, using both lattice Boltzmann simulations and experiments, we show that the key factor determining the direction of motion is the wettability difference of the droplet on the solid surface and on the lubricant film. The bidirectional motion is shown for various combinations of droplets and lubricants, as well as for different forms of topographical gradients.

preprint2020arXiv

Control of superselectivity by crowding in three-dimensional hosts

Motivated by the fine compositional control observed in membraneless droplet organelles in cells, we investigate how a sharp binding-unbinding transition can occur between multivalent client molecules and receptors embedded in a porous three-dimensional structure. In contrast to similar superselective binding previously observed at surfaces, we have identified that a key effect in a three-dimensional environment is that the presence of inert crowding agents can significantly enhance or even introduce superselectivity. In essence, molecular crowding initially suppresses binding via an entropic penalty, but the clients can then more easily form many bonds simultaneously. We demonstrate the robustness of the superselective behavior with respect to client valency, linker length and binding interactions in Monte Carlo simulations of an archetypal lattice polymer model.

preprint2020arXiv

Critical pressure asymmetry in the enclosed fluid diode

Joint physically and chemically pattered surfaces can provide efficient and passive manipulation of fluid flow. The ability of many of these surfaces to allow only unidirectional flow mean they are often referred to as fluid diodes. Synthetic analogues of these are enabling technologies from sustainable water collection via fog harvesting, to improved wound dressings. One key fluid diode geometry features a pore sandwiched between two absorbent substrates, an important design for applications which require liquid capture while preventing back-flow. However, the enclosed pore is particularly challenging to design as an effective fluid diode, due to the need for both a low Laplace pressure for liquid entering the pore, and a high Laplace pressure to liquid leaving. Here, we calculate the Laplace pressure for fluid travelling in both directions on a range of conical pore designs with a chemical gradient. We show that this chemical gradient is in general required to achieve the largest critical pressure differences between incoming and outgoing liquids. Finally, we discuss the optimisation strategy to maximise this critical pressure asymmetry.

preprint2020arXiv

Factors Controlling the Pinning Force of Liquid Droplets on Liquid Infused Surfaces

Liquid infused surfaces with partially wetting lubricants have recently been exploited for numerous intriguing applications, such as for droplet manipulation, droplet collection and spontaneous motion. When partially wetting lubricants are used, the pinning force is a key factor that can strongly affect droplet mobility. Here, we derive an analytical prediction for contact angle hysteresis {in the limit where the meniscus size is much smaller than the droplet}, and numerically study how it is controlled by the solid fraction, the lubricant wetting angles, and the various fluid surface tensions. We further relate the contact angle hysteresis and the pinning force experienced by a droplet on a liquid infused surface, and our predictions for the critical sliding angles are consistent with existing experimental observations. Finally, we discuss why a droplet on a liquid infused surface with partially wetting lubricants typically experiences stronger pinning compared to a droplet on a classical superhydrophobic surface.

preprint2020arXiv

Modelling double emulsion formation in planar flow-focusing microchannels

Double emulsion formation in a hierarchical flow-focusing channel is systematically investigated using a free energy ternary lattice Boltzmann model. A three dimensional formation regime diagram is constructed based on the capillary numbers of the inner ($Ca_i$), middle ($Ca_m$) and outer ($Ca_o$) phase fluids. The results show that the formation diagram can be classified into periodic two-step region, periodic one-step region, and non-periodic region. By varying $Ca_i$ and $Ca_m$ in the two-step formation region, different morphologies are obtained, including the regular double emulsions, decussate regimes with one or two alternate empty droplets, and structures with multiple inner droplets contained in the continuous middle phase thread. Bidisperse behaviors are also frequently encountered in the two-step formation region. In the periodic one-step formation region, scaling laws are proposed for the double emulsion size and for the size ratio between the inner droplet and the overall double emulsion. Furthermore, we show that the interfacial tension ratio can greatly change the morphologies of the obtained emulsion droplets, and the channel geometry plays an important role in changing the formation regimes and the double emulsion sizes. In particular, narrowing the side inlets or the distance between the two side inlets promotes the conversion from the two-step formation regime to the one-step formation regime.

preprint2020arXiv

Phase transitions on non-uniformly curved surfaces: Coupling between phase and location

For particles confined to two dimensions, any curvature of the surface affects the structural, kinetic and thermodynamic properties of the system. If the curvature is non-uniform, an even richer range of behaviours can emerge. Using a combination of bespoke Monte Carlo, molecular dynamics and basin-hopping methods, we show that the stable states of attractive colloids confined to non-uniformly curved surfaces are distinguished not only by the phase of matter but also by their location on the surface. Consequently, the transitions between these states involve cooperative migration of the entire colloidal assembly. We demonstrate these phenomena on toroidal and sinusoidal surfaces for model colloids with different ranges of interactions as described by the Morse potential. In all cases, the behaviour can be rationalised in terms of three universal considerations: cluster perimeter, stress, and the packing of next-nearest neighbours.

preprint2020arXiv

Self-assembly of small molecules at hydrophobic interfaces using group effect

Although common in nature, the self-assembly of small molecules at sold-liquid interfaces is difficult to control in artificial systems. The high mobility of dissolved small molecules limits their residence at the interface, typically restricting the self-assembly to systems under confinement or with mobile tethers between the molecules and the surface. Small hydrogen-bonding molecules can overcome these issues by exploiting group-effect stabilization to achieve non-tethered self-assembly at hydrophobic interfaces. Significantly, the weak molecular interactions with the solid makes it possible to influence the interfacial hydrogen bond network, potentially creating a wide variety of supramolecular structures. Here we investigate the nanoscale details of water and alcohols mixtures self-assembling at the interface with graphite through group effect. We explore the interplay between inter-molecular and surface interactions by adding small amounts of foreign molecules able to interfere with the hydrogen bond network and systematically varying the length of the alcohol hydrocarbon chain. The resulting supramolecular structures forming at room temperature are then examined using atomic force microscopy with insights from computer simulations. We show that the group-based self-assembly approach investigated here is general and can be reproduced on other substrates such as molybdenum disulphide and graphene oxide, potentially making it relevant for a wide variety of systems.

preprint2020arXiv

The Effect of Ageing on the Structure and Properties of Model Liquid Infused Surfaces

Liquid infused surfaces (LIS) exhibit unique properties that make them ideal candidates for a wide range of applications, from anti-fouling and anti-icing coatings to self-healing surfaces and controlled wetting. However, when exposed to realistic environmental conditions, LIS tend to age and progressively lose their desirable properties, potentially compromising their application. The associated ageing mechanisms are still poorly understood, and results reflecting real-life applications are scarce. Here we track the ageing of model LIS composed of glass surfaces functionalized with hydrophobic nanoparticles and infused with silicone oil. The LIS are fully submerged in aqueous solutions and exposed to acoustic pressure waves for set time intervals. The ageing is monitored by periodic measurements of the LIS' wetting properties. We also track changes to the LIS' nanoscale structure. We find that the LIS rapidly lose their slippery properties due to a combination of oil loss, smoothing of the nanoporous functional layer and substrate degradation when directly exposed to the solution. The oil loss is consistent with water microdroplets entering the oil layer and displacing oil away from the surface. These mechanisms are general and could play a role in the ageing of most LIS.

preprint2019arXiv

Modelling ternary fluids in contact with elastic membranes

We present a thermodynamically consistent model of a ternary fluid interacting with elastic membranes. Following a free-energy modelling approach and taking into account the thermodynamics laws, we derive the equations governing the ternary fluid flow and dynamics of the membranes. We also provide the numerical framework for simulating such fluid-structure interaction problems. It is based on the lattice Boltzmann method, employed for resolving the evolution equations of the ternary fluid in an Eulerian description, coupled to the immersed boundary method, allowing for the membrane equations of motion to be solved in a Lagrangian system. The configuration of an elastic capsule placed at a fluid-fluid interface is considered for validation purposes. Systematic simulations are performed for a detailed comparison with reference numerical results obtained by Surface Evolver, and the Galilean invariance of the proposed model is also proven. The proposed approach is versatile, and a wide range of geometries can be simulated. To demonstrate this, the problem of a capillary bridge formed between two deformable capsules is investigated here.

preprint2016arXiv

Apparent Contact Angle and Contact Angle Hysteresis on Liquid Infused Surfaces

We theoretically investigate the apparent contact angle and contact angle hysteresis of a droplet placed on a liquid infused surface. We show that the apparent contact angle is not uniquely defined by material parameters, but also has a strong dependence on the relative size between the droplet and its surrounding wetting ridge formed by the infusing liquid. We derive a closed form expression for the contact angle in the limit of vanishing wetting ridge, and compute the correction for small but finite ridge, which corresponds to an effective line tension term. We also predict contact angle hysteresis on liquid infused surfaces generated by the pinning of the contact lines by the surface corrugations. Our analytical expressions for both the apparent contact angle and contact angle hysteresis can be interpreted as `weighted sums' between the contact angles of the infusing liquid relative to the droplet and surrounding gas phases, where the weighting coefficients are given by ratios of the fluid surface tensions.

preprint2016arXiv

Kinetic Transition Networks for the Thomson Problem and Smale's 7th Problem

The Thomson Problem, arrangement of identical charges on the surface of a sphere, has found many applications in physics, chemistry and biology. Here we show that the energy landscape of the Thomson Problem for $N$ particles with $N=132, 135, 138, 141, 144, 147$ and $150$ is single funnelled, characteristic of a structure-seeking organisation where the global minimum is easily accessible. Algorithmically constructing starting points close to the global minimum of such a potential with spherical constraints is one of Smale's 18 unsolved problems in mathematics for the 21st century because it is important in the solution of univariate and bivariate random polynomial equations. By analysing the kinetic transition networks, we show that a randomly chosen minimum is in fact always `close' to the global minimum in terms of the number of transition states that separate them, a characteristic of small world networks.

preprint2016arXiv

Phase Separation on Bicontinuous Cubic Membranes: Symmetry Breaking, Re-entrant and Domain Facetting

We study the phase separation of binary lipid mixtures that form bicontinuous cubic phases. The competition between non-uniform Gaussian membrane curvature and line tension leads to a very rich phase diagram, where we observe symmetry breaking of the membrane morphologies and reentrant phenomena due to the formation of bridges between segregated domains. Upon increasing the line tension contribution, we also find facetting of lipid domains that we explain using a simple argument based on the symmetry of the underlying surface and topology.

preprint2016arXiv

Ternary Free Energy Lattice Boltzmann Model with Tunable Surface Tensions and Contact Angles

We present a new ternary free energy lattice Boltzmann model. The distinguishing feature of our model is that we are able to analytically derive and independently vary all fluid-fluid surface tensions and the solid surface contact angles. We carry out a number of benchmark tests: (i) double emulsions and liquid lenses to validate the surface tensions, (ii) ternary fluids in contact with a square well to compare the contact angles against analytical predictions, and (iii) ternary phase separation to verify that the multicomponent fluid dynamics is accurately captured. Additionally we also describe how the model here presented here can be extended to include an arbitrary number of fluid components.

preprint2015arXiv

Free energy pathways of a Multistable Liquid Crystal Device

The planar bistable device [Tsakonas \textit{et al., Appl. Phys. Lett.}, 2007, {\textbf{90}}, 111913] is known to have two distinct classes of stable equilibria: the diagonal and rotated solutions. We model this device within the two-dimensional Landau-de Gennes theory, with a surface potential and without any external fields. We systematically compute a special class of transition pathways, referred to as minimum energy pathways, between the stable equilibria that provide new information about how the equilibria are connected in the Landau-de Gennes free energy landscape. These transition pathways exhibit an intermediate transition state, which is a saddle point of the Landau-de Gennes free energy. We numerically compute the structural details of the transition states, the optimal transition pathways and the free energy barriers between the equilibria, as a function of the surface anchoring strength. For strong anchoring, the transition pathways are mediated by defects whereas we get defect-free transition pathways for moderate and weak anchoring. In the weak anchoring limit, we recover a cusp catastrophe situation for which the rotated state acts as a transition state connecting two different diagonal states.

preprint2015arXiv

Moving contact line dynamics: from diffuse to sharp interfaces

We reconcile two scaling laws that have been proposed in the literature for the slip length associated with a moving contact line in diffuse interface models, by demonstrating each to apply in a different regime of the ratio of the microscopic interfacial width $l$ and the macroscopic diffusive length $l_D= (Mη)^{1/2}$, where $η$ is the fluid viscosity and $M$ the mobility governing intermolecular diffusion. For small $l_D/l$ we find a diffuse interface regime in which the slip length scales as $ξ\sim(l_Dl)^{1/2}$. For larger $l_D/l>1$ we find a sharp interface regime in which the slip length depends only on the diffusive length, $ξ\sim l_D \sim (Mη)^{1/2}$, and therefore only on the macroscopic variables $η$ and $M$, independent of the microscopic interfacial width $l$. We also give evidence that modifying the microscopic interfacial terms in the model's free energy functional appears to affect the value of the slip length only the diffuse interface regime, consistent with the slip length depending only on macroscopic variables in the sharp interface regime. Finally, we demonstrate the dependence of the dynamic contact angle on the capillary number to be in excellent agreement with the theoretical prediction of \cite{Cox1986}, provided we allow the slip length to be rescaled by a dimensionless prefactor. This prefactor appears to converge to unity in the sharp interface limit, but is smaller in the diffuse interface limit. The excellent agreement of results obtained using three independent numerical methods, across several decades of the relevant dimensionless variables, demonstrates our findings to be free of numerical artifacts.

preprint2015arXiv

Surveying the Free Energy Landscapes of Continuum Models: Application to Soft Matter Systems

A variety of methods are developed for characterising the free energy landscapes of continuum, Landau-type free energy models. Using morphologies of lipid vesicles and a multistable liquid crystal device as examples, I show that the methods allow systematic study of not only the most relevant minimum energy configurations, but also the transition pathways between any two minima, as well as their corresponding energy barriers and transition state configurations. A global view of the free energy landscapes can then be visualized using either a disconnectivity graph or a network representation. Different forms of free energy functionals and boundary conditions can be readily implemented, thus allowing these tools to be utilised for a broad range of problems.

preprint2014arXiv

Design principles for Bernal spirals and helices with tunable pitch

Using the framework of potential energy landscape theory, we describe two in silico designs for self-assembling helical colloidal superstructures based upon dipolar dumbbells and Janus-type building blocks, respectively. Helical superstructures with controllable pitch length are obtained using external magnetic field driven assembly of asymmetric dumbbells involving screened electrostatic as well as magnetic dipolar interactions. The pitch of the helix is tuned by modulating the Debye screening length over an experimentally accessible range. The second design is based on building blocks composed of rigidly linked spheres with short-range anisotropic interactions, which are predicted to self-assemble into Bernal spirals. These spirals are quite flexible, and longer helices undergo rearrangements via cooperative, hinge-like moves, in agreement with experiment.

preprint2013arXiv

Defect Motifs for Constant Mean Curvature Surfaces

The energy landscapes of electrostatically charged particles embedded on constant mean curvature surfaces are analysed for a wide range of system size, curvature, and interaction potentials. The surfaces are taken to be rigid, and the basin-hopping method is used to locate the putative global minimum structures. The defect motifs favoured by potential energy agree with experimental observations for colloidal systems: extended defects (scars and pleats) for weakly positive and negative Gaussian curvatures, and isolated defects for strongly negative Gaussian curvatures. Near the phase boundary between these regimes the two motifs are in strong competition, as evidenced from the appearance of distinct funnels in the potential energy landscape. We also report a novel defect motif consisting of pentagon pairs.

preprint2012arXiv

Anisotropic wetting and de-wetting of drops on substrates patterned with polygonal posts

We present results showing how water drops, produced by ink-jet printing, spread on surfaces patterned with lattices of diamond or triangular posts. Considering post widths typically ~7 m and lattice spacings between 15-40 m, we observe drop shapes with 3,4 and 6-fold symmetry, depending on both the symmetry of the lattice and the shape of the posts. This is a result of the different mechanisms of interface pinning and depinning which depend on the direction of the contact line motion with respect to the post shape. Lattice Boltzmann simulations are used to describe these mechanisms in detail for triangular posts. We also follow the motion of the contact line as the drops evaporate showing that they tend to return to their original shape. To explain this we show that the easy direction for movement is the same for spreading and drying drops. We compare the behaviour of small drops with that of larger drops created by jetting several drops at the same position. We find that the contact line motion is unexpectedly insensitive to drop volume, even when a spherical cap of fluid forms above the posts. The findings are relevant to microfluidic applications and to the control of drop shapes in ink-jet printing.

preprint2012arXiv

Non-isomorphic nucleation pathways arising from morphological transitions of liquid channels

Motivated by unexpected morphologies of the emerging liquid phase (channels, bulges, droplets) at the edge of thin, melting alkane terraces, we propose a new heterogeneous nucleation pathway. The competition between bulk and interfacial energies and the boundary conditions determine the growth and shape of the liquid phase at the edge of the solid alkane terraces. Calculations and experiments reveal a "pre-critical" shape transition (channel-to-bulges) of the liquid before reaching its critical volume along a putative shape-conserving path. Bulk liquid emerges from the new shape, and depending on the degree of supersaturation, the new pathway may have two, one, or zero energy barriers. The findings are broadly relevant for many heterogeneous nucleation processes because the novel pathway is induced by common, widespread surface topologies (scratches, steps, etc.).

preprint2012arXiv

Wetting-induced budding of vesicles in contact with several aqueous phases

Osmotic deflation of vesicles enclosing two liquid phases can lead to bulging of one of the phases from the vesicle body. This budding process is preceded by a complete to partial wetting transition of one of the liquid phases on the membrane and depends on the membrane tensions and the tension of the interface between the enclosed liquid phases. These tensions dominate in different morphology regimes, the crossover of which initiates the budding process. In addition, the degree of budding can be controlled by aspiration via micropipettes. We also demonstrate that the budding direction can be reversed if there are two external phases in contact with the vesicle.

preprint2011arXiv

Droplet-Induced Budding Transitions of Membranes

Motivated by recent experiments on biomimetic membranes exposed to several aqueous phases, we theoretically study the morphology of a membrane in contact with a liquid droplet formed via aqueous phase separation. We concentrate on membranes with negligible spontaneous curvature. At small droplet volume, bending energy dominates and the droplet is only partially wrapped by the membrane. At large volume, this configuration can become unstable and undergo a discontinuous transition to a state, in which the droplet is (almost) completely wrapped by the membrane. A morphology diagram, showing the parameter region where such budding transition occurs, is constructed as a function of the membrane tension and the intrinsic contact angle of the liquid with the membrane. The effects of spontaneous curvature are discussed qualitatively.

preprint2010arXiv

Equilibrium morphologies and effective spring constants of capillary bridges

We theoretically study the behaviour of a liquid bridge formed between a pair of rigid and parallel plates. The plates are smooth, they may either be homogeneous or decorated by circular patches of more hydrophilic domains, and they are generally not identical. We calculate the mechanical equilibrium distance of the liquid bridge as a function of liquid volume, contact angle and radius of the chemical domain. We show that a liquid bridge can be an equilibrium configuration as long as the sum of the contact angles at the two walls is larger than $180^\circ$. When comparisons are possible, our results agree well with recent analytical and molecular dynamics simulation results. We also derive the effective spring constant of the liquid bridge as it is perturbed from its equilibrium distance. The spring constant diverges when the sum of the contact angles is $180^\circ$ and is finite otherwise. The value of the spring constant decreases with increasing contact angle and volume, and the rate, at which it decreases, depends strongly on the properties of the two plates.