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
|
# problem from
# PHASE-SPACE FINITE ELEMENT METHODS APPLIED TO THE FIRST-ORDER FORM OF THE TRANSPORT EQUATION
# WILLIAM R. MARTIN,CARL E. YEHNERT, LEONARD LORENCE and JAMES J. DUDERSTAD
# Annals of Nuclear Energy. Vol. 8. pp. 633 to 646, 1981
DEFAULT_ARGUMENT_VALUE 1 quarter
DEFAULT_ARGUMENT_VALUE 2 8
PROBLEM neutron_sn 2D SN $2 ALLOW_UNRESOLVED_BCS MESH $0-$1.msh
c = 1/3
Sigma_t = 0.75
Sigma_s = c*Sigma_t
BC mirror mirror
BC vacuum vacuum
MATERIAL src Sigma_t1=Sigma_t Sigma_s1.1=Sigma_s S1=1
MATERIAL bulk Sigma_t1=Sigma_t Sigma_s1.1=Sigma_s S1=0
sn_alpha = 0.15
SOLVE_PROBLEM
PRINT %.1f 10*phi1(0,1.875)
# PRINT "\# memory= " memory()
WRITE_RESULTS FILE $0-$1-$2.msh
# profile(x) = 10*phi1(x,1.875)
# PRINT_FUNCTION profile MIN 0 MAX 2 NSTEPS 50
|
