Search

Leaves shed by deciduous trees contain 40\% of the annually sequestered carbon, and include nutrients vital to the expansion and health of forest ecosystems. To achieve this, leaves must fall quickly to land near the parent tree -- otherwise, they are lost to the wind, like pollen or gliding seeds. However, the link between leaf shape and sedimentation speed remains unclear. To gauge the relative performance of extant leaves, we developed an automated sedimentation apparatus (ASAP) capable of performing $\sim100$ free fall experiments per day on biomimetic paper leaves. The majority of 25 representative leaves settle at rates similar to our control (a circular disc). Strikingly, the Arabidopsid mutant asymmetric leaves1 (as1) fell 15\% slower than the wild type. Applying the as1-digital mutation to deciduous tree leaves revealed a similar speed reduction. Data correlating shape and settling across a broad range of natural, mutated, and artificial leaves support thefast-leaf-hypothesis: Deciduous leaves are symmetric and relatively unlobed in part because this maximizes their settling speed and concomitant nutrient retention.

In this study, we investigate the dynamics of particles overcoming the hydrodynamic barrier posed by a thin fluid film to achieve contact in finite time, a phenomenon critical in various natural and engineered processes such as enzyme docking, catalysis, and vesicular transport. Using the framework of lubrication theory, which posits that drag force scales inversely with the film thickness, we explore how particles can achieve finite-time contact despite theoretical predictions of infinite time under constant force. We conduct experiments where a spherical particle settles under gravity and magnetic attraction, the latter introducing a spatially varying force. Our findings reveal that a spatially varying force significantly alters the settling trajectory, enabling finite-time contact. The results are supported by a simple model that links hydrodynamic drag and the impact of spatially varying forces. Finally, we illustrate that forces can be inferred from kinematic observations. In the future, this may provide insights into biological and microscale systems where direct force measurements are challenging. Our study demonstrates that varying forces can be harnessed to overcome lubrication barriers, offering potential applications in designing self-assembly systems and improving surface interaction processes.

Soft intertwined channel systems are frequently found in fluid flow networks in nature. The passage geometry of these systems can deform due to fluid flow, which can cause the relationship between flow rate and pressure drop to deviate from Hagen-Poiseuille's linear law. Although fluid-structure interactions in single deformable channels have been extensively studied, such as in Starling's resistor and its variations, the flow transport capacity of an intertwined channel with multiple self-intersections (a "hydraulic knot"), is still an open question. We present experiments and theory on soft hydraulic knots formed by interlinked microfluidic devices comprising two intersecting channels separated by a thin elastomeric membrane. Our experiments show flow-pressure relationships similar to flow limitation, where the limiting flow rate depends on the knot configuration. To explain our observations, we develop a mathematical model based on lubrication theory coupled with tension-dominated membrane deflections that compares favorably to our experimental data. Finally, we present two potential hydraulic knot applications for microfluidic flow rectification and attenuation.

Exploring fluid-structure interactions is essential for understanding the physical principle underlying flow control in biological and man-made systems. Traditionally, we assume that the geometry is known, and from it, the solution to the coupled elastohydrodynamic problem is determined. Solving the inverse problem -- finding the geometry that leads to a desired flow -- has received comparatively less attention. Here, we present a strategy for solving inverse hydroelastic problems. Specifically, we compute the shape of a soft channel that yields a desired flow-rate versus pressure-drop relationship. The analysis is based on low-Reynolds-number hydrodynamics and linear elasticity. We demonstrate its usefulness in understanding intercellular transport in plants and the design of check valves. The sensitivity of the algorithm to fabrication errors and other limitations are discussed.

Fluid flow networks are ubiquitous and can be found in a broad range of contexts, from human-made systems such as water supply networks to living systems like animal and plant vasculature. In many cases, the elements forming these networks exhibit a highly non-linear pressure-flow relationship. Although we understand how these elements work individually, their collective behavior remains poorly understood. In this work, we combine experiments, theory, and numerical simulations to understand the main mechanisms underlying the collective behavior of soft flow networks with elements that exhibit negative differential resistance. Strikingly, our theoretical analysis and experiments reveal that a minimal network of nonlinear resistors, which we have termed a `fluidic memristor', displays history-dependent resistance. This new class of element can be understood as a collection of hysteresis loops that allows this fluidic system to store information. Our work provides insights that may inform new applications of fluid flow networks in soft materials science, biomedical settings, and soft robotics, and may also motivate new understanding of the flow networks involved in animal and plant physiology.

We describe imbibition in real and artificial plant seeds, using a combination of experiments and theory. In both systems, our experiments demonstrate that liquid permeates the substrate at a rate which decreases gradually over time. Tomographic imaging of soy seeds is used to confirmed this by observation of the permeating liquid using an iodine stain. To rationalize the experimental data, we propose a model based on capillary action which predicts the temporal evolution of the radius of the wet front and the seed mass. The depth of the wetting front initially evolves as $t^{1/2}$ in accord with the Lucas-Washburn law. At later times, when the sphere is almost completely filled, the front radius scales as $(1-t/t_{\text{max}})^{1/2}$ where $t_{\text{max}}$ is the time required to complete imbibition. The data obtained on both natural and artificial seeds collapse onto a single curve that agrees well with our model, suggesting that capillary phenomena contribute to moisture uptake in soy seeds.

Vascular plants rely on differences of osmotic pressure to export sugars from regions of synthesis (mature leaves) to sugar sinks (roots, fruits). In this process, known as M\"unch pressure flow, the loading of sugars from photosynthetic cells to the export conduit (the phloem) is crucial, as it sets the pressure head necessary to power long-distance transport. Whereas most herbaceous plants use active mechanisms to increase phloem concentration above that of the photosynthetic cells, in most tree species, for which transport distances are largest, loading seems to occur via passive symplastic diffusion from the mesophyll to the phloem. Here, we use a synthetic microfluidic model of a passive loader to explore the nonlinear dynamics that arise during export and determine the ability of passive loading to drive long-distance transport. We first demonstrate that in our device, phloem concentration is set by the balance between the resistances to diffusive loading from the source and convective export through the phloem. Convection-limited export corresponds to classical models of M\"unch transport, where phloem concentration is close to that of the source; in contrast, diffusion-limited export leads to small phloem concentrations and weak scaling of flow rates with the hydraulic resistance. We then show that the effective regime of convection-limited export is predominant in plants with large transport resistances and low xylem pressures. Moreover, hydrostatic pressures developed in our synthetic passive loader can reach botanically relevant values as high as 10 bars. We conclude that passive loading is sufficient to drive long-distance transport in large plants, and that trees are well suited to take full advantage of passive phloem loading strategies.

Fluid-structure interactions are ubiquitous in nature and technology. However, the systems are often so complex that numerical simulations or ad hoc assumptions must be used to gain insight into the details of the complex interactions between the fluid and solid mechanics. In this paper, we present experiments and theory on viscous flow in a simple bioinspired soft valve which illustrate essential features of interactions between hydrodynamic and elastic forces at low Reynolds numbers. The setup comprises a sphere connected to a spring located inside a tapering cylindrical channel. The spring is aligned with the central axis of the channel and a pressure drop is applied across the sphere, thus forcing the liquid through the narrow gap between the sphere and the channel walls. The sphere's equilibrium position is determined by a balance between spring and hydrodynamic forces. Since the gap thickness changes with the sphere's position, the system has a pressure-dependent hydraulic resistance. This leads to a non-linear relation between applied pressure and flow rate: flow initially increases with pressure, but decreases when the pressure exceeds a certain critical value as the gap closes. To rationalize these observations, we propose a mathematical model that reduced the complexity of the flow to a two-dimensional lubrication approximation. A closed-form expression for the pressure-drop/flow rate is obtained which reveals that the flow rate $Q$ depends on the pressure drop $\Delta p$, sphere radius $a$, gap thickness $h_0$, and viscosity $\eta$ as $Q\sim \eta^{-1} a^{1/2}h_0^{5/2}\left(\Delta p_c-\Delta p\right)^{5/2}\Delta p$, where the critical pressure $\Delta p_c$ scales with the spring constant $k$ and sphere radius $a$ as $\Delta p_c\sim k a^{-2}$. These predictions compared favorably to the results of our experiments with no free parameters., Comment: 10 pages, 5 figures

Stimuli-responsive materials that modify their shape in response to changes in environmental conditions -- such as solute concentration, temperature, pH, and stress -- are widespread in nature and technology. Applications include micro- and nanoporous materials used in filtration and flow control. The physiochemical mechanisms that induce internal volume modifications have been widely studies. The coupling between induced volume changes and solute transport through porous materials, however, is not well understood. Here, we consider advective and diffusive transport through a small channel linking two large reservoirs. A section of stimulus-responsive material regulates the channel permeability, which is a function of the local solute concentration. We derive an exact solution to the coupled transport problem and demonstrate the existence of a flow regime in which the steady state is reached via a damped oscillation around the equilibrium concentration value. Finally, the feasibility of an experimental observation of the phenomena is discussed. Please note that this version of the paper has not been formally peer reviewed, revised or accepted by a journal.

We study the viscous dissipation in pipe flows in long channels with porous or semipermeable walls, taking into account both the dissipation in the bulk of the channel and in the pores. We give simple closed form expressions for the dissipation in terms of the axially varying flow rate $Q(x)$ and the pressure $p(x)$, generalizing the well known expression $\dot W=Q\,\Delta p$ for the case of impenetrable walls with constant $Q$ and a pressure difference $\Delta p$ between the ends of the pipe. When the pressure $p_0$ outside the pipe is constant, the result is the straightforward generalization $\dot W=\Delta \left[(p-p_0) \,Q\right]$. Finally, applications to osmotic flows are considered., Comment: 6 pages, 2 figures

The phloem vascular system facilitates transport of energy-rich sugar and signaling molecules in plants, thus permitting long range communication within the organism and growth of non-photosynthesizing organs such as roots and fruits. The flow is driven by osmotic pressure, generated by differences in sugar concentration between distal parts of the plant. The phloem is an intricate distribution system, and many questions about its regulation and structural diversity remain unanswered. Here, we investigate the phloem structure in the simplest possible geometry: a linear leaf, found, for example, in the needles of conifer trees. We measure the phloem structure in four tree species representing a diverse set of habitats and needle sizes, from 1 cm (\textit{Picea omorika}) to 35 cm (\textit{Pinus palustris}). We show that the phloem shares common traits across these four species and find that the size of its conductive elements obeys a power law. We present a minimal model that accounts for these common traits and takes into account the transport strategy and natural constraints. This minimal model predicts a power law phloem distribution consistent with transport energy minimization, suggesting that energetics are more important than translocation speed at the leaf level., Comment: 4 figures, supplemental information

We consider the slow flow of a viscous incompressible liquid in a channel of constant but arbitrary cross section shape, driven by non-uniform suction or injection through the porous channel walls. A similarity transformation reduces the Navier-Stokes equations to a set of coupled equations for the velocity potential in two dimensions. When the channel aspect ratio and Reynolds number are both small, the problem reduces to solving the biharmonic equation with constant forcing in two dimensions. With the relevant boundary conditions, determining the velocity field in a porous channels is thus equivalent to solving for the vertical displacement of a simply suspended thin plate under uniform load. This allows us to provide analytic solutions for flow in porous channels whose cross-section is e.g. a rectangle or an equilateral triangle, and provides a general framework for the extension of Berman flow (Journal of Applied Physics 24(9), p. 1232, 1953) to three dimensions.

We perform experimental investigations of osmotically driven flows in artificial microchannels by studying the dynamics and structure of the front of a sugar solution traveling in 200 um wide and 50-200 um deep microchannels. We find that the sugar front travels with constant speed, and that this speed is proportional to the concentration of the sugar solution and inversely proportional to the depth of the channel. We propose a theoretical model, which, in the limit of low axial flow resistance, predicts that the sugar front indeed should travel with a constant velocity. The model also predicts an inverse relationship between the depth of the channel and the speed and a linear relation between the sugar concentration and the speed. We thus find good agreement between the experimental results and the predictions of the model. Our motivation for studying osmotically driven flows is that they are believed to be responsible for the translocation of sugar in plants through the phloem sieve element cells. Also, we suggest that osmotic elements can act as integrated pumps with no movable parts in lab-on-a-chip systems., Comment: 7 Pages, submitted to Lab-on-a-chip on October 24th, 2008

We report a novel and spectacular instability of a fluid surface in a rotating system. In a flow driven by rotating the bottom plate of a partially filled, stationary cylindrical container, the shape of the free surface can spontaneously break the axial symmetry and assume the form of a polygon rotating rigidly with a speed different from that of the plate. With water we have observed polygons with up to 6 corners. It has been known for many years that such flows are prone to symmetry breaking, but apparently the polygonal surface shapes have never been observed. The creation of rotating internal waves in a similar setup was observed for much lower rotation rates, where the free surface remains essentially flat. We speculate that the instability is caused by the strong azimuthal shear due to the stationary walls and that it is triggered by minute wobbling of the rotating plate. The slight asymmetry induces a tendency for mode-locking between the plate and the polygon, where the polygon rotates by one corner for each complete rotation of the plate.