MCNP - z-coordinates of cylindrical mesh >0?

In summary, the conversation discusses setting up a cylindrical coordinate in a reactor channel using MCNP6.2. The user encountered a fatal error message when attempting to set up a cylindrical FMESH with the origin on the z-axis at the bottom of the channel. The error message stated that z-coordinates of cylindrical mesh tally must be greater than zero, and the user asks if there is a workaround for this. The issue was eventually solved by adjusting the value of jmesh to be positive in the cylindrical mesh. The conversation ends with a suggestion to choose the axis direction in the desired direction for the desired effect.
  • #1
19matthew89
47
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
  • #2
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
  • #3
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.
 

FAQ: MCNP - z-coordinates of cylindrical mesh >0?

What is MCNP and how does it handle cylindrical coordinates?

MCNP (Monte Carlo N-Particle) is a software package for simulating nuclear processes, including neutron, photon, electron, or coupled neutron/photon/electron transport. It uses a variety of coordinate systems, including cylindrical coordinates, to define the geometry and mesh tallies. In cylindrical coordinates, the z-coordinate typically represents the height or axial position along the cylinder's axis.

Why do the z-coordinates in cylindrical mesh tally need to be greater than zero?

In MCNP, the z-coordinates in a cylindrical mesh tally need to be greater than zero to ensure that the mesh cells are well-defined and non-overlapping. This constraint helps in accurately modeling the physical geometry and ensures that the simulation results are reliable and meaningful.

How can I define a cylindrical mesh tally in MCNP?

To define a cylindrical mesh tally in MCNP, you need to specify the radial (r), angular (theta), and axial (z) coordinates. The z-coordinates must be specified in increasing order and should be greater than zero. You can define the mesh using the "FMESH" card, where you provide the necessary parameters for the cylindrical coordinate system.

What issues might arise if z-coordinates are set to zero or negative in MCNP?

If the z-coordinates are set to zero or negative in MCNP, it can lead to undefined or overlapping mesh cells, resulting in errors or inaccurate simulation results. The software might not properly interpret the geometry, leading to potential issues in the transport calculations and tally results.

How can I troubleshoot problems with cylindrical mesh tallies in MCNP?

To troubleshoot problems with cylindrical mesh tallies in MCNP, ensure that all z-coordinates are greater than zero and correctly ordered. Check the input file for any syntax errors or inconsistencies. Using visualization tools to inspect the mesh geometry can also help identify and resolve issues. Additionally, reviewing the MCNP user manual and consulting with experienced users or forums can provide further guidance.

Back
Top