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

Bohua Sun; Xiao-Lin Guo

doi:10.20944/preprints202206.0262.v2

Preprint

Article

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

This version is not peer-reviewed.

Bohua Sun^*,

Xiao-Lin Guo

Bohua Sun^*,

Xiao-Lin Guo

This version is not peer-reviewed.

Downloads

487

Views

367

Comments

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:

Dandelion

;

pappus

;

flexible filament

;

wind-dispersal

;

aerodynamic shape

;

drag

;

Reynolds number

;

scaling laws

Subject:

Physical Sciences - Fluids and Plasmas Physics

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 \sim U^{2 + V},

(1)

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] \approx \frac{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} \frac{d^{2} κ}{d s^{2}} + \frac{1}{2} J_{c} κ^{3} = \frac{1}{2} ρ U^{2} d_{c} cos θ

. If introducing

\bar{s} = s / L_{c}

and

\bar{κ} = d θ / d \bar{s} = L_{c} κ

, The ODE and corresponding boundary conditions can be nondimensionalized as to the form:

\frac{d^{2} \bar{κ}}{d {\bar{s}}^{2}} + \frac{1}{2} {\bar{κ}}^{3} = η^{2} cos θ, {\bar{κ}}_{\bar{s} = \frac{1}{2}} = {(\frac{d \bar{κ}}{d \bar{s}})}_{\bar{s} = \frac{1}{2}} = 0,

(2)

where

η = {[\frac{\frac{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] = \frac{1}{2} ρ U^{2}

is the fluid pressure jump across the filament.

The aerodynamic shape of the fibre can be reconstructed by relations:

\frac{d \bar{x}}{d \bar{s}} = cos θ (s)

and

\frac{d \bar{y}}{d \bar{s}} = sin θ (s)

, where

\bar{x} = x / L_{c}

\bar{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}) .

(3)

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}} = \frac{D_{f}}{\frac{1}{2} ρ U^{2} A_{c}}

, Similarly, for

J_{c}

, we have

Π_{J} = \frac{J_{c}}{\frac{1}{2} ρ U^{2} A_{c}^{2}}

and for

ν

, we have

Π_{ν} = \frac{ν}{ρ U \sqrt{A_{c}}}

From Buckingham

Π

theorem [14,15,16,17], the Eq.3 can be equivalency expressed as

Π_{D_{f}} = f (Π_{J}, Π_{ν})

, namely

D_{f} = \frac{1}{2} ρ U^{2} A f (\frac{J_{c}}{\frac{1}{2} ρ U^{2} A_{c}^{2}}, Π_{ν}) .

(4)

For a single characteristic filament, we can propose an approximate

f (\frac{J_{c}}{\frac{1}{2} ρ U^{2} A_{c}^{2}}, Π_{ν}) \approx C_{0} {(\frac{J_{c}}{\frac{1}{2} ρ U^{2} A_{c}^{2}})}^{α} Π_{ν}^{β}

, therefore

D_{f} \approx C_{0} \frac{1}{2} ρ U^{2} A_{c} {(\frac{J_{c}}{\frac{1}{2} ρ U^{2} A_{c}^{2}})}^{α} Π_{ν}^{β},

(5)

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

η \leq 1

, we have

α = 0, β = 0

, for

η \geq 1

, we have

α = \frac{1}{3}, β = 0

, therefore the drag of single filament is given by

D_{f} = \{\begin{matrix} C_{1} \frac{1}{2} ρ L_{c} d_{c} U^{2}, & (η \leq 1), \\ C {[{(\frac{1}{2} ρ)}^{2} J L_{c} d_{c}]}^{1 / 3} U^{4 / 3}, & (η \geq 1) . \end{matrix}

(6)

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

η \geq 1

, and

V = 0

in case of

η \leq 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} = \{\begin{matrix} N C_{1} \frac{1}{2} ρ L_{c} d_{c} U^{2}, & (η \leq 1), \\ N C {[{(\frac{1}{2} ρ)}^{2} J_{c} L_{c} d_{c}]}^{1 / 3} U^{4 / 3}, & (η \geq 1) . \end{matrix}

(7)

Thus we can get the total drag coefficients of the dandelion as follows

C_{D} = \frac{D}{\frac{1}{2} ρ U^{2} N L_{c} d_{c}} = \{\begin{matrix} C_{1}, & (η \leq 1), \\ C {[\frac{J}{\frac{1}{2} ρ U^{2} {(L_{c} d_{c})}^{2}}]}^{1 / 3}, & (η \geq 1) . \end{matrix}

(8)

in which, the case of

η \leq 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 \approx N D_{f}

, which gives the terminal velocity U of the dandelion as follows

U = C_{U} {(m g)}^{3 / 4},

(9)

where the coefficient

C_{U} = {(N C)}^{- 3 / 4} {(\frac{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 = \frac{μ 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 {(\frac{2 ρ J_{c}}{μ^{2} d_{c}^{2}})}^{1 / 3} R e^{- 2 / 3}, (η \geq 1) .

(10)

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} \sim 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 {(\frac{π ρ E d_{c}^{2}}{32 μ^{2}})}^{1 / 3} R e^{- 2 / 3}, (η \geq 1) .

(11)

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} \sim 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} \sim 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

C. Cummins, et al. A separated vortex ring underlies the flight of the dandelion. Nature 562, 414-418 (2018). [CrossRef]
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]
D. F. Greene, The role of abscission in long-distance seed dispersal by the wind. Ecology 86, 3105-3110 (2005). [CrossRef]
D, F. Greene and E. A. Johnson, The aerodynamics of plumed seeds. Funct. Ecol. 4, 117-125(1990). [CrossRef]
S. Vogel, Life in Moving Fluids: The Physical Biology of Flow 2nd edn (Princeton University Press, 2020).
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]
M. C. Andersen, Diaspore morphology and seed dispersal in several wind-dispersed Asteraceae. Am. J. Bot. 80, 487-492(1993). [CrossRef]
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]
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).
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).
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).
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]
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]
E. Buckingham, On physically similar systems: illustration of the use of dimensional equations. Phys. Rev. 1914, 4:345-376. [CrossRef]
P.W. Bridgman, Dimensional Analysis. (Yale University Press, New Haven, 1922).
B.H. Sun, Dimensional Analysis and Lie Group. (China High Education Press, Bejing, 2016).
P.G. de Gennes, Scaling Concepts in Polymer Physics. (Connel Unversity Press, 1979).

Figure 1. Dandelion wind dispersal.

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} = \frac{\sum_{i = 1}^{N} L_{i}}{N}

and a characteristic diameter

d_{c} = \frac{\sum_{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} = \frac{\sum_{i = 1}^{N} L_{i}}{N}

and a characteristic diameter

d_{c} = \frac{\sum_{i = 1}^{N} d_{i}}{N}

, where N is total number of filaments.

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 4. The aerodynamics shape

(\bar{x}, \bar{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

(\bar{x}, \bar{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 5. Terminal velocity vs.

m g

Figure 5. Terminal velocity vs.

m g

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

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

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.

© 2023 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).

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

Subscription

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

MDPI Initiatives

Important Links

Choose an area of interest and we will send you notifications of new preprints at your preferred frequency.

Disclaimer

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

Abstract

Keywords:

Subject:

1. Introduction

2. Aerodynamic Shape of the Dandelion

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

4. Conclusions and Perspectives

Data Availability Statement

Acknowledgments

References

MDPI Initiatives

Important Links

Subscribe