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 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...
Back
Top