计算地球流体力学若干新进展

本文简要介绍了有关计算地球流体力学问题研究的若干新进展,其中包括显式平方守恒格式,经济算法,辛算子法,泥沙冲积和三角洲发育的数值模拟等,从中也可以看到计算地球流体力学的发展是具有重要理论意义和实用价值的。     

高精度差分求解气动方程的几个问题

本文探讨了发展高精度格式的必要性，研究了高精度格式与网格以及熵增条件的关系，并发展一种半离散化的空间为三阶精度的格式。从此半离散化的格式出发，可建立多步显式格式及隐式格式。模型问题的计算表明，该三阶精度的格式具有好的精度，且激波附近基本上没有波动。     

求解多维欧拉方程的二阶旋转输运格式

提出了一个基于旋转近似Riemann求解器的二阶精度迎风型有限体积方法.不同于"网格相关"(Gridaligned)的有限体积方法或者维数分裂的有限差分方法,本格式在求解Riemann问题时不依赖网格的方向,而是沿具有一定的物理意义的两个方向(称为迎风方向).我们发现当迎风方向取为控制体界面两侧速度差矢量方向(及与之正交的方向)时,该格式能够完全消除基于Riemann求解器的通量差分裂格式存在的激波不稳定或者所谓"红斑"(carbuncle)现象.为了提高格式的时间和空间精度,我们通过在控制体界面处求解线化的广义Riemann问题的方法体现输运过程对通量计算的影响.这种方案,使得我们有可能以此为基础,构造真正多维的有限体积型迎风格式.     

AUSM+-up格式在无网格算法中的推广

将AUSM -up格式推广运用到了无网格算法中.在点云离散的基础上,运用曲面拟合和AUSM -up格式求得数值通量;采用四步龙格-库塔方法进行时间推进,并引入当地时间步长和残值光顺等加速收敛措施.为了验证本文方法的计算精度和鲁棒性,对NACA0012翼型跨声速流动、某三段翼型低速绕流、平行错位NACA0012双翼超声速、高超声速流动进行了数值模拟.     

不同壁温及差分格式下超燃冲压发动机的仿真

针对圆截面超燃冲压发动机相关实验进行数值仿真,采用了不同的壁面温度条件和差分格式,并探究其对仿真结果的影响规律。结果表明采用常用的绝热壁面边界条件会导致激波边界层交互位置处的回流区变大,进而导致流道收缩,壁面压力升高。采用高精度格式不仅不能更为精确地计算回流区的大小,反而会使回流区面积变得更大。通过改变壁面边界条件,采用300 K的等温壁面,可以使激波边界层交互位置处的回流区变得极小,使仿真结果与实验结果十分吻合。计算结果表明,针对当前构型与实验条件,盲目提高数值格式的精度并不一定会得到更好的结果,相反可能会使仿真与实验更加偏离,而适当修改壁面温度条件,即使数值格式精度较低,仍然可以得到很好的仿真结果。最后,针对与实验结果吻合的仿真结果,分析了圆截面超燃冲压发动机的流场特性,重点研究了凹腔剪切层及其质量交换特性。     

抗扰动复合非线性伺服控制器的设计与应用

针对典型的电机伺服系统,提出了一种鲁棒复合非线性伺服控制器的离散域设计方案。把系统的扰动和不确定性归结为一个斜坡信号(其变化率恒定),设计一个降维线性扩展状态观测器,对系统未测量状态和未知扰动加以估计。把设计的控制律应用于永磁同步伺服电机,先在MATLAB上进行仿真分析,随后基于TMS320F28335DSC进行试验测试。结果表明系统在各种类型扰动作用下,对目标位置都能实现快速、平稳和准确的跟踪,具有较好的鲁棒性。     

带压缩性修正的离散涡方法在结冰数值模拟中的应用

离散涡方法（DVM）是一种无网格的涡运动算法,适用于解决易产生分离流的非定常问题,将其应用于结冰过程中的流场求解,在有效模拟分离流动的同时能避免冰形尖角对网格质量的影响。但离散涡方法基于不可压 N-S 方程,无法应用于预测可压缩流动下的结冰过程。本文在离散涡方法的基础上添加普朗特—格劳尔特压缩性修正,进行基于离散涡方法的可压缩流动下的数值模拟,并将其应用于翼型结冰预测;对流场分布、结冰冰形和结冰模型计算过程等仿真结果进行对比和分析。结果表明：引入压缩性修正后的离散涡方法能较好地模拟可压缩流动,与实验值相比,基于该方法得到的结冰数值模拟结果符合良好,对结冰数值模拟在工程上的应用提供了一定的参考。     

矩的显式积分算法研究与应用

矩的求解通常被用于求解有限元、体积、惯性矩等问题中.基于矩的叠加性,首先给出了在三维空间中计算域的离散方式,并推导了矩的显式积分公式,随后将其推广到n维空间中,该表达式易于在计算机上实现;设计了矩的并行计算算法,并通过Fortran和Python混编的方式,实现了矩的并行计算;对多重精度下的样例数据给出了一个算例,实现了零阶矩和二阶矩的计算,并和串行算法、逐次降维算法作出比较,进行了效率分析和误差分析.结果显示,矩的显式积分并行计算算法易于程序实现,并且在效率上高于串行算法,能够很容易推广到高维空间,该算法具有高度可并行性,误差主要来自计算域离散.      

高精度激波捕捉格式的性能分析

在数值模拟激波-湍流干扰问题时,为了精确模拟湍流脉动的小尺度结构,数值格式应尽量减小自身的数值耗散.但是为了稳定、无震荡地捕捉激波,数值格式又必须具有一定程度的数值耗散.基于上述这2种相互矛盾的要求,采用无黏泰勒格林涡、可压缩湍流等经典算例对目前广受好评的几种高精度格式进行了系统地对比.一系列的计算分析表明,当前工程中广泛应用的高精度格式虽然都能一定程度地满足上述2种要求,但与理想目标仍有较大的差距.