Buckling of a beam extruded into highly viscous fluid
Abstract
Inspired by microscopic paramecies which use trichocyst extrusion to propel themselves away from thermal aggressions, we propose a macroscopic experiment to study the stability of a slender beam extruded in a highly viscous fluid. Piano wires were extruded axially at constant speed in a tank filled with corn syrup. The force necessary to extrude the wire was measured to increase linearly at first until the compressive viscous force causes the wire to buckle. A numerical model, coupling a lengthening elastica formulation with resistive force theory, predicts a similar behaviour. The model is used to study the dynamics at large time when the beam is highly deformed. It is found that at large time, a large deformation regime exists in which the force necessary to extrude the beam at constant speed becomes constant and length-independent. With a proper dimensional analysis, the beam can be shown to buckle at a critical length based on the extrusion speed, the bending rigidity and the dynamic viscosity of the fluid. Hypothesising that the trichocysts of paramecies must be sized to maximise their thrust per unit volume as well as avoid buckling instabilities, we predict that their bending rigidity must be about . The verification of this prediction is left for future work.
pacs:
87.16.Qp 47.63.Gd 46.70.De 46.32.+xI Introduction
Self propulsion in micro-organisms is mostly achieved through the beating of flagella or cilia Bray (2000). A widely studied model organism is unicellular Paramecium which swims by beating the cilia covering its body, giving rise to elegant metachronal waves. Although the cell generally employs cilia-based swimming to feed and move, it has been observed recently that when perturbed enough, paramecium can escape with a “jumping” gait Hamel et al. (2011). To avoid a severe thermal aggression, it can extrude trichocysts in the direction of the aggression to burst away rapidly. These tricocysts are rod-like organelles bound to the plasma membrane. Upon activation, these organelles grow about 10-fold in 3 milliseconds before detaching from the animal Sperling et al. (1987); Hamel et al. (2011). The drag force created by these extruding rods pushing against the viscous fluid generates thrust in the opposite direction and enables the Paramecium to accelerate itself up to .
The function of trichocysts is still a subject of investigation Plattner (2002); Buonanno et al. (2013). Their extrusion is known to be a defense mechanism against predators Harumoto and Miyake (1991); Buonanno et al. (2013), and it was shown that explosive local release of trichocysts leads to a cell push back as an escape mechanism Knoll et al. (1991). With their model aggression in the form of local laser heating, Hamel et al. Hamel et al. (2011) convincingly demonstrated through microscope image time sequences, cell tracking and simple fluid mechanics calculations that paramecium achieves emergency propulsion by trichocyst extrusion. It comes as a new development in the understanding of micro-organism locomotion Lauga (2011), and it raises many questions.
Once fully extended, the trichocysts of paramecium are in length and just 200 nm in radius. Considering that trichocysts are left behind after use, they must be regenerated in a process that requires several hours. It is reasonable to expect that these organelles should be optimised to minimise the amount of material that must be left behind to achieve propulsion. Hamel et al. Hamel et al. (2011) suggest that trichocysts are well optimised for low Reynolds number viscous flow: “Because the drag force acting on the trichocysts scales linearly with their length but depends weakly on their width, it is optimal for Paramecium to make the thinnest and longest possible trichocysts.” From a purely fluid-mechanics point of view, this argument holds true. However, trichocysts are very slender and are loaded in axial compression during deployment. Would buckling be a limit to the propulsion burst attainable by trichocyst extrusion? How does buckling occur on an extruding beam?
The extrusion of trichocysts bears resemblance to what is known in dynamics as the inverse spaghetti problem Mansfield and Simmonds (1987); Païdoussis (2003). This problem is encountered in the ejection of paper from a copy machine Mansfield and Simmonds (1987); Hong et al. (2007), morphing aircrafts with deploying wings Zhang et al. (2013), deployment of spacecraft structures such as solar arrays or antennas Downer and Park (1993), extrusion fabrication processes, and robotic manipulators Païdoussis (2003). In these different applications, inertia, gravitational, and aero/hydrodynamic forces can come into play. Most relevant to the problem at hand is the extrusion of a beam in a dense fluid at large Reynolds number Taleb and Misra (1981); Gosselin et al. (2007). Viscous friction acting against the extrusion can compress the beam and make it unstable, leading to large transverse vibrations. However, in these linear stability studies Taleb and Misra (1981); Gosselin et al. (2007) significant inertia is considered leading to a vibratory dynamics significantly different from the overdamped one expected in the present study at low Reynolds number.
Experimental results on the inverse spaghetti problem are few and limited to some work on paper ejection Hong et al. (2007). Research has been conducted almost exclusively through theoretical models. The inverse spaghetti problem is non-autonomous (i.e., time appears explicitly in the governing equations), and its solution may be sought by time-integration techniques. The dynamics of a deploying Euler-Bernoulli beam has been considered under small displacement assumption Taleb and Misra (1981); Gosselin et al. (2007), with non-linear models coupling axial and transverse deformations Park et al. (2013); Zhang et al. (2013), and also using elastica models for arbitrary large amplitude Mansfield and Simmonds (1987); Downer and Park (1993); Stolte and Benson (1993). In the studies considering the deployment of a beam in a stagnant fluid Taleb and Misra (1981); Gosselin et al. (2007), the fluid forces on the beam are modelled using a local approach combining an analytically calculated inviscid added mass with semi-empirical pressure drag and skin friction formulations based on the equations of Taylor Taylor (1952). More generally, this simplified approach has been proven succesful for studying the dynamics of beams in large Reynolds number axial flow (Païdoussis (2003, Chapters 8 and 9)).
At low Reynolds number, resistive-force theory offers a conceptually similar local drag force model for Stokes flow Gray and Hancock (1955); Lauga and Powers (2009); Powers (2010). In the limit of an asymptotically thin filament, resistive-force theory allows evaluating the hydrodynamic force per unit length acting at a point along the filament independently of the rest of the filament. Combined with theoretical modelling, macroscopic scale experimental models of the swimming of micro-organisms have been employed to study the propulsive efficiency of actuated passive rods Koehler and Powers (2000); Kim et al. (2003); S. Tony and Hosoi (2006); Qian et al. (2008); Coq et al. (2008); Rodenborn et al. (2013). These macro-scale experiments allow elegant demonstration of the physics of micro-organisms as well as a validation of the theoretical fluid-structure interactions models.
In this paper, we present a macro scale experiment inspired by the extrusion of trichocysts by Paramecium. We study the stability of slender steel rods extruded in corn syrup. The high viscosity of the fluid makes the Reynolds number very small and the inertial effects negligible. We also propose a theoretical model coupling the 2D formulation of a lengthening elastica with the resistive-force theory. This combination of experiments and theory allows understanding how flexibility influences the propulsive force which can be generated by beam extrusion in a viscous fluid as well as investigate its post-buckling dynamics.
Ii Experiment
Steel piano-wires of diameter mm and length mm were extruded vertically into a cubic aquarium with sides of mm filled of highly viscous corn syrup (Fig. 1). The bending modulus of the wires, , was determined by hanging weights to their centre and measuring the deflection. The viscosity of corn syrup was found to be fairly constant for sheer rates of to in a rheology test conducted with a rotational rheometer. The density of corn syrup was measured to be . The wire is pushed by a MIG welder wire-drive assembly consisting of a DC motor that drives two counter rotating wheels (USA Weld, Elk Grove Village, IL). A 3D printed funnel with a 0.5 mm outlet guides the wire vertically down to the fluid surface. The distance of the outlet of the tip to the motor wheels is approximately mm, so that the total extrusion length . The vertical force, , necessary to extrude the wire is measured using a Tedea-Huntleigh 1022 load cell (Vishay Precision Group, Malvern, PA) connected to a 10-bit U3-LV data-acquisition-system (LabJack, Lakewood, CO). A mirror is placed next to the aquarium so that a camera can capture orthogonal views simultaneously.
For each run, the wire is loaded in the feeding mechanism with its tip just above the fluid surface. When the power supply is switched on, the wire is extruded at a constant velocity ranging in equidistant steps from up to . At least 5 tests were conducted for each of the 7 extruding velocity.
Based on the total extruded length, depending on the extruding velocity, the Reynolds number is found to be . Based on the diameter of the wire, then the value is bound by . The two corresponding values of Reynolds number in the case of real trichocysts are and based on the numbers in Ref. Hamel et al. (2011). The Reynolds similitude is satisfactory and we expect the fluid flow in the experiment to be Stokesian.
A representative time sequence of the wire extrusion at a fast speed is shown in Fig. 2 (a). See Supplemental Material at [URL will be inserted by publisher] for the raw video footage of this extrusion. At first, the wire comes out of the nozzle straight. On the frame, the wire starts deforming significantly. In the subsequent frames, a distinctive S-mode of deformation becomes more and more pronounced. Note that from inspection of the front and mirror views (not shown), the deformation is largely two-dimensional.
Force measurements obtained from averaging 5 experimental runs for each extrusion velocity are shown in Fig. 3 (a). For the two slowest extrusion speeds, the curves display linear increase in the extrusion force with time until the extrusion is stopped. For velocities and higher, the extrusion force increases linearly initially until it reaches a maximum value and decreases thereafter. The faster the extrusion speed is, the steeper the initial rise in extrusion force. Also, note that for the fastest extrusion speed , the initial trend does not appear linear because the acceleration of the feeder motor becomes non-negligible.
From both perpendicular video recordings, image processing in Matlab is used to find two functions describing the instantaneous position of the wire in the and planes, respectively. It is based on the minimum light intensity for every value along the vertical. Each of the two generated functions is then least-square fitted on an expansion of 10 cantilevered beam modes. The curvature of these fitted functions in both planes is used to compute the elastic energy of the beam from its clamped end to its free end , i.e., Landau and Lifshitz (1970). The evolution with time of the bending energy computed this way and averaged over 5 experimental runs is shown in Fig. 3 (b) for the 7 extrusion velocities. For all extrusion velocities, the bending energy is small initially then rises suddenly, indicating a loss of stability by buckling. For faster extrusion speeds, this instability occurs earlier and the rise in bending energy is faster.
To make the considered quantities dimensionless, we introduce a ‘bending length’ Alben et al. (2002), which captures the competition between the wire rigidity and longitudinal viscous fluid forcing:
(1) |
where is the dimensional longitudinal friction coefficient. It is expected that this coefficient vary as Powers (2010), and thus should vary during the extrusion. However, it can be shown that this logarithmic variation has minimal influence on the estimation of the force in the problem at hand
The experimental data of Fig. 3 are plotted in dimensionless form in Fig. 4. All the force measurements collapse onto a master curve in Fig. 4 (a). For time , the dimensionless extrusion force increases linearly, and the bending energy in Fig. 4 (b) is small indicating that the wire is straight. At dimensionless time between 6 and 8, the extrusion force peaks and the bending energy starts increasing.
Iii Theoretical Model
For a deeper insight and to consider longer extrusions, we propose a theoretical model coupling a lengtening elastica with a fluid force formulation based on resistive force theory.
Consider a beam lengthening at a constant speed such that at any instant , its length is given by . The deformation of the beam is described by the local angle or curvature . The centreline of the beam is traced by the curvilinear coordinate defined from its clamped end to its free end . The bending and compressive forces in the beam are balanced by the fluid viscous forces
(2) |
where is the bending rigidity, the axial compressive force, the fluid forces acting on the beam, the unit tangent vector to the beam, and the unit vector normal to the beam.
For a very slender beam at low Reynolds number, the fluid forces acting on the beam can be approximated using the resistive force theory Gray and Hancock (1955); Lauga and Powers (2009); Powers (2010). The viscous force per unit length acting at one point on the beam is given by
(3) |
where and are respectively the local tangential and normal velocity components of the beam, and where and are the friction coefficients. For a cylinder in an infinite fluid domain, these friction coefficients can be evaluated using analytical expressions of slender-body hydrodynamics Keller and Rubinow (1976); Gray and Hancock (1955). Here, however, we evaluate experimentally (described below Eq. 1) and make use of the analytical relationship of , true for cylinders of infinite aspect ratio.
A material point on the beam moves along the centreline at the extruding velocity since we assume the beam to be inextensible, i.e., . With the material derivative, the tangential and normal velocity components can be reconstructed from the Eulerian velocity components
(4) | |||||
(5) |
where and are the eulerian coordinates aligned with the tangential and normal directions when the beam is undeformed.
Solution of the beam Eq. 2 with the fluid force of Eq. 3 is obtained by integrating the spatial domain with finite differences (function bvp4c in Matlab using default adaptive mesh size). The Eulerian coordinates are integrated in time using the backwards Euler scheme using time steps of . The model has only three parameters: we fix , and impose an initial perturbation from the straight equilibrium position in the form of a uniform curvature at time =1.
For a small initial uniform curvature, , the computed shape of the beam as it is extruded is shown in Fig. 2 (b). The similarity with the experimental measurements is appreciable.
The variation with time of the dimensionless extrusion force and bending energy predicted by the model are shown in Fig. 4 (a) and (b), respectively. The dynamics of the beam is simulated for a small and a large uniform initial curvature. For , the extrusion force increases linearly at first, then reaches a maximum force of 8.7 at time 8.8 before decreasing and oscillating lightly about a constant value of 3.7 for times larger than 13. The bending energy is small initially and raises abruptly when the extrusion force peaks. For the initial curvature , the behaviour is the same, but a lower force peak of 6.3 is reached earlier at time 6.8. The model agrees well with the experiments: in both, an initial linear force increase is interrupted when the length reaches a critical value of . It is informative to compare this prediction of the dynamical non-linear model with linear static theory. A cantilevered beam under a static distributed compressive axial load (i.e. a standing beam subjected to its own weight) has a critical length beyond which it becomes linearly unstable Timoshenko and Gere (1963). This discrepancy is attributable to the fact that upon buckling, perturbations have a finite growth rate and it takes time for the beam to achieve significant deformation when moving in a viscous fluid. Moreover, as the beam keeps lengthening, the bending energy already accumulated gets redistributed in more material.
The non-linear model is used to study the dynamics at large time. Snapshots of the beam shape at incrementing time steps are shown in the inset in Fig. 5. See Supplemental Material at [URL will be inserted by publisher] for a video animation of this extrusion. At large time, the beam becomes more and more compact. Its deformation shares visual similarities with the inflexional elastica solutions for the case of a rod subjected to terminal forces alone, without couple Love (1944); Timoshenko and Gere (1963). The time history of the beam curvature at the base (extrusion point), shown in Fig. 6, indicates an oscillatory variation. Similarly to the folding of a falling viscous sheet Ribe (2003) or a viscous thread flowing along a less viscous liquid Cubaud and Mason (2006), the buckling of a beam extruded in a viscous liquid gives rise to an oscillatory motion. The frequency of oscillation is proportional to the number of bending lengths extruded per second, i.e., . Although buckling is usually viewed as a static instability, here it gives rise to a dynamic oscillatory motion.
The contour plot shown in Fig 5 demonstrates the evolution of the beam curvature as the extrusion time and thus the beam length increase. In this figure, dashed lines represent the loci of inflection (zero curvature) points, and the regions bounded by these lines are lobes of either positive or negative curvature. It is seen that, after some initial transient, curvature is generated close to the extrusion point and then it moves along the beam as a diminishing-amplitude travelling wave. In other words, after a lobe is generated, it moves along the beam while its curvature reduces and finally becomes flat (destroyed) at the beam tip. In terms of speed, the lobes accelerate from a small velocity near the base () to an almost constant velocity for most their journey and then decelerate as approaching the tip.
For times larger than and up to the largest simulated time of 70, the extrusion force is constant with only small fluctuations. This means that once the beam is buckled, assuming that its material still behaves elastically, the force required to extrude this forever lengthening beam is constant. This counterintuitive result can be obtained by dimensional analysis following a reasoning similar to Gosselin et al. (2010): at large time, when the beam is buckled and completely deformed (see inset in Fig. 5), its dimension becomes irrelevant. Thus the only four quantities remaining can be combined into a single dimensionless number required to describe this physical problem at large time, i.e., or equivalently . For a given beam with a fixed bending length and a fixed bending rigidity, the extrusion force must be constant, irrespective of the instantaneous length of the beam.
It is insightful to compare the behaviour observed in this version of the inverse spaghetti problem at low Reynolds number with the vibrational dynamics simulated in the linear stability analysis of a beam deployed in a dense fluid where inertia is significant Taleb and Misra (1981); Gosselin et al. (2007). In both problems, viscous forces can lead to a buckling instability. We have described how this instability influences the extrusion force in the low Reynolds case in Fig. 4, however it is not clear how this influence is felt in the problem with inertia because the linear stability analysis employed in these studies Taleb and Misra (1981); Gosselin et al. (2007) does not offer this insight. Inertia greatly influences the dynamics of the deploying beam, making it vibrational whereas here the large viscosity leads to an overdamped motion at small time. Interestingly though, an oscillating motion is observed at large time when the extruding beam is highly deformed (see Fig. 6). Another key difference between both problems is that the model with inertia depends on many different parameters and coefficients, some of which are difficult to evaluate. Here, while it is understood that we made simplifying assumptions regarding the fluid forces, there is basically only the initial conditions which can be varied. And even those do not influence the dynamics tremendously (see Fig. 4). When comparing with the complex dynamics of other variants of the inverse spaghetti problem Mansfield and Simmonds (1987); Païdoussis (2003); Hong et al. (2007); Zhang et al. (2013); Downer and Park (1993); Taleb and Misra (1981); Gosselin et al. (2007), the relative simplicity of the dynamics observed for the extruding beam at low Reynolds number is in itself a meaningful result. It is this relative simplicity which allows making predictions about the mechanical nature of trichocysts in the following section.
Iv Trichocyst Extrusion
Lastly, we consider how this study applies to trichocyst extrusion. To use material as efficiently as possible because tricocysts must be regenerated after use, we hypothesise that paramecies have evolved trichocysts sized to be as slender as possible yet respecting the critical length of . From Hamel et al. (2011) we can evaluate , and . However, to the best of our knowledge, no measurement of the bending rigidity of trichocysts exists in the literature. This makes it impossible to check if trichocysts do respect the critical length. However, if indeed they do, we can predict that their bending rigidity should be about . To put this number into perspective, the bending rigidity of Paramecium cilia is Hill et al. (2010). If our prediction is correct, trichocysts must be 5 times more rigid in bending than cilia.
Inherent to this calculation of the bending rigidity of trichocysts is the approximation of the fluid forces using resistive force theory where we consider a rod moving in an infinite domain of fluid a rest. It is probable that the presence and the motion of the paramecium cell, and the fact that many trichocysts are extruded in parallel should result in trichocysts perceiving a more complex flow than is considered here. The importance of these effects is hard to estimate without involved CFD analyses. Considering that this resistive force approach was shown to correctly estimate the thrust force of trichocyst extrusion on a swimming paramecium Hamel et al. (2011), we can expect it to approximate reasonnably well the internal loads carried by the trichocyst which leads to its instability.
V Conclusion
In conclusion, a lengthening elastica model coupled with resistive force theory successfully predicts the time-dependent buckling phenomenon observed when a piano wire is extruded in corn syrup. This combination of theory and experiment allowed us to highlight the following aspects of the dynamics of a beam extruded in highly viscous fluid: buckling limits the amount of thrust achievable through beam extrusion; a post-buckling large-deformation regime exists where the force required to extrude a beam at constant speed is constant and independent of the instantaneous length of the beam; in the post-buckling regime, a self-excited oscillatory motion is generated, sending travelling waves from the root to the tip of the beam; the bending length is the governing length-scale of the problem which allows collapsing the experimental measurements of the extrusion force and which scales the problem at large time.
Moreover, because we hypothesise that trichocysts must be optimised to be as slender as possible while avoiding large amplitude buckling during their ejection, we predict that their bending rigidity should be about . Future work should test this prediction, possibly using an atomic force microscope. Confirmation of this prediction would be an additional clue that tricocysts evolved for propulsion.
Acknowledgements
The authors are thankful to Sébastien Berger and Ryan Mathews for their design work on preliminary versions of the experimental set-up and to Emmanuel de Langre for enlightening discussions. The financial support of the Natural Sciences and Engineering Research Council of Canada and the Fonds de recherche du Québec - Nature et technologies is acknowledged.
Footnotes
- thanks: Corresponding author. Electronic mail: frederick.gosselin@polymtl.ca.
- From slender body hydrodynamics, it is expected that the viscous force on the cylinder obeys a function of the form . This function can be approximated with a linear fit, which over a range , gives a correlation coefficient . For our application, an approximation of the viscous force varying linearly with is accurate enough.
References
- D. Bray, Cell movements (Taylor & Francis, 2000).
- A. Hamel, C. Fisch, L. Combettes, P. Dupuis-Williams, and C. N. Baroud, PNAS 108, 7290 (2011).
- L. Sperling, A. Tardieu, and T. Gulik-Krzywicki, The Journal of cell biology 105, 1649 (1987).
- H. Plattner, Bioessays 24, 649 (2002).
- F. Buonanno, T. Harumoto, and C. Ortenzi, Zoological science 30, 255 (2013).
- T. Harumoto and A. Miyake, Journal of Experimental Zoology 260, 84 (1991).
- G. Knoll, B. Haacke-Bell, and H. Plattner, European journal of protistology 27, 381 (1991).
- E. Lauga, Proceedings of the National Academy of Sciences 108, 7655 (2011).
- L. Mansfield and J. Simmonds, Journal of applied mechanics 54, 147 (1987).
- M. P. Païdoussis, Fluidï¿½ï¿½-Structure Interactions: Slender Structures and Axial Flow, Vol. 2 (Elsevier Academic Press, London, 2003).
- S.-K. Hong, J.-G. Jee, Y. H. Jang, N.-C. Park, and Y.-P. Park, Microsystem technologies 13, 1141 (2007).
- W. Zhang, L. Sun, X. Yang, and P. Jia, Journal of Sound and Vibration 332, 6785 (2013).
- J. D. Downer and K. Park, AIAA journal 31, 339 (1993).
- I. Taleb and A. Misra, Journal of Hydronautics 15, 62 (1981).
- F. Gosselin, M. Païdoussis, and A. Misra, Journal of sound and vibration 299, 123 (2007).
- S. Park, H. H. Yoo, and J. Chung, AIAA journal 51, 686 (2013).
- J. Stolte and R. C. Benson, Journal of vibration and acoustics 115, 308 (1993).
- G. Taylor, Proceedings of the Royal Society of London. Series A. Mathematical and Physical Sciences 214, 158 (1952).
- J. Gray and G. Hancock, Journal of Experimental Biology 32, 802 (1955).
- E. Lauga and T. R. Powers, Rep. Prog. Phys. 72 (2009).
- T. R. Powers, Reviews of Modern Physics 82, 1607 (2010).
- S. A. Koehler and T. R. Powers, Phys. Rev. Lett. 85, 4827 (2000).
- M. Kim, J. C. Bird, A. J. V. Parys, K. S. Breuer, and T. R. Powers, PNAS 100, 15481 (2003).
- E. L. S. Tony and A. E. Hosoi, Phys. Fluids 18 (2006).
- B. Qian, T. R. Powers, and K. S. Breuer, Phys. Rev. Lett. 100 (2008).
- N. Coq, O. du Roure, J. Mathelot, D. Bartolo, and M. Fermigier, Phys. Fluids 20 (2008).
- B. Rodenborn, C.-H. Chen, H. L. Swinney, B. Liu, and H. Zhang, Proceedings of the National Academy of Sciences 110, E338 (2013).
- L. D. Landau and E. M. Lifshitz, Theory of elasticity (Pergamon Press, 1970).
- S. Alben, M. Shelley, and J. Zhang, Nature 420, 479 (2002).
- From slender body hydrodynamics, it is expected that the viscous force on the cylinder obeys a function of the form . This function can be approximated with a linear fit, which over a range , gives a correlation coefficient . For our application, an approximation of the viscous force varying linearly with is accurate enough.
- J. B. Keller and S. I. Rubinow, J. Fluid Mech. 75, 705 (1976).
- S. P. Timoshenko and J. M. Gere, Theory of elastic stability, Second edition (McGraw-Hill, 1963).
- A. Love, A treatise on the mathematical theory of elasticity (Dover, 1944).
- N. M. Ribe, Phys. Rev. E 68, 036305 (2003).
- T. Cubaud and T. G. Mason, Phys. Rev. Lett. 96, 114501 (2006).
- F. Gosselin, E. de Langre, and B. A. Machado-Almeida, Journal of Fluid Mechanics 650, 319 (2010).
- D. B. Hill, V. Swaminathan, A. Estes, J. Cribb, E. T. O’Brien, C. W. Davis, and R. Superfine, Biophysical journal 98, 57 (2010).