توضیح مفصل راجعبه چشمه ها - complete describtion about Sources

چشمه ها در کد mcnp

ساده ترین نوع چشمه که میتوانیم بنویسیم یک چشمه نقطه ای است بدون هیچ گونه ای توزیع - توزیع در موقعیت / توزیع در انرژی - مثال

نکته : نوع ذره را زیاد ملاک قرار نداده ایم و فرض کردیم نوترون باشد بسته به موضوع میتوانید مثال را به سبک خود دست خوش تغییر قرار دهید

SDEF par=1 pos= 0 0 0 erg=1

این یک چشمه همسانگرد با انرژی 1 مگا الکترون ولت در مکان x=0 y=0 z=0 است

نکته : این چشمه نقطه ای انرژی خود را مثل یک خورشید در تمام جهات تابش میکند

اگر در انرژی ها توزیع داشته بایشم مثلا چشمه دو طیف انرژی 1.2 و 3.2 مگا الکتون ولت تابش میکند در همان موقعیت مبدا مختصات 0 و 0 و 0

در این حالت مینویسیم 

SDEF par=1 pos= 0 0 0 erg=d1

SI1 L 1.2 3.2

SP1 0.5 0.5

توجه داشته باشید در این قسمت اول فرض کردم که توزیع انرژی چشمه گسسته باشد یعنی L 

خط بعد از تعریف کلی که با SI1 شروع میشه مربوط به اطلاعات چشمه است که در این بخش گفتم دو چشمه دارم با توضیع گسسته و با انرژی های 1.2 و 3.2 مگاالکترون ولت

خط بعد مربوط به بیان احتمال رخداد هریک از انرژی هاست یعنی طبق مثال چشمه بالا انرژی 1.2 را با احتمال 50% یا همان 0.5 و 3.2 را همچنین با همین احتمال ساطع میکند

نکته : در این مثال هر دو توزیع که ذکر شد توزیع گسسته هستند در ادامه با توزیع های بیشتری روبرو خواهید شد

حال فرض میکنیم که چهار توزیع انرژی گسسته داریم با احتمال های 10% و 20% و 30% و 40% با انرژی های بترتیب 0.5 و 0.3 و 0.1 و 2.5 مگاالکترون ولت - این مثال همان مثال mcnp primer است - که در تصویر زیر میبینیم

برای مشاهده تصویر بزرگتر اینجا را کلیک کنید

click here for larger picture

همانطور که در تصویر میبینید این چهار طیف انرژی از یک نقطه بصورت همسانگرد در تمام جهات با احتمالات ذکر شده منتشر میشوند

حالا میخواهیم مثال دیگه ای رو با طیف انرژی پیوسته تعریف کنیم نکته اینجاست که اینبار باید بجای L از حرف H استفاده کنید که معرف طیف انرژی پیوسته است

در نهایت یک همچین فرمی رو خواهیم داشت :

SDEF pos= 0 0 0 par=1 erg=d1

SI1 H 0.1 0.3 0.5 1.0 2.5

SP1 0.2 0.4 0.3 0.2

یا اینطور تفسیرش کنیم که : 

SDEF pos= 0 0 0 par=1 erg=d1

SI1 H E1 E2 E3 E4 E5

SP1 P1 P2 P3 P4 P5

-----------------------

SP1 0 P2 P3 P4 P5

p1=0

طبق تصویر و انچه در نظرمان بود احتمال رخداد یا گسیل طیف انرژی از 0 تا E1 صفر است

به زبان دیگر :

از 0 تا e1 احتمال 0 است

از e1 تا e2 احتمال p2 است

از e2 تا e3 احتمال p3 است

از e3 تا e4 احتمال p4 است

از e4 تا e5 احتمال p5 است

یعنی : 

نکته در چند مثال بالا وقتی وقتی نوشتیم 

SP1 P1 P2 P3 P4 P5

در واقع این طور هم میتوانستیم بنویسیم که چون حالت پیش فرض بود ساده ترش کردیم

SP1 D P1 P2 P3 P4 P5

البته بجای D حروف C و V که در زیر امده اند را نیز میتوانید قرار دهید : ص 350/898 mcnp4c manual

نکته : برای پارامتر اطلاعات چشمه یا source information یا همون SI که در مثال های بالا دو نوع توزیع L , H را براشون اوردیم در همین راستا میتوانید بجای L یا H از A , S هم بقرار زیر استفاده کنید 

اگه توزیع مکان داشته باشید یعنی در چند نقطه تابش صورت میگیره و در انرژی توزیعی نداشته باشید به این طریق عمل میکنیم :

SDEF par=1 erg=1 pos=d1

SI1 L -10 0 0 10 0 0

SP1 0.75 0.25

در مثال هر دو چشمه نقطه ای روی محور x ها قرار دارند در مختصات 10 و 0 و 0 و  ... با احتمال 25% و دیگری 75%

حالا اگر چشمه شما یه چشمه پیوسته باشه اما نه با توزیع پله ای یا همون هیستوگرامی بلکه توزیع با تابع ریاضی داشته باشه در این شرایط میبایستی مطابق زیر عمل کنید :

SDEF pos= 0 0 0 par=1 erg=d1

SP1 -2 0.5

طیف انرژی که اینجا اوردیم یه طیف ماکسولیه که در 0.5 مگاالکترون ولت پیک داره

همچنین میتوانید از سایر توابع توزیع ریاضی که در زیر امده اند نیز استفاده کنید :

click here for larger picture

برای مشاهده تصویر بزرگتر اینجا کلیک کنید

نکته توجه داشته باشید که : 

برای مشاهده تصویر بزرگتر کلیک کنید

click here for larger picture

refrence : mcnp4c manual page 353/898

اگر دو چشمه داشته باشیم - دو چشمه نقطه ای - که یکی توزیع پیوسته داشته باشد یا همان هیستوگرامی و دیگری توزیع گسسته داشته باشد در این حالت میتوانیم بگوییم که کارت SDEF و بخصوص پارامتر erg ما تابعی است از مکان d2 یا fpos یا function of position- توجه داشته باشید که مکان ما نیز توزیع دارد d1

SDEF par=2 pos=d1 erg=fpos d2

در ادامه مینویسیم :

برای مشاهده تصویر بزرگتر اینجا کلیک کنید

for larger picture click here