Large eddy simulation using a dynamic mixing length subgrid-scale model

A novel dynamic mixing length (DML) subgrid-scale model for large eddy simulations is proposed in this work to improve the cutoff length of the Smagorinsky model. The characteristic mixing length (or the characteristic wave number) is dynamically estimated for the subgrid-scale fluctuation of turbulence by the cutoff wave-number, kc, and the dissipation wave-number, kd. The dissipation wave number is derived from the kinetic energy spectrum equation and the dissipation spectrum equation. To prove the promise of the DML model, this model is used to simulate the lid-driven cubical cavity with max-velocity-based Reynolds numbers 8850 and 12,000, the channel flows with friction-velocity-based Reynolds numbers 180, 395, 590, and 950, and the turbulent flow past a square cylinder at the higher Reynolds number 21,400, respectively, compared with the Smagorinsky model and Germano et al.'s dynamic Smagorinsky model. Different numerical experiments with different Reynolds numbers show that the DML model can be used in simulations of flows with a wide range of Reynolds numbers without the occurrence of singular values. The DML model can alleviate the dissipation of the Smagorinsky model without the loss of its robustness. The DML model shows some advantages over Germano et al.'s dynamic Smagorinsky model in its high stability and simplicity of calculation because the coefficient of the DML model always stays positive. The characteristic mixing length in the DML model reflects the subgrid-scale fluctuation of turbulence in nature and thus the characteristic mixing length has a spatial and temporal distribution in turbulent flow. Copyright (C) 2011 John Wiley & Sons, Ltd.