Preprint
Article

Aerodynamic Shape, Scaling Law of Drag and Terminal Speed of the Dandelion

This version is not peer-reviewed.

Submitted:

11 April 2023

Posted:

12 April 2023

You are already at the latest version

A peer-reviewed article of this preprint also exists.

Abstract
The common dandelion uses a bundle of drag-enhancing bristles (the pappus) that enables seed dispersal over formidable distances; however, the scaling laws of aerodynamic drag underpinning pappus-mediated flight remains unresolved. In this paper, we study the aerodynamic shape of dandelion, derive the scaling law of resistance, determine the Vogel exponent. In particular, we find that the total drag coefficient is proportional to the -2/3 power of the dandelion pappus Reynolds number, and obtain the terminal velocity of the dandelion seed under gravitation field.
Keywords: 
;  ;  ;  ;  ;  ;  ;  

1. Introduction

The seeds of plants on Earth, such as dandelion, leave the mother and fly with the wind during the ripening season, and it is a very interesting phenomenon that plants use natural wind to spread seeds. After hundreds of millions of years of natural evolution, different seeds have evolved various but unique structures and thus have their own drag-enhancing pneumatic/aerodynamic behaviors. In order for seeds to spread over long distances, their structure is generally a disc structure composed of flexible filaments. During the process of seed flight, due to the action of hydrodynamic pressure, the flexible filament will deform to form an aerodynamic shape, thereby reducing the flow drag [1,2,3,4,5,6,7,8,9,10,11,12,13].
The common dandelion uses a bundle of drag-enhancing bristles (the pappus) that helps to keep their seeds aloft (as shown in Figure 1). This passive flight mechanism is highly effective, enabling seed dispersal over formidable distances and decreasing its terminal velocity; however, the physics underpinning pappus mediated flight had not been understood until the discovery of the separated vortex ring attached to the pappus [1].
Inspired by the amazing aerodynamics features of the dandelion, very recently, Lyer et al.[9] demonstrated wind-dispersal of battery-free wireless sensing devices.
Whether it is the experimental study of Cummins, et al.[1], or the sensing devices development of Lyer et al.[9], all need to understand the pneumatic/aerodynamic shape and drag behavior of dandelion structures. Regarding the aerodynamic drag of dandelion, their study only gave experimental scatter plots of drag versus Reynolds numbers, and did not obtain a universal scale law; As for the pneumatic/aerodynamic shape of the dandelion, there have been no relevant theoretical and experimental research reports so far.
It is not difficult to understand that the drag of dandelion is closely related to its pneumatic/aerodynamic shape. During the movement of the dandelion, its filaments are pneumatically deformed, so that the dandelion is deformed as a whole. This least drag shape is the aerodynamic shape of a dandelion.
By bending and twisting under fluid loading, on the one hand, dandelion reduce their projected area perpendicular to the flow, and, on the other hand, they also become more streamlined [5,6]. Through these two mechanisms of reconfiguration, the drag load that dandelions must support does not grow with the square of the velocity of the flow they are subjected to – as it would on a rigid bluff body–but rather more slowly. How much slower is described with the Vogel exponent V [5] such that
Drag U 2 + V ,
where U is the flow velocity. However, the exponent V for dandelion has not been determined yet [1,5,6,9].
This article will examine the aerodynamic shape and drag of dandelion. According to the structural characteristics of dandelion, assuming that all filaments are the same, the elastic deformation of each filament is first studied. To simplify, the filament is seen as a free elastic cantilever beam with a fixed bipartisan in the middle, and hydrodynamic pressure acts on the filament. In this way, the aerodynamic shape of the elastic filament is obtained. We regard the overall aerodynamic resistance of the dandelion as the sum of the drag of each filament, and obtain the relationship between the aerodynamic drag coefficient and the dandelion Reynolds number.

2. Aerodynamic Shape of the Dandelion

Dandelions have evolved mechanisms to use wind for seed dispersal over a wide area [1,2,3,4,5] including creating lightweight diaspores with plumose or comose structures that act as drag-enhancing parachutes [7,8,9]. A bundle of flexible filaments bend in a wind in order to absorb deformation energy and reduce the drag force of the wind on the dandelion as shown in Figure 1. The flexible filaments immersed in a flowing medium will adjust its shape to counteract flow drag to minimizes drag by reducing, or delaying, the turbulence in boundary layer of the flow nearest the moving body.
Similar to a flexible fibre modelling in [1,9,13], the dandelion pappus (a bundle of filaments) structure is regarded as a porous disk composed of many flexible filaments as shown in Figure 2.
Consider a flexible filament in the flow medium as shown in Figure 3. The filament is modelled as a thin, inextensible elastic beam loaded by the difference in fluid pressure p between its upstream and downstream sides. The flow characteristic velocity is U and no wake is considered, thus the flexible filament acted fluid pressure [ p ] 1 2 ρ U 2 uniformly.
The Euler-Bernoulli beam under fluid dynamic pressure was formulated by Alben et al. [10,11], Sun and Guo [13], and can be used to dandelion flexible filament. According Sun and Guo [13], leads to a single ordinary differential equation (ODE): J c d 2 κ d s 2 + 1 2 J c κ 3 = 1 2 ρ U 2 d c cos θ . If introducing s ¯ = s / L c and κ ¯ = d θ / d s ¯ = L c κ , The ODE and corresponding boundary conditions can be nondimensionalized as to the form:
d 2 κ ¯ d s ¯ 2 + 1 2 κ ¯ 3 = η 2 cos θ , κ ¯ s ¯ = 1 2 = d κ ¯ d s ¯ s ¯ = 1 2 = 0 ,
where η = 1 2 ρ U 2 L c 2 d c ( J c / L c ) 1 / 2 is called Alben-Shelly-Zhang number [10,13], d is the diameter of the filament, J c = E I c is the bending rigidity, the Young modulus is E, area moment of inertia is I c , [ p ] = 1 2 ρ U 2 is the fluid pressure jump across the filament.
The aerodynamic shape of the fibre can be reconstructed by relations: d x ¯ d s ¯ = cos θ ( s ) and d y ¯ d s ¯ = sin θ ( s ) , where x ¯ = x / L c , y ¯ = y / L c , and shown in Figure 4.

3. Scaling Law of Drag and Terminal Speed of the Dandelion

Similar to the analysis of a flexible fibre in a flowing medium [13], we can also deal with this problem by dimensional method. There are 6 quantities in this problem, namely D , ρ , J c , ν , U , A c , where ν is kinematic viscosity and A = L c d c . The quantity dimensions are listed in the following Table 1:
The drag D f can be expressed as the function of quantities ρ , E , ν , U , A c , namely,
D f = f ( ρ , E , ν , U , A c ) .
In the above relation, there are 6 quantities, two of them are dimensionless. Since only 3 dimensional basis L,M,T are used, so according to dimensional analysis of Buckingham [14], Eq. 3 produces 3 dimensionless quantities Π . The first one is: Π D = D ρ a U b A c = 0 M 0 T 0 ,, we have Π D f = D f 1 2 ρ U 2 A c , Similarly, for J c , we have Π J = J c 1 2 ρ U 2 A c 2 and for ν , we have Π ν = ν ρ U A c .
From Buckingham Π theorem [14,15,16,17], the Eq.3 can be equivalency expressed as Π D f = f ( Π J , Π ν ) , namely
D f = 1 2 ρ U 2 A f ( J c 1 2 ρ U 2 A c 2 , Π ν ) .
For a single characteristic filament, we can propose an approximate f ( J c 1 2 ρ U 2 A c 2 , Π ν ) C 0 J c 1 2 ρ U 2 A c 2 α Π ν β , therefore
D f C 0 1 2 ρ U 2 A c J c 1 2 ρ U 2 A c 2 α Π ν β ,
where C 0 is a constant, α and β are exponents to be determined by experiments.
From Alben, Shelley and Zhang [10,11], Sun and Guo [13], for η 1 , we have α = 0 , β = 0 , for η 1 , we have α = 1 3 , β = 0 , therefore the drag of single filament is given by
D f = C 1 1 2 ρ L c d c U 2 , ( η 1 ) , C ( 1 2 ρ ) 2 J L c d c 1 / 3 U 4 / 3 , ( η 1 ) .
where C 1 , C are a constants.
Comparing with the Vogel law in Eq.1, we can determine the Vogel exponent as follows: V = 2 / 3 in case of η 1 , and V = 0 in case of η 1 .
Since the pappus of the dandelion is consist of N filaments, the total drag of the dandelion is the summation of each filament, therefore we have the total drag of the dandelion
D = N D f = N C 1 1 2 ρ L c d c U 2 , ( η 1 ) , N C ( 1 2 ρ ) 2 J c L c d c 1 / 3 U 4 / 3 , ( η 1 ) .
Thus we can get the total drag coefficients of the dandelion as follows
C D = D 1 2 ρ U 2 N L c d c = C 1 , ( η 1 ) , C J 1 2 ρ U 2 ( L c d c ) 2 1 / 3 , ( η 1 ) .
in which, the case of η 1 is corresponding to the rigid filament. we will not going to use C 1 .
With the total drag in Eq.7, if we assume the weight of total dandelion seed is m g , then we have balance relation m g = D N D f , which gives the terminal velocity U of the dandelion as follows
U = C U ( m g ) 3 / 4 ,
where the coefficient C U = ( N C ) 3 / 4 ( 1 4 ρ 2 J c L c d c ) 1 / 4 . The analytical expression of the terminal velocity has not been seen in literature and whose log profile is depicted in Figure 5.
The Reynolds number is a non-dimensional parameter characterizing the relative importance of inertial to viscous forces in a fluid. The flow through and around the pappus involves two different Reynolds numbers: that of the entire pappus ( R e = ρ U L c / μ , in which U is the velocity of the seed, L c is the characteristic diameter of the pappus and μ the dynamics viscosity of the fluid) and that of an individual filament ( R e f = ρ U d c / μ ). The study of Cummins et al.[1] revealed that the pappus of a dandelion benefits from a ’wall effect’ at low R e f = ρ U d c / μ . Neighbouring filaments interact strongly with one another because of the thick boundary layer around each filament, which causes a considerable reduction in air flow through the pappus.
From the Reynolds numbers of the entire pappus R e = ρ U L c / μ , we have U = μ R e ρ L c . Replacing the velocity U in the 2nd expression of Eq.8, we have the total drag coefficient in terms of entire pappus Reynolds number as follows
C D = C 2 ρ J c μ 2 d c 2 1 / 3 R e 2 / 3 , ( η 1 ) .
The expression reveals that the total drag of the dandelion is proportional to 2 / 3 power law of the pappus Reynolds number, namely C D R e 2 / 3 .
If assume all filaments are circular solid cross-section, then the area moment of inertia I can be calculated as I = π d c 4 / 64 , inserting it to Eq.10, we have
C D = C π ρ E d c 2 32 μ 2 1 / 3 R e 2 / 3 , ( η 1 ) .
To verify our scaling law in Eq.10, we depicted comparisons with Cummins et al. [1] and Lyer, et al. [9] in Figure 6 and Figure 7, respectively.
With the help of our scaling law in Eq.10 and data fitting, we obtain C D = f ( R e ) approximate analytical relations for Cummins et al. [1]. The solid lines are from our formula Eq.10, for blue solid line C D = 350 R e 2 / 3 and for red solid line C D = 250 R e 2 / 3 as shown in Figure 6.
In the same way, we obtain C D = f ( R e ) approximate analytical relations for Lyer et al. [9]. The solid lines are from our formula Eq.10, for green solid line C D = 380 R e 2 / 3 , for blue solid line C D = 295 R e 2 / 3 and for red solid line C D = 250 R e 2 / 3 as shown in Figure 7.

4. Conclusions and Perspectives

Our study obtained not only the aerodynamic/aerodynamic shape of the filament, but also the universal drag scaling laws of the dandelion, as well as the terminal velocity of the dandelion seed under gravitation field.
To the best of the authors’ knowledge, this is the first detailed study of dandelion’s drag in the context of the dimensional analysis. The investigation showns that the drag of the dandelion is proportional to 2 / 3 power law of the pappus Reynolds number, namely C D R e 2 / 3 , which reveals that the softer filament the smaller the drag, which is the secret of reducing drag by fine hairs [18].
For insight perspectives, due to the generality of the scale law we obtain using dimensional analysis, we can even bravely predict that the resistance coefficients of all structures with flexible filaments obey the same law as C D R e 2 / 3 .

Data Availability Statement

The data that support the findings of this study are available from the corresponding author upon reasonable request.

Acknowledgments

This work was supported by Xi’an University of Architecture and Technology (Grant No. 002/2040221134). The authors wish to thank Prof. Hai-Hang Cui from Xi’an University of Architecture and Technology for bringing the dandelion to our attention.

References

  1. C. Cummins, et al. A separated vortex ring underlies the flight of the dandelion. Nature 562, 414-418 (2018). [CrossRef]
  2. D. Lentink, W.B. D. Lentink, W.B. Dickson, J. L. van Leeuwen and M. H. Dickinson, Leading-edge vorticeselevate lift of autorotating plant seeds. Science 324, 1438-1440 (2009). [CrossRef]
  3. D. F. Greene, The role of abscission in long-distance seed dispersal by the wind. Ecology 86, 3105-3110 (2005). [CrossRef]
  4. D, F. Greene and E. A. Johnson, The aerodynamics of plumed seeds. Funct. Ecol. 4, 117-125(1990). [CrossRef]
  5. S. Vogel, Life in Moving Fluids: The Physical Biology of Flow 2nd edn (Princeton University Press, 2020).
  6. F. Gosselin, E. F. Gosselin, E. de Langre and B. A. Machado-Almeida, Drag reduction of flexible plates by reconfiguration, J. Fluid Mech. 650, 319-341(2010). [CrossRef]
  7. M. C. Andersen, Diaspore morphology and seed dispersal in several wind-dispersed Asteraceae. Am. J. Bot. 80, 487-492(1993). [CrossRef]
  8. V. Casseau, G. V. Casseau, G. De Croon, D. Izzo and C. Pandolfi, Morphologic and aerodynamic considerations regarding the plumed seeds of tragopogon pratensis and their implications for seed dispersal. PLoS ONE 10, e0125040 (2015). [CrossRef]
  9. V. Lyer, H. V. Lyer, H. Gaensbauer, T.L. Daniel and S. Gollakota, Wind dispersal of battery-free wireless devices, Nature, 603, 427-433(2022).
  10. S. Alben, M. S. Alben, M. Shelley and J. Zhang, Drag reduction through self-similar bending of a flexible body, Nature 420, 479-481 (2002).
  11. S. Alben, M. S. Alben, M. Shelley and J. Zhang, How flexibility induces streamlining in a two-dimensional flow, Phy. Fluids, 16(5):1694-1713 (2004).
  12. S. Gao, S. S. Gao, S. Pan, H.C. Wang and X.L. Tian, Shape Deformation and Drag Variation of a Coupled Rigid-Flexible System in a Flowing Soap Film, Phy. Rev. Lett. 125, 034502 (2020). [CrossRef]
  13. B.H. Sun and X.L. Guo, Aerodynamic shape and drag scaling law of a flexible fibre in a flowing medium. Theo.Appl.Mech.Lett. (2022) 100397. 100397. [CrossRef]
  14. E. Buckingham, On physically similar systems: illustration of the use of dimensional equations. Phys. Rev. 1914, 4:345-376. [CrossRef]
  15. P.W. Bridgman, Dimensional Analysis. (Yale University Press, New Haven, 1922).
  16. B.H. Sun, Dimensional Analysis and Lie Group. (China High Education Press, Bejing, 2016).
  17. P.G. de Gennes, Scaling Concepts in Polymer Physics. (Connel Unversity Press, 1979).
Figure 1. Dandelion wind dispersal.
Figure 1. Dandelion wind dispersal.
Preprints 70862 g001
Figure 2. The 3D structure pappus is modelled as a planner disk that is consist of filaments with different lengths l i and diameters d i . Transforming the disk to a model with a characteristic length L c = i = 1 N L i N and a characteristic diameter d c = i = 1 N d i N , where N is total number of filaments.
Figure 2. The 3D structure pappus is modelled as a planner disk that is consist of filaments with different lengths l i and diameters d i . Transforming the disk to a model with a characteristic length L c = i = 1 N L i N and a characteristic diameter d c = i = 1 N d i N , where N is total number of filaments.
Preprints 70862 g002
Figure 3. The filament is supported by a thin stainless-steel rod, which is clamped at one end. Fluid drag force acting on the filament deflects this support slightly downwards.
Figure 3. The filament is supported by a thin stainless-steel rod, which is clamped at one end. Fluid drag force acting on the filament deflects this support slightly downwards.
Preprints 70862 g003
Figure 4. The aerodynamics shape ( x ¯ , y ¯ ) for different parameter η = 6 n , n = 1 , 2 , 3 . . . , 8 , the aerodynamic shape is symmetric to the y axis. The aerodynamic shape is a revolution of 2D curve in this figure.
Figure 4. The aerodynamics shape ( x ¯ , y ¯ ) for different parameter η = 6 n , n = 1 , 2 , 3 . . . , 8 , the aerodynamic shape is symmetric to the y axis. The aerodynamic shape is a revolution of 2D curve in this figure.
Preprints 70862 g004
Figure 5. Terminal velocity vs. m g .
Figure 5. Terminal velocity vs. m g .
Preprints 70862 g005
Figure 6. Comparisons with experimental data provided Cummins et al. [1]. The drag coefficient C D for natural (red filled circle) and artificially weighted/clipped (black filled circle) dandelion seeds as a function of R e .
Figure 6. Comparisons with experimental data provided Cummins et al. [1]. The drag coefficient C D for natural (red filled circle) and artificially weighted/clipped (black filled circle) dandelion seeds as a function of R e .
Preprints 70862 g006
Figure 7. Comparisons with experimental data provided Lyer et al. [9]. The drag coefficient C D for disk with film thickness 7.5 μ m (red filled circle), for disk with film thickness 12 μ m (blue filled circle) and for disk with film thickness 25 μ m (green filled circle) dandelion seeds as a function of R e .
Figure 7. Comparisons with experimental data provided Lyer et al. [9]. The drag coefficient C D for disk with film thickness 7.5 μ m (red filled circle), for disk with film thickness 12 μ m (blue filled circle) and for disk with film thickness 25 μ m (green filled circle) dandelion seeds as a function of R e .
Preprints 70862 g007
Table 1. Dimensions of physical quantities.
Table 1. Dimensions of physical quantities.
Variables Notation Dimensions
Filament’s Drag D f M T 2
Mass density ρ M 3
Rigidity J c M 3 T 2
Kinematic viscosity ν 2 t 1
Flow velocity U L 1
Area( = L c d c ) A c 2
The dimensional basis used is length (), mass (M) and time (T).
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.
Copyright: This open access article is published under a Creative Commons CC BY 4.0 license, which permit the free download, distribution, and reuse, provided that the author and preprint are cited in any reuse.

Downloads

487

Views

367

Comments

1

Subscription

Notify me about updates to this article or when a peer-reviewed version is published.

Email

Prerpints.org logo

Preprints.org is a free preprint server supported by MDPI in Basel, Switzerland.

Subscribe

© 2025 MDPI (Basel, Switzerland) unless otherwise stated