Conservative integral form of the incompressible Navier-Stokes equations for a rapidly pitching airfoil

This study provides a simple moving-grid scheme which is based on a modified conservative form of the incompressible Navier-Stokes equations for flow around a moving rigid body. The modified integral form is conservative and seeks the solution of the absolute velocity. This approach is different from previous conservative differential forms [1-3] whose reference frame is not inertial. Keeping the reference frame being inertial results in simpler mathematical derivation to the governing equation which includes one dyadic product of velocity vectors in the convective term, whereas the previous [2,3] needs to obtain the time derivative with respect to non-inertial frames causing an additional dyadic product in the convective term. The scheme is implemented in a second-order accurate Navier-Stokes solver and maintains the order of the accuracy. After this verification, the scheme is validated for a pitching airfoil with very high frequencies. The simulation results match very well with the experimental results [4,5], including vorticity fields and a net thrust force. This airfoil simulation also provides detailed vortical structures near the trailing edge and time-evolving aerodynamic forces that are used to investigate the mechanism of the thrust force generation and the effects of the trailing edge shape. The developed moving-grid scheme demonstrates its validity for a rapid oscillating motion.

[1]  R. Ramamurti,et al.  Simulation of Flow About Flapping Airfoils Using Finite Element Incompressible Flow Solver , 2001 .

[2]  P. J. Zwart,et al.  The Integrated Space-Time Finite Volume Method and Its Application to Moving Boundary Problems , 2000 .

[3]  A. Hedenström A general law for animal locomotion? , 2004, Trends in ecology & evolution.

[4]  Anthony J. Calise,et al.  Flow Simulation of a Controlled Airfoil With Synthetic Jet Actuators , 2009 .

[5]  Joel H. Ferziger,et al.  Computational methods for fluid dynamics , 1996 .

[6]  A. D. Gosman,et al.  Two calculation procedures for steady, three-dimensional flows with recirculation , 1973 .

[7]  Haecheon Choi,et al.  Immersed boundary method for flow around an arbitrarily moving body , 2006, J. Comput. Phys..

[8]  D. Gottlieb,et al.  Time-stable boundary conditions for finite-difference schemes solving hyperbolic systems: methodology and application to high-order compact schemes , 1994 .

[9]  L. K. Taylor,et al.  Strong Conservative Form of the Incompressible Navier—Stokes Equations in a Rotating Frame with a Solution Procedure , 1996 .

[10]  P. Moin,et al.  Effects of the Computational Time Step on Numerical Solutions of Turbulent Flow , 1994 .

[11]  Miguel R. Visbal,et al.  High-Fidelity Simulation of Transitional Flows Past a Plunging Airfoil , 2009 .

[12]  Thomas Bewley,et al.  On the contravariant form of the Navier-Stokes equations in time-dependent curvilinear coordinate systems , 2004 .

[13]  P. Moin,et al.  A numerical method for large-eddy simulation in complex geometries , 2004 .

[14]  Keiji Kawachi,et al.  Regular Article: A Numerical Study of Undulatory Swimming , 1999 .

[15]  Dongjoo Kim,et al.  A Second-Order Time-Accurate Finite Volume Method for Unsteady Incompressible Flow on Hybrid Unstructured Grids , 2000 .

[16]  Manoochehr Koochesfahani,et al.  MTV measurements of the vortical field in the wake of an airfoil oscillating at high reduced frequency , 2009, Journal of Fluid Mechanics.

[17]  M. Koochesfahani Vortical patterns in the wake of an oscillating airfoil , 1987 .

[18]  P. J. Zwart,et al.  Regular Article: The Integrated Space-Time Finite Volume Method and Its Application to Moving Boundary Problems , 1999 .

[19]  J. Whitelaw,et al.  Flow-around airfoils with blunt, round, and sharp trailing edges , 1988 .

[20]  M. Damodaran,et al.  Computational Study of Unsteady Low-Reynolds-Number Airfoil Aerodynamics Using Moving Overlapping Meshes , 2008 .

[21]  D. Spalding,et al.  A calculation procedure for heat, mass and momentum transfer in three-dimensional parabolic flows , 1972 .

[22]  William H. Rae,et al.  Low-Speed Wind Tunnel Testing , 1966 .