High order accurate and low dissipation method for unsteady compressible viscous flow computation on helicopter rotor in forward flight

An improved fifth-order weighted essentially non-oscillatory (WENO-Z) scheme combined with the moving overset grid technique has been developed to compute unsteady compressible viscous flows on the helicopter rotor in forward flight. In order to enforce periodic rotation and pitching of the rotor and relative motion between rotor blades, the moving overset grid technique is extended, where a special judgement standard is presented near the odd surface of the blade grid during search donor cells by using the Inverse Map method. The WENO-Z scheme is adopted for reconstructing left and right state values with the Roe Riemann solver updating the inviscid fluxes and compared with the monotone upwind scheme for scalar conservation laws (MUSCL) and the classical WENO scheme. Since the WENO schemes require a six point stencil to build the fifth-order flux, the method of three layers of fringes for hole boundaries and artificial external boundaries is proposed to carry out flow information exchange between chimera grids. The time advance on the unsteady solution is performed by the full implicit dual time stepping method with Newton type LU-SGS subiteration, where the solutions of pseudo steady computation are as the initial fields of the unsteady flow computation. Numerical results on non-variable pitch rotor and periodic variable pitch rotor in forward flight reveal that the approach can effectively capture vortex wake with low dissipation and reach periodic solutions very soon.

[1]  Wai-Sun Don,et al.  An improved weighted essentially non-oscillatory scheme for hyperbolic conservation laws , 2008, J. Comput. Phys..

[2]  Lakshmi N. Sankar,et al.  Hybrid Navier-Stokes/Free-Wake Method for Modeling Blade-Vortex Interactions , 2010 .

[3]  Brian E. Wake,et al.  Recent rotor wake simulation and modeling studies at United Technologies Corporation , 2000 .

[4]  H. Lomax,et al.  Thin-layer approximation and algebraic model for separated turbulent flows , 1978 .

[5]  B. V. Leer,et al.  Towards the Ultimate Conservative Difference Scheme , 1997 .

[6]  A. Jameson,et al.  Lower-upper Symmetric-Gauss-Seidel method for the Euler and Navier-Stokes equations , 1988 .

[7]  A. Jameson Time dependent calculations using multigrid, with applications to unsteady flows past airfoils and wings , 1991 .

[8]  Robert L. Meakin,et al.  On automating domain connectivity for overset grids , 1995 .

[9]  Jasim Ahmad,et al.  Helicopter rotor blade computation in unsteady flows using moving overset grids , 1996 .

[10]  D. Whitfield,et al.  Numerical simulation of a helicopter rotor in hover and forward flight , 1995 .

[11]  Juan J. Alonso,et al.  Time-accurate simulation of helicopter rotor flows including aeroelastic effects , 1997 .

[12]  W. J. McCroskey Vortex wakes of rotorcraft , 1995 .

[13]  Robert E. Harris,et al.  An Efficient Adaptive Cartesian Vorticity Transport Solver for Rotorcraft Flowfield Analysis , 2010 .

[14]  A. T. Conlisk,et al.  Modern helicopter rotor aerodynamics , 2001 .

[15]  Robert Meakin A new method for establishing intergrid communication among systems of overset grids , 1991 .

[16]  Siegfried Wagner,et al.  Tip Vortex Conservation on a Helicopter Main Rotor Using Vortex-Adapted Chimera Grids , 2007 .

[17]  Chengjian He,et al.  Modeling Rotor Wake Dynamics with Viscous Vortex Particle Method , 2008 .

[18]  Berend G. van der Wall,et al.  Chimera simulations of multibladed rotors in high-speed forward flight with weak fluid-structure-coupling , 2005 .

[19]  Implicit multigrid algorithm for the Navier-Stokes equations , 1991 .

[20]  F. X. Caradonna,et al.  Experimental and Analytical Studies of a Model Helicopter Rotor in Hover , 1980 .

[21]  S. Osher,et al.  Weighted essentially non-oscillatory schemes , 1994 .

[22]  Wai-Sun Don,et al.  High order weighted essentially non-oscillatory WENO-Z schemes for hyperbolic conservation laws , 2011, J. Comput. Phys..

[23]  J. M. Powers,et al.  Mapped weighted essentially non-oscillatory schemes: Achieving optimal order near critical points , 2005 .

[24]  Soohyung Park,et al.  Euler and Navier-Stokes Simulations of Helicopter Rotor Blade in Forward Flight Using an Overlapped Grid Solver , 2009 .

[25]  Wai-Sun Don,et al.  Accuracy of the weighted essentially non-oscillatory conservative finite difference schemes , 2013, J. Comput. Phys..

[26]  V. P. Kolgan,et al.  Application of the principle of minimizing the derivative to the construction of finite-difference schemes for computing discontinuous solutions of gas dynamics , 2011, J. Comput. Phys..

[27]  C. Angelopoulos High resolution schemes for hyperbolic conservation laws , 1992 .

[28]  Rohit Jain,et al.  Capability Enhancements in Version 3 of the Helios High-Fidelity Rotorcraft Simulation Code , 2012 .

[29]  Li Xu,et al.  Rotor wake capture improvement based on high-order spatially accurate schemes and chimera grids , 2011 .

[30]  Chi-Wang Shu,et al.  Efficient Implementation of Weighted ENO Schemes , 1995 .

[31]  Prabhu Ramachandran,et al.  Approximate Riemann solvers for the Godunov SPH (GSPH) , 2014, J. Comput. Phys..

[32]  Lakshmi N. Sankar,et al.  Application of Higher Order Spatially Accurate Schemes to Rotors in Hover , 2006 .

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

[34]  Lakshmi N. Sankar,et al.  High-Order Essentially Nonoscillatory Schemes for Rotary-Wing Wake Computations , 2004 .

[35]  Eleuterio F. Toro,et al.  WENO schemes based on upwind and centred TVD fluxes , 2005 .