离子输运蒙特卡罗模拟的步长选取

介绍带电离子在物质中输运模拟的详细历史法,开发模拟程序ITR(ions transport and reaction),同时实现抽取离子飞行步长的基本方法和脉冲近似法,计算结果与CORTEO及TRIM程序的结果比较符合很好.讨论两种步长选取方法,对脉冲近似法的最大步长进行限制,分析CORTEO程序与TRIM程序计算结果差别的主要原因.ITR同时具有TRIM和CORTEO的优点,采用脉冲近似法减少需要计算的核碰撞次数,通过索引(Index)方法对散射角的计算进行加速,计算效率大幅提高.     

ECR微波等离子体源离子输运的蒙特卡罗模拟

采用蒙特卡罗方法，模拟了ＥＣＲ微波等离子体源中离子离开放电室后历经中性区，鞘层区，最后被加负偏压的工件表面吸收的全过程。研究了鞘层外界处的衔接问题，得到了自洽电势分布和鞘层区不同位置处的速度分布，能量分布与角分布。同时，讨论了磁场对鞘势的影响。     

离子原子碰撞过程中反冲离子飞行时间谱的蒙特卡罗模拟研究

基于蒙特卡罗方法,建立了程序模拟离子与原子碰撞中的反冲离子飞行时间谱.模拟入射离子束的空间密度分布,原子束密度的空间分布,靶原子的初速度分布, 根据反应几率产生不同电荷态反冲离子,并求出反冲离子飞行时间后,对每次模拟事件根据飞行时间进行累计,就得到反冲离子的飞行时间谱.模拟结果与实验(Ar12+-Ar) 测量到的TOF谱进行了分析比较,并进行了定性讨论     

盐离子对蛋白质带电特性的影响

本文应用分子理论,研究盐离子对蛋白质带电特性的影响,理论模型考虑蛋白质与阴离子的结合作用。研究发现,由于蛋白质与阴离子的结合,距离蛋白质表面附近处的阴离子被吸附在了蛋白质表面,在距离蛋白质表面附近区域,阴离子分布较少。通过计算体系中的静电势,我们发现,在距离蛋白质表面附近,静电势呈现了较大的负值,带正电荷的阳离子感受到静电吸引,会出现在距离蛋白质表面附近的区域,这会使得在距离蛋白质表面附近的区域,阳离子数目增多。这样,在不同阴离子浓度、以及阴离子与蛋白质不同结合能条件下,阴离子会在不同程度上影响蛋白质的带电特性、影响体系中的静电特性。通过考察不同结合能条件下,蛋白质表面电荷面密度随阴离子浓度的变化关系还发现,较大的结合能会使得阴离子与蛋白质结合增快,蛋白质表面会呈现从正电荷态向负电荷态的转变。理论结果符合实验观测,由此表明,盐离子与蛋白质的结合导致蛋白质表面带电特性的改变,是盐离子影响蛋白质带电特性的本质。     

微通道板离子壁垒膜粒子阻透特性的蒙特卡罗模拟

在三代微光像管中,微通道板(MCP)输入面上覆盖一层超薄离子壁垒膜(IBF),目的是保护光电阴极,延长像管使用寿命.为了深入研究离子壁垒膜的特性,本文对Al2O3和SiO2两种离子壁垒膜的粒子阻透能力进行了蒙特卡罗模拟,结果表明:5 nm厚Al2O3和SiO2离子壁垒膜的死电压分别为230 ～240 V和220～230...     

ECR微波等离子体离子输运的数值模拟

建立了ECR微波等离子体源离子输运的平板和圆柱模型,对离子历经的空间区域的输运过程进行了数值研究。采用Monte Carlo(M-C)方法模拟了存在外磁场情况下,离子离开开放电室后历经中性区、鞘层区、最后被加负偏压的工作表面吸收的全过程,考虑了离子与中性粒子的电荷交换碰撞玫弱性散性,统一处理了中性区和鞘层区电势的衔接,采用曲线拟合,电势自洽迭代方法把中性区和鞘层区衔接起来,得到了光滑自治的电势分布曲线和鞘层区不同位置处的速度分布、能量分布及角分布。     

气体放电阴极鞘层氩离子的自洽蒙特卡罗模拟

建立了气体放电阴极鞘层离子的自洽蒙特卡罗模拟模型,考虑了离子与中性原子的电荷交换碰撞和弹性散射,用到了精确依赖于离子能量的电荷交换和动量输运截面。模拟了氩离子在阴极鞘层中的运动,得到了不同气压下自治电场分布,离子的能量分布和角分布,发现:离子由鞘层边界向阴极运动过程中,离子能量分布的高能部分逐渐增大,角分布向小角度部分压缩,鞘层中的强电场加速和聚焦了离子;在鞘层边界附近的电场呈非线性。     

基于蒙特卡罗模拟的射束硬化校正方法

传统的射束硬化校正方法, 通常需要针对每一种材料测量该材料对射线的吸收曲线. 由于吸收曲线对实验条件有很大的依赖性, 每当改变X光机电压或者被测工件的材料等条件时,需要重新测量吸收曲线才能完成硬化校正过程. 这种方法费事费时. 本文提出了基于蒙特卡罗模拟计算物质吸收曲线的硬化校正方法. 实验中, 分别用本方法和传统的硬化校正方法对铝工件进行硬化校正, 经过比较, 确认本方法是有效的. 然后用该方法对不同材料(铝、铁和铜)的工件进行校正. 实验结果表明, 本方法能有效消除各种材料工件图像中的硬化伪影, 是快速的和切实可行的.     

蒙特卡罗方法对各向异性介质辐射特性的模拟

本文提出一种计算含微粒半透明介质辐射特性的蒙特卡罗方法，考虑了微粒群的独立多次散射，直接用米氏相函数考虑微粒的各向异性性质，计算了含微粒层的双向反射率与双向透射率。与实验数据和离散座标法的计算结果比较，吻合良好。     

相对论高次谐波的蒙特卡罗模拟

根据相对论条件下的经典电磁辐射理论,采用蒙特卡罗计算方法对类氢离子在强激光场下辐射的高次谐波进行了模拟计算,获得了相对论领域内激光参数（频率、脉冲波形、强度等）对高次谐波影响的各种计算结果,并对结果进行了理论分析。     

相对论高次谐波的蒙特卡罗模拟

 根据相对论条件下的经典电磁辐射理论,采用蒙特卡罗计算方法对类氢离子在强激光场下辐射的高次谐波进行了模拟计算,获得了相对论领域内激光参数（频率、脉冲波形、强度等）对高次谐波影响的各种计算结果,并对结果进行了理论分析。     

蒙特卡罗模拟中相关变量随机数序列的产生方法

蒙特卡罗模拟有时需要对多维相关随机变量进行模拟抽样．本文介绍基于Choesky因子线性变换一非线性变换产生具有指定边缘分布和相关系数的多维相关随机变量抽样序列的一般方法,给出一种简单易行的高效数值实现途径和一些模拟结果．模拟结果表明,该方法产生的各随机变量抽样序列间具有预期要求的相关性,并能通过指定边缘分布Kolmogorov—Smimov非参数假设检验．对该方法应用中的一些限制问题进行了讨论．     

接触热阻的蒙特卡罗法模拟

 用蒙特卡罗随机模拟法计算了固体接触界面间的接触热阻。取服从高斯分布的随机数作为粗糙峰高度，分别对每个粗糙峰计算出单点接触热阻，再并联所有单点接触热阻得到总的接触热阻。修正了“积分法”的理论结果且与发表的文献数据相符。     

热辐射输运问题的高效蒙特卡罗模拟方法

惯性约束聚变研究中,热辐射光子在介质中的输运以及热辐射光子与介质的相互作用是重要研究课题,蒙特卡罗方法是该类问题的重要研究手段之一.隐式蒙特卡罗方法虽然能正确地模拟热辐射在介质中的输运过程,但当模拟重介质(材料的吸收系数大)问题时,该方法花费的计算时间将变得很长,导致模拟效率很低.本文以离散扩散蒙特卡罗方法为基础,开发了"离散扩散蒙特卡罗方法辐射输运模拟程序",可以较好地解决重介质区的计算效率问题,但是离散扩散蒙卡罗方法在模拟轻介质区时精度不够高.辐射输运问题中通常既有轻介质也有重介质,为了能同时解决蒙特卡罗方法模拟的效率和精度问题,本文研究了离散扩散蒙特卡罗方法与隐式蒙特卡罗方法相结合的模拟方法,并提出了新的扩散区与输运区界面处理方法,研制了混合蒙特卡罗方法的辐射输运模拟程序.典型辐射输运问题模拟显示:在模拟重介质问题时,该程序能大幅缩短模拟时间,且能取得与隐式蒙特卡罗方法一致的结果;在模拟轻重介质均存在的问题时,与隐式蒙特卡罗方法相比,混合蒙特卡罗方法的模拟精度与其相当且计算效率同样能够得到显著提升.     

混响室的概率统计分析方法及其蒙特卡罗模拟

提出了一种研究混响室问题的理论分析方法,从理想混响室的离散概率模型出发,采用蒙特卡罗模拟来获得室内电磁场量及其相关物理量的概率统计特性。运用该方法获得了混响室内若干场量（均为归一化值）的概率密度函数、室内一段均匀无耗单导体传输线的负载平均耦合截面及负载电流响应模值的概率分布,所得结果与已有文献给出的结果（包括由解析公式计算或测量所得的结果）吻合良好。由于蒙特卡罗方法的普适性,采用该方法可以获得任意待分析参量的全部概率统计特性（只要样本数足够大）。同时该方法无需针对具体混响室进行,只要在其测试区域能较好地与理想混响室的基本概率模型相吻合,就能给出较精确的计算结果,且其计算量远远小于基于全波分析的蒙特卡罗方法。     

混响室的概率统计分析方法及其蒙特卡罗模拟

 提出了一种研究混响室问题的理论分析方法,从理想混响室的离散概率模型出发,采用蒙特卡罗模拟来获得室内电磁场量及其相关物理量的概率统计特性。运用该方法获得了混响室内若干场量（均为归一化值）的概率密度函数、室内一段均匀无耗单导体传输线的负载平均耦合截面及负载电流响应模值的概率分布,所得结果与已有文献给出的结果（包括由解析公式计算或测量所得的结果）吻合良好。由于蒙特卡罗方法的普适性,采用该方法可以获得任意待分析参量的全部概率统计特性（只要样本数足够大）。同时该方法无需针对具体混响室进行,只要在其测试区域能较好地与理想混响室的基本概率模型相吻合,就能给出较精确的计算结果,且其计算量远远小于基于全波分析的蒙特卡罗方法。     

直接模拟蒙特卡罗法对连续流体传热和流动的模拟

直接模拟蒙特卡罗法在稀薄气体模拟中获得了成功的应用。然而,在处理连续介质传热和流动问题时,模拟的速度极大地制约了该方法的应用。尤其对于大系统的连续流动,该方法的收敛速度几乎无法实现。鉴于此,本文通过引入超粒子,对直接模拟蒙特卡罗方法进行改进,并将其应用于连续介质传热和流动模拟之中。