نوع مقاله : مقاله پژوهشی
نویسندگان
1
پژوهشکده راکتور و ایمنی هسته‌ای، پژوهشگاه علوم و فنون هسته‌ای، سازمان انرژی اتمی ایران، صندوق پستی: 1339-14155، تهران، ایران
2
پژوهشکده کاربرد پرتوها، پژوهشگاه علوم و فنون هسته‌ای، سازمان انرژی اتمی ایران، صندوق پستی: 3486-11365، تهران، ایران
چکیده
در این پژوهش، یک روش عددی به‌منظور مدل‌سازی جریان جوشش زیر سرد در یک کانال عمودی با استفاده از مدل شار رانشی ارایه شده است. سیستم معادلات غیرخطی با روشی کاملاً ضمنی و با استفاده از روش نیوتن- کریلف مستقل از ژاکوبین (JFNK) حل شده‌اند. به منظور بهبود کارایی روش JFNK و پایداری عددی، پیش شرطساز نیمه‌ضمنی بر پایه فیزیک مسئله پیاده‌سازی شده است. اخیراً پژوهشگران برای حل معادلات تنک، بزرگ و کاهش هزینه محاسبات، به‌طور گستردهای از روش JFNK در تدوین نرم‌افزارهای ترموهیدرولیک هسته‌ای استفاده نموده‌اند. به‌منظور اعتبارسنجی، نتایج به‌دست آمده با داده‌های آزمایشگاهی، نتایج نرم‌افزار RELAP5 و نیز نتایج حاصل از شبیه‌سازی با روش نیمه‌ضمنی مقایسه شده و توافق مناسبی به‌دست آمد. مقایسه نتایج نشان می‌دهد که روش JFNK با پیش شرطساز نیمه‌ضمنی، همگرایی را حداقل به مقدار 50% سرعت بخشیده و هم‌چنین قابلیت پیش‌بینی کسر بخار با میانگین خطای مطلق کم‌تر از 34/7% در گستره قابل ملاحظه‌ای از فشار را داراست.
کلیدواژه‌ها
عنوان مقاله [English]
Simulation of two-phase flow drift flux model in a vertical channel using Newton–Krylov method
نویسندگان [English]
-
H. Esmaili
1
-
H. Kazeminejad
2
-
H. Khalafi
1
-
S.M. Mirvakili
1
1
Nuclear Science and Technology Research Institute (NSTRI), AEOI, P.O. Box 11365-3486, Tehran, Iran
2
Radiation Application Research School, Nuclear Science and Technology Research Institute , ,AEOI, P.O.Box:11365-3486, Tehran Iran
چکیده [English]
In the present work, a numerical method is proposed in order to model the subcooled boiling flow in a vertical channel using the Drift Flux Model. The system of nonlinear equations is solved with the fully-implicit scheme using the Jacobian-free Newton–Krylov (JFNK) method. In order to improve the efficiency of the JFNK method and its numerical stability, a semi-implicit physics based preconditioning (PBP) is implemented. Recently, the JFNK method has been widely used to solve large and sparse system of equations. To validate the proposed method, the results were compared with the experimental data, the results of modeling by using the RELAP5 code, and the available numerical results in the literature. It was found that the results corresponding to the present work have a good agreement with those of the other mentioned methods. Also, it was found that the convergence rate of the JFNK method with the PBP is at least 50% higher than the JFNK method, and the void fraction mean absolute percentage error (MAPE) is less than 7.34% over a wide range of flow and pressure
کلیدواژه‌ها [English]
-
Two-phase flow
-
Drift Flux Model
-
JFNK method
2. M. Ishii, T. Hibiki, Thermo-Fluid Dynamics of Two-Phase Flow, Springer (2006).
3. D. Liles, W. H. Reed, A Semi-implicit Method for Two phase Fluid Dynamics, J. Comput. Phys., 26, 390-407 (1978).
4. Y. G. Lee, G. C. Park, TAPINS: a thermal-hydraulic system code for transient analysis of a fully-passive integral PWR, Nucl. Eng. Technol. 45,172-185 (2013).
5. S. Talebi, H. Kazeminejad, H. Davilu, A numerical technique for analysis of transient two-phase flow in a vertical tube using the drift flux model, Nucl. Eng. Des. 242, 316–322 (2012).
8. C. H. Frepoli, J. Mahaffy, K. Ohkawa, Notes on the implementation of a fully implicit numerical scheme for a two-phase three-field flow model, Nucl. Eng. Des. 225 (2003) 191–217.
9. D.A. Knoll, D.E. Keyes, Jacobian-free Newton-Krylov methods: a survey of approaches and applications, J. Comput. Phys. 193 (2004) 357–397.
11. H. Elman, et al., Block preconditioners based on approximate commutators, SIAM J. Sci. Comput. 27 (2006) 1651-1668.
12. A. Hajizadeh, H. Kazeminejad, S. Talebi, A new numerical method for solution of boiling flow using combination of SIMPLE and Jacobian-free Newton-Krylov algorithms, Progr. Nucl. Energy, 95 (2017) 48-60.
13. L. Zou, H. Zhao, H. Zhang, Numerical implementation, verification and validation of two-phase flow four-equation drift flux model with
14. V. A. Mousseau, Implicitly balanced solution of the two-phase flow equations coupled to nonlinear heat conduction, J. Comp. Phys. 200 (2004) 104–132.
15. H. Elman, et al., A taxonomy and comparison of parallel block multi-level preconditioners for the incompressible Naviere Stokes equations, J. Comput. Phys. 227 (2008) 1790-1808.
16. J. Reisner, et al., An efficient physics-based preconditioner for the fully implicit solution of small-scale thermally driven atmospheric flows, J. Comput. Phys. 189 (2003) 30–44.
17. L. Hu, et al., JFNK method with a physics-based preconditioner for the fully implicit solution of one-dimensional drift-flux model in boiling two-phase flow, App. Therm. Eng. 116 (2017) 610–622.
18. Y. Saad, Iterative Methods for Sparse Linear Systems, second ed., (2003).
19. T. Hibiki, M. Ishii, One-dimensional drift-flux model and constitutive equations for relative motion between phases in various two-phase flow regimes, Int. J. Heat Mass Transf. 46 (2003) 4935–4948.
20. W. Wagner, A. Pruss, The IAPWS formulation 1995 for thermodynamic properties of ordinary water for general and science use, J. Phys. Chem. Ref. Data 31 (2002) 387–535.
23. Y. Saad, M. H. Schultz, GMRES: a generalized minimal residual algorithm for solving non-symetric linear systems, SIAM J. Sci. Stat. Comp. 7 (1986) 856-569.
24. G. G. Bartolomei, V. M. Chanturiya, Experimental study of true void fraction when boiling subcooled water in vertical tubes, Therm. Eng. 14 (1967) 123–128.
25. G. G. Bartolomej, et al., An experimental investigation of the true volumetric vapor content with subcooled boiling tubes, Therm. Eng. 29 (1982) 132–1
)