Development of iterative techniques for the solution of unsteady compressible viscous flows

The development of efficient iterative solution methods for the numerical solution of two- and three-dimensional compressible Navier-Stokes equations is discussed. Iterative time marching methods have several advantages over classical multi-step explicit time marching schemes, and non-iterative implicit time marching schemes. Iterative schemes have better stability characteristics than non-iterative explicit and implicit schemes. In this work, another approach based on the classical conjugate gradient method, known as the Generalized Minimum Residual (GMRES) algorithm is investigated. The GMRES algorithm has been used in the past by a number of researchers for solving steady viscous and inviscid flow problems. Here, we investigate the suitability of this algorithm for solving the system of non-linear equations that arise in unsteady Navier-Stokes solvers at each time step

Development of iterative techniques for the solution of unsteady compressible viscous flows

Efficient iterative solution methods are being developed for the numerical solution of two- and three-dimensional compressible Navier-Stokes equations. Iterative time marching methods have several advantages over classical multi-step explicit time marching schemes, and non-iterative implicit time marching schemes. Iterative schemes have better stability characteristics than non-iterative explicit and implicit schemes. Thus, the extra work required by iterative schemes can also be designed to perform efficiently on current and future generation scalable, missively parallel machines. An obvious candidate for iteratively solving the system of coupled nonlinear algebraic equations arising in CFD applications is the Newton method. Newton's method was implemented in existing finite difference and finite volume methods. Depending on the complexity of the problem, the number of Newton iterations needed per step to solve the discretized system of equations can, however, vary dramatically from a few to several hundred. Another popular approach based on the classical conjugate gradient method, known as the GMRES (Generalized Minimum Residual) algorithm is investigated. The GMRES algorithm was used in the past by a number of researchers for solving steady viscous and inviscid flow problems with considerable success. Here, the suitability of this algorithm is investigated for solving the system of nonlinear equations that arise in unsteady Navier-Stokes solvers at each time step. Unlike the Newton method which attempts to drive the error in the solution at each and every node down to zero, the GMRES algorithm only seeks to minimize the L2 norm of the error. In the GMRES algorithm the changes in the flow properties from one time step to the next are assumed to be the sum of a set of orthogonal vectors. By choosing the number of vectors to a reasonably small value N (between 5 and 20) the work required for advancing the solution from one time step to the next may be kept to (N+1) times that of a noniterative scheme. Many of the operations required by the GMRES algorithm such as matrix-vector multiplies, matrix additions and subtractions can all be vectorized and parallelized efficiently

Application of Extended Messinger Models to Complex Geometries

Since, ice accretion can significantly degrade the performance and the stability of an airborne vehicle, it is imperative to be able to model it accurately. While ice accretion studies have been performed on airplane wings and helicopter blades in abundance, there are few that attempt to model the process on more complex geometries such as fuselages. This paper proposes a methodology that extends an existing in-house Extended Messinger solver to complex geometries by introducing the capability to work with unstructured grids and carry out spatial surface streamwise marching. For the work presented here commercial solvers such as STAR-CCM+ and ANSYS Fluent are used for the flow field and droplet dispersed phase computations. The ice accretion is carried out using an in-house icing solver called GT-ICE. The predictions by GT-ICE are compared to available experimental data, or to predictions by other solvers such as LEWICE and STAR-CCM+. Three different cases with varying levels of complexity are presented. The first case considered is a commercial transport airfoil, followed by a three-dimensional MS(1)-317 swept wing. Finally, ice accretion calculations performed on a Robin fuselage have been discussed. Good agreement with experimental data, where applicable, is observed. Differences between the ice accretion predictions by different solvers have been discussed