The numerical continuation of the Reynolds Averaged Navier-Stokes equations involves computation of the solution of large non-linear systems, resulting in high computational cost. Three different methods are implemented to reduce the computational cost of these calculations. Firstly, the turbulence model is decoupled from the main flow equations allowing two smaller systems to be solved. The Recursive Projection Method is implemented to allow existing time integration techniques to be used in the computation of unstable solutions. Finally, bifurcation tracking methods are presented that allow the bifurcation point to be tracked with two parameters without the need for multiple continuation runs. The decoupling method works up to parameter values close to the bifurcation point, where flow separation occurs and the cross-derivative terms in the Jacobian become more important. The Recursive Projection Method shows potential in extending the applicability of time integration methods by allowing them to converge onto unstable solutions. The bifurcation trackers successfully follow the Hopf bifurcation in multiple parameters allowing loci of bifurcations to be mapped in two parameter space. Results are shown for aerofoils at a range of transonic conditions.
[1]
Cyrus K. Aidun,et al.
A direct method for computation of simple bifurcations
,
1995
.
[2]
Peter Wriggers,et al.
A general procedure for the direct computation of turning and bifurcation points
,
1990
.
[3]
Andreas Griewank,et al.
The Calculation of Hopf Points by a Direct Method
,
1983
.
[4]
Edward N. Tinoco,et al.
Abridged Summary of the Third AIAA Computational Fluid Dynamics Drag Prediction Workshop
,
2008
.
[5]
Willy Govaerts.
Stable solvers and block elimination for bordered systems
,
1991
.
[6]
David W. Zingg,et al.
Newton-Krylov Algorithm with a Loosely Coupled Turbulence Model for Aerodynamic Flows
,
2007
.
[7]
P. Spalart.
A One-Equation Turbulence Model for Aerodynamic Flows
,
1992
.
[8]
Louis A. Romero,et al.
Bifurcation Tracking Algorithms and Software for Large Scale Applications
,
2005,
Int. J. Bifurc. Chaos.
[9]
Gautam M. Shroff,et al.
Stabilization of unstable procedures: the recursive projection method
,
1993
.
[10]
Christopher L. Rumsey,et al.
Prediction of high lift: review of present CFD capability
,
2002
.