MCNP cone source definition

Click For Summary
SUMMARY

The discussion focuses on defining a cone source using the SDEF command in MCNP, specifically addressing the SI, SP, and SB distributions. The user aims to generate particles at equal probabilities within a cone defined by a half angle of 25.8 degrees, utilizing the cosine value of 0.90 in the SI card. The provided SDEF parameters include X=0, Y=0, Z=0, ERG=14.0, and VEC=1 0 0, with a suggestion to sample cosine values between 0.9 and 1 to accurately represent the cone geometry. The user expresses uncertainty regarding the necessity of the -1 0 pairing in the SI card.

PREREQUISITES
  • Familiarity with MCNP (Monte Carlo N-Particle Transport Code)
  • Understanding of SDEF command syntax and parameters
  • Knowledge of particle distribution functions
  • Basic concepts of geometric probability in particle transport simulations
NEXT STEPS
  • Review MCNP documentation on SDEF command usage
  • Explore examples of cone source definitions in MCNP
  • Learn about the implications of SI, SP, and SB distributions in particle simulations
  • Investigate the significance of cosine values in defining angular distributions
USEFUL FOR

Researchers and practitioners in radiation transport simulations, particularly those using MCNP for modeling cone sources and angular distributions in particle physics.

MadGander
Messages
28
Reaction score
1
TL;DR
Need help defining a cone collimated source in MCNP
Hi folks,

I'm attempting to define an SDEF cone source but am getting tripped up in the SI/SP/SB distributions. I need all particles to be generated at angles within the cone, with all angles in the cone having equal probability. I feel like this should be relatively simple to define since the probability is 1 inside the cone and 0 outside. Any help is appreciated.
 
Engineering news on Phys.org
This document is helpful. The plot doesn't make it look like a cone section of an isotropic source, it looks flat, but that is because it is showing a log plot.
Code:
SDEF X=0 Y=0 Z=0 ERG=14.0 VEC= 1 0 0 DIR=D1
SI1 -1 0.90 1
SP1 0 0.0 1
That is with a half angle of 25.8 degrees (the cosine of which is 0.90, in the SI card). It's a clever use of the histogram up to a value feature, sampling the range of cosine values between 0.9 and 1 to make the cone. I don't know why there is a -1 0 pairing, this doesn't seem to be needed. It's probably worth sticking with the full example anyway.
 

Similar threads

  • · Replies 4 ·
Replies
4
Views
3K
  • · Replies 1 ·
Replies
1
Views
2K
  • · Replies 0 ·
Replies
0
Views
1K
Replies
1
Views
2K
  • · Replies 2 ·
Replies
2
Views
3K
  • · Replies 3 ·
Replies
3
Views
3K
  • · Replies 1 ·
Replies
1
Views
3K
  • · Replies 3 ·
Replies
3
Views
4K
Replies
7
Views
3K
Replies
1
Views
3K