ش | ی | د | س | چ | پ | ج |
1 | 2 | 3 | 4 | 5 | ||
6 | 7 | 8 | 9 | 10 | 11 | 12 |
13 | 14 | 15 | 16 | 17 | 18 | 19 |
20 | 21 | 22 | 23 | 24 | 25 | 26 |
27 | 28 | 29 | 30 | 31 |
روش کار ما اینطور است که : چون با mcnp4c یا mcnpx راحتیم برنامه را در اون پکیج ها مینویسیم و وارسی هندسی را با این پکیج انجام میدهیم اما کار های زیادی را میتونه انجام بده و فوقالاده ست
طبق شکل زیر
مسیر زیر را دنبال کنید
File ->Open->select your mcnp file
فایل mcnp خودتان را که در نسخه 4c یا x نوشته اید از داخل پوشه های mcnp4c یا mcnpx باز کنید یا open کنید
طبق انچه که در تصویر بالا میبینید میتوانید صفحه مشاهده تان را از yz به xy یا هر نمای دیگری تغییر دهید همچنین توانایی بزرگنمایی و کوچکنمایی را هم دارید
نکته : اگر فایل را باز کردید و چیزی ندیدید نقاط مشخص شده در تصویر بالا را تغییر دهید چون اگر برنامه خطا یا هشدار داشته باشد در بدو باز شدن پیام دریافت خواهید کرد
دانلود نسخه Mcnp Visual Editor
پس از Extract کردن توسط نرم افزار Winrar یا Winzip با کلیک کردن روی ایکون زیر میتوانید از ان استفاده کنید نیازی به نصب ندارید
دانلود Winrar نسخه فارسی 32bit , 64bit
دانلود Winrar نسخه انگلیسی 32bit , 64bit
تصویر محیط نرم افزار
مشکل مساله شما حل شد
نکته : شعاع کره را باید از 150 بیشتر میگرفتید چون کم بود چند سلول بینام درست میشد و اون خطای کذایی رو میداد
همچنین سلول 102 را به این شکل میبایستی تعریف میکردیم
مساله شما خروجی میدهد
+ پلوت
مش بندی هم میشه
شکل هندسی مساله
مساله مشبندی شده
پلوت کانتور مساله
نکته راستای x ثابت است
امیدوارم منظورم از مشبندی رو تو تصویر زیر بتونید بگیرید
هر مربع 1 سانت در 1 سانته
توجه کنید :
هدف از مشبندی :
در mcnp4c برای شبکه کردن یک محیط یا تکرار یک سلول مسیر دشواری رو میبایستی طی میکردیم و سر اخر با تالی که مینوشتیم جزییاتی که در شبکه اتفاق میافته را بررسی میکردیم
در mcnpx این کار ساده تر شده و با چند خط دستور شما قدرت بالاتری را در اختیار خواهید داشت
این جزییاتی که گفتیم شامل
شار عبوری از هر سطح از شبکه
تعداد یا جریان عبوری از هر سطح از شبکه
و ....
پیشنهاد میکنم که manual mcnpx قسمت مشبندی را مطالعه کنید
- مش بندی استوانه ای
مثال :
توضیح :
cora11 برای تعیین شعاع استوانه است
corb11 در واقع محور z یا محور استوانه است
corc11 زاویه صفر تا 360 درجه است که این سه مختصه با هم یک استوانه کامل را میسازند
همان طور که برای مشبندی مستطیلی گفتیم برای این نوع مشبندی هم همان ها صدق میکنند
- برای پلوت کانتور گرفتن یکی از مولفه ها را ثابت گرفتم اینجا من زاویه صفر درجه را فقط در نظر گرفتم و r , Z هرکدام متحول میشوند
وقتی برنامه را در mcnpx نوشتیم یه فایل mdat به ما میدهد + خروجی و... که با gridconv فایل mdat را باید تبدیل کنیم به فایلی مناسب برای برنامه tecplot
من این کار را انجام دادم و برای برنامه بالا خروجی به شکل زیر در امد
نکته :
این مش بندی که انجام دادیم در زاویه 0 درجه و پلوت 2 بعدی rz بوده و برای شار نوشتم که در تصویر متن برنامه نوشته شده track که ان را به flux تغییر دادم شما میتوانید با تغییر کلمات کلیدی وهمچنین تغییر پارامتر های شعاع و زاویه و... شکل های 2 بعدی و 3 بعدی خوبی ترسیم کنید.
در برنامه تغییراتی دادم و این بار pedep را حساب کردم شکل زیر
شار را نیز حساب کردم شکل زیر
در 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
چشمه را با 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% یا یک باشد.
ادامه مطلب ...کار با این کد ساده است
ابتدا برنامه را از قسمت دانلود ها دانلود کرده و نصب کنید
سپس برنامه را اجرا کنید
در ابتدا این صفحه ظاهر میشود
طبق شکل گزینه load the simulator را انتخاب کنید
بقیه مطلب را در قسمت ادامه مطالب ببینید
دانلود کد VMC یا Visual Monte carlo Program
از این کد بیشتر برای دزیمتری استفاده میشه
توضیحات بیشتر را در سایت این کد بخوانید
لینک دانلود نرم افزار
سرور 1
نمونه ای از کار هایی که با این کد انجام شده
شکل های هندسی زیر را ترسیم کنید ؟
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
جواب : دانلود کنید
برای رسم یک مکعب میبایستی 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 قطع میکند
حالا نوبت نوشتن سلول ها ست
ما در کل در این مساله دو سلول بیشتر نداریم یکی خود مکعب و هر چه داخلش هست و دیگری فضای بیرون مکعب هر چه هست
نکته : جلوتر برای جلو گیری از ثبت اندرکنش های ضاعد سعی میکنیم که مکعب را در یک کره فرض کنیم حالا با این کار یک سطح و سلول جدید اضافه میشه ولی داده ها ی بهتری در زمان کمتر بدست میاریم
ادامه را در قسمت ادامه مطالب مشاهده کنید
ادامه مطلب ...