اشنایی با کد مونت کارلو SIMIND

The SIMIND Monte Carlo program

دانلود فایل اصلی

دانلود فانتوم ها البته در خود فایل اصلی هست

دانلود اموزش

انواع کد های مونت کارلو و مقایسه انها کار تحقیقاتی مفیدیه

برای دریافت اطلاعات بیشتر به سایت زیر سربزنید

http://www.radfys.lu.se/simind/index.asp

نحوه نصب این کد :

یک دایرکتوری در درایو C درست کنید به این نام 

c:\simind

سپس فایل اصلی را در این پوشه Extract کنید با Winrar یا Winzip

وارد Control Panel شده و طبق شکل روی System کلیک کنید

تب Advanced را انتخاب کنید

Envirement Variable را انتخاب کنید

New را انتخاب کنید

سپس Ok کنید

Ok کنید

سپس میتوانید با این کد کار کنید

نحوه مقایسه دو پلوت یا تعداد بیشتری از پلوت ها با Origin 8

در پست اول از Origin گفتم که چطور یک ترسیم را از یک سری بسیار بلند بالا از داده ها ترسیم کنید

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

من برای مثال 4 فایل اماده دارم که میخوام اونها رو پلوت کنم

محیط فعلی برنامه

پس از کامنت گذاری خواهیم داشت

پلوت خروجی های بدست امده از mcnp4c یا mcnpx در Oprigin 8

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

حالت مورد نظر ما خروجی ها باید در نهایت این شکلی شوند

نکته : خروجی ها ستونی و در ردیف های منظم باید باشند

سپس نرم افزار Origin را باز کنید

نکته : فایل txt اماده شده را درگ کنید در پنجره Origin 8 سپس مراحل زیر را دنبال کنید

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 گرفتم شما برای بلای نیم ساعت امتحان کنید