Bifurcation Tracking for High Reynolds Number Flow Around an Airfoil

High Reynolds number flows are typical for many applications including those found in aerospace. In these conditions nonlinearities arise which can, under certain conditions, result in instabilities of the flow. The accurate prediction of these instabilities is vital to enhance understanding and aid in the design process. The stability boundary can be traced by following the path of a bifurcation as two parameters are varied using a direct bifurcation tracking method. Historically, these methods have been applied to small-scale systems and only more recently have been used for large systems as found in Computational Fluid Dynamics. However, these have all been concerned with flows that are inviscid. We show how direct bifurcation tracking methods can be applied efficiently to high Reynolds number flows around an airfoil. This has been demonstrated through the presentation of a number of test cases using both flow and geometrical parameters.

[1]  Ken Badcock,et al.  Direct Aeroelastic Bifurcation Analysis of a Symmetric Wing Based on Euler Equations , 2005 .

[2]  P. Spalart A One-Equation Turbulence Model for Aerodynamic Flows , 1992 .

[3]  Richard B. Lehoucq,et al.  Large‐scale eigenvalue calculations for stability analysis of steady flows on massively parallel computers , 2001 .

[4]  Louis A. Romero,et al.  Bifurcation Tracking Algorithms and Software for Large Scale Applications , 2005, Int. J. Bifurc. Chaos.

[5]  Werner C. Rheinboldt,et al.  Algorithm 596: a program for a locally parameterized , 1983, TOMS.

[6]  Willy Govaerts,et al.  MATCONT: A MATLAB package for numerical bifurcation analysis of ODEs , 2003, TOMS.

[7]  K. Badcock,et al.  Fast prediction of transonic aeroelastic stability and limit cycles , 2007 .

[8]  Thomas F. Fairgrieve,et al.  AUTO 2000 : CONTINUATION AND BIFURCATION SOFTWARE FOR ORDINARY DIFFERENTIAL EQUATIONS (with HomCont) , 1997 .

[9]  Christopher L. Rumsey,et al.  Prediction of high lift: review of present CFD capability , 2002 .

[10]  R. Lehoucq,et al.  Computational bifurcation and stability studies of the 8: 1 thermal cavity problem , 2002 .

[11]  Heinrich Voss,et al.  A new justification of the Jacobi–Davidson method for large eigenproblems , 2007 .

[12]  Daniella E. Raveh,et al.  Reynolds-Averaged Navier-Stokes Study of the Shock-Buffet Instability Mechanism , 2012 .

[13]  Timothy A. Davis,et al.  Algorithm 832: UMFPACK V4.3---an unsymmetric-pattern multifrontal method , 2004, TOMS.

[14]  Dirk Roose,et al.  An Adaptive Newton-Picard Algorithm with Subspace Iteration for Computing Periodic Solutions , 1998, SIAM J. Sci. Comput..

[15]  Willy Govaerts,et al.  Stable solvers and block elimination for bordered systems , 1991 .

[16]  Timothy A. Davis,et al.  A combined unifrontal/multifrontal method for unsymmetric sparse matrices , 1999, TOMS.

[17]  Edward N. Tinoco,et al.  Abridged Summary of the Third AIAA Computational Fluid Dynamics Drag Prediction Workshop , 2008 .

[18]  Fred Wubs,et al.  Numerical Bifurcation Methods and their Application to Fluid Dynamics: Analysis beyond Simulation , 2012 .

[19]  Ken Badcock,et al.  Hopf Bifurcation Calculations for a Symmetric Airfoil in Transonic Flow , 2004 .

[20]  Chao Yang,et al.  ARPACK users' guide - solution of large-scale eigenvalue problems with implicitly restarted Arnoldi methods , 1998, Software, environments, tools.

[21]  Willy Govaerts,et al.  Numerical computation of bifurcations in large equilibrium systems in matlab , 2014, J. Comput. Appl. Math..

[22]  Peter Wriggers,et al.  A general procedure for the direct computation of turning and bifurcation points , 1990 .

[23]  Andreas Griewank,et al.  The Calculation of Hopf Points by a Direct Method , 1983 .

[24]  Alan R. Champneys,et al.  Numerical continuation of high Reynolds number external flows , 2012 .

[25]  C. WALES,et al.  An Initial Study of the Flow around an aerofoil at High Reynolds numbers using Continuation , 2012, Int. J. Bifurc. Chaos.

[26]  A. Spence,et al.  The numerical analysis of bifurcation problems with application to fluid mechanics , 2000, Acta Numerica.

[27]  Timothy A. Davis,et al.  A column pre-ordering strategy for the unsymmetric-pattern multifrontal method , 2004, TOMS.

[28]  Timothy A. Davis,et al.  An Unsymmetric-pattern Multifrontal Method for Sparse Lu Factorization , 1993 .

[29]  Cyrus K. Aidun,et al.  A direct method for computation of simple bifurcations , 1995 .