基于SPH方法的表面张力修正算法及其应用

针对传统SPEI方法中基于CSF模型的表面张力算法,在计算边界、尖角等粒子缺失部位的曲率时存在偏差较大,且粒子秩序较差,对大变形问题表面张力计算精度较低的问题,在Morris提出的表面张力SPH方法基础上,通过引入CSPM方法对边界法向的计算和曲率的计算进行修正,得到了表面张力修正方程组.应用本文方法模拟了水溶液中初始...     

计算地球物理流体力学中的非线性计算不稳定问题

一、序论§1 引言计算地球流体力学中的非定常问题(包括数值天气预报、气候数值模拟、海流数值模拟和风暴潮数值预报等)大多是非线性时变偏微分方程的求解问题,无论是用有限差分法、有限元法还是用谱展开法都可能出现非线性计算不稳定。就以数值天气预报和大气环流数值模拟问题来说,一般是对一组复杂的非线性偏微分方程的初、边值问题数值求解,譬     

钢压杆稳定设计的直接计算法

本文根据钢结构设计理论及规范,提出了钢压杆整体稳定设计的直接计算法,既不需试算或迭代计算,又不依赖图表,便于机算设计     

改进的物理融合神经网络在瑞利?泰勒不稳定性问题中的应用

基于相场法的物理融合神经网络PF-PINNs被成功用于两相流动的建模, 为两相流动的高精度直接数值模拟提供了全新的技术手段. 相场法作为一种新兴的界面捕捉方法, 其引入确保了界面的质量守恒, 显著提高了相界面的捕捉精度; 但是相场法中高阶导数的存在也降低了神经网络的训练速度. 为了提升计算训练过程的效率, 本文在PF-PINNS框架下, 参考深度混合残差方法MIM, 将化学能作为辅助变量以及神经网络的输出之一, 并修改了物理约束项的形式, 使辅助变量与相分数的关系式由硬约束转为了软约束. 上述两点改进显著降低了自动微分过程中计算图的规模, 节约了求导过程中的计算开销. 同时, 为了评估建立的PF-PINNS在雷诺数较高、计算量较大的场景中的建模能力, 本文将瑞利?泰勒RT不稳定性问题作为验证算例. 与高精度谱元法的定性与定量对比结果表明, 改进PF-PINNs有能力捕捉到两相界面的强非线性演化过程, 且计算精度接近传统算法, 计算结果符合物理规律. 改进前后的对比结果表明, 深度混合残差方法能够显著降低PF-PINNS的训练用时. 本文所述方法是进一步提升神经网络训练速度的重要参考资料, 并为探索高精度智能建模方法提供了全新的见解.      

弹塑性数值流形方法在边坡稳定分析中的应用

针对现有数值流形方法只能进行弹性计算的不足,建立了一个能够反映完整岩块弹塑性变形特征的本构模型,并借助VC++开发了内置该本构模型的弹塑性数值流形程序。利用该程序模拟了含单节理岩样的室内压缩试验,分析得到了其强度和变形特性,计算结果符合实际的物理现象,表明程序是正确有效的。考虑到数值流形方法本身能够有效模拟材料的不连续变形,新增的弹塑性分析功能又可以反映岩石的强度特性,将弹塑性数值流形程序应用于某含有不连续面的岩石边坡的稳定性分析。并结合锚杆单元的使用,对比分析了不同锚固方案的加固效果。程序提供的变形、应力等计算结果表明：预应力锚杆不仅可以防止不连续面发生剪切破坏,增强坡体的稳定性,限制塑性变形的发展;而且可以使不连续软弱层面对岩体变形的消极影响得以减弱,起到提高岩体整体性的作用。     

混合层中的流向涡与Rayleigh离心不稳定

用有限差分方法求解三维Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程和连续性方程，对时间发展的平面混合层中流向涡的产生原因进行了分析．将Ｒａｙｌｅｉｇｈ的轴对称无粘离心不稳定理论推广应用于分析混合层的二维基本流，并导出一无量纲量Ｒａｙ＝－（ｒ／νθ）νθ／ｒ，其中νθ为一流体质点相对于平均速度的速度，ｒ为该质点迹线的曲率半径．当Ｒａｙ＞１．０时就会发生离心不稳定．采用这一判别式后发现混合层中展向涡的周围，尤其是在辫带区，的确存在离心不稳定区域，而过去的实验结果也表明三维不稳定产生于展向涡之间的辫带区．因此有理由认为，除非雷诺数特别小，离心不稳定是流向涡产生和发展的主要原因     

平板湍流边界层底层的不稳定波

湍流边界层底层的相干结构及猝发现象对决定边界层的特性有很重要的意义,但直至现在对相干结构的性质及起因仍不清楚。在本文中,对湍流边界层的底层进行了非线性稳定性分析,发现该处确可产生不稳定波,可能用以说明相干结构的起因。     

无网格局部Petrov-Galerkin方法在弹塑性断裂力学问题中的应用

采用无网格局部Petroy-Galerkin方法来分析弹塑性断裂力学问题.这种无网格方法采用移动最小二乘法(MLS)来构造近似试函数和采用Heaviside函数作为加权残值法中的权函数,由于近似函数不满足KroneckerDelta条件,因此采用直接插值法来施加本质边界条件.如果不考虑体力,所形成的整体刚度矩阵只包含局部边界积分,而不包含局部域积分和奇异积分.采用增量Newton-Raphson迭代法来求解弹塑性增量形式的局部Petrov-Galerkin方程.数值算例结果表明,该文方法对于弹塑性断裂力学问题的求解是可行的和有效的,并且所得到的结果具有较好的精度.     

液体在不均质双层油层中不稳定一维渗流问题

液体在不均质双层油层中的不稳定渗流问题可以化为解抛物线型微分方程组。本文用Laplace变换解所提问题,得出了一系列解答。这些解答所阐明的不均质双层油层中压力分布的一些性质,对于多油层的弹性驱动油田开发以及油田试井都有一定的意义。同时,这类问题的已有结果是本文解答的特殊情况。     

物理参数直接识别法对于分布系统的应用

本文将物理参数直接识别法用于均匀分布系统等效离散模型的物理参数识别，数值模拟结果表明，识别得到的离散模型的动力特性较好地接近分布系统的动力特性，离散系统物理参数识别的有关结论对分布系统仍然成立，数值模拟与实验结果还表明，均匀分布系统时域响应测量点应等距离布置。     

离散时间系统Ляпунов函数的构造与稳定性判定

本文旨在以矩阵的相似变换理论为基础,以电子计算机为工具,提供一种无需求解离散?矩阵方程,而直接由系统矩阵来构造离散时间系统?函数和判定系统稳定性的方法.     

INSTABILITY OF COUPLED THERMO-SOLUTE CAPILLARY CONVECTION IN LIQUID BRIDGE AND CONTROL BY ROTATING MAGNETIC FIELD

浮区法因具有无坩埚接触污染的生长优点而成为生长高完整性和高均匀性单晶材料的重要技术.但熔体中存在的毛细对流会给浮区法晶体生长带来极大挑战,这是由于对流的不稳定会导致晶体微观瑕疵的产生和宏观条纹等缺陷的形成.为了提高浮区法生长单晶材料的品质,研究浮区法晶体生长中毛细对流特性及如何控制其不稳定性显得尤为重要.本文采用数值模拟的方法对半浮区液桥内Si_xGe_1-x体系中存在的热质毛细对流展开研究并施加旋转磁场对其进行控制.结果表明：纯溶质毛细对流表现为二维轴对称模式,温度场主要由热扩散作用决定,而浓度场则由对流和溶质扩散共同支配;纯热毛细对流呈现三维稳态非轴对称流动,浓度分布与熔体内热毛细对流的流向密切相关,等温线在对流较大的区域发生弯曲;耦合溶质与热毛细对流则为三维周期性旋转振荡流.施加旋转磁场后,熔体周向速度沿径向向外增大,熔体内浓度场和流场均呈现二维轴对称分布.     

表面张力液柱中粘性引起的模分裂

本文证明粘性会使得表面张力所维持的无限长液柱运动模发生分裂(或称分叉),结果是每个轴对称模一次分裂为无穷个。这种现象类似于原生物理学中的 Zeemen 和 Stark效应。无粘条件下的不稳定模分裂后还出现可数无穷个负势能的稳定模,这是关于连续系统模稳定性的 Hare 和 Chandrasekhar 势能判据的一个反例,并解释其物理意义。     

无网格径向点插值法在拓扑优化中的应用

论文基于无网格径向点插值法(RPIM)对连续体结构进行拓扑优化设计.以高斯点的相对密度为设计变量,以结构的柔度最小化为目标函数,利用带惩罚的各向同性固体材料模型(SIMP)和优化准则法,建立了设计变量的迭代格式.利用灵敏度过滤技术有效地消除了点状棋盘格现象.给出了相应的计算流程,并用Fortran程序语言实现其算法.算例证明,应用无网格径向点插值法能够有效地对连续体结构进行拓扑优化设计.     

圆环旋转黏性液体射流空间不稳定性研究

利用线性稳定性理论, 进行了液体黏性对不同旋转强度下圆环旋转液体射流 空间不稳定性影响的研究. 在推导出的三维扰动下具有固体涡核型旋转速度分布的圆环旋转 黏性液体射流色散方程的基础上, 针对中低速射流, 进行了类反对称模式与类对称模式下圆 环旋转黏性液体射流的空间不稳定性分析. 研究结果表明, 对于旋转强度较大的圆环旋转液 体射流, 液体黏性的增加, 不利于射流的破碎; 随着液体黏性的增加, 射流的特征频率和最 不稳定波数减小. 然而, 对于旋转强度较小的圆环旋转液体射流, 液体黏性的增加, 有利于 射流的破碎; 随着液体黏性的增加, 类反对称模式下射流特征频率先减小后增大, 类对称模 式下射流特征频率增大; 随着液体黏性的增加, 类反对称模式下射流最不稳定波数先减小后 增大, 类对称模式下射流最不稳定波数增大.     

幂律流体在单根弯曲毛细管中的毛细上升

基于分形理论和数值模拟的方法, 给出了幂律流体在单根弯曲毛细管的上升高 度和上升累积质量随时间的变化关系曲线. 研究结果表明: 在上升初期阶段, 重力因素可以 忽略; 但随着时间增加, 重力因素的影响越来越大; 幂律流体上升的最大高度只与毛细管直 径$\lambda$、幂律流体密度$\rho $有关, 与毛细管弯曲程度$\tau _0 $和幂指数$n$无关; 幂指数$n$越小, 上升初期上升速度越快; 迂曲度分形维数$D_T $越大, 平衡时吸入的幂律流体质量也越大.     

APPLICATION OF THE PSEUDOSPECTRAL METHOD TO THERMOHYDRODYNAMIC LUBRICATION

The pseudospectral method is used for the first time to solve the thermohydrodynamic lubrication equations for a slider bearing. The orthogonal polynomials used in the series expansions are Lagrangian interpola nts derived from a Legendre basis. Exponential convergence to exact solutions is demonstrated and favourable comparisons with previous work are made.