RE: [fluka-discuss]: USRYIELD two-way vs. one-way

From: <me_at_marychin.org>
Date: Wed, 12 Mar 2014 23:13:09 -0700 (PDT)

Hi Brandon,

As pointed out earlier, the following are not something we should expect FLUKA
to digest:
- numbers 100.0000.0 and 50.000.0
- USRBDX scoring from angle zero to angle zero in zero bin

In addition to the above which FLAIR flashes in red, I modified as follow, since
you still haven't divulged what you've been trying to score:

$ grep ^USR brandon.inp

USRYIELD 2374. NEUTRON 21. BOD1 BOD2 1.0Nu35
USRYIELD 2.5 1.0E-14 1. 180. 0.0 1006. &
USRYIELD 2374. PROTON 22. BOD1 BOD2 1.0Pr35
USRYIELD 2.5 1E-14 1. 180. 0.0 1006. &
USRYIELD 2374. NEUTRON 23. BOD2 BOD1 1.0Nu35b
USRYIELD 2.5 1.0E-14 1. 0.5 -0.5 1006. &
USRYIELD 2374. PROTON 24. BOD2 BOD1 1.0Pr35b
USRYIELD 2.5 1.0E-14 1. 1.5 0.5 1006. &
USRBDX -1.0 NEUTRON 25. BOD1 BOD2 N1c
USRBDX 2.5 1E-14 1. 1. &
USRBDX -1. NEUTRON 26. BOD2 BOD1 N1f
USRBDX 2.5 1E-14 1. 1. &
USRBDX 9. NEUTRON 27. BOD1 BOD2 N2c
USRBDX 2.5 1E-14 1. 1. &
USRBDX 9. NEUTRON 28. BOD2 BOD1 N2f
USRBDX 2.5 1E-14 1. 1. &

The change of fortran output units is only for the convenience which follows:
Results: there is no symmetry in one-way scoring whether in USRYIELD or USRBDX.
The only symmetry is in USRBDX two-way scoring:

$ ls brandon001_fort.2[1-4] | while read f; do echo -n $f; tail -n 1 $f ; done
brandon001_fort.21 2.9876E-03
brandon001_fort.22 3.2540E-02
brandon001_fort.23 0.0000E+00
brandon001_fort.24 1.8812E+00

$ ls brandon001_fort.2[5-8] | while read f; do awk 'NR==19{print FILENAME,$0}'
$f ; done
brandon001_fort.25 1.1459E-03
brandon001_fort.26 6.3662E-05
brandon001_fort.27 6.0479E-04
brandon001_fort.28 6.0479E-04

There is indeed mis-positioning of spaces in the above .inp excerpt I pasted
inline, which is good to prevent further cut-and-paste at this point ;)

Mis-positioning of spaces, however, was obviously not an issue when I
recommended the three steps for solving your problem.

:) mary



> On 11 March 2014 at 09:31 "REDDELL, BRANDON D. (JSC-EV511)"
> <brandon.d.reddell_at_nasa.gov> wrote:
>
>
> Hi Mary,
>
>
>
> I think perhaps the input CARD items below (in the email) got shifted in the
> email. In a separate email, I sent an attached input file that runs in FLUKA.
> I am trying to score the proton and neutron flux across a plane in the solid.
> I use USRYIELD all the time and it works usually. For the USRYIELD cards,
> there is not angular binning. I am obtaining the differential flux in energy
> of protons or neutrons. The second differential is over charge. The simulation
> works fine and the results look good compared with other codes. The only
> problem is that forward and backward fluxes are exactly the same.
>
>
>
> As for the USRBDX cards, I added those to the input file for comparison to
> see if they give the same for 1-way and 2-way scoring. They do as well. Have
> you had a chance to look at that file?
>
>
>
> Brandon
>
>
>
>
>
> From: me_at_marychin.org [mailto:me_at_marychin.org]
> Sent: Monday, March 10, 2014 11:33 PM
> To: fluka-discuss_at_fluka.org; REDDELL, BRANDON D. (JSC-EV511)
> Subject: RE: [fluka-discuss]: USRYIELD two-way vs. one-way
>
>
>
> Brandon,
>
> Both of your input files appear to contain multiple errors, including
> numerical formats and the request for zero angular binning from omega zero to
> omega zero. What are you trying to score? Have you checked with FLAIR? Did
> FLUKA actually run?
>
> > >
> > On 09 March 2014 at 22:12 "REDDELL, BRANDON D. (JSC-EV511)"
> > <brandon.d.reddell_at_nasa.gov <mailto:brandon.d.reddell_at_nasa.gov> > wrote:
> >
> > My expectation was that switching the upstream and downsteam assignments
> > will give forward or backward USRYIELD distributions. I have a simple slab,
> > a proton beam normally incident and I am looking at the differential flux of
> > neutrons crossing the surface. I have the following USRYIELD CARDS for the
> > forward flux of neutrons and the backward flux of neutrons between two
> > adjacent bodies. The two files are exactly the same? I haven’t seen this
> > before in previous versions of FLUKA, so I am wondering if I am doing
> > something wrong. My source is a beam of protons normally incident on a slab
> > surface. I must be overlooking something very simple here?
> >
> >
> >
> > RPP BLKH -500000.0 500000.0 -500000.0 500000.0 -500000.0 100.0000.0
> >
> > RPP SPCE -50000.0 50000.0 -50000.0 50000.0 -50000.0 50.000.0
> >
> > RPP TARG -5000.0 500.0 -5000.0 5000.0 0.0 14.9
> >
> > XYP P1 12.967766
> >
> > XYP P2 14.820304
> >
> > END
> >
> > BLKH 5 +BLKH -SPCE
> >
> > SPCE 5 +SPCE -TARG
> >
> > BOD1 5 +TARG +P1
> >
> > BOD2 5 +TARG +P2 -P1
> >
> > BOD3 5 +TARG -P2
> >
> > END
> >
> > GEOEND
> >
> > *
> >
> > ASSIGNMA BLCKHOLE BLKH
> >
> > ASSIGNMA VACUUM SPCE
> >
> > ASSIGNMA ALUMINUM BOD1
> >
> > ASSIGNMA ALUMINUM BOD2
> >
> > ASSIGNMA ALUMINUM BOD3
> >
> > * Forward
> >
> > USRYIELD 2174 NEUTRON -50.0 BOD1 BOD2 1.0Nu35
> >
> > USRYIELD 2.5 1.0D-5 50.0 0.5 -0.5 1006&
> >
> > * Backward
> >
> > USRYIELD 2174 NEUTRON -60.0 BOD2 BOD1 1.0Nu35
> >
> > USRYIELD 2.5 1.0D-5 50.0 0.5 -0.5 1006&
> >
> > *
> >
> > From: me_at_marychin.org <mailto:me_at_marychin.org> [mailto:me_at_marychin.org]
> > Sent: Sunday, March 09, 2014 8:46 PM
> > To: fluka-discuss_at_fluka.org <mailto:fluka-discuss_at_fluka.org> ; REDDELL,
> > BRANDON D. (JSC-EV511)
> > Subject: Re: [fluka-discuss]: USRYIELD two-way vs. one-way
> >
> >
> >
> > Hi Brandon,
> >
> > Swapping the upstream and downstream regions generally would not produce
> > the same results. Is there any corresponding symmetry of some sort in your
> > set up?
> >
> > Cheers,
> >
> > mary
> >
> > > > >
> > > On 08 March 2014 at 05:35 "REDDELL, BRANDON D. (JSC-EV511)"
> > > <brandon.d.reddell_at_nasa.gov <mailto:brandon.d.reddell_at_nasa.gov> > wrote:
> > >
> > > Hello,
> > >
> > >
> > >
> > > I recently performed a simulation and used the USRYIELD card to get the
> > > double differential spectrum crossing a boundary defined by reg1 and reg2.
> > > I swapped the upstream and downstream assignments in WHAT(4) and WHAT(5)
> > > expecting different results to get both forward and backward estimates. I
> > > got the same results in both cases. Can someone tell me if this is the
> > > expected result? In the past, I was under the impression that it would
> > > give me forward or backward double differential scoring and not just the
> > > total.
> > >
> > >
> > >
> > > I am aware that I can use USRBDX to do this same problem.
> > >
> > >
> > >
> > > Thanks,
> > >
> > >
> > >
> > > Brandon
> > >
> > > > >
> >
> >
> >
> > >
>
>
>
Received on Thu Mar 13 2014 - 08:24:55 CET

This archive was generated by hypermail 2.3.0 : Thu Mar 13 2014 - 08:24:57 CET