Research ArticleAPPLIED SCIENCES AND ENGINEERING

Electrostatically driven fog collection using space charge injection

See allHide authors and affiliations

Science Advances  08 Jun 2018:
Vol. 4, no. 6, eaao5323
DOI: 10.1126/sciadv.aao5323

Abstract

Fog collection can be a sustainable solution to water scarcity in many regions around the world. Most proposed collectors are meshes that rely on inertial collision for droplet capture and are inherently limited by aerodynamics. We propose a new approach in which we introduce electrical forces that can overcome aerodynamic drag forces. Using an ion emitter, we introduce a space charge into the fog to impart a net charge to the incoming fog droplets and direct them toward a collector using an imposed electric field. We experimentally measure the collection efficiency on single wires, two-wire systems, and meshes and propose a physical model to quantify it. We identify the regimes of optimal collection and provide insights into designing effective fog harvesting systems.

INTRODUCTION

There are more than 1.1 billion people who lack access to safe drinking water worldwide, according to the World Water Council (1). Many remote drought-prone coastal areas have little or no rain and prohibitively expensive water transportation costs but have dense fog that occurs on a regular basis (2). Fog harvesting is a promising solution to provide clean water in these regions. Areas prone to dense fog formation are usually close to oceans, where fog clouds form over the water and are then transported by the wind (3). Fog is composed of tiny droplets with a typical diameter of 10 μm. Researchers have designed various artificial fog harvesting systems, some of which mimic natural fog collection mechanisms in animals and plants (411), and have successfully implemented small-scale fog collectors (1214).

The most common design for fog collectors is a mesh that stands perpendicular to the fog-laden wind (13, 1517). Upon impact, droplets stick to the mesh and grow as they coalesce with other incoming droplets. When they reach a critical size, they shed by gravity into a container. Meshes are used because they cause a smaller deviation of the incoming flow streamlines by letting air pass through their openings. Nevertheless, field studies have shown that these meshes typically have very low efficiencies of around 1 to 2% (17).

The main mechanisms that limit the efficiency of mesh-based fog collectors are the shedding rate and the aerodynamic deviation of fog droplets. The shedding rate can decrease the efficiency by two mechanisms. Water can clog the mesh openings, making the mesh locally act as a plate. Re-entrainment of the droplets due to wind drag may also occur before these droplets are collected. Researchers have developed coatings to improve the shedding rate, but the overall efficiency of these meshes remained low, around 10% in laboratory setups, suggesting that the main limitation is not the shedding rate (7, 9, 17).

The main limitation is the aerodynamic deviation of fog droplets, which can occur on two scales: the mesh (aerodynamic efficiency ηa) and individual mesh wires (deposition efficiency ηd). The overall collection efficiency is η = ηaηd. ηa is the number of droplets directed toward the mesh wires divided by the number of droplets that were directed toward the collector far upstream (15). It depends on the shading coefficient SC (surface fraction of the wires), and it has been shown that an SC around 55% gives a maximum ηa (16, 17).

A more significant bottleneck in the fog collection process is the deviation of the droplets around the individual wires. ηd is defined as the ratio of captured droplets to those initially directed toward the wire. The flow through the mesh wires can be modeled as a flow past a cylinder (18). Close to the cylinder of radius Rc, in a region of characteristic size Rc, air streamlines start deviating and the flow separates, as schematically shown in Fig. 1A. The ability of droplets to follow streamlines is characterized by the Stokes number, which is the ratio of the droplet inertia to the drag force (18). Embedded Image, where Rd is the droplet radius, ρw is the water density, U is the air speed, and ηg is the air viscosity. At low Stokes number, the droplet trajectories follow the streamlines closely, and few droplets are collected. Figure 1B and movie S1 show the air streamlines and droplet trajectories around a cylindrical wire for St = 0.05. At high Stokes number, drag forces do not affect the trajectories, and droplets directed toward the cylinder continue along their trajectories and collide with the cylinder. An empirical formula has been established: Embedded Image (17). However, large Stokes numbers require very fine meshes, which are difficult to fabricate and lack structural integrity. Hence, low deposition efficiencies remain a significant challenge in fog collection.

Fig. 1 Trajectories of fog droplets around a cylinder with and without the application of corona discharge.

(A and B) Schematic of air streamlines and droplet trajectories and photograph of droplet trajectories in the absence of an electric field. The inset shows the velocities of the wind Embedded Image and the particle Embedded Image and the drag force acting on the droplet. The bright ring is the edge of the cylinder. (C and D) Schematic of air streamlines, electric field lines, and droplet trajectories and photograph under corona discharge. Droplets closely follow the electric field lines in this case. The inset shows the additional electric force acting on a droplet. The cylinders in (B) and (D) have a diameter of 1.88 mm.

Here, we propose to overcome the fundamental aerodynamic limitation of streamline deviation around the collector wires by introducing an additional electric force that will overcome the aerodynamic drag force and propel the fog droplets toward the collector. Active control of droplets with electromagnetic and other fields has recently received significant attention (1921). Here, we develop a new approach to enhancing fog collection, inspired by the principles of electrostatic precipitators (2224). The mechanism is schematically shown in Fig. 1C. We use an ion emitter to inject a net charge into the droplets (25) and use an electric field to direct them toward a collector. The electric field lines go from the emitter to the grounded collector, as shown in Fig. 1C. When the electrical forces are much larger than the air drag forces, the fog droplets follow the field lines and deposit on both sides of the collector wire. We demonstrate this concept experimentally in Fig. 1D and movie S2, which shows that droplets are collected all over the wire and that some droplets that were not initially directed toward the wire are also captured. Intuitively, this suggests that the deposition efficiency, which was previously defined as the fraction of droplets directed toward the wire that are collected, can become greater than one as we will quantitatively show in this paper. The physics of the collection on a single wire will be generalized to a system of two wires and to the more general case of a mesh.

RESULTS AND DISCUSSION

Collection on a single cylindrical wire

The experimental setup is shown schematically in Fig. 2A. A sharp electrode that serves as an ion emitter is placed at a distance L from the horizontal cylindrical wire of radius Rc, which serves as a collector. LRcRd (droplet radius). The collector is electrically grounded, while a high voltage V is applied to the emitter electrode, thus establishing an electric field between them. When V is above a critical value Vc (26), corona discharge occurs, the air surrounding the emitter is ionized, and a plasma region forms. Electrons are accelerated and have enough energy to ionize the air atoms when they collide with them. A chain reaction occurs, with every collision creating additional electrons and ions. After a collision, the electric field pulls electrons and ions in opposite directions, preventing recombination. At a certain distance from the emitter, the electric field can no longer provide enough energy to the electrons to sustain the reaction. In this region, ions travel freely in the air toward the opposite electrode, collide with fog droplets, and are captured by them. Thus, the droplets acquire a net charge q of the same polarity as V (22, 27, 28).

Fig. 2 Mechanism of droplet collection on a cylindrical wire.

(A) Schematic of simplified experimental setup and droplet trajectories. (B) Schematic of the acceleration phase undergone by droplets. The electric field, the initial and terminal velocities, as well as the forces acting on a droplet are shown. (C) Added velocity as a function of V2. A linear fit of the data (R2 = 0.94) gives a slope of 0.006 m/s per kV2. The gray area is where the voltage is not high enough to induce corona discharge. The error bars reflect the SD over four measurements. (D) Schematic of the cross section of the collection phase near the cylinder. Streamlines, field lines, and trajectories of the droplets are shown. (E) Nondimensional collection area as a function of V2 for four different wind speeds. The gray area is where there is no corona discharge.

As L >> Rc, apart from small regions around the emitter and the collector, the electric field lines are essentially parallel in the central region and the electric field Embedded Image. We assume in the following analysis that the electric field is not disrupted by the presence of fog, and we neglect interdroplet interactions.

In the central region, fog droplets accelerate (Fig. 2B). The droplets enter the region with the wind velocity U0, and the electric force Fe = qE accelerates them. As their velocity increases, they are subject to a drag force Fd from the air. When Fe = Fd, the droplet reaches a terminal velocity Uf.

The equation of motion for the droplet isEmbedded Image(1)where Embedded Image is the droplet’s velocity, m is its mass, and Embedded Image is the air velocity.

The acceleration time scale Embedded Image is one to two orders of magnitude smaller than the time of flight Embedded Image. Thus, the droplets reach their terminal velocityEmbedded Image

Figure 2C shows experimental measurements of the added velocity UfU0 as a function of V2. The gray region of the plot corresponds to the case V < Vc, where no charge injection occurs. In this regime, the electric force acting on the drops arises because of dielectrophoresis, whose magnitude ~ Embedded Image (29, 30) is six orders smaller than that of drag force and, hence, no added velocity is observed. Upon the onset of corona discharge, droplets are charged and the added velocity is proportional to V2 since both q and E are proportional to V.

As droplets are conductive, the charge is localized at the surface, and the surface charge per unit area can be estimated as σ ~ ϵ0E (31).

Droplets gain charge when ions attach to them but cannot lose charge to the air (insulator), so the final charge of a droplet is determined by the electric field that it encounters. We can then estimate the surface charge as Embedded Image, and more precise calculations give Embedded Image (see the Supplementary Materials) (22).

This calculation is valid here because the time of flight τf ~ 10− 2s is much larger than the charge relaxation time in the droplet Embedded Image, where σw is the conductivity of deionized water (2.10− 4S/m). Thereby, fog droplets have enough time to acquire their maximum charge before they reach the collector (31).

Thus, the acquired charge is given byEmbedded Image(2)and the added velocity isEmbedded Image(3)

The linear fit to the data in Fig. 2C gives a slope of 0.006. Using the above formula in our case, we find that UfU0 = 0.008 V2, which matches well with the experimental values, considering the uncertainties in measuring the parameters.

After the acceleration phase, the fog droplets approach the collector wire. In this region (Fig. 2D), both the streamlines and the electric field lines deviate from the parallel configuration. The electric field near the collector scales as Embedded Image (32). Near the cylinder, electric forces dominate, but further from the cylinder, the electrical forces are weaker and unable to overcome inertia of the droplets. Therefore, there will be a distance beyond which droplets will not be collected for a given intensity of the applied voltage. Assuming an infinitely long cylinder, the problem reduces to a two-dimensional flow. We then define the area of collection A as the cross-sectional area of the flow of incoming droplets that are collected per unit depth of the cylinder. The deposition efficiency is then Embedded Image, where A0 is the area of the cylinder projected in the direction of flow, equal to 2Rc. We experimentally measured the area of collection for different voltages and wind speeds by imaging the region around the cylinder. Figure 2E presents the nondimensional area A/A0 as a function of V2 for different values of U0. We find that A/A0 increases with V2 and decreases with wind speed (fig. S1).

To rationalize these observations, we nondimensionalize Eq. 1 by dividing all distances by Rc, velocities by Uf, and time by Embedded Image, the characteristic time in the collector region (23, 3337) Embedded Image(4)Ke is the ratio of electric and viscous forces, and we will call it the electrical number. Apart from these nondimensional numbers, all the terms of the equation are of order 1. Around the cylinder, E scales as V/Rc and we can use the expressions of q (Eq. 2) and Uf (Eq. 3) to obtainEmbedded Imagewhere we define Embedded Image, a nondimensional number expressing the ratio between the added velocity due to the electric field and the wind velocity.

In the low Stokes number limit (relatively large wires), which is the more practical case, Eq. 4 is further simplified intoEmbedded Image

The electrical number should then govern the physics of the problem, in particular the deposition efficiency. In Fig. 3, we plot A/A0 as a function of Ke and find that the data for various voltages and wind speeds collapse into a linear master curve. The single-wire deposition efficiency is given by

Embedded Image(5)where c can be determined experimentally.

Fig. 3 Dependence of the deposition efficiency on the electrical number.

The data correspond to five values of the wind speed and five values of the voltage. The colors represent different wind speeds (red, 0.55 m/s; blue, 0.6 m/s; green, 1 m/s; yellow, 1.65 m/s; purple, 3.3 m/s). The solid line is a linear fit (R2 = 0.92), with a slope of 0.26.

This linear behavior is expected as Ke represents the relative amplitude of the driving force causing the collection.

We explore the limiting cases of Ke. For low U* (low voltages or high wind speeds), the term containing V2 in the denominator of Ke is negligible, Ke scales linearly as V2, and hence, the deposition efficiency increases linearly with V2. However, for high U*, U0 becomes negligible in the denominator and Ke tends toward a constant (Embedded Image). Any increase in voltage beyond this point does not lead to more collection of droplets. We will call this the voltage saturation. Voltage saturation does not change the linear behavior of the collection rate but limits the maximum reachable Ke. Voltage saturation sets in as Ke approaches this maximum value. Ke can be written as Embedded Image. Complete voltage saturation occurs when U* >> 1, but its effects start to be visible when U* becomes of order 1. U* ranged here from ~0.05 to ~1.2, and voltage saturation is observed experimentally in the bending of the experimental curves in Fig. 2E for high voltages and low wind speeds, which correspond to the highest U*. Qualitatively, this limitation in the reachable collection rates arises because, as we increase the voltage, we have a higher electric force that both attracts the droplets and increases their speed, thereby increasing the drag force. Eventually, at high voltages, these two effects balance each other, and the collection cannot be enhanced any further.

Collection on two parallel cylindrical wires

To extend the single-wire model to an array of wires, we now investigate fog collection on a system of two parallel cylindrical wires. We assume that, when the wires are far from each other, air streamlines, electric field lines, and droplet trajectories around each cylinder are not affected by the presence of the other cylinder (Fig. 4, A and B). Thus, the cylinders behave as two single wires. However, when they are close enough, electric field lines and droplet trajectories are influenced by the neighboring wire. When the wires are too close, because of geometric confinement, the droplets above the plane of symmetry are propelled toward the upper wire, while the droplets below go to the lower wire. The deposition efficiency is limited by the system’s geometry because of the competition between the two wires over the same droplets, as shown schematically in Fig. 4C. We define two areas of collection, Ain and Aout, accounting for the projected area of the flow of incoming droplets that are collected in the regions between the wires and outside them, respectively. The corresponding deposition efficiencies ηd, in and ηd, out are defined as the ratios of these areas to their respective collection area A0 = Rc.

Fig. 4 Mechanism of droplet collection on two parallel cylindrical wires.

(A) Schematic of droplet trajectories with two distant wires. The collection area of each single wire Asw and the distance D are shown. (B) Schematic of droplet trajectories in a two-wire system with spacing saturation. The parameters Ain, Aout, and D are shown. The white arrows show the single-wire areas of collection Asw. (C) Photograph of the droplet trajectories for two distant wires. The wire diameters are 1.88 mm. The distance between them is 10 mm. The applied voltage is 10 kV. (D) Photograph of the droplet trajectories in a spacing saturation case. The wire diameters are 1.88 mm. The distance between them is 6 mm. The applied voltage is 14 kV. (E) Embedded Image as a function of Ke for three different wire distances. The gray region covers theoretically inaccessible values. The vertical dashed lines represent the predicted spacing saturation values for D* = 1.7 and 4.2.

Figure 4D is a photograph showing the two cylinders and the droplet trajectories, in a case where D is small enough that most of the droplets between the wires are collected. When this happens (small D), a simple geometrical analysis shows that Ain reaches a saturation value that is equal to Embedded Image, and Embedded Image, where Embedded Image. We will call this limitation spacing saturation.

We hypothesize that Embedded Image and Embedded Image grow with the single wire law until spacing saturation is reached, at which point Ain will plateau, while Aout will still follow the single wire law. Knowing that saturation happens when Embedded Image, we can use the single wire law (Eq. 5) to predict the electrical number Kesat (or equivalently the voltage) or the critical distance Embedded Image at which spacing saturation happens.

We experimentally measured Ain and Aout for different Ke and D*. The results for Ain are reported in Fig. 4E (Aout data are reported in fig. S2). We observe that Embedded Image increases with Ke and plateaus at 1 for low values of D*. Before the plateau, Embedded Image increases linearly with Ke, with a slope equal to what was measured in the single-wire case (fig. S2). The zone of the graph representing Embedded Image is inaccessible because of the competition over a finite number of droplets that leads to saturation. The vertical dashed lines represent the expected saturation Kesat for the cases of D* = 1.7 and 4.2 from our model, and Embedded Image plateaus beyond these lines for the corresponding curves.

For D* = 11, spacing saturation is not observed here and is expected to occur at a much higher Ke. However, in our case, such a high Ke cannot be reached because of the voltage saturation, which will then become the limiting factor here.

Collection on meshes

Having rationalized the behavior of two parallel wires, we studied the macroscopic collection on the ultimate system of interest: meshes. We conducted similar experiments on 5 cm × 5 cm meshes, and we measured the mass of the collected water after 5 min of exposure. We used five meshes that had the same wire radius, but the spacing D between the wires increased from mesh 1 to mesh 5 (details in Materials and Methods). We quantified the efficiency of the collection process as the ratio of the collected mass to the total mass of water directed toward the mesh Embedded Image. As stated before, η = ηaηdηs, where Embedded Image is the deposition efficiency. de Dios Rivera have studied the aerodynamic efficiency previously (16), and it is given by Embedded Image, where SC is the shading coefficient. While wettability can improve the shedding efficiency ηs in passive collectors by facilitating shedding and preventing clogging and re-entrainment (17), its effect here is minimal. Re-entrainment was not observed in any of our experiments. Some local clogging was observed. In passive inertial collectors, clogging of pores leads to a local deviation of the flow and loss of fog droplets that pass through the neighboring openings. However, droplets here are charged and attracted to the collector, so they can still be collected even if they are deviated and pass through another opening. Therefore, the influence of local clogging is minimal and we took ηs ~ 1 for all cases.

It should be noted that mtotal was taken as the quantity of fog generated by our fog machine and does not take into account evaporation and losses of fog before reaching the mesh. Thereby, the reported efficiency is a lower bound for the actual efficiency.

The collection on meshes with corona discharge is drastically increased. As shown in Fig. 5 (A and B) and movies S3 and S4, when there is no corona discharge, there is hardly any collection on the mesh, and the collection beaker remains empty. In contrast, when a corona discharge is established, the mesh is covered with water within a few seconds, and the collection beaker (30 ml) is filled within 30 min. Longer experiments of up to 10 hours were also performed and showed a steady performance (fig. S3)

Fig. 5 Fog collection on meshes.

(A) Snapshots of meshes at different time intervals of fog exposure. In the first row, a 15-kV voltage was applied, while there was no electric field in the second row. Red dye was added to the dispersed fog for visualization purposes. (B) Photographs showing the collection mesh and the storage beaker for collected water after 30 min of exposure. The case with high voltage resulted in the collection of 30 ml of water, while only three droplets were collected without electric field. Complete video of collection is shown in movies S3 to S4. (C) Mass of the collected water as a function of Ke for different meshes. The vertical dashed lines represent the predicted onset of spacing saturation from the two-wire model for meshes 1 to 3. (D) Deposition efficiency of the five meshes as a function of Ke. The colors represent different meshes according to the color code of (C).

Figure 5C shows the mass of collected water after 5 min of exposure for different Ke values. The mass m increases with Ke for all meshes. We notice that the mesh with the highest collection is not the one with the highest wire density; mesh 2 has fewer wires but still collects more water because of a higher aerodynamic efficiency. Finally, the collection of meshes 1 and 2 starts to slow down at high Ke, which is due to spacing saturation.

Finally, we plotted the deposition efficiency Embedded Image for the cases below voltage saturation as a function of Ke in Fig. 5D. We see that the data collapse around a linear curve, which shows that our models for single and double wires can be extended to meshes. The slope is lower than that for the case of a single wire, which may be due to evaporation and leaks as well as the nonuniform distribution of fog across the mesh. Still, ηd reaches values over 2 and can therefore compensate for other inefficiencies and lead to a significantly high overall efficiency for the collection system.

While the operating voltages that we are using are on the order of 10 kV, the current is low and on the order of 0.01 mA. The power consumption in our experimental setup is thus ~0.1 W for a 5 cm × 5 cm mesh, which translates to about 40 W/m2. Given the output of our experimental system (>60 g/hour), the energy consumption is on the order of 2 kWh/m3, which is less than current reverse-osmosis desalination consumption (~3 to 5 kWh/m3). In addition, the power consumption is due to the formation and transport of ions. It does not depend on the presence of fog. In particular, one of the major applications of this study is the recovery of water from cooling tower plumes in power plants, where fog is much more concentrated and is on the order of 5000 liters/m2 per day. In that case, the energy consumption is estimated to be ~0.2 kWh/m3, which is one to two orders of magnitude lower than the energy consumption to produce clean water from conventional sources.

In practice, for a given mesh, we can use our scaling law to predict the deposition efficiency as a function of ambient conditions. We have shown that four nondimensional numbers govern the problem: St, Ke, U*, and D*. ηd depends on the first two numbers and, in the case of low St, ηd is proportional to Ke. U* and D* predict two important limitations, which are the voltage and spacing saturations. We can adjust the voltage according to wind speed and droplet size to maintain a constant throughput of harvested fog water. When designing an active collection mesh, our model also offers design guidelines for the wire size, wire spacing, and operating voltage that should be chosen so as to remain below the voltage saturation and the spacing saturation in typical operating conditions to maximize collection while conserving energy. This method may be combined with other techniques, such as surface treatment of the collector, to slightly enhance the efficiency. These results can be used to design efficient fog harvesters in drought-prone areas to collect water for drinking, irrigation, and afforestation. It is also possible to apply this method to power plant cooling towers to recover water that is lost in the plumes. Fog removal systems to increase visibility on roads and airports may also be designed.

MATERIALS AND METHODS

Experimental setup and procedure

Samples were placed 4 cm away from the outlet, perpendicular to the axis, of two concentric cylinders (inner diameters, 6.3 and 5 cm) from which came a uniform stream of fog. The fog consisted of a cloud of air-suspended water droplets of an average radius of 3.5 μm, generated using an ultrasonic humidifier (Air-O-Swiss AOS 7146) delivering a volume rate of up to 0.1 liters/hour. The average droplet size was measured with high-magnification imaging of fog droplets. Fog was generated directly into the smaller cylinder through an orifice. Fog was released upstream of the emitter. The fog droplets, the emitter, and the collector were all at the same temperature (room temperature), and the relative humidity in fog conditions was close to 100%. These are also the typical conditions in practical applications. A speed-tunable fan (Thermaltake Mobile Fan II External USB Cooling Fan) was placed at the inlet of the larger cylinder to create the airflow that would convect the fog toward the collection area. A honeycomb flow straightener (Saxon Computers 120mm Honeycomb Airflow Straightener) was placed after the fan to ensure that the wind velocity was uniform through the area of the cylinder, thus reproducing real-fog conditions. The outlet velocity was measured with an anemometer (Testo 405 Hot Wire Thermo-Anemometer) and was spatially uniform within a 15% interval. Corona discharge was produced by placing a sharp metallic needle inside the cylinders, its tip coinciding with the outlet of the smaller cylinder. The needle was connected to a high-voltage generator (Spellman SL600) delivering voltages from 0 to −25 kV. Corona discharge was observed to start at a voltage around −7 kV. The onset of corona was detected by the sudden increase in current as ions were generated and traveled between the two electrodes and by the sudden change in the flow pattern of fog. In all experiments, the collector was connected to the ground, its voltage at 0 V. All experiments were performed under ambient temperature and humidity conditions. The collection area was determined from measuring the radius at which the limit between trajectories of collected drops and noncollected drops occurs, in photographs like Fig. 1D. We used five different photographs for each case to do the measurements, and we reported the average area.

Wires and meshes

In single- and two-wire experiments, cylindrical needles, made of stainless steel, with a length of 4 cm and a diameter of 1.88 mm, were used as collectors.

Five-centimeter square meshes were used for collection tests. They were purchased from McMaster-Carr (Corrosion-Resistant Type 304 Stainless Steel Wire Cloth), and their individual characteristics are given in Table 1.

Table 1 Characteristics of meshes.
View this table:

SUPPLEMENTARY MATERIALS

Supplementary material for this article is available at http://advances.sciencemag.org/cgi/content/full/4/6/eaao5323/DC1

Supplementary Text

fig. S1. Nondimensional collection area as a function of the inverse of wind speed for five different voltages.

fig. S2. Nondimensional collection area for two wires as a function of Ke for three different wire distances.

fig. S3. Continuous 10-hour fog collection experiment.

movie S1. Droplet trajectories in the absence of an electric field.

movie S2. Droplet trajectories with corona discharge.

movie S3. Ninety-second video of fog collection on meshes with and without corona discharge.

movie S4. Thirty-minute video of fog collection on meshes with and without corona discharge.

This is an open-access article distributed under the terms of the Creative Commons Attribution-NonCommercial license, which permits use, distribution, and reproduction in any medium, so long as the resultant use is not for commercial advantage and provided the original work is properly cited.

REFERENCES AND NOTES

Acknowledgments: We thank S. R. Mahmoudi for his contributions to the initial brainstorming for this work. Funding: We are grateful for support from the Massachusetts Institute of Technology (MIT) Tata Center for Technology and Design. Author contributions: M.D. and K.K.V. conceived the project. M.D. and K.K.V. analyzed the data and wrote the paper. M.D. carried out the experiments. K.K.V. supervised the work. Competing interests: M.D. and K.K.V. are inventors on a patent application related to this work filed by the MIT (no. PCT/US2016/054230, filed 28 September 2016). The other authors declare no competing interests. Data and materials availability: All data needed to evaluate the conclusions in the paper are present in the paper and/or the Supplementary Materials. Additional data related to this paper may be requested from the authors.
View Abstract

Navigate This Article