Flow solver: Proust / TurbFlow

The equations solved are the 3D, unsteady, compressible, Reynolds aver­aged (RANS) or spatially filtered (LES), Navier-Stokes equations cast in the absolute frame where the laminar viscosity is assumed constant or calculated by the Sutherland’s law.

1 Spatial discretization. The space discretization is based on a MUSCL finite volume formulation with moving structured meshes, which uti­lizes vertex variable storage. The convective flixes are evaluated using a 3rd order upwind scheme (Van Leer’s Flux Vector Splitting with the Hanel cor­rection, Roe’s Approximate Riemann Solver, or Liou’s Advection Upwind Splitting Method), or 4th order centered scheme (for LES). An hybrid method combining the advantages of the central scheme in subsonic regions with the properties of the upwind scheme through discontinuities has been introduced to reduce the numerical losses in very low Mach number regions. The viscous terms are computed by a second order centered scheme. The resulting semi­discrete scheme is integrated in time using an explicit five steps Runge-Kutta time marching algorithm.

2 Boundary conditions. Compatibility relations are used to take into account physical boundary conditions. The outgoing characteristics are retained, since these provide information from inside the domain. The incom­ing characteristics, on the other hand, are replaced by physical boundary con­ditions, i. e. total pressure, total temperature and ft>w angles for a subsonic inlet, static pressure for a subsonic outlet, zero normal velocity component for a slip wall and zero velocity and heat fhx for an adiabatic wall. Ghost cells for which the equations are not solved are built around the domain to simulate ge­ometrical boundary conditions, like periodicity and symmetry. Non refhctive boundary conditions are implemented by retaining the equations associated to the incoming characteristics, in which the wave velocity is fixed to zero to prohibit propagation directed into the computational domain.

3 Turbulence Models. For the RANS approach, turbulence is taken into account by a k — и model. Two transport equations are implemented, gov­erning the turbulent energy k and the dissipation u. The evaluation of the Reynolds tensor and of the turbulent viscosity are carried out by different tur­bulence models: the linear model of Wilcox, 1993b, the low-Reynolds model

Figure 1. VKI turbine. Left: Mach numbers contour’s. Right: Comparison of density iso­contours [Upper : linear k-w model; Middle : experimental data (ONERA); Lower : non-linear k-w model]

of Wilcox, 1993a, the non-linear model of Shih et al., 1995 and the non-linear model of Craft et al., 1996.

For the LES, the subgrid scales are represented by a viscosity, computed using the auto-adaptive model of Casalino et al., 2003. This formulation enables an effective evaluation of the subgrid-scale viscosity, even for complex geome­tries.