McStas logo McStas - A neutron ray-trace simulation package ILL; PSI; ESS Niels Bohr Institute DTU Physics NEXMAP

McStas

About McStas
 Conditions of use
 Authors/Contacts
 Project funding
 Screenshots

Download
 Components
 Linux Install (deb/rpm)
 Mac OS X Install
 Unix Install (src code)
 Windows Install
 Other Downloads (share)

Mailing list

Search web/mailinglist

Documentation
 McStas manual
 FAQ
 Known problems
 Publications
 C Compilers
 Other
 Tools
 Tutorial

Workshops/conferences

Developments

Links

Report bugs

Git

McStas Ubuntu live-dvd


McStas: ESS_butterfly Component

[ Identification | Description | Input parameters | Output parameters | Links ]

The ESS_butterfly Component

ESS butterfly moderator, 2016 revision

Identification

  • Author: Peter Willendrup and Esben Klinkby
  • Origin: DTU
  • Date: August-September 2016
  • Version:(Unknown)

Description

ESS butterfly moderator with automatic choice of coordinate system, with origin
placed at relevant "Moderator Focus Coordinate System" depending on sector location.

To select beamport N 5 simply use 
 
 COMPONENT Source = ESS_butterfly(sector="N",beamline=5,Lmin=0.1,Lmax=20,dist=2,
                                  cold_frac=0.5, yheight=0.03,focus_xw=0.1, focus_yh=0.1)

Geometry
The geometry corresponds correctly to the latest release of the butterfly moderator,
including changes warranted by the ESS CCB in July 2016, the so called BF1 type moderator.
A set of official release documents are available with this component, see the benchmarking
website mentioned below.

Brilliances, geometry adapted from earlier BF2 design
The geometry and brightness data implemented in the McStas ESS source component ESS_butterfly.comp, 
are released as an updated component library for McStas 2.3, as well as a stand alone archive for
use with earlier versions of McStas.

The following features are worth highlighting: 
  • The brightness data are still based on last years MCNP calculations, based on the Butterfly 2 geometry. As a result, the spatial variation of the brightness across the moderator face should be considered to have an uncertainty of the order of 10%. Detailed information on the reasoning behind the change to the Butterfly 1 geometry can be found in [1] and detailed information on horizontal spatial brightness variation can be found in [2]. The spectral shape has been checked and has not changed significantly.
  • A scaling factor has been introduced to in order to account for the decrease in brightness since 2015. To accommodate the influence of the changed geometry, this scaling factor has been applied independently for the cold and thermal contributions and is beamline dependent. It is adjusted to agree with the spectrally-integrated 6cm width data shown in [1],Figure 3.
  • To allow future user adjustments of brilliance, the scalar parameters c_performance and t_performance have been implemented. For now, we recommend to keep these at their default value of 1.0.
  • The geometry has been updated to correspond within about 2 mm to the geometry described in [1]. This has been done by ensuring that the position and apparent width of the moderators correspond to [1],Figure 2, which has been derived from current MCNP butterfly 1 model.
  • The beamport is now defined directly by its sector and number (e.g. 'W' and '5'), rather than giving the angle, as before. [1],Figure 5 shows the geometry of the moderator2, beamport insert and beamline axis for beamline W5. Since the underlying data is still from last years MCNP run, when the brightness was calculated at 10-degree intervals, this means that the spectral curve for the nearest beamport on the grid 5,15,25,35,45,55 degrees is used. The use of this grid has no effect on the accuracy of the geometry or brilliance because of the above- mentioned beamline-dependent adjustments to the brilliance and geometry. See the website [3] for details.
As before, the beamports all originate at the focal point of the sector. The beamline will in almost all cases be horizontally tilted in order to view the cold or thermal moderator, which should be done using an Arm component.

We expect to release an MCNP-event-based source model later in 2016, and possibly also new set of brilliance functions for ESS_butterfly.comp. These are expected to include more realistic brilliances in terms of variation across sectors and potentially also performance losses due to engineering reality. Engineering reality An ad-hoc method for future implementation of "engineering reality" is included, use the "c_performance/t_performance" parameters to down-scale performance uniformly across all wavelengths. References:

  1. Release document "Update to ESS Moderators, latest version"
  2. Release document "Description and performance of the new baseline ESS moderators, latest version"
  3. http://ess_butterfly.mcstas.org/ benchmarking website with comparative McStas-MCNP figures
  4. html-based, interactive 3D model of moderators and monolith, as seen from beamline N4.
  5. Source code for ESS_butterfly.comp at GitHub.

Input parameters

Parameters in boldface are required; the others are optional.
Name Unit Description Default
sector str Defines the 'sector' of your instrument position. Valid values are "N","S","E" and "W" "N"
beamline 1 Defines the 'beamline number' of your instrument position. Valid values are 1..10 or 1..11 depending on sector 1
yheight m Defines the moderator height. Valid values are 0.03 m and 0.06 m 0.03
cold_frac 1 Defines the statistical fraction of events emitted from the cold part of the moderator 0.5
target_index 1 Relative index of component to focus at, e.g. next is +1 this is used to compute 'dist' automatically. 0
dist m Distance from origin to focusing rectangle; at (0,0,dist) - alternatively use target_index 0
focus_xw m Width of focusing rectangle 0
focus_yh m Height of focusing rectangle 0
c_performance 1 Cold brilliance scalar performance multiplicator c_performance > 0 1
t_performance 1 Thermal brilliance scalar performance multiplicator t_performance > 0 1
Lmin AA Minimum wavelength simulated  
Lmax AA Maximum wavelength simulated  
tmax_multiplier 1 Defined maximum emission time at moderator, tmax= tmax_multiplier * ESS_PULSE_DURATION. 3
n_pulses 1 Number of pulses simulated. 0 and 1 creates one pulse. 1
acc_power MW Accelerator power in MW 5
tfocus_dist m Position of time focusing window along z axis 0
tfocus_time s Time position of time focusing window 0
tfocus_width s Time width of time focusing window 0

Output parameters

Name Unit Description Default
cx  
cz  
sign_bl_angle  
orientation_angle  
jmax  
C1_x  
C1_z  
C2_x  
C2_z  
C3_x  
C3_z  
T1_x  
T1_z  
T2_x  
T2_z  
T3_x  
T3_z  
rC1_x  
rC1_z  
rC2_x  
rC2_z  
rC3_x  
rC3_z  
rT1_x  
rT1_z  
rT2_x  
rT2_z  
rT3_x  
rT3_z  
tx  
ty  
tz  
r11  
r12  
r21  
r22  
xf  
yf  
zf  
w_mult  
w_stat  
w_geom  
w_focus  
w_tfocus  
w_geom_c  
w_geom_t  
tx  
ty  
tz  
dt  
lambda  
l_range  
k  
v  
r  
dx  
dy  
dz  
internal_angle  
x0  
z0  
cos_beamport_angle  
sin_beamport_angle  
cos_thermal  
cos_cold  
Lmin_sampled  
Lmax_sampled  

Links


[ Identification | Description | Input parameters | Output parameters | Links ]

Generated automatically by McDoc, Peter Willendrup <peter.willendrup@risoe.dk> / Mon May 29 08:29:49 2017


Last Modified: Monday, 29-May-2017 10:29:49 CEST
Search website mailinglist archive GitHub repos