Application of MORPHY program in lead-cooled fast reactor

Jinzhou Li; Tengfei Zhang; Donghao He; Qingquan Pan; Xiaojing Liu

doi:10.11884/HPLPB202436.230357

Journals >High Power Laser and Particle Beams >Volume 36 >Issue 6 >Page 066001 > Article

High Power Laser and Particle Beams
Vol. 36, Issue 6, 066001 (2024)

Application of MORPHY program in lead-cooled fast reactor

Jinzhou Li, Tengfei Zhang, Donghao He^*, Qingquan Pan, and Xiaojing Liu

Author Affiliations

School of Nuclear Science and Engineering, Shanghai Jiao Tong University, Shanghai 200240, China

show less

DOI: 10.11884/HPLPB202436.230357 Cite this Article

Jinzhou Li, Tengfei Zhang, Donghao He, Qingquan Pan, Xiaojing Liu. Application of MORPHY program in lead-cooled fast reactor[J]. High Power Laser and Particle Beams, 2024, 36(6): 066001 Copy Citation Text

show less

Fig. 1. ELSY assembly arrangement

Download full size | View in the Article

Fig. 2. Different types of grid division

Download full size | View in the Article

Fig. 3. 2D full reactor model

Download full size | View in the Article

Fig. 4. Single assembly model of different kinds of assemblies

Download full size | View in the Article

Fig. 5. Normalized flux distribution of 2D model

Download full size | View in the Article

Fig. 6. Normalized flux deviations distribution of 2D model

Download full size | View in the Article

Fig. 7. 3D model with all control rods (CRs) removed

Download full size | View in the Article

Fig. 8. Normalized flux deviations distribution of all CRs removed model

Download full size | View in the Article

Fig. 9. 3D model with CRs inserted halfway

Download full size | View in the Article

Fig. 10. Normalized flux distribution of different insertion modes

Download full size | View in the Article

Fig. 11. Normalized flux deviation distribution of four inserted-halfway models

Download full size | View in the Article

Fig. 12. CR fully-inserted model

Download full size | View in the Article

Fig. 13. Normalized flux distribution of different full insertion modes

Download full size | View in the Article

Fig. 14. The normalized flux deviation distribution of different full insertion modes

Download full size | View in the Article

fuel type	mass density/(g·cm⁻³)	enrichment /%		mass fraction of Pu-U/%
fuel type	mass density/(g·cm⁻³)	²³⁵U	²³⁹Pu	mass fraction of Pu-U/%
inner fuel	11.0049	0.4038	56.87	14.6
middle fuel	11.0092	0.4039	56.87	15.5
outer fuel	11.0236	0.4037	56.87	18.5

Table 1. Material parameters of various fuel assemblies

method	k_eff (MORPHY)	k_eff (Serpent)	deviation/10⁻⁵
no correction	1.02503	1.00790$ \pm $0.00009	+1713
transport correction	1.00392	1.00790$ \pm $0.00009	−398
(n,xn) correction	1.00656	1.00790$ \pm $0.00009	−134

Table 2. The k_eff calculation results of different correction methods were added

Model	$ {S}_{N} $	k_eff		k_eff deviation/$ {10}^{-5} $	MORPHYcalculation time/(CPU·h)
Model	$ {S}_{N} $	MORPHY	Serpent	k_eff deviation/$ {10}^{-5} $	MORPHYcalculation time/(CPU·h)
	$ {S}_{2} $	1.00656	1.00790$ \pm $0.00009	−134	5
Serpent(2D)	$ {S}_{4} $	1.00797	1.00790$ \pm $0.00009	+7	12
	$ {S}_{6} $	1.00813	1.00790$ \pm $0.00009	+23	27
	$ {S}_{2} $	1.00476	1.00790$ \pm $0.00009	−314	5
Serpent(S)	$ {S}_{4} $	1.00531	1.00790$ \pm $0.00009	−259	12
	$ {S}_{6} $	1.00528	1.00790$ \pm $0.00009	−262	27

Table 3. The k_eff calculation results of different S_N expansion orders

model	k_eff calculation results		deviation/10⁻⁵
model	k_eff (MORPHY)	k_eff (Serpent)	deviation/10⁻⁵
all CRs removed	0.99572	0.99594$ \pm $0.00008	−22

Table 4. The k_eff calculation results of 3D model with all CRs removed

model	k_eff calculation results		deviation/10⁻⁵
model	MORPHY	Serpent	deviation/10⁻⁵
all CRs inserted halfway	0.97691	0.97729$ \pm $0.00008	−38
only CR A inserted halfway	0.99507	0.99486$ \pm $0.00007	+21
only CR C inserted halfway	0.99516	0.99514$ \pm $0.00007	+2
CRs A and B inserted halfway	0.99438	0.99454$ \pm $0.00007	−16

Table 5. The k_eff calculation results of 3D full reactor model with CRs inserted halfway

model	reactivity of MORPHY	reactivity of Serpent	deviation/10⁻⁵
all CRs insert halfway	1881	1865	+16
only CR A insert halfway	65	108	−43
only CR C insert halfway	56	80	−24
CRs A and B insert halfway	134	140	−6

Table 6. The calculation results of reactivity change of 3D full reactor with CRs inserted halfway

model	min negative deviation/%	max positive deviation/%	mean absolute deviation/%
only A inserted halfway	−6.57	2.97	1.60
only C inserted halfway	−8.87	3.19	1.70
CRs A and B inserted halfway	−9.73	3.11	1.76

Table 7. Flux deviation of 3D full reactor with CRs inserted halfway

model	k_eff calculation results
model	MORPHY	Serpent	deviation/10⁻⁵
all CRs fully insert	0.96911	0.96888$ \pm $0.00008	+23
only CR A fully insert	0.99431	0.99460$ \pm $0.00007	−29
only CR C fully insert	0.99440	0.99470$ \pm $0.00007	−30
CR A and B fully insert	0.99405	0.99429$ \pm $0.00007	−24

Table 8. The k_eff calculation results of three dimensional full reactor model with CRs fully inserted

model	CR worth calculation results
model	MORPHY/10⁻⁵	Serpent/10⁻⁵	deviation/10⁻⁵
all CRs fully insert	2661	2706	−45
only CR A fully insert	141	134	+7
only CR C fully insert	132	124	+8
CR A and B fully insert	167	165	+2

Table 9. The calculation results of CR worth for CR assembly

model	min negative deviation/%	max positive deviation/%	mean absolute deviation/%
all CRs fully inserted	−7.68	1.62	1.81
only CR A fully inserted	−5.39	2.84	1.62
only CR C fully inserted	−5.68	3.43	1.78
CR A and B fully inserted	−5.85	2.11	1.37

Table 10. The flux calculation deviation of 3D full reactor with CRs fully inserted

Jinzhou Li, Tengfei Zhang, Donghao He, Qingquan Pan, Xiaojing Liu. Application of MORPHY program in lead-cooled fast reactor[J]. High Power Laser and Particle Beams, 2024, 36(6): 066001

Download Citation

Tools

Save the article for my favorites

Paper Information

微信扫一扫：分享

微信扫一扫：分享