blob: e6960fc263b6e3a19b3b099bff488ea0d7b4112c (
plain)
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
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
|
*--------------------------------------------------------
* Procédure p9_Data
* Base de données grappe de combustible, supercellule et coeur CANDU3
* Appel :
* p9_Data :: <<proc>> [>>XXXXX<<]
* avec en entree
* S procAPPEL = procédure d'appel pour les données de sortie
* et en sortie
* XXXXX = données associées à la procédure d'appel
* Auteur : Emmanuel Saint Aubin
* D&Dv5 Ivan Trancart - 28-11-2024 (maj 05-05-2025)
*--------------------------------------------------------
MODULE END: ;
******
STRING procAPPEL ;
:: >>procAPPEL<< ;
***
REAL Pitch RinTF RoutTF RinTC RoutTC RinAs RoutAs RinBs RoutBs Lgrap
RC0s RC1s RC2s RC3s AC0s AC1s AC2s AC3s RinAc RoutAc RinBc
RC0c RC1c RC2c RC3c AC0c AC1c AC2c AC3c Rcyl
r1s r2s r3s r4s r5s r6s r7s r8s r1c r2c r3c r4c r5c r6c r7c r8c ;
REAL fwt dens denss densf Pgrap Bdev eF eK eA alp Pth B1g dB1
Tcalo Tmod Tfuel Pcalo Pmod Peau Dforce Dgaz Dcala Dgaine Dstnlss
Dzrii Dcadm Dunat Dth Ddupic devpitch PLIMC PLIMG VfuelS VfuelC
RoutBc Bref BrG2 BrPL BrHP Pg n0s n1s n2s n3s n0c n1c n2c n3c ;
INTEGER N0s N1s N2s N3s N0c N1c N2c N3c nang nangs nangf ;
***
INTEGER s1 s2 s3 s4 ns4 ns8 ;
STRING MeM MeP MeD ;
******
*** GRAPPES - GÉOMÉTRIE 2D (+ devpitch)
EVALUATE Pitch Lgrap Rcyl devpitch := 28.575 49.53 16.12171 7.0 ;
EVALUATE RinTF RoutTF RinTC RoutTC := 5.168878 5.6032 6.4478 6.587482 ;
EVALUATE r1s r2s r3s := 0.7221626 2.160324 3.600681 ;
EVALUATE RinAs RoutAs RinBs RoutBs := 0.6122 0.654 0.6122 0.654 ;
EVALUATE RC0s RC1s RC2s RC3s AC0s AC1s AC2s AC3s :=
0.0 1.4885 2.8755 4.3305 0.0 0.0 0.261799 0.0 ;
EVALUATE N0s N1s N2s N3s := 1 6 12 18 ;
EVALUATE n0s n1s n2s n3s := 1.0 6.0 12.0 18.0 ;
EVALUATE r1c r2c r3c := 0.865 2.4525 3.7295 ;
EVALUATE RinAc RoutAc RinBc RoutBc := 0.631 0.675 0.536 0.575 ;
EVALUATE RC0c RC1c RC2c RC3c AC0c AC1c AC2c AC3c :=
0.0 1.73 3.075 4.384 0.0 0.0 0.224399 0.0 ;
EVALUATE N0c N1c N2c N3c := 1 7 14 21 ;
EVALUATE n0c n1c n2c n3c := 1.0 7.0 14.0 21.0 ;
EVALUATE VfuelS := 37.0 $Pi_R RinAs RinAs Lgrap * * * * ;
EVALUATE VfuelC := 8.0 RinAc RinAc * * 35.0 RinBc RinBc * * +
$Pi_R Lgrap * * ;
******
*** TRACKING - TRANSPORT 2D
EVALUATE nangs nangf := 7 11 ;
EVALUATE denss densf := 10.0 20.0 ;
******
*** TRACKING - TRANSPORT 3D
EVALUATE dens := 5.0 ;
EVALUATE nang := 6 ;
******
*** GEOMÉTRIE CANDU 3D
EVALUATE MeM MeP MeD := "MCFD" "PRIM" "DUAL" ;
EVALUATE s1 s2 s3 s4 := 1 2 3 4 ;
EVALUATE ns4 ns8 := 4 8 ;
EVALUATE fwt := 19.2 ; ! kg
******
*** FLUX - DIFFUSION
EVALUATE eF := 1.0E-6 ;
EVALUATE eK := 1.0E-4 ;
EVALUATE eA := 1.0E-3 ;
EVALUATE alp := 0.65 ; !<1
******
*** PARAMÈTRES IMPOSÉS
EVALUATE Pth := 1390.0 ; !MW thermique CANDU3
EVALUATE Pg := 615.0 ; !kW thermique A VERIFIER CANDU3
EVALUATE Bdev := 4000.0 ; !MWj/T
******
*** TAUX DE COMBUSTION (GUESS à partir des données officielles)
EVALUATE PLIMC := 6700.0 ; !kW
EVALUATE PLIMG := 860.0 ; !kW
EVALUATE Bref := 7498.708 ; !MWj/T
EVALUATE BrG2 := Pth 15.531 fwt * / 1000. * ; !MWj/T
EVALUATE BrPL := Pth 15.498 fwt * / 1000. * ; !MWj/T
EVALUATE BrHP := Pth 14.244 fwt * / 1000. * ; !MWj/T
******
*** MÉLANGES CONSTANTS
EVALUATE Tcalo Tmod Tfuel := 560.7 345.7 941.3 ;
EVALUATE Pcalo Pmod Peau := 99.3 99.92 0.0 ;
EVALUATE Dforce Dgaz Dcala Dgaine Dstnlss
Dzrii Dcadm Dunat Dth Ddupic :=
6.57 0.0014 6.44 6.44 7.8894
6.4003 1.16 10.4371 10.0 10.4 ;
******
! Sortie
***
IF procAPPEL "GEO2D-STANDARD" = THEN
:: <<r1s>> <<r2s>> <<r3s>>
<<RinTF>> <<RoutTF>> <<RinTC>> <<RoutTC>> <<Rcyl>>
<<RC0s>> <<RC1s>> <<RC2s>> <<RC3s>> <<AC0s>> <<AC1s>>
<<AC2s>> <<AC3s>> <<N0s>> <<N1s>> <<N2s>> <<N3s>>
<<RinAs>> <<RoutAs>> <<RinBs>> <<RoutBs>> <<Pitch>>
<<nangs>> <<denss>> <<nangf>> <<densf>> ;
ELSEIF procAPPEL "GEO2D-CANFLEX" = THEN
:: <<r1c>> <<r2c>> <<r3c>>
<<RinTF>> <<RoutTF>> <<RinTC>> <<RoutTC>> <<Rcyl>>
<<RC0c>> <<RC1c>> <<RC2c>> <<RC3c>> <<AC0c>> <<AC1c>>
<<AC2c>> <<AC3c>> <<N0c>> <<N1c>> <<N2c>> <<N3c>>
<<RinAc>> <<RoutAc>> <<RinBc>> <<RoutBc>> <<Pitch>>
<<nangs>> <<denss>> <<nangf>> <<densf>> ;
ELSEIF procAPPEL "GEOCELL-STANDARD" = THEN
:: <<r1s>> <<r2s>> <<r3s>>
<<RinTF>> <<RoutTF>> <<RinTC>> <<RoutTC>> <<Rcyl>>
<<RC0s>> <<RC1s>> <<RC2s>> <<RC3s>> <<AC0s>> <<AC1s>>
<<AC2s>> <<AC3s>> <<N0s>> <<N1s>> <<N2s>> <<N3s>>
<<RinAs>> <<RoutAs>> <<RinBs>> <<RoutBs>> <<nangs>> <<denss>>
<<nangf>> <<densf>> <<Pitch>> ;
ELSEIF procAPPEL "GEOCELL-CANFLEX" = THEN
:: <<r1c>> <<r2c>> <<r3c>>
<<RinTF>> <<RoutTF>> <<RinTC>> <<RoutTC>> <<Rcyl>>
<<RC0c>> <<RC1c>> <<RC2c>> <<RC3c>> <<AC0c>> <<AC1c>>
<<AC2c>> <<AC3c>> <<N0c>> <<N1c>> <<N2c>> <<N3c>>
<<RinAc>> <<RoutAc>> <<RinBc>> <<RoutBc>> <<nangs>> <<denss>>
<<nangf>> <<densf>> <<Pitch>> ;
ELSEIF procAPPEL "GEOSCELL-STANDARD" = THEN
:: <<r1s>> <<r2s>> <<r3s>>
<<RinTF>> <<RoutTF>> <<RinTC>> <<RoutTC>> <<Rcyl>>
<<RC0s>> <<RC1s>> <<RC2s>> <<RC3s>> <<AC0s>> <<AC1s>>
<<AC2s>> <<AC3s>> <<N0s>> <<N1s>> <<N2s>> <<N3s>>
<<RinAs>> <<RoutAs>> <<RinBs>> <<RoutBs>> <<Lgrap>> <<Pitch>> ;
ELSEIF procAPPEL "GEOSCELL-CANFLEX" = THEN
:: <<r1c>> <<r2c>> <<r3c>>
<<RinTF>> <<RoutTF>> <<RinTC>> <<RoutTC>> <<Rcyl>>
<<RC0c>> <<RC1c>> <<RC2c>> <<RC3c>> <<AC0c>> <<AC1c>>
<<AC2c>> <<AC3c>> <<N0c>> <<N1c>> <<N2c>> <<N3c>>
<<RinAc>> <<RoutAc>> <<RinBc>> <<RoutBc>> <<Lgrap>> <<Pitch>> ;
ELSEIF procAPPEL "MIXS-STANDARD" = THEN
:: <<Tcalo>> <<Tmod>> <<Tfuel>> <<Pcalo>> <<Pmod>> <<Peau>>
<<Dforce>> <<Dgaz>> <<Dcala>> <<Dgaine>> <<Dstnlss>>
<<Dzrii>> <<Dcadm>> <<Dunat>> <<Dth>> <<Ddupic>> <<Pg>>
<<RinAs>> <<RinAs>> <<n0s>> <<n1s>> <<n2s>> <<n3s>> <<Lgrap>> ;
ELSEIF procAPPEL "MIXS-CANFLEX" = THEN
:: <<Tcalo>> <<Tmod>> <<Tfuel>> <<Pcalo>> <<Pmod>> <<Peau>>
<<Dforce>> <<Dgaz>> <<Dcala>> <<Dgaine>> <<Dstnlss>>
<<Dzrii>> <<Dcadm>> <<Dunat>> <<Dth>> <<Ddupic>> <<Pg>>
<<RinAc>> <<RinBc>> <<n0c>> <<n1c>> <<n2c>> <<n3c>> <<Lgrap>> ;
ELSEIF procAPPEL "GEO3D-STANDARD" = THEN
:: <<nang>> <<dens>> <<Pitch>> <<Lgrap>> <<devpitch>>
<<r1s>> <<r2s>> <<r3s>> <<RinTF>> <<RoutTC>> ;
ELSEIF procAPPEL "GEO3D-CANFLEX" = THEN
:: <<nang>> <<dens>> <<Pitch>> <<Lgrap>> <<devpitch>>
<<r1c>> <<r2c>> <<r3c>> <<RinTF>> <<RoutTC>> ;
ELSEIF procAPPEL "EVOL2D" = THEN :: <<Bdev>> ;
ELSEIF procAPPEL "INSTAGE" = THEN :: <<eF>> <<Pth>> <<alp>> ;
ELSEIF procAPPEL "TAVG-MX" = THEN :: <<Bref>> <<PLIMC>> <<PLIMG>> ;
ELSEIF procAPPEL "TAVG" = THEN :: <<Bref>> <<PLIMC>> <<PLIMG>> ;
ELSEIF procAPPEL "TAVG-G2" = THEN :: <<BrG2>> <<PLIMC>> <<PLIMG>> ;
ELSEIF procAPPEL "TAVG-PL" = THEN :: <<BrPL>> <<PLIMC>> <<PLIMG>> ;
ELSEIF procAPPEL "TAVG-HP" = THEN :: <<BrHP>> <<PLIMC>> <<PLIMG>> ;
ELSEIF procAPPEL "RADIAL" = THEN :: <<Bref>> <<PLIMC>> <<PLIMG>> ;
ELSEIF procAPPEL "INSTREFUEL" = THEN :: <<eF>> <<Pth>> <<alp>> ;
ELSEIF procAPPEL "REFUELDIRECT" = THEN :: <<eF>> <<Pth>> <<alp>> ;
ELSEIF procAPPEL "REACTIVITY" = THEN :: <<eF>> <<Pth>> ;
ELSEIF procAPPEL "CALFLX" = THEN :: <<eF>> <<Pth>> ;
ELSEIF procAPPEL "FLUXAX" = THEN :: <<eA>> ;
ELSEIF procAPPEL "EQUILIB" = THEN :: <<eK>> <<eA>> <<eF>> <<Pth>> ;
***
ELSEIF procAPPEL "CANDU6-M1" = THEN :: <<MeM>> <<s1>> <<s1>> ;
ELSEIF procAPPEL "CANDU6-P2" = THEN :: <<MeP>> <<s2>> <<s4>> ;
ELSEIF procAPPEL "CANDU6-M2" = THEN :: <<MeM>> <<s2>> <<s4>> ;
ELSEIF procAPPEL "CANDU6-D2" = THEN :: <<MeD>> <<s2>> <<s4>> ;
ELSEIF procAPPEL "CANDU6-M3" = THEN :: <<MeM>> <<s3>> <<s4>> ;
ENDIF ;
***
END: ;
|
