STRUCTURED PDMS USED AS ACTIVE ELEMENT FOR A BIOMIMETICS INSPIRED FLUID TRANSPORTER

Overview

A pneumatically actuated fluid transporter to transport fluids along surfaces is introduced. The biomimetic approach is based on the transportation principle of fluids by cilia or comb-row arrays due to the generation of metachronal waves. Rows of PDMS flaps which mimic comb-rows are asymmetrically positioned on flexible membranes. Each membrane is deflected by applying a defined pressure profile to achieve a metachronal wave on the surface. The simulations of the membrane behavior as well as a description of the concept by applying metachronal waves to the artificial comb row arrays are presented. The proof-of-concept shows fluid transport of up to 64 µm/s near the flap tips.

Keywords:
PDMS, microfluidics, fluid transport, particle transport, biomimetics, cilia

Introduction

In microfluidic devices a couple of methods exist to transport particles in closed channels, e.g. pressure driven flow [1], electrophoresis [2], dielectrophoresis [3], or surface acoustic wave devices [4]. Particle or fluid transportation in open channels may be required in some cases. Here, the above described methods may fail due to the large distance or even absence of a cover plate.

In nature, many ciliated surfaces are known, e.g. in the respiratory tract or in the fallopian tube. Cilia are also used as propelling mechanism for microorganisms, e.g. Mesodinium rubrum, Paramecium, Volvox, Chlamydomonas reinhardtii, Ctenophore. These cilia arrays often move in a synchronized manner with a certain phase shift between adjoined cilia to optimize the transportation or self-propulsion. The biological actuation principle can vary over the different species.

In the biomimetic approach to mimic biological cilia, different techniques of forming artificial cilia and methods to stimulate those were introduced in the last decade, such as electrostatic forces [5, 6], magnetic field [7–11], light-excitation [12], chemical [13–15], or mechanical forces [16]. However, the fabrication processses of all these approaches suffer from high costs and complexity. On the other hand, soft-lithography and micro molding of polydimethyl-siloxane (PDMS) are established as techniques for the production of low-cost microfluidic devices. They can ideally be used for fabrication of flexible membranes and flaps.

The objective of this paper is the investigation of the propelling mechanism of a pneumatically working artificial comb row, which can be easily fabricated using standard low-cost soft-lithography and micromolding techniques. The concept is shown in Figure 1 and is already published in [17]. Rows of cilia-like flaps are positioned on bendable membranes out of the center of regularly separated cavities. Each flap row can be deflected separately by an induced pneumatic pressure from the back side into the cavity, which then leads to a bending of the supporting membrane. Due to the high aspect ratio of the flaps the membrane deflection results in a considerable large transversal motion of the flap tips. In addition, as a reaction to the pressure profile the outward motion may differ in the temporal profile from the inward motion of the membrane, an asymmetry in the temporal motion cycle between the forward and backwards directed stroke is imposed. The bending cycle between adjacent channels can be adjusted so that it replicates a metachronal wave.

Design Concept

Layout details of the membrane and the flap structure are shown in Figure 1. The width of the supporting structure was set to 400 µm and the width of the cavity to 600 µm, respectively. Typically, 500 µm high (h) and 50 µm wide flaps were used as artificial comb rows that extend into the fluid from the flexible surface wall. To move the fluid effectively and not only in the vicinity of the membrane the design must show some kind of asymmetry, therefore the flap is positioned out of center. This first asymmetry enables the flap to transport the fluid in the desired direction.

A second asymmetry is induced by the asymmetric movement of the flap. By applying the metachronal wave, the movement of each flap is strongly coupled to its neighbors. This coupling enforces the stroke of the flap by changing the beat path. The membranes are stretched and compressed by the neighboring membranes. When a membrane is moving upwards the pressure in the cavities in front are pushing the flap forward. This increases the absolute height of the flap in the fluid during the effective stroke.

On the other hand, the contracted membranes behind the recovering membrane are pulling the flap back and down. The behavior is shown in Figure 3 which is composed of superposed pictures of a beating single flap. The colored track represents the path of one edge of the flap tip. The total height difference on the 342 µm long path is 85 µm which corresponds to 25% relative to the flap height.

Furthermore, a third type of asymmetry corresponds to the applied pressure profile under the membranes, which is also advantageous for the fluid transport: The pressure behavior depends on several parameters: magnitude and direction of the pressure change, fluid impedance and resistance of the connecting tubes between valve and connector. For illustration, Figure 4 shows the pressure profiles for four applied pressures measured in one cavity for a 500 mm long tube. The black curve refers to the valve position and reflects the switching behavior of the valve.

After switching the valve from high to low pressure level the pressure in the cavity decreases until it reaches the lower level. A higher pressure difference corresponds to a longer decay time for reaching the lower level, e.g. a pressure step from +0.75 bar to −0.75 bar takes 90 ms, whereas a step from +0.2 bar to –0.2 bar takes only 30 ms. The opposite pressure change from low to high level is significantly faster. Here, a pressure change from –0.75 bar to +0.75 bar takes 45 ms, whereas a step from +0.2 bar to –0.2 bar takes only 20 ms.

The characteristic pressure delay behavior between the applied pressure and the pressure in the cavity under the membrane can be used for generating an asymmetric behavior of the flap movement. A fast positive pressure difference induces the energy faster into the system then the slower negative pressure difference does. Following the notation of biology, a positive switch (from low to high) defines the effective stroke of the flaps, whereas a negative switch (from high to low) represents the slow recovery stroke. This anisotropy is assumed to support the net transport in direction of the effective stroke. To create a travelling wave mechanism in the flap array pressures were switched from +0.6 bar to –0.6 bar to all of the twenty membranes consecutively, as demonstrated in a simulation model (Figure 5).

Simulations

Simulations were carried out using ANSYS 12 Structural program based on a 2D model with high deformation. The material data for PDMS (poly(dimethylsiloxane), Sylgard 184, Dow Corning, Midland, USA), which serves as suitable material for the membrane and flaps, was defined as an isotropic material model. The material parameters were generated by tensile testing [19]. Optimum membrane bending conditions were extracted for a 600 µm wide and 100 µm thick PDMS membrane. The supporting structure between two membranes was optimized to 400 µm. The van Mises stress distribution of the PDMS membrane is shown in Figure 6, when the membrane is deflected by a pressure of 40 kPa. With respect to the soft membrane material, the bending deforms the membrane shape and increases the membrane area. These deformations appear in both bending directions of the membrane and are responsible for a shift of the tip deflection.

The asymmetric flap position on the membrane results in an asymmetrical movement of the flap. Figure 6 shows the simulated movement of the flap by setting a minimum and maximum pressure difference of 40 kPa and -30 kPa compared to ambient pressure. The flap moves 480 µm in x-direction and 75 µm in z-direction. This behavior generates a lateral movement in the fluid. A similar characteristic movement was published by Khaderi et al. [20].

 

Results

At first, the real flap movement were analyzed and compared to the simulation results. As mentioned before, the coupling of neighboring membranes enforces the stroke of the flap by increasing the absolute height of the flap in the fluid during the effective stroke. This behavior was simulated in Figure 6 and is verified in Figure 3 experimentally. Both, simulation and experiment showed the same difference in forward and backward stroke. In the simulation the tip movement of 480 µm in x-direction and 75 µm in z-direction was calculated. This is a relative difference of 16%.

This behavior is crucial to transport fluids in the regime of laminar flow, because in general a reversible movement in a laminar flow regime cannot produce an effective transport of the fluid. In a rough estimation the Reynolds number Re can be calculated to:

where r is the density (1060 kg/m³), L the height of the fluid (6 mm) and m the viscosity (0.001 Pa*s), respectively. The characteristic velocity v is hard to determine, but can be estimated to be between the maximum velocity of flap tip in the power stroke of 7600 µm/s and the total flow velocity of the fluid of 40 µm/s. Using this velocity range, Re varies from 0.25 to 48. Both values are in the range of laminar flow. This indicates the importance that only asymmetries in the stroke result in fluid transport, which is also known from natural cilia and comb-rows [21, 22].

To perform fluid transport in the system the transporter device was connected to a pneumatic triggering module and integrated into a microfluidic chamber. A glycol-water mixture of 3:10 was used as fluid to ensure levitating particles. Fluorescence-marked 60 µm polystyrene particles were used for experiments. Figure 7 shows a snap-shot of the proof-of-concept. Several particles are shown in the flap array. The different diameters correspond to the different vertical positions the particles have cor-responding to the focus depth of the objective lens (focus is adjusted on top of the flap). The trajectories of some particles are shown. The particle transportation direction was from the left to the right. Depending on the relative position between particle and flap, a net transport velocity was determined to be 40 µm/s.

Depending on the vertical distance of the particle from the membrane surface the particle velocity varies considerably. Particle tracks in vicinity to the flap tips at y ≈ h showed often straight lines (lilac particle, 1) with velocities up to 64 µm/s, whereas particles farther above the tip at y ≈ 2h showed less flow speed (light blue particle, 2) of 22 to 31 µm/s. This behavior indicates a rather homogenous fluid stream with high velocity next to the flap tips. The net transport velocity was determined from all tracks to approximately 40 µm/s. Similar fluid speeds were reported using discontinued flap arrays with a single flap length of 0.9 mm and a gap of 0.1 µm between two flaps [17]. In addition, the same order of flow speeds was calculated in numerical simulations of flow generated by the cilia on Ctenophores [23] using flaps of the same scale as in our experiments.

Conclusion

The active movement of comb row like flaps induces a fluid flow. The flaps are actuated through pneumatic pressure applied to the backside of a membrane. Simulations are carried out to optimize the design. The asymmetric position of the flaps on the membrane results in an asymmetric movement of the flaps. Mean velocities of 40 µm/s were obtained by applying metachronal waves to the flap array.

In future devices the transport efficiency must be optimized. The device dimensions must be significantly reduced to be able to integrate it in microsystems, e.g. for microfluidic applications.

Acknowledgement

Funding was given by the German Federal Ministry of Education and Research under grant #16SV5341 (PaTra).

Uwe Schnakenberg, Dr.-Ing.

Institute of Materials in Electrical Engineering 1

RWTH Aachen University

Sommerfeldstraße 24, 52074 Aachen, Germany

E-mail: schnakenberg@iwe1.rwth-aachen.de

Phone: +49 241 80 2784