Numerical study on the turbulent mixing in channel with Large Eddy Simulation (LES) using spectral element method

Abstract Turbulent mixing is an important thermal hydraulic phenomenon in the reactor core rod bundles and it leads to strong momentum and energy transfer among adjacent subchannels in fuel assembly. Generally, turbulent mixing coefficient is usually referred to describe turbulent effect, which commonly sets to a constant or calculated utilizing Reynolds number dependent fitting correlations based on experimental data. In this paper, a high-scalable and high-performance spectral element method combining with high passed filtered Large Eddy Simulation (LES) model has been proposed to simulate the turbulent flow through parallel wall channel and square channel with a cylindrical rod for single-phase flow. In the parallel wall channel case, mesh sensitivity analysis and model validation were accomplished by comparing the simulation data with DNS references utilizing key parameters such as velocity profile and stresses in the near wall region. Meanwhile, for the square channel with a cylindrical rod case, instantaneous lateral velocity monitored in the gap center has been compared with the experimental data with different Reynolds numbers and geometric conditions. From the simulation approach, both the absolute oscillation amplitude and the oscillation frequency increase with the pitch diameter ratio. Finally, the simulation results of turbulent mixing coefficient dependent on Reynolds number were validated with both experimental and theoretic references to demonstrate the feasibility of spectral element method combining with LES model. The root mean square (RMS) value of lateral fluctuating velocity was adopted to reflect the effective mixing capability, predicting the turbulent mixing phenomena with a reasonable degree of accuracy from the specific calculation strategy. Fast Fourier Transform (FFT) also performed and got an agreement with reference data.

[1]  P. Fischer,et al.  High-Order Methods for Incompressible Fluid Flow , 2002 .

[2]  C. Lee,et al.  Effect of gap width on turbulent mixing of parallel flow in a square channel with a cylindrical rod , 2013 .

[3]  Elia Merzari,et al.  Accurate prediction of the wall shear stress in rod bundles with the spectral element method at high Reynolds numbers , 2014 .

[4]  Elia Merzari,et al.  Numerical simulation of the flow in wire-wrapped pin bundles: effect of pin-wire contact modeling. , 2012 .

[5]  Javier Jiménez,et al.  Reynolds number effects on the Reynolds-stress budgets in turbulent channels , 2008 .

[6]  M. L. Corradini,et al.  Potential Recovery Actions from a Severe Accident in a PWR: MELCOR Analysis of a Station Blackout Scenario , 2018, Nuclear Technology.

[7]  S. Qiu,et al.  DEM-CFD simulation of helium flow characteristics in randomly packed bed for fusion reactors , 2018, Progress in Nuclear Energy.

[8]  Jan Vierendeels,et al.  Benchmark exercise for fluid flow simulations in a liquid metal fast reactor fuel assembly , 2016 .

[9]  G. Su,et al.  Study on the coolant mixing phenomenon in a 45° T junction based on the thermal-mechanical coupling method , 2018, Applied Thermal Engineering.

[10]  E. Merzari,et al.  Flow-induced vibration analysis of a helical coil steam generator experiment using large eddy simulation , 2017 .

[11]  Stavros Tavoularis,et al.  Rod bundle vortex networks, gap vortex streets, and gap instability: A nomenclature and some comments on available methodologies , 2011 .

[12]  Leonhard Kleiser,et al.  High-pass filtered eddy-viscosity models for large-eddy simulations of transitional and turbulent flow , 2005 .

[13]  Dirk Lucas,et al.  Computational modelling of flash boiling flows: A literature survey , 2017 .

[14]  G. Su,et al.  Thermal hydraulic and stress coupling analysis for AP1000 Pressurized Thermal Shock (PTS) study under SBLOCA scenario , 2017 .

[15]  Elia Merzari,et al.  Spectral Element Methods for Liquid Metal Reactors Applications , 2017, ArXiv.

[16]  Zhijian Zhang,et al.  CFD simulation for the optimal design and utilization of experiment to research the flow process in PWR , 2016 .

[17]  D. Bestion,et al.  The difficult challenge of a two-phase CFD modelling for all flow regimes , 2014 .

[18]  Simon P. Walker,et al.  CFD simulation of turbulent flow in a rod bundle with spacer grids (MATIS-H) using STAR-CCM+ , 2014 .

[19]  S. V. Möller,et al.  On phenomena of turbulent flow through rod bundles , 1991 .

[20]  E. Merzari,et al.  Anisotropic Turbulence and Coherent Structures in Eccentric Annular Channels , 2009 .

[21]  P. Fischer,et al.  A novel numerical treatment of the near-wall regions in the k−ω class of RANS models , 2018 .

[22]  Nicolás Silin,et al.  Experimental study on the Reynolds number dependence of turbulent mixing in a rod bundle , 2006 .

[23]  Einar M. Rønquist,et al.  An Operator-integration-factor splitting method for time-dependent problems: Application to incompressible fluid flow , 1990 .

[24]  Elia Merzari,et al.  Large-scale large eddy simulation of nuclear reactor flows: Issues and perspectives , 2017 .

[25]  J. Chen,et al.  CFD investigation on thermal-hydraulic behaviors of a wire-wrapped fuel subassembly for sodium-cooled fast reactor , 2018 .

[26]  G. Su,et al.  CFD simulation of secondary side fluid flow and heat transfer of the passive residual heat removal heat exchanger , 2018, Nuclear Engineering and Design.

[27]  S. Tavoularis,et al.  Measurements of turbulent flow in a large-scale model of a 37-rod bundle , 2018, Nuclear Engineering and Design.

[28]  P. Fischer,et al.  Petascale algorithms for reactor hydrodynamics , 2008 .

[29]  Erwin Laure,et al.  Evaluation of Parallel Communication Models in Nekbone, a Nek5000 Mini-Application , 2015, 2015 IEEE International Conference on Cluster Computing.