Large eddy simulation of turbulent flow using the parallel computational fluid dynamics code GASFLOW-MPI

Abstract GASFLOW-MPI is a widely used scalable computational fluid dynamics numerical tool to simulate the fluid turbulence behavior, combustion dynamics, and other related thermal–hydraulic phenomena in nuclear power plant containment. An efficient scalable linear solver for the large-scale pressure equation is one of the key issues to ensure the computational efficiency of GASFLOW-MPI. Several advanced Krylov subspace methods and scalable preconditioning methods are compared and analyzed to improve the computational performance. With the help of the powerful computational capability, the large eddy simulation turbulent model is used to resolve more detailed turbulent behaviors. A backward-facing step flow is performed to study the free shear layer, the recirculation region, and the boundary layer, which is widespread in many scientific and engineering applications. Numerical results are compared with the experimental data in the literature and the direct numerical simulation results by GASFLOW-MPI. Both time-averaged velocity profile and turbulent intensity are well consistent with the experimental data and direct numerical simulation result. Furthermore, the frequency spectrum is presented and a –5/3 energy decay is observed for a wide range of frequencies, satisfying the turbulent energy spectrum theory. Parallel scaling tests are also implemented on the KIT/IKET cluster and a linear scaling is realized for GASFLOW-MPI.

[1]  Y. Saad,et al.  Iterative solution of linear systems in the 20th century , 2000 .

[2]  R. V. Westphal,et al.  Effect of initial conditions on turbulent reattachment downstream of a backward-facing step , 1984 .

[3]  Luca Ammirabile,et al.  Simulation of helium release in the Battelle Model Containment facility using OpenFOAM , 2013 .

[4]  Y. Saad,et al.  GMRES: a generalized minimal residual algorithm for solving nonsymmetric linear systems , 1986 .

[5]  Yang Zhiyin,et al.  Large-eddy simulation: Past, present and the future , 2015 .

[6]  P. Royl,et al.  Three-dimensional all-speed CFD code for safety analysis of nuclear reactor containment: Status of GASFLOW parallelization, model development, validation and application , 2016 .

[7]  Y. Dubief,et al.  On coherent-vortex identification in turbulence , 2000 .

[8]  E.M.J. Komen,et al.  Validation of the CFX4 CFD code for containment thermal-hydraulics , 2008 .

[9]  Yousef Saad,et al.  Iterative methods for sparse linear systems , 2003 .

[10]  U. Piomelli,et al.  Turbulent structures in accelerating boundary layers , 2000 .

[11]  J. P. Magnaud,et al.  The TONUS CFD code for hydrogen risk analysis: Physical models, numerical schemes and validation matrix , 2008 .

[12]  S. W. Hong,et al.  Spray effect on the behavior of hydrogen during severe accidents by a loss-of-coolant in the APR1400 containment , 2006 .

[13]  G. Grötzbach,et al.  Direct numerical and large eddy simulations in nuclear applications , 1999 .

[14]  Tong Wang,et al.  An Adaptive Control Strategy for Proper Mesh Distribution in Large Eddy Simulation , 2010 .

[15]  James J. Sienicki,et al.  Hydrogen Mixing Analyses for a VVER Containment , 2002 .

[16]  Toshio Kobayashi,et al.  Large eddy simulation of backward-facing step flow , 1992 .

[17]  David M. Driver,et al.  Backward-facing step measurements at low Reynolds number, Re(sub h)=5000 , 1994 .

[18]  M. Freitag,et al.  THAI test facility for experimental research on hydrogen and fission product behaviour in light water reactor containments , 2015 .

[20]  M. Hestenes,et al.  Methods of conjugate gradients for solving linear systems , 1952 .

[21]  M. Saunders,et al.  Solution of Sparse Indefinite Systems of Linear Equations , 1975 .

[22]  P. Moin,et al.  Direct numerical simulation of turbulent flow over a backward-facing step , 1997, Journal of Fluid Mechanics.

[23]  P. Royl,et al.  Analysis of steam and hydrogen distributions with PAR mitigation in NPP containments , 2000 .

[24]  Jürgen Eyink,et al.  Computational validation of the EPR™ combustible gas control system , 2012 .

[25]  P. Spalart.  Direct simulation of a turbulent boundary layer up to Rθ = 1410 , 1988, Journal of Fluid Mechanics.

[26]  J. R. Travis,et al.  Numerical analysis of hydrogen risk mitigation measures for support of ITER licensing , 2010 .