MCNP6.2 - Combination of transformations

In summary, the conversation discusses difficulties with combining multiple transformations, particularly rotations. There is confusion about which reference frame is used for the rotation angle when applying consecutive transformations. An example involving a tilted channel and successive transformations is given to clarify the issue. The correct order of transformations is also mentioned, with the rotational part being applied first and then the translation. Ultimately, it is determined that the reference for a transformation will be the new origin created by the previous transformation.
  • #1
19matthew89
47
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
  • #2
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
 
  • Like
Likes Alex A

What is MCNP6.2 and how is it used in scientific research?

MCNP6.2 is a Monte Carlo N-Particle transport code used in nuclear engineering and radiation physics research. It is used to simulate the transport of particles such as neutrons, photons, and electrons through various materials and geometries, allowing scientists to study their interactions and effects.

What are "transformations" in MCNP6.2 and why are they important?

In MCNP6.2, transformations refer to the mathematical operations used to transform a particle's position and direction as it moves through a simulation. These transformations are crucial in accurately modeling particle interactions and determining their behavior in different materials and geometries.

How does MCNP6.2 handle combinations of transformations?

MCNP6.2 uses a combination of rotations, translations, and reflections to transform particles as they move through a simulation. These transformations are applied sequentially, with each one building upon the previous one to accurately track the particle's position and direction.

What are some common applications of MCNP6.2 and its combination of transformations?

MCNP6.2 is commonly used in nuclear reactor design, radiation shielding analysis, and medical physics research. Its combination of transformations allows scientists to model complex systems and study the effects of radiation on different materials and environments.

Are there any limitations to MCNP6.2 and its use of transformations?

While MCNP6.2 is a powerful tool for simulating particle transport, it does have some limitations. For example, it may not accurately model particles with very high energies or extremely small wavelengths. Additionally, the accuracy of the simulation depends on the accuracy of the input data and the assumptions made by the user.

Similar threads

Replies
4
Views
2K
  • Nuclear Engineering
Replies
1
Views
1K
Replies
1
Views
1K
  • Linear and Abstract Algebra
Replies
8
Views
2K
Replies
4
Views
730
  • Linear and Abstract Algebra
Replies
5
Views
2K
Replies
5
Views
1K
  • Classical Physics
Replies
3
Views
2K
  • General Math
Replies
2
Views
4K
  • Special and General Relativity
Replies
7
Views
1K
Back
Top