New mesh plotter capabilities

This category replaces the missing input manual
User avatar
Jaakko Leppänen
Site Admin
Posts: 2387
Joined: Thu Mar 18, 2010 10:43 pm
Security question 2: 0
Location: Espoo, Finland
Contact:

New mesh plotter capabilities

Post by Jaakko Leppänen » Fri Feb 03, 2012 4:39 pm

The reaction rate mesh plotter works similar to Serpent 1:

Code: Select all

mesh <or> <nx> <ny> [ <sym> <x0> <x1> <y0> <y1> <z0> <z1> ]
(see description of the input parameters in Chapter 10 in Serpent manual).

In addition to the old thermal flux / fission rate plots, there are other types that can be used by setting the first parameter > 3. The syntax is then slightly different, and the types are:

4 - thermal flux / fission rate plot
5 - collision points
6 - collided weight
7 - gamma heating
8 - detector value

For types 4 - 7 the syntax is:

Code: Select all

mesh <type> <colmap> <or> <nx> <ny> [ <sym> <x0> <x1> <y0> <y1> <z0> <z1> ]
where <colmap> is the color map used for plotting the distribution. Options are:

1 - hot
2 - cold
3 - hot/cold
4 - jet
5 - grayscale

Colormap 3 is to be used with plots that are associated with two parameters (type 4 = the old thermal flux / fission rate plot).

For type 8 the syntax is:

Code: Select all

mesh <type> <colmap> <det> <or> <nx> <ny> [ <sym> <x0> <x1> <y0> <y1> <z0> <z1> ]
where <det> is the detector name. The plot illustrates the distribution of the detector result, taking into account the response function and the integration domain.
- Jaakko

User avatar
Jaakko Leppänen
Site Admin
Posts: 2387
Joined: Thu Mar 18, 2010 10:43 pm
Security question 2: 0
Location: Espoo, Finland
Contact:

Re: New mesh plotter capabilities

Post by Jaakko Leppänen » Fri Mar 02, 2012 4:42 pm

The efficiency of the delta-tracking rejection sampling routine was added as a mesh plot type in version 2.1.3. The type numbers are:

12 - DT efficiency for neutrons
13 - DT efficiency for gammas
- Jaakko

User avatar
Jaakko Leppänen
Site Admin
Posts: 2387
Joined: Thu Mar 18, 2010 10:43 pm
Security question 2: 0
Location: Espoo, Finland
Contact:

Re: New mesh plotter capabilities

Post by Jaakko Leppänen » Wed May 09, 2012 5:33 pm

Density distribution for the multi-physics interface was added in update 2.1.5. The type number is 11. This capability will probably later be replaced by improved geometry plotter routine.
- Jaakko

User avatar
Jaakko Leppänen
Site Admin
Posts: 2387
Joined: Thu Mar 18, 2010 10:43 pm
Security question 2: 0
Location: Espoo, Finland
Contact:

Re: New mesh plotter capabilities

Post by Jaakko Leppänen » Tue Jun 05, 2012 12:33 pm

I updated the Serpent website and added a few new examples in the plot gallery:

http://montecarlo.vtt.fi/gallery/atr.htm
- Jaakko

MDaeubler
Posts: 3
Joined: Tue Jul 16, 2013 5:54 pm
Security question 1: No
Security question 2: 92

Re: New mesh plotter capabilities

Post by MDaeubler » Wed Nov 27, 2013 7:31 pm

Would it be possible to add the temperature distribution to mesh plotter for the multi-physics interface? It would make it easier to check one's input in the ifc files.

Thanks a lot in advance.

User avatar
Jaakko Leppänen
Site Admin
Posts: 2387
Joined: Thu Mar 18, 2010 10:43 pm
Security question 2: 0
Location: Espoo, Finland
Contact:

Re: New mesh plotter capabilities

Post by Jaakko Leppänen » Thu Nov 28, 2013 10:32 am

Try type 10 without the colormap entry:

Code: Select all

mesh 10 <or> <nx> <ny> [ <sym> <x0> <x1> <y0> <y1> <z0> <z1> ]
You should also see the density and temperature distributions read from the interface in the geometry plot (temperature overrides density).
- Jaakko

User avatar
markdehart
Posts: 40
Joined: Fri May 17, 2013 6:10 pm
Security question 1: No
Security question 2: 92
Location: Idaho National Laboratory

Re: New mesh plotter capabilities

Post by markdehart » Thu Jan 09, 2014 1:32 am

I've been playing with the mesh plotting recently for INL's TREAT reactor. I can't seem to figure out how to locate the plot on the axis of interest. I can certainly pick the axis, and the upper and lower coordinates for each of the other bounds. But the specification you give here doesn't seem to work.
mesh <type> <colmap> <or> <nx> <ny> [ <sym> <x0> <x1> <y0> <y1> <z0> <z1> ]
For or=3 (xy plot) I can specify x0 x1 y0 y1 but not z0 z1. I.e.,
mesh 3 1000 1000 0 -150.0 150.0 -150.0 150.0 70 70
generates the error
Input error in parameter "mesh" on line 575 in file "sss2_m8cal_zr_nodet":

Error in mesh boundaries
I thought maybe I just need to specify z0 (since z0 and z1 should be the same). The code does accept this and runs, but the plots appear to be generated at z=0.0. Leaving this value off completely, i.e,
mesh 3 1000 1000 0 -150.0 150.0 -150.0 150.0
will generate the same plot.

This also happens with the new mesh plotter specs.

Running sss2, Version 2.1.14 (June 28, 2013)

(It still generates very nice pictures though.)
Attachments
sss2_m8cal_zr_nodet_mesh1.jpg
sss2_m8cal_zr_nodet_mesh1.jpg (192.38 KiB) Viewed 20344 times

User avatar
Jaakko Leppänen
Site Admin
Posts: 2387
Joined: Thu Mar 18, 2010 10:43 pm
Security question 2: 0
Location: Espoo, Finland
Contact:

Re: New mesh plotter capabilities

Post by Jaakko Leppänen » Thu Jan 09, 2014 12:06 pm

Unlike the geometry plotter, which is 2-dimensional, the mesh plotter is based on a flux tally that needs a finite volume for calculating the integral reaction rate. If you set z0=z1, the volume of integration is reduced to zero, which is why you get the error message. Putting some distance between z0 and z1 should fix the problem, but even so, getting sufficient statistics in a thin slab may require running a lot of neutron histories.
- Jaakko

User avatar
markdehart
Posts: 40
Joined: Fri May 17, 2013 6:10 pm
Security question 1: No
Security question 2: 92
Location: Idaho National Laboratory

Re: New mesh plotter capabilities

Post by markdehart » Thu Jan 09, 2014 6:25 pm

Oooo. Okay, now that I know this everything makes more sense. And when I was omitting the pair, it likely used the min and max ranges, which explains another oddity I was seeing in a plot (large fission energy deposition occurring in an air channel).

Thanks!

User avatar
Jaakko Leppänen
Site Admin
Posts: 2387
Joined: Thu Mar 18, 2010 10:43 pm
Security question 2: 0
Location: Espoo, Finland
Contact:

Re: New mesh plotter capabilities

Post by Jaakko Leppänen » Thu Jan 09, 2014 11:43 pm

Yes, that's right, the limits are set to geometry min and max by default. The mesh plotter works best in 2D geometries. In 3D the distributions are flattened, which sometimes makes interpreting the plots a bit more complicated.
- Jaakko

Post Reply