基于优化有限差分和混合吸收边界条件的三维VTI介质声波和弹性波数值模拟

传统有限差分系数是通过泰勒级数展开求取的,这样导致所计算的频散曲线在大波数区域会产生较强的数值误差.针对二阶空间偏导数的显式有限差分离散,本文发展了一种新的优化差分系数方法:首先将泰勒级数展开与多点采样方法结合应用于空间频散关系,基于最大范数建立直观有效的优化目标函数,采用Remez算法求解该目标函数,从而获得最优化差分系数.利用优化有限差分方法求解三维垂直对称轴横向各向同性(VTI)介质中的声波和弹性波方程.另外,本文将二维混合吸收边界条件推广到三维VTI介质中,用于吸收人工截断边界反射;基于各向异性特征,合理调整了边界区域的速度值来提高吸收效果.考虑到三维情况下计算效率的问题,本文波场外推过程中采用图形处理器(GPU)取代传统的中央处理器(CPU).数值精度分析表明,相比较于传统的泰勒级数展开方法,优化有限差分方法在大波数区域对频散误差的压制效果更明显.在三维均匀和修改的Hess VTI模型中的数值模拟实验证明了本文方法具有更高的精度与效率,混合吸收边界条件在三维VTI介质中具有良好的边界吸收效果.

小离地高度下翼型气动特性计算模型的研究

在地面效应影响下,翼型的气动性能会发生很大变化.研究过程中发现离地高度较小,翼型与地面的距离与壁面边界层厚度处于同一数量级时,使用高雷诺数湍流模型和标准壁面函数计算所得的升力系数会严重偏离实验值.通过改进湍流模型和壁面边界,即用低雷诺数湍流模型结合混合壁面边界的处理方法解决了以上问题.计算结果与实验值的比较说明,该方法在离地高度较小时是十分有效和可靠的,而随着离地高度的增大,该方法的计算误差也将随之增大.所以在计算翼型的近地气动特性时,应先比较翼型表面与地面之间的间距和壁面边界层厚度,再确定采用相应的计算模型来提高计算的准确度.

污水土地处理系统重塑土渗透性及混合比的确定方法

污水土地处理系统的渗透性是能否实现污水处理设计要求的关键。本文采用地下水数值模拟方法刻画了污水土地处理系统复杂的结构和边界条件，准确地得到了能够满足污水处理设计要求的污水土地处理系统的渗透系数，并通过不同砂-土混合比例的重塑土渗透试验，研究了能够满足污水处理设计要求的重塑土的砂-土混合比，以为充分利用原场地黏土资源提供依据，也为类似研究和设计提供新的思路和方法。

非定常预处理方法在倾转旋翼飞行器悬停状态气动干扰模拟中的应用

针对倾转旋翼飞行器的多种速度尺度并存的复杂流动问题,发展了能够统一求解同时包含不可压缩流动和可压缩流动的非定常预处理方法和一种考虑了流场物理信息传播方向的预处理远场边界条件。结合适用于倾转旋翼飞行器的复杂外形及大位移问题的非结构混合网格技术和重叠插值技术,建立了倾转旋翼飞行器的数值模拟方法。在此基础上,首先通过悬停算例验证了数值方法的有效性,确定了选择预处理截断参数的原则。然后探究了倾转旋翼飞行器直升机模式悬停状态的尾迹涡结构演化过程,研究了"喷泉效应"和"地面效应"对倾转旋翼气动性能的影响。计算结果表明,旋翼尾迹涡的耗散在不同方位角存在较大差异。"地面效应"抑制了旋翼尾迹的发展,提高了旋翼整个方位角区域的升力系数,而不仅限于机翼上方方位角区域。"地面效应"作用强于"喷泉效应"作用,使旋翼相对于孤立旋翼拉力系数略微增大了5.26%,但同时整个倾转旋翼飞行器的总的向上的力只有孤立旋翼拉力的82.46%。结论表明预处理方法和预处理远场边界条件能够较好地模拟倾转旋翼飞行器悬停状态的干扰流场,兼顾了数值计算的耗散性和稳定性。气动干扰对倾转旋翼飞行器整体气动性能有重大影响。

基于无散平滑的4D flow MRI流场优化

目的时间解析的三维相位差核磁共振图像速度编码技术(4D flow MRI)在大血管临床诊断方面有重要应用价值。但是其测量的速度场通常包含许多噪音和坏点,需要利用后处理方法来优化并进行血流动力学的分析。方法提出一种带壁面边界条件的无散平滑(DFSwt)方法来处理4D flow MRI的速度场数据,增强具有物理约束的速度场质量。