HyperFLOW软件非结构网格亚跨声速湍流模拟的验证与确认

计算流体力学(computational fluid dynamics,CFD)数值模拟在航空航天等领域发挥越来越重要的作用,然而CFD数值模拟结果的可信度仍然需要通过不断地验证与确认来提高.本文给出了从制造解精度测试、简单到复杂外形湍流模拟网格收敛性研究等三个方面开展CFD软件验证与确认的方法,并对自主研发的CFD软件平台HyperFLOW在非结构网格上模拟亚跨声速湍流问题的能力进行了验证与确认.首先通过基于Euler方程和标量扩散方程的制造解精度测试,分别验证了HyperFLOW在非结构网格上对Euler方程和黏性项的求解精度,结果表明其能够在任意非结构网格上达到设计的二阶精度. 其次,通过NASATurbulence Modeling Resource中的湍流平板、二维翼型近尾迹流动、二维Bump等几个典型的亚声速湍流算例的网格收敛性研究,量化考察了数值结果的观测精度阶和网格收敛性指数,并与国外知名CFD解算器CFL3D,FUN3D的计算结果进行了对比,验证了HyperFLOW对简单湍流问题的模拟能力,且具有良好的网格收敛性和计算精度(阶). 最后,通过NASA CommonResearchModel标模定升力系数的网格收敛性研究和升阻极曲线预测,验证了软件在复杂外形亚跨声速湍流流动数值模拟中也具有良好的可信度.      

砂土类材料离散元细观参数的自动识别

准确的细观参数是采用离散元进行合理模拟的前提,基于最优化理论,本文提出了一种砂土类材料离散元细观参数的自动标定方法:采用PFC离散元软件模拟砂土的双轴压缩试验,以使PFC数值模拟的应力–应变曲线与某标准砂土室内双轴压缩实验结果的误差最小为优化目标,采用模式搜索法实现砂土类材料离散元细观参数的自动标定。本文提出的离散元细观参数自动标定方法具有识别效率快,精度高的优点。     

具有恒定冲击载荷的梯度泡沫金属材料设计

多胞材料可通过大变形大量地吸收冲击能量，引入密度梯度可进一步提高其耐撞性。梯度多胞材料的宏观力学响应对材料密度分布极为敏感，不同类型的细观构型的影响也极为不同。已有的研究工作主要局限在对给定的密度梯度分析其动态响应，较少对耐撞性设计方法进行研究。本文针对梯度闭孔泡沫金属材料，基于非线性塑性冲击波模型发展了耐撞性反向设计方法，以维持冲击物受载恒定为目标，运用级数法获得了简化模型和渐近解。利用变胞元尺寸法构建了连续梯度变化的三维Voronoi细观有限元模型，并利用ABAQUS/Explicit有限元软件对理论设计进行数值验证。结果表明，反向设计理论简化模型的渐近解对于梯度闭孔泡沫金属材料的耐撞性设计是有效的，所提出的耐撞性设计方法在控制冲击吸能过程和冲击物受载方面具有指导意义。     

基于DIC的爆炸加载下脆性材料裂纹扩展规律的试验研究

通过数字图像相关法（DIC），应用PMMA对爆炸加载条件下脆性材料的裂纹扩展规律进行了试验研究。基于对称性试验模型，实现了裂纹尖端位置和应变场信息的同步记录。以此为基础，通过对比分析获知，主应变场应变值最大点不能作为裂纹尖端的判断依据。并以动态裂纹扩展速度为参量，应用断裂动力学和最小二乘牛顿迭代法，计算出了考虑惯性效应的Ⅰ-Ⅱ混合型裂纹的应力强度因子:K_Ⅰ和K_Ⅱ值会随着裂纹扩展方向改变而发生突变；K_Ⅰ最大值为2.63 MPa·m1/2，最小值为0.89 MPa·m1/2；其整体变化趋势表明，爆炸加载条件下脆性材料裂纹扩展随能量积聚和释放呈循环阶梯式递减发展。     

考虑空化效应的螺旋槽机械密封液膜动力学特性研究

考虑液膜空化效应的影响，研究螺旋槽液体润滑机械密封的动力学特性. 基于液体润滑理论和小扰动法，建立了考虑液膜空化的密封微扰膜压控制方程，采用有限单元法对端面液膜三自由度微扰下的液膜刚度和阻尼系数进行了数值求解，分析了不同参数对液膜密封动力系数的影响. 螺旋槽深度在10 μm左右、槽坝比在0.75左右、槽宽比在0.4左右，螺旋角在9°左右时液膜具有最大的轴向和角向刚度系数. 螺旋槽深度在5 μm左右、槽宽比在0.6左右、螺旋角在20°左右时，两角向交叉阻尼绝对值最大. 初始偏角的存在使密封压力呈现非对称性，从而使两角向动力系数绝对值不再相等. 液膜轴向刚度k_zz在数量级上远大于其余液膜刚度值，液膜轴向阻尼d_zz、角向阻尼d_αα和d_ββ远大于液膜其余阻尼值且随着转速和间隙的增大而减小.      

硅基底复合自组装膜的制备及其微摩擦磨损性能研究

采用共吸附法，在硅基底表面制备3-氨基丙基三乙氧基硅烷(APS)和十二烷基三甲氧基硅烷(WD-10)复合自组装膜. 通过分子动力学模拟不同温度与混合分子在不同比例下的混合体系界面结合能；依据模拟结果，采用正交试验法设计试验方案制备9种不同条件下的自组装膜；采用原子力显微镜、接触角测定仪以及X射线光电子能谱仪对自组装膜的表面形貌、湿润性能和化学成分进行表征分析；利用微摩擦测试仪对自组装膜的微摩擦磨损性能进行性能测试. 结果表明:混合分子成功组装到羟基化硅基底表面，并且当组装温度为25 ℃，组装时间为4 h，组装溶液的pH为6时，自组装膜的质量较好；制备的复合自组装膜由于引起了边界润滑效应，有效减小了试件表面的摩擦磨损，且两种混合分子比例为1:1时自组装膜的减摩特性最佳.      

坐标点分布对三维坐标变换参数估计不确定度的影响（英文）

相机坐标系和目标坐标系中公共特征点的测量精度和相对几何分布对参数估计具有重要影响,直接影响视觉导航精度。由于坐标点的测量精度依赖于测量系统的精度,而公共点的分布选取更容易调整,因此仅讨论公共特征点相对几何分布对平移旋转参数估计精度的影响。首先,采用非线性模型和加权整体最小二乘算法求解变换参数,并采用蒙特卡罗方法归纳总结了不同的特征点分布对参数估计精度的影响。同时,为了工程化应用,将坐标点分布的分散度进行量化作为其权重,提出一种重构协因数矩阵的方法。最后,基于实际测量数据的对比试验表明,在没有先验知识的情况下,所提出的方法可以改善参数的估计精度,即变换参数的后验标准差至少减小6.7%。     

不可压粘性流体浸入边界直接力法及软件实现

浸入边界法通过在N-S方程中施加体积力模拟不可滑移固壁边界及动边界，避免生成复杂贴体网格及动网格，极大地节省了网格建模时间及动网格计算消耗。本文提出一种新型附加体积力简化计算方法，将简化附加体积力以源项形式嵌入动量方程迭代中，通过用户自定义函数对CFD软件FLUENT二次开发，实现了浸入边界法和通用流体力学求解器的耦合计算。通过静止圆柱和动圆柱绕流数值模拟进行了验证，并探讨了插值函数对计算精度的影响。研究表明，通过引入浸入边界模型，能够提高计算效率，并实现结构网格背景下复杂边界和动边界的高效建模。     

基于Hopf-Cole变换的Burgers方程初边值问题的Sinc-Galerkin法

利用Sinc-Galerkin法数值求解Burgers方程的初边值问题。首先,用Hopf-Cole变换将二阶非线性的Burgers方程变换为二阶线性方程,同时把第一类边界条件变为第二类边界条件。时间上的导数采用θ加权格式离散,空间导数采用Sinc-Galerkin法离散,端点处分别引入权函数处理变换后的第二类边界条件。最后,通过数值算例验证了Sinc-Galerkin法的指数收敛性,与精确解相比,本文构造的数值格式精度高,能够有效捕捉激波等物理现象。     

含参广义集值强向量平衡问题的稳定性

本文借助集合极限的性质和弱f-性假设证明了含参广义集值强向量平衡问题解集映射的下半连续性,其方法不同于最近文献(Zhao,2016和Meng,2018).此外,建立了含参广义集值强向量平衡问题解集连通性的充分条件,并举例验证了所得结果的正确性.本文得结果推广和改进了已有文献(Gong,2008,Xu,2009,Chen,2010,Xu,2013和Zhao,2013)中相应结果.