<content> <p>NEUTRON TRANSPORT INPUT TO YIELD CALCULATIONS</p> </content>

06fc258d-d2ad-4c5f-85dd-0ff4b7c5738e 69 Reg Curry 4 5 false false false 0 3 decimal

&[DATE] &[TIME] - &[FILENAME]

NUCLEAR FISSION YIELD CALCULATIONS

TSTART 1 time :

ALPHA AND LAMDA MODEL

DEFINE SOME UNITS

CHOOSE TAMPER MATERIAL; "DU" OR "WC"

Tamper DU :

M.Tamper 231.0 kg * :

CHOOSE CORE MATERIAL: "U" OR "Pu"

CORE Pu :

M.Core 6.3 kg * :

COMPRESSION RATIO: \[REGION[9aa76a8d31f14602b1b72effb33e911a]]\

CR 2.5 :

AMU 1.66 1027 -^* kg * :

1 eV * 1.602 1019 -^* J *

1 MeV * 1.602 1013 -^* J *

Avogadro's Number

N.A 6.022 1023^* 1 mol / *

Ton.TNT 4.2 109^* J * :

KT 1000 Ton.TNT * :

CORE DEFINITION

HEU 235U PERCENT: \[REGION[ee5d80608d224c02b46b1fbc905fb865]]\

HEU% 100 :

HEUMASS M.Core :

M.core HEUMASS HEU% 100 / * : 6.300 kg *

GRAM MOLECULAR WEIGHT OF 235U: \[REGION[030ac39937c445e5971bc1bfa63e3a99]]\

U235.gmw 235 g * :

MolesU235.kg kg U235.gmw / : mol 4.255 1 mol / *

NumberAtomsU235.kg MolesU235.kg N.A * : 2.563 1024^* 1 mol / *

DATA ELEMENT r (gm/cm3) lf (cm) lt (cm) n U235 1.871 101^* 1.689 101^* 3.596 2.637 Pu239 1.560 101^* 1.414 101^* 4.108 3.172 DU 1.895 101^* 0.000 4.342 0.000 WC 1.480 101^* 0.000 3.159 0.000 55 mat :

ρ DATA 22 el g cm 3^/ * : g cm 3^/ 18.710

λ.f DATA 23 el g cm 3^/ * : g cm 3^/ 16.890

λ.t DATA 24 el g cm 3^/ * : g cm 3^/ 3.596

ρ.DU DATA 42 el g cm 3^/ * : g cm 3^/ 18.950

λ.TampDU DATA 44 el g cm 3^/ * : g cm 3^/ 4.342

ρ.WC DATA 52 el g cm 3^/ * : g cm 3^/ 14.800

λ.TampWC DATA 54 el g cm 3^/ * : g cm 3^/ 3.159

ρ.pu DATA 32 el g cm 3^/ * : g cm 3^/ 15.600

λ.puf DATA 33 el g cm 3^/ * : g cm 3^/ 14.140

λ.put DATA 34 el g cm 3^/ * : g cm 3^/ 4.108

ν DATA 25 el : 2.637

ν.pu DATA 35 el : 3.172

NOW DEFINE λs IN TERMS OF NUCLEAR CROSS SECTION AND NUCLEI NUMBER DENSITY

A 235.04 g mol / * :

n ρ N.A * A / : cm 3 -^4.794 1022^*

σ.f 1.235 barn * :

λ.f 1 n σ.f * / : 0.169 m *

σ.el 4.566 barn * :

σ.t σ.f σ.el + :

σ.t barn 5.801

λ.t 1 n σ.t * / : cm 3.596

A.WC 195.85 g mol / * :

n.WC ρ.WC N.A * A.WC / : cm 3 -^4.551 1022^*

σ.elWC 6.857 barn * :

λ.TampWC 1 n.WC σ.elWC * / : cm 3.205

A.DU 238.05 g mol / * :

n.DU ρ.DU N.A * A.DU / : cm 3 -^4.79 1022^*

σ.elDU 4.804 barn * :

λ.TampDU 1 n.DU σ.elDU * / : cm 4.34

σ.put σ.puel σ.puf + :

A.pu 239.05 g mol / * :

σ.puel 4.394 barn * :

σ.puf 1.80 barn * :

ν CORE U ≡ ν ν.pu if :

ρ.core CORE U ≡ ρ ρ.pu if :

ρ.tamper Tamper WC ≡ ρ.WC ρ.DU if :

R.c 3 M.core * 4 π * ρ.core * / 3 nthroot : cm 4.585

R.T 3 M.Tamper * 4 π * ρ.tamper * / R.c 3^+ 3 nthroot : cm 14.433

ρ.core g cm 3^/ 15.600

ρ.tamper g cm 3^/ 18.950

M.core 6.300 kg *

M.Tamper 231.000 kg *

NOW LET THE RADII OF THE CORE AND TAMPER EXPAND BY Δr \[REGION[40f85b2fef844511aa77e4d06a791dd3]]\

maxk 50 :

k 1 maxk range :

Δr 0.01 cm * :

Rexp R.c k 1 - (Δr * + :

Rexp max cm 5.08

RTexp R.T k Δr * + :

RTexp max cm 14.93

ρexp 3 M.core * 4 π * / Rexp (3 -^* vectorize :

nexp CORE U ≡ ρexp N.A * A / ρexp N.A * A.pu / if :

σt CORE U ≡ σ.t σ.put if :

σf CORE U ≡ σ.f σ.puf if :

ρexp 1 el g cm 3^/ 15.600

ρTexp 3 M.Tamper * 4 π * RTexp 3^Rexp 3^- (* / vectorize :

nTexp Tamper WC ≡ ρTexp N.A * A.WC / ρTexp N.A * A.DU / if :

σ.tamper Tamper WC ≡ σ.elWC σ.elDU if :

σ.tamper barn 4.804

ρTexp 1 el g cm 3^/ 18.909

λf 1 σf nexp * / vectorize :

λt 1 σt nexp * / vectorize :

λT 1 σ.tamper nTexp * / vectorize :

Λ λT λt / vectorize :

λT max λT min / 1.097

Rexp Rexp CR 3 nthroot / :

RTexp RTexp CR 3 nthroot / :

λT λT CR / :

λf λf CR / :

λt λt CR / :

α f1 e 23 α * λT λf * / sqrt * Rexp RTexp - (*^vectorize :

α f2 3 ν α - 1 - (* λf λt * / sqrt Rexp * 3 ν α - 1 - (* λf λt * / sqrt Rexp * cot * 1 - Λ 3 α * λT λf * / sqrt Rexp * 1 - (* - RTexp 23 / λT * 3 α * λT λf * / sqrt RTexp * 1 - (* + / vectorize :

α f3 3 ν α - 1 - (* λf λt * / sqrt Rexp * 3 ν α - 1 - (* λf λt * / sqrt Rexp * cot * 1 - Λ 3 α * λT λf * / sqrt Rexp * 1 + (* + RTexp 23 / λT * 3 α * λT λf * / sqrt RTexp * 1 + (* - / vectorize :

α F α f1 α f2 * α f3 - vectorize :

Ans α F α solve 2 el vectorize :

α Ans :

r Rexp :

slope α maxk el α 1 el - r maxk el cm / r 1 el cm / - / : 0.719 -

αfit slope r r 1 el - cm / (* α 1 el + :

alpha α 1 el slope 21 mat :

alpha 0.797 0.719 - 21 mat

x r cm / :

λslope λf maxk el cm / λf 1 el cm / - r maxk el cm / r 1 el cm / - / :

λfit λslope r r 1 el - cm / (* λf 1 el cm / + :

α 0.797 0.791 0.785 0.779 0.773 0.768 0.762 0.756 0.751 91 mat

x 3.379 3.386 3.393 3.401 3.408 3.415 3.423 3.430 3.437 91 mat

αfit 0.797 0.791 0.786 0.781 0.775 0.770 0.765 0.759 0.754 91 mat

λf cm 5.65 5.69 5.73 5.77 5.80 5.84 5.88 5.92 5.96 91 mat

lambda λf 1 el m / λslope 21 mat :

lambda 0.057 5.577 21 mat

plot x α augment :

Fit x αfit augment :

Lambda x λf cm / augment :

λFit x λfit augment :

λf cm / max 7.668

α max 0.797

x max 3.74

λf cm / min 5.655

α min 0.537

x min 3.38

plot Fit 2 1 sys Lambda λFit 2 1 sys

MassCore HEUMASS :

CR 2.50

R.c cm 4.59

CORE Pu

M.core 6.30 kg *

M.Tamper 231.00 kg *

Tamper DU

Δr 1.00 104 -^* m *

A g 235.04 1 mol / *

E.f 180 MeV * :

M.total MassCore M.Tamper + :

M.total 237.30 kg *

N.o 1 cm 3 -^* :

γ 13 / :

E 1 el 1010 -^KT * :

v 1 el 0 m s / * :

r 1 el R.c CR 3 nthroot / :

V.core 1 el 43 / π * r 1 el 3^* : cm 3^161.54 11 mat

r 1 el cm 3.379

LambdaFit λf 1 el m / λslope 21 mat : 0.06 5.58 21 mat

NeutronFissionEnergy 2 MeV * :

m.n 1.67493 1027 -^* kg *

α 1 el 0.80

Nvel 2 NeutronFissionEnergy * m.n / sqrt :

τ 1 el LambdaFit 1 el m * Nvel / : 2.89 109 -^* s * 11 mat

alpha 0.797 0.719 - 21 mat

LambdaFit 0.057 5.577 21 mat

mint 0.1 μs * :

npoints 8000 :

Δt τ 1 el α 1 el / 100 / : μs 3.63 105 -^*

τ 1 el α 1 el / μs 3.629 103 -^*

j 1 npoints range :

t j el mint j Δt * + :

t max μs 0.39

r j el r 1 el Δr j * + :

r 1 el 0.03 m *

r max 0.83 m *

α j el alpha 2 el r j el r 1 el - cm / (* alpha 1 el + :

alpha 1 el 0.797

FITλ j el LambdaFit 2 el r j el r 1 el - (* LambdaFit 1 el m * + :

alpha 2 el 0.72 -

αf 1 el alpha 1 el :

j 2 npoints range r j el r j 1 - el v j 1 - el Δt * + eval : αf j el alpha 2 el r j 1 - el r 1 el - cm / (* alpha 1 el + eval : τ j el LambdaFit 2 el r j 1 - el r 1 el - (* LambdaFit 1 el m * + Nvel / eval : v j el v j 1 - el 4 π * r j 1 - el (2^* γ * E j 1 - el * V.core j 1 - el M.total * / Δt * + eval : V.core j el 43 / π * r j 1 - el (3^* eval : E j el E j 1 - el N.o V.core j 1 - el * E.f * τ j 1 - el / (e αf j 1 - el τ j 1 - el / (t j 1 - el *^* Δt * + eval : 61 line for

αf 1 el 0.7965

r 1 el cm 3.389

v 1 el 0.00

τ 1 el μs 2.891 103 -^*

V.core 1 el cm 3^161.538

E 1 el KT 1.00 1010 -^*

\[REGION[98d7c3041f874c179dd974343704a8c1]]\

CORE Pu

Yield t μs / E KT / augment :

E length 8000.00

Yield max 20.23

t length 8000.00

Ngrowth t μs / αf augment :

αf max 0.7965

αf length 8000.00

t max 3.90 107 -^* s *

t length 8000.00

Ymax E max :

1 time TSTART - 26.38 s *

Yield Ngrowth 2 1 sys

Ymax KT 20.234

CORE Pu

Tamper DU

M.core 6.300 kg *

M.Tamper 231.000 kg *

CR 2.500