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

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

基于多尺度方法的平动圆柱贮箱航天器刚——液耦合动力学研究

以充液航天器为工程背景,借助多尺度方法研究刚--液耦合动力学系统非线性动力学特性.利用多维模态方法,将描述横向外激励下圆柱贮箱中液体非线性晃动的自由边界问题转换为液体模态系数相互耦合的有限维非线性常微分方程组.推导液体晃动产生的作用于贮箱壁的晃动力和晃动力矩的解析表达式,进而建立航天器刚体部分平动和液体晃动耦合的非线性动力学方程组.应用多尺度方法对刚--液耦合系统的动力学特性进行解析分析,通过固有频率的特征方程求解耦合系统固有频率,推导外激励频率接近耦合系统第一阶固有频率时液体晃动稳态解的幅值频率响应方程.结合数值方法,研究了液体晃动稳态解的幅值频率响应曲线和激励--幅值响应曲线.结果表明, 随充液比变化,液体晃动稳态解的幅值频率响应曲线会发生软、硬弹簧特性转换现象和"跳跃"现象;幅值频率响应曲线的软、硬弹簧特性转换点受重力加速度和弹簧刚度系数影响;以上所得研究结果表明,考虑非线性效应时的刚--液耦合系统动力学特性与传统的线性系统模型所显示的动力学特性具有本质区别.本文的研究工作对进一步分析充液航天器刚--液耦合非线性动力学特性具有重要参考价值.      

不同加载速率下砂岩弯曲破坏的细观机理

岩石细观破裂形貌是岩石破坏机制的重要反映，为研究不同加载速率对砂岩弯曲破坏的影响，通过三点弯曲实验和扫描电镜方法，对某煤矿关键层砂岩弯曲破断裂纹细观形态以及裂纹的自相似性进行了研究。选取6个不同加载速率对岩样进行三点弯曲实验，观察其宏观断裂情况，并利用扫描电镜对弯曲断裂面表面裂纹细观结构进行观察，并拍摄不同倍数下的扫描电镜图片。对图片进行图像处理后得到砂岩弯曲断裂破坏细观裂纹信息，并计算得到微裂纹的分形盒维数值。结果显示:随着加载速率的提高，砂岩穿晶断裂的比例也随之升高，裂纹分形维数亦随着加载速率的增大而增加，同时，分形维数还与弯曲断裂破坏荷载和抗弯强度成正比。可见，加载速率对断裂方式有一定的影响，且加载速率越大断裂所需的破坏能越大，裂纹分布越广，表明开采速度与岩爆等岩体动力灾变有密切关系。     

基于双轴转台的捷联惯性导航系统8位置系统级标定方法(英文)

惯性测量单元(IMU)是捷联惯性导航系统的核心设备,IMU的精度直接影响捷联惯性导航系统的精度。标定精度和速度是评价IMU标定方法的主要指标。传统标定方法是先进行温度补偿,再利用高精度三轴转台进行位置和速率实验,共需耗费约4天时间。虽然能够保证IMU的标定精度,但该方法复杂耗时。为改善复杂耗时的传统标定方法并保证标定精度,以便于批量生产,提出了基于双轴转台的8位置系统级标定方法。该方法基于27维状态向量的卡尔曼滤波器,以导航速度误差为观测量。该方法基于带温箱双轴转台设计,每个温度点需40 min进行标定。仿真和试验结果表明,该方法能够标定21个IMU误差参数,且该标定方法与传统标定方法相比标定时间更短,导航精度更高。     

非Hermitian正定线性方程组的外推的HSS迭代方法

为了高效地求解大型稀疏非Hermitian正定线性方程组，在白中治、Golub和Ng提出的Hermitian和反Hermitian分裂（HSS）迭代法的基础上，通过引入新的参数并结合迭代法的松弛技术，对HSS迭代方法进行加速，提出了一种新的外推的HSS迭代方法（EHSS），并研究了该方法的收敛性.数值例子表明：通过参数值的选择，新方法比HSS方法具有更快的收敛速度和更少的迭代次数，选择了合适的参数值后，可以提高HSS方法的收敛效率.     

三次Gauss和及二项指数和的混合均值

我们用三角和的性质研究一类三次Gauss和与两项指数和混合均值的计算问题,并给出一个精确的计算公式.     

Euler-Bernoulli梁的高阶二次摄动解及收敛性讨论

首次用解析的方式给出了Euler-Bernoulli梁后屈曲与非线性弯曲问题的高阶二次摄动解答.假定梁的中线不可伸长，用精确曲率公式与能量变分原理导出了非线性Euler-Bernoulli梁的模型.通过与精确解或高阶摄动解的比较，讨论了二次摄动解答的收敛性及适用域.得到主要结论如下:低阶摄动解适用于描述梁的初始后屈曲阶段及初始非线性弯曲阶段；更高阶次的摄动解适用于描述梁的深度后屈曲以及深度非线性弯曲.从这个意义上去说，该文不仅仅指出某些文献上的部分结果不精确是由于摄动解答超出了其特定的适用域，并且还进一步发展与完善了二次摄动法.     

一维正方准晶中裂纹和刚性线夹杂的封闭形式解（英文）

通过引入广义复变函数方法,研究含裂纹和刚性体夹杂物的反平面模型的一维正方准晶问题.对于一维正方准晶,考虑周期平面为(x_1, x_2),含有宏观裂纹或刚性线夹杂,具有准周期x_3方向的原子结构存在相位位移,本文重点研究相位位移对相关物理量的影响.利用广义复变函数方法,将这两个模型简化为Riemann-Hilbert问题,得到反平面的声子场与相位场的封闭解.同时求得声子场和相位场的应力强度因子的显式解,这在断裂力学和工程领域具有广泛的应用价值.结果表明,反平面情形下,含裂纹和刚性体夹杂物的声子和相位的应力强度因子,与声子场和相位场的耦合无关.     

含分布时滞的时滞微分系统多步龙格-库塔方法的时滞相关稳定性

本文研究了一类含分布时滞的时滞微分系统的多步龙格-库塔方法的稳定性.基于辐角原理，本文给出了多步龙格-库塔方法弱时滞相关稳定性的充分条件，并通过数值算例验证了理论结果的有效性.