DNS 1-3 Computational Details: Difference between revisions
No edit summary |
No edit summary |
||
| Line 17: | Line 17: | ||
terms in the budget equations are computed, in particular if there are differences to the proposed | terms in the budget equations are computed, in particular if there are differences to the proposed | ||
approach in Introduction. | approach in Introduction. | ||
==References== | |||
#<div id="1">'''Ohlsson, J., Schlatter, P., Fischer P.F. and Henningson, D.S. (2009):''' DNS of three-dimensional separation in turbulent diffuser flows. In ''Advances in Turbulence XII'', Proceedings of the 12th EUROMECH European Turbulence Conference, Marburg. Springer Proceedings in Physics, Vol. 132, ISBN 978-3-642-03084-0</div> | |||
<br/> | <br/> | ||
---- | ---- | ||
Revision as of 10:21, 13 February 2021
Computational Details
Computational approach
Alya is a parallel multi-physics/multiscale simulation code developed at the Barcelona Supercomputing Centre to run efficiently on high-performance computing environments. The convective term is discretized using a Galerkin finite element (FEM) scheme recently proposed [11,13], which conserves linear and angular momentum, and kinetic energy at the discrete level. Both second- and third-order spatial discretizations are used.
Neither upwinding nor any equivalent momentum stabilization is employed. In order to use equalorder elements, numerical dissipation is introduced only for the pressure stabilization via a fractional step scheme [12], which is similar to approaches for pressure-velocity coupling in unstructured, collocated finite-volume codes [13]. The set of equations is integrated in time using a third order Runge-Kutta explicit method combined with an eigenvalue-based time-step estimator [14]. This approach has been shown to be significantly less dissipative than the traditional stabilized FEM approach [15]. Thus, is an optimal methodology for high-fidelity simulations of complex flows as the ones required in the present project.
Spatial and temporal resolution, grids
A mesh resulting in approximately 250 million degrees of freedom (DoF). With a stretched grid, the maximum grid resolution in the duct centre is reported at Δz + = 11.6, Δy + = 13.2 and Δx += 19.5. Correspondingly, the wall resolution (in terms of the first grid point) is reported as z += 0.074, y + = 0.37 in the spanwise and normal directions, respectively, is used. This resolution was deemed sufficient to compute the flow in the diffuser and is based on the previous work of Ohlsson et al. (2010). For the temporal resolution a third order explicit Runge Kutta method using a dynamic time stepping with CFL below 0.9 has been used.
Computation of statistical quantities
Describe how the averages and correlations are obtained from the instantaneous results and how terms in the budget equations are computed, in particular if there are differences to the proposed approach in Introduction.
References
- Ohlsson, J., Schlatter, P., Fischer P.F. and Henningson, D.S. (2009): DNS of three-dimensional separation in turbulent diffuser flows. In Advances in Turbulence XII, Proceedings of the 12th EUROMECH European Turbulence Conference, Marburg. Springer Proceedings in Physics, Vol. 132, ISBN 978-3-642-03084-0
Contributed by: Oriol Lehmkuhl, Arnau Miro — BSC
© copyright ERCOFTAC 2024