Minimum-dissipation model for large-eddy simulation using symmetry-preserving discretization in OpenFOAM

J Sun,R.W.C.P Verstappen

2023-11-03

Abstract:The minimum-dissipation model is applied to channel flow up to $Re_\tau = 2000$, flow past a circular cylinder at $Re=3900$, and flow over periodic hills at $Re=10595$. Numerical simulations were performed in OpenFOAM which is based on the finite volume methods. We used both symmetry-preserving and standard second-order accurate discretization methods in OpenFOAM on structured meshes. The results are compared to DNS and experimental data. The results of channel flow demonstrate a static QR model performs equally well as the dynamic models while reducing the computational cost. The model constant of $C=0.024$ gives the most accurate prediction, and the contribution of the sub-grid model decreases with the increase of the mesh resolution and becomes very small (less than 0.2 molecular viscosity) if a fine mesh is used. Furthermore, the QR model is able to predict the mean and rms velocity accurately up to $Re_\tau = 2000$ without a wall damping function. The symmetry-preserving discretization outperforms the standard OpenFOAM discretization at $Re_\tau=1000$. The results for the flow over a cylinder show that the mean velocity, drag coefficient, and lift coefficient are in good agreement with the experimental data and the central difference schemes conjugated with the QR model predict better results. The various comparisons carried out for flows over periodic hills demonstrate the need to use central difference schemes in OpenFOAM in combination with the minimum dissipation model. The best model constant is again $C=0.024$. The single wind turbine simulation shows that the QR model is capable of predicting accurate results in complex rotating scenarios.

Fluid Dynamics,Mathematical Physics

What problem does this paper attempt to address?

The problem that this paper attempts to solve is how to effectively describe the sub - grid scale (SGS) turbulent fluxes in large - eddy simulation (LES), especially by applying the minimum - dissipation model to reduce the computational cost while maintaining the accuracy of the model. Specifically, the paper focuses on the following aspects:

1. **Application of the minimum - dissipation model**:
- Research on the performance of static and dynamic minimum - dissipation models in channel flows at different Reynolds numbers.
- Evaluate the influence of the model constant \(C\) on the model performance, especially its performance at a high Reynolds number \(Re_{\tau}=2000\).

2. **Comparison of numerical discretization methods**:
- Compare the performance of the symmetry - preserving discretization method and the standard OpenFOAM discretization method under different flow conditions.
- In particular, evaluate the advantages of the symmetry - preserving discretization method in channel flows at \(Re_{\tau} = 1000\).

3. **Application under complex geometric conditions**:
- Research on the applicability of the minimum - dissipation model in flow past a circular cylinder and flow over periodic hills.
- Evaluate the influence of different numerical schemes (such as central difference, upwind - biased central difference, etc.) on the flow behavior.

4. **Selection of model constants**:
- Determine the optimal model constant \(C\) so that the minimum - dissipation model can provide accurate results under different flow conditions.
- For example, in channel flows, the model constant \(C = 0.024\) is considered optimal.

Through these studies, the paper aims to verify the effectiveness and reliability of the minimum - dissipation model under different flow conditions and explore its potential in practical engineering applications.