Skip to main content

Thank you for visiting You are using a browser version with limited support for CSS. To obtain the best experience, we recommend you use a more up to date browser (or turn off compatibility mode in Internet Explorer). In the meantime, to ensure continued support, we are displaying the site without styles and JavaScript.

Cellular fluidics


The natural world provides many examples of multiphase transport and reaction processes that have been optimized by evolution. These phenomena take place at multiple length and time scales and typically include gas–liquid–solid interfaces and capillary phenomena in porous media1,2. Many biological and living systems have evolved to optimize fluidic transport. However, living things are exceptionally complex and very difficult to replicate3,4,5, and human-made microfluidic devices (which are typically planar and enclosed) are highly limited for multiphase process engineering6,7,8. Here we introduce the concept of cellular fluidics: a platform of unit-cell-based, three-dimensional structures—enabled by emerging 3D printing methods9,10—for the deterministic control of multiphase flow, transport and reaction processes. We show that flow in these structures can be ‘programmed’ through architected design of cell type, size and relative density. We demonstrate gas–liquid transport processes such as transpiration and absorption, using evaporative cooling and CO2 capture as examples. We design and demonstrate preferential liquid and gas transport pathways in three-dimensional cellular fluidic devices with capillary-driven and actively pumped liquid flow, and present examples of selective metallization of pre-programmed patterns. Our results show that the design and fabrication of architected cellular materials, coupled with analytical and numerical predictions of steady-state and dynamic behaviour of multiphase interfaces, provide deterministic control of fluidic transport in three dimensions. Cellular fluidics may transform the design space for spatial and temporal control of multiphase transport and reaction processes.

Access options

Rent or Buy article

Get time limited or full article access on ReadCube.


All prices are NET prices.

Fig. 1: Cellular fluidics: a nature-inspired, unit-cell-based platform.
Fig. 2: Architected design for predictable capillary rise.
Fig. 3: Numerical simulations and high-speed videography reveal dynamics.
Fig. 4: Transpiration in cellular fluidics.
Fig. 5: Gas absorption in cellular fluidics.
Fig. 6: Preferential fluid pathways for continuous flow and selective patterning.

Data availability

Source data for calculated porosities and interfacial areas of different unit cells are included with the paper. Reasonable requests for additional data that support the findings of this study will be honoured upon securing release approval from the Lawrence Livermore National Laboratory and the US Department of Energy. Source data are provided with this paper.

Code availability

Reasonable requests for the code that support the findings of this study will be honoured upon securing release approval from the Lawrence Livermore National Laboratory and the US Department of Energy.


  1. 1.

    Kirkham, M. B. in Principles of Soil and Plant Water Relations 2nd edn (ed. Kirkham, M. B.) Ch. 15, 243–266 (Academic, 2014).

  2. 2.

    Grigoryan, B. et al. Multivascular networks and functional intravascular topologies within biocompatible hydrogels. Science 364, 458–464 (2019).

    ADS  CAS  PubMed  PubMed Central  Google Scholar 

  3. 3.

    Yaraghi, N. A. & Kisailus, D. Biomimetic structural materials: inspiration from design and assembly. Annu. Rev. Phys. Chem. 69, 23–57 (2018).

    ADS  CAS  PubMed  Google Scholar 

  4. 4.

    Huang, W. et al. Multiscale toughening mechanisms in biological materials and bioinspired designs. Adv. Mater. 31, 1901561 (2019).

    CAS  Google Scholar 

  5. 5.

    Torres, A. M. et al. Bone-inspired microarchitectures achieve enhanced fatigue life. Proc. Natl Acad. Sci. USA 116, 24457–24462 (2019).

    ADS  CAS  PubMed  PubMed Central  Google Scholar 

  6. 6.

    Whitesides, G. M. The origins and the future of microfluidics. Nature 442, 368–373 (2006).

    ADS  CAS  Google Scholar 

  7. 7.

    Shui, L., Pennathur, S., Eijkel, J. C. & van den Berg, A. Multiphase flow in lab on chip devices: a real tool for the future? Lab Chip 8, 1010–1014 (2008).

    CAS  PubMed  Google Scholar 

  8. 8.

    Pennathur, S. Flow control in microfluidics: are the workhorse flows adequate? Lab Chip 8, 383–387 (2008).

    CAS  PubMed  Google Scholar 

  9. 9.

    Zheng, X. et al. Ultralight, ultrastiff mechanical metamaterials. Science 344, 1373–1377 (2014).

    ADS  CAS  PubMed  Google Scholar 

  10. 10.

    Zheng, X. et al. Multiscale metallic metamaterials. Nat. Mater. 15, 1100–1106 (2016).

    ADS  CAS  PubMed  Google Scholar 

  11. 11.

    Piorek, B. D. et al. Free-surface microfluidic control of surface-enhanced Raman spectroscopy for the optimized detection of airborne molecules. Proc. Natl Acad. Sci. USA 104, 18898–18901 (2007).

    ADS  CAS  PubMed  PubMed Central  Google Scholar 

  12. 12.

    Delamarche, E. & Kaigala, G. V. (eds) Open-Space Microfluidics: Concepts, Implementations, Applications (John Wiley & Sons, 2018).

  13. 13.

    Olanrewaju, A., Beaugrand, M., Yafia, M. & Juncker, D. Capillary microfluidics in microchannels: from microfluidic networks to capillaric circuits. Lab Chip 18, 2323–2347 (2018).

    CAS  PubMed  Google Scholar 

  14. 14.

    Martinez, A. W., Phillips, S. T. & Whitesides, G. M. Three-dimensional microfluidic devices fabricated in layered paper and tape. Proc. Natl Acad. Sci. USA 105, 19606–19611 (2008).

    ADS  CAS  PubMed  PubMed Central  Google Scholar 

  15. 15.

    Martinez, A. W., Phillips, S. T., Whitesides, G. M. & Carrilho, E. Diagnostics for the developing world: microfluidic paper-based analytical devices. Anal. Chem. 82, 3–10 (2010).

    CAS  PubMed  Google Scholar 

  16. 16.

    Kotz, F. et al. Fabrication of arbitrary three-dimensional suspended hollow microstructures in transparent fused silica glass. Nat. Commun. 10, 1439 (2019).

    ADS  PubMed  PubMed Central  Google Scholar 

  17. 17.

    Schaedler, T. A. et al. Ultralight metallic microlattices. Science 334, 962–965 (2011).

    ADS  CAS  PubMed  Google Scholar 

  18. 18.

    Cui, H., Hensleigh, R., Chen, H. & Zheng, X. Additive manufacturing and size-dependent mechanical properties of three-dimensional microarchitected, high-temperature ceramic metamaterials. J. Mater. Res. 33, 360–371 (2018).

    ADS  CAS  Google Scholar 

  19. 19.

    Jackson, J. A. et al. Field responsive mechanical metamaterials. Sci. Adv. 4, eaau6419 (2018).

    ADS  CAS  PubMed  PubMed Central  Google Scholar 

  20. 20.

    Wang, Q. et al. Lightweight mechanical metamaterials with tunable negative thermal expansion. Phys. Rev. Lett. 117, 175901 (2016).

    ADS  PubMed  Google Scholar 

  21. 21.

    Farzinazar, S., Schaedler, T., Valdevit, L. & Lee, J. Thermal transport in hollow metallic microlattices. APL Mater. 7, 101108 (2019).

    ADS  Google Scholar 

  22. 22.

    Xu, H. & Pasini, D. Structurally efficient three-dimensional metamaterials with controllable thermal expansion. Sci. Rep. 6, 34924 (2016).

    ADS  CAS  PubMed  PubMed Central  Google Scholar 

  23. 23.

    Dudukovic, N. A. et al. Predicting nanoparticle suspension viscoelasticity for multimaterial 3D printing of silica–titania glass. ACS Appl. Nano Mater. 1, 4038–4044 (2018).

    CAS  Google Scholar 

  24. 24.

    Dylla-Spears, R. et al. 3D printed gradient index glass optics. Sci. Adv. 6, eabc7429 (2020).

    ADS  CAS  PubMed  PubMed Central  Google Scholar 

  25. 25.

    Zhu, C. et al. Highly compressible 3D periodic graphene aerogel microlattices. Nat. Commun. 6, 6962 (2015).

    ADS  CAS  PubMed  Google Scholar 

  26. 26.

    Golobic, A. M. et al. Active mixing of reactive materials for 3D printing. Adv. Eng. Mater. 21, 1900147 (2019).

    Google Scholar 

  27. 27.

    Xia, X. et al. Electrochemically reconfigurable architected materials. Nature 573, 205–213 (2019).

    ADS  CAS  PubMed  Google Scholar 

  28. 28.

    Hensleigh, R. et al. Charge-programmed three-dimensional printing for multi-material electronic devices. Nat. Electron. 3, 216–224 (2020).

    CAS  Google Scholar 

  29. 29.

    Zhu, C. et al. Toward digitally controlled catalyst architectures: hierarchical nanoporous gold via 3D printing. Sci. Adv. 4, eaas9459 (2018).

    ADS  CAS  PubMed  PubMed Central  Google Scholar 

  30. 30.

    Fink, K. D., Kolodziejska, J. A., Jacobsen, A. J. & Roper, C. S. Fluid dynamics of flow through microscale lattice structures formed from self-propagating photopolymer waveguides. AIChE J. 57, 2636–2646 (2011).

    CAS  Google Scholar 

  31. 31.

    Li, C. et al. Bioinspired inner microstructured tube controlled capillary rise. Proc. Natl Acad. Sci. USA 116, 12704–12709 (2019).

    ADS  CAS  PubMed  PubMed Central  Google Scholar 

  32. 32.

    Ragelle, H. et al. Surface tension-assisted additive manufacturing. Nat. Commun. 9, 1184 (2018).

    ADS  PubMed  PubMed Central  Google Scholar 

  33. 33.

    Esarte, J., Blanco, J. M., Bernardini, A. & Sancibrián, R. Performance assessment of a three-dimensional printed porous media produced by selective laser melting technology for the optimization of loop heat pipe wicks. Appl. Sci. 9, 2905 (2019).

    CAS  Google Scholar 

  34. 34.

    Rapp, B. E. Microfluidics: Modeling, Mechanics and Mathematics (Elsevier, 2016).

  35. 35.

    Watts, S., Arrighi, W., Kudo, J., Tortorelli, D. A. & White, D. A. Simple, accurate surrogate models of the elastic response of three-dimensional open truss micro-architectures with applications to multiscale topology design. Struct. Multidiscipl. Optim. 60, 1887–1920 (2019).

    MathSciNet  Google Scholar 

  36. 36.

    Nguyen, P. Q., Dorval Courchesne, N.-M., Duraj-Thatte, A., Praveschotinunt, P. & Joshi, N. S. Engineered living materials: prospects and challenges for using biological systems to direct the assembly of smart materials. Adv. Mater. 30, 1704847 (2018).

    Google Scholar 

  37. 37.

    Huang, G., Zhu, Y., Liao, Z., Xu, R. & Jiang, P.-X. Biomimetic self-pumping transpiration cooling for additive manufactured porous module with tree-like micro-channel. Int. J. Heat Mass Transf. 131, 403–410 (2019).

    Google Scholar 

  38. 38.

    Hanks, D. F. et al. High heat flux evaporation of low surface tension liquids from nanoporous membranes. ACS Appl. Mater. Interfaces 12, 7232–7238 (2020).

    CAS  PubMed  Google Scholar 

  39. 39.

    Vericella, J. J. et al. Encapsulated liquid sorbents for carbon dioxide capture. Nat. Commun. 6, 6124 (2015).

    ADS  CAS  PubMed  Google Scholar 

  40. 40.

    Tsunazawa, Y., Yokoyama, T. & Nishiyama, N. An experimental study on the rate and mechanism of capillary rise in sandstone. Prog. Earth Planet. Sci. 3, 8 (2016).

    ADS  Google Scholar 

  41. 41.

    So, J.-H. et al. Reversibly deformable and mechanically tunable fluidic antennas. Adv. Funct. Mater. 19, 3632–3637 (2009).

    CAS  Google Scholar 

  42. 42.

    Vafaei, S., Tuck, C., Wildman, R. & Ashcroft, I. Spreading of the nanofluid triple line in ink jet printed electronics tracks. Addit. Manuf. 11, 77–84 (2016).

    CAS  Google Scholar 

  43. 43.

    Tumbleston, J. R. et al. Continuous liquid interface production of 3D objects. Science 347, 1349–1352 (2015).

    ADS  CAS  PubMed  Google Scholar 

  44. 44.

    Saha, S. K. et al. Scalable submicrometer additive manufacturing. Science 366, 105–109 (2019).

    ADS  CAS  PubMed  Google Scholar 

  45. 45.

    Kelly, B. E. et al. Volumetric additive manufacturing via tomographic reconstruction. Science 363, 1075–1079 (2019).

    ADS  CAS  PubMed  Google Scholar 

  46. 46.

    Walker, D. A., Hedrick, J. L. & Mirkin, C. A. Rapid, large-volume, thermally controlled 3D printing using a mobile liquid interface. Science 366, 360–364 (2019).

    ADS  CAS  PubMed  PubMed Central  Google Scholar 

  47. 47.

    Suss, M. E. et al. Capacitive desalination with flow-through electrodes. Energy Environ. Sci. 5, 9511–9519 (2012).

    CAS  Google Scholar 

  48. 48.

    Wu, L. et al. Highly efficient three-dimensional solar evaporator for high salinity desalination by localized crystallization. Nat. Commun. 11, 521 (2020).

    ADS  PubMed  PubMed Central  Google Scholar 

  49. 49.

    Dinh, C.-T. et al. CO2 electroreduction to ethylene via hydroxide-mediated copper catalysis at an abrupt interface. Science 360, 783–787 (2018).

    CAS  PubMed  PubMed Central  Google Scholar 

  50. 50.

    Weislogel, M. M. et al. The capillary flow experiments aboard the International Space Station: Status. Acta Astronaut. 65, 861–869 (2009).

    ADS  Google Scholar 

  51. 51.

    Stewart, A. Woman tree roots digital design. Thingiverse (2017).

  52. 52.

    Kinstlinger, I. S. & Miller, J. S. 3D-printed fluidic networks as vasculature for engineered tissue. Lab Chip 16, 2025–2043 (2016).

    CAS  PubMed  Google Scholar 

  53. 53.

    Rosàs-Canyelles, E., Modzelewski, A. J., Geldert, A., He, L. & Herr, A. E. Assessing heterogeneity among single embryos and single blastomeres using open microfluidic design. Sci. Adv. 6, eaay1751 (2020).

    ADS  PubMed  PubMed Central  Google Scholar 

  54. 54.

    Nguyen, D. T. et al. 3D-printed transparent glass. Adv. Mater. 29, 1701181 (2017).

    Google Scholar 

  55. 55.

    Kotz, F. et al. Three-dimensional printing of transparent fused silica glass. Nature 544, 337–339 (2017).

    ADS  CAS  PubMed  Google Scholar 

  56. 56.

    Helmer, D. & Rapp, B. E. Divide and print. Nat. Mater. 19, 131–133 (2020).

    ADS  CAS  PubMed  Google Scholar 

  57. 57.

    Gao, L. & McCarthy, T. J. Contact angle hysteresis explained. Langmuir 22, 6234–6237 (2006).

    CAS  PubMed  Google Scholar 

  58. 58.

    Rodríguez-Valverde, M. Á. & Miranda, M. T. Derivation of Jurin’s law revisited. Eur. J. Phys. 32, 49–54 (2011).

    MATH  Google Scholar 

  59. 59.

    Messner, M. C. Optimal lattice-structured materials. J. Mech. Phys. Solids 96, 162–183 (2016).

    ADS  MathSciNet  Google Scholar 

  60. 60.

    NASA Johnson. Capillary flow experiments on Space Station. (2014).

  61. 61.

    de Gennes, P. G. Wetting: statics and dynamics. Rev. Mod. Phys. 57, 827–863 (1985).

    ADS  MathSciNet  Google Scholar 

  62. 62.

    Leger, L. & Joanny, J. F. Liquid spreading. Rep. Prog. Phys. 55, 431–486 (1992).

    ADS  CAS  Google Scholar 

  63. 63.

    Hubbe, M. A., Gardner, D. J. & Shen, W. Contact angles and wettability of cellulosic surfaces: a review of proposed mechanisms and test strategies. BioResources 10, 8657–8749 (2015).

    CAS  Google Scholar 

  64. 64.

    Shusteff, M. et al. One-step volumetric additive manufacturing of complex polymer structures. Sci. Adv. 3, eaao5496 (2017).

    PubMed  PubMed Central  Google Scholar 

Download references


We are grateful to the late J. J. Vericella for his contributions to the CO2 capture technologies that enabled this work. We thank R. Goldsberry and J. Long for support with figure illustrations and artwork; D. Macknelly for providing guidance on lattice generation; C. Ye and J. Mancini for assistance with high-speed video experiments; and J. Knipe, S. Pang, X. Xia, V. Beck and C. Spadaccini for discussions. This work was performed under the auspices of the US Department of Energy by the Lawrence Livermore National Laboratory under contract DE-AC52-07NA27344 and LDRD 19-SI-005 with information management number LLNL-JRNL-811671.

Author information




N.A.D. designed structures, devices and experiments, performed 3D printing, capillary rise tests and analytical calculations, recorded images and videos, and wrote the manuscript. E.J.F. designed and fabricated devices, performed capillary rise, transpiration, and active gas and liquid flow experiments, recorded images and videos, carried out thermal and Schlieren imaging and analysis, and colorimetric analysis of CO2 absorption. H.B.G. fabricated structures, performed capillary rise, active flow, transpiration and CO2 absorption experiments and analysis, performed selective metallization, recorded images and videos. J.R.D. designed structures, carried out numerical simulations and analysis, performed high-speed videography experiments and analysis. M.C.R. performed CO2 absorption experiments and analysis. B.D.M. designed and assembled the Schlieren imaging set-up and participated in data collection. J.T.D. fabricated structures and carried out electroless plating. S.E.B. proposed and supervised the project. E.B.D. conceived the idea, proposed and supervised the project, and wrote the manuscript. All authors participated in discussions of the research and revisions of the manuscript.

Corresponding author

Correspondence to Eric B. Duoss.

Ethics declarations

Competing interests

A patent application has been filed by the Lawrence Livermore National Security LLC on the basis of the cellular fluidics concept described here (US Patent Application Serial No. 16/549,543), on which N.A.D., S.E.B., V. Beck, S. Chandrasekaran, J.R.D., E.B.D., J. Feaster, J. Knipe, J. Mancini, J. Oakdale, F. Qian and M. Worsley are listed as inventors.

Additional information

Peer review information Nature thanks Sung Kang, Abraham D. Stroock and the other, anonymous, reviewer(s) for their contribution to the peer review of this work. Peer reviewer reports are available.

Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Extended data figures and tables

Extended Data Fig. 1 Effects of cellular architecture on fluid transport.

a, Unit cell types used in this study. b, Porosity as a function of relative strut diameter D* can be used to determine the liquid-holding capacity of a given cell type. c, The cell design can be optimized for a desired liquid–solid and gas–liquid interfacial area by choosing an optimal D*. d, Illustration of the repeating structural motif in BCC arrays. Owing to the symmetry of the BCC unit cell, the structure defined by the green control volume (0, 0, 0) is also repeated in the magenta control volume, offset by a half-cell distance (L/2, L/2, L/2). e, Geometric considerations for a simplified analysis of capillary rise in BCC structures. Assuming the force equilibrium always occurs at the central node where Fadh is at a minimum (Fig. 2a), the contact perimeter can be approximated by a circle circumscribing the elliptical cross-section. The capillary height is explicitly solved for from the force equilibrium at the node. For an array of cells, the contact perimeter in the control volume can be approximated by the dashed lines, resulting in a greater contact perimeter and higher capillary rise.

Source data

Extended Data Fig. 2 Numerical simulations of surface-tension-driven and actively pumped liquid flow.

a, Simulations of capillary rise in an array of BCC cells. Similar trends in the spatiotemporal evolution of the gas–liquid and liquid–solid interfaces are observed as for a single column of cells (Fig. 3c). The liquid fronts at the outer corners accelerate and decelerate in an alternating manner. The overall velocity decreases over time. See Supplementary Video 1. b, Simulations of flow through a channel of simple cubic structures of different strut diameters were carried out on the outlined domain of 1.5 mm cells using quarter symmetry. The gas–liquid interface evolution and pressure map along the central cross-section plane are shown at time steps of 2, 10 and 20 ms. The white circles represent the cross-sections of the cell struts. For a cell with 500-μm struts, the gas–liquid interface initially advances, then recedes. c, Flow rate variation over time for cubic cells with different strut diameters. The negative magnitude of the flow rate (black line) corresponds to the receding interface in a cell with 500 μm. d, Average pressure variation over time. In cells with 500-μm struts, the average pressure becomes positive, indicating that the interface will not propagate in the desired direction and suggesting a potential leak scenario. In cells with 600-μm struts, the pressure hovers around zero, indicating that this strut thickness can be considered as the minimum critical value, and that slightly thicker struts are needed to ensure flow in the desired direction. e, The maximum pressure at the interface across the first cell (0 to L) corresponds to the swelling of the liquid in the cell before it propagates into the subsequent cell. Exceeding the Laplace pressure of the advancing interface indicates a leak can occur (Supplementary Fig. 6).

Extended Data Fig. 3 Gas–liquid transport under different environmental conditions.

a, Illustration of water evaporation for cellular fluidics with and without a liquid reservoir and filter paper (left). Evaporative mass loss of the test structures at different relative humidity (middle). When supplied with a liquid reservoir, the evaporation rate of cellular fluidics is comparable to filter paper (right). b, Flow cell used for measurements of humidification and pressure drop for gas flow across wetted devices (left). Comparable levels of humidification are achieved between cellular fluidic devices and filter paper (middle); Data are mean ± s.d.; n = 6. However, significant pressure drops were observed across filter paper, whereas no appreciable pressure drop could be detected for the cellular fluidics. Gas-driven liquid ejection was observed for both the filter paper (11–15 ml min−1) and a cellular fluidic device (>25 ml min−1) with no gas pores (right). Incorporating open (non-liquid-filled) gas pores allowed robust operation with no liquid ejection at all tested flow rates. c, Flow cell used for measurements of pressure drop of liquid flow across a uniform isotruss lattice (left). The data agree well with calculations from the Ergun equation, implying that pressure drop due to inertial effects is significant at higher flow rates (middle). Pressure drops calculated for lattices of different cell types at a liquid velocity of 10 cm s−1 (right). The circle symbol in the right plot corresponds to the circle symbol data point in the middle plot.

Extended Data Fig. 4 Effect of unit cell coordination number on CO2 absorption rate.

a, Images (top) of the hierarchical lattice filled with a solution of Na2CO3 and a pH-sensitive colour indicator that changes from blue to yellow upon CO2 saturation. The rate of absorption in the lattice is significantly faster than a liquid pool of the same sorbent volume. b, Hierarchical structure 750 s after exposure to CO2, with analysed regions of interest outlined. Scale bar, 5 mm. c, Colorimetric map depicting the reaction rate during the first 500 s. The extent of CO2 absorption is estimated by the colour change from blue to yellow (the sum of green and red channels minus the blue channel for each pixel). Higher reaction rates occur at cells with a lower coordination number—that is, more gas–liquid interfaces (blue arrows)—compared to cells with a higher coordination number (black arrows). The enlarged images of the outlined regions show high rates of reaction occurring at interfaces in the centre of each unit cell. Although it is expected that the reaction will occur more rapidly at struts and corners, where the liquid layer is the thinnest and has the highest surface-area-to-volume ratio, the thin liquid layer results in high reflectivity and low levels of pH indicator, which limits the measurement accuracy in these regions. The cell at the intersection with a higher coordination number shows a slower reaction rate than adjacent cells. d, The extent of reaction progression for cells selected in b. Cells with a coordination number 5 (black) show lower levels of reaction than cells with coordination number 2 (blue and pink).

Extended Data Fig. 5 Actively driven flow in channels of simple cubic cells with uniform and graded density structure.

a, In the uniform density structure, the pressure drop across the channel causes the liquid to halt its propagation in the desired direction and break the surface tension at a gas–liquid interface facing an undesirable direction, which results in a leak. b, In the graded density structure, the retaining capillary force along the desired flow direction better matches pressure drop across the channel, which results in uniform and leak-free flow through the entire structure. Scale bars, 5 mm.

Extended Data Fig. 6 Capillary rise in Nx × Ny × Nz arrays of isotruss and simple cubic cells.

a, For isotruss cells, similar to BCC cells (Fig. 2), capillary rise of Nx × Ny × Nz arrays is higher compared to a single column of the same relative density. b, The same effect is not observed in arrays of simple cubic cells, where there is no strut structure that is internal to the cubic cell. Scale bars, 5 mm.

Extended Data Fig. 7 Transpiration experiments comparing cellular fluidic structures with filter paper.

Transpiration was observed by the water level drop in the reservoir and the concentration of green dye at the tip of the structures. a, A cellular fluidic structure and filter paper with an equivalent surface area exhibit similar evaporation rates, plotted as the decrease in liquid level in the reservoir over time. b, Filter paper was cut to match the shape of the cellular fluidic structure, but the filter paper collapsed under its own weight after 90 min of use. Conversely, the cellular fluidic component is a mechanically robust structure that exhibits no sign of collapse or deformation. c, Transpiration from the stochastic tree structure and filter paper holding equivalent volumes of liquid. Scale bars, 5 mm.

Extended Data Table 1 Dimensionless parameters derived for cellular fluidics
Extended Data Table 2 Parameter values of cellular fluidic structures

Supplementary information

Supplementary Information

This file contains Supplementary Methods (analytical calculations) and Supplementary Figures (unmodified images from Figs. 2d, 2f and 6c, evaporation front analysis from Fig. 4c, simulated velocity from Fig. 3c and thermal imaging data from Fig. 4b without smoothing, microscopy of micro-roughness effects, illustration of filter paper testing cell).

Peer Review File

Supplementary Video 1

Capillary rise – simulations and experiments. Simulation (left) and high-speed video (right) of single cell filling and capillary rise in a 1×1×Nz cell column shown in Fig. 3. High-speed video of capillary rise in a column. Simulation of capillary rise in a Nx×Ny×Nz array (Extended Data Fig. 2a). High-speed video of capillary rise in an array.

Supplementary Video 2

Multiphase processes in cellular fluidics. Video of transpiration of green-dyed ethanol through a stochastic tree structure. Infrared videos corresponding to data in Figs. 4b and 4c. Schlieren imaging used to calculate evaporative flux shown in Fig. 4d. Temporal color change of the hierarchical structure is used to estimate CO2 absorption in Extended Data Fig. 4.

Supplementary Video 3

Active flow. Videos of active flow through 1D and 2D cellular fluidic structures, and demonstration of programmed fluidic and selectively metallized pathways in 3D.

Source data

Rights and permissions

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Dudukovic, N.A., Fong, E.J., Gemeda, H.B. et al. Cellular fluidics. Nature 595, 58–65 (2021).

Download citation


By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.


Quick links