Abstract
This work investigates a high‐order numerical method which is suitable for performing large‐eddy simulations, particularly those containing wall‐bounded regions which are considered on stretched curvilinear meshes. Spatial derivatives are represented by a sixth‐order compact approximation that is used in conjunction with a tenth‐order non‐dispersive filter. The scheme employs a time‐implicit approximately factored finite‐difference algorithm, and applies Newton‐like subiterations to achieve second‐order temporal and sixth‐order spatial accuracy. Both the Smagorinsky and dynamic subgrid‐scale stress models are incorporated in the computations, and are used for comparison along with simulations where no model is employed. Details of the method are summarized, and a series of classic validating computations are performed. These include the decay of compressible isotropic turbulence, turbulent channel flow, and the subsonic flow past a circular cylinder. For each of these cases, it was found that the method was robust and provided an accurate means of describing the flowfield, based upon comparisons with previous existing numerical results and experimental data. Published in 2003 by John Wiley & Sons, Ltd.