一种基于湍动能方程的转捩判定方法

转捩现象是阻碍阻力高精度求解的主要问题之一. Menter 和Langtry 所提出的γ-θ转捩模型通过引入涡量雷诺数和间歇因子输运方程来驱动转捩,但是其中很多经验公式的理论立足点有待商榷. 驱使层流转变到湍流依赖的仍然是平均速度的一阶和二阶相关量,它们组合构成了湍动能方程的耗散尺度. 在湍动能方程中做合适的耗散平衡后,仅仅依靠湍动能方程可以有效地捕捉转捩现象. 采用自然转捩和旁路转捩测试算例进行了验证,结果证明该方法与试验值匹配较好,具有一定的工程实用价值.      

结构可靠性仿真方法研究简

可靠性仿真是对复杂结构系统进行可靠性分析计算最为有效的方法,引起了研究者越来越广泛的关注. 综合目前结构可靠性仿真方法研究现状,介绍了蒙特卡罗（Monte-Carlo,MC）法、极限状态方程重构法、随机有限元法（stochastic finite element method,SFEM）进行结构可靠性仿真计算的特点及发展现状,归纳了固体火箭发动机（solid rocket motor,SRM）结构可靠性的特点及当前进行仿真的方法;并在此基础上分析了目前结构可靠性仿真方法及对SRM 结构可靠性进行仿真分析计算时存在的不足及需要进一步研究的内容.     

新型涡轮流量计的理论模型简

基于动量力矩定理和流体边界层理论,建立了一新型涡轮流量计的理论模型. 利用该理论模型可以分析叶轮几何参数对流量计计量性能的影响. 以一DN50 涡轮流量计为例,利用本文所建立的理论模型对该流量计的仪表系数进行了计算,并在一体积管标准装置对该流量计的仪表系数进行标定. 结果表明,计算值和实验值较为吻合,误差在±;3.5% 以内,验证了模型的有效性,从而为仪表的结构优化设计提供了理论依据.     

基于有限变形理论的数值流形方法研究简

原有数值流形方法通过积累每一时步的小变形而得到结构最终的大变形,然而,当结构发生大变形、大转动时往往产生较大计算误差. 针对该问题,从动量守恒方程以及应力边界条件的积分弱形式出发,引入流形方法的插值函数,建立了基于有限变形理论的数值流形方法. 通过对比改进前后流形方法的计算迭代格式,指出了原有流形方法计算大变形问题时的误差来源. 最后,通过大变形悬臂梁和旋转块体算例对有限变形流形方法进行了验证. 数值结果表明,改进后的流形方法能够很好地处理大变形大转动问题,消除了转动所带来的计算误差,其计算结果与解析解及ABAQUS 软件求得的数值解相吻合.     

绕振荡水翼流动及其转捩特性的数值计算研究简

通过对比标准k-ω SST 湍流模型和基于标准k-ω SST 湍流模型修正的γ-Re_θ 转捩湍流模型对绕振荡NACA66 水翼流动的数值计算结果与实验结果,对水翼振荡过程的水动力特性和流场结构变化进行了分析研究. 结果表明：与标准k-ω SST 湍流模型的数值计算结果相比,基于标准k-ω SST 湍流模型修正的γ-Re_θ 转捩湍流模型能有效预测绕振荡翼型流场结构和水动力特性,捕捉流场边界层发生的流动分离和转捩现象;绕振荡水翼的流动过程可分为5 个特征阶段,当来流攻角较小时,在水翼前缘发生层流向湍流的转捩现象,水翼动力特征曲线出现变化拐点;随着来流攻角的增大,顺时针尾缘涡逐渐形成并向水翼前缘发展;当攻角较大时,前缘涡分离导致动力失速,水翼的动力特征曲线出现大幅波动;水翼处于顺时针向下旋转阶段,绕水翼的流动状态逐渐由湍流过渡为层流.     

基于能量耗散的岩石脆性评价方法

脆性评价是页岩气钻井以及压裂工艺设计的重要环节.然而,现今对于脆性的定义十分模糊,且脆性评价的方法没有成熟的理论作为支撑.本文以岩石应力应变关系为基础,采用能量法来进行储层岩石脆性评价.根据断裂力学裂纹尖端的非线性场理论,以前端区作为尖端特殊的屏蔽构型,来描述储层岩石在断裂过程中能量的耗散. 通过计算本征内聚力与前端区稳态阻力比值来表征储层岩石的脆性指数. 研究表明,能量耗散相对越大,岩石的脆性指数越低,其塑性特征表现得越明显. 基于能量耗散的储层岩石脆性评价是分析岩石力学特性的重要内容,也为压裂选层提供了重要依据.      

一种适用于超音速边界层的湍流转捩模式

建立一种合理反映扰动模态和可压缩性影响的新型k-ω-γ转捩模式.其主要特点为:(1)假设脉动动能k 由湍流脉动部分和非湍流脉动部分组成,且在模化后者时采用了稳定性分析的结果; (2)在间歇因子γ方程的源项中,构造了具有``自动判断转捩起始位置'功能的函数; (3)通过构造新型的物面法向长度尺度,保证了模式中所有的表达式均由当地变量构成,可以方便地应用于现代CFD程序之中. 该模式在亚音速、超音速和高超音速条件下的边界层流动中进行了验证. 计算结果表明,该模式可应用于较宽马赫数范围内的自然转捩以及旁路转捩过程,所具有的捕捉流动转捩的性能优于国际上的现有模式.      

稀薄流到连续流的气体运动论模型方程算法研究

通过引入碰撞松弛参数和当地平衡态分布函数对BGK模型方程进行修正，确定含流态控制参数可描述不同流域气体流动特性的气体分子速度分布函数的简化控制方程。发展和应用离散速度坐标法于气体分子速度空间，利用一套在物理空间和时间上连续而速度空间离散的分布函数来代替原分布函数对速度空间的连续依赖性。基于非定常时间分裂数值计算方法和无波动、无自由参数的NND耗散差分格式，建立直接求解气体分子速度分布函数的气体运动论有限差分数值方法。推广应用改进的Gauss-Hermite无穷积分法和华罗庚－王元提出的以单和逼近重积分的黄金分割数论积分方法等，对离散速度空间进行宏观取矩获取物理空间各点的气体流动参数，由此发展一套从稀薄流到连续流各流域统一的气体运动论数值算法。通过对不同Knudsen数下一维激波管问题、二维圆柱绕流和三维球体绕流的初步数值实验表明文中发展的数值算法是可行的。     

基于卡门尺度和滤波方法的SST方程改进

基于滤波方法和卡门尺度对原始剪切应力输运(shear stress transport, SST)湍流模型进行了改进,提出了一种卡门尺度修正的滤波SST 方法. 湍流多尺度效应必须在分离流场模拟中给予反映,该方法减弱了雷诺平均(Reynolds averaged Navier-Stokes, RANS)方法时间平均特性对于流场脉动量的压迫作用,在流场中引入了大涡模拟(large eddy simulation, LES)方法的亚格子模型,形成一种新型的脱体涡模拟方法(detached eddy simulation,DES)方法;同时,为了降低原始DES方法在网格加密过程中产生网格诱发的雷诺应力损耗,利用卡门尺度对滤波因子进行修正. 平板边界层算例中,卡门尺度对于RANS方法的跟随性远远强于DES方法,在边界层内的速度型和RANS方法吻合很好,而DES方法在加密过程中速度型的鲁棒性较差,说明卡门尺度在有效地保护了边界层内使用RANS求解,降低速度型偏离对数率现象的产生;HGR-01翼型算例证明BY-SST方法可以有效的避免网格诱导分离现象的产生;证明BY-SST方法在分离流动中的精度高于DES方法.     

关于具有剪应力的普通湍流理论发展

本文叙述了周培源1945年发表具有剪应力的普通湍流理论以后,各种湍流模式理论的主要发展概况.      

考虑损伤的内变量黏弹-黏塑性本构方程

基于Rice 不可逆内变量热力学框架,在约束构型空间中讨论材料的蠕变损伤问题. 通过给定具体的余能密度函数和内变量演化方程推导出考虑损伤的内变量黏弹-黏塑性本构方程. 通过模型相似材料单轴蠕变加卸载试验对一维情况下的本构方程进行参数辨识和模型验证,本构方程能很好地描述黏弹性变形和各蠕变阶段.不同的蠕变阶段具有不同的能量耗散特点. 受应力扰动后,不考虑损伤的材料系统能自发趋于热力学平衡态或稳定态. 在考虑损伤的整个蠕变过程中,材料系统先趋于平衡态再背离平衡态发展. 能量耗散率可作为材料系统热力学状态偏离平衡态的测度;能量耗散率的时间导数可用于表征系统的演化趋势;两者的域内积分值可作为结构长期稳定性的评价指标.      

A NEW INTERFACE CAPTURING METHOD BASED ON DOUBLE INTERFACE FUNCTIONS

基于代数重构思想,发展了一种新的双界面函数重构方法,并采用双正弦函数构造了双正弦界面重构方法（double sine interface capturing,DSINC）.为验证不同界面函数对界面捕捉效果的影响,用数值方法求解了可压缩五方程模型,其中对流项的离散采用五阶WENO（weighted essentially non-oscillatory method）格式,时间积分采用三阶Runge——Kutta方法,通量计算分别考虑了HLL和HLLC方法,而状态方程采用Mie-Grüneisen状态方程.在数值计算中,在界面附近,采用DSINC来获得体积分数的重构,而在远离界面的区域采用WENO格式来获得高阶插值状态.相比采用单界面函数的方法,如双曲正切界面重构方法（tangent of hyperbola for interfacecapturing,THINC）,DSINC方法同样具有界面重构算法简单,在程序中添加方便等特点,两者区别在于,DSINC方法在重构过程中未知函数更易于求解,而无需求解复杂的非线性超越方程,这就使其具有易于向多维扩展的能力.一些典型的两相流动问题,如圆形水柱对流问题,两相三波点问题和激波——界面不稳定性问题等被用作不同界面函数对界面捕捉效果的影响对比.对比分析发现,DSINC与THINC在界面捕捉效果上大致保持一致,并在计算中表现出了较好的稳定性.双界面函数重构思想可以为多相流动界面的代数重构提供了一种新的思路.     

AN EXPERIMENTAL METHOD BASED ON ENERGY RELEASE RATE OF INTERFACIAL FRACTURE

基于对断裂力学常用实验方法的研究,结合界面断裂问题的特殊性,以断裂力学为理论基础,通过能量释放率建立了界面断裂测量的实验分析方法,并且利用文献中的实验数据进行了验证,取得了良好的一致性.该方法通过测量试验件的载荷位移关系,利用裂纹扩展过程中的能量变化关系得到该裂纹长度下的临界能量释放率;在此基础上,根据试件的阻抗能量曲线预测结构的最大承载能力.该方法以能量释放率为理论基础,为界面裂纹的强度分析提供了合理的手段,基于能量角度建立的实验分析方法也具有良好的实用性和适用性.     

基于步态切换的欠驱动双足机器人控制方法

由于高维、非线性、欠驱动等特点, 3-D双足机器人的稳定性控制依然是一个研究难点. 一些传统的控制方法, 如基于事件的反馈控制方法和PD控制方法, 抗扰动能力较弱, 鲁棒性较差. 通过观察, 人类受到外部扰动影响时, 会通过调整步态重新获得稳定性,相较之下仅依靠一个步态获得的稳定性是有限的. 受此启发, 本文针对上述问题提出一种基于步态切换的欠驱动3-D双足机器人控制方法. 首先, 以能耗最少为优化目标, 通过非线性优化方法预先设计多组不同步长、步速的步态作为参考步态, 以构建一个步态库; 然后, 通过综合考虑步态切换过程中的稳定性与能效, 建立了多目标步态切换函数; 最后, 将该步态切换函数作为优化目标, 并求解该最小化问题获得下一步的参考步态, 从而实现步态切换, 达到使用步态库?多轨迹方法来提高鲁棒性的目的. 在仿真实验中运用该步态切换控制方法, 欠驱动3-D双足机器人可实现相对高度在[-20, 20] mm内随机变化的不平整地面上行走, 而仅采用单步态控制策略则无法克服这样的外部扰动, 从而说明了基于步态切换的欠驱动双足机器人控制方法的有效性.