Conservation properties of the mass, momentum, and kinetic energy equations for incompressible flow are specified as analytical requirements for a proper set of discrete equations. Existing finite difference schemes in regular and staggered grid systems are checked for violations of the conservation requirements and a few important discrepancies are pointed out. In particular, it is found that none of the existing higher order schemes for a staggered mesh system simultaneously conserve mass, momentum, and kinetic energy. This deficiency is corrected through the derivation of a general family of fully conservative higher order accurate finite difference schemes for staggered grid systems. Finite difference schemes in a collocated grid system are also analyzed, and a violation of kinetic energy conservation is revealed. The predicted conservation properties are demonstrated numerically in simulations of inviscid white noise, performed in a two-dimensional periodic domain. The proposed fourth order schemes in a staggered grid system are generalized for the case of a non-uniform mesh, and the resulting scheme is used to perform large eddy simulations of turbulent channel flow.

This article reviews state-of-the-art adaptive, multiresolution wavelet methodologies for modeling and simulation of turbulent flows with various examples. Different numerical methods for solving the Navier-Stokes equations in adaptive wavelet bases are described. We summarize coherent vortex extraction methodologies, which utilize the efficient wavelet decomposition of turbulent flows into space-scale contributions, and present a hierarchy of wavelet-based turbulence models. Perspectives for modeling and computing industrially relevant flows are also given.

Large-eddy simulation (LES) with regular explicit filtering is investigated. The filtered-scale stress due to the explicit filtering is here partially reconstructed using the tensor-diffusivity model: It provides for backscatter along the stretching direction(s), and for global dissipation, both also attributes of the exact filtered-scale stress. The necessary LES truncations (grid and numerical method) are responsible for an additional subgrid-scale stress. A natural mixed model is then the tensor-diffusivity model supplemented by a dynamic Smagorinsky term. This model is reviewed, together with useful connections to other models, and is tested against direct numerical simulation (DNS) of turbulent isotropic decay starting with ${\mathrm{Re}}_{\lambda }\text{=90}$ (thus moderate Reynolds number): LES started from a ${256}^3$ DNS truncated to ${64}^3$ and Gaussian filtered. The tensor-diffusivity part is first tested alone; the mixed model is tested next. Diagnostics include energy decay, enstrophy decay, and energy spectra. After an initial transient of the dynamic procedure (observed with all models), the mixed model is found to produce good results. However, despite expectations based on favorable a priori tests, the results are similar to those obtained when using the dynamic Smagorinsky model alone in LES without explicit filtering. Nevertheless, the dynamic mixed model appears as a good compromise between partial reconstruction of the filtered-scale stress and modeling of the truncations effects (incomplete reconstruction and subgrid-scale effects). More challenging ${48}^3$ LES are also done: Again, the results of both approaches are found to be similar. The dynamic mixed model is also tested on the turbulent channel flow at ${\mathrm{Re}}_{\tau }\text{=395.}$ The tensor-diffusivity part must be damped close to the wall in order to avoid instabilities. Diagnostics are mean profiles of velocity, stress, dissipation, and reconstructed Reynolds stresses. The velocity profile obtained using the damped dynamic mixed model is slightly better than that obtained using the dynamic Smagorinsky model without explicit filtering. The damping used so far is however crude, and this calls for further work.