基于弱可压缩方法的水中声流数值分析

为研究水中的非线性声波传播问题,本文主要开展了基于弱可压缩方程的理论分析及水中声波动及声流模拟工作。首先,采用正压流体的密度与压力函数,结合粘性流体动力学方程,对方程组采用特征变量的波动分析,证明了该方程能描述水中的涡旋和声波传递。然后,本文结合时空守恒元解元算法求解该方程组,并对一些典型气动声学问题进行了计算和对比研究：通过计算分析了二维顶盖半圆腔驱动流,本文计算结果与其他研究结论进行对比分析,证实了该方法可以准确模拟出流体的粘性;通过对平面声波在不同粘性流体中传播进行了数值模拟,声波衰减趋势与理论解完全一致;最后,本文还通过对Rayleigh声流问题进行非定常数值模拟,得到的模拟结果与其他研究学者结论完全一致。这说明,采用二维弱可压缩方程组及时空守恒元解元方法,能够模拟水中的声波传播及耗散过程。     

基于初始压力间断的冲击波模拟方法

冲击波是引起结构破坏的首要因素,抗爆结构性能分析通常有两种方法,一是把冲击波以压力函数的形式直接加于结构上;二是通过分别建立炸药与空气的网格模拟得到冲击波,然后将冲击波荷载以节点速度或加速度的形式耦合到结构结点上进行计算.前一种方法基于许多假定,不够真实,而爆炸模拟存在的问题较多,主要是炸药网格单元较小,由此考虑炸药单元与空气或结构单元网格的协调性会导致单元较多,计算量较大,峰值压力和作用时间的可调整范围有限.为解决上述问题提出用初始压力间断的原理直接生成冲击波的方法,并对决定冲击波峰值压力和作用时间的因素进行了分析.     

多孔介质孔隙尺度下不可压缩流体流动特性SPH模拟

基于孔隙尺度下的多孔介质模型,使用光滑粒子动力学(SPH)数值模拟方法研究多孔介质中不可压缩流体流动现象,分析流体在孔隙尺度下的流动特性.首先通过模拟经典的Poiseuille流和Couette流验证方法的正确性.随后,针对流体在分别由圆柱和方柱为骨架构造的规则多孔介质内的流动情况,同时使用SPH方法和有限元法(FEM)进行模拟,所得的结果十分吻合.最后利用SPH方法对流体在基于贝雷岩心获得的多孔介质模型和以同尺度的圆柱为基元构造的一种人工随机多孔介质模型内的流动进行模拟,证明平均流体速度和单位质量体积力具有很好的线性关系.     

SPH方法后处理的新途径

重点研究并挖掘对化合物分子三维立体结构可视化显示软件RasMol的功能。在用实心球体粒子集合显示连续体,并提出了一套比较完整的基于Rasmol软件的光滑粒子流体动力学（SmoothedParticleHydrodynamics,简称SPH）后处理的新途径。为sPH和其它无网格方法的后处理问题提供可行的方法和新途径。     

用于模拟弱可压缩Navier-Stokes方程的格子气模型

给出压缩因子的形式．通过引入鬼场，消除压缩因子，得到了平衡态分布函数的一般表达式及可压缩的等温的Navier－Stokes方程，同时给出了自映射条件所要求的平衡态分布函数的范围．     

弱可压缩移动粒子半隐式方法对溃坝溢流问题的数值模拟

无网格粒子法(拉格朗日法)是计算流体力学领域最新一代的计算方法。由于不需要生成网格,无网格粒子法能够模拟任意形状的界面及其变形和碎裂,尤其适用于计算如溃坝溢流及波浪破碎等具有复杂自由表面运动的流体力学问题。其中,光滑粒子流体动力学方法(SPH)和移动粒子半隐式法(MPS)是无网格计算方法的典型代表。采用基于MPS法的修正弱可压缩移动粒子半隐式法(weakly compressible corrected moving particle semi-implicit method)构建了用于求解自由表面流问题的数值计算模型,并用于模拟干床面与湿床面上的溃坝溢流问题。计算结果与实验结果的对比分析表明,采用基于WCCMPS无网格法建立的数值计算模型对干、湿床面上溃坝溢流的溢流形态、水柱高度以及波前速度等均给出了良好的计算结果。     

多重赋值语句的翻译方法

本文介绍一种翻译多重赋值语句的方法,用所提出的方法对Pascal S编译程序中翻译赋值语句的assignment过程作了修改。给出了新的assignment过程代码。最后列举一两个编译实例,说明修改后的Pascal S编译程序能成功地编译多重赋值语句。     

确定厚壁圆筒初始屈服压力的一种实验方法

基于厚壁圆筒的弹性失效准则所确定的圆筒的初始屈服压力与材料的屈服极限的关系，设计了一种测定厚壁圆筒初始屈服压力的实验方法，并测得了一组实验数据。通过对该实验数据的分析得到了圆筒的初始屈服压力。该实验值与理论值误差较小，表明了该实验方法具有较好的可靠性。     

高速碰撞问题的SPH方法模拟

给出了全应力张量空间中的光滑粒子流体动力学方法插值公式，利用该方法编制程序对高速碰撞问题进行了数据模拟，给出了正碰撞和斜碰撞2个算例，得到了合理的模拟结果。表明该方法可以运用于冲击力学问题的数值模拟计算。     

非线性赋值集分析的拓扑方法

本书作者给出了许多非线性赋值集分析的拓扑方法领域研究的原创性结果,如经济数学中的不动点定理、次数定理、KKM原理、变分不等式理论、Nash均衡点,最优化中的Pareto最优值及其在最佳逼近理论和偏微分方程边值问题中的应用,都是赋值集分析的拓扑方法领域的创新性成果。另外,     

一类多元多项式的有效赋值方法

一类多元多项式的有效赋值方法张晓鹏（西北大学数学系，西安７１００６９；作者，男，３２岁，讲师）矩形域和立方体上定义的多元多项式是逼近论和计算机辅助设计等方面应用很广泛的函数形式［１．２］，因而人们很关注它的有效赋值问题．本文对矩形域上的Ｂｅｒｎｓｔｅ...     

低雷诺数流动问题的SPH数值模拟及与FPM方法的比较

采用SPH方法对低雷诺数流动问题进行了数值模拟,讨论了初始光滑长度以及核函数影响域大小对结果的影响。典型Poiseuille流和Couette流的模拟结果表明,SPH方法能够很好地模拟低雷诺数流动。并比较了SPH方法和FPM方法的精度和计算效率。     

基于Duffing振子的弱信号初始相位同步算法

针对传统混沌相位同步方法复杂度高,预警准确率低的问题,提出了一种基于耦合振子的混沌同步算法。首先给出了混沌相位同步模型及双耦合振子检测方程;其次,对耦合混沌模型的非线性动力学特性进行了深入分析,并对系统模型的抗噪性和稳定性进行了分析;最后,采用基于电力载波信号的仿真实验证明了算法的有效性。     

排污权初始权分配的两种方法及其效率比较

在分析环境容量的基础上,提出了排污权初始分配的两种方式,建立了两种分配方式的数学模型,分析了两种分配方式下企业购置排污权的行为和排污治污行为并讨论了这两种分配方式效率之间的关系,证明了在不考虑市场交易的交易成本的条件下,排污权初始分配的市场分配方式可以实现计划分配方式的整体最大效益.     

修正变光滑长度SPH方法及其应用

提出了修正的变光滑长度SPH方法及其算法实现.与传统变光滑长度SPH法相比,新方程组基于对称形式核函数近似,并从根本上考虑了变光滑长度影响.方程组中,密度演化方程与变光滑长度方程隐式关联;而SPH动量方程和能量方程在Springel提出的全守恒SPH方程组基础上,改进其分散核近似形式为对称核近似形式得到.由于新方程组密度演化方程和变光滑长度方程隐式耦合,采用迭代求解密度演化方程和变光滑长度方程,显式求解SPH动量方程和能量方程的办法,迭代过程增加的计算量相对很少.给出了2个一维激波管算例和二维 Sedov算例来验证方法的有效性.结果表明,新方法很好地解决了传统SPH法中变光滑长度影响,特别是在模拟二维 Sedov问题时新方法能得到比Springel方法更准确的压强峰值位置,中心压强也更接近理论值.特别适合模拟爆炸与冲击、大变形大扭曲等密度和光滑长度变化剧烈的问题.     

基于SPH方法的碎片穿甲数值模拟

基于显示中心差分格式应用光滑粒子流体力学方法(SPH)和有限元(FEM)两种方法研究在三种不同接触方式下,钢制方块碎片对质地均匀的铝板侵彻过程中的力学过程。计算结果获得金属碎片在侵彻过程中粒子云的扩散分布变化情况。对在不同角度下侵彻过程中的破坏规律作出了一些总结。所得结论为进一步深入研究侵彻及防护问题以及引信设计的参数优化提供一定的依据。     

压力校直下压量计算方法的比较分析

校直下压量是压力校直的主要工艺参数。校直下压量有三种计算方法：基于弹塑性理论基础的解析计算方法,有限元方法和经验公式的计算方法。通过实例计算结果与实验数据比较,证实了计算方法的实用性,对开展校直工艺理论研究和开发设计自动校直机具有一定的指导意义。     

解决SPH方法应力不稳定问题的改进方法

文章研究了光滑粒子流体动力学(smoothed particle hydrodynamics,SPH)方法中的应力不稳定问题,将改进的核函数和改进的动量方程的离散形式相结合,提出一种能同时改善压应力和拉应力不稳定性的方法。针对应力不稳定产生的原因,即应力不稳定性与应力状态和核函数有关,首先提出了一种改进的核函数来解决流体中的压应力不稳定问题,并将其应用到计算机图形学的流体模拟中,对二维和三维溃坝进行仿真实验,并与其他方法进行比较,结果表明,该方法的稳定性和模拟效率都有明显提高。为了同时改善压应力和拉应力的不稳定性,对动量方程的离散形式进行改进,把拉伸状态下的吸引力转化为排斥力,使得排斥力随粒子间距减小而增大,避免了原先吸引力造成的粒子聚集现象,有效改善了拉应力不稳定性。二维液滴在受到表面张力作用下形变的仿真实验效果验证了文中方法能使液滴粒子分布更均匀,稳定性大大提高。