241Am-Be中子参考辐射的蒙特卡罗模拟研究

研究用于校准场所中子剂量监测仪表的241Am-Be中子参考辐射场计量特性。采用蒙特卡罗方法模拟了空气自由中子参考辐射(FRNR),GB/T 14055规定的最小尺寸中子参考辐射(SRNR)和实际中子参考辐射(ARNR)中不同检验点处中子周围剂量当量率、散射中子占比和能谱分布特征。研究结果表明,空气对FRNR中的剂量率和能谱分布影响小,近似为理想中子参考辐射;采用5%含硼聚乙烯作屏蔽的最小尺寸SRNR可减少热中子,降低散射中子占比,影锥法不适用于小尺寸中子参考辐射中对散射中子的修正;ARNR中的散射中子更少、占比更低,影锥法所得散射中子占比与理论值基本一致。     

均匀磁场下碳离子笔形束的剂量变化分析及位置修正方法

分析外加均匀磁场对于碳离子笔形束剂量分布的影响,并考虑修正这种影响,为磁共振成像引导碳离子放射治疗的临床应用提供指导。本文利用蒙特卡罗方法模拟计算了不同能量碳离子笔形束在不同强度磁场下的剂量分布情况,发现垂直于碳离子束入射方向的均匀磁场对于碳离子笔形束射程缩短的影响很小,磁场对碳离子束的主要影响是引起束流横向偏转,特别是碳离子束布拉格峰位置的横向侧移。横向侧移程度与碳离子束的能量和磁场强度相关,根据模拟结果,得到了一个计算碳离子束布拉格峰在磁场中相对横向偏转的方程,并提出一种校正外加磁场引起的碳离子束布拉格峰横移的角度修正方法。这些结果可用于指导磁共振图像引导碳离子放射治疗计划系统的研发。     

80.5 MeV/u碳离子诱发铜靶的放射性剩余产物测量

高能重带电粒子能直接穿透靶原子核外电子层,与原子核发生直接碰撞,发生散裂反应,产生一系列具有放射性的剩余产物核.重带电粒子诱发靶材放射性剩余核与辐射防护和人员安全有着密切联系,当前,大部分剩余核产额主要依靠蒙特卡罗粒子输运程序进行模拟计算,其准确程度亟需通过实验测量进行准确评估.本文利用能量为80.5 MeV/u的(12)^C⁶⁺粒子对薄铜靶开展了辐照实验与伽玛射线测量,结合伽玛谱学分析方法,得出了辐照产生的18种放射性剩余产物的初始活度和产生截面值,并与PHITS模拟结果进行对比.结果表明,PHITS模拟程序对放射性剩余核种类的估计具有较高可靠性,在其绝对产额方面,与实验测量仍具有较大偏差.     

基于蒙特卡罗法的水下激光光幕探测性能研究

为了提高水下激光光幕的探测性能, 根据水下光束传播规律, 构建了水下激光光幕探测模型, 基于蒙特卡罗模拟方法对水下激光光幕探测性能进行研究。基于水下激光光幕探测模型, 利用MATLAB软件进行仿真, 分析海水衰减系数、初始功率及传输距离对水下激光光幕传输的影响。仿真结果表明:海水的衰减系数越小, 水下激光光幕传输率受到的影响越小。海水参数和传输距离一定时, 随着初始功率的增加, 只会影响到达探测端的最终功率, 但对传输率影响不大。当海水衰减系数一定, 传输距离为1 m时, 其传输率约为15%且变化稳定; 当传输距离增加到30 m时, 传输率在5%以下。     

质子辐照硼硅酸盐玻璃盖片的物理效应分析

作为航天器电源系统的重要组成部分,太阳电池需要更高的转换效率和可靠性以及更长的使用寿命。通过在太阳电池表面覆盖抗辐照玻璃盖片,可以增强太阳电池对粒子辐射的防护,延长太阳电池的服役寿命,使航天器获得可靠的能源供应。硼硅酸盐玻璃就是一种理想的太阳电池玻璃盖片材料。采用蒙特卡罗方法,结合SRIM软件模拟研究质子辐照硼硅酸盐玻璃的损伤物理机理。基于粒子与物质相互作用的理论以及基本公式,通过分析不同入射能量的质子在硼硅酸盐玻璃中的阻止本领、电离能损、位移能损、空位的产生情况,对辐照损伤的物理机制进行研究。结果表明:能量为30~120 keV的质子辐照损伤主要发生在硼硅酸盐玻璃表面;质子沉积、空位分布等均为Bragg峰型分布;电离能损是能量损失的主要部分,随入射能量的增加而增大,导致电子的电离和激发;位移能损在玻璃内部随能量降低而增大,导致硼、氧和硅等空位缺陷的产生;电离效应和缺陷的产生是硼硅酸盐玻璃色心形成的重要原因。     

用人工神经网络方法测量生物组织光学参数γ

提出了一种利用亚扩散空间分辨漫反射预测生物组织约化散射系数μ_s'和相函数参量γ的人工神经网络方法 .采用蒙特卡罗方法得到光子经生物组织漫反射的数据样本,利用这些数据样本训练反向传播神经网络,用于从亚扩散散射光中预测γ的信息.为了解决同时预测μ_s'和γ两个参数时会产生较大误差的问题,分段数据训练两个BP网络,依次识别μ_s'和γ.研究发现3.64l_(th)(l_(th)表示平均输运自由程)是对γ的不敏感点,可用该点附近的数据样本训练网络用于预测μ_s',用2l_(th)范围内的数据样本训练网络用于预测γ.蒙特卡罗仿真结果表明,在1.3≤γ≤1.9范围内,预测结果与真实值的相对均方根误差在1%以内.与现有的测量方法相比,所提的人工神经网络方法更加简单,且提高了预测精度.     

碳离子在不同材料叶片的多叶光栅上产生的次级粒子研究

在碳离子放射治疗中，碳离子束与治疗头设备和患者身体相互作用产生的次级粒子可以到达患者体内的许多区域，在产生的次级粒子中以中子和$\gamma $射线的产额为最大。在不影响束流配送功能的情况下，减少碳离子放疗中产生的次级中子和$\gamma $射线对于降低放疗后出现的正常组织并发症及二次肿瘤风险有着非常重要的意义。本文利用蒙特卡罗(Monte Carlo)方法模拟计算了被动式束流配送系统下，400 MeV/u碳离子束照射到由不同材料叶片构成的多叶光栅(MLC)形成典型的10 cm×10 cm方形射野时，在水模体中产生的次级中子和$\gamma $射线所沉积的剂量及空间分布等。模拟结果显示:碳离子束通过MLC形成射野后在水模体中产生的次级中子主要分布在水模体的入射端，次级$\gamma $射线较为均匀的分布在整个水模体内，且较多分布在具有展宽Bragg峰(SOBP)射野在水模体中贯穿时的坪区。对于MLC叶片材料的选择，则需根据实际情况对叶片厚度以及次级粒子当量剂量的要求来确定。本文通过模拟研究不同MLC叶片材料产生次级粒子的情况，为被动式束流配送系统中MLC叶片及其他元件的材料选择提供了科学依据。     

商品和金融摇摆期权定价模型及其数值研究

采用具有均值回复特点的三叉树模型研究上摇一下摇特征或带有惩罚条款的复杂商品摇摆期权,并运用森林树法进行数值计算.进而分析最大可执行次数、每次执行最大可执行量、惩罚系数对摇摆期权价格的影响.此外还运用多次最优停时研究连续时间框架下的金融摇摆期权定价问题.为了避免高维诅咒,选择改进的最小二乘蒙特卡罗方法完成数值计算,并分析标的资产价格及最大可执行次数对期权价格的影响.数值分析结果表明,当实施次数为一次时,摇摆期权价格与相应的美式期权价格重合,从而在一定程度上验证了模型的正确性.     

迭代加权质心法机理及多星定位误差特性研究

提高星敏感器星点定位精度是高精度星敏感器关键技术之一。常规质心法和高斯拟合法定位最小均方误差在0.1 pixel左右,难以满足未来高精度定位需求。在自适应光学领域,为研究Shack-Hartmann传感器波前定位而提出的迭代加权质心算法(IWCOG),其定位误差优于前两种算法,但目前该算法的定位机理及多星定位误差特性未知。从Meanshift理论角度推出该算法定位机制,证明算法的收敛性;每个仿真采用蒙特卡罗方法随机生成10000个星点样本,通过分析星点参数、信噪比与定位误差的关系研究IWCOG算法在星敏感器多星定位的特点。模拟结果显示,在信噪比R=40下,IWCOG法星点提取误差低于0.02 pixel的样本占44%,低于0.04 pixel的占72%,远高于质心法的0.8%和1.9%,基本满足甚高精度星敏感器的定位要求。     

反应堆蒙特卡罗临界模拟中均匀裂变源算法的改进

在反应堆pin-by-pin精细建模及蒙特卡罗模拟计算研究中, 由于不同栅元的功率密度差异较大, 导致蒙特卡罗方法临界计算的样本在不同栅元之间的分配不均衡, 由此引起栅元内的各种计数的统计误差差异较大. 为使大部分栅元内计数的统计误差降至一个合理的水平, 单纯增加总样本已不是一个高效的解决方法. 通过在特定临界计算迭代算法的基础上改进并实现均匀裂变源算法的思想, 对大亚湾压水堆pin-by-pin模型取得了具有较高效率的数值结果. 本工作为具有自主知识产权的蒙特卡罗粒子输运模拟软件JMCT最终达到反应堆pin-by-pin模型(包括一系列国际基准模型)的模拟性能要求提供了一个有效的工具.