skip to main content


Title: Effects of body shape on aerodynamic performance and wake structures in snake-like gliding
Flying snakes are the only snakes on Earth capable of aerial gliding, taking advantage of fluid dynamic principles to leap from point to point among the trees. During their gliding, the locomotion of aerial undulation is observed. We hypothesize that this locomotion and its associated unsteady vortex dynamics are critical to their aerodynamic performance. However, there is a lack of detailed three-dimensional flow field information around the snake body in gliding due to the difficulties in experimental flow visualizations of live animals. In this study, a computation fluid dynamics (CFD) study has been conducted to study the fluid dynamics of a snake-like gliding. A mathematical equation describing the horizontal undulation motion was applied for constructing snake-like 3D computational models and a series of flow simulations were conducted. An immersed-boundary-method (IBM)-based direct numerical simulation (DNS) flow solver along with adaptive mesh refinement (AMR) was used in the simulation. Specifically, different head positions, corresponding to different horizontal wave shapes and their effect on aerodynamic performance, flow field and wake structures behind the body will be studied. In addition, the dynamic undulating motion is introduced in the model and a CFD simulation is also conducted. Results from this study are expected to bring a step stone to understanding snake-inspired locomotion.  more » « less
Award ID(s):
2027534
NSF-PAR ID:
10470900
Author(s) / Creator(s):
Publisher / Repository:
AIAA
Date Published:
Format(s):
Medium: X
Location:
Orlando, FL
Sponsoring Org:
National Science Foundation
More Like this
  1. This paper numerically studies the flow dynamics of aerial undulation of a snake-like model, which is adapted from the kinematics of the flying snake (Chrysopelea) undergoing a gliding process. The model applies aerial undulation periodically in a horizontal plane where a range of angle of attack (AOA) is assigned to model the real gliding motion. The flow is simulated using an immersed-boundary-method-based incompressible flow solver. Local mesh refinement mesh blocks are implemented to ensure the grid resolutions around the moving body. Results show that the undulating body produces the maximum lift at 45° of AOA. Vortex dynamics analysis has revealed a series of vortex structures including leading-edge vortices (LEV), trailing-edge vortices, and tip vortices around the body. Changes in other key parameters including the undulation frequency and Reynolds number are also found to affect the aerodynamics of the studied snake-like model, where increasing of undulation frequency enhances vortex steadiness and increasing of Reynolds number enhances lift production due to the strengthened LEVs. This study represents the first study of both the aerodynamics of the whole body of the snake as well as its undulatory motion, providing a new basis for investigating the mechanics of elongated flexible flyers.

     
    more » « less
  2. null (Ed.)
    ABSTRACT Flying snakes flatten their body to form a roughly triangular cross-sectional shape, enabling lift production and horizontal acceleration. While gliding, they also assume an S-shaped posture, which could promote aerodynamic interactions between the fore and the aft body. Such interactions have been studied experimentally; however, very coarse models of the snake's cross-sectional shape were used, and the effects were measured only for the downstream model. In this study, the aerodynamic interactions resulting from the snake's posture were approximated using two-dimensional anatomically accurate airfoils positioned in tandem to mimic the snake's geometry during flight. Load cells were used to measure the lift and drag forces, and flow field data were obtained using digital particle image velocimetry (DPIV). The results showed a strong dependence of the aerodynamic performance on the tandem arrangement, with the lift coefficients being generally more influenced than the drag coefficients. Flow field data revealed that the tandem arrangement modified the separated flow and the wake size, and enhanced the lift in cases in which the wake vortices formed closer to the models, producing suction on the dorsal surface. The downforce created by the flow separation from the ventral surface of the models at 0 deg angle of attack was another significant factor contributing to lift production. A number of cases showing large variations of aerodynamic performance included configurations close to the most probable posture of airborne flying snakes, suggesting that small postural variations could be used to control the glide trajectory. 
    more » « less
  3. Abstract

    Flow resistance in mountain streams is important for assessing flooding hazard and quantifying sediment transport and bedrock incision in upland landscapes. In such settings, flow resistance is sensitive to grain‐scale roughness, which has traditionally been characterized by particle size distributions derived from laborious point counts of streambed sediment. Developing a general framework for rapid quantification of resistance in mountain streams is still a challenge. Here we present a semi‐automated workflow that combines millimeter‐ to centimeter‐scale structure‐from‐motion (SfM) photogrammetry surveys of bed topography and computational fluid dynamics (CFD) simulations to better evaluate surface roughness and rapidly quantify flow resistance in mountain streams. The workflow was applied to three field sites of gravel, cobble, and boulder‐bedded channels with a wide range of grain size, sorting, and shape. Large‐eddy simulations with body‐fitted meshes generated from SfM photogrammetry‐derived surfaces were performed to quantify flow resistance. The analysis of bed microtopography using a second‐order structure function identified three scaling regimes that corresponded to important roughness length scales and surface complexity contributing to flow resistance. The standard deviationσzof detrended streambed elevation normalized by water depth, as a proxy for the vertical roughness length scale, emerges as the primary control on flow resistance and is furthermore tied to the characteristic length scale of rough surface‐generated vortices. Horizontal length scales and surface complexity are secondary controls on flow resistance. A new resistance predictor linking water depth and vertical roughness scale, i.e. H/σz, is proposed based on the comparison betweenσzand the characteristic length scale of vortex shedding. In addition, representing streambeds using digital elevation models (DEM) is appropriate for well‐sorted streambeds, but not for poorly sorted ones under shallow and medium flow depth conditions due to the missing local overhanging features captured by fully 3D meshes which modulate local pressure gradient and thus bulk flow separation and pressure distribution. An appraisal of the mesh resolution effect on flow resistance shows that the SfM photogrammetry data resolution and the optimal CFD mesh size should be about 1/7 to 1/14 of the standard deviation of bed elevation. © 2019 John Wiley & Sons, Ltd.

     
    more » « less
  4. The rising global trend to reduce dependence on fossil fuels has provided significant motivation toward the development of alternative energy conversion methods and new technologies to improve their efficiency. Recently, oscillating energy harvesters have shown promise as highly efficient and scalable turbines, which can be implemented in areas where traditional energy extraction and conversion are either unfeasible or cost prohibitive. Although such devices are quickly gaining popularity, there remain a number of hurdles in the understanding of their underlying fluid dynamics phenomena. The ability to achieve high efficiency power output from oscillating airfoil energy harvesters requires exploitation of the complexities of the event of dynamic stall. During dynamic stall, the oncoming flow separates at the leading edge of the airfoil to form leading ledge vortex (LEV) structures. While it is well known that LEVs play a significant role in aerodynamic force generation in unsteady animal flight (e.g. insects and birds), there is still a need to further understand their spatiotemporal evolution in order to design more effective energy harvesting enhancement mechanisms. In this work, we conduct extensive experimental investigations to shed-light on the flow physics of a heaving and pitching airfoil energy harvester operating at reduced frequencies of k = fc=U1 = 0.06-0.18, pitching amplitude of 0 = 75 and heaving amplitude of h0 = 0:6c. The experimental work involves the use of two-component particle image velocimetry (PIV) measurements conducted in a wind tunnel facility at Oregon State University. Velocity fields obtained from the PIV measurements are analyzed qualitatively and quantitatively to provide a description of the dynamics of LEVs and other flow structures that may be present during dynamic stall. Due to the difficulties of accurately measuring aerodynamic forces in highly unsteady flows in wind tunnels, a reduced-order model based on the vortex-impulse theory is proposed for estimating the aerodynamic loadings and power output using flow field data. The reduced-order model is shown to be dominated by two terms that have a clear physical interpretation: (i) the time rate of change of the impulse of vortical structures and (ii) the Kutta-Joukowski force which indirectly represents the history effect of vortex shedding in the far wake. Furthermore, the effects of a bio-inspired flow control mechanism based on deforming airfoil surfaces on the flow dynamics and energy harvesting performance are investigated. The results show that the aerodynamic loadings, and hence power output, are highly dependent on the formation, growth rate, trajectory and detachment of the LEV. It is shown that the energy harvesting efficiency increases with increasing reduced frequency, peaking at 25% when k = 0.14, agreeing very well with published numerical results. At this optimal reduced frequency, the time scales of the LEV evolution and airfoil kinematics are matched, resulting in highly correlated aerodynamic load generation and airfoil motion. When operating at k > 0:14, it is shown that the aerodynamic moment and airfoil pitching motion become negatively correlated and as a result, the energy harvesting performance is deteriorated. Furthermore, by using a deforming airfoil surface at the leading and trailing edges, the peak energy harvesting efficiency is shown to increase by approximately 17% and 25% relative to the rigid airfoil, respectively. The performance enhancement is associated with enhanced aerodynamic forces for both the deforming leading and trailing edges. In addition, The deforming trailing edge airfoil is shown to enhance the correlation between the aerodynamic moment and pitching motion at higher reduced frequencies, resulting in a peak efficiency at k = 0:18 as opposed to k = 0:14 for the rigid airfoil. 
    more » « less
  5. Road accidents caused by heavy rain have become a frightening issue in recent years requiring investigation. In this regard, an aerodynamic comparative and experimental rain study is carried out to observe the flow phenomena change around a generic ground vehicle (Ahmed Body at a scale) and the utility truck. In this paper, a Discrete Phase Model (DPM) based computational methodology is used to estimate the effect of rain on aerodynamic performance. First, an experimental rain study of the Ahmed body at a scale that is representative of a car or light truck was conducted at the Wall of Wind (WOW) large-scale testing facility using force measurement equipment. In addition, the experiment allowed drag, lift, and side-force coefficients to be measured at yaw angles up to 55 degrees. Next, experimental results are presented for the Ahmed Body back angle of 35 degrees, then compared to validate the computational model for ground vehicle aerodynamics. Afterwards, we investigated the effect of heavy rainfall (LWC = 30 g/m3) on the external aerodynamics of the utility truck with the morphing boom equipment using the validated computational fluid dynamics method, and the external flow is presented using a computer visualization. Finally, force & moment coefficients and velocity distributions around the utility truck are computed for each case, and the results are compared. Keywords: Experimental Wind-Driven Rain Wind Tunnel Testing, Heavy Rainfall, The Ahmed Body, Utility Truck, Morphing Boom Equipment, Discrete Phase Model (DPM), Automotive Aerodynamics, Computational Fluid Dynamics (CFD) 
    more » « less