Solution of the Unsteady Euler Equations Using an Implicit Dual-Time Method

An unfactored implicit time-marching method for the solution of the unsteady two-dimensional Euler equations on deforming grids is described. The present work is placed into a multiblock framework and e ts into the development of a generally applicable parallel multiblock e ow solver. The convective terms are discretized using an upwind total variation diminishing scheme, whereas the unsteady governing equations are discretized using an implicit dual-time approach. The large sparse linear system arising from the implicit time discretization at each pseudotime step is solved efe ciently by using a conjugate-gradient-type method with a preconditioning based on a block incomplete lower-upper factorization. Results are shown for a series of pitching airfoil test cases selected from the AGARD aeroelastic cone gurations for the NACA 0012 airfoil. Comparisons with experimental data and previous published results are presented. The efe ciency of the method is demonstrated by looking at the effect of a number of numerical parameters, such as the conjugate gradient tolerance and the size of the global time step and by carrying out a grid ree nement study. Finally, a demonstration test case forthe Williamsairfoil (Williams, B. R., “ An Exact Test Case for the Plane Potential Flow About Two Adjacent Lifting Aerofoils,” National Physical Lab., Aeronautical Research Council, Research Memorandum 3717, London, 1973 )with an oscillating e ap is presented, highlighting the capability of the grid deformation technique.