MCNP - z-coordinates of cylindrical mesh >0?

AI Thread Summary
In MCNP6.2, cylindrical mesh tallies require z-coordinates to be greater than zero, leading to errors when attempting to set the origin at z<0. The manual's description can be misleading, as it refers to the cylindrical coordinate system rather than the original reference system. To resolve this, users must ensure that the jmesh values are positive, positioning the origin at the bottom of the mesh. A workaround involves orienting the axis in the negative z direction while selecting intervals in the positive direction. This adjustment allows for effective modeling within the cylindrical mesh framework.
19matthew89
Messages
46
Reaction score
12
TL;DR Summary
z-coordinates of cylindrical mesh tally must be greater than zero?
Hi everyone,

I am using MCNP6.2 and trying to set up a cylindrical coordinate in a reactor channel. The origin as the midplane of the channel.

In my attempt of setting up a cylindrical FMESH with the origin on the z-axis at the bottom of the channel (so z<0) I got this fatal error message: "z-coordinates of cylindrical mesh tally must be greater than zero"!!!
Is it really like that? Do the cylindrical meshes really only work the z>0 halfspace?
Is there a workaround (like a transformation) that could be used? As it says that z-coordinates have to be positive it doesn't even work putting the origin at the top of the channel.

Thanks a lot in advance.
 
Engineering news on Phys.org
EDIT: Problem solved

I am writing here as I solved the issue and it might be useful for someone else.

In the cylindrical mesh, the value of jmesh have to be positive! The description in the manual is misleading as it says:

"Locations of the coarse mesh points in the y direction for rectangular geometry or in the z direction for cylindrical geometry"

so you might expect it's the position of the coarse meh point expressed with respect to the original system of reference but actually no! It refers to the cylindrical system of reference of the mesh itself and since the origin has to be considered at the bottom of the mesh, the z can consequently only be positive.

I hope it helps in case.
 
  • Like
  • Informative
Likes Grelbr42, berkeman and Alex A
If I am not mistaken, you can get the effect you want by choosing the axis to be in the direction you need. So if you point your axis in the -ve z direction, then choose the intervals in the +ve sense of that vector, you get what you need.
 
Hello, I'm currently trying to compare theoretical results with an MCNP simulation. I'm using two discrete sets of data, intensity (probability) and linear attenuation coefficient, both functions of energy, to produce an attenuated energy spectrum after x-rays have passed through a thin layer of lead. I've been running through the calculations and I'm getting a higher average attenuated energy (~74 keV) than initial average energy (~33 keV). My guess is I'm doing something wrong somewhere...
Back
Top