The approximate deconvolution model (ADM) for the large-eddy simulation of incompressible flows is detailed and applied to turbulent channel flow. With this approach an approximation of the unfiltered solution is obtained by repeated filtering. Given a good approximation of the unfiltered solution, the nonlinear terms of the filtered Navier–Stokes equations can be computed directly. The effect of nonrepresented scales is modeled by a relaxation regularization involving a secondary filter operation. Large-eddy simulations are performed for incompressible channel flow at Reynolds numbers based on the friction velocity and the channel half-width of Reτ=180 and Reτ=590. Both simulations compare well with direct numerical simulation (DNS) data and show a significant improvement over results obtained with classical subgrid scale models such as the standard or the dynamic Smagorinsky model. The computational cost of ADM is lower than that of dynamic models or the velocity estimation model.

High-resolution simulations are presented of particle-driven gravity currents in the lock-exchange configuration. The study concentrates on dilute flows with small density differences between particle-laden and clear fluid. Moreover, particles are considered which have negligible inertia, and which are much smaller than the smallest length scales of the buoyancy-induced fluid motion. For the mathematical description of the particulate phase a Eulerian approach is employed with a transport equation for the local particlenumber density. The governing equations are integrated numerically with a high-order mixed spectral/ spectral-element technique. In the analysis of the results, special emphasis is placed on the sedimentation of particles and the influence of particle settling on the flow dynamics. Time-dependent sedimentation profiles at the channel floor are presented which agree closely with available experimental data. A detailed study is conducted of the balance between the various components of the energy budget of the flow, i.e. the potential and kinetic energy, and the dissipative losses. Furthermore, the simulation results, along with a modified Shields criterion, are used to show that resuspension of sediment back into the particle-driven current is unlikely to occur in the cases considered. Two-dimensional (2D) and three-dimensional (3D) computations are compared which reveal that, for the present configuration, a 2D model can predict reliably the flow development at early times. However, concerning the long-time evolution of the flow, more substantial differences exist between 2D and 3D simulations. Ó

The gravitational settling of an initially random suspension of small solid particles in homogeneous turbulence is investigated numerically. The simulations are based on a pseudospectral method to solve the fluid equations combined with a Lagrangian point-particle model for the particulate phase ͑Eulerian-Lagrangian approach͒. The focus is on the enhancement of the mean particle settling velocity in a turbulent carrier fluid, as compared to the settling velocity of a single particle in quiescent fluid. Results are presented for both one-way coupling, when the fluid flow is not affected by the presence of the particles, and two-way coupling, when the particles exert a feedback force on the fluid. The first case serves primarily for validation purposes. In the case with two-way coupling, it is shown that the effect of the particles on the carrier fluid involves an additional increase in their mean settling velocity compared to one-way coupling. The underlying physical mechanism is analyzed, revealing that the settling velocity enhancement depends on the particle loading, the Reynolds number, and the dimensionless Stokes settling velocity if the particle Stokes number is about unity. Also, for particle volume fractions ⌽ v տ 10 −5 , a turbulence modification is observed. Furthermore, a direct comparison with recent experimental studies by Aliseda et al. ͓J. Fluid Mech. 468, 77 ͑2002͔͒ and Yang and Shy ͓J. Fluid Mech. 526, 171 ͑2005͔͒ is performed for a microscale Reynolds number Re Ϸ 75 of the turbulent carrier flow.

No abstract

A formulation of the approximate deconvolution model (ADM) for the large-eddy simulation (LES) of compressible flows in complex geometries is detailed. The model is applied to supersonic compression ramp flow where shock-turbulence interaction occurs. With the ADM approach an approximation to the unfiltered solution is obtained from the filtered solution by a series expansion involving repeated filtering. Given a sufficiently good approximation of the unfiltered solution at a time instant, the flux terms of the underlying filtered transport equations can be computed directly, avoiding the need to explicitly compute subgrid-scale closures. The effect of nonrepresented scales is modeled by a relaxation regularization involving a secondary filter operation and a dynamically estimated relaxation parameter. Results of the large-eddy simulation of the turbulent supersonic boundary layer along a compression ramp compare well with filtered DNS data. The filtered shock solution is correctly predicted by the ADM procedure, demonstrating that turbulent and nonturbulent subgrid-scales are properly modeled. We found that a computationally expensive shock-capturing technique was not necessary for stable integration. As a consequence, the computational effort for simulations with ADM is approximately as large as for a coarse-grid DNS with a hybrid compact-upwind-ENO scheme, since the additional computational cost for the subgrid-scale model is more than compensated due to the fact that in the LES flux-derivatives can be computed by linear central finite differences on the entire domain.