1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
|
##############################################################################
# #
# Description : PyGan scritp for VVER pin cell simulation with DRAGON5 #
# Author : L. Fede #
# Date : 2023 #
# #
##############################################################################
#
import lifo
import lcm
import cle2000
import numpy as np
from vverGEO import *
from vverMIX import *
from vverTRACK import *
from vverBU import *
from assertS import *
# Parameters of calculus
assembly='V1_S1'
Library='DLIB_J3p1p1'
Multicompo=1
Lifo3=lifo.new()
#################################
# ----------------------- vverGEO
#################################
pyGEOM = vverGEO("GEOM")
#################################
# ----------------------- vverMIX
#################################
pyMIX = vverMIX("LIBRARY",Library)
###################################
# ----------------------- vverTRACK
###################################
[pyUOX_TBH,pyFIG,pyTRACK,pyTF_EXC] = vverTRACK("UOX_TBH","FIG1.ps","TRACK","TF_EXC",pyGEOM)
#################################
# ------------------ BU iteration
#################################
StepList = lcm.new('LCM','burnup_steps')
StepList['ListBU'] = np.array([0.0, 15.0, 30.0, 50.0], dtype='f')
StepList['ListAutop'] = np.array([15.0, 30.0, 50.0], dtype='f')
StepList['ListCompo'] = np.array([0.0 , 30.0, 50.0], dtype='f')
StepList.close() # close without erasing
pyCOMPO = vverBU("COMPO",pyMIX,pyTRACK,pyTF_EXC,StepList,Multicompo)
Keff_conv = pyCOMPO['EDIBU']['MIXTURES'][0]['CALCULATIONS'][0]['K-INFINITY']
print("testVVER: step 1 converged K-effective=", Keff_conv)
Microlib = pyCOMPO['EDIBU']['MIXTURES'][0]['CALCULATIONS'][0]
Microlib.lib()
assertS(Microlib,'K-EFFECTIVE', 0, 1.269919)
Keff_conv = pyCOMPO['EDIBU']['MIXTURES'][0]['CALCULATIONS'][2]['K-INFINITY']
print("testVVER: step 3 converged K-effective=", Keff_conv)
Microlib = pyCOMPO['EDIBU']['MIXTURES'][0]['CALCULATIONS'][2]
Microlib.lib()
assertS(Microlib,'K-EFFECTIVE', 0, 1.271938)
print("test testVVER completed")
|
