Numerical simulation of flapping wings using a panel method and a high‐order Navier–Stokes solver

SUMMARY The design of efficient flapping wings for human engineered micro aerial vehicles (MAVs) has long been an elusive goal, in part because of the large size of the design space. One strategy for overcoming this difficulty is to use a multifidelity simulation strategy that appropriately balances computation time and accuracy. We compare two models with different geometric and physical fidelity. The low-fidelity model is an inviscid doublet lattice method with infinitely thin lifting surfaces. The high-fidelity model is a high-order accurate discontinuous Galerkin Navier–Stokes solver, which uses an accurate representation of the flapping wing geometry. To compare the performance of the two methods, we consider a model flapping wing with an elliptical planform and an analytically prescribed spanwise wing twist, at size scales relevant to MAVs. Our results show that in many cases, including those with mild separation, low-fidelity simulations can accurately predict integrated forces, provide insight into the flow structure, indicate regions of likely separation, and shed light on design–relevant quantities. But for problems with significant levels of separation, higher-fidelity methods are required to capture the details of the flow field. Inevitably high-fidelity simulations are needed to establish the limits of validity of the lower fidelity simulations.Copyright © 2012 John Wiley & Sons, Ltd.

[1]  Ahmed F. Ghoniem,et al.  Three-Dimensional Vortex Simulation of Time Dependent Incompressible Internal Viscous Flows , 1997 .

[2]  Jacob K. White,et al.  A combined pFFT‐multipole tree code, unsteady panel method with vortex particle wakes , 2007 .

[3]  Per-Olof Persson,et al.  A High Order Discontinuous Galerkin Method for Fluid-Structure Interaction , 2007 .

[4]  A. Leonard Computing Three-Dimensional Incompressible Flows with Vortex Elements , 1985 .

[5]  J. Rayner A vortex theory of animal flight. Part 1. The vortex wake of a hovering animal , 1979, Journal of Fluid Mechanics.

[6]  Per-Olof Persson,et al.  Multifidelity Approaches for the Computational Analysis and Design of Eective Flapping Wing Vehicles , 2008 .

[7]  Max F. Platzer,et al.  Thrust Generation due to Airfoil Flapping , 1996 .

[8]  R. Alexander Diagonally implicit runge-kutta methods for stiff odes , 1977 .

[9]  Y. Saad,et al.  GMRES: a generalized minimal residual algorithm for solving nonsymmetric linear systems , 1986 .

[10]  A. U.S.,et al.  Curved Mesh Generation and Mesh Refinement using Lagrangian Solid Mechanics , 2009 .

[11]  Joseba Murua,et al.  Structural and Aerodynamic Models in Nonlinear Flight Dynamics of Very Flexible Aircraft , 2010 .

[12]  J. D. Delaurier,et al.  An aerodynamic model for flapping-wing flight , 1993, The Aeronautical Journal (1968).

[13]  Sheila E. Widnall,et al.  Energetics and optimum motion of oscillating lifting surfaces of finite span , 1986, Journal of Fluid Mechanics.

[14]  Thomas J. Mueller,et al.  Euler Solutions for a Finite-Span Flapping Wing , 2001 .

[15]  Jacob K. White,et al.  A precorrected-FFT method for electrostatic analysis of complicated 3-D structures , 1997, IEEE Trans. Comput. Aided Des. Integr. Circuits Syst..

[16]  Per-Olof Persson,et al.  Implicit Large Eddy Simulation of transition to turbulence at low Reynolds numbers using a Discontinuous Galerkin method , 2011 .

[17]  Jaime Peraire,et al.  Discontinuous Galerkin Solution of the Navier-Stokes Equations on Deformable Domains , 2007 .

[18]  Gareth Jones,et al.  Bird Flight Performance. A Practical Calculation Manual, C.J. Pennycuick. Oxford University Press, Oxford (1989), x, +153. Price £25 , 1990 .

[19]  L. Greengard The Rapid Evaluation of Potential Fields in Particle Systems , 1988 .

[20]  Chi-Wang Shu,et al.  Runge–Kutta Discontinuous Galerkin Methods for Convection-Dominated Problems , 2001, J. Sci. Comput..

[21]  N. J.,et al.  Distributions of sources and normal dipoles over a quadrilateral panel , 2004 .

[22]  Thomas J. Mueller,et al.  Fixed and Flapping Wing Aerodynamics for Micro Air Vehicle Applications , 2001 .

[23]  Per-Olof Persson,et al.  RANS Solutions Using High Order Discontinuous Galerkin Methods , 2007 .

[24]  Hao Liu,et al.  Flapping Wings and Aerodynamic Lift: The Role of Leading-Edge Vortices , 2007 .

[25]  Bret Stanford,et al.  Analytical Sensitivity Analysis of an Unsteady Vortex Lattice Method for Flapping Wing Optimization , 2009 .

[26]  James D. DeLaurier,et al.  A Nonlinear Aeroelastic Model for the Study of Flapping Wing Flight , 2008 .

[27]  Kenneth C. Hall,et al.  Power Requirements for Large-Amplitude Flapping Flight , 1997 .

[28]  D. Crighton The Kutta Condition in Unsteady Flow , 1985 .

[29]  P. Roe Approximate Riemann Solvers, Parameter Vectors, and Difference Schemes , 1997 .

[30]  Per-Olof Persson,et al.  Newton-GMRES Preconditioning for Discontinuous Galerkin Discretizations of the Navier--Stokes Equations , 2008, SIAM J. Sci. Comput..

[31]  J. N. Newman Distributions of sources and normal dipoles over a quadrilateral panel , 1986 .

[32]  K. Morganyand,et al.  Unstructured Grid Finite Element Methods for Fluid Mechanics , 1997 .

[33]  C. E. Lan,et al.  The unsteady quasi-vortex-lattice method with applications to animal propulsion , 1979, Journal of Fluid Mechanics.

[34]  W. Rodden,et al.  A doublet-lattice method for calculating lift distributions on oscillating surfaces in subsonic flows. , 1969 .

[35]  L. Morino,et al.  Subsonic Potential Aerodynamics for Complex Configurations: A General Theory , 1974 .

[36]  Joseph Katz,et al.  Aerodynamic study of a flapping-wing micro-UAV , 1999 .

[37]  Z. Jane Wang,et al.  DISSECTING INSECT FLIGHT , 2005 .

[38]  Per-Olof Persson,et al.  The Compact Discontinuous Galerkin (CDG) Method for Elliptic Problems , 2007, SIAM J. Sci. Comput..

[39]  J. Katz,et al.  Low-Speed Aerodynamics , 1991 .

[40]  Per-Olof Persson,et al.  A Computational Framework for Fluid Structure Interaction in Biologically Inspired Flapping Flight , 2007 .

[41]  Per-Olof Persson,et al.  A Simple Mesh Generator in MATLAB , 2004, SIAM Rev..

[42]  J. Hess,et al.  The problem of three-dimensional lifting potential flow and its solution by means of surface singularity distribution☆ , 1974 .

[43]  Andrew W. Appel,et al.  An Efficient Program for Many-Body Simulation , 1983 .

[44]  Per-Olof Persson,et al.  Curved mesh generation and mesh refinement using Lagrangian solid mechanics , 2008 .

[45]  Lyle N. Long,et al.  Object-oriented unsteady vortex lattice method for flapping flight , 2004 .

[46]  Smith,et al.  The advantages of an unsteady panel method in modelling the aerodynamic forces on rigid flapping wings , 1996, The Journal of experimental biology.

[47]  Per-Olof Persson,et al.  An Efficient Low Memory Implicit DG Algorithm for Time Dependent Problems , 2006 .

[48]  Per-Olof Persson,et al.  Scalable parallel Newton-Krylov solvers for discontinuous Galerkin discretizations , 2009 .