Design and Development of Neutronics and Thermal Hydraulics Modeling Code for ACP1000 Nuclear Reactor Dynamics in LabVIEW


  • Arshad Habib Malik Department of Basic Training, Chashma Centre of Nuclear Training, Pakistan Atomic Energy Commission, Chashma, Pakistan
  • Feroza Arshad Department of Information System, Karachi Nuclear Power Generating Station, Pakistan Atomic Energy Commission, Karachi, Pakistan
  • Aftab Ahmad Memon Department of Telecommunication Engineering, Mehran University of Engineering and Technology, Jamshoro, Pakistan



Reactor Neutronics, Thermal Hydraulics, Linear Optimization, Nonlinear Optimization, ACP1000, Nuclear Power Plant


An advanced neutronics and thermal hydraulics nuclear code, called GNTHACP code, is designed and developed in LabVIEW as Graphical Neutronics and Thermal Hydraulics toolkit for 1100 MWe Advanced Chinese PWR (ACP-1000) nuclear power plant.  The reactor neutronics model is developed using a nonlinear point reactor kinetics model, while the reactor thermal hydraulics model is developed based on nonlinear fuel and coolant temperature dynamics. The heart of the GNTHACP code is the control rod reactivity model. Control rod reactivity banks are comprised of four power compensation banks G1, G2, N1, N2 and one temperature compensation bank R. The reactivity control configuration of these banks is highly nonlinear, complex and challenging in nature. The control rod reactivity model as a function of G1, G2, N1, N2 and R banks is optimized using two distinct techniques. The control rod reactivity model is optimized using Simplex Linear Programming (SLP) technique under constraints of reactor power as safety limit and control rod speed as maximum speed limit in LabVIEW. The control rod reactivity model is also optimized and investigated using nonlinear Sequential Quadratic Programming (SQP) technique under same constraints in LabVIEW. All the models are integrated and the state-of-the-art virtual instruments (VIs) are designed for cost function optimization, configuring models and calibration of model parameters in LabVIEW. The integrated model as graphical coupled neutronics and thermal hydraulics modeling code is optimized and validated against the Final Safety Analysis Report (FSAR) and different parameters of interest are investigated and proved within design limits as reported with CORCA and CORTH benchmark nuclear codes. The proposed code is stable, highly efficient and accurate as compared to other nuclear codes.


M. Johnson, S. Lucas, and P. Tsvetkov. Modeling of reactor kinetics and dynamics. Idaho National Laboratory, Report Idaho Falls, Idaho 83415, U.S. Department of Energy, Canada (2010).

A. Laureau, D. Heuer, E.M. Lucotte, P.R. Rubiolo, M. Allibert, and M. Aufiero. Transient coupled calculations of the molten salt fast reactor using transient fission matrix approach. Nuclear Engineering and Design 316 (01): 112-124 (2017).

W.K. Lam. Pressurized water reactor simulator. IAEA Workshop Material, Second Edition, IAEA, Vienna (2005).

A.S. Mollah. Education tool for simulation of safety and transient analysis of a pressurized water reactor. International Journal of Integrated Sciences & Technology 03: 01-10 (2018).

S.U.E. Hakim, A. Abimanyu, and Sutanto. Simulator design of Kartini reactor based on LabVIEW. Journal forum Nukir 12 (1): 29-41 (2018).

S.M.A. Ibrahim. Thermal hydraulic simulations of a PWR nuclear power plant. International Journal of Safety and Security 013 (01): 31-520 (2019).

L.A. Macedo, W.M. Torres, G. Sabundjian, D.A. Andrade, A.B. Junior, P.E. Umbehaun, T.N. Conti, R.N. Mesquita, P.H.F. Masotti, and G. Angelo. Development of a LabVIEW web-based simulator for RELAP. International Nuclear Atlantic Conference, Brazil: 01-13 (2011).

J. Park, J. Jang, H. Kim, J. Choe, D. Yun, P. Zhang, A. Cherezov, and D. Lee. Three-Dimensional nodal diffusion code for pressurized water reactor core analysis. Energies 13(23): 6324 (2020).

A.H, Malik, A.A. Memon, and F. Arshad. Advanced multi-modeling of PWR dynamics and deep learning based computational tool in SIMULINK and LabVIEW. Proceedings of Pakistan Academy of Sciences-A: Physical and Computational Sciences 59(01): 71-81 (2022).

A. Rais, D. Siefman, G. Girardin, M. Hursin, and A. Pautz. Methods and models for the coupled neutronics and thermal-hydraulics analysis of the CROCUS reactor at EFPL. Science and Technology of Nuclear Installations 2015: 01-10 (2015).

F.P. Ribeiro, and J. Su. Coupled neutronics-thermal hydraulics analysis in a coolant sub-channel of a PWR using CFD techniques. International Nuclear Atlantic Conference, 22-27 October: 01-25 (2017).

S.M.H. Mousakazemi. Control of a pressurized light-water nuclear reactor two-point kinetics model with performance index-oriented PSO. Nuclear Engineering and Technology 53 (08): 2556-2563 (2021).

A. Tanaka, K. Futahashi, K. Takanabe, and C. Kurimura. Development of a 3-D simulation analysis system for a PWR control rod drive mechanism. International Journal of Pressure Vessels and Piping 85: 655-661 (2008).

C. Wang, C. Yan, J. Wang, C. Tian, and S. Yu. Parametric optimization of steam cycle in PWR nuclear power plant using improved genetic-simplex algorithm. Applied Thermal Engineering 125: 830-845 (2017).

P.T. Boggs, and J.W. Tolle. Sequential quadratic programming for large scale nonlinear optimization. Journal of Computational and Applied Mathematics 124: 123-137 (2000).

M. Torabi, A. Lashkari, S.F. Maoudi, and S. Bagheri. Neutronic analysis of control rod effect on safety parameters in Tehran Research Reactor. Nuclear Engineering and Technology 50(07): 1017-1023 (2018).




How to Cite

Arshad Habib Malik, Feroza Arshad, & Aftab Ahmad Memon. (2023). Design and Development of Neutronics and Thermal Hydraulics Modeling Code for ACP1000 Nuclear Reactor Dynamics in LabVIEW. Proceedings of the Pakistan Academy of Sciences: A. Physical and Computational Sciences, 60(2), 55–62.



Research Articles