Simulation of flapping wings subjected to gusty inflow

2019 ◽  
Vol 123 (1266) ◽  
pp. 1170-1192 ◽  
Author(s):  
M. M. De ◽  
J. S. Mathur ◽  
S. Vengadesan

ABSTRACTOrnithopters and entomopters should be insensitive to the gusty environment during outdoor operations. Hence, it becomes imperative to understand their behaviour under the influence of gust for ensuring stable flight. In light of this, the present numerical study focused on understanding the aerodynamics of flapping wings with five different planform shapes under the influence of a spatiotemporally varying frontal gust. 3D, unsteady, laminar, and incompressible Navier-Stokes equations were solved using finite volume formulation. A canonical case of asymmetric 1 degree of freedom (DoF) flapping kinematics was considered. Horizontal and vertical force patterns in constant and gusty inflow conditions were numerically computed and compared. Findings were analyzed quantitatively by comparing the differences in the instantaneous force patterns, ordinal scoring approach, and phase space plots. Qualitative comparisons were made based on plots of vortex structures and surface pressure contours for constant and gusty inflow conditions for wings with different planform shapes. Spanwise Lagrangian Coherent Structures (LCS) of all the five wings were also compared. Studies revealed that the elliptical wing exhibited low sensitivity and inverse semi-elliptical wing exhibited high sensitivity to the gusty inflow. Rectangular, triangular and semi-elliptical shaped wings were moderately sensitive to the gusty inflow. This finding, within the limitations of the flapping kinematics and simulation conditions considered for the present study, supported the fact that many natural flyers like forest raptors, non-migratory passerines, pheasants, and partridges have adopted elliptical wing planform for efficient flight.

Processes ◽  
2021 ◽  
Vol 9 (5) ◽  
pp. 841
Author(s):  
Yuzhen Jin ◽  
Huang Zhou ◽  
Linhang Zhu ◽  
Zeqing Li

A three-dimensional numerical study of a single droplet splashing vertically on a liquid film is presented. The numerical method is based on the finite volume method (FVM) of Navier–Stokes equations coupled with the volume of fluid (VOF) method, and the adaptive local mesh refinement technology is adopted. It enables the liquid–gas interface to be tracked more accurately, and to be less computationally expensive. The relationship between the diameter of the free rim, the height of the crown with different numbers of collision Weber, and the thickness of the liquid film is explored. The results indicate that the crown height increases as the Weber number increases, and the diameter of the crown rim is inversely proportional to the collision Weber number. It can also be concluded that the dimensionless height of the crown decreases with the increase in the thickness of the dimensionless liquid film, which has little effect on the diameter of the crown rim during its growth.


1989 ◽  
Vol 209 ◽  
pp. 285-308 ◽  
Author(s):  
R. J. Bodonyi ◽  
W. J. C. Welch ◽  
P. W. Duck ◽  
M. Tadjfar

A numerical study of the generation of Tollmien-Schlichting (T–S) waves due to the interaction between a small free-stream disturbance and a small localized variation of the surface geometry has been carried out using both finite–difference and spectral methods. The nonlinear steady flow is of the viscous–inviscid interactive type while the unsteady disturbed flow is assumed to be governed by the Navier–Stokes equations linearized about this flow. Numerical solutions illustrate the growth or decay of the T–S waves generated by the interaction between the free-stream disturbance and the surface distortion, depending on the value of the scaled Strouhal number. An important result of this receptivity problem is the numerical determination of the amplitude of the T–S waves.


1996 ◽  
Vol 118 (2) ◽  
pp. 335-340 ◽  
Author(s):  
Mohamed Selmi

This paper is concerned with the solution of the 3-D-Navier-Stokes equations describing the steady motion of a viscous fluid inside a partially filled spinning and coning cylinder. The cylinder contains either a single fluid of volume less than that of the cylinder or a central rod and a single fluid of combined volume (volume of the rod plus volume of the fluid) equal to that of the cylinder. The cylinder rotates about its axis at the spin rate ω and rotates about an axis that passes through its center of mass at the coning rate Ω. In practical applications, as in the analysis and design of liquid-filled projectiles, the parameter ε = τ sin θ, where τ = Ω/ω and θ is the angle between spin axis and coning axis, is small. As a result, linearization of the Navier-Stokes equations with this parameter is possible. Here, the full and linearized Navier-Stokes equations are solved by a spectral collocation method to investigate the nonlinear effects on the moments caused by the motion of the fluid inside the cylinder. In this regard, it has been found that nonlinear effects are negligible for τ ≈ 0.1, which is of practical interest to the design of liquid-filled projectiles, and the solution of the linearized Navier-Stokes equations is adequate for such a case. However, as τ increases, nonlinear effects increase, and become significant as ε surpasses about 0.1. In such a case, the nonlinear problem must be solved. Complete details on how to solve such a problem is presented.


Author(s):  
David Gross ◽  
Yann Roux ◽  
Benjamin Rousse ◽  
François Pétrié ◽  
Ludovic Assier ◽  
...  

The problem of Vortex-Induced Vibrations (VIV) on spool and jumper geometries is known to present several drawbacks when approached with conventional engineering tools used in the study of VIV on risers. Current recommended practices can lead to over-conservatism that the industry needs to quantify and minimize within notably cost reduction objectives. Within this purpose, the paper will present a brief critical review of the Industry standards and more particularly focus on both experimental and Computational Fluid Dynamic (CFD) approaches. Both qualitative and quantitative comparisons between basin tests and CFD results for a 2D ‘M-shape’ spool model will be detailed. The results presented here are part of a larger experimental and numerical campaign which considered a number of current velocities, heading and geometry configurations. The vibratory response of the model will be investigated for one of the current velocities and compared with the results obtained through recommended practices (e.g. Shear7 and DNV guidelines). The strategy used by the software K-FSI to solve the fluid-structure interaction (FSI) problem is a partitioned coupling solver between fluid solver (FINE™/Marine) and structural solvers (ARA). FINE™/Marine solves the Reynolds-Averaged Navier-Stokes Equations in a conservative way via the finite volume method and can work on structured or unstructured meshes with arbitrary polyhedrons, while ARA is a nonlinear finite element solver with a large displacement formulation. The experiments were conducted in the BGO FIRST facility located in La Seyne sur Mer, France. Particular attention was paid towards the model design, fabrication, instrumentation and characterization, to ensure an excellent agreement between the structural numerical model and the actual physical model. This included the use of a material with low structural damping, the performance of stiffness and decay tests in air and in still water, plus the rationalization of the instrumentation to be able to capture the response with the minimum flow perturbation or interaction due to instrumentation.


1998 ◽  
Vol 14 (1) ◽  
pp. 23-29
Author(s):  
Robert R. Hwang ◽  
Sheng-Yuh Jaw

ABSTRACTThis paper presents a numerical study on turbulent vortex shedding flows past a square cylinder. The 2D unsteady periodic shedding motion was resolved in the calculation and the superimposed turbulent fluctuations were simulated with a second-order Reynolds-stress closure model. The calculations were carried out by solving numerically the fully elliptic ensemble-averaged Navier-Stokes equations coupled with the turbulence model equations together with the two-layer approach in the treatment of the near-wall region. The performance of the computations was evaluated by comparing the numerical results with data from available experiments. Results indicate that the present study gives good agreement in the shedding frequency and mean drag as well as in some phase profiles of the mean velocity.


2014 ◽  
Vol 752 ◽  
pp. 602-625 ◽  
Author(s):  
Kengo Deguchi ◽  
Philip Hall

AbstractOur concern in this paper is with high-Reynolds-number nonlinear equilibrium solutions of the Navier–Stokes equations for boundary-layer flows. Here we consider the asymptotic suction boundary layer (ASBL) which we take as a prototype parallel boundary layer. Solutions of the equations of motion are obtained using a homotopy continuation from two known types of solutions for plane Couette flow. At high Reynolds numbers, it is shown that the first type of solution takes the form of a vortex–wave interaction (VWI) state, see Hall & Smith (J. Fluid Mech., vol. 227, 1991, pp. 641–666), and is located in the main part of the boundary layer. On the other hand, here the second type is found to support an equilibrium solution of the unit-Reynolds-number Navier–Stokes equations in a layer located a distance of $\def \xmlpi #1{}\def \mathsfbi #1{\boldsymbol {\mathsf {#1}}}\let \le =\leqslant \let \leq =\leqslant \let \ge =\geqslant \let \geq =\geqslant \def \Pr {\mathit {Pr}}\def \Fr {\mathit {Fr}}\def \Rey {\mathit {Re}}O(\ln \mathit{Re})$ from the wall. Here $\mathit{Re}$ is the Reynolds number based on the free-stream speed and the unperturbed boundary-layer thickness. The streaky field produced by the interaction grows exponentially below the layer and takes its maximum size within the unperturbed boundary layer. The results suggest the possibility of two distinct types of streaky coherent structures existing, possibly simultaneously, in disturbed boundary layers.


2013 ◽  
Vol 397-400 ◽  
pp. 783-788
Author(s):  
Xing Wei Zhang ◽  
Chao Wang ◽  
Hang Liu

This paper investigates the aerodynamic forces of several plunging wing models by means of computational fluid dynamics. A finite volume method was used to solve the two-dimensional unsteady incompressible Navier-Stokes equations. The forces and power efficiency have been calculated and compared between sets of different models. Current work found that the nonsymmetrical moving can enhance the lift and thrust forces. The numerical results also prove that the flexible wing model can be use to improve the efficiency and reduce the input. Additionally, a new conceptual model for flapping wing mechanism with active deformation and adaptive nonsymmetrical driving motion is proposed base on the numerical results.


Author(s):  
Fayçal Hammami ◽  
Nader Ben Cheikh ◽  
Brahim Ben Beya

This paper deals with the numerical study of bifurcations in a two-sided lid driven cavity flow. The flow is generated by moving the upper wall to the right while moving the left wall downwards. Numerical simulations are performed by solving the unsteady two dimensional Navier-Stokes equations using the finite volume method and multigrid acceleration. In this problem, the ratio of the height to the width of the cavity are ranged from H/L = 0.25 to 1.5. The code for this cavity is presented using rectangular cavity with the grids 144 × 36, 144 × 72, 144 × 104, 144 × 136, 144 × 176 and 144 × 216. Numerous comparisons with the results available in the literature are given. Very good agreements are found between current numerical results and published numerical results. Various velocity ratios ranged in 0.01≤ α ≤ 0.99 at a fixed aspect ratios (A = 0.5, 0.75, 1.25 and 1.5) were considered. It is observed that the transition to the unsteady regime follows the classical scheme of a Hopf bifurcation. The stability analysis depending on the aspect ratio, velocity ratios α and the Reynolds number when transition phenomenon occurs is considered in this paper.


Sign in / Sign up

Export Citation Format

Share Document