MCNP6.2 - Combination of transformations

  • Thread starter Thread starter 19matthew89
  • Start date Start date
  • Tags Tags
    Mcnp6 Rotations
AI Thread Summary
The discussion centers on the complexities of combining multiple transformations in MCNP6.2, particularly regarding the reference frames used for successive rotations. The user seeks clarification on whether the second rotation's reference frame is the original or the one defined by the first transformation. It is concluded that the reference for subsequent transformations is indeed the new origin established by the first transformation. Additionally, the order of transformations is clarified, revealing that the rotational component must be applied before the translation, contrary to the user's initial assumption. The final transformation card is defined to ensure the correct axial segmentation of the tilted channel.
19matthew89
Messages
46
Reaction score
12
TL;DR Summary
It's unclear which reference frame is used for consecutive transformations and so how to combine TR
Hi everyone.

I am struggling understanding how to combine more than one transformations, especially rotations. This stems mainly form the fact that it's unclear to me what reference frame is used to define the transformations angle if two consecutive transformations are applied. If I have a rotation TR1 and then another rotation TR2, will the reference for rotation TR2 be the original reference frame or the tilted one defined by TR1.
I give the example of interest.
I have a channel which is tilted and translated with respect to center (translation doesn't matter, but let's say of (0 0 1)).
I want to have axial segmentation of this tilted channel to have tallies axially.

The model (too complicated to be changed so I'll stick to how it defines the channel) defines the surface of the tilted and translated channel in two steps

First a transformation TR1 of all the surfaces (let's assume phi=30°, theta=10° and no psi, i.e. psi=0°) is defined and applied in the definition of the surfaces.
E.g. TR1 reads
Code:
TR1= 0                  0              0   
         0.85286867     -0.5         0.15038294
         0.49240396     0.8660254    0.08682363
        -0.17364727      0.0       0.98480791

and for instance a plane 11 (originally parallel to XY plane) is defined via
Code:
11 PZ 1 10.0
.

Then another successive transformation for the cell defining the channel is applied with TRCL, and this rotation is about the z axis (which one? the original one?) and of -60° so it reads:
Code:
TRCL= 0.0         0.0         1.0
       0.5       0.8660254     0.0
      -0.8660254    0.5        0.0
        0.0         0.0        1.0
Now my question. In order to have a series of planes orthogonal to the axis of the channel, do I need to define another transformation, let' say TR3, given by the rotational part of TRCL_rot and TR1 (i.e TRCL_rot*TR1), and then apply that to the planes, such as
Code:
11 PZ 3 10.0
.
?

Or the last one is not necessary because the rotation of TRCL will occur around the new z axis defined by TR1 and so the orthogonality of the planes is already guaranteed by the first transformation TR1?

Thanks in advance

P.S. I have a related question about the orders of the transformation. Is the rotational part applied first and then the result being translated, or first translated and then rotated? I ask because in the latter case, if the orginal reference frame is considered, rotation phi of 180° after a translation is equivalent to bringing the item on the other side of the origin.
 
Engineering news on Phys.org
Ok,
I solved it and found out that if a transformation is applied the reference for the following transformation will be the new origin, given by the first transformation.

About the order I solved it and it's actually the opposite of what I was expecting.
For the example given above, indeed a new transformation TR3 has to be defined with: the right translation vector, but teh order of the rotation is TR1*TRCL_rot, with TRCL_rot the rotational part of the TRCL card defined.

So the right transformation card reads:

Code:
TR3        0         0         1 
       0.85944703  0.48860593  0.15038294
       -0.50379801  0.85944704  0.08682363
       -0.08682363 -0.15038295  0.98480791
 
Hello everyone, I am currently working on a burnup calculation for a fuel assembly with repeated geometric structures using MCNP6. I have defined two materials (Material 1 and Material 2) which are actually the same material but located in different positions. However, after running the calculation with the BURN card, I am encountering an issue where all burnup information(power fraction(Initial input is 1,but output file is 0), burnup, mass, etc.) for Material 2 is zero, while Material 1...
Hi everyone, I'm a complete beginner with MCNP and trying to learn how to perform burnup calculations. Right now, I'm feeling a bit lost and not sure where to start. I found the OECD-NEA Burnup Credit Calculational Criticality Benchmark (Phase I-B) and was wondering if anyone has worked through this specific benchmark using MCNP6? If so, would you be willing to share your MCNP input file for it? Seeing an actual working example would be incredibly helpful for my learning. I'd be really...

Similar threads

Back
Top