Artificial and Eddy Viscosity in Large Eddy Simulation Part 1: Temporal and Spatial Schemes
Fuente:
arXiv
Salvato in:
| Autori principali: | , |
|---|---|
| Natura: | Preprint |
| Pubblicazione: |
2025
|
| Soggetti: | |
| Accesso online: | |
| Tags: |
Aggiungi Tag
Nessun Tag, puoi essere il primo ad aggiungerne!!
|
| _version_ | 1866913712362749952 |
|---|---|
| author | Sun, Jing Verstappen, Roel |
| author_facet | Sun, Jing Verstappen, Roel |
| contents | We propose a novel method to quantify artificial dissipation in large eddy simulation. Here, artificial dissipation is defined as the residual of the discrete turbulent kinetic energy (TKE) equation. This method is applied to turbulent channel flow at Reynolds number 180 and 550 using a minimum-dissipation model within the OpenFOAM framework, employing both symmetry-preserving and standard OpenFOAM discretizations. Our analysis reveals that artificial viscosity can both produce and dissipate TKE. To quantify the interaction between sub-grid-scale model contributions and numerical error, we introduce viscous activity parameters, which also allow for evaluating the balance between molecular and non-molecular viscosity. Examining temporal discretization schemes, we find that all methods can deliver accurate flow predictions if an appropriate time step is used. For very small time steps, Forward-Euler is more efficient and accurate than Crank-Nicolson, though a better balance of accuracy and computational efficiency for larger time steps, making it preferable for large-scale applications. On collocated meshes, Van Kan's pressure correction method predicts better results than Chorin's. Our spatial discretization study shows that the symmetry-preserving scheme outperforms standard OpenFOAM schemes, as excessive artificial errors from spatial discretization can suppress the LES model. For the symmetry-preserving scheme, both eddy and artificial viscosities vary similarly, with eddy viscosity being dominant. The role of artificial viscosity changes with mesh refinement: on coarse meshes, it mainly dissipates TKE, while on finer meshes, it can produce TKE, counteracting eddy viscosity. The most accurate Reynolds shear stress predictions occur when non-molecular dissipation is nearly zero. For accurate mean flow kinetic energy, larger non-molecular dissipation is needed on finer meshes. |
| format | Preprint |
| id |
arxiv_https___arxiv_org_abs_2502_08154 |
| institution | arXiv |
| publishDate | 2025 |
| record_format | arxiv |
| spellingShingle | Artificial and Eddy Viscosity in Large Eddy Simulation Part 1: Temporal and Spatial Schemes Sun, Jing Verstappen, Roel Fluid Dynamics We propose a novel method to quantify artificial dissipation in large eddy simulation. Here, artificial dissipation is defined as the residual of the discrete turbulent kinetic energy (TKE) equation. This method is applied to turbulent channel flow at Reynolds number 180 and 550 using a minimum-dissipation model within the OpenFOAM framework, employing both symmetry-preserving and standard OpenFOAM discretizations. Our analysis reveals that artificial viscosity can both produce and dissipate TKE. To quantify the interaction between sub-grid-scale model contributions and numerical error, we introduce viscous activity parameters, which also allow for evaluating the balance between molecular and non-molecular viscosity. Examining temporal discretization schemes, we find that all methods can deliver accurate flow predictions if an appropriate time step is used. For very small time steps, Forward-Euler is more efficient and accurate than Crank-Nicolson, though a better balance of accuracy and computational efficiency for larger time steps, making it preferable for large-scale applications. On collocated meshes, Van Kan's pressure correction method predicts better results than Chorin's. Our spatial discretization study shows that the symmetry-preserving scheme outperforms standard OpenFOAM schemes, as excessive artificial errors from spatial discretization can suppress the LES model. For the symmetry-preserving scheme, both eddy and artificial viscosities vary similarly, with eddy viscosity being dominant. The role of artificial viscosity changes with mesh refinement: on coarse meshes, it mainly dissipates TKE, while on finer meshes, it can produce TKE, counteracting eddy viscosity. The most accurate Reynolds shear stress predictions occur when non-molecular dissipation is nearly zero. For accurate mean flow kinetic energy, larger non-molecular dissipation is needed on finer meshes. |
| title | Artificial and Eddy Viscosity in Large Eddy Simulation Part 1: Temporal and Spatial Schemes |
| topic | Fluid Dynamics |
| url | https://arxiv.org/abs/2502.08154 |