Directional Fluidity of Dense Emulsion Activated by Transverse Wedge-Shaped Microroughness

TL;DR

This paper shows how wedge-shaped microroughness can enhance the flow of emulsions in microfluidic channels by activating directional fluidization.

Contribution

The study introduces wedge-shaped microroughness as a novel method to control and enhance emulsion flow in microfluidic systems.

Findings

Wedge-shaped microroughness increases emulsion flow in the direction of the climbing ramp.

A 30% gain in volumetric flow rate is observed compared to the opposite direction.

The effect depends on the magnitude of the applied pressure gradient.

Abstract

The handling and fluidization of amorphous soft solids, such as emulsions, foams, or gels, is crucial in many technological processes. This is generally achieved by applying mechanical stress that overcomes a critical threshold, known as yield stress, below which these systems behave as elastic solids. However, the interaction with the walls can facilitate the transition from solid to fluid by activating rearrangements of the fluid constituents close to the wall, resulting in increased fluidity of the system up to distances greater than the spatial scale of the rearrangements. We address the impact of wedge-shaped microroughness on activating the fluidization of emulsion droplets in pressure-driven flow through microfluidic channels. We realize the micro wedges by maskless photolithography to texture one wall of the channel and measure the velocity profiles for flow directed accordingly and against the increasing ramp of the wedge-shaped grooves. We report the enhancement of the emulsion flow in the direction of the climbing ramp of the wedge activated by increasing the magnitude of the pressure gradient. A gain for the volumetric flow rate is registered with respect to the opposite direction as being to 30%, depending on the pressure drop.

Figures

Full text

1. Introduction

The flow of dense emulsions and pastes is essential for their processing and technological use [1,2,3,4,5]. They behave as elastic solids at low stresses and flow only when subjected to stresses exceeding a critical value, known as yield stress Their ability to flow depends on many factors, including rheology [6,7,8,9], the confinement within the capillary or microfluidic channels [10], and the interaction with the walls of the channels [11,12,13,14,15,16,17]. While the roughness of the channel has been extensively considered to increase the fluidity of dense emulsions, directional effects introduced by specific geometries are still poorly explored. It has been shown that the introduction of herringbone-shaped microroughness in pressure-driven flow introduces an activation of emulsion flow in the form of the banding of the flow in specific regions of the channel cross section, the activation being at a maximum on the tip of the herringbone grooves, when the emulsion flows in the direction of grooves narrowing [18]. Recently, 2D numerical simulations on triangular-shaped roughness reported similar activation when the flow is directed to the upward slope of the wedge-shaped obstacles, compared to the direction oriented against the vertical wall of the wedge [19]. However, no experimental confirmation has been presented so far for this wedge-like geometry.

This work aims to fill this gap by studying the velocity profiles and stress distribution of dense emulsions within a microfluidic channel, the bottom of which is textured by a periodic arrangement of wedge-shaped microgrooves, aligned perpendicularly to the flow direction. To address the role of increasing or decreasing wedges, the flow direction is alternatively inverted.

In standard (mask-based) photolithography, fabrication of gradient structures, such as modular wedges or pyramids, requires multiple masks and iterations, which is challenging, expensive, and time consuming [20,21,22].

In the last decade, maskless photolithography has progressively been established based on the direct projection of light onto the photoresistor using a computer-controlled Spatial Light Modulator (SLM) in place of the traditional, high-resolution photomask [23]. Grayscale photolithography has been widely adopted in the manufacturing of microelectromechanical systems (MEMSs), microlens arrays, Fresnel lenses, and molds for soft lithography [24,25,26]. Historically, grayscale photolithography required complex sets of physical hard masks that were applied sequentially to expose specific regions to different light levels. However, this approach had several limitations, including low spatial resolution, high costs, long fabrication times, poor process control, and a lack of flexibility, all of which hindered its broader adoption [27,28]. In contrast, digital designs allow for the precise control of light exposure according to the desired dose and target regions. This method enables faster prototyping, easy design modification, and precise control over grayscale variations to produce gradients. A digital mask allows continuous redesign and reexposure, while grayscale modulation enables varying doses for smooth transitions in structure height. This innovative technique is highly precise, versatile, and efficient, offering superior control and smoother microstructures, making it ideal for advanced microfabrication tasks [29,30,31,32]. Compared to traditional exposure through photomasks, maskless photolithography significantly simplifies and accelerates the fabrication of structures with height gradients.

The manuscript is organized as follows. In Section 2, we describe the procedure to realize the wedge-shaped microroughness by using maskless photolithography, preparing the emulsions, and controlling the emulsion flow. In Section 3, we address the fluid dynamic behavior of the textured channel by reporting the velocity profiles of the emulsion flow at different pressure loads in opposite directions with respect to the wedge roughness. Moreover, we compare the velocity profiles and the volumetric flow rates in textured and flat channels. In Section 4, we rationalize the results and review some relevant perspectives. In Appendix A, we show a comprehensive characterization of the realized structures, while Appendix B describes the experimental setup used to measure the velocity profiles.

2. Materials and Methods

Microfabrication. The microfluidic textures are fabricated using maskless photolithography with positive liquid MAP-1275G photoresist (Micro Resist Technology GmbH, Berlin, Germany) spun on a 75 × 25 mm^2^ glass slide [33]. The photolithographic setup is the Tabletop Micro Maskless Aligner $[eqn]$ -MLA (Heidelberg Instrument, Heidelberg, Germany), which uses a 365 nm LED of 6 W power to engrave programmable images on the photoresistor through a Spatial Light Modulator [34,35]. The mask is a digital drawing that can be easily redesigned and reexposed as needed. In addition, the direct-write process allows the laser to be modulated according to grayscale variations, as shown in Figure 1. Each grayscale value in the mask corresponds to a different dose, enabling the automatic representation of gradients [36].

Starting from the spin-coating process, 1 mL of the photoresistor is deposited on the substrate, which is spun at 3000 rpm for a uniform layer thickness $[eqn]$ of approximately 7.5 μm. Subsequently, the substrate is baked at 95 °C for 1 min to harden the resistor and improve adhesion. The selected digital mask is exposed to a dose of 180 mJ/cm^2^ using μ-MLA. The process is automatically repeated to cover an area 1 cm in width and 5 cm in length. The controlled roughness is applied exclusively to one wall of the channel, which is textured with wedges that span the entire channel width and are evenly distributed along its length. The specifications for the microfabrication process of wedges and their characteristics can be found in Appendix A. Finally, the exposed substrate is developed for 90 s within the developer solution ma-D 532/S (Micro Resist Technology GmbH, Germany).

The channel was then assembled using as the top wall a microscope slide previously drilled to apply the tubes for inlets and outlets. To close the channel, an H = 120 μm thick dry photoresist film (WBR-2000, DuPont, Wilmington, DE, USA) was used, specifically leveraging its excellent adhesion properties to surfaces [37,38]. The dry photoresist acts both as a glue to bond the top and bottom glass slides and as a sidewall of constant height. The dry photoresist was cut into a channel shape using a knife plotter (Craft Robo CC200-20, Graphtec, Yokohama, Japan) and then laminated to the glass plate, removing trapped air bubbles. To complete the bond procedure, the sample was placed under a hydraulic press for 5 min at 80 °C with a pressure of 2 bar. With this step, the cover plate was pre-bonded to the cover glass. Finally, to fix the different layers and polymerize the photoresist, the sample was exposed to UV light (i-line 365 nm) for 60 s. We characterized these microstructures using a Tencor P-17 Stylus Profiler from KLA Corporation (Milpitas, CA, USA).

The sequence of wedge-shaped textures is arranged perpendicularly to the x direction of the flow and has heights of 4.5 μm, somewhat similar to the mean droplet size of the emulsion. A hydrophilic coating was applied to the inner walls of the microfluidic channel to homogenize the wetting properties of the surfaces and to avoid the adhesion and coalescence of oil droplets on the walls. The coating was achieved by treating the closed microfluidic chip with an oxygen plasma to activate the inner surfaces for 1 min at 300 W and an oxygen flow of 35 sccm ( $[eqn]$ ). The plasma treatment was performed using the ‘Femto’ plasma system from Diener Electronic, GmbH, Germany. Then, a 2.5% (w/w) PVP (Polyvinylpyrrolidone K90, AppliChem, GmbH, Darmstadt, Germany) water solution flowed within the channel for at least 2 h with a flow rate of 1.5 $[eqn]$ using a syringe pump [39] After treatment, the channel was rinsed by flushing 5 mL of MilliQ water (Millipore Sigma, Burlington, MA, USA) at a flow rate of 2.5 $[eqn]$ to remove excess PVP and dried with nitrogen flow at a pressure of 1 bar.

Emulsion preparation. The concentrated emulsion is prepared by dispersing silicone oil (polydimethyl siloxane, 1000 mPas) in a glycerine/water mixture (54% w/w) stabilized with 1% w/w tetradecyl trimethylammonium bromide (TTAB). The surfactant concentration is set to be high enough to prevent droplet coalescence but low enough to avoid flocculation by depletion [10]. The continuous phase also contains rhodamine B-labeled dye microparticles of 0.2 μm (FluoSpheres, ThermoFisher, Waltham, MA, USA), which act as fluorescent tracers. The volume fraction occupied by the oil phase relative to the total volume of the emulsion is set to $[eqn]$ . This value ensures that the emulsion is stable for months and unyielding so that the roughness-induced fluidization can be addressed [16,18,40]. At the same time, it is low enough not to require a high stress, which may cause losses in the microfluidic channel.

The emulsion is prepared using a custom power homogenizer consisting of a commercial drill operating at a stabilized 24 V DC voltage and specific blades adapted to a commercial mortar (model Z247472, Merck Group, Darmstadt, Germany), which were designed and 3D-printed for this purpose. The continuous phase is placed in the mortar, and the oil phase is added at a rate of approximately 100 μL/min using a syringe pump (11 Elite Module Pico Single, Harvard Apparatus, Holliston, MA, USA). The mixture is subjected to constant and intense shearing until the oil volume fraction reaches the desired value. The voltage applied is 11.2 ± 0.1 V, corresponding to a mixing rate of 30 Hz.

To measure the droplet size distribution, the emulsion was diluted in an aqueous solution containing 1% TTAB, which helps to maintain the silicone oil in the dispersed phase. A few drops of this solution were placed between two microscope coverslips. The images were then captured using a bright-field optical microscope with a 100× magnification objective, as shown in Figure 2. The mean particle diameter of the emulsion is $[eqn]$ = 4 μm, with a coefficient of variation (CV) of 33.5%, derived from a standard deviation ( $[eqn]$ ) of 1.34 μm. The polydispersity index (PDI) is 1.5.

Flow Control. The emulsion flow was pressure-controlled using a Microfluidic Flow Control System (MFCS series, Fluigent, Paris, France). Stable pressure loads in the range of (1000–7000) ± 5 mbar were applied between the inlet and the outlet of the microfluidic channel. The pressure difference was alternated to reverse the flow direction along the inclined side of the wedge-shaped grooves or in the opposite direction, referred to as forward and backward, respectively (see Figure 3a,b).

Velocimetry. The velocity profiles of the emulsion flow were measured by tracking fluorescent microtracers within a 36 × 14 μm^2^ region of interest (see Figure A2c) [40,41]. The microtracers were illuminated using a 532 nm diode-pumped solid-state (DPSS) laser (Roithner Lasertechnik GmbH, Vienna, Austria) and imaged through an inverted motorized microscope (Eclipse Ti-E, Nikon, Tokyo, Japan) in an epifluorescence configuration. The detection system consisted of a scientific CMOS camera (ORCA Fusion, Hamamatsu Photonics, Hamamatsu, Japan). A scheme of the entire experimental setup is provided in Appendix B. Images were acquired as z-stacks with 6 μm steps, and only trajectories crossing the optical field within the depth of focus were considered. The velocity $[eqn]$ was determined by averaging tracer displacements at the channel center. Data analysis was performed using a custom routine based on the TrackMate 7 algorithm [42,43,44].

Geometry of the experiment. The grooves are designed so that the height h of the wedges is comparable to the mean droplet size $[eqn]$ , while their length l and the groove separation g are about 10 times greater, with $[eqn]$ . All the main dimensions of the channel are summarized in Table 1. In particular, $[eqn]$ and $[eqn]$ . This configuration results in a tiny directional parameter, characterized by a smooth tilt of the wedges, with trains of several droplets free to travel over the flat wall before encountering a wedge-shaped groove (see Figure 3b).

3. Results

In Figure 4, we report the velocity profiles of the emulsion at the volume fraction $[eqn]$ = 0.76 flowing in the patterned channel at different pressure drops. Under low-pressure drop conditions, the emulsion is flowing at constant velocity as a rigid elastic solid (plug flow) regardless of the flow direction.

As the pressure drop increases, the emulsion becomes progressively more fluidized, as reported by an increase in the velocity with the distance from the channel walls (shear flow). However, this increase is not the same for both forward and backward flow directions, while the magnitude of applied stress is unchanged. The effect of directionality emerges, leading to a noticeable gap between the forward and backward profiles. This gap widens as the pressure gradient increases, reflecting the growing impact of the directional roughness on the emulsion’s behavior. To evaluate the effectiveness of the wedge-shaped microroughness in enhancing the velocity profile, we systematically examined the emulsion flow in a microfluidic channel featuring two flat (unstructured) walls, at the same pressure drops applied in the textured channel. Both channels have the same height H = 120 μm, so that the magnitude of the pressure gradient is the same for both channels. This comparison highlights the role of directional roughness in the fluidization of the emulsion. As expected, in the flat–flat channel, the velocity profiles on forward and backward flows are the same within the experimental error bars (see Figure 5). In parallel, the presence of the wall texturing promotes more fluidization than the flat channel, regardless of the specific direction of the emulsion flow. By integrating the velocity profile over the channel cross section Σ, we obtain the volumetric flow rate $[eqn]$ . For the textured channel, we distinguish between the flow rates $[eqn]$ and $[eqn]$ , taken in the forward and backward directions, respectively. In the flat channel, the flow rate $[eqn]$ is obtained by averaging the velocity profiles in both directions. Clearly, the flow rate increases as the pressure drop increases. For the flat channel, it is the same within the error bars, regardless of the flow direction, as expected. In contrast, for the textured channel, the flow rate does not increase with pressure by the same amount in both directions, the increase being more consistent in the forward direction (see Figure 6a). To quantify the efficiency of the roughness design in increasing the emulsion transport at constant applied stress, we introduce the flow gain $[eqn]$ as

$$[eqn]$$

With $[eqn]$ being the flow gains of the forward (backward) flow with respect to the flat channel, respectively; $[eqn]$ being the flow gain of the forward flow with respect to backward, provided that all flows are compared at the same pressure load with that different pressure loads (Figure 6b).

As shown in Figure 6b, when the emulsion flows in the direction of the raising ramp (forward) of the wedge-shaped grooves, the flow rate gains about $[eqn]$ compared to the opposite direction, for all pressure drops explored but the lowest. Indeed, for $[eqn]$ P = 800 mbar, the emulsion is only at the beginning of the roughness fluidization process. Conversely, with respect to the flat channel, the gain in the forward direction is maximum and increases until about $[eqn]$ as the pressure drop increases, while it is similar to that computed considering forward and backward directions and is rather independent of the pressure, except for the lowest one. To investigate the stress profile within the channel, we make use of a specific protocol developed and validated by numerical simulations [19]. In fact, due to the presence of different boundary conditions at the channel walls (flat and rough), the stress profile within the channel is not directly accessible by simply knowing the pressure drop and the length of the channel. The protocol is rooted in the calibration between the slip velocity, i.e., the velocity of the emulsion flow at the contact with the solid wall and the shear stress, within a channel featuring both walls flat, and the same height of the flat–rough channel. The corner stone of the protocol consists of the validation that allows assigning the stress at the flat wall in a flat–rough microfluidic channel once the relation between stress and slip velocity in a flat–flat channel is measured for the same emulsion (see Figure 7a).

The mechanical-balance condition $[eqn]$ implies a linear stress profile $[eqn]$ as a function of the vertical position z, where the slope is given by the pressure gradient $[eqn]$ . The stress profiles are different for the opposite flow directions, corresponding to different sides of the wedged grooves encountered by the emulsion droplets (see Figure 7b). Notably, this implies that the rheological response of the emulsion is directional in the sense that it is different for the different configurations of the wedge-shaped roughness, here obtained by reversing the emulsion flow (see Figure 7c). This result is similar to that reported either for channels textured by herringbone-shaped roughness or numerical wedge-shaped roughness, although the details of the experimental wedge-grooves do not match those of the simulations [18,19]. Furthermore, this aligns with the report of a different rheological response when boundary wettability is changed and an increase in shear stress is observed for a hydrophobic wall compared to a hydrophilic one [12,15]. In this case, that is, in the presence of wedge-shaped grooves, the rheological response is not determined by a change in the wettability pattern but is due to a change in flow direction with respect to a fixed suitable physical texture of the wall roughness.

4. Discussion and Conclusions

From a micromechanical perspective, the roughness-induced fluidization of dense emulsion is attributed to an increase in the number of plastic rearrangements [10,16,45]. Droplets interacting with the textured wall undergo more variations in speed and direction compared to the flat wall, resulting in a scrambling effect that promotes interactions with neighboring droplets and facilitates flow. This is marked by some key changes in the velocity profiles, such as the increase in the velocity, the growth in the shear flow regions close to walls, and the corresponding thinning of the plug flow region in the center of the channel height, as reported by Figure 4 and Figure 5. The directional design has shown that the roughness geometry plays a key role. For herringbone-shaped roughness, a flow banding along the channel cross section has been reported, with flow boost at the herringbone tip [18]. For wedge-shaped roughness, numerical simulations have elucidated the lift-off of emulsion droplets climbing on the wedges as the driving mechanism that promotes directional fluidization [19]. Here, we confirm this view by showing a significant gain in the volumetric flow rate along the direction that promotes the droplets’ lift-off, compared to the opposite direction or the flat wall. The directional fluidization takes place at a high enough droplet concentration and pressure. This is consistent with the reporting within different geometries for the directional roughness [18]. In fact, at a low volume fraction, the emulsion is fluid-like for vanishing pressure drops. Therefore, it is rather insensitive to the specific wall texturing. Conversely, for dense emulsions, the pressure load should overcome a threshold to trigger fluidization and, even more, to activate directional fluidization. Indeed, plug (or almost plug) flows are independent of or less sensitive to directional details of the roughness texturing. Our results confirm that the gain in the flow rate increases with increasing pressure. In the limit of a very high pressure load, a recovery to a fluid system is expected and also reported [18]. However, in the present study, we could only address a limited range of pressure drop since $[eqn]$ mbar is the highest pressure before the channels start leaking.

Concluding, to investigate the directional fluidization of emulsions flowing through a microfluidic channel with one wall textured by wedge-shaped grooves, we successfully employed grayscale maskless photolithography to fabricate structures with a height gradient. The optimization of the exposure parameters for the photoresist enabled the fabrication of grooves with a height of approximately 4 μm, demonstrating the precision and effectiveness of this process for microfluidic applications. At the fluid-dynamic level, the wedge-shaped design promotes a directional fluidization of dense emulsion, being more effective in the direction that follows the rising side of the wedge. By running a protocol specifically validated to measure the stress profile in channels with different boundary conditions on opposite walls, we obtained the stress distribution in the channel and reported the presence of two flow curves $[eqn]$ according to the two flow directions, forward and backward with respect to the wedge side [19]. This reflects the selective increase or decrease in the emulsion flow rate promoted by directional roughness, regardless of the specific geometry considered. Considering different combinations $[eqn]$ of the wedge-shaped microgrooves, it is natural to ask which one would produce the optimal flow throughput in a microfluidic channel at a given pressure gradient. The resulting space of parameters would be very demanding for practical realization, and the help of numerical simulations would be invaluable to individuate the optimal set.

The work presented here has significant implications for several emerging processes such as bioprinting, where emulsions are used as scaffold structures to support cell growth in efforts to replicate living tissues [46,47,48]. While microfluidics enables precise control over fluid flows, the printability of various bioinks is often compromised by the increased stress during extrusion, which can damage cells and negatively affect print quality. Our findings reveal a notable improvement in fluid flow under the same pressure gradient when compared to unstructured channels. This enhancement holds promise for bioprinting applications as it enables higher fluid throughput. Consequently, this reduces the time cells are exposed to damaging forces, such as those within the nozzle. By minimizing this exposure, we can improve cell viability and ensure better material quality post-printing [49]. Furthermore, incorporating directional roughness into the printing channels could further optimize fluid flow and cell arrangement.