QUOTE(barvi7 @ 19.7.2015, 12:50)

Можно "надеяться уверенно",

что спектр будет очень теплым (мягше чем на воде обычной).
А вот и вопрос.
Я поступил сейчас тупо и решительно. Взял двухзонную ячейку, в топливную зону впихнул все 12 твэлов (оболочки и пр. не учитывал), во второй зоне D2O, радиус вычисляется из питча.
Фактически при этом не учитывается ещё и эффект гетерогенности.
Питчей взял два - 11 см и 22 см (в два раза больший).
26 групп, БНАБ-78.
На графике группповой спектр нейтронов в первой (топливной) зоне.
Сплошная линия - питч 11 см (то, что предлагается).
Прерывистая - питч 22 см (то есть, удвоенный питч 11*2).
26-ая группа (ноль по энергии) справа.
Видно, что когда питч увеличивается по сравнению с заложенным 11 см, спектр становится
более тепловым.

P.S. В цифрах
(не нормировано, в отличие от графика):
шаг 11 смCODE
KINF=1.8475594
I=1 ROT=0.0250454
I=2 ROT=0.1393075
I=3 ROT=0.3070543
I=4 ROT=0.4629312
I=5 ROT=0.4054835
I=6 ROT=0.3911151
I=7 ROT=0.2998576
I=8 ROT=0.2494684
I=9 ROT=0.2354656
I=10 ROT=0.2164327
I=11 ROT=0.2048320
I=12 ROT=0.1917865
I=13 ROT=0.1760358
I=14 ROT=0.1572866
I=15 ROT=0.1321001
I=16 ROT=0.1081010
I=17 ROT=0.0888489
I=18 ROT=0.0712913
I=19 ROT=0.0649059
I=20 ROT=0.0555877
I=21 ROT=0.0563577
I=22 ROT=0.0735893
I=23 ROT=0.0476844
I=24 ROT=0.0270787
I=25 ROT=0.0101784
I=26 ROT=0.0490287
шаг 22 смCODE
KINF=1.9955770
I=1 ROT=0.0214112
I=2 ROT=0.1193509
I=3 ROT=0.2577135
I=4 ROT=0.3833526
I=5 ROT=0.3135297
I=6 ROT=0.2648007
I=7 ROT=0.1559485
I=8 ROT=0.0972259
I=9 ROT=0.0745310
I=10 ROT=0.0604625
I=11 ROT=0.0551910
I=12 ROT=0.0512230
I=13 ROT=0.0478352
I=14 ROT=0.0437835
I=15 ROT=0.0379979
I=16 ROT=0.0324462
I=17 ROT=0.0280504
I=18 ROT=0.0238280
I=19 ROT=0.0230388
I=20 ROT=0.0209934
I=21 ROT=0.0226356
I=22 ROT=0.0310982
I=23 ROT=0.0212260
I=24 ROT=0.0128963
I=25 ROT=0.0052372
I=26 ROT=0.1042525
P.P.S. Как считал (интуитивно всё должно быть понятно).
CODE
Integer Nf=12,i
Real dt,S,pitch,Sd,mut,Denw,rot,Avogadro,Denf,enr
Data Cell
Driver=Std
Zones=2
Isotopes=4
Rad=1,2
Names=U235,U238,O,D
Den=8:0
EndData
dt=0.65 # диаметр топливной таблетки
Nf=12 # число твэлов
pitch=11 # шаг под ключ ТВС
* Радиусы
* Площадь твэлов
S=Pi*dt*dt/4*Nf
DC.R[1]=sqrt(S/Pi)
* Площадь тяжёлой воды
Sd=sqrt(3)/2*pitch*pitch
DC.R[2]=sqrt((S+Sd)/Pi)
* Плотности
Avogadro=0.6023
* Топливо
Denf=9.5 # плотность топлива
enr=3.67 # обогащение в процентах
mut=Weight(U238)*(1-enr/100)+Weight(U235)*enr/100+Weight(O)*2
rot=Avogadro*Denf/mut
DC.Den[1,U235]=rot*(1-enr/100)
DC.Den[1,U238]=rot*enr/100
DC.Den[1,O]=rot*2
* D2O
Denw=1 # плотность тяжёлой воды
mut=Weight(D)*2+Weight(O)
rot=Avogadro*Denw/mut
DC.Den[2,D]=2*rot
DC.Den[2,O]=rot
Solve
Put Kinf
For i=1 To 26 By 1
rot=flux["i",1]
put i, rot
EndFor
Drop dt,Nf,S,pitch,Sd,mut,Denw, rot,Avogadro,Denf,enr,i
stop