FLUKA: RESNUCLEi
Previous Index Next
RESNUCLEi
Scores stopping nuclei on a region basis.
WHAT(1)
: type of products to be scored
= 1.0 : spallation products (all inelastic interactions except
those induced by neutrons below the threshold for
multigroup treatment)
= 2.0 : low energy products, i.e. those produced by neutrons
below the threshold for multigroup treatment (provided
the information is available, see Note 1).
= 3.0 : all residual nuclei are scored (if available, see above)
<= 0.0 : resets the default (= 1.0)
Default
= 1.0
WHAT(2)
= logical output unit
> 0.0 : formatted data are written on the WHAT(2)
unit
< 0.0 : unformatted data are written on the |WHAT(2)
| unit
Values of |WHAT(2)
| < 21 should be avoided (with the
exception of WHAT(1)
= +11).
= 0.0 : resets the default = 11.0
Default
= 11.0 (standard output unit)
WHAT(3)
= Maximum atomic number Z of the residual nuclei distribution
Default
: according to the Z of the element(s) of the
material assigned to the scoring region
WHAT(4)
= Maximum M = N - Z - (NMZ)_min of the residual nuclei
distribution (see Notes 2 and 3 below)
Default
: according to the A, Z of the element(s) of the
material assigned to the scoring region.
WHAT(5)
= scoring region number or name
= -1: all regions (see Note 10)
Default
= 1.0
WHAT(6)
= volume of the region in cm**3 (or, more in general, a
normalization factor by which the scoring shall be divided)
Default
= 1.0
SDUM
= any character string identifying the detector
(max. 10 characters)
Notes:
1) Elements or isotopes for which the FLUKA low-energy neutron cross
sections contain information on the production of residual nuclei
are indicated by "Y" in the "RN" column of the Table in 10} where
the components of the neutron cross section library are listed.
The same information can be obtained by setting the printing flag
in the LOW-NEUT option (WHAT(4)
> 0.0). If such data are available
for a given nuclide, the following message is printed on standard
output:
(RESIDUAL NUCLEI INFORMATIONS AVAILABLE)
2) To minimise storage, nuclei are indexed by Z (with Z_min = 1) and
NMZ = N - Z (with (NMZ)_min = -5). The parameter M is defined as
M = NMZ - (NMZ)_min: therefore M_min = 1. The following relations
can also be useful:
N - Z = M + (NMZ)_min N = M + Z + (NMZ)_min
3) In the case of heavy ion projectiles the default NMZ, based
on the region material, is not necessarily sufficient to score all
the residual nuclei, which could include possible ion fragments.
4) In order to achieve reasonable results for residual nuclei
production the new evaporation module must be activated
(it is the default) and heavy fragment evaporation should
also be activated (it is not the default because of the related
large CPU penalty). Coalescence should also be activated (see
option PHYSICS for all these settings). The old evaporation is
still available, mostly because of historical reasons, but it does
not produce meaningful results for residuals. The new evaporation,
available since 1997, is far more sophisticated in this respect,
while differences in emitted particle spectra and multiplicities
are weak.
5) Starting with Fluka2006.3 protons are scored, together with 2-H,
3-H, 3-He, 4-He, at the end of their path, if transported (see
option IONTRANS). This is a change with respect to previous
versions where protons were not scored.
6) All residual nuclei are scored when they have been fully de-excited
down to their ground or isomeric state.
7) Radioactive decay of residual nuclei can be performed by FLUKA in
the same run (see commands RADDECAY, DCYSCORE, DCYTIMES and
IRRPROFIle) or can be done off-line by a user-written code (see for
instance the program USRSUWEV available with the normal FLUKA
distribution). If command IRRPROFI has been issued, RESNUCLEi
results provided by detectors associated to a cooling time index
via DCYSCORE will be expressed in Bq.
8) An example on how to read RESNUCLEi unformatted output is given
below. An explanation of the meaning of the different variables is
given in the comments at the beginning of the program. The program
lists the Z and A of the produced nuclei, followed by the
corresponding amount per unit volume.
9) A more complex program USRSUW, which allows to compute also
standard deviations over several runs, is available with the normal
FLUKA code distribution in directory $FLUPRO/flutil.
A special version of the same program, USRSUWEV, provides in
addition a calculation of induced activity and of its evolution
in time.
10) Setting WHAT(5)
= -1 will provide the sum of the residual nuclei in
all regions, divided by the value set by the user for WHAT(6)
.
PROGRAM RDRESN
PARAMETER ( MXRSNC = 400 )
CHARACTER*10 TIURSN
CHARACTER RUNTIT*80, RUNTIM*32, FILNAM*80
DIMENSION TIURSN(MXRSNC), ITURSN(MXRSNC), NRURSN(MXRSNC),
& VURSNC(MXRSNC), IMRHGH(MXRSNC), IZRHGH(MXRSNC)
DIMENSION RNDATA(MXRSNC,100,260)
WRITE(*,*)' Type the name of the input file:'
READ (*,'(A)') FILNAM
LQ = INDEX(FILNAM,' ') - 1
OPEN (UNIT=1, FILE=FILNAM, STATUS='OLD', FORM='UNFORMATTED')
OPEN (UNIT=2, FILE=FILNAM(1:LQ)//'.txt', STATUS='UNKNOWN')
READ (1) RUNTIT, RUNTIM, WEIPRI, NCASE
WRITE(2,100) RUNTIT, RUNTIM, NCASE, WEIPRI
DO 1 IRN = 1, MXRSNC
READ (1, END=1000) NRN, TIURSN(NRN), ITURSN(NRN),
& NRURSN(NRN), VURSNC(NRN), IMRHGH(NRN), IZRHGH(NRN), K
IF (ABS(ITURSN(NRN)) .LE. 1) THEN
WRITE (2,200) NRN, TIURSN(NRN), NRURSN(NRN),
& VURSNC(NRN), IZRHGH(NRN), IMRHGH(NRN) + K, K + 1
ELSE IF (ABS(ITURSN(NRN)) .LE. 2) THEN
WRITE (2,300) NRN, TIURSN(NRN), NRURSN(NRN),
& VURSNC(NRN), IZRHGH(NRN), IMRHGH(NRN) + K, K + 1
ELSE
WRITE (2,400) NRN, TIURSN(NRN), NRURSN(NRN),
& VURSNC(NRN), IZRHGH(NRN), IMRHGH(NRN) + K, K + 1
END IF
WRITE(2,'(/,A)') ' Z A Residual nuclei'
WRITE(2,'(A,/)') ' per cm**3 per primary'
READ (1) ((RNDATA(NRN,I,J), I=1,IZRHGH(NRN)), J=1,IMRHGH(NRN))
DO 2 I = 1, IZRHGH(NRN)
DO 3 J = 1, IMRHGH(NRN)
IF(RNDATA(NRN,I,J) .GT. 0.)
& WRITE(2,'(2I4,1P, G15.6)') I, J+K+2*I, RNDATA(NRN,I,J)
3 CONTINUE
2 CONTINUE
1 CONTINUE
1000 CONTINUE
100 FORMAT(/,1X,'*****',2X,A80,2X,'*****',/,/,10X,A32,/,/,
& 10X,'Total number of particles followed ',I9,', for a ',
& 'total weight of ',1P,E15.8,/)
200 FORMAT (/,3X,'Res. nuclei n. ',I3,' "',A10,
& '" , "high" energy products, region n. ',I5,
& /,6X,'detector volume: ',1P,E11.4,' cm**3',/
& 6X,'Max. Z: ',I3,', Max. N-Z: ',I3,' Min. N-Z:',I3)
300 FORMAT (/,3X,'Res. nuclei n. ',I3,' "',A10,
& '" , "low" energy products, region n. ',I5,
& /,6X,'detector volume: ',1P,E11.4,' cm**3',/
& 6X,'Max. Z: ',I3,', Max. N-Z: ',I3,' Min. N-Z:',I3)
400 FORMAT (/,3X,'Res. nuclei n. ',I3,' "',A10,
& '" , all products, region n. ',I5,
& /,6X,'detector volume: ',1P,E11.4,' cm**3',/
& 6X,'Max. Z: ',I3,', Max. N-Z: ',I3,' Min. N-Z:',I3)
END
Example:
MATERIAL 26.0 0.0 7.87 11. 0.0 0. IRON
MATERIAL 30.0 0.0 7.133 12. 0.0 0. ZINC
ASSIGNMAT 11.0 6.0 9.0 0.0 ! Four Fe slabs
ASSIGNMAT 12.0 10.0 0.0 0.0 ! Zn vessel
IONTRANS -2.0
PHYSICS 2.0 0.0 0.0 0.0 0.0 0. EVAPORAT
RESNUCLEI 3.0 22.0 0.0 0.0 6.0 0. FirstFe
RESNUCLEI 1.0 23.0 0.0 0.0 10.0 0. Znvessel
Previous Index Next
