پاسخ خانم فهیمه - 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 ارجاع میدم )

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

موفق باشید

پاسخ اقای فریدی - Answer to Mr Faridi

سلام
با راهنمایی های خوبتون و استفاده از منوال صفحه 356/898 تونستم  چشمه حجمی یکنواخت تولید کنم ولی وقتی plot source رو اجرا میکنم ذرات چشمه در خارج از حجم سلول هم تولید میشن که احتمالا در نتایج تاثیرکذار هستند.
چطور میتونم این مشکل رو رفع کنم
اگه لازم میدونید یه ایمیل بزنید تا برنامه رو براتون بفرستم
با سلام 

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

[e-mail]

یه نکته دیگه هم هست چون درست خاطرم نیست چشمه شما رادیو ایزوتوپ باید باشه و رو این حساب روشی رو که گفتم یک پای مساله حساب کنید و مهمتر اینکه در گروه گوگل تیم mcnpx پس از اینکه وارد شدید دنبال فردی بنام Biagio Zaffora بگردید یه پست داره تحت عنوان spectre gamma by tally F8 پس از پیدا کردن این پست یه سری ردوبدل شدن اطلاعات راجعبه اکتیویته و ... رو میبینید متن زیر گزیده ایست از پاسخ زافورا زبانش ایتالیایی است اما با مترجم گوگل برگردونیدش به انگلیسی اونوقت راحت میشه فهمید :

Voici un extrait d'un travail que j'ai fait il y a quelque temps

Le fichier d’entré du spectre de Cs-137 est donné en Annexe 5. Observons les commandes nécessaires à la création du spectre

C  Physical Data

MODE        P

IMP:P        1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 0 1

SDEF        CEL=17  RAD=D1  EXT=D2 AXS= 0 0 1  ERG=D3  PAR=2  WGT=73592000

F8:P          3

FT8  GEB  .0003493  .001110908  0

SI1            0  5

SI2            6 8

SI3   L       .031817 .032194 .036350 .037340 .66160

SP3 D       19.3  35.6  10.5  2.5 852

E8             0 1E-5 1E-3 16380I 1.9488

Les nouvelles fonctions WGT, FT8 GEB, SIn L et SPn D apparaissent. Expliquons-en le sens

Dans cet exemple on a voulu créer le spectre du Cs-137 avec les pics à 31,817 keV, 32,194 keV, 36,35 keV, 37,34 keV et 661,66 keV d’intensités respectives 1,93 %, 3,56 %, 1,05 %, 0,25 % et 85,2 %. Pour ce faire on a associé à la fonction ERG la variable Dx (D3 dans cet exemple) laquelle fait référence à la source d’information SI3 contenant les énergies du spectre de Cs-137 exprimées en MeV. La lettre L indique qu’il s’agit d’un spectre discret avec les énergies indiquées sur la même ligne de commande

Egalement, on associe à la variable D3 une source de probabilité SP3 que, dans le cas de la spectrométrie gamma, fournit la probabilité d’émission de la source à l’énergie correspondante. Dans le cas présent, par exemple, en correspondance de l’énergie 0,031817 MeV on a une source de probabilité donnée par

SP(E)=A*I(E)=1000*0,0193=19,3

La lettre D suivante la fonction SPn indique simplement qu’il s’agit d’une distribution de probabilité

Etroitement liée aux deux commandes qu’on vient d’analyser est la commande WGT laquelle permet de prendre en compte le temps de comptage du spectre simulé et l’activité totale. La formule mathématique de WGT est la suivante

WGT = (ΣΣ A*I)*t 

avec t exprimé en seconds

Dans le cas du seul Cs-137 WGT est calculé ainsi

 

WGT = (A*I1 + A*I2 + ... + A*I5) = (1000*0,0193 + 1000*0.0356 + ... + 1000*0,852) * 80000= 73592000

 

Dans le cas de radionucléides avec activités defferentes la formule reste la même mais tu auras des activités A qui changent 

Cordialement,

 

Biagio Zaffora

ادرس ایمیلش اینه : [e-mail] پیشنهاد میکنم یه مشورت ازش بگیر سابقه بیشتری داره تو این کد

پاسخ اقای مجید موسوی - Answer to Majid musavi

سلام
من visual editor که شما تو سایتتون گذاشته بودید رو دانلود کردم ولی کتابخانه موجود تو فولدر  BIN رو نمی خونه . و وقتی  read cross section  رو کلیک می کنم خطای cant find xsdir رو می ده . در صورتی که  xsdir_old تو فولدر هست.

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

خسته نباشید شما input file خودتونو open نکردید بعلاوه خودتون که میگین cross section این فایل حجمش خیلی بیشتر از 150 کیلوبایت xdir-old هست باز هم میگم این نرم افزار یه واسط گرافیکی متنیه سورس mcnp نیست بلکه از mcnp استفاده میکنه چه 4c چه x . توجه داشته باشید برای مشاهده سطح مقطع ها حتما باید material  تعریف کرده باشید . 

پاسخ اقای فریدی - Answer to Mr Faridi

سلام و خداقوت
وقتی یک چشمه حجمی با توزیع یکنواخت و تک انرژی به صورت سلول محاط در کره تعریف می کنم و sp رو با پارامتر 2    21-  معرفی می کنم این خطا رو میده:
the sampling efficiency in source cell  10 is too low
ولی وقتی به عدد دیگری مثل 100   0  تبدیل می کنم درست میشه.
علت چیه و چه اعدادی رو طبق چه قاعده ای باید انتخاب کنیم؟ سپاس

سلام

شما کارت توزیع احتمال چشمه تون رو بر حثب تابع نوشتین ص 353/898 منوال 4C توضیح داده

اگر برای چشمه راستای تابش تعریف کردین بجای a باید مقدار 1 قرار بدین اگر توزیع زاویه ای تعریف کردین با Rad باید بجای a عدد 2 را بنویسید و در هر صورت با این کارت شما با این ساختار دارید مشخص میکنید که توضیع شار چشمه چطوری باشه 21- هم مشخص کننده تابعی است که توضیع بر حثب اون انجام میشه

سلول 10 همون سلولی هست که چشمه رو در اون قرار دادید ؟ انرژی چشمه چقدره ؟ پارامتر هایی که برای چشمه معرفی کردید کدام ها هستند ؟ Rad ? Dir ? این خطایی که داده داره میگه که بهره چشمه شما کمه با این تغییری که شما دادین یعنی تعویض اعداد قبلی با 0 و 100 دیگه نمیشه گفت توزیع احتمالاتی شار چشمه شما چطوریه ؟ قبول دارید ؟ چشمه شما نوترونیه / فوتونیه یا الکترونی ؟

پاسخ اقای اوسط - Answer to Mr Ousat

با سلام
سوال 1 - ورودی کد ام سی ان پی فانتوم rando را از کجا میتوانم تهیه کنم؟
زیاد سرچ کردم و همه ی فانتوم ها را پیدا کردم اما این یک مورد را نتوانستم
با تشکر

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

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

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

سوال 2 - آیا برنامه ی scan2mcnp را دارید؟

قیمت این برنامه 750 $ میباشد
در صورت ممکن اگر آنرا پیدا کردید برای ما هم مفید خواهد بود

نه این برنامه رو ندارم . چشم اگر پیدا کردم حتما میزارم رو سایت

موفق باشید

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

سلام و خسته نباشید. من تازه کار با mcnp رو شروع کردم .ژئومتری مساله ام رو تا حدی انجام دادم اما با پلات گرفتن مشکل دارم.4c خط در میون ارور میده.x هم این ارور رو می ده. xlib:  unable to connect to server display 
لطف می کنید راهنمایی بدید من چه کار کنم.

سلام

من نمیدونم شما میخواین چی رو پلات کنید هندسه رو اگه میخواین وارسی کنین از visual editor استفاده کنید نسخه 4C هم این قابلیت رو داره تا بخش های مختلف هندسه رو نمایش بده اگه خطا داره پس حتما جایی رو اشتباه دارین اما نسخه x برای دیدن هندسه در این نسخه باید Xming رو نصب و همزمان اجرا کنید در غیر این صورت همون خطایی رو نشون میده که ذکر کردین

از Visual editor استفاده کنید در وقت صرفه جویی میشه

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

موفق باشید