构造可压流体的格子玻尔兹曼模型是计算流体领域的难点,现有的模型,要么不能稳定地模拟高马赫流,要么运算效率偏低。本文提供一个可压流体的高效格子玻尔兹曼模型。本工作借用Kataoka和Tsutahara[Physical Review E 69,035701(R)(2004)...
详细信息
构造可压流体的格子玻尔兹曼模型是计算流体领域的难点,现有的模型,要么不能稳定地模拟高马赫流,要么运算效率偏低。本文提供一个可压流体的高效格子玻尔兹曼模型。本工作借用Kataoka和Tsutahara[Physical Review E 69,035701(R)(2004)]的二维十六速离散速度模型,引入合理的有限差分格式、色散项和耗散项。色散项有效地减小间断处的振荡,提高数值精度。耗散项使新模型更容易满足Von Neumann稳定性条件。新模型成功模拟了几个一维Riemann问题,以及冲击波反射问题。这些数值试验表明新模型在低速和高速流系统中均适用,且可以达到较高的精确性和稳定性。
讨论非结构网格上温度扩散方程的能流计算方法.应用有限点方法(Finite Point Method,简称FPM)导出基于有限点两点公式和三点公式的能流计算公式,该公式适用于任意多边形及非匹配网格等非结构网格;给出网格角点温度新的计算公式.数值试...
详细信息
讨论非结构网格上温度扩散方程的能流计算方法.应用有限点方法(Finite Point Method,简称FPM)导出基于有限点两点公式和三点公式的能流计算公式,该公式适用于任意多边形及非匹配网格等非结构网格;给出网格角点温度新的计算公式.数值试验表明:基于两点公式的离散解和基于三点公式的离散解均具有平方阶的收敛速度;基于三点公式的离散解的精度总优于基于两点公式的离散解.
用并行化的MCNP(Monte Carlo N-particle)程序模拟了美国堪萨斯州立大学1977的光子散射实验基准问题。针对不加顶盖屏蔽层的情况,对问题的几何和材料进行了详细的描述,模拟了地面、水泥柱壳屏蔽体和小车的散射。由于问题并不是完全轴对...
详细信息
用并行化的MCNP(Monte Carlo N-particle)程序模拟了美国堪萨斯州立大学1977的光子散射实验基准问题。针对不加顶盖屏蔽层的情况,对问题的几何和材料进行了详细的描述,模拟了地面、水泥柱壳屏蔽体和小车的散射。由于问题并不是完全轴对称的,因而采用的是点探测器估计,用100个处理器模拟了1亿个样本,给出了计算结果和实验值的对比,模拟结果和试验结果存在一些差别,但趋势相同。同时,用指向概率法的MCCO程序与MCNP进行了比较。模拟中还考虑了加重要性和不加重要性对计算效率的影响。
根据Snyder解析模型建立了一个含有肿瘤的BNCT(硼中子俘获治疗)-4 mm网格模型,应用MCNP(Monte Carlo N-particle)程序进行模拟计算,并对结果进行了物理分析;使用新研制的MCDB(Monte Carlodosimetry in brain)程序对Snyder加肿瘤模型进...
详细信息
根据Snyder解析模型建立了一个含有肿瘤的BNCT(硼中子俘获治疗)-4 mm网格模型,应用MCNP(Monte Carlo N-particle)程序进行模拟计算,并对结果进行了物理分析;使用新研制的MCDB(Monte Carlodosimetry in brain)程序对Snyder加肿瘤模型进行了模拟,得到和MCNP程序几乎完全一致的结果,验证了MCDB程序的正确性,由于MCDB采用了适合均匀网格的快速粒子径迹算法,因而提高了计算速度,较MCNP节省了36%的模拟时间。
暂无评论