• NUCLEAR TECHNIQUES
  • Vol. 47, Issue 7, 070605 (2024)
Lixuan ZHANG1, Guangliang CHEN1,*, Zhaofei TIAN1, Xinli YIN1..., Hao QIAN1 and Dabin SUN2|Show fewer author(s)
Author Affiliations
  • 1College of Nuclear Science and Technology, Harbin Engineering University, Harbin 150001, China
  • 2Wuhan Second Ship Design and Research Institute, Wuhan 430205, China
  • show less
    DOI: 10.11889/j.0253-3219.2024.hjs.47.070605 Cite this Article
    Lixuan ZHANG, Guangliang CHEN, Zhaofei TIAN, Xinli YIN, Hao QIAN, Dabin SUN. CFD computational analysis techniques for core feature component domain based on momentum source and detailed porous media[J]. NUCLEAR TECHNIQUES, 2024, 47(7): 070605 Copy Citation Text show less
    Schematic diagram of rod bundle channel domain identification (a) Layout of mixing vanes, (b) Illustration of the ray method, (c) Axial diagram of the rod bundle channel, (d) Structured mesh for modeling, (e) Results of domain identification for vane zone, (f) Results of domain identification for spacer zone
    Fig. 1. Schematic diagram of rod bundle channel domain identification (a) Layout of mixing vanes, (b) Illustration of the ray method, (c) Axial diagram of the rod bundle channel, (d) Structured mesh for modeling, (e) Results of domain identification for vane zone, (f) Results of domain identification for spacer zone
    General momentum source scheme and pressure distribution in the vane zone (a) Illustration of the momentum source in the windward side of vane, (b) Pressure distribution in the vane zone, (c) Body-fitted mesh, (d) Structured mesh, (e) Nearest neighbor interpolation with inverse distance weighting correction
    Fig. 2. General momentum source scheme and pressure distribution in the vane zone (a) Illustration of the momentum source in the windward side of vane, (b) Pressure distribution in the vane zone, (c) Body-fitted mesh, (d) Structured mesh, (e) Nearest neighbor interpolation with inverse distance weighting correction
    Flowchart of integrated simulation scheme for detailed porous media and general momentum sources
    Fig. 3. Flowchart of integrated simulation scheme for detailed porous media and general momentum sources
    Independence validation for body-fitted mesh and modeling scheme
    Fig. 4. Independence validation for body-fitted mesh and modeling scheme
    Partial results of mesh sensitivity analysis
    Fig. 5. Partial results of mesh sensitivity analysis
    Validation of simulation results of body-fitted mesh (a) z=17 mm, x=6.3 mm, (b) z=17 mm, x=6.3 mm, (c) z=30 mm, x=6.3 mm, (d) z=90 mm, x=6.3 mm
    Fig. 6. Validation of simulation results of body-fitted mesh (a) z=17 mm, x=6.3 mm, (b) z=17 mm, x=6.3 mm, (c) z=30 mm, x=6.3 mm, (d) z=90 mm, x=6.3 mm
    Flow contour of detailed porous media simulation spacer zone (a) Body-fitted mesh calculation results at z=-2.5Dh, (b) Detailed porous media calculation results at z=-2.5Dh, (c) Body-fitted mesh calculation results at x=0.5 P, (d) Detailed porous media calculation results at x=0.5 P
    Fig. 7. Flow contour of detailed porous media simulation spacer zone (a) Body-fitted mesh calculation results at z=-2.5Dh, (b) Detailed porous media calculation results at z=-2.5Dh, (c) Body-fitted mesh calculation results at x=0.5 P, (d) Detailed porous media calculation results at x=0.5 P
    Simulation of mixing vane zone with different schemes (a) Body-fitted mesh calculation results, (b) Calculation results of global momentum source model, (c) Calculation results of local momentum source model
    Fig. 8. Simulation of mixing vane zone with different schemes (a) Body-fitted mesh calculation results, (b) Calculation results of global momentum source model, (c) Calculation results of local momentum source model
    Comparison of cross-flow at the outlet of mixing vanes with different schemes
    Fig. 9. Comparison of cross-flow at the outlet of mixing vanes with different schemes
    Simulation of axial flow distribution under different schemes at 1 m·s-1 operating condition (a) z=30 mm, x=-6.3 mm, (b) z=90 mm, x=-6.3 mm
    Fig. 10. Simulation of axial flow distribution under different schemes at 1 m·s-1 operating condition (a) z=30 mm, x=-6.3 mm, (b) z=90 mm, x=-6.3 mm
    Simulation results for different schemes at at an average axial flow velocity of 7 m·s-1 (a) Axial velocity distribution of body-fitted mesh at x = 0.5Dh, (b) Axial velocity distribution of local momentum source model at x = 0.5Dh, (c) Axial velocity distribution of global momentum source model at x = 0.5Dh, (d) Axial velocity distribution of body-fitted mesh at z = 1Dh, (e) Axial velocity distribution of local momentum source model at z = 1Dh, (f) Axial velocity distribution of global momentum source model at z = 1Dh, (g) Lateral flow distribution of body-fitted mesh, (h) Lateral flow distribution of local momentum source model, (i) Lateral flow distribution of global momentum source model
    Fig. 11. Simulation results for different schemes at at an average axial flow velocity of 7 m·s-1 (a) Axial velocity distribution of body-fitted mesh at x = 0.5Dh, (b) Axial velocity distribution of local momentum source model at x = 0.5Dh, (c) Axial velocity distribution of global momentum source model at x = 0.5Dh, (d) Axial velocity distribution of body-fitted mesh at z = 1Dh, (e) Axial velocity distribution of local momentum source model at z = 1Dh, (f) Axial velocity distribution of global momentum source model at z = 1Dh, (g) Lateral flow distribution of body-fitted mesh, (h) Lateral flow distribution of local momentum source model, (i) Lateral flow distribution of global momentum source model
    (a) The left image shows the simulation results of body-fitted mesh, (b) while the right image shows integrated simulation scheme for detailed porous media and general momentum sources
    Fig. 12. (a) The left image shows the simulation results of body-fitted mesh, (b) while the right image shows integrated simulation scheme for detailed porous media and general momentum sources
    Diagram of turbulent structure and cross-flow distribution at z=0.5Dh (a) Calculation results of body-fitted mesh, (b) Calculation results of global momentum source model
    Fig. 13. Diagram of turbulent structure and cross-flow distribution at z=0.5Dh (a) Calculation results of body-fitted mesh, (b) Calculation results of global momentum source model
    Diagram of temperature distributions of body-fitted mesh at z = 1 Dh (a), global momentum source model at z = 1 Dh (b)
    Fig. 14. Diagram of temperature distributions of body-fitted mesh at z = 1 Dh (a), global momentum source model at z = 1 Dh (b)
    Circumferential nusselt number distribution of central fuel rod
    Fig. 15. Circumferential nusselt number distribution of central fuel rod
    Axial pressure drop curve in rod bundle channels
    Fig. 16. Axial pressure drop curve in rod bundle channels
    方案 Scheme

    贴体网格

    Body-fitted mesh

    多孔介质和全域动量源联合仿真方案

    Detailed porous media and global momentum source

    网格类型 Mesh type

    结构化网格与非结构化网格

    Hybird mesh (structured mesh and unstructured mesh)

    结构化网格Structured mesh
    网格数量 Cell number / 1065.20.6

    域识别耗时

    Domain identification time / s

    22

    提取贴体网格的动量源耗时

    Extracted source terms

    from body-fitted mesh time / s

    150

    CFD 求解耗时

    CFD solving time / s

    917183
    总耗时 Total / s917355
    Table 1. Simulation efficiency of different schemes
    Lixuan ZHANG, Guangliang CHEN, Zhaofei TIAN, Xinli YIN, Hao QIAN, Dabin SUN. CFD computational analysis techniques for core feature component domain based on momentum source and detailed porous media[J]. NUCLEAR TECHNIQUES, 2024, 47(7): 070605
    Download Citation