再入飞行器驻点流光辐射特性的理论建模

本文用直接模拟蒙特卡罗法，对再入体驻点区域的流场、光辐射特性进行了研究．针对头部锥体驻点区（包括激波层），建立了流场中的化学反应模型，中性粒子和等离子体的热运动、粒子间碰撞的理论模型，碰撞激发辐射和化学发光辐射的物理模型．给出了样本粒子运动方向和速率、粒子碰撞对的选择等蒙特卡罗法概率模型和计算方法．     

近红外光谱LASSO特征选择方法及其聚类分析应用研究

近红外光谱技术是一种通过分析样本的特征光谱数据,实现定性或定量分析的无损检测方法,特征数据的完整性和代表性决定了所建模型的性能,而现有分析方法只能实现光谱子区间特征筛选,导致分析模型稳定性差、且难以再优化。为实现近红外光谱区间高维数特征提取,有效提高近红外光谱定性分析模型的精度和稳定性,提出一种基于最小绝对收缩和选择算法(LASSO)的光谱特征筛选方法,并以我国特色高值外贸产品云南松茸为分析对象进行聚类应用研究,讨论了该方法对于高维光谱特征筛选的有效性、分析对比了LASSO筛选特征变量及主元分析(PCA)降维算法所建松茸真伪甄别及食用菌分类模型的预测精度及稳定性。通过调研发现,云南产鲜松茸因其独特外形易于分辨,而片状的干松茸失去其独有的外形特征,导致国内干松茸掺假事件屡禁不止。选取云南产松茸、杏鲍菇、老人头、姬松茸四种干样共166样本数据进行分析,采用光谱范围为900～1 700 nm的NIRQuest512型近红外光谱仪获得166×512维原始光谱数据,剔除异常数据后采用标准正态变换对光谱数据进行预处理。在此基础上,利用LASSO筛选出全光谱区间的特征变量,再使用Kennard-Stone法并结合典型线性(KNN)和非线性建模(BP)算法,构建松茸真伪甄别模型和食用菌分类模型,对两种模型进行盲样测试,并分析了LASSO与PCA算法的不同点,最后使用蒙特卡罗方法检测两种模型的稳定性。实验结果表明基于LASSO光谱特征选择的松茸真伪甄别模型和食用菌分类模型预测精度和稳定性均高于PCA方法,其中基于原始光谱数据所建真伪甄别模型的预测准确率为69.57% (BP)和60.87% (KNN),食用菌分类模型准确率为67.39% (BP)和65.22% (KNN),基于LASSO特征筛选的真伪甄别模型预测准确率分别达到100% (BP)和78.26% (KNN),食用菌分类模型预测准确率分别达到89.13% (BP)和80.43% (KNN),对两种模型进行10次蒙特卡罗实验,其结果平均值分别为99.93%和97.22%,由此可知,与PCA等数据降维算法相比,LASSO可实现全光谱区间的光谱特征选择和数据降维,有效地提高了近红外定性分析模型的预测性能,为近红外分析提供了一种新的特征筛选方法。     

上证50ETF期权定价有效性的研究:基于B-S-M模型和蒙特卡罗模拟

上证50ETF期权作为中国资本市场上股票期权的第一个试点产品,其定价问题尤为重要。本文分别运用B-S-M期权定价模型和蒙特卡罗模拟方法对其定价进行实证研究,分析结果表明:1)IGARCH模型比传统的GARCH模型更能较好地拟合上证50ETF的波动率;2)当模拟次数为1000时,蒙特卡罗方法的效率一致地高于B-S-M模型,并且除了对偶变量技术的拟蒙特卡罗其他模型的精确度也都高于B-S-M模型;3)B-S-M模型和蒙特卡罗模拟方法都可以较为准确地、有效地模拟出上证50ETF期权价格。这些研究将为今后期权定价模型的发展和完善提供必要的参考和指引。     

一种针对弱场搜寻态冷极性分子的可控静电表面囚禁方案

提出一种使用带电金属环和六个球电极和一个外加偏置电场实现对冷极性分子静电囚禁的新方案.计算装载和囚禁时空间电场分布.囚禁中心距离芯片表面的高度可以通过外电场和环形电极所加电压来操控.蒙特卡罗模拟表明对于中心速度为15 m·s^-1的ND₃分子束,装载效率可以达到70%,得到冷分子的温度大约为45 mK.当继续增加偏置电场强度时,单阱分裂为对称的两个阱.如果同时改变球形电极上所加电压,得到不对称的两个阱,可以借此来调节两个阱中所囚禁的冷分子数目的比例.为了易于理解,用蒙特卡罗方法模拟了装载、囚禁、分裂冷分子波包的动力学过程.     

γ射线探测器响应函数的JMCT模拟计算

采用自主研发的JMCT软件,模拟计算了射线在NaI闪烁探测器中的响应函数。JMCT模拟计算了探测过程中所有的、完整的光子-电子耦合输运过程。其模拟计算的射线能量沉积谱、探测器响应函数,与MCNP的计算结果符合得非常好,与Berger解析法计算的结果也基本符合,表明JMCT软件可以在探测器响应函数方面得到满意的结果,JMCT将在实验核物理、核分析技术应用等方面发挥越来越重要的作用。     

蒙特卡罗模拟分析电子束发射度对照射量空间分布影响

加速器靶前正前方1 m处的照射量是衡量加速器光源辐射能力的重要物理参量。电子束轰击高原子序数靶产生的X射线空间分布具有很强的前冲性,照射量空间分布与电子束的发射度密切相关。采用高斯函数对不同电子束发射度条件下入射电子的空间分布和角度分布进行建模。应用蒙特卡罗方法对电子束打靶的轫致辐射过程进行模拟,分析电子束发射度对照射量的影响。同时还对整靶结构和叠靶结构下的轫致辐射光源照射量进行计算比较。结果表明,电子束的发散角是影响轫致辐射光源照射量的主要因素。与采用整靶结构相比,采用叠靶结构所获得的照射量空间分布基本一致,在正前方小角度范围内（0~4）的照射量有2%~3%的降低。     

辐射成像系统模拟软件NUCRPD的研发与应用

在辐射成像系统的设计研发中,经常需要计算系统的各项物理性能指标。为了提升产品研发的效率,同方威视联合清华大学,基于蒙特卡罗程序包Geant4,研发了国内首套具有自主知识产权的辐射成像系统模拟软件NucRPD(NUCTECH Radiography Performance Design Tools),其能够快速、精确地模拟系统的各项性能指标。该软件对常用的一些辐射成像系统进行了参数化建模,用户通过修改软件界面上的少量参数,就可以快速建立各类辐射成像系统的几何体、源项、物理模型和统计量,然后在服务器上以并行计算方式完成模拟计算,能在较短的时间内模拟出系统的性能指标,并给出直观的图形帮助用户深入理解模拟结果。NucRPD的模拟结果经过了大量的实验验证,其剂量场分布和物理指标等模拟结果和实验结果符合得很好。NucRPD已经应用于同方威视辐射成像系统产品的设计研发,在产品的物理指标优化和辐射防护优化中发挥了重要作用。     

JAEA 10A���Դ�ŵ��е缫����ģ��

采用自主研发的全三维PIC/MCC软件对JAEA 10A多峰负氢离子源进行系统仿真,分析了电极布局对H?体积产生效率的影响.模拟结果显示:放电功率一定的情况下,增加电极数目,高能电子导致的横向漂移相对减弱,进而得到相对均匀的H?束;但电子的相对利用率越低,导致H?体积产生效率较低.电极越靠近过滤磁区,高能电子导致的横向漂移相对增强,H?体积产生效率较低.改变电极数量,H?体积产生效率和空间均匀性不可兼得,电极远离磁过滤区,可实现双赢.     

外力驱动作用下高分子链在表面吸附性质的计算机模拟

采用退火法模拟研究受外力F驱动的高分子链在吸引表面的吸附特性.通过高分子链的平均表面接触数〈M〉与温度T之间的关系计算临界吸附温度T_c,并发现T_c随着F的增加而减小;进而通过高分子链的均方回转半径分析外力驱动作用对高分子链构象的影响,并从回转半径极小值或者垂直外力方向的y和z分量的变化交叉校验临界吸附点T_c.模拟计算了处于吸附状态的高分子链随着外力F的增加是否会发生吸附状态到脱附状态的相变以及发生相变所需施加的外力是否由温度所决定.模拟结果表明:两种不同温度下高分子链的吸附性质和构象性质受外力驱动作用而产生不同现象,在温度区间T*_cTT_c时会发生脱附现象,而在TT*_c时不会发生脱附现象.     

材料二次电子产额对腔体双边二次电子倍增的影响

采用蒙特卡罗抽样与粒子模拟相结合的方法,数值研究了材料二次电子产额对腔体双边二次电子倍增瞬态演化及饱和特性的影响.研究发现:随着材料二次电子产额的增加,二次电子增长率以及稳态二次电子数目和振幅均呈现增加的趋势,放电电流起振时间逐步缩短,稳态电流幅值以及放电功率平均值和振幅值均呈现逐步增加并趋于饱和的规律,沉积功率波形延时以及脉宽呈现逐步增加并趋于饱和的趋势.粒子模拟给出了高/低二次电子产额情况下的电子相空间分布、电荷密度分布、平均碰撞能量、平均二次电子产额、二次电子数目和放电电流的细致物理图像.模拟结果表明:高二次电子产额材料,饱和时更倾向趋于单边二次电子倍增类型分布;低二次电子产额材料的二次电子倍增饱和特性由空间电荷场的"去群聚"效应和"反场"效应同时决定,而高二次电子产额材料的二次电子倍增饱和特性则主要是由发射面附近的强空间电荷场"反场"效应决定的.