Wall-modeled large-eddy simulation based on spectral-element discretization

Timofey Mukha,Philipp Schlatter
2024-04-08
Abstract:This article analyses the simulation methodology for wall-modeled large-eddy simulations using solvers based on the spectral-element method (SEM). To that end, algebraic wall modeling is implemented in the popular SEM solver Nek5000. It is combined with explicit subgrid-scale (SGS) modeling, which is shown to perform better than the high-frequency filtering traditionally used with the SEM. In particular, the Vreman model exhibits a good balance in terms stabilizing the simulations, yet retaining good resolution of the turbulent scales. Some difficulties associated with SEM simulations on relatively coarse grids are also revealed: jumps in derivatives across element boundaries, lack of convergence for weakly formulated boundary conditions, and the necessity for the SGS model as a damper for high-frequency modes. In spite of these, state-of-the-art accuracy is achieved for turbulent channel flow and flat-plate turbulent boundary layer flow cases, proving the SEM to be a an excellent numerical framework for massively-parallel high-order WMLES.
Fluid Dynamics,Numerical Analysis
What problem does this paper attempt to address?
This paper aims to solve several key problems in Wall - Modeled Large - Eddy Simulation (WMLES), especially when using a solver based on the Spectral - Element Method (SEM) for simulation. Specifically, the paper focuses on the following aspects: 1. **Improving prediction accuracy**: The paper evaluates the accuracy of SEM - based WMLES in predicting turbulent channel flows and zero - pressure - gradient turbulent boundary - layer flows over a flat plate. By comparing with Direct Numerical Simulation (DNS) data, the effectiveness of the method is verified. 2. **Improving sub - grid - scale models**: The paper explores the performance of different Subgrid - Scale (SGS) models in WMLES, especially the Vreman model and the Sigma model. These models can stabilize the simulation while maintaining a good resolution of the turbulence scale. 3. **Dealing with challenges brought by numerical methods**: The paper points out some difficulties encountered when using SEM for simulation on relatively coarse grids, such as derivative jumps on element boundaries, convergence problems of weak - form boundary conditions, and the necessity of the SGS model as a high - frequency mode damper. 4. **Exploring different wall - boundary conditions**: Besides the traditional Neumann boundary condition, the paper also tests a boundary - condition method based on sub - grid viscosity. This new method can control the direction of wall - shear stress more precisely and can switch flexibly between wall - resolved and wall - modeled LES. 5. **Open - source implementation**: The paper provides an open - source WMLES code implementation so that other researchers can reproduce and further develop these methods. ### Specific problems and solutions - **Problem 1: How to obtain high - precision WMLES results on a coarse grid?** - **Solution**: By using advanced SGS models (such as the Vreman and Sigma models) and combining appropriate wall - modeling methods, high - precision results can be obtained on a coarse grid. The paper shows that even on a relatively coarse grid, an accuracy comparable to DNS data can be achieved. - **Problem 2: How to deal with the numerical stability problem of the spectral - element method near the wall?** - **Solution**: The paper proposes using a boundary - condition method based on sub - grid viscosity. This method can effectively reduce numerical instability and error. In addition, by adjusting the parameters of the SGS model, the resolution of the simulation can be improved while ensuring stability. - **Problem 3: How to evaluate the performance of different SGS models in WMLES?** - **Solution**: By comparing the performance of different SGS models in turbulent channel flows and flat - plate boundary - layer flows, the paper concludes that the Vreman model and the Sigma model have a good balance in WMLES, being able to stabilize the simulation and maintain good resolution. ### Conclusion Through the above research, the paper proves the effectiveness and accuracy of SEM - based WMLES in dealing with complex turbulence problems. Although there are some numerical challenges on coarse grids, by reasonably selecting SGS models and wall - boundary conditions, the stability and precision of the simulation can be significantly improved. These research results provide important references and guidance for future high - precision WMLES.