Krylov methods for the incompressible Navier-Stokes equations
Journal of Computational Physics
A high-resolution hybrid compact-ENO scheme for shock-turbulence interaction problems
Journal of Computational Physics
Low-storage, explicit Runge-Kutta schemes for the compressible Navier-Stokes equations
Applied Numerical Mathematics
A Krylov--Schur Algorithm for Large Eigenproblems
SIAM Journal on Matrix Analysis and Applications
Iterative Methods for Sparse Linear Systems
Iterative Methods for Sparse Linear Systems
Jacobian-free Newton-Krylov methods: a survey of approaches and applications
Journal of Computational Physics
SLEPc: A scalable and flexible toolkit for the solution of eigenvalue problems
ACM Transactions on Mathematical Software (TOMS) - Special issue on the Advanced CompuTational Software (ACTS) Collection
Predicting the onset of flow unsteadiness based on global instability
Journal of Computational Physics
Journal of Computational Physics
Journal of Computational Physics
Hi-index | 31.45 |
The direct and adjoint operators play an undeniably important role in a vast number of theoretical and practical studies that range from linear stability to flow control and nonlinear optimization. Based on an existing nonlinear flow solver, the design of efficient and straightforward procedures to access these operators is thus highly desirable. In the case of compressible solvers, the use of high-order numerical schemes combined with complicated governing equations makes the derivation of efficient procedures a challenging and often tedious undertaking. In this work, a novel technique for the evaluation of the direct and adjoint operators directly from compressible flow solvers is presented and extended to include nonlinear differentiation schemes and turbulence models. The application to the incompressible counterpart is also discussed. The presented method requires minimal additional programming effort and automatically takes into account subsequent modifications in the governing equations and boundary conditions. The introduced methodology is demonstrated on existing numerical codes, and direct and adjoint global modes are calculated for three typical flow configurations. Implementation issues and the performance measures are also discussed. The proposed algorithm presents an easy-to-implement and efficient technique to extract valuable information for the quantitative analysis of complex flows.