STURM 4 progress meeting Clamart 13 Oct 2009

STURM 4 -progress meeting, Clamart 13 Oct 2009 EDF – U o Man collaborations on Unstructured meshes for DNS and LES * 2 conf. papers presented by Y. Addad (Uo. M) and R. Howard (EDF) * Int. Symp on Convective Heat and Mass Transfer in Sustainable Energy, April 26 – 1 May, 2009 Hammamet, Tunisia * France-Japan joint Seminar on Thermal fatigue, 5 th to 6 th October 2009, Tokyo, Japan 1

International Symposium on Convective Heat and Mass Transfer in Sustainable Energy April 26 – 1 May, 2009 Yasmine Hammamet, Tunisia THE EFFECT OF DIFFERENTS GRIDS AND DIFFERENT LES MODELS ON TURBULENCE STATISTICS IN CHANNEL FLOWS Richard Howard and Yacine Addad

DIFFERENTS GRIDS AND DIFFERENT LES MODELS IN CHANNEL FLOWS Stretched uniform hybrid Box Reynolds (u tau) 180 Code Saturne http///www. code-saturne. org Collocalised, Cell centered, 2 nd order central differences 3

DIFFERENTS GRIDS AND DIFFERENT LES MODELS IN CHANNEL FLOWS Production Dissipation Subgrid model contributions 4 Viscous Pressure Turbulent

DIFFERENTS GRIDS AND DIFFERENT LES MODELS IN CHANNEL FLOWS Second derivatives Gradients of turbulent viscosity Products of velocity gradients 5

DIFFERENTS GRIDS AND DIFFERENT LES MODELS IN CHANNEL FLOWS Stretched 6

DIFFERENTS GRIDS AND DIFFERENT LES MODELS IN CHANNEL FLOWS Stretched WALE 7 Term I (second derivatives) Term II (gradients of turbulent viscosity) Term III (products of gradients) Smagorinsky

DIFFERENTS GRIDS AND DIFFERENT LES MODELS IN CHANNEL FLOWS Uniform 8

DIFFERENTS GRIDS AND DIFFERENT LES MODELS IN CHANNEL FLOWS Hybrid 9

DIFFERENTS GRIDS AND DIFFERENT LES MODELS IN CHANNEL FLOWS • The dominant term is Term. III (products of velocity gradients which includes a kind subgrid dissipation term) • The WALE model dissipates the resolved fluctuations in v and w more than the other models. • The hybrid grid show numerical noise particularly at the channel centreline. The sloping faces of the non structured grid act to transfer the local mean velocity field between each velocity component. 10

Thermal mixing in T-Junction: Hexa Grid, STAR V 4 LES 3. 1 Dh Grid Cells= 2. 56 M. 3 Dc LES RUNS: 13 Dc SGS Model: Smago. Pr. SGS=0. 9 (Default). URANS RUNS: Non-linear k- model Suga et al. 2006. Numerical schemes: 11 2 nd order in space. 2 nd order in time.

Conclusions and future work - URANS model: - Fails to capture the complex features of the flow. - Most information needed for thermal fatigue studies is lost. - Not able to capture high frequency events. - Can this be improved or other approaches tested (example DES) ? ? - LES with Unstructured grids and Professional Software: - Second order accuracy seems OK. - Mesh is extremely important, adapted to (1/10) large eddy scale ? => need unstructured mesh - Work in progress: How about tetrahedral and polyhedral cells, different SGS models? ? Acknowledgements: 12 UK research council project “Keeping the Nuclear Option Open” (KNOO) FR ANR-CIS STURM 4