Universitat Politècnica de Catalunya. Departament de Física, Universitat Politècnica de Catalunya. Doctorat en Física Computacional i Aplicada, Universitat Politècnica de Catalunya. Departament de Màquines i Motors Tèrmics, Barcelona Supercomputing Center, Universitat Politècnica de Catalunya. TUAREG - Turbulence and Aerodynamics in Mechanical and Aerospace Engineering Research Group, Miró Jané, Arnau, Eiximeno Franch, Benet, Rodríguez Pérez, Ivette María, Lehmkuhl Barba, Oriol, Universitat Politècnica de Catalunya. Departament de Física, Universitat Politècnica de Catalunya. Doctorat en Física Computacional i Aplicada, Universitat Politècnica de Catalunya. Departament de Màquines i Motors Tèrmics, Barcelona Supercomputing Center, Universitat Politècnica de Catalunya. TUAREG - Turbulence and Aerodynamics in Mechanical and Aerospace Engineering Research Group, Miró Jané, Arnau, Eiximeno Franch, Benet, Rodríguez Pérez, Ivette María, and Lehmkuhl Barba, Oriol
A direct numerical simulation of a three-dimensional diffuser at Reynolds number Re = 10,000 (based on inlet bulk velocity) has been performed using a low-dissipation finite element code. The geometry chosen for this work is the Stanford diffuser, introduced by Cherry et al. (Int. J. Heat Fluid Flow 29:803–811, 2008). Results have been exhaustively compared with the published data with a quite good agreement. Additionally, further turbulent statistics have been provided such as the Reynolds stresses or the turbulent kinetic energy. A proper orthogonal decomposition and a dynamic mode decomposition analyses of the main flow variables have been performed to identify the main characteristics of the large-scale motions. A combined, self-induced movement of the large-scales has been found to originate in the top-right expansion corner with two clear features. A low-frequency diagonal cross-stream travelling wave first reported by Malm et al. (J. Fluid Mech. 699:320–351, 2012), has been clearly identified in the spatial modes of the stream-wise velocity components and the pressure, associated with the narrow band frequency of St¿[0.083,0.01] . This movement is caused by the geometrical expansion of the diffuser in the cross-stream direction. A second low-frequency trait has been identified associated with the persisting secondary flows and acting as a back and forth global accelerating-decelerating motion located on the straight area of the diffuser, with associated frequencies of St<0.005 . The smallest frequency observed in this work has been St=0.0013 . This low-frequency observed in the Stanford diffuser points out the need for longer simulations in order to obtain further turbulent statistics., The research leading to this work has been partially funded by the European Project NextSim which has received funding from the European High-Performance Computing Joint Undertaking (JU) under grant agreement No 956104 and co-founded by the Spanish Agencia Estatal de Investigacion (AEI) under grant agreement PCI2021-121962. Benet Eiximeno also acknowledges the financial support by the Ministerio de Economía y Competitividad, Secretaría de Estado de Investigación, Desarrollo e Innovación, Spain (Refs: PID2020-116937RB-C21 and PID2020-116937RB-C22). Oriol Lehmkuhl has been partially supported by a Ramon y Cajal postdoctoral contract (Ref: RYC2018-025949-I). He also acknowledges the support of the European Project HiFi-TURB which has received funding from the European Union’s Horizon 2020 research and innovation programme under grant agreement No 814837. We also acknowledge the Barcelona Supercomputing Center for awarding us access to the MareNostrum IV machine based in Barcelona, Spain. The authors acknowledge the support of Departament de Recerca i Universitats de la Generalitat de Catalunya to the Research Group Large-scale Computational Fluid Dynamics (Code: 2021 SGR 00902)., Peer Reviewed, Postprint (published version)