二维柱坐标系辐射输运方程保球对称的离散纵标格式

 离散纵标格式是计算辐射输运方程的常用格式之一. 但是, 传统的离散纵标格式求解二维柱坐标系辐射输运方程模拟一维球对称问题时, 会出现明显的非对称现象, 球对称性被破坏. 针对该问题, 本文分析了传统离散纵标格式不能够保持球对称性的原因, 提出了空间基于柱坐标系、方向基于球坐标系的辐射输运方程, 并对该方程设计了新的离散纵标格式, 从理论上证明了当空间网格取球对称剖分时该离散格式能够保持一维球对称性的充分必要条件. 通过对真空球区域辐射输运、与物质耦合辐射输运等球对称算例的数值模拟, 验证了该格式的保球对称性, 球对称误差能够达到机器精度. 非对称辐射驱动模型以及非对称网格剖分条件下的数值模拟等算例也取得了较好的结果.     

一类非线性Burgers方程组的基于Hopf-Cole变换的直接间断Galerkin有限元方法

<正>1引言Burgers方程可以作为描述许多物理现象的数学模型,如交通流、激波、扰流问题和连续的随机过程.它还可以用于检验数值方法的效率.由于其具有较广的实用范围,一些学者对其近似解进行了较多的研究.如Adomian分解方法、混合有限差分和边界元方法、样条有限元方法、精确显式有限差分方法、Douglas有限差分格式,直接变分方法和变分迭代方法被用于Burgers方程近似解的研究~([1-13]).Hopf-Cole变换~([14,15])是研究Burgers方程较好的分析工具,利用它可以获得Burgers方程一些精确解.近年来,人们意识到该变换也是一个很好的数值工具并利用其得到了一     

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

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

燃烧系统的离散Boltzmann建模与模拟研究进展

燃烧系统的诸多模拟依托于流体建模, 离散Boltzmann方法(discrete Boltzmann method, DBM) 是近年来发展起来的一种新的流体介观建模方法. 本文简要评述DBM发展的两个方向——Navier-Stokes等偏微分方程的数值逼近解法和复杂系统的微介观动理学建模. 主要介绍在燃烧系统模拟方面DBM已有的工作、新近的思路、与传统流体建模的异同以及近期的研究成果. 本文重点传递的信息为: 作为复杂系统微介观动理学建模出现的DBM在模拟过程中同时给出“流动”及其相伴随的、关系最密切的那部分“热动”非平衡效应; 它为燃烧等复杂系统中各类非平衡行为的描述、非平衡信息的提取、非平衡程度的度量提供了一种简洁、有效的方法; 它所提供的热动非平衡测量量有两类: 一类是直接比较分布函数和平衡态分布函数的动理学矩关系得到的, 一类是来自于Chapman-Enskog多尺度分析给出的热传导和黏性项. 基于第二类DBM, 可以实现(燃烧等)一大类复杂流体系统的多尺度物理建模.     

爆轰波在炸药-金属界面上的折射分析

针对爆轰波在炸药-金属界面上折射时由实验获得的金属折射冲击波压力与经典爆轰波极曲线理论预测的压力存在显著差异这一问题, 本文展开了进一步的理论和数值模拟分析研究. 首先通过分析指出经典爆轰波极曲线理论的缺陷, 并对爆轰波极曲线理论进行了改进, 改进爆轰波极曲线理论给出了炸药爆轰波折射类型以及折射冲击作用点处的压力值. 然后发展了一个基于次特征理论来数值求解爆轰反应流动控制方程的二阶中心型Lagrange方法, 并数值模拟了一个典型的炸药爆轰波折射实验. 改进爆轰波极曲线理论和数值模拟分析结果表明, 爆轰波折射类型有三种:反射冲击波的正规折射、带Mach反射的非正规折射、无反射波的正规折射, 并且金属折射冲击波压力值随入射角增大而单调减小.     

复杂网络视角研究反应系统

 网络科学作为统计物理与复杂性科学的一个重要分支, 在近十几年中取得了长足的发展。巴拉巴希(A-L. Barabási)等人提出的“无标度”模型和瓦茨(D. Watts)等人提出的“小世界”模型用简洁优雅的规则开创性地描述了现实世界的复杂系统。网络的思想源于图论这一研究顶点和边所组成图形的数学分支, 它将系统中的个体看作节点, 个体间的关系看作边, 将节点连边的数量称为度值。“无标度”模型认为, 新加入的节点在连边时服从“偏好依附”规律, 即更加倾向于连接度值大的节点, 最终形成的网络所有节点的度值呈P(k)=k^-a形式的幂律分布(k为度值, a为幂律分布指数, 常见网络的a值在2到3之间), 这种分布下k的取值可以相差悬殊, 无法用一个值标度k的范围, 故模型被称为“无标度”。“小世界”模型解释了著名的“六度分离”现象(图1), 发现现实世界中许多系统介于规则网络与随机网络之间, 具有较短的平均最短路径(两个节点之间的边数)和较高的聚类系数(节点的两个邻居相互也是邻居的概率)。一个远方的朋友, 或是一个出国的亲戚就能大大扩展社交范围;人与人之间往往只需要通过几个人就能联系上, 相对于几十亿的人口, 我们所处的世界真不大。这些研究工作将图论与统计物理、计算机、生物乃至社会学等诸多学科进行了交叉, 将科学的疆域拓展到了新的广度与深度。     

PBX炸药粘弹性统计微裂纹本构模型的改进

基于描述PBX炸药本构关系的粘弹性统计微裂纹模型(Visco-SCRAM),通过修正其微裂纹方向矢量的取值范围、初始缺陷尺寸和体量本构关系,得到模型的改进形式。将改进的模型嵌入有限元程序中,开展PBX炸药JO-9159的平板撞击实验数值模拟研究,讨论了微裂纹数密度参数对计算结果的影响。结合实测曲线,对模型改进前、后的计算结果进行对比和分析。结果表明,改进的模型能够更合理地描述低速冲击下PBX炸药的损伤和破坏。     

透明导电薄膜ITO对浅面浮雕VCSEL的影响

对于顶面出光的浅面浮雕VCSEL结构,有源区的电流密度分布的不均性制约着单模稳定性的提高。为此,提出了一种新型结构:氧化铟锡透明导电薄膜(ITO)浅面浮雕VCSEL。该结构不仅能够增大高阶模式的阈值增益,还能够提高基模的增益,实现基模对高阶模式的稳定抑制。研究了ITO的厚度对阈值增益的影响及ITO对VCSEL有源区电流密度分布的影响。研究结果表明:在ITO的厚度为半波长的整数倍时,基模对高阶模式的限制作用最强;ITO通过改善VCSEL有源区的电流密度分布,达到了增大基模的增益和降低高阶模式增益的目的,同时还降低了串联电阻和外电压。     

核壳纳米结构壳中Yb离子对稀土上转换发光温度特性的影响

采用油酸盐法分别制备出均匀的上转换发光裸核纳米粒子及其包覆具有不同Yb³⁺浓度掺杂的NaYF₄:Yb³⁺,Er³⁺核壳纳米结构的上转换纳米粒子。在不同温度下(90~450 K),研究分析了在壳中掺杂不同浓度Yb³⁺的NaYF₄:Yb³⁺,Er³⁺@NaYF₄:x%Yb³⁺核壳纳米体系的上转换发光特性。结果表明:在NaYF₄:Yb³⁺,Er³⁺上转换体系中,惰性壳中的525 nm(²H_11/2 →⁴I_15/2)发射峰呈现出与活性壳中不一样的趋势。壳层中掺杂的Yb³⁺通过声子对纳米粒子内部发光与表面及外界之间的相互作用起到了重要的"桥"连作用。     

学习于敏院士的科学精神

 2015年1月9日,习近平主席亲自授予于敏院士国家最高科学技术奖,表彰他在国防科研上作出杰出成就,为国家安全和科技进步作出的卓越贡献。