zbMATH — the first resource for mathematics

Enabling off-design linearised aerodynamics analysis using Krylov subspace recycling technique. (English) Zbl 1390.76267
Summary: The major computational challenge, when using frequency domain linearised computational fluid dynamics in the analysis of aeroelastic problems such as aircraft flutter, gust response or shock buffet, are the excessive memory and CPU time requirements to solve the large sparse linear systems of equations. To address these issues found with the generalised minimal residual linear equation solver, the generalised conjugate residual solver with deflated restarting is adopted here in which an invariant Krylov subspace is recycled both between restarts when solving a single linear frequency domain problem and for a sequence of equations when varying the system matrix and forcing terms. The proposed method is implemented in an industrial code and applied to three test cases including the forced excitation and buffet onset of a pitch-plunge aerofoil, a realistic passenger aircraft in inviscid transonic flow and a generic half wing-body model at a pre-buffet condition. The memory requirements for the problems investigated are reduced by up to an order of magnitude, while the CPU times are reduced by up to a factor of three.

76G25 General aerodynamics and subsonic flows
65F10 Iterative numerical methods for linear systems
74F10 Fluid-solid interactions (including aero- and hydro-elasticity, porosity, etc.)
76-04 Software, source code, etc. for problems pertaining to fluid mechanics
76H05 Transonic flows
Full Text: DOI
[1] Thormann, R.; Widhalm, M., Linear-frequency-domain predictions of dynamic-response data for viscous transonic flows, AIAA J, 51, 11, 2540-2557, (2013)
[2] Timme, S.; Marques, S.; Badcock, K., Transonic aeroelastic stability analysis using a Kriging-based Schur complement formulation, AIAA J, 49, 6, 1202-1213, (2011)
[3] Da Ronch, A.; Badcock, K. J.; Wang, Y.; Wynn, A.; Palacios, R. N., Nonlinear model reduction for flexible aircraft control design, AIAA Paper 2012-4404, (2012)
[4] Timme, S.; Badcock, K.; Da Ronch, A., Linear reduced order modelling for gust response analysis using the DLR-TAU code, Proc IFASD, (2013)
[5] Crouch, J.; Garbaruk, A.; Magidov, D., Predicting the onset of flow unsteadiness based on global instability, J Comput Phys, 224, 2, 924-940, (2007) · Zbl 1123.76018
[6] Iorio, M.; González, L.; Ferrer, E., Direct and adjoint global stability analysis of turbulent transonic flows over a NACA0012 profile, Int J Numer Methods Fluids, 76, 3, 147-168, (2014)
[7] Sartor, F.; Mettot, C.; Sipp, D., Stability, receptivity, and sensitivity analyses of buffeting transonic flow over a profile, AIAA J, 50, 4, 880-890, (2014)
[8] Nitzsche J, A numerical study on aerodynamic resonance in transonic separated flow. Proc IFASD-2009-126.
[9] Iovnovich, M.; Raveh, D. E., Reynolds-averaged Navier-Stokes study of the shock-buffet instability mechanism, AIAA J, 50, 4, 880-890, (2012)
[10] Thormann, R.; Nitzsche, J.; Widhalm, M., Time-linearized simulation of unsteady transonic flows with shock-induced separation, 6th European Congress on Computational Methods in Applied Sciences and Engineering, (2012)
[11] Cai, X.-C.; Gropp, W. D.; Keyes, D. E.; Tidriri, M. D., Newton-Krylov-Schwarz methods in CFD, Numerical methods for the Navier-Stokes equations, 17-30, (1994), Springer · Zbl 0876.76059
[12] Gropp, W.; Keyes, D.; Mcinnes, L. C.; Tidriri, M., Globalized Newton-Krylov-Schwarz algorithms and software for parallel implicit CFD, Int J High Perform Comput Appl, 14, 2, 102-136, (2000)
[13] Osusky, M.; Zingg, D. W., Parallel Newton-Krylov-Schur flow solver for the Navier-Stokes equations, AIAA J, 51, 12, 2833-2851, (2013)
[14] Xu, S.; Radford, D.; Meyer, M.; Müller, J.-D., Stabilisation of discrete steady adjoint solvers, J Comput Phys, 299, 175-195, (2015) · Zbl 1351.76131
[15] Dwight, R., Efficiency improvements of RANS-based analysis and optimization using implicit and adjoint methods on unstructured grids, (2006), University of Manchester, Ph.D. thesis
[16] Rüde, U.; Schmid, W., Inverse multigrid correction for generalized eigenvalue computations, (1995), Technical Report, Universität Augsburg
[17] Saad, Y.; Schultz, M. H., GMRES: A generalized minimal residual algorithm for solving nonsymmetric linear systems, SIAM J Sci Stat Comput, 7, 3, 856-869, (1986) · Zbl 0599.65018
[18] Eisenstat, S.; Elman, H.; Schultz, M., Variational iterative methods for nonsymmetric systems of linear equations, SIAM J Numer Anal, 20, 2, 345-357, (1983) · Zbl 0524.65019
[19] Parks, M. L.; de Sturler, E.; Mackey, G.; Johnson, D. D.; Maiti, S., Recycling Krylov subspaces for sequences of linear systems, SIAM J Sci Comput, 28, 5, 1651-1674, (2006) · Zbl 1123.65022
[20] de Sturler, E., Nested Krylov methods based on GCR, J Comput Appl Math, 67, 1, 15-41, (1996) · Zbl 0854.65026
[21] de Sturler, E., Truncation strategies for optimal Krylov subspace methods, SIAM J Numer Anal, 36, 3, 864-889, (1999) · Zbl 0960.65031
[22] Morgan, R. B., GMRES with deflated restarting, SIAM J Sci Comput, 24, 1, 20-37, (2002) · Zbl 1018.65042
[23] Kilmer, M. E.; de Sturler, E., Recycling subspace information for diffuse optical tomography, SIAM J Sci Comput, 27, 6, 2140-2166, (2006) · Zbl 1103.65036
[24] Liesen, J.; Tichỳ, P., Convergence analysis of Krylov subspace methods, GAMM-Mitteilungen, 27, 2, 153-173, (2004) · Zbl 1071.65041
[25] Titley-Peloquin, D.; Pestana, J.; Wathen, A. J., GMRES convergence bounds that depend on the right-hand-side vector, IMA J Numer Anal, 34, 2, 462-479, (2014) · Zbl 1302.65083
[26] McCracken, A. J.; Timme, S.; Badcock, K. J., Accelerating convergence of the CFD linear frequency domain method by a preconditioned linear solver, 6th European Congress on Computational Methods in Applied Sciences and Engineering, (2012)
[27] Thormann, R.; Widhalm, M., Forced motion simulations using a linear frequency domain solver for a generic transport aircraft, Proc IFASD, (2013)
This reference list is based on information provided by the publisher or from digital mathematics libraries. Its items are heuristically matched to zbMATH identifiers and may contain data conversion errors. It attempts to reflect the references listed in the original paper as accurately as possible without claiming the completeness or perfect precision of the matching.