<div>Hello colleagues,<br><br> I've set up a residual stress instrument.The simulation result suggests(Neutron count=10<font size="1">9<font size="3">)</font></font>: Neutrons after Fe(211) diffraction are I I_err N=[29927.1 31233.5 3];And neutrons reaching the detector through aperture are I I_err N=[0 0 0]. Both I and N are too low.<br> I raised horizontal divergences of the first and second collimator but I and N remain too low. What should I do in order to get high I and N on detector?<br> <br> Here is the code:<br>/*******************************************************************************<br>* McStas instrument definition URL=http://www.mcstas.org<br>*<br>* Instrument: residual stress instrument <br>*<br>* %Identification<br>* Written by: Tom<br>* Date:10 Jan 2008<br>* Origin: institute<br>* Release: McStas 1.6<br>* Version: 0.2<br>* %INSTRUMENT_SITE:thesis<br>*<br>* residual stress instrument <br>*<br>* %Description<br>* try to optimize parameters of residual stress instrument. default geometry is <br>* from blueprints of residual stress instruments at CARR,FRM-II(Germany),E3* (Germany) * and SALSA(France) <br>* Example: mcrun residual_stress.instr <RV=-1><br>*<br>* %Parameters<br>* dist: [m] focusing distance from source to SlitIn<br>* GapM2_G: [m] distance between M2Exit and guide entrance<br>* xw: [m] width of neutron guide entrance<br>* yh: [m] height of neutron guide entrance<br>* l: [m] length of straight neutron guide<br>* m: [1] m-value of material1.5~3<br>* CH1: ['] horizontal divergence of the first collimator<br>* len_C1 [m] length of the first collimator<br>* L1: [m] distance between exit of first collimator and <br>* monochromator,0.01~0.3m<br>* A1: [deg] rotation of mono<br>* DM: [Angs] d-spacing of monochromator<br>* NH: [columns] number of slabs horizontal <br>* NV: [rows] number of slabs vertical <br>* mosaich: ['] Horisontal mosaic FWHM <br>* mosaicv: ['] Vertical mosaic FWHM <br>* RV: [m] Monochromator vertical curvature, 0 for flat, -1 for automatic<br>* setting(then RV=2*L2*sin(DEG2RAD*A1) )<br>* RH: [m] Monochromator horizontal curvature, 0 for flat, -1 for automatic<br>* setting(then RH = 2*L2/sin(DEG2RAD*A1))<br>* L2: [m] distance between monochromator and sample,1.9~2.2<br>* CH2: ['] horizontal divergence of the second collimator<br>* DS: [Angs] d-spacing of sample:default Fe(hkl) <br>* L3: [m] distance between sample and detector <br>* %Link<br>* A reference/HTML link for more information<br>*<br>* %End<br>*******************************************************************************/<br>/* Change name of instrument and input parameters with default values */<br>DEFINE INSTRUMENT residual_stress(<br> dist=2, GapM2_G=0, xw=0.06, yh=0.11, l=6, m=2, CH1=10,<br> L1=0.3, A1=45, DM=1.0886, len_C1=0.6, NH=1, NV=15, mosaich=30,mosaicv=30,<br> RV=-1, RH=0, L2=1.9,CH2=30,DS=1.1782,L3=1.0) <br>/* The DECLARE section allows us to declare variables or small */<br>/* functions in C syntax. These may be used in the whole instrument. */<br><br>DECLARE<br>%{<br>double M2Exit_mono; /* distance between M2Exit and monochromator */<br>double lambda; /*wavelength after diffraction*/<br>double monow; /*width of mono*/ <br>double A2; /*half of sample take off angle*/<br>%}<br><br>/* The INITIALIZE section is executed when the simulation starts */<br>/* (C code). You may use them as component parameter values. */<br>INITIALIZE<br>%{<br>M2Exit_mono=GapM2_G+l+len_C1+L1;<br>lambda=2*DM*sin(DEG2RAD*A1);<br>if (RV<0) RV=2*L2*sin(DEG2RAD*A1);<br>if (RH<0) RH = 2*L2/sin(DEG2RAD*A1);<br>monow=xw/sin(DEG2RAD*A1);<br>A2=asin(lambda/(2*DS))*RAD2DEG;<br>printf("\nlambda=%f\tRV=%f\tRH=%f\tmonow=%f\tA2=%f\n",lambda,RV,RH,monow,A2); <br>%}<br><br>/* Here comes the TRACE section, where the actual */<br>/* instrument is defined as a sequence of components. */<br>TRACE<br>/* Progress_bar is an Arm displaying simulation progress.*/<br>COMPONENT Origin = Progress_bar()<br> AT (0,0,0) ABSOLUTE<br><br>/********************* neutron source********************************/<br>COMPONENT Source = Source_Maxwell_3(<br> height = 0.75, width =0.52988 , <br> l_low =0.2,l_high =7, <br> dist = dist, <br> xw = 0.06, yh = 0.11,<br> T1 = 320,I1 =9.87E12) <br> AT (0, 0, 0) RELATIVE Origin<br><br>/************channel(Al and Fe) inside reactor, length:4.049m**********/<br>COMPONENT SlitIn = Slit(<br> width = 0.06, height = 0.11)<br> AT (0, 0, 0 ) RELATIVE Source<br><br>COMPONENT M2Exit = Slit(<br> width = 0.06, height = 0.11)<br> AT (0, 0, 4.049 ) RELATIVE SlitIn <br><br><br> /*************************straight_guide*******************/<br>COMPONENT straight_guide = Guide(<br> w1 = xw, h1 = yh, w2 = xw, h2 = yh, l =l, m = m)<br> AT (0, 0, GapM2_G) RELATIVE M2Exit<br><br><br>/**************************first collimator**************************/<br>COMPONENT collimator1 = Collimator_linear(<br> xmin = -xw/2, xmax = xw/2, ymin = -yh/2, ymax = yh/2,<br> len = len_C1, divergence = CH1)<br> AT (0, 0, l) RELATIVE straight_guide<br><br><br>/*************************mono_craddle***************************/<br>COMPONENT mono_craddle = Arm(<br> )<br> AT (0, 0, M2Exit_mono) RELATIVE M2Exit<br> ROTATED (0, A1, 0) RELATIVE Origin<br><br>/***************************monochromator Ge(511)****************************/<br>COMPONENT monochromator = Monochromator_curved(<br> width=monow, height=0.15, NH=NH, NV=NV,reflect = "0", transmit = "0",<br> gap = 0.0005, mosaich = mosaich, mosaicv = mosaicv, r0 = 0.7, <br> t0 = 1.0, RV = RV, RH = RH, DM = DM)<br> AT (0, 0, 0) RELATIVE mono_craddle<br><br>/*********************mono_out*******************/<br>COMPONENT mono_out = Arm(<br> )<br> AT (0, 0, 0) RELATIVE mono_craddle<br> ROTATED (0, 2*A1, 0) RELATIVE Origin<br><br>/********************detectors before second collimator**************************/<br>COMPONENT LDbeforeC2 = L_monitor(<br> nchan = 200, filename = "LDbeforeC2", xwidth = 0.15,<br> yheight = 0.15, Lmin = 0, Lmax = lambda+0.1)<br> AT (0, 0, 0.8) RELATIVE mono_out<br><br>COMPONENT PSDbeforeC2 = PSD_monitor(<br> filename = "PSDbeforeC2", xwidth = 0.15, yheight = 0.15)<br> AT (0, 0, 0.8) RELATIVE mono_out<br><br>/***********************second collimator***************************/<br> COMPONENT collimator2 = Collimator_linear(<br> xmin = -xw/2, xmax = xw/2, ymin = -yh/2, ymax = yh/2, len = 0.4,<br> divergence = CH2)<br> AT (0, 0, 0.8) RELATIVE mono_out<br><br>/**********************detector after second collimator********************/<br>COMPONENT PSDafterC2 = PSD_monitor(<br> filename = "PSDafterC2", xwidth = 0.06, yheight = 0.11)<br> AT (0, 0, 0.4) RELATIVE collimator2<br><br><br>/**********************aperture before sample*************************/<br> COMPONENT slit_before = Slit(<br> width = 0.008, height = 0.008)<br> AT (0, 0, 1.88) RELATIVE mono_out<br><br>/**********************detector behind Slit_before***********************/<br>COMPONENT PSDafterSlit_before = PSD_monitor(<br> filename = "PSDafterSlit_before", xwidth = 0.01, yheight = 0.01)<br> AT (0, 0, 0.01) RELATIVE slit_before<br><br><br>/************************Fe(221) sample****************************/<br>COMPONENT Sample = PowderN(<br> reflections = "Fe.laz", radius = 0.004, yheight = 0.8)<br> AT (0, 0, L2) RELATIVE mono_out<br>EXTEND<br>%{<br>if(!SCATTERED) ABSORB; /*perfect beamstop*/<br>%}<br><br>/************************sample_out************************/<br>COMPONENT sample_out = Arm(<br> )<br> AT (0, 0, 0) RELATIVE Sample<br> ROTATED (0, -2*A2, 0) RELATIVE Sample<br><br>/************************detector before slit_after********************/<br>COMPONENT PSDbeforeSlit_after = PSD_monitor(<br> filename = "PSDbeforeSlit_after", xwidth = 0.008,<br> yheight = 0.050)<br> AT (0, 0, 0.01) RELATIVE sample_out<br><br>/************************aperture after sample***************************/<br>COMPONENT slit_after = Slit(<br> width = 0.008, height = 0.050)<br> AT (0, 0, 0.02) RELATIVE sample_out<br><br>/************************detector*********************************/<br>COMPONENT PSDafterSample = PSDlin_monitor(<br> nx = 200, filename = "PSDafterSample", xwidth = 0.2,<br> yheight = 0.2)<br> AT (0, 0, L3) RELATIVE sample_out<br><br> <br>/* This section is executed when the simulation ends (C code). Other */<br>/* optional sections are : SAVE */<br>FINALLY<br>%{<br>%}<br>/* The END token marks the instrument definition end */<br>END<br><br>Regards,<br>Tom<br><br><br> <br> <br></div><div> </div><div> </div><br><!-- footer --><br>
<hr>
<a style="font-size:12px;line-height:15px; color:blue; text-decoration:underline;" href="http://popme.163.com/link/004572_0721_2489.html">想免费升级vip邮箱吗?</a><font style="font-size:12px;line-height:15px;">快来参加喝汇源游网易活动吧!</font>