Learning MCNP4C and MCNPX

پاسخ خانم رازقندی - Answers


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

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

محمدزاده

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

با سلام

من یک استوانه را mesh بندی کردم و از mesh tally نوع 3 استفاده کرده ام . که بعد از اجرای برنامه در خروجی داده ها به ازای z های خاص در r های  مختلف ergdeposit را داده حالا من می خواهم در خروجی به ازای rهای خاص در z های متفاوت بدهد چون می خواهم مثلا به ازای یک r مشخص تغییرات را در راستای z رسم کنم ،باید چه کار کنم؟


سلام

در حالت اول به ازای ظ های خاص پلات کانتور گرفتید یعنی

z=constant

teta and radius = variable

شعاع و زاویه متغییر بودند (برای اینها توزیع در نظر گرفتید)

z ها توزیعی نخواهند داشت

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

حالا برای شعاع خاص میخواهید پلات بگیرید مثل حالت قبل یکی از پارامتر ها که در اینجا شعاع است باید ثابت در نظر گرفته شود پس برای شعاع توزیعی ننویسید

برای z ها و زاویه ها باید توزیع بنویسید

منظورم از توزیع بازه بندی کردن اونهاست

مثلا

شروع 100i پایان


=========== پرسش و پاسخ جدید ===========

سلام  
ممنون که جواب دادید.اما من جوابم رو کامل نگرفتم.
برنامه زیر را نوشتم:

tmesh
cmesh3 total
cora3 0 100i 5
corb3 0 60i 6
corc3 360
endmd

در خروجی بر حسبz های یکسان rهای متفاوت مرتب شده من می خواهم خروجی برنامه در r های یکسان و zهای متفاوت مرتب شود
باید چه کار کنم.


با سلام

به این صورت بنویسید

tmesh
cmesh3 total
cora3 0 100i 5
corb3 0 60i 6
corc3 0 360
endmd

میبینید برای کور۳ توضیع تعریف نکردم (یعنی این پارت مد نظرم نیست اما دوتای بالا باید پارتیشن بندی بشن که قبلا شما انجام دادید)

حالا برای شرایط دیگه میتونید با این الگو عمل کنید .

پاسخ اقای مظاهری - Your answer

سلام

اقای مظاهری متن بلند بالایی نوشتید طبق خواسته شما از این پس سعی میکنم به میل شما بفرستم اما من هر وقت قصد اضافه کردن مطلبی دارم یا به صندوق پستی وبلاگ سر میزنم (اگه در اوقات فراغت از درس باشم که زیاد پیش میاد) سعی میکنم جواب ها رو همین جا پاسخ بدم حالا باید تدبیر دیگه ای رو اتخاذ کنم چون این شکلی گفتمان گونه پیش نمیره حتما باید یک فوروم در نظر گرفته بشه که تو فکرش هستم یا من یا یکی دیگه فرقی نمیکنه

اما در مورد سوالا من هیچ وقت ناراحت نشدم . اصلا ناراحتی نداره . احتمالا شما از نحوه پاسخ گویی من چنین برداشتید .

این عالیه که در کارشناسی این کد رو شروع کردید و باید بگم من و سایر دوستان نیز در حال یادگیری بیشتر هستیم چون همون نسخه 4c رو هنوز کسی کامل بلد نیست اما میتونید کتاب هایی در مورد این کد منتشر شده رو مطالعه کنید اما اونها هم از روی منوال هستند (من جزوه ای ندارم از منوال و کمک دو نفر از استاد های جوان منو بیشتر راه انداخت بقیه اش رو خودم تمرین کردم)

در مورد پراکنده گویی در وبلاگ باید بگم خیلی وقت ها دوستان اینو میگن ولی هر مطلب با عنوان خاص در شاخه خاصی قرار گرفته اما در بسیاری موقعیت ها اینو نمیشه رعایت کرد مثل استفاده از دستورات مشبندی که در جای خودش نیست ظرف یکی دو ماه ایند حتما درستش میکنم

در مورد گفته اساتید مورد اشاره باید بگم درسته منوال فقط یه کمکه و باید مواردی که با این کد کار شده رو تست کرد (مقایسه تجربه با شبیه سازی ها) البته در خیلی از موارد مثلا همین شبه سازی که من انجام میدم اصل دستگاه رو که نداریم باید از چندین مقاله اعداد رو بگیریم بعد اجرا کنیم یا .... این کار هاست که بیس خوبی برای فعالیت هی مهمتر بعدیه

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

در مورد نسخه های 2.4 و 2.6 هم چون از یک خانواده هستن تا الان یکجا گفتم بعدا مطالب جدا پست میکنم

در مورد بند اخر چوب کاری نکنید قبل از من 2 تا وبلاگ در حد چند پست وجود داشته پس من اولی نیستم قبل از اینها هم همگی در گروه های گوگل پرسش ها و پاسخ هایی داشتیم

نکته افتخار همه ما پرفسور علیرضا حقیقت هستند که بانی mcnp-polimi هستند و ارزوی من اینه که یه روز ملاقاتشون کنم

کد هایی مثل phits , mars هم هستند که تابحال سایتی که اونها رو اموزش بده ندیدم (فارسی). جا برای کار روی اونها هم هست بخصوص فیتس که خیلی از نظر ساختاری با mcnp نزدیکه و امکاناتش در شبیه سازی متحرک منو متعجب کرده و ازش خوشم اومده خلاصه پاسخ منم زیاد شد

همچنین کلی داده پردازی حاصل از کارهای خودم دارم که چند وقت دیگه بصورت هم پرسش پاسخ هم اموزش + روش های پردازش خروجی سریعتر و ... رو قرار میدم البته استفاده از مطلب رو قبلا گفتم مابقی رو باید توضیح بدم

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

iranhpc.com

امیدوارم موفق باشید

پاسخ به سوالات - Answers

سلام اول سوال اقای مظاهری رو جواب بدم چون پاسخ ایشون دلیل کاهش فعالیت منه


سلام.دوست عزیز مدتی است که به پیامهای ما جواب نمیدهی!
این رسم جوانمردی نبودا!
آقا شما کتاب آشکارسازی نول به زبان فارسی را دارید؟


اقای مظاهری من در حال تایپ و نهایی کردن پایان نامه ام هستم حق بدید گاهی اوقات نمیرسم

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

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


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

خانم ایدا

سلام پیام شما رو دریافت کردم

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

خانم هاشمپور

سلام پیام شما را نیز دریافت کردم

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

خانم ریاحی

سلام . مشکل شما با ادرس میلی که فرستادم حل شد ایشون جواب دادند ؟

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

اقای مظاهری

سوال قبلی رو در مورد راکتور اب سنگین

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

حوضه اطلاعاتم در مورد راکتور خیلی کمه ولی فرصت داشته باشم بیشتر مطالعه میکنم

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


پاسخ خانم هشمپور - Answer to Mrs Hashempoor


سلام

ممنون از راهنمایتون من دارم رو کسر جذبی بتا کار می کنم و از تالیF8*استفاده می کنم سلول من تنها یه بیضی و چشمه من هم تو کل حجم سلول گسترده شده و با توجه به اینکه mcnp چشمه بیضی تعریف نکرده من اونو داخل یه استوانه قرار دادم می خواستم بدونم من باید تو قسمت sdef از ccc استفاده کنم یا از cel نمونه ای از برناممو واستون می فرستم


سلام خانم هشمپور

برنامه رو به این ادرس میل بفرستید فایل (txt) :

physics.scholars@gmail.com

پاسخ به سوالات شما - your answers


خانم هاشمپور


واسه عضویت در این سایت و پرسیدن سوالاتم درباره کد mcnp باید چکار کنم ممنون میشم راهنماییم کنید


سلام

میتونید در قسمت تماس با من در بالای صفحه اول سوالاتتونو بپرسید

این سرویس وبلاگه و اپشن فعلا عضویت مهم نیست

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

که بصورت عضو پذیر و فوروم دار باشه تا در تالار پرسش ها و پاسخ ها رد و بدل بشه. 


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


پاسخ اقایان محمدی و مظاهری - your Answers

اقای محمدی

سلام

اقای محمدی من چشمه های بحرانی رو کار نکردم

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

پیشنهاد میکنم این رو در گروه های گوگل بپرسید امیدوارم پاسخ بگیرید


اقای حسین مظاهری

سلام

این کتاب رو خودم ندارم بنظر سخت میشه گیر اورد . سعی میکنم با همین عنوان اگه بین کتابام داشتم براتون بفرستم

موفق باشید


دانلود کتاب املیو سگره - Download Book

Finally i found Emilio segre's Book Nuclei and Particles , you can download it from bellow link

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


Download Emilio Segre , Nuclei and Particles 1st Edition

Size 85 MB , Format PDF


Source Link : Click Here

پاسخ به سوال اقایان مظاهری یزدانی و دادگر - your answers

اقای یزدانی

سلام

تمرین 2-2

این فایلایی که برا دانلود گذاشتی ایراد دارن!


سلام اقای یزدانی

منظورتون از ایراد چیه ؟ این لینک ها همگی سالم هستند چند بار چک شدند

مگه شما تا بحال mcnp کار نکردید ؟

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

اقای حسین مظاهری

سلام

خواهشمندم سولات دکتری مهندسی هسته ای گرایش راکتور حاوی زبان و استعداد تحصیلی و دروس تخصصی را برای دانلود بگذارید یا حداقل به ایمیلم بفرستید.


سلام اقای مظاهری

سوالات رو از سایت مدرسان شریف میتونید بگیرید البته نمیدونم الان هم سوالات رو برای دانلود داره یا نه ( همین سوالات رو به قیمت 50 هزار تومن بیرون میفروختند الان هم بهای اونها کم نیست خدا اخر عاقبت ما رو بخیر کنه ) .

سوالاتی رو که قرار دادم پرتوپزشکی و کاربرد پرتو ها و فیزیک هستند .


موفق باشید

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

اقای رضا دادگر

سلام

اطاعت لینک ها رو قرار میدم


موفق باشید

مشبندی - Meshing


قبل از هر صحبتی شما رو به صفحه ۱۶۶ / ۴۲۷ mcnpx نسخه 2.4 ارجاع میدهم .

تعریف مشبندی :

همانطور که میدانید پیاده سازی یک هندسه از سیستم مورد نظر وقت گیر همراه با خطا و دشواریهای خود را دارد حال اگر به جزییات بیشتری احتیاج داشتیم (محاسبه شار در بخشی خاص یا ... به هر ترتیب باید سطح و سلول تعریف شود پس ... ) میبایستی دوباره هندسه قبلی را دستخوش تغییر قرار دهیم و در این حالت کار دشوارتر میشود اما در نهایت با قربانی کردن زمان برنامه را مینویسیم . برای اینکه ساختار برنامه اولیه یا بهتر توصیف کنم اسکلت اولیه مساله دستخوش تغییر قرار نگیرد تغییرات مرحله دوم را که میبایست اعمال شوند را بصورت مجازی با عکلگر های مشبندی وارد میکنیم

در این مرحله هیچ کاه سطح یا سلولی تعریف نمیکنیم فقط به کد میفهمانیم که مساله ما را با استفاده از این دستگاه مختصات (دکارتی یا استوانه ای یا کروی که باید تقارن را نیز در نظر داشته باشید مثلا برای یک کره از مشبندی مستطیلی یا همان دستگاه دکارتی استفاده نکنید) شبکه بندی کن (شبکه های کوچک)

مشبندی نوع اول : track averaged mesh tally

از این مشبندی برای محاسبه : شار / انرژی توزیع شده / دوز / و ... استفاده میکنیم

پرکاربرد ها رو در بالا اشاره کردم

tmesh                                                                          
rmesh11:p pedep                                                                 
cora11 -25 25                                                                  
corb11 -25 100i 25                                                             
corc11 0 40i 50
endmd

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

همیشه tmesh شروع کار و endmd پایان کار است

rmesh یعنی مشبندی مستطیلی

cora=x

corb=y

corc=z

cmesh یعنی مشبندی استوانه ای

smesh یعنی مشبندی کروی

عدد جلوی rmesh باید 1 / 11 / 21 / 31 / ... عدد اخر باید 1 باشد که معرف مشبندی نوع اول است به مثال توجه کنید !

در بالا مختصه x هیچ تغییری نمیکند تنها مختصه های y , z را قسمت بندی کرده ایم به عبارت دیگر بصورت دوبعدی مساله را بررسی میکنیم و تغییرات انرژی توضیع شده را در راستای y تحلیل خواهیم کرد تصویر زیر یک مثال از این نوع است :

نکته : شماره مشتالی هایی که مینویسید نباید با تالی های اصلی مساله یکی باشد . (به پانویس ص 145 منوال 2.4 دقت کنید )

تصویر بالا پلات کانتور مشبندی مستطیلی است که میزان انرژی توزیع شده فوتون را در اب نشان میدهد

البته عملگر pedep معادل F6 است منتها با یک تفاوت که F6 انرژی توزیع شده در سلول به ازای واحد جرم سلول است و pedep به ازای واحد مشبندی شبکه . این دو اختلاف اندکی دارند که با انتخاب درست دستگاه مختصات و اندازه واحد های شبکه به حدقل میرسند . البته هر دو مقدار قابل تایید اند .

مشبندی نوع دوم : مشبندی چشمه

برای حفاظ سازی به کاربر اجازه میدهد تا نشتی ها رو مشخص کند (منوال 2.6 صفحه 198/551)

tmesh
rmesh2 p trans
cora2 -25 25 
corb2 -25 100i 25
corc2 0 50i 50
endmd

مساله قبل را برای یافتن نشتی فوتونی با همان ساختار بستیم

شما میتوانید برای ذرات مختلف این کار را انجام دهید

نکته : شماره مش تالی باید 2 یا 12 یا 22 یا ... عدد اخر به 2 ختم شود

مشبندی نوع سوم : Energy Deposition Mesh Tally

نتایج این مشتالی با تالی F6+ یکسان است حتما میدانید که F6+ انرژی توزیع شده به ازای واحد حجم را بدست میدهد . (صفحه ۱۹۹/ ۵۵۱ منوال ۲.۶)

tmesh
rmesh3 de/dx
cora3 cora11 -25 25
corb3 corb11 -25 50i 25
corc3 corc11 0 50i 50
endmd

کلمات کلیدی که در این مشبندی استفاده میشوند :

dE/dx توان توقف ذره تابیده شده در محیط

recol میزان طیف بازتابش شده در محیط

و ... سایر عملگر ها را میتوانید در منوال مطالعه کنید

نکته : شماره مشتالی کماکان طبق روالهای ذکر شده در بالا این بار عدد اخر باید به 3 ختم شود


خروجی این مشتالی ها بصورت فایل mdat ذخیره میشود که باید با gridconv تبدیل کنید به ورودی مناسب برای برنامه هایی مثل tecplot (من از این استفاده میکنم) و gnuplot و دو تا برنامه دیگه که در طی مراحل تبدیل اسم این برنامه ها رو مینویسه . قبلا نحوه تبدیل این فایل ها رو توضیح دادم امید وارم این توضیحات مفید واقع بشه هرچند چند بار اصلاح شد و حجم متنش رو به حداقل رسوندم تا خسته کننده نشه . موفق باشید


PageRank