Learning MCNP4C and MCNPX

Learning MCNP4C and MCNPX

شماره تماس 09358768859
Learning MCNP4C and MCNPX

Learning MCNP4C and MCNPX

شماره تماس 09358768859

Mcnp Visual Editor قسمت اول

روش کار ما اینطور است که : چون با mcnp4c یا mcnpx راحتیم برنامه را در اون پکیج ها مینویسیم و وارسی هندسی را با این پکیج انجام میدهیم اما کار های زیادی را میتونه انجام بده و فوقالاده ست

طبق شکل زیر 

مسیر زیر را دنبال کنید

File ->Open->select your mcnp file

فایل mcnp خودتان را که در نسخه 4c یا x نوشته اید از داخل پوشه های mcnp4c یا mcnpx باز کنید یا open کنید

طبق انچه که در تصویر بالا میبینید میتوانید صفحه مشاهده تان را از yz به xy یا هر نمای دیگری تغییر دهید همچنین توانایی بزرگنمایی و کوچکنمایی را هم دارید 

نکته : اگر فایل را باز کردید و چیزی ندیدید نقاط مشخص شده در تصویر بالا را تغییر دهید چون اگر برنامه خطا یا هشدار داشته باشد در بدو باز شدن پیام دریافت خواهید کرد 





دانلود های ضروری 3

دانلود نسخه Mcnp Visual Editor

دانلود

پس از Extract کردن توسط نرم افزار Winrar یا Winzip با کلیک کردن روی ایکون زیر میتوانید از ان استفاده کنید نیازی به نصب ندارید

دانلود Winrar نسخه فارسی 32bit , 64bit

دانلود Winrar نسخه انگلیسی 32bit , 64bit

تصویر محیط نرم افزار


پاسخ به سوال شما

مشکل مساله شما حل شد

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

همچنین سلول 102 را به این شکل میبایستی تعریف میکردیم

مساله شما خروجی میدهد

+ پلوت

مش بندی هم میشه

شکل هندسی مساله

مساله مشبندی شده

پلوت کانتور مساله

نکته راستای x ثابت است 

امیدوارم منظورم از مشبندی رو تو تصویر زیر بتونید بگیرید

هر مربع 1 سانت در 1 سانته

توجه کنید :

هدف از مشبندی :

در mcnp4c برای شبکه کردن یک محیط یا تکرار یک سلول مسیر دشواری رو میبایستی طی میکردیم و سر اخر با تالی که مینوشتیم جزییاتی که در شبکه اتفاق میافته را بررسی میکردیم

در mcnpx این کار ساده تر شده و با چند خط دستور شما قدرت بالاتری را در اختیار خواهید داشت

این جزییاتی که گفتیم شامل 

شار عبوری از هر سطح از شبکه

تعداد یا جریان عبوری از هر سطح از شبکه

و ....

پیشنهاد میکنم که manual mcnpx قسمت مشبندی را مطالعه کنید


مش بندی در mcnpx قسمت دوم مش بندی استوانه ای

 - مش بندی استوانه ای

مثال :


توضیح :

cora11 برای تعیین شعاع استوانه است

corb11 در واقع محور z یا محور استوانه است

corc11 زاویه صفر تا 360 درجه است که این سه مختصه با هم یک استوانه کامل را میسازند

همان طور که برای مشبندی مستطیلی گفتیم برای این نوع مشبندی هم همان ها صدق میکنند

- برای پلوت کانتور گرفتن یکی از مولفه ها را ثابت گرفتم اینجا من زاویه صفر درجه را فقط در نظر گرفتم و r , Z هرکدام متحول میشوند

وقتی برنامه را در mcnpx نوشتیم یه فایل mdat به ما میدهد + خروجی و... که با gridconv فایل mdat را باید تبدیل کنیم به فایلی مناسب برای برنامه tecplot

من این کار را انجام دادم و برای برنامه بالا خروجی به شکل زیر در امد

نکته :

این مش بندی که انجام دادیم در زاویه 0 درجه و پلوت 2 بعدی rz بوده و برای شار نوشتم که در تصویر متن برنامه نوشته شده track که ان را به flux تغییر دادم شما میتوانید با تغییر کلمات کلیدی وهمچنین تغییر پارامتر های شعاع و زاویه و... شکل های 2 بعدی و 3 بعدی خوبی ترسیم کنید.

در برنامه تغییراتی دادم و این بار pedep را حساب کردم  شکل زیر

شار را نیز حساب کردم شکل زیر

مش بندی در mcnpx قسمت اول مش بندی مستطیلی

در mcnpx از نظر هندسی سه نوع مشبندی داریم  . به بخش 8 از کتاب manual of mcnpx که در بخش دانلود ها وجود دارد مراجعه کنید و شکل بسیط تری از مش ها را خواهید یافت


1 - مستطیلی

مثال :


1    pz -0.005

2    pz  0.005

3    py -0.075

4    py  0.075

5    px -0.075

6    px  0.17

7    pz  -0.105

8    pz  0.035

9    pz  0.03816

10   px  0.825 

11   so  1.5

12   px  0.02  


m1    28063  1

m2    13027  1

m3    22000  1

mode  e

sdef   erg=0.0669 cel=1  pos=0.095 0 0  x=d1 y=d2 z=d3  par=3

si1  0.02  0.17 

sp1  0  1

si2  -0.075  0.075

sp2  0  1

si3  -0.005  0.005

sp3  0  1

+f8:e 5

شروع مش بندی مستطیلی --------------------------- c

tmesh

rmesh11:e traks

 cora11 0 30i 0.9

corb11 -0.005 5i 0.005

corc11 0.035 10i 0.03816

endmd

پایان مش بندی مستطیلی ---------------------------- c

ctme 30

توضیح : cora بیانگر محور x است و دستور مش 11 را دنبال خواهد کرد که مربوط به کرک مسیر الکترون هاست در تالی 8 و این محور را از 0 تا 0.9 به 30 قسمت مساوی تقسیم کرده است

corb بیانگر محور y هاست درست مانند بالا از 0.005- تا 0.005 به 5 قسمت تقسیم بندی شده

corc هم محور z ها را مشخص میکند

نکته : برای پلوت گرفتن بصورت کانتوری با tecplot بهتر است که یکی از این محور ها را ثابت فرض کنید مثلا

 cora11 0 0.9

یعنی قسمت بندیش نکنید
اجرای برنامه :
نکته : من این برنامه را برای 5 دقیقه Run گرفتم شما برای بلای نیم ساعت امتحان کنید
ادامه مطلب ...

معرفی چشمه و جزیات مرتبط با ان

چشمه را با SDEF مشخص میکنیم سپس نوع ذره PAR=1 یا 2 یا 3 که 1 نوترون ها 2 فوتون ها 3 الکترون ها هستند بعد موقعیت چشمه را با POS=x y z مشخص خواهید کرد سر اخر انرژی چشمه را تعیین میکنییم با ERG=1 یا هر انرژی که مطلوب مساله شماست در واحد MEV مینویسید

مثال :

SDEF Par=2 Pos=0 0 -14 Erg=0.14

چشمه بالا یک چشمه فوتونی با انرژی 140 کیلو الکترون ولت و در موقعیت مشخص شده در دستگاه دکارتی است

نکته چشمه بالا یک چشمه همسانگرد است و نقطه ای و مانند یک خورشید طیف فوتونی یکدست 0.14 خود را به همه جهت میتاباند

=========================

توزیع انرژی در چشمه ها :

1 = توزیع گسسته

2 = توزیع پیوسته پله ای یا هیستوگرامی

3 = توزیع پیوسته با تابع ریاضی

=========================

-> 1 - این نوع توزیع انرژی در چشمه ها را به این شکل مینویسیم 

SDEF Par=2 Erg=d1 Pos=0 0 0

SI1 L E1 E2

SP1  P1 P2

درمجموع این نوع چشمه ها از فرم کلی زیر پیروی میکنند

Erg=dn

SIn L E1 E2 ... Em

SPn   P1 P2 ... Pm

همانطور که میبینید چهار نوع انرژی 0.3 0.5 1 2.5 mev بترتیب با احتمال های 20% 10% 30% , 40% مشخص شده است

نکته مجموع احتمالات باید 100% یا یک باشد.

ادامه مطلب ...

اشنایی با کد VMC یا Visual Monte Carlo

کار با این کد ساده است

ابتدا برنامه را از قسمت دانلود ها دانلود کرده و نصب کنید

سپس برنامه را اجرا کنید

در ابتدا این صفحه ظاهر میشود

طبق شکل گزینه load the simulator را انتخاب کنید


بقیه مطلب را در قسمت ادامه مطالب ببینید


ادامه مطلب ...

تمرین 2 - 2 Exercise

شکل های هندسی زیر را ترسیم کنید ؟

Try to draw this examples with mcnp

1 - فرض کنید تمام سطوح کروی باشند

think that all of the surfaces are spherical

 جواب : دانلود کنید با برنامه Textpad یا Notepad میتوانید فایل را باز کنید

answer click here

2 -

یک استوانه است که قطاع بندی شده

this is a cylinder that we cut it

شکل سه بعدی این استوانه

3d view of above example

جواب : دانلود کنید

answer click here

3 -

this is a collimating system for neutrons i wrote it at the first weeks that i start mcnp

this was for bnct from this article

یک کولیماتور فوقالعده تو در تو و مقداری پیچیده برای نوترون ها نوشته بودم میتونین تغییرش بدین

این نوع کولیماتور ها رو سعی کردم از روی کولیماتور مقاله مربوط به BNCT برنامشو بنویسم

شکل دو بعدی

شکل سه بعدی همین تمرین بصورت Wire Frame

3d view of above example in wire frame mode


شکل سه بعدی بصورت Solid

in solid mod

جواب : دانلود کنید

answer click here

4 - این شکل را ترسیم کنید ؟

try to write it

شکل سه بعدی این تمرین (2 مخروط تو در تو و 2 استوانه تو در تو)

3d view of above example

جواب : دانلود کنید

answer click here

5 - دو استوانه تو در تو ترسیم کنید

draw 2 cylinder one in another one

شکل سه بعدی

3d view of above example

جواب : دانلود کنید

answer click here

ادامه مطلب ...

چگونه یک مکعب رسم کنیم

برای رسم یک مکعب میبایستی 6 وجه تعریف کنیم

مبدا مختصات را مرکز مکعب فرض میکنیم و سپس شروع به شماره گزاری سصوح مینماییم

سطح شماره 1

نکته : برای پرهیز از سردرگمی سطوح را با اعداد 1 تا 100 نامگذاری میکنم و سلول ها را از شماره 101 شروع میکنم

سطح 1 : فرض میکنیم که در صفحه ZY باشد و محور x را در 2- قطع کند

مینویسیم

سطح 2 : فرض میکنیم که در همان صفحه ZY قرار دارد و این بار محور x را در 2 قطع کند

سطح 3 : فرض میکنیم که در صفحه zx قرار دارد و محور y را در 2- قطع کرده

سطح 4 : فرض میکنیم که در همان صفحه zx قرار گرفته و محور y را در 2 قطع کرده است

سطح 5 : در صفحه xy قرار میگیرد و محور z را در 2- قطع میکند

سطح 6 : درست در همان صفحه قرار دارد و محور z را در 2 قطع میکند

حالا نوبت نوشتن سلول ها ست

ما در کل در این مساله دو سلول بیشتر نداریم یکی خود مکعب و هر چه داخلش هست و دیگری فضای بیرون مکعب هر چه هست

نکته : جلوتر برای جلو گیری از ثبت اندرکنش های ضاعد سعی میکنیم که مکعب را در یک کره فرض کنیم حالا با این کار یک سطح و سلول جدید اضافه میشه ولی داده ها ی بهتری در زمان کمتر بدست میاریم

ادامه را در قسمت ادامه مطالب مشاهده کنید

ادامه مطلب ...