Parallelization of the Implicit RPLUS Algorithm

The multiblock reacting Navier-Stokes flow solver RPLUS2D was modified for parallel implementation. Results for non-reacting flow calculations of this code indicate parallelization efficiencies greater than 84% are possible for a typical test problem. Results tend to improve as the size of the problem increases. The convergence rate of the scheme is degraded slightly when additional artificial block boundaries are included for the purpose of parallelization. However, this degradation virtually disappears if the solution is converged near to machine zero. Recommendations are made for further code improvements to increase efficiency, correct bugs in the original version, and study decomposition effectiveness

Implementation of a two-equation k-omega turbulence model in NPARC

The implementation of a two-equation k-omega turbulence model into the NPARC flow solver is described. Motivation for the selection of this model is given, major code modifications are outlined, new imputs to the code are described, and results are presented for several validation cases: an incompressible flow over a smooth flat plate, a subsonic diffuser flow, and a shock-induced separated flow. Comparison of results with the k-epsilon model indicate that the k-omega model predicts simple flows equally well whereas, for adverse pressure gradient flows, the k-omega model outperforms the other turbulence models in NPARC

A Discontinuous Galerkin Chimera scheme

The Chimera overset method is a powerful technique for modeling fluid flow associated with complex engineering problems using structured meshes. The use of structured meshes has enabled engineers to employ a number of high-order schemes, such as the WENO and compact differencing schemes. However, the large stencil associated with these schemes can significantly complicate the inter-grid communication scheme and hole cutting procedures. This paper demonstrates a methodology for using the Discontinuous Galerkin (DG) scheme with Chimera overset meshes. The small stencil of the DG scheme makes it particularly suitable for Chimera meshes as it simplifies the inter-grid communication scheme as well as hole cutting procedures. The DG-Chimera scheme does not require a donor interpolation method with a large stencil because the DG scheme represents the solution as cell local polynomials. The DG-Chimera method also does not require the use of fringe points to maintain the interior stencil across inter-grid boundaries. Thus, inter-grid communication can be established as long as the receiving boundary is enclosed by or abuts the donor mesh. This makes the inter-grid communication procedure applicable to both Chimera and zonal meshes. Details of the DG-Chimera scheme are presented, and the method is demonstrated on a set of two-dimensional inviscid flow problems

Assessment of Computational Fluid Dynamics (CFD) Models for Shock Boundary-Layer Interaction

A workshop on the computational fluid dynamics (CFD) prediction of shock boundary-layer interactions (SBLIs) was held at the 48th AIAA Aerospace Sciences Meeting. As part of the workshop numerous CFD analysts submitted solutions to four experimentally measured SBLIs. This paper describes the assessment of the CFD predictions. The assessment includes an uncertainty analysis of the experimental data, the definition of an error metric and the application of that metric to the CFD solutions. The CFD solutions provided very similar levels of error and in general it was difficult to discern clear trends in the data. For the Reynolds Averaged Navier-Stokes methods the choice of turbulence model appeared to be the largest factor in solution accuracy. Large-eddy simulation methods produced error levels similar to RANS methods but provided superior predictions of normal stresses

On the Use of Optimization Techniques for Turbulence Model Calibration

This article discusses the use of numerical optimization procedures to aid in the calibration of turbulence model coefficients. Such methods would increase the rigor and repeatability of the calibration procedure by requiring clearly defined and objective optimization metrics, and could be used to identify unique combinations of coefficient values for specific flow problems. The approach is applied to the re-calibration of an explicit algebraic Reynolds stress model for the incompressible planar mixing layer using the Nelder-Mead simplex algorithm and a micro-genetic algorithm with minimally imposed constraints. Three composite fitness functions, each based upon the error in the mixing layer growth rate and the normal and shear components of the Reynolds stresses, are investigated. The results demonstrate a significant improvement in the target objectives through the adjustment of three pressure-strain coefficients. Adjustments of additional coefficients provide little further benefit. Issues regarding the effectiveness of the fitness functions and the efficiency of the optimization algorithms are also discussed