26.10.2014 Views

JAEA-Data/Code 2007-004 - Welcome to Research Group for ...

JAEA-Data/Code 2007-004 - Welcome to Research Group for ...

JAEA-Data/Code 2007-004 - Welcome to Research Group for ...

SHOW MORE
SHOW LESS

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

EPSI6 = 0.0 (not used)<br />

Card-003-4 Miscellaneous data (6E12.5)<br />

XMIS1 External extrapolated boundary constant C g<br />

> 0 The constant <strong>for</strong> all extrapolated boundaries (See NUAC11-16) <strong>for</strong> all<br />

groups<br />

< 0 This is the <strong>to</strong>tal number of energy groups and Card-003-5 are <strong>to</strong> follow<br />

this card, which give the group-dependent extrapolated boundary<br />

constants <strong>for</strong> problem boundaries <strong>for</strong> all energies.<br />

= 0 The code will use the built-in value <strong>for</strong> all extrapolated boundaries<br />

<br />

Note:<br />

The extrapolated boundary constant C g of energy group g is defined as<br />

C<br />

g<br />

⎡−<br />

D<br />

= ⎢<br />

⎢⎣<br />

ϕ g<br />

g<br />

∂ϕ<br />

g<br />

⎤<br />

⎥<br />

∂x<br />

⎥⎦<br />

X =<br />

X 0<br />

(x) ϕ g<br />

and as the extrapolated length d g is<br />

defined by<br />

∂ϕ g<br />

∂x<br />

x 0<br />

1 ⎡ −1<br />

∂ϕ<br />

g<br />

⎤<br />

≡ ⎢ ⎥<br />

d g ⎢⎣<br />

ϕ g ∂x<br />

⎥⎦<br />

X = X<br />

0<br />

x0<br />

d g<br />

x<br />

thus the relation C g = D g /d g holds, where D g is diffusion coefficient.<br />

As the diffusion coefficient D g is defined as<br />

1 λtr,<br />

g<br />

D g<br />

= = ,<br />

3Σ<br />

3<br />

tr,<br />

g<br />

the constant C g has the relation below;<br />

C = λ / 3d<br />

g<br />

tr,<br />

g<br />

g<br />

If the extrapolation length of the Milne’s problem<br />

d g 0.7105λtr,<br />

g<br />

= is applied,<br />

the extrapolation constant C g becomes the defaulted value (0.4692) which does not<br />

depend on neutron energy. It is the reason <strong>to</strong> apply XMIS1=0.0 (defaulted value <strong>for</strong><br />

all groups). If a sufficiently large value of extrapolation constant is applied,<br />

d g =D g /C g tends <strong>to</strong> zero so that the zero flux at the boundary is assumed.<br />

XMIS2<br />

The extrapolation constant <strong>for</strong> internal black absorber boundary specified by<br />

106

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!