ALE方法中一种新的二阶保界守恒重映算法

在用拉格朗日格式求解流体力学问题时,随着时间的推进,计算网格会扭曲变形,影响格式精度,甚至导致计算中断。因此,要在网格变形较大时进行网格重分和物理量重映,以保证网格质量和格式精度。针对间断有限元方法求解流体力学问题的二阶拉格朗日格式,给出了一种守恒重映算法。该重映算法包括两步:第一步是用已有重映方法计算新网格上的单元平均值,并用相应修补算法对单元平均值进行调整,保证单元平均值的保界性;第二步是由已得到的新单元平均值重构出新网格上分片一次多项式,再使用Van Leer限制器对新网格上的梯度进行限制,使之不出现新的极值。最后用数值算例表明了该重映算法的保界性和二阶收敛性。     

基于人群搜索算法的Dykas汽轮机叶栅型线重构

针对汽轮机叶片智能优化设计中的叶片参数化建模与型线重构问题,首先基于三次Bezier曲线推导了叶片型线的坐标方程,根据叶片基本几何参数实现了汽轮机叶片型线的参数化表达。在构建的参数化模型的基础上,根据已知的叶型坐标点数据,采用人群搜索算法(SOA)重构实际汽轮机叶片型线,并通过叶型重构算例对SOA和粒子群算法(PSO)进行了对比分析。结果表明：基于Bezier曲线构建的参数化型线光滑性好质量高,所提方法可以任意修改叶型几何参数并再生模型,实用性较强且效率较高。基于SOA优化算法的叶片型线重构方法可以准确重构Dykas汽轮机叶片型线,且SOA重构叶型的方法收敛速度快、收敛结果更稳定,为汽轮机叶片优化设计奠定基础。     

非线性系统中多传感器目标跟踪融合算法研究

 研究了在非线性系统中 ,基于转换坐标卡尔曼滤波器的多传感器目标跟踪融合算法。通过分析得出 :在非线性系统的多传感器目标跟踪中 ,基于转换坐标卡尔曼滤波器 ( CMKF)的分布融合估计基本可以重构中心融合估计。仿真实验也证明了此结论。由此可见分布的 CMKFA是非线性系统中较优的分布融合算法     

涡扇发动机气路传感器故障诊断

为了实现对某涡扇发动机传感器故障的在线诊断,提出并设计了1种基于在线贯序极端学习机的故障诊断算法。其核心思想是在定位某传感器故障后,在线建立针对该故障传感器"预学习"的信号重构算法,解决多故障混叠问题。在线信号重构算法以泛化能力指标为判定条件,利用选择策略对算法网络权值进行选择性更新,提高了故障诊断系统的实时性。以某型涡扇发动机为对象开展了传感器故障诊断与重构仿真,结果表明:该算法能够对发动机单、双传感器故障进行准确地诊断与信号重构,且具有良好的实时性。     

面向再入目标跟踪的估计与辨识联合优化算法

准确的弹道系数辨识和精确的目标状态估计是再入目标高精度跟踪与高可靠识别的关键。一方面,状态估计的误差会造成模型参数(弹道系数)的辨识风险;另一方面,模型参数的辨识偏差又会导致模型失配从而降低目标状态的估计精度。因此,需要实现再入目标的状态估计和参数辨识的联合优化。针对再入目标弹道系数未知情形,提出了一种基于期望最大化(EM)框架并采用粒子滤波(PF)平滑器实现的PF-EM联合优化算法。在E步基于粒子平滑器得到目标状态的后验平滑估计,M步采用数值优化算法更新上一次迭代的弹道系数,通过E步和M步的不断迭代,以保证状态估计和弹道系数辨识的一致性。算法仿真对比表明:所提算法的状态估计和参数辨识精度均优于传统的状态增广算法。     

基于比拟理论计算圆柱/翼型的气动噪声

采用计算流体力学(CFD)与"声比拟"相结合的方法,通过两个步骤模拟流动的声学远场。第一步,在包含所有声源的近场区域内,通过求解URANS方程获得控制面处的非定常流场参数;第二步,采用基于可穿透数据面的FW-H方法模拟声学远场。该方法与传统的半经验方法相比,计算精确,易于工程实现,并且可以计算非线性噪声。根据气动声学计算积分公式,构造了高阶的时间离散格式,采用高斯积分公式计算圆柱/翼型的气动噪声。在时域和频域上计算了观测点处的声压及声压级随时间步数的变化情况,其结果与国外实验结果对比取得了较好的一致性,此外研究了翼型尾流部分界面对计算精度的影响。     

基于QSS的自适应多步校正算法及其在柔性航天器动力学中的应用

量化状态系统（QSS）算法是一种基于状态变量离散化的数值积分方法,该方法与基于时间离散的传统方法显著不同,QSS通过计算状态变量每次跃迁所需的时间来推进下一步的积分。在求解非刚性常微分方程时,QSS算法比传统算法更具优势,但它不适合求解刚性问题。为此提出一种基于QSS的自适应多步校正算法（AMCQSS）,该算法以QSS为基础、结合隐式多步法思想,在计算过程中可以自适应选择二步法或三步法,以有效提高求解刚性问题的精度和效率。通过对柔性航天器动力学的仿真求解,验证了算法的可行性。将该算法与ODE23tb、ODE15s、ODE45以及QSS等算法进行对比,结果表明AMCQSS算法既能保证求解的效率及精度,又具有较好的收敛性和稳定性。     

超声速翼型气动优化设计

首先分析了几何外形和相对厚度对超声速翼型气动特性的影响。基于遗传算法(GA)和气动力快速工程算法,对于相对厚度为3.5%的多边形翼型进行优化设计,多边形翼型的优化外形趋于四边形,最大厚度点后移到翼型弦线的60%左右,随着迎角或者马赫数增大下翼面会变薄,上翼面变厚,最大厚度点相应稍有后移。对于相对厚度为4%的双圆弧翼型,采用两步优化设计方法,第1步优化结合基于B样条的类别形状函数变换(CST)参数化方法与小波分解方法,实现几何外形的局部控制与光顺处理,并且采用本征正交分解(POD)代理模型降低优化过程中流场计算的工作量;第2步优化采用基于Navier-Stokes方程的最速下降法(SDA),修正第1步优化中代理模型和小波光顺引入的误差;优化设计得到的翼型近似为四边形,其相对厚度最大点后移到翼型弦线的60%~65%处,升阻比可以提高7%。     

二维非结构网格的一个TVD型有限体积方法

考虑双曲型守恒律方程,对二维非结构三角形网格给出一种 TVD型有限体积方法,主要思想是在一阶单调格式的基础上,在每一个单元上对变量作单调线性重构函数,时间离散采用二阶 Runge-Kutta方法。通过计算分析了该方法的精度,对平面激波反射和空穴流动的计算结果表明该方法是成功的。     

运载火箭动力故障下的自主救援轨道规划

针对运载火箭上升段出现动力下降故障的问题,提出了一种基于二阶锥规划方法的在线救援轨道规划策略与飞行轨迹重构算法。根据故障大小设置不同的终端条件,规划合适的救援轨道;同时为了方便求解,设计了新的控制量,然后采用逐次凸化、线性化和梯形离散等方法,将原始规划问题转化为序列迭代求解的二阶锥规划问题。仿真结果表明,所提策略及算法能够满足动力系统故障后的在线救援轨道规划和飞行轨迹重构需求。     

基于Euler／N—S方程的跨音速非线性静气动弹性问题研究

在C-H网格的基础上,采用Jameson的中心差分有限体积法求解Euler/N-S方程,采用结构影响系数法计算结构的弹性变形,用三角元面积加权法和常体积转换法(CVT)实现流固耦合,计算了弹性后掠机翼在跨音速状态的静气动弹性问题     

A Preconditioned Multigrid Method for Efficient Simulation of Three-dimensional Compressible and Incompressible Flows

将矩阵预处理方法与多重网格方法结合,发展了适用于机翼及复杂外形气动数值优化设计的高效数值模拟方法和程序。在不改变定常流动解的前提下,对 Navier-Stokes 方程的时间导数项实施 Choi-Merkle 矩阵预处理,从而改善可压缩控制方程在低速情况下的系统刚性。Choi-Merkle 预处理最初用于二维低马赫数粘性流动数值模拟,本文将其推广应用到计算三维流动,并针对提高跨音速流动计算收敛性进行了相应的改进。预处理后控制方程采用中心格式有限体积进行空间离散,Runge-Kutta 混合多步方法进行显式时间推进求解,并应用 FAS 多重网格方法加速收敛。采用所发展的方法和程序数值模拟了绕 M6 机翼和 M100 机翼的低速和跨音速流动(马赫数从0.01 到 0.839)以及绕 F4 翼身组合体流动。计算结果与实验值吻合良好,验证了所发展的方法的正确性。算例研究表明所发展方法较大幅度地提高了三维流动数值模拟的效率,且同时适用于可压和不可压流动计算,在机翼及复杂外形气动优化设计方法具有应用前景     

使用高阶逆风通量差分裂格式的悬停旋翼流场数值模拟

为减少Jameson二阶中心差分有限体积法导致的旋翼尾迹的数值耗散,将三阶逆风格式(MUSCL)与通量差分裂方法相结合,建立了一个三维雷诺平均N-S方程数值模拟悬停状态旋翼流场的方法。为了充分考虑旋翼尾迹对流场的影响,采用周期性边界条件和由动量理论导出的远场边界条件。为进一步减少尾迹数值耗散和便于添加上述边界条件,采用了嵌套网格方法。然后,进行了算例计算,给出了桨叶表面的压强分布,与可得到的试验数据及二阶中心差分方法的计算结果进行了对比,并针对几种不同桨尖形状的旋翼悬停流场进行计算,数值结果显示:后掠桨尖可减弱超临界流动。此外,还计算和分析了旋翼下方不同位置上的涡量分布。计算结果表明,本文方法能够有效地减少旋翼尾迹的数值耗散。      

一种改进的非定常激波装配算法

基于非结构动网格技术和格心型有限体积方法,提出一种改进的非定常激波装配算法,进一步拓展了其在包含有运动激波的非定常流场的应用范围。首先,针对激波在直/曲壁面传播这类问题,分别建立了壁面间断节点的运动模型;其次,为保证激波在大范围运动时装配阵面不产生失真,基于Bézier曲线拟合方法实现了间断节点分布的自动重构;接着,通过嵌入局部网格自动重构模块,提高了算法的计算效率和自动化程度;最后,对于激波相交点的运动,设计了一种根据位移推算速度的方法进行装配。数值算例表明,所提算法能够有效地处理激波传播问题,相比激波捕捉方法可以提取更多的流场信息,同时可以获得流场间断更加直观清晰的图谱。     

利用多状态试车数据修正发动机部件特性

针对经典的航空发动机部件特性修正因子法的不足,提出一种于多状态试车数据的修正方法.利用发动机部件匹配模型方程以及试车参数计算值和实测值的差值组成目标方程组,采用改进的免疫粒子群算法优化求解获得部件特性修正系数.以多个不同状态下的修正系数为基础,利用滑动最小二乘方法拟合修正系数曲面,进而得到修正的部件特性图.数值试验表明:综合利用不同状态下发动机试车数据进行部件特性修正,克服了传统方法的不足.提高了修正部件特性在整个工作范围的精度.      

三维可压和不可压缩流动数值模拟的一种预处理多重网格方法研究(英文)

将矩阵预处理方法与多重网格方法结合,发展了适用于机翼及复杂外形气动数值优化设计的高效数值模拟方法和程序。在不改变定常流动解的前提下,对 Navier-Stokes 方程的时间导数项实施 Choi-Merkle 矩阵预处理,从而改善可压缩控制方程在低速情况下的系统刚性。Choi-Merkle 预处理最初用于二维低马赫数粘性流动数值模拟,本文将其推广应用到计算三维流动,并针对提高跨音速流动计算收敛性进行了相应的改进。预处理后控制方程采用中心格式有限体积进行空间离散,Runge-Kutta 混合多步方法进行显式时间推进求解,并应用 FAS 多重网格方法加速收敛。采用所发展的方法和程序数值模拟了绕 M6 机翼和 M100 机翼的低速和跨音速流动(马赫数从0.01 到 0.839)以及绕 F4 翼身组合体流动。计算结果与实验值吻合良好,验证了所发展的方法的正确性。算例研究表明所发展方法较大幅度地提高了三维流动数值模拟的效率,且同时适用于可压和不可压流动计算,在机翼及复杂外形气动优化设计方法具有应用前景     

基于Kinect三维重建的地下井室可视化方法研究

针对现有地下井室病害探测与维护方法的不足,提出了一种基于Kinect三维重建的地下井室可视化方法,以实现其三维可视化探测与维护.还提出了一种基于多项式曲面拟合的Kinect深度测量误差修正方法,利用联合双边滤波算法对深度图像数据进行预处理;结合SIFT特征匹配和改进的RANSAC算法获取相邻点云间的初始位姿,并利用基于邻域特征的ICP算法进一步实现不同视角点云的精确配准,从而获取全局一致的稠密三维点云.最后,在三维稠密点云的基础上进行曲面重建和纹理贴图,以实现地下井室真实三维重建.实验结果表明:所提方法可有效修正Kinect深度相机的深度测量误差,在0.5～4.5m的测距范围内,其三维重建精度可达2cm;在4.5～7m的测距范围内,精度也可以保持在4.5cm以内.所提重建方法可实现地下井室真实场景的三维可视化,为地下井室的探测和维护提供了技术支持.     

一种基于N—S方程的CFD／CSD耦合计算方法

对于气动弹性问题,在跨音速和大变形情况下,模态叠加法失效,为此需要探讨全场求解耦合方程的新方法。利用耦合计算流体动力学和计算结构动力学（CFD／CSD）方法,数值模拟了机翼的气弹响应问题。气动力的求解采用三维Navier—Stokes方程的有限体积法,实现了一种多块结构网格的动网格方法生成气动网格,采用有限元方法计算应...     

高速列车进入带缓冲结构隧道的压力变化研究（I）

给出了列车穿越带有缓冲结构的隧道压力变化的三维粘性流场数值模拟过程，控制方程为三维粘性、可压缩、非定常流的N．S方程，空间离散采用了中心有限体积法格式，时间采用预处理二阶精度多步后差分格式进行离散，对列车与隧道之间的相对运动采用移动网格技术处理。对不同的缓冲结构缓解隧道内瞬变压力及压力梯度的作用进行了研究。研究结果表明，缓冲结构的设置能够有效地降低隧道内的压力和压力梯度的最大值，其原因在于缓冲结构延长了压缩波压力上升的时间，降低列车突入隧道时所形成的最大压力梯度；另一方面由于压缩波在缓冲结构和列车、隧道之间多次的反射，也降低了压力峰值。     

高速列车进入带缓冲结构隧道的压力变化研究(Ⅰ)

给出了列车穿越带有缓冲结构的隧道压力变化的三维粘性流场数值模拟过程,控制方程为三维粘性、可压缩、非定常流的N-S方程,空间离散采用了中心有限体积法格式,时间采用预处理二阶精度多步后差分格式进行离散,对列车与隧道之间的相对运动采用移动网格技术处理。对不同的缓冲结构缓解隧道内瞬变压力及压力梯度的作用进行了研究。研究结果表明,缓冲结构的设置能够有效地降低隧道内的压力和压力梯度的最大值,其原因在于缓冲结构延长了压缩波压力上升的时间,降低列车突入隧道时所形成的最大压力梯度;另一方面由于压缩波在缓冲结构和列车、隧道之间多次的反射,也降低了压力峰值。