模拟颗粒凝并过程的快速Monte Carlo方法

对常规异权值Monte Carlo(MC)方法进行改进,基于强核函数思想,通过颗粒群单重遍历即可求得强核函数最大值,采用接受-拒绝法随机搜寻凝并对,并利用搜寻过程中拒绝和接受的所有凝并对的信息来估计凝并事件的等待时间(时间步长),从而避免颗粒群的双重遍历,以提高MC的效率.对典型工况的模拟结果显示该快速方法计算代价仅为O(N_s),能够显著提高计算效率,同时保持足够的计算精度,较好地协调计算代价与计算精度之间的矛盾.     

利用液体夹心法测量LiF单晶高压热导率

提出了液体夹心法热导率测量技术, 采用CHBr3液体作为夹心材料实现了样品/窗口界面的理想接触, 并将动载荷作用下的夹心法高压热导率实验测量压力下限拓展至40 GPa, 为绝缘介电晶体高温高压热导率测量提供了技术支持. 实验利用平面碰撞和DPS测试技术, 结合液体夹心法实测了LiF单晶高压热导率数据, 对现有热导率理论模型进行了研究和探讨, 结果显示, 在γ/γ0=(ρ0/ρ)2时, 修正后的Roufosse理论公式与实验数据符合较好, 这一研究结果为非透明材料冲击波温度测量中的热传导修正提供了实验数据和理论模型。     

基于表面阻抗边界条件的等离子体薄涂层电磁散射的时域有限差分分析

基于表面阻抗边界条件时域有限差分(FDTD)方法研究了一维斜入射情况下非磁化等离子体薄涂层涂敷金属材料的电磁散射特性, 该方法忽略对薄层背景材料进行网格剖分, 大大减少了计算量. 首先推导了理想导体涂敷等离子体薄涂层的表面阻抗频域表达式, 然后代入边界条件并变换到时域, 再用分段线性递推卷积方法将时域表达式离散得到FDTD迭代式. 编程计算了垂直及斜入射情形下的平行极化和垂直极化反射系数, 通过验证算例与解析解对比, 结果表明该方法的准确性和有效性. 最后利用该方法分析了不同入射角对反射系数的影响.     

修正的BBM方程新的精确解

应用进一步修正的简单方程法对修正的 Benjamin -Bona -Mahoney (mBBM )方程进行求解,给出了mBBM方程新的精确类孤波解,取定某些参数值,便可得到精确孤波解.这种方法也可用于寻找其它常系数以及变系数非线性发展方程(组)的精确解,具有一定的普适性.     

Novel loop-like solitons for generalized Vakhnenko equation

A nontraveling wave solution of a generalized Vakhnenko equation arising from the high-frequent wave motion in a relaxing medium is derived via the extended Riccati mapping method. The solution includes an arbitrary function of an independent variable. Based on the solution, two hyperbolic functions are chosen to construct new solitons. Novel single-loop-like and double-loop-like solitons are found for the equation.     

光滑粒子流体动力学方法固壁处理的一种新型排斥力模型

与传统网格法相比, 光滑粒子流体动力学方法不能直接施加壁面边界条件, 这就限制了该方法在工程中的应用.为此, 本文基于Galerkin加权余量法并结合传统排斥力方法, 推导出一种新的排斥力公式来施加壁面边界条件.该方法不含未知参数, 能在不减小边界粒子尺寸的情形下有效地防止流体粒子穿透壁面, 同时可避免邻近边界的流体粒子的速度及压力振荡. 分别通过静止液柱算例、液柱坍塌算例、容器中液体静止算例及溃坝算 例来验证本文方法的有效性, 并与传统边界处理方法进行对比, 结果表明: 本文方法克服了传统方法存在的缺陷, 是一种有效的固壁边界处理方法. 关键词： 光滑粒子流体动力学法 固壁边界 排斥力 加权余量法     

梯度计算的集合变分方案及其在大气Ekman层湍流系数反演中的应用

大气Ekman层湍流系数的准确计算对数值天气预报和污染物扩散计算有着重要的意义. 将集合计算和变分法结合起来, 提出了目标泛函梯度计算的集合变分方案, 并根据正演模式的线性化情况提出了两种计算流程. 利用这种集合变分梯度算法及两种流程对Ekman层湍流系数进行了反演试验, 结果表明这种算法实施简单、方便, 根据观测资料能够比较准确地反演湍流系数值. 关键词： 梯度计算 集合变分 湍流系数 反演     

改进的节块展开方法求解圆柱几何对流扩散方程

为高效求解球床高温气冷堆物理-热工耦合问题,发展改进节块展开法求解圆柱几何下的对流扩散方程.针对圆柱几何和对流扩散方程的特殊性,采用三阶多项式和指数函数作为r向横向积分方程的展开函数,在节块展开法的框架下高效求解对流扩散方程.数值验证表明,改进的节块展开方法具有固有的迎风特性,在使用粗网节块时依然能保持稳定性和较高的计算精度.     

CTRW模型等待时间分布的随机抽样方法

研究等待时间分布为截断幂律分布的连续时间随机行走(CTRW)模型,分析比较截断幂律分布的可能随机抽样方法,并提出一套可以对其进行精确抽样的乘分布舍选算法.数值结果表明:当幂律指数α > 0.2时,该算法的抽样效率在80%以上.利用该算法可以成功得出CTRW模型描述的三种反常输运过程的标度关系.     

A Pseudo-Stokes Mesh Motion Algorithm

This work presents a moving mesh methodology based on the solution of a pseudo flow problem. The mesh motion is modeled as a pseudo Stokes problem solved by an explicit finite element projection method. The mesh quality requirements are satisfied by employing a null divergent velocity condition. This methodology is applied to triangular unstructured meshes and compared to well known approaches such as the ones based on diffusion and pseudo structural problems. One of the test cases is an airfoil with a fully meshed domain. A specific rotation velocity is imposed as the airfoil boundary condition. The other test is a set of two cylinders that move toward each other. A mesh quality criterion is employed to identify critically distorted elements and to evaluate the performance of each mesh motion approach. The results obtained for each test case show that the pseudo-flow methodology produces satisfactory meshes during the moving process.