Rotating Translating STL Bunny

This example demonstrates the use of time-dependent transformations in the context of STL-objects. Time dependent transformations are also applicable to other Serpent objects.

Input

The main input file:

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% --- Geometry and material information %
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

% --- Box containing bunny:

surf 1 cyl 0.0 0.0 17.0 -0.1 27

cell 1 0 fill 1  -1 % Fill box with mesh based universe
cell 2 0 outside  1 % Outside of the box is outside

% --- Create the background universe (filled with void)

surf 3 inf
cell 4 2 void -3

% --- Create the STL-bunny

solid 2 1 2                        % type, universe, bacground universe
10 4 5 4 3 2                       % search mesh split, depth and sizes
1 1E-5                             % ray test, merge radius
body fuel fuel fuel                % solid name, cell name, material name
file fuel "bunny.stl" 1.00 0 0 0   % body name, file, scaling factor, origin

% --- Composition for bunny:

mat fuel  4.7984E-02 rgb 222 222 222
92235.03c 4.4994E-02
92238.03c 2.4984E-03
92234.03c 4.9184E-04

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% --- Add up and down translations using acceleration %
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

% --- Upwards acceleration of 1e13 cm/s^2 from 0.0 s to 0.9e-6 s

transa u 1 tlim 0 0.9e-6      3 0 0  1e13

% --- Downwards acceleration of 1e13 cm/s^2 from 0.9e-6 s to 2.7e-6 s
%     First for slowing down and then for acceleration.

transa u 1 tlim 0.9e-6 2.7e-6 3 0 0 -1e13

% --- Upwards acceleration of 1e13 cm/s^2 from 2.7e-6 s to 3.6e-6 s
%     For slowing down

transa u 1 tlim 2.7e-6 3.6e-6 3 0 0  1e13

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% --- Add rotations using acceleration %
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

% --- Angular acceleration to the negative mathematical direction
%     around axis defined by vector starting from 0.0 0.0 0.0 with i, j, k
%     components 0.0 0.0 1.0. Magnitude of angular acceleration is
%     1.11111111111111e+14 degrees/s^2 from 0.0 s to 1.8e-6 s
%     (results in rotation of 180 degrees in 1.8 e-6 s)
%     NOTE: The signs are inverted for Serpent 2.1.29 (positive sign corresponds to
%     negative rotation)

transa u 1 tlim 0 1.8e-6      3 rot 0.0 0.0 0.0 0.0 0.0 1.0 1.11111111111111e+14

% --- Angular acceleration to the positive mathematical direction
%     around axis defined by vector starting from 0.0 0.0 0.0 with i, j, k
%     components 0.0 0.0 1.0. Magnitude of angular acceleration is
%     1.11111111111111e+14 degrees/s^2 from 1.8e-6 s to 3.6e-6 s

transa u 1 tlim 1.8e-6 3.6e-6 3 rot 0.0 0.0 0.0 0.0 0.0 1.0 -1.11111111111111e+14

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% --- Run & plot parameters: %
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

% --- Link cross section libraries

set acelib "/home/vvvillehe/XSdata/xsdir0K"

% --- Link dynamic source (this way the geometry will be plotted for each time interval)

set dynsrc "./source" 1

% --- Define time-interval structure to be used with the simulation
%     180 equally wide time-intervals between 0 s and 3.6e-6 s

tme simutime 2 180 0.0 3.6e-6

% --- Simulate a small number of neutrons just to get the plots

set nps 5000 10 simutime

% --- Increase the precursor and neutron buffer sizes for the dynamic simulation

set pbuf 40
set nbuf 40

% --- Create plots for each plot plane

plot 1 500 399
plot 2 500 399
plot 3 500 500

Additional input-files:

• bunny.stl
• The dynamic source files (source.*) need to be created using a separate criticality source simulation.

Running

NOTE: Since the delayed neutron precursors are not rotated and/or translated alongside with the geometry, some of them will fall outside of the bunny during some parts of the simulation and be emitted from the void-part of the geometry instead. This will lead to a number of warning messages such as

***** Tue Jun 20 14:58:13 2017 (seed = 1497959831)
Warning message from function SamplePointDelnu:

Precursor is not in material?

 x = -2.443193E+00, y = -4.201429E+00, z = 9.463751E+00
 Cell name 4

Output

The following animation can be compiled from the 180 plots rendered for each of the three plotting planes:

Resulting movement of the bunny. Compiled from 180 individual frames plotted during the simulation time-intervals.

Notes

• The transformations are updated at time-interval boundaries, which means that the geometry is static for each time-interval.
• The positions of neutrons and precursors are not updated as the transformations are updated.
• In external source simulations (no delayed neutrons, no coupled calculation) the geometry plots are only produced for the first time interval, although the transformations are applied similarly.