Xiang Yang

Integral wall model for large eddy simulations of wall-bounded turbulent flows

A new approach for wall modeling in Large-Eddy-Simulations (LES) is proposed and tested in various applications. To properly include near-wall physics while preserving the basic economy of equilibrium-type wall models, we adopt the classical integral method of von Karman and Pohlhausen (VKP). A velocity profile with various parameters is proposed as an alternative to numerical integration of the boundary layer equations in the near-wall zone. The profile contains a viscous or roughness sublayer and a logarithmic layer with an additional linear term that can account for inertial and pressure gradient effects. Similar to the VKP method, the assumed velocity profile coefficients are determined from appropriate matching conditions and physical constraints. The proposed integral wall-modeled LES (iWMLES) method is tested in the context of a pseudo-spectral code for fully developed channel flow with a dynamic Lagrangian subgrid model as well as in a finite-difference LES code including the immersed boundary method and the dynamic Vreman eddy-viscosity model. Test cases include a fully developed half-channel at various Reynolds numbers, a fully developed channel flow with unresolved roughness, a standard developing turbulent boundary layer flows over smooth plates at various Reynolds numbers, over plates with unresolved roughness, and a case with resolved roughness elements consisting of an array of wall-mounted cubes. The comparisons with data show that the proposed iWMLES method provides accurate predictions of near-wall velocity profiles in LES while, similarly to equilibrium wall models, its cost remains independent of Reynolds number and is thus significantly lower compared to existing zonal or hybrid wall models. A sample application to flow over a surface with truncated cones (representing idealized barnacle-like roughness elements) is also presented, which illustrates effects of subgrid scale roughness when combined with resolved roughness elements.

A Web services accessible database of turbulent channel flow and its use for testing a new integral wall model for LES

abstract The output from a direct numerical simulation (DNS) of turbulent channel flow at Reτ ≈ 1000 is used to construct a publicly and Web services accessible, spatio-temporal database for this flow. The simulated channel has a size of 8πh × 2h × 3πh, where h is the channel half-height. Data are stored at 2048 × 512 × 1536 spatial grid points for a total of 4000 time samples every 5 time steps of the DNS. These cover an entire channel flow-through time, i.e. the time it takes to traverse the entire channel length 8πh at the mean velocity of the bulk flow. Users can access the database through an interface that is based on the Web services model and perform numerical experiments on the slightly over 100 terabytes (TB) DNS data on their remote platforms, such as laptops or local desktops. Additional technical details about the pressure calculation, database interpolation, and differentiation tools are provided in several appendices. As a sample application of the channel flow database, we use it to conduct an a-priori test of a recently introduced integral wall model for large eddy simulation of wall-bounded turbulent flow. The results are compared with those of the equilibrium wall model, showing the strengths of the integral wall model as compared to the equilibrium model.

Recycling inflow method for simulations of spatially evolving turbulent boundary layers over rough surfaces

ABSTRACTThe technique by Lund et al. to generate turbulent inflow for simulations of developing boundary layers over smooth flat plates is extended to the case of surfaces with roughness elements. In the Lund et al. method, turbulent velocities on a sampling plane are rescaled and recycled back to the inlet as inflow boundary condition. To rescale mean and fluctuating velocities, appropriate length scales need be identified and for smooth surfaces, the viscous scale lν = ν/uτ (where ν is the kinematic viscosity and uτ is the friction velocity) is employed for the inner layer. Different from smooth surfaces, in rough wall boundary layers the length scale of the inner layer, i.e. the roughness sub-layer scale ld, must be determined by the geometric details of the surface roughness elements and the flow around them. In the proposed approach, it is determined by diagnosing dispersive stresses that quantify the spatial inhomogeneity caused by the roughness elements in the flow. The scale ld is used for rescali...

Aerodynamic Properties of Rough Surfaces with High Aspect-Ratio Roughness Elements: Effect of Aspect Ratio and Arrangements

We examine the effect of varying roughness-element aspect ratio on the mean velocity distributions of turbulent flow over arrays of rectangular-prism-shaped elements. Large-eddy simulations (LES) in conjunction with a sharp-interface immersed boundary method are used to simulate spatially-growing turbulent boundary layers over these rough surfaces. Arrays of aligned and staggered rectangular roughness elements with aspect ratio >1 are considered. First the temporally- and spatially-averaged velocity profiles are used to illustrate the aspect-ratio effects. For aligned prisms, the roughness length (

On the Mean Flow Behaviour in the Presence of Regional-Scale Surface Roughness Heterogeneity

z_\mathrm{o}

Statistics of turbulence in the energy-containing range of Taylor-Couette compared to canonical wall-bounded flows

zo) and the friction velocity (

Large eddy simulations and parameterisation of roughness element orientation and flow direction effects in rough wall boundary layers

u_*