激光脉冲在各向异性散射介质内的瞬态热效应

本文考察了激光脉冲在吸收、发射、各向异性散射介质内引起的瞬态热效应。将激光脉冲辐射在介质内的传递过程分为两个子过程：发射－衰减－反射过程和吸收－散射过程。用光线踪迹法结合节点分析导出辐射传递系数和辐射源项,用控制容积法解瞬态能量方程。检验结果表明,本文的计算方法准确。在此基础上,考察了散射特性、初始温度对激光脉冲响应的影响。     

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

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

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

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

层状各向异性固体媒质的声反射和声透射

本文导出了联系任意对称的各向异性固体层上下界面力学量（位移和应力）的传递矩阵，从而提出了层状固体媒质声反射和声透射的一般模型，作为例子，本文数值计算了金属铝中双层正交对称固体（其间界面呈刚性联接或滑移联接）对声波的反射和透射系数。     

双向蒙特卡罗法模拟参与性介质中的辐射换热

双向蒙特卡罗法继承了正向蒙特卡罗法和反向蒙特卡罗法两种方法的优势,不仅可以解决各种各向异性散射问题,还能简易直接地解决辐射能量的方向性和定位置问题。为了验证双向蒙特卡罗法模拟热辐射传递,与文献中其它方法计算结果进行了比较。最后,利用双向蒙特卡罗法模拟计算圆柱形参与性介质中的辐射换热。     

结构各向异性对传输平均自由程的影响研究

本文研究了由定向排布的椭球形粒子组成的无序介质中微纳结构的各向异性对辐射传输平均自由程的影响机理。采用离散偶极子算法(DDA)计算单个椭球形粒子的各向异性辐射特性,为了得到多重散射下的各向异性介质的辐射特性,采用随机行走方法研究了传输平均自由程随粒子浓度、长径比、分布角度的变化关系,并与各向同性理论的计算结果进行了对比,证明了现有理论的不足。最后,利用辐射传递方程的扩散理论得到了椭球形粒子各向异性介质的表面反射能量分布图。结果表明,微纳结构的各向异性对辐射传输平均自由程和表面辐射能量分布有较大影响,对于各向异性较强的介质,在理论计算中必须加以考虑。     

漫反射各向异性散射介质内的温度热流场

本文采用射线踪迹结合节点分析法和谱带模型,研究了漫反射不透明边界下吸收、发射、各向异性散射介质内的热辐射传递过程。考虑介质辐射能的入射和散射方向,导出漫反射、不透明边界、各向异性散射介质的辐射传递系数。在辐射平衡的情况下,考察了表面发射率和散射反照率对介质内辐射热流和温度场的影响。研究表明,介质不透明边界处存在温度跃迁现象,而且,内界面发射率越大,相应界面温度跃迁越小。     

具有复杂界面的三维弹性液滴动力学模拟

具有复杂界面特性(如界面弹性、界面抗弯特性等)的弹性液滴(如血液细胞)的动力学行为对理解血流等液滴群流动极其重要。本文采用界面追踪法(Front tracking method)捕捉界面演变,界面有限元考虑界面弹性和抗弯特性,有限差分法求解控制方程,发展了研究复杂界面液滴动力学的直接数值模拟方法。通过模拟球形弹性液滴的变形运动验证了该方法的准确性,进而研究了线性剪切流场中弹性液滴的坦克履带式、摇摆式和翻滚式运动。     

基于区域分解和逐级光线分裂的杂散光分析

针对蒙特卡罗法运算量大的问题,通过引入区域分解和逐级光线分裂技术,减少杂散光线与表面交点求解时所涉及的表面数及跟踪光线数,降低计算量.区域分解时,考虑子区域与表面数的匹配关系;并基于杂散光传播过程的模拟信息,逐级确定光线分裂倍数.以大衰减比抑制结构和多传感器系统为例,进行杂散光分析.结果表明,对复杂系统杂散光传播的蒙特卡罗法模拟,应用区域分解和逐级光线分裂技术,可以在保持计算精度的前提下,有效地提高计算效率.     

Vlasov辐射器反射特性研究

 利用矩量法导出Vlasov辐射器端口处各种模式模系数（振幅）的计算方法，并利用此方法通过数值计算研究了斜切形Vlasov辐射器的反射特性，得到了辐射器端口反射系数和斜切角度及直径波长比之间的关系。结果表明斜切形辐射器具有很小的反射系数，并从反射特性上为辐射器结构参数的选择提供了依据。     

一个模拟光子照射人体的蒙特卡罗程序

本文分绍了一个用于辐射防护和医疗目的的蒙特卡罗程序。它能计算出受到光子照射时人体20种器官(组织)的吸收剂量和人体有效剂量当量。程序中使用了一种特殊的方法以加快计算结果的收敛。本文所介绍的方法也可用于其他复杂几何系统的剂量计算。     

高温堆乏燃料贮罐中子屏蔽性能计算

球床模块式高温气冷堆采用包覆颗粒球形燃料元件,在反应堆运行过程中,不断排出的乏燃料球将被装入乏燃料贮罐。乏燃料贮罐应选取适当的材料和厚度,对光子和中子进行有效屏蔽,使罐外的剂量率满足相应的限值要求。为此,使用张弛长度法和蒙特卡罗模拟法研究乏燃料贮罐的中子屏蔽性能。屏蔽材料为铁和含硼聚乙烯,计算了铁和不同B4C含量聚乙烯的屏蔽性能,并给出了乏燃料贮罐装满乏燃料球后,乏燃料球自吸收对贮罐外剂量率的影响。两种方法计算结果吻合很好,可以为实际工程中的屏蔽设计提供参考意见。     

Monte Carlo方法在单粒子翻转数值模拟中的应用

描述了在单粒子翻转数值模拟中的一种有效方法-Monte Carlo方法,并对该方法中的粒子输运过程和相关的随机抽样进行了描述.对14MeV的中子从存储器硅片表面随机入射引起的单粒子翻转进行了计算和分析,同时计算了Monte Carlo方法引起的误差.     

晶柱粘连对CsI:Na转换屏分辨特性的影响

本文用MonteCarlo方法研究了CsI:Na转换屏中晶柱之间的粘连对其空间分辨特性的影响,给出了不同粘连程度下可见光在转换屏中的点扩展函数及相应的MTF曲线.通过几种不同粘连系数下模拟结果的比较可见,要获得一个高分辨的转换屏应尽可能减小晶柱之间的粘连.对于一个实际的转换屏其粘连系数至少应在40%以下,最好控制在20%以内.     

电子加速器屏蔽室空间散射的蒙特卡罗计算

为解决电子直线加速器屏蔽室的空间散射计算相关问题,采用蒙特卡罗程序MCNP4B进行模拟,使用分步计算的方法,结合分裂和轮盘赌技巧,计算出9MeV驻波直线加速器运行时的空间散射剂量率大小.并实际测量了加速器运行时屏蔽室周围距地面上方1m处的空间散射剂量分布.将计算得到的结果、计算空间散射的经验公式结果和实际测量的结果三者进行了比较分析.结果表明:蒙特卡罗方法计算值比较符合实际测量值,经验公式计算值则存在一定的偏差.     

DRESOR法对平行入射辐射问题的研究

本采用一种基于蒙特卡洛法(Monte Carlo Method,MCM)求解辐射传递方程(Radiative Transfer Equation, RTE)的快捷、有效的方法-DRESOR法(Distributions of Ratios of Energy Scattered Or Reflected)在一维充满吸收、各向同性散射介质平行平板中,外部有平行入射条件下,求解计算空间点的辐射强度沿空间方向角的分布,而不需要辐射平衡和在空间位置坐标和方向角度坐标上同时离散辐射传递方程进行迭代求解。     

路径积分Monte Carlo方法与非线性随机振动

给出了一种求解非线性随机振动问题的新方法-基于路径积分表述的Monte Carlo方法。受白噪声激励的非线性系统的各响应统计值被表成路径积分形式,并采用Monte Carlo方法进行计算。讨论了应用方面的问题,并计算了两个实例。     

DRESOR法对瞬态辐射传递问题的研究

用基于蒙特卡洛法(Monte Carlo Method,MCM)的DRESOR法(Distributions of Ratios of Energy Scattered by the medium Or Reflected by the boundary surface)求解入射辐射经过介质散射、壁面反射传递后辐射强度随时间变化的瞬态辐射传递方程(Transient RadiativeTransfer Equation,TRTE)问题。通过在系统内计算一单位瞬态入射辐射对介质的DRESOR数分布,就能计算任意时间内入射辐射在系统内时间响应特性,这样有效提高数值方法处理瞬态辐射问题的通用性。并且能够获得高方向分辨率的辐射强度随时间变化的结果,这是目前大多数数值处理方法比较难做到的,显示出了DRESOR法处理瞬态入射辐射问题的能力．     

求解介质内热辐射传递的双向统计蒙特卡罗法

基于热辐射传输的光路可逆性原理,提出了求解介质内热辐射传递的双向统计蒙特卡罗法(BSMC法)。该方法采用等温等权抽样,利用能束传输路径的可逆关系进行辐射传递的双向统计计算,充分利用了能束跟踪的计算信息。以二维矩形区域内吸收性介质的热辐射传递为例,介绍了BSMC的求解过程,分析了其计算误差。通过数值模拟,从辐射传递因子计算结果的倒易性满足程度与辐射平衡温度场两方面,将BSMC法与传统的蒙特卡罗法(TMC法)进行了比较。结果表明,在相同的计算量下,BSMC法比TMC有更高的模拟精度。     

蒙特卡罗粒子输运问题的全局降方差方法

精确获得中子通量的时间-空间-能量精细分布对于一类非定常粒子输运蒙特卡罗模拟至关重要.以零方差理论为基础,采取化整为零的策略,把蒙特卡罗方法与离散纵标方法相耦合,用近似重要函数指导蒙特卡罗模拟,给出一种实现蒙特卡罗全局降方差的计算方案.数值算例表明,该计算方案获得了全局降方差的效果.