首页期刊导航|计算物理
期刊信息/Journal information
计算物理
计算物理

沈隆钧

双月刊

1001-246X

jswl@iapcm.ac.cn

010-62014411-2647、2292、2171

100094

北京市海淀区丰豪东路2号

计算物理/Journal Chinese Journal of Computational PhysicsCSCD北大核心CSTPCDEI
查看更多>>本刊是由中国科学技术协会主管、中国核学会主办、北京应用物理与计算数学研究所承办的计算物理学科综合性学术刊物，其宗旨是充分反映我国计算物理学科前沿的研究成果，促进国内外学术交流，推进计算物理学的发展。计算物理是物理学与科学计算相结合的交叉学科。在物理学相关领域内，凡从事物理建模和计算方法研究，并应用先进计算机获得计算结果的工作，均属投稿范围。本刊特别欢迎在国民经济建设和高科技领域内，运用高性能计算方法解决涉及物理学问题的有创造性的研究成果。
正式出版
收录年代

    基于Ginzburg-Landau方法的快速结构复振幅展开晶体相场模型

    王昆陈军王裴祝文军...
    547-558页
    查看更多>>摘要：利用 Ginzburg-Landau模型思路推广了传统复振幅扩展晶体相场(APFC)模型的发展思路,提出一种描述不同晶体结构的简单有效方法,即快速结构 APFC模型.以方形相和矩形相为例,系统地确定了快速结构 APFC 模型中与结构相关的模型参数,并通过数值算例检验了该方法的有效性.特别是在处理矩形相的过程中,研究发现这种方法不仅可以解决矩形相结构稳定性的问题,还可以描述矩形相和正交层状相之间的结构相变,证明了该模型具有描述多结构相变的能力.最后,通过模拟研究经典的圆形晶粒转动-收缩问题,检验了该模型对物理规律的准确预测能力,并揭示了不同晶体对称性对晶粒转动-收缩规律的影响.

    晶体相场Ginzburg-Landau模型位错晶界晶体各向异性

      原文链接:
    • 万方数据
    • 维普

    非均质炸药的Eshelby夹杂理论相场方法

    徐云龙瑶向美珍陈军...
    559-568页
    查看更多>>摘要：考虑到夹杂-基体界面微结构复杂失效过程对非均质 PBX(高聚物黏结炸药)材料损伤起始具有主导作用,本文提出基于 Eshelby 夹杂理论的相场损伤模型,对非均质 PBX 炸药损伤形核演化进行数值模拟.相场能量由弹性能和夹杂-基体相互作用能构成,结合 Mori-Tanaka方法推导不同级配下 PBX 炸药等效力学模量,相场内变量的变化直接反映了非线性脱黏效应下的损伤分布.采用 Eshelby 夹杂理论相场方法计算高模量比非均质 PBX炸药圆形和多边形颗粒夹杂典型结构,分析加载条件、颗粒形状、体积占比、模型参数对夹杂-基体界面脱黏机制的影响.数值结果表明:夹杂-基体界面微结构演化加速了界面脱黏与宏观损伤的形成,与实验观测是一致的.

    非均质PBX夹杂损伤演化相场模型Eshelby理论

      原文链接:
    • 万方数据
    • 维普

    弹塑性问题壁热误差分析及抑制方法

    李肖沈智军郭虹平方俊...
    569-581页
    查看更多>>摘要：使用单元中心型拉氏 Godunov 方法,研究理想弹塑性流体的数值格式和减少壁热误差的方法.给出Godunov格式的黏性修正方程,描述初值间断情况下黏性激波的形成和传播过程,分析修正方程的黏性行为和壁热误差的关系.在此基础上,提出一个 HLLC 型近似黎曼解法器.在该解法器中,引入一种自适应热通量黏性项用来抑制界面处的内能和密度壁热误差;提出一个额外接触速度用来抑制偏应力的过热现象.

    弹塑性流体黎曼解法器HLLC壁热

      原文链接:
    • 万方数据
    • 维普

    液态重金属反应堆燃料棒及控制棒设计

    徐昌恒潘晖何明涛赵常有...
    582-588页
    查看更多>>摘要：旨在从中子学角度分析液态重金属快堆燃料棒和控制棒的物理特性,以支撑其结构设计及材料选型.使用快中子反应堆分析程序 SARAX对燃料棒的液铀比选择、径向功率分布、包壳材料选型、反射层材料选型,控制棒径向结构设计、吸收体材料选型、慢化结构设计进行计算分析.计算结果表明:燃料棒栅元 Kinf 与液铀比呈线性关系,液铀比的选择更多需要考虑载热因素.燃料棒包壳优先选择铁素体/马氏体钢材料.活性段轴向两端反射层材料可选择钢类材料.基于 Pu迁移现象分析结果,当前物理计算程序中均未耦合 Pu 迁移现象,计算得到的功率分布相对均匀,导致燃料中心温度计算结果偏不保守.控制棒采用细棒设计能有效降低自屏效应,但会降低空间利用率.快堆常用 B4C作为吸收体材料.控制棒吸收体外包裹慢化体能使中子能谱软化,B-10 中子吸收截面大幅上升,提高控制棒价值.同一外径下,纯吸收体芯块的棒价值小于使用慢化体(ZrH2)包裹的芯块.

    液态重金属反应堆燃料棒设计控制棒设计

      原文链接:
    • 万方数据
    • 维普

    基于G-TF/SF技术的地面HEMP环境研究

    杨超谢海燕翟鑫洋乔海亮...
    589-595页
    查看更多>>摘要：对于半空间问题的时域有限差分方法(FDTD)解决方案,基于传统总场散射场(TF/SF)技术的平面波加载方法无法直接用于地面电磁散射问题中,本文将广义总场散射场(G-TF/SF)方法推广至三维地面的高空电磁脉冲(HEMP)环境计算,四周采用卷积完全匹配层(CPML)吸收边界条件.G-TF/SF 方法可有效消除数值计算中截取无限大地面产生的边缘效应,提高地面电磁环境的计算精度,且 G-TF/SF边界只需要加载入射场,避免了反射场与透射场的计算.

    时域有限差分广义总场散射场高空电磁脉冲环境地面

      原文链接:
    • 万方数据
    • 维普

    燃料棒接触问题的新型自适应预处理方法

    刘威震刘振海辛勇叶舒愉...
    596-606页
    查看更多>>摘要：通过分析接触前后雅可比矩阵的稀疏结构,构造基于接触压力判断准则的新型自适应预处理方法.在接触压力为零时使用完全的雅可比矩阵作为 Preconditioned Jacobian-free Newton-Krylov(PJFNK)方法的预处理矩阵,在接触压力非零时使用雅可比矩阵的块对角作为预处理矩阵.数值试验表明:本文的自适应预处理方法能够有效提升模拟效率,求解总时间与 Multiphysics Object-Oriented Simulation Environment(MOOSE)中现有方法相比减少了15.2%～33.1%.

    MOOSEPJFNK燃料棒力学接触自适应预处理

      原文链接:
    • 万方数据
    • 维普

    基于界面重构的三维弹塑性MMALE计算方法

    郭少冬周海兵熊俊
    607-618页
    查看更多>>摘要：针对弹塑性流体特点,通过弹塑性混合网格封闭模型、三维大变形网格上界面重构以及保物理特性弹塑性量重映方法等技术,给出一种适应弹塑性流体的三维 MMALE(multi-material arbitrary Lagragian-Eulerian)算法.首先,在混合网格封闭模型方面,针对传统体模量加权模型应用于弹塑性流体导致的压力不平衡问题,通过引入松弛机制和自适应松弛系数,实现压力平衡,同时该封闭模型无须迭代,避免了在迭代过程可能遇到的不收敛问题,提升了健壮性;在界面重构方面,针对三维多面体网格由于变形及三维特殊性导致的守恒性问题,通过引入逐级表面三角化算法,实现了适应三维任意多面体网格上的严格守恒的界面重构.在弹塑性量重映方面,针对偏应力按分量重映后导致的张量性质破坏,应变能不守恒的问题,通过引入对应力张量不变量的重映,实现重映后应力不破坏屈服准则,并保持应变能的守恒性.最后,采用 MMALE算法对典型算例进行数值模拟,验证了其正确性和处理大变形问题的健壮性.

    弹塑性多介质界面重构MMALE

      原文链接:
    • 万方数据
    • 维普

    含裂缝线弹性问题的压缩预条件共轭梯度算法

    刘杏康陈星玎余云龙
    619-629页
    查看更多>>摘要：基于几何型扩展有限元离散,研究含静态裂缝线弹性问题的高效压缩预条件算法.不仅构造适用于含裂缝线弹性问题的压缩子空间矩阵,而且给出压缩点的选点原则.为进一步提高计算效率,将该压缩技巧以乘性的方式与"裂尖型"区域分解预条件子相结合,提出一类高效的自适应压缩预条件共轭梯度算法,该算法能同时消去迭代求解中的高频误差和低频误差,数值实验验证了算法的有效性.

    扩展有限元方法压缩技巧区域分解预条件子静态裂缝问题

      原文链接:
    • 万方数据
    • 维普

    一种稳健的水力压裂裂缝平面识别算法

    林子愉刘月田裴雪皓樊平天...
    630-642页
    查看更多>>摘要：提出一种稳健的水力压裂裂缝平面识别算法——抽样投影算法(RANSAC-MP),通过随机抽样方法弱化无关破裂事件导致的离群噪声,并提出最大投影平面拟合算法以减小环境噪声的影响,同时结合了 RANSAC 算法的抗噪性和投影方法的降维效果.以实际直井微地震数据为例进行裂缝平面拟合,结果表明:RANSAC-MP 算法在多重噪声影响下表现出更强的稳健性和更高的计算精度,当压裂仅形成单一裂缝时,该算法可以直接处理原始数据.

    多重噪声影响抽样投影(RANSAC-MP)裂缝平面识别水力压裂微地震数据点稳健性

      原文链接:
    • 万方数据
    • 维普

    反演分布函数的多块网格LBM局部加密方法

    魏如普丁鹏
    643-650页
    查看更多>>摘要：针对格子玻尔兹曼方法(LBM)局部加密复杂的耦合计算,提出从宏观物理量反演得到粒子分布函数以避免进行分布函数转换的思路.通过重写 Boltzmann 方程推导得到分布函数内部关于 Knudsen 数的关系式,给出分布函数的反演计算方法,以此为基础设计新的多块网格加密算法,并将该算法应用于算例中验证其性能.结果表明:该算法能够较好地展现不同 Reynolds数下流场的数值特征,流线图连续完整,得到的流场中涡旋的中心坐标与文献数据相符,可以用于流场模拟计算,为局部加密方法的发展提供指导.

    格子玻尔兹曼方法多块网格分布函数反演局部加密

      原文链接:
    • 万方数据
    • 维普