基于移动粒子半隐式法的表面张力模拟

采用移动粒子半隐式法(MPS)模拟了受表面张力影响的自由面流动。表面张力的计算采取了一种较适合于MPS方法的表面自由能模型。方形液滴振荡和射流断裂的模拟结果分别与理论分析和试验结果一致,同时进行了三维射流注水模拟,从而验证了MPS方法结合该表面张力模型可以有效、方便地进行自由面流动中表面张力现象的模拟。     

平面叶栅跨音流场气动特性的数值研究

通过对模型方程的分析，给出了一种新的隐格式构造思想。将它运用到关通量分裂格式中，可得到无近似因子分解、无矩阵运算的高效二阶精度隐式矢通量分裂差分格式，并用来直接求解时间平均Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程组。数值计算标明：该方法具有精度高、稳定性好、计算量少、收敛快等优点，在平面叶栅跨音流场的计算中，较好地捕获了激波，与实验比较，结果令人满意。     

推进－迭代算法计算效率的研究

本文对求解三维定常超音速动性流场的一次空间推进，在每一个推进站沿伪时间层局部迭代的推进－迭代算法作了进一步的研究．在每一推进站（侧向平面）沿伪时间层局部迭代时，给出了四种不同的隐式迭代方法，即沿侧面两个方向（法向和周向）全用隐式；法向隐式而周向采用Ｇａｕｓｓ－Ｓｉｌｄｌｅ来回扫描迭代；法向隐式而周向显式及以系数矩阵谱半径代替系数矩阵的简化标量隐式算法．用这四种算法模拟了三维球锥黏性绕流，给出了四种不同算法的计算效率和收敛特性比较．     

不可压缩流动的移动粒子半隐式法数值模拟

介绍了一种新的移动粒子半隐式方法。该方法采用Lagrange描述,用一系列离散粒子代替流体,通过半隐式求解和时间推进法来进行计算,弥补了传统网格法对复杂形状、大变形、高速撞击等情况下网格划分和重构过于繁琐的缺点,在模拟不可压缩流动问题中得到了较好的效果。     

岩土弹塑性有限元分析的隐式积分弹性刚度法

本文提出了用于岩土弹塑性有限元分析的隐式积分弹性刚度算法。该算法既具有隐式积分法精度好，效率高，无条件稳定等优点，也具有弹性刚度矩阵正定、对称的特点，更重要的是它避免了传统切线刚度法在处理岩土非相关联塑性流动和屈服面“角角”所遇到的非对称性和奇异性计算问题。通过算例分析了该算法的精度、效率。     

显式模拟类橡胶材料Mullins效应滞回圈

通过显式、直接的方法提出一个多轴可压缩应变能函数,用来模拟类橡胶材料在加载——卸载作用下,由于Mullins效应而产生的应力——应变滞回圈. 本文的创新点在于将表征能量耗散的变量引入到应变能函数.新的弹性势具有以下两个特点:第一,在加载情况下,新引入的变量不会对弹性势产生任何影响,因此,只要给出合适的形函数显式表达,3个基准实验,包括单轴拉伸和压缩,等双轴拉伸和压缩,以及平面应变,都可精确模拟;第二,新引入的变量在卸载情况下将被激活.在不同的卸载应力下,变量将发生改变,从而影响弹性势,使其最终产生不同的应力——应变关系卸载曲线,与对应的加载曲线共同构成应力——应变滞回圈.通过对Mullins效应实验数据进行分析和研究,得出了卸载形函数在不同卸载应力下变化的规律,并预测不同卸载应力下的应力——应变关系.最后,我们将得到精确匹配实验数据的数值模拟结果,从而证明本文方法不仅可以精确匹配至少3个基准实验,还可以模拟和预测类橡胶材料在加载——卸载作用下由于Mullins效应而产生的滞回圈.      

无网格粒子法壁面边界特性研究

针对无网格粒子法下的离散壁面边界,通过基于粒子数密度的机理分析,证明了离散数值边界与真实边界之间存在差异。分别在二维和三维尺度下,通过计算并分析典型情况下近壁面零散流体粒子的运动轨迹,验证了该差异的存在并展示了壁面对流动产生影响的方式和相关现象。结果表明,数值壁面边界未能完全准确描述真实物理边界,壁面对流体的影响区域并非完全非光滑平顺,这对近壁面流动尤其是低速零散流体粒子会产生一定影响。     

ON THE EXPECTED DISCOUNTED PENALTY FUNCTION IN A MARKOV-DEPENDENT RISK MODEL WITH CONSTANT DIVIDEND BARRIER

This article considers a Markov-dependent risk model with a constant dividend barrier. A system of integro-differential equations with boundary conditions satisfied by the expected discounted penalty function, with given initial environment state, is derived and solved. Explicit formulas for the discounted penalty function are obtained when the initial surplus is zero or when all the claim amount distributions are from rational family. In two state model, numerical illustrations with exponential claim amounts are given.     

基于碎点法的动态断裂分析

工程中的冲击防护结构在撞击、爆炸等冲击载荷下可能发生动态断裂并最终破坏, 抑制结构的动态断裂是提升结构防护能力的重要手段, 为此需要准确预测结构在动态载荷下的断裂行为. 数值仿真是预测动态断裂的重要手段, 然而当前工程中常用的有限元法在模拟断裂方面仍存在网格畸变和难以显式引入裂纹等问题. 碎点法是近年来提出的一种不连续型伽辽金弱形式无网格方法, 适合模拟断裂问题, 本文提出一种显式动力学格式的碎点法并将该方法应用于动态断裂分析. 一方面, 碎点法参考弱形式无网格类方法, 将求解域离散为空间中的节点和子域, 并基于支持域内的节点群构造子域的位移试函数, 因此该方法的子域具有抵抗畸变的能力. 另一方面, 碎点法参考间断伽辽金有限元法, 使用分片连续的位移试函数, 并引入内部界面数值通量修正保证方法的一致性和稳定性, 因此该方法易于在结构中显式引入裂纹. 本文首先介绍碎点法的核心思想和离散形式, 接着推导了动力学碎点法弱形式动量方程, 然后建立了碎点法的显式动力学求解格式, 最后通过算例验证动力学碎点法预测应力波传播和动态断裂行为的能力.      

Three-stage Stiffly Accurate Runge-Kutta Methods for Stiff Stochastic Differential Equations

In this paper we discuss diagonally implicit and semi-implicit methods based on the three-stage stiffly accurate Runge-Kutta methods for solving Stratonovich stochastic differential equations（SDEs）.Two methods,a three-stage stiffly accurate semi-implicit（SASI3） method and a three-stage stiffly accurate diagonally implicit （SADI3） method,are constructed in this paper.In particular,the truncated random variable is used in the implicit method.The stability properties and numerical results show the effectiveness of these methods in the pathwise approximation of stiff SDEs.