利用大涡模拟方法和Lagrangian方法研究了腐蚀产物(颗粒)在核反应堆冷却回路管道中的沉积现象(Reynolds数2.5×105). 利用Lagrangian方法跟踪颗粒轨迹, 颗粒运动方程中包含牵引力、升力、浮力和重力. 流体和颗粒是单相耦合, 计算中忽略颗粒--颗粒间碰撞. 流场计算结果与实验数据以及直接数值模拟计算结果吻合良好. 计算结果显示, 小颗粒趋于沉积在管道中心处, 而大颗粒趋于沉积在管道边缘, 这种趋势随颗粒粒径增大而加强. 在近管道底部区域, 颗粒沉积分布密度随颗粒粒径的增大而提高, 小颗粒在底部分布均匀, 而大颗粒在流向速度低速区聚集.
The deposition of corrosion products (particles) in cooling circuits of water-cooled nuclear reactors has been investigated using large eddy simulation and Lagrangian method (Reynolds number 2.5×105). A particle equation of motion including Stokes drag, lift, buoyancy and gravitational forces is used for particle trajectory analysis. The fluid-particle effect is considered and the particle-particle interact is ignored in this work. Results obtained from the fluid field calculation showed good agreement with experimental data and the predictions of direct numerical simulations. The particle size, drag force, shear-induced lift force and gravity all affected the particle deposition process. The small sized particles tended to deposite near the duct center while large sized particles tended to deposite near the duct edge, which become more obvious with increasing particle size. Close to the bottom of the duct, the particle number density increased with particle size increasing, and a high concentration of large particles appeared in the region with flow velocities lower than the mean, while small particles distribute evenly throughout the flow.
[1] Comley G C W. Prog Nucl Energy, 1985; 16: 41
[2] Murray A P. Nucl Technol, 1986; 74: 324
[3] Hirschberg G, Baradlai P, Varga K, Myburg G, Schunk J, Tilky P, Stoddart P. J Nucl Mater, 1999; 265: 273
[4] Aizawa M, Chiba Y, Hosokawa H, Ohsumi K, Uchida S,
Ishizawa N. J Nucl Sci Technol, 2002; 39: 1041
[5] Lin C C. Prog Nucl Energy, 2009; 51: 207
[6] Germano M. J Fluid Mech, 1992; 238: 325
[7] Piomelli U, Liu J. Phys Fluids, 1995; 7: 839
[8] di Mare F, Jones W P. Int J Heat Fluid Flow, 2003; 24: 606
[9] Rhie C M, Chow W L. AIAA J, 1983; 21: 525
[10] Choi H, Moin P. J Comput Phys, 1994; 13: 1
[11] Yao J, Zhao Y L, Hu G L, Fan J R, Cen K F. Aerosol Sci Technol, 2009; 43: 174
[12] Saffman P G. J Fluid Mech, 1965; 22: 385
[13] Saffman P G. J Fluid Mech, 1968; 31: 624
[14] Mei P. Int J Multiphas Flow, 1992; 18: 145
[15] Yao J, Zhang Y, Wang C H, Matsusaka S, Masuda H. Ind Eng Chem Res, 2004; 43: 7181
[16] Yao J, Zhang Y, Wang C H, Liang C Y. AIChE J, 2006; 52: 3775
[17] Armenio V, Fiorotto V. Phys Fluids, 2001; 13: 2437
[18] Gessner F B, Emery A F. J Fluid Eng, 1981; 103: 445
[19] Gessner F B, Po J K, Emery A F. In: Durst F, Launder B E, Schmidt F W, eds., Measurements of Developing Turbulent Flow in a Square Duct. Berlin: Springer Verlag, 1979: 279
[20] Pan Y, Banerjee S. Phys Fluids, 1996; 8: 2733
[21] Pedinotti S, Mariotti G, Banerjee S. Int J Multiphas Flow, 1992; 18: 927
[22] Eaton J K, Fessler J R. Int J Multiphas Flow, 1994; 20: 169
[23] Zhang H F, Ahmadi G. J Fluid Mech, 2000; 406: 55
