采用离散坐标法计算一维平行入射辐射问题

针对离散坐标法难以应用于平行入射辐射的问题，提出了一种简单有效的方法；在原有的离散坐标系列中，直接添加一个与平行入射辐射的方向重合的离散方向，并将该离散方向的加权取得极其微小（合如1）10^-200)，这一新的离散方向对于离散坐标系列的各阶段矩方程没有任何影响，因此可以任意地设置，这样无论平行入射辐射来自哪个方向，都可以被准确地捕捉到，并采用一个一维大气辐射问题对本方法进行了检验，计算结果与精确解吻合。     

蒙特卡罗方法用于研究堆的屏蔽计算

为解决研究堆的屏蔽计算问题 ,讨论了应用 MonteCarlo方法的 MCNP程序进行反应堆的屏蔽计算 ,引入了结合抽样技术的 MCNP分步计算的方法。该方法特别适用于计算物理模型比较大 ,屏蔽层比较厚 ,粒子的穿透距离比较长的问题。在这一类的问题中 ,单步计算根本无法计算出结果 ,利用分步计算和抽样技术则能在较短的时间内以较高的精度 (统计误差小于 10 % )得出结果。结合抽样技术的MCNP分步计算方法能够很好地解决深穿透问题     

用离散坐标法计算任意方向辐射强度(Ⅱ):验证

为验证提出的离散坐标算法，采用了一个带内热源的平行平板测试问题和两个二维吸收—各向同性散射测试问题．数值试验表明本方法能够以较少的离散方向和较短的计算时间来获得具有比较精确的任意方向上的辐射强度．     

一种可随意调节的旋转对称离散坐标格式

提出了一种既可随意布置离散方向 ,又能满足 90°旋转对称要求的新的离散坐标格式 ,其基本思路是 :先将每个卦限中的单位球表面划分为三个相等的部分 ,再在其中一部分内根据需要任意地布置离散方向 ,然后根据对称性得到其他两部分内的设置 .数值试验表明本格式具有良好的精度     

采用蒙特卡罗法分析离散坐标法的假散射

采用蒙特卡罗法获得了离散坐标方程的不含任何假散射的高精度解，其基本思路是：当采用一个离散坐标格式来计算一个漫射表面的热辐射时，就意味着用有限个离散方向去“代表”2π立体空间上的无穷多个方向，因此，不妨假定存在着这样一个虚拟表现：该表面的确只沿着该离散坐标格式的离散方向上发射热辐射，并且在每个方向上的热流也完全遵循离散坐标方程，然后将蒙特卡罗法应用于此虚拟表现，由此所获得解即为此有面的高精度解，不言而喻它也等效于离散坐标方程的高精度解，在此基础上，分析了离散坐标法的假散射对计算结果的影响。     

用离散坐标法计算任意方向辐射强度(Ⅰ):原理

针对目前所提出的各种离散坐标格式中的离散方向大都是固定的，不能任意调节，因而无法满足一些工程辐射问题的需要的问题，在“带无穷小加权离散坐标格式”的基础上，提出了一种离散坐标方法来计算任意方向上的辐射强度．同时，提出了带微小控制角有限体积法的概念．     

用于离散坐标法的一种局部分解菱形格式

提出了一种二阶精度的局部分解菱形空间差分格式，按照辐射在控制容积中的传播方向，将控制容积分解为三个部分，然后将辐射传递方程分别在各部分积分，同时沿辐射的传播方向应用菱形格式来推导出各个部分的下游边界强度，整个控制容积下游界面强度等于这三个部分的下游界强度的算术平均值，其他数值步骤与标准菱形差分格式完全相同，数值试验表明本方法完全消除了数值解的振荡，同时还可以避免负强度的发生。     

贴体坐标系下离散坐标法计算尾喷焰辐射特性

在计算尾喷焰燃气组分、辐射物性的基础上，用贴体坐标系下的离散坐标法(DOM)，研究了不同工况下尾喷焰的红外光谱辐射特性随飞行高度的变化特征，分析了介质的辐射物性对尾喷焰辐射特性的影响．研究结果表明：尾喷焰的辐射物性对尾喷焰辐射传输计算有很大的影响；在所考虑的工况下，与喷焰轴向全立的方向上光谱辐射最强。     

坐标法计算面积在储量计算中的应用

介绍了储量计算中几种面积计算的方法 ,阐述了坐标法计算面积的原理及其在矿井储量计算中的应用以及该方法计算储量程序的设计要点     

坐标法计算面积在诸量计算中的应用

介绍了诸量计算中几种面积计算的方法，阐述了坐标计算面积的原理及其在矿井储量计算中的应用以及该方法计算储量程序的设计要点。     

WSGGM结合离散坐标法研究非灰气体辐射换热问题

推导了描述非灰气体辐射特性的灰气体加权和模型(WSGGM)与离散坐标法结合的公式,通过数值模拟以及与段法的比较,表明在WSGGM参数相同的情况下,所推导的方法与段法的结果相吻合,这证明将WSGGM引入离散坐标法能成功地用于模拟非灰气体辐射问题·不同WSGGM参数对计算结果和计算时间都有影响,可根据具体场合和精度要求选择使用WSGGM参数·该法容易实现辐射场与流场的耦合求解,因此是研究非灰气体辐射换热的一种实用的方法·     

三维离散纵标中子输运稳态及瞬态动力学计算

基于著名输运程序DOT4．2，开发了三维离散纵标中子输运稳态及瞬态动力学计算程序．对时间变量采用直接的无条件稳定且具有高精度的全隐式向后差分格式进行离散处理，将瞬态输运方程转化为各个离散时间步上的固定源型稳态中子输运方程进行求解．通过对TAKEDA基准题及输运瞬态基准题的校验计算表明：该三维稳态程序能准确地给出有效增殖系数以及各区的区域平均中子通量密度，其相对误差分别小于0．1％和4．0%；对于缓发超临界和瞬发超临界问题，三维瞬态程序的计算结果均与参考值吻合良好．本三维程序可以为复杂反应堆堆芯的临界计算和瞬态动力学特性分析提供一个有效而准确的分析工具和手段．     

三维中子输运方程的非结构网格离散纵标数值解法

从一阶三维中子输运方程出发，对方向变量采用离散纵标方法展开，得到一系列关于空间变量的偏微分方程，从而避免了二阶方程由于分母上存在截面，不能准确描述内含真空介质的问题．对这些关于空间变量的方程采用最小二乘有限元方法进行离散，形成的刚度矩阵是对称的，因此可以采用快速迭代方法求解．据此编制了三维中子输运方程的非结构网格离散纵标计算程序，并采用三棱柱元素和四面体元素剖分对一系列基准问题做了验算．计算结果表明，该方法能用于非结构网格，并具有较高的计算精度，对多数问题，有效增值系数的误差都小于0．3％，通量误差都小于3．0％．     

快中子发生器屏蔽的Monte－Carlo估计

本文介绍的是一个在微机上实现的计算快中子屏蔽的Ｍｏｎｔｅ－Ｃａｒｌｏ程序．采用历次飞行法对结果进行统计估计，对中子发生器周围的屏蔽作了简单模拟，结果表明采用的物理模型是正确的．该程序改进后也可用于多次散射的计算，或模拟其他中子源的屏蔽，故有一定的实用价值．     

离散坐标法积分方案的选取

以辐射平衡条件下一维平行平板介质层的辐射换热问题为例,介绍了求解辐射传递方程的离散坐标法,对影响离散坐标法计算精度与速度的角度积分方案的选取进行了研究.分析了离散坐标法常用的SN积分方案、Fiveland等权值积分方案FN及高斯积分方案GN对计算结果准确性的影响.研究表明:积分方案近似阶数越高,计算结果的精度越高,但随着积分方案阶数的增加,所需的计算量也增大;同一近似阶数的各种积分方案的准确性各不相同,其中FN积分方案在相同阶数下的精度最高,SN积分方案次之,GN积分方案的精度最差.在近似阶数足够高时上述三种积分方案均可达到较高的计算精度.     

散体单元法确定计算时步模拟刚度法

 散体单元法DEM的核心内容是运用所谓的中心差分法,实际上使用了迭代计算,在计算中确定合理的计算时步及其法向刚度是研究问题非常重要的方面,时步的选取与被研究颗粒的刚度、质量、阻尼,以及它们之间的匹配关系密切相关。目前,确定刚度的一般做法是实测,但实测只能对特定条件进行,且需要研究人员具有较丰富的经验,才能保证计算结果接近实际。本文通过质量-弹簧系统,对DEM计算过程进行理论与数值模拟讨论,给出确定计算时步和法向刚度的选取方法,对同类计算具有一定的参考意义。     

深部巷道围岩蠕变分析与离散元模拟

将顶管法应用于深部煤矿地质条件下深部岩体施工中,会遇到因岩石流变等大变形导致顶管机的卡停甚至抱死,因此充分考虑岩石的流变特性,对顶管法在高地应力区的应用十分重要。以张集矿北区1413A高抽巷为例,利用离散元PFC3D(particle flow code 3D)程序对巷道蠕变特性进行了系统的研究。首先,调用伯格斯模型,通过虚拟蠕变试验对巷道所处粉砂岩层的微观参数进行标定,并与室内试验数据对比,从而验证离散元方法用于岩石蠕变问题分析的可行性。然后,建立巷道的平面应变离散元模型,设置多个测量球,得到了巷道围岩的应力应变在时间和空间上的分布情况,并与理论解对比,发现两者吻合较好,可为后续巷道围岩蠕变的深入研究提供有益的参考。     

同步辐射装置辐射屏蔽计算的探讨

从第三代同步辐射装置的特点出发 ,提出局部屏蔽的重要性 ,并以光束线前端为例 ,用 M- C方法模拟了 1Ge V的电子在其中的剂量分布 .     

图算角变力矩法

力矩分配法,是工程界常用的手算结构的力学分析方法.以力矩表达结构结点角位移的量,称为角变力矩.以角变力矩为未知量,以计算图为对象利用力矩分配法概念求解超静定结构的方法,称为图算角变力矩法,即边看图边计算角变力矩的方法.图算角变力矩法,克服了力矩分配法直接待求量个数多、计算的渐近性和有截断误差的近似性等缺点,发扬了力矩分配法计算程式规范,省脑力、思路清晰等优点,成为一种直接精确的力学分析方法,为建立超静定结构普遍统一的力学分析方法奠定了基础.     

基于动态元件离散化的柔性直流电网故障电流计算方法

故障电流计算是柔性直流电网设计与规划、故障类型判断的重要环节,开展相关问题的深入研究尤为重要。基于电路理论和数值积分原理,针对合理简化后的电路提出了一种基于动态元件离散化的故障电流数值计算方法。首先,通过电路理论分析故障放电机理,对电网进行合理的简化;然后,通过后向欧拉法对动态元件进行离散化处理,建立简化后电路的离散模型;最后,基于简化后电路的离散模型建立改进节点方程,求得各时刻的故障电流值。通过将所提方法的计算结果与在PSCAD/EMTDC中搭建的柔性直流电网仿真模型的结果进行对比分析,验证了所提方法的可行性与通用性。