پاسخ اقای مولایی - Answer to alireza molaei

سلام
من  mcnp vised را از همین وبلاگ دانلود کردم. مثالهای  سایته mcnpvised.com  را باش run کردم ،مشکلی نداشت
ولی هر کار میکنم قسمت plot particle track  جواب نمیده ، با اینکه خود مثال  آموزشیو   خط به خط اجرا میکنم.
میخواسم بپرسم شما تا حالا تونستین از   particle display  استفاده کنین؟

لینک مثالها
http://www.mcnpvised.com/sample_exercises/sample_exercises.htm

اقای مولایی سلام

برنامه ریجستر نشده هست ولی میتونید خیلی کار ها انجام بدید فقط ترک ذرات در سه بعد رو - بصورت سه بعدی - غیر فعاله همین ایرادش فقط همینه . تصویر زیر ترک الکترون در فانتوم اب هست از نمای بالا

من فکر میکنم مسیری رو که طی کردین یه جاهایی رو نادیده گرفتین

بگذریم این مثال ساده رو داشته باش

وقتی در visual editor به قسمت particle tracks برین پنجره زیر را خواهید دید

شماره یک - در این اپشن شما نوع نمایش ترک ها رو تعیین میکنید :

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

شماره دو - باید تیک جلوی این ذره که در این مثال نوترون است زده بشه در غیر این صورت هیچ ترکی رو نخواهید داشت چون با این تیک رد نوترون ها در محیط دنبال میشه

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

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

شماره 4 : میتونید حد انرژی تعریف کنید برای ترک ها

شماره 5 : اگه میخواهید ترک ها ترسیم شوند این تیک رو بزنید در غیر این صورت ذرات ثبت شده بصورت نقطه هایی روی صفحه مشهاده خواهند شد که باز هم میتونید اندازه این نقاط رو کوچک و بزرگ کنید

شماره 6 : برای انرژی های پایین و انرژی های بالا و ترک ها رنگ جدا گانه انتخاب کنید

شماره 7 : ضخامت خطوط ترک ها را میتوانید تغییر دهید

در نهایت پلات دو بعدی تالی که نوشتم اینطور در میاد

اطلاعیه

Some articles about contamination and Dosimetry in photon beams

A MonteCarlostudyonneutronandelectroncontaminationofan unflattened18-MVphotonbeam

Download

A review of dosimetry studies on external-beam radiation treatment with respect to second cancer induction_X George Xu,...

Download

Calculation of electron contamination doses produced using blockingtrays for 6 MV X-rays_Martin J. Butsona

Download

Correcting for electron contamination at dose maximum in photon beams_D. W. O. Rogers

Download

Depth for dose calibration in high energy photon beams

Download

Dosimetric Characteristics of Clinical Photon Beams_Jatinder R Palta

Download

Monte Carlo estimation of electron contamination in a 18 MV clinical photon beam_Mesbahi

Download

reducing electron contamination_LiRogers-1994

Download

Simulation and measurement of air generated electron contamination in radiotherapy_Martin J. Butson

Download

پاسخ اقا/خانم Karrabi

با عرض تشکر از راهنماییهای شما من هنوز تازه یادگیریر
mcnpرا شروع کردم و به خیلی چیزها واقف نیستم در تمرین های خودتان(ex 77) زیر چرا  عددی که برای دانسیته جرمی در کارت سلول  برای سلولهانوشته اید 1 است در حالیکه در کارت ماده مواد متفاوتی را  در هر سلول تعریف کرده اید

سلام

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

فرکشن های جرمی

فرکشن های اتمی

فرکشن های جرمی رو با منفی مینویسند و فرکشن های اتمی رو با عدد مثبت

میتونید این فرکشن ها رو تو کتاب زیر پیدا کنید

Compendium of Material Composition Data for Radiation Transport Modeling_R.G. Williams III,C.J. Gesh,R.T. Pagh_April 2006_Prepared for the U.S. Department of Energy under Contract DE-AC05-76RL01830

Download it

امید وارم منظورم رو فهمیده باشید اگر نه یاداوری کنید با مثال پست کنم

پاسخ خانم فهیمه - Answer to Mrs Fahimeh

سلام.
چه طور می شه در یک سلول شش وجهی که در مرکزش یک استوانه است داخلش رو از استوانه های دیگه ای که کاملا مشابه هستند پر کرد؟

سلام

سلول شش وجهی و استوانه داخلش رو تعریف کنید

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

مثال :

فرض میکنم سلول شش وجهی شماره ۱۰۱

سلول استوانه ای شماره ۱۰۲ باشه و در مرکز شش وجهی باشه حالا سلول بعدی رو اینطور مینویسم

103 Like 102 But Trn(مختصات x y z)

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

این نوع دستور به دستور Like But معروف است برای مطالعه بیشتر به منوال رجوع کنید

پاسخ اقای اوسط - Answer of Mr Aosat

با سلام
در شبیه سازی سر شتاب دهنده ی پرتو درمانی در یک پایان نامه ذکر شده بود که بیم الکترونی با پهنای باریکه fwhm= 2 میلی متر به هدف برخورد میکند
بنده نتوانستم در ام-سی-ان-پی این سورس را شبیه سازی کنم.در واقع منظور این جمله را اصلا نمیفهمم
در صورت ممکن بنده را راهنمایی بفرمایید .

اقا سلام 

من با این اپشن کار نکردم ولی

FWHM یعنی Full Width Half Maximum و برمیگرده به نحوه ی تابش طیف از چشمه با توجه به انچه که شما ذکر کردید منظور از متن پایان نامه این بوده که چشمه مورد نظر از تابع توزیع گوسی که مقدار نیم پهنای ماکزیمم ان 2 است طبعیت میکند و برای تعریف چنین حالتی میبایستی از کارت FTn که به Special Treatment Cards مربوط میشه با اپشن GEB کار کنید

به ص 185 / 898 منوال 4C و همچنین ص 398/898 از همین منوال رجوع کنید

پاسخ اقا/خانم S Karrabi

سلام

سوال اول در مورد برنامه زیر بود و دلیل خطا :

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

fatal error.     1 tally volumes or areas were not input nor calculated

سلول 30 رو بصورت زیر بنویسید :

و اما در مورد فایل دومی که فرستادین مرز بلاک ها رعایت نشده بود بطوری که کارت سطح را بصورت دوبلاکه نوشته بودید کارت داده ها سه تکه نوشته شده بود - البته این فواصل رو وقتی فایل txt شما رو باز کردم دیدم - مطمئنم به این مورد اشراف کامل دارید اما یه نکته سعی کنید از عملگر # استفاده نکنید یا کمتر استفاده کنید هندسه سیستمی که طرح ریزی کرده اید انقدر پیچیده نیست که نتوان سطوح سلول ها رو براحتی نوشت اینو هم به این خاطر گفتم که هنوز سر استفاده از این عملگر توافقی بین کد نویس ها صورت نگرفته برای راحتی کار ازش استفاده میکنند ولی بعضی جاهای حساس کار دست ادم میده بگذریم 

نوع چشمه ای که نوشتید چیه ؟ فوتونیه / نوترونیه / چیه ؟ از روی mode p میشه حدس زد فوتونیه و سلول 4 که پر از اب هست هم به عنوان چشمه انتخاب شده 

با این اوصاف چشمه شما isotropic volumetric source بحساب میاد پیشنهاد میکنم به منوال an mcnp primer  ص 17/45 مراجعه کنید

+ چرا imp:n  ؟ انرژی چشمه فوتونی اگه تقریبا 6 مگا الکترون ولت باشه میتونه نوترون ازاد تولید کنه - ایا چشمه ای که مد نظر دارین به این انرژی میرسه ؟ 

فایل اصلاح شده رو میتونید از این لینک دانلود کنید 

سوال شما : فانتوم چیست ؟ فانتوم از اسمش پیداست شبح یا بهتره بگیم بدل یه موجود زنده مثلا یه عضو از بدن را با کد طراحی میکنیم تمامی تلاشمون رو روی درست طراحی کردن از نظر هندسی و حجم و غیره انجام میدهیم و در نهایت چیزی که درست شده فانتوم اون عضو از بدنه حالا فانتوم های معروف یه تعدادی هستند مثل فش و مش و اشنایدر و ... بگذریم البته میتونیم بجای فانتوم عضو بدن یه حجم رو در نظر بگیریم و درونش اب بریزیم به این فانتوم ها water phantom میگن و تالی هایی که از فانتوم اب گرفته میشه بخصوص در پزشکی همخوانی خوبی با داده های واقعی دارند بخاطر همین اکثرا دستگاه هایی رو که طراحی میکنند روی فانتوم اب ازمایش میکنند و بعد روی فانتوم بدن و پس از اتمام کار شبیه سازی کار ها بصورت عملی انجام میشه تا نتایج با تجربه مقایسه شوند 

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

در موضوعات به بخش mcnpx مراجعه کنید 

پاسخ حسین - Answer to hossein

سلام. من چندتا مشکل با mcnpx  دارم. اول اینکه بیشتر کسانی که در ایران روی این کد کار می کنند روی ترابرد نوترون و فوتون متمرکز شده اند. من دوست دارم روی ترابرد الکترون و پروتون هم کار کنم ولی ظاهرا کسی نیست کمکی بهم بکنه. شما کسی رو یا رفرنسی رو میشناسید به من معرفی کنید؟ دوم اینکه وقتی الکترون رو ترابرد می کنم ارور میده که مواد تعریف شده رسانا هستندconductive .چه جوری تعریف مواد با ایزوتوپ هاشونو انتخاب کنم یا مثلا برای پروتون چه طوری سطح مقطع ها رو از کتابخونه اش انتخاب کنم؟ مورد سوم اینکه وقتی کد رو RUN می کنم با یه کامپیوتر core I 7 بیش از 2 میلیارد ذره نمی تونه ران کنه و خطا هم زیاده؟مشکل بعدی اینه که در یه قسمتی از هندسه ی من تالی صفر میشه و من نمی دونم مشکل از کجاست؟میشه راهنماییم کنید.

سلام

اول اینکه میشه گفت درسته تو ایران بیشتر علاقه مندی کاربران این کد به نوترون معطوف میشه مثل کار روی راکتور / حفاظ سازی / پزشکی مثلا bnct ولی در هر اندرکنشی که فکر کنید کم کمش طیف گاما یا ایکس یا در مواردی میشه گفت الودگی فوتونی و الودگی الکترونی و سایر الودگی ها رو میتونید مشاهده کنید . برای مثال میگم من هم روی ترابرد فوتون و الکترون دارم کار میکنم ولی اینو میدونم که وقتی انرژی طیف فوتونی من از 6 مگاالکترون ولت بیشتر بشه احتمال اینکه در محیط مفروض الودگی نوترونی داشته باشم رفته رفته با افزایش انرژی طیف بالا میره پس من نوترون ها رو هم نباید فراموش کنم - به بیان دیگه هر کسی که کد کار میکنه باید محدوده ها رو بدونه تا به بیراهه کشیده نشه -

همچنین اگه مایل هستین روی پروتون کار کنید باید بدونید که کد mcnp برای ترابرد این ذره از مدل استفاده میکنه - مدل های متفاوتی ارائه شده برای محدوده های خاص - من تا بحال با مدل کار نکردم ولی میتونید از منوال اطلاعات بیشتری رو بدست بیارید .

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

مرجع تصویر :

mcnpx transport in the tabular regime _ H.grady Hughes

Hadronic Shower Simulation Workshop

Fermilab, 6-8 September 2006

اگه فایل Pdf رو خواستید با من تماس بگیرید [e-mail]

و اما ادامه خیر کسی رو نمیشناسم با این کد روی پروتون ها کار کنه بیشتر میرن سراغ Geant و Phits یا Mars البته دوتای اخری رو جایی ندیدم ازشون صحبتی بشه فکر نمیکنم اینجا کسی این کد ها رو داشته باشه . توصیه میکنم از منوال استفاده کنید + گروه mcnpx گوگل هم میتونه مشکل گشا باشه

مورد دومی رو که گفتید درست متوجه نشدم وقتی ترابرد الکترون رو بررسی میکنی برنامه fatal error میده ؟ من فکر میکنم باید هشدار باشه اگه warning بده میتونی با نوشتن cond=1 در انتهای متریالی که رساناست و cond=0 برای متریال نارسانا از شر هشدار خلاس شی اما اگه خطا میده باید برنامه رو بررسی کنی رو هوا نمیشه صحبت کرد ونظر داد

درمورد پروتونها پیشنهاد میکنم به منوال رجوع کنید

مورد سوم اینکه وقتی کد رو RUN می کنم با یه کامپیوتر core I 7 بیش از 2 میلیارد ذره نمی تونه ران کنه و خطا هم زیاده؟مشکل بعدی اینه که در یه قسمتی از هندسه ی من تالی صفر میشه و من نمی دونم مشکل از کجاست؟

حالا چرا تعداد ذرات براتون مهمه ؟ برنامه شما باید از نظر 10 چک اماری حکم yes رو اتخاذ کنه

با تعداد ذرات پایین شروع کنید ببینید خطا چنده اگه روند رو به کاهش داشت ارام ارام تعداد رو اضافه کنید تا خطا رو به زیر 0.05 یا 5% برسونید البته ممکنه با این کار Fom پایین بیاد یا نوسان داشته باشه همچنین slop دست خوش تغییر بشه باید تعداد ذرات را تا انجا اضافه کنید که از تمام چک های اماری yes بگیرید البته در کار با ذرات سنگین روش های کاهش واریانس هم نقش دارند و باید این روش ها رو هم چاشنی کنید

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

موفق باشید

علی اصغر محمدزاده

[e-mail]

پاسخ اقای دادگر - Answer to Mr Dardgar

اقای دادگر سلام

زحمتی نیست ترویج علم رحمته و اما پیرامون سوال های شما :

1 - خیر نرم افزار FastRad رو هنوز نتونستم مفتش رو پیدا کنم فکر نمیکنم تو ایران کسی داشته باشه

2 - باید فایل Catia رو تبدیل کنید به ورودی مناسب برای اتوکد برای اشنایی بیشتر به ادرس های زیر سر بزنید - ( من روی سیستمم Catia یا Solid Work ندارم اگه داشتم حتما تست میکردم - خلاصه از این بابت شما رو به لینک 1 و لینک 2 ارجاع میدم )

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

موفق باشید

Dpm Download link

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

دانلود DPM 11 

Download DPM 11