大气水平不均匀条件下电波折射修正法——辐射计法

对于大气水平不均匀性较大的地区,用59型探空仪测量大气剖面,然后用大气球面分层法计算电波折射误差,其精度不会太高。本文提出了用微波辐射计直接测量电波传播路径上大气辐射亮度温度而得到电波折射修正量的方法——辐射计法。由于辐射计法是直接在电波传播路径上进行测量,直接反应大气的空间结构和时变特性,并且不需要施放探空仪,具有方便、灵活、实时、全天候和修正精度高等优点,因此是一种非常实用的高精度电波折射修正方法。     

大气水平不均匀条件下电波折射修正法——辐射计法

对于大气水平不均匀性较大的地区,用５９型探空仪测量大气剖面,然后用大气球面分层法计算电波折射误差,其精度不会太高。本文提出了用微波辐射计直接测量电波传播路径上大气射亮度温度而得到电波折射正量的方法－辐射计法。由于辐射计法是直接在电波传播路径上进行测量,直接反应大气的空间结构和时变特性,并且不需要施放探空仪,具有方便、灵活、实时、全天候和修正精度高等优点,因此是一种非常实用的高精度电波折射修正方法。     

无线传感器网络覆盖新算法研究

目前无线传感器网络对时间空间连续分布的物理量进行监测时存在着节点分布不合理，覆盖效率低，网络计算量大等问题。采用流体力学中数值分析的动态嵌套网格技术，设计了一种新的基于动态嵌套网格技术的无线传感器网络覆盖算法。最后对该算法进行仿真实验，结果表明该算法在相同监测精度时所得到的覆盖效率、覆盖均匀性、能效指标等比传统的覆盖算法有所优化。     

基于域分解的双模式PE

高频近似算法的误差比抛物方程(PE)大,叠加PE计算不能反映同频干扰,对流层电波传播计算域内存在不同的电波传播机制。文章提出了域分解双模式PE,将计算域分解为海洋或地形上的阻抗边界域和空间传播域,分别采用DMFT和SSFT求解。交叠区域实现了电波在两域间的互相传播,并避免了边界截断产生反射。与功率叠加PE相比,反映了多源电波间的矢量叠加,并降低了一半的计算量。在标准大气单源、蒸发波导双源条件下的数值试验验证了方法的正确性和有效性。     

基于混合网格的声爆/气动一体化设计方法研究

探索了基于混合网格的近场预测以及基于波形参数法的远场预测相结合的声爆预测方法,综合利用了结构化网格计算效率高以及非结构网格对复杂外形适应性好的优点,在保证声爆预测精度的基础上提高了计算效率。在此基础上,建立了基于声爆预测方法、响应面模型和Pareto遗传算法的声爆/气动一体化设计方法。分析表明,发动机短舱对超声速升阻特性和声爆特征都有较大影响,利用声爆/气动一体化设计方法对发动机短舱位置进行了一体化设计研究。结果表明,建立的响应面模型的精度完全满足要求,根据不同的设计偏向选择了三种优化方案,优化后方案的性能较原始方案都有不同程度的提高。     

广义Burgers方程声爆传播模型高阶格式离散

声爆预测技术是研制新一代环保型超声速民机的关键技术之一。针对当前基于广义Burgers方程的声爆远场传播模型的数值求解格式精度较低的情况，开展了高阶格式离散方法研究。通过分析该模型中各效应项的物理性质，应用合适的高阶精度格式对各项分别进行离散求解，并采用经典几何声学射线法计算声爆传播路径，实现了对地面声爆波形的精准预测。通过对美国F-5E机型的声爆飞行试验和第二届国际声爆预测研讨会的典型算例进行数值模拟，验证了本文高阶离散方法的可靠性。进一步的结果分析表明，采用高阶离散格式的预测方法相比传统二阶精度预测方法具有明显的低耗散特性，在同等网格数目上可获得更高分辨率的计算结果，且其网格收敛性远优于传统二阶精度格式，同时在计算效率上也有一定优势；在Burgers方程的各效应项中，非线性项的影响更为明显，因而采用高阶离散格式时其优势更为突出；热黏性吸收项对数值结果影响较小，采用高阶离散格式对计算精度的提升效果也不明显，实际计算中仍可以采用传统离散格式，甚至忽略该项的贡献。     

用强耦合RANS方法模拟旋翼悬停流场

在旋转坐标系下,将Spalart-Allmaras(S-A)一方程湍流模型和Reynolds-averaged NavierStokes(RANS)方程耦合成一个新的RANS方程,并发展了基于多块重叠网格的强耦合RANS求解方法,用于直升机旋翼悬停流场的数值模拟.为了提高计算效率,针对多重网格方法在多块重叠网格上实施的困难,提出了一种基于重叠网格的多重网格实施方法.通过对Caradonna-Tung(C-T)和ONERA 7A旋翼悬停算例验证了发展的强耦合RANS方法和基于重叠网格的多重网格实施方法的有效性.研究结果表明:发展的基于重叠网格的多重网格方法有较高的计算效率,3层网格的加速比约为7.7;强耦合RANS法的计算精度明显高于传统的松耦合RANS方法,特别是在与阻力相关性能参数的预测中,强耦合RANS方法的预测结果更加精确.     

基于KZK方程的声爆频域预测法

建立了一种快速预测声爆传播特性的频域方法,基于传统Khohklov-Zabolotskaya-Kuznetsov (KZK)方程,描述声爆沿激波波阵面法线方向的传播。为了验证模型正确性,以NASA TD N-161中试件C为对象,首先基于不同网格,利用计算流体力学(CFD)方法得到超音速流场;将CFD结果作傅立叶变换后代入模型方程,快速求解声爆传播特性。预测结果与NASA实验结果符合很好,研究表明：预测方法能够捕捉声爆的非线性传播,声衍射项使得声场在远场趋于轴对称分布,远距离传播后声能量集中于低频分量。     

非线性壁板颤振计算的子循环预测校正方法研究

为了解决流固耦合问题中弱耦合计算方法存在的滞后误差引起计算发散的问题, 将子循环方法与预测校正方法相结合, 使用有限元方法求解壁板动力响应, 使用Fluent求解受扰动的流场, 使用基于拉伸弹簧与扭转弹簧的动态网格方法更新流场网格, 分析了弹性壁板在超声速气流作用下的非线性动力响应, 结果表明此方法能准确预测壁板颤振的响应特性, 计算稳定性优于弱耦合方法, 可以保持长时间计算结果的稳定性.      

基于Khokhlov-Zabolotskaya-Kuznetsov方程的声爆频域预测法

建立了一种快速预测声爆传播特性的频域方法,基于传统Khokhlov-Zabolotskaya-Kuznetsov(KZK)方程,描述声爆沿激波波阵面法线方向的传播.为了验证模型正确性,以NASA TD N-161中试件C为对象,首先基于不同网格,利用计算流体动力学(CFD)方法得到超声速流场;将CFD结果作傅里叶变换后代入模型方程,快速求解声爆传播特性.预测结果与NASA实验结果符合很好,研究表明:预测方法能够捕捉声爆的非线性传播,声衍射项使得声场在远场趋于轴对称分布,远距离传播后声能量集中于低频分量.      

不规则地形条件下双向DMFT电波传播特性算法研究

 在传统中心差分离散混合傅里叶变换(DMFT)抛物方程模型基础上,针对传统算法计算动态阻抗边界条件时的数值振荡及忽略后向传播与散射造成计算误差增大的问题,提出基于前-后向混合差分方法的递归双向DMFT(TW-DMFT)模型,以一阶前-后向差分方程代替传统的中心差分方程拟合阻抗边界条件,以递归双向传播模型改进传统单向预测模型.仿真分析了双向DMFT模型计算不同边界及媒介条件下单刃峰和实际地形的电波传播分布特性,与传统算法及几何绕射理论(GTD)进行比较,新算法提高稳定性及精度的同时弥补了传统算法的不足,更具有普适性.     

耦合高效配平策略的旋翼气动特性分析方法

建立了一个耦合高效配平策略和计算流体力学（CFD）技术的旋翼气动特性分析方法.采用CFD求解器进行雅克比矩阵和旋翼流场的计算，建立了一个包含直接配平法的旋翼气动特性分析平台；引入差量法思想对配平策略进行简化，发展了一种差量配平法，解决了直接配平法中雅克比矩阵计算效率低下的问题，显著缩短了计算时间.旋翼流场求解中控制方程采用Navier-Stokes方程，空间离散格式选取2阶迎风Roe格式，时间推进采用隐式lower-upper symmetric Gauss-Seidel (LU-SGS)格式，湍流模型为Baldwin-Lomax (B-L)模型.分别采用建立的基于直接配平方法和差量配平方法的气动特性分析平台对AH-1G和Helishape 7A直升机的旋翼进行了配平计算和研究，对比分析了两种配平方法的计算效率和精度及背景网格的网格量对配平结果的影响.结果表明:与直接配平法相比，耦合差量配平法的气动特性分析方法计算圈数更少，可节省至少60%的计算时间；使用疏化的背景网格虽然不利于捕捉流场细节，但是可以在兼顾配平结果准确性的同时，降低配平计算代价，提高旋翼配平和流场计算的效率.      

非结构弹性体动网格方法应用策略及改进

针对3种代表性的经典弹性体动网格方法，从网格最大变形能力、质量评价指标和空间分布趋势3个方面综合分析了各方法的网格变形特性，重点总结了各方法的适用范围.结合工程实际，研究了特定问题的弹性体动网格应用策略.最后，整合各经典弹性体动网格方法的优势提出了改进方法，并通过NACA 0006翼型的超声速俯仰振荡加以验证.研究结果表明:基于面积的弹性体法适用于加密型网格，而基于距离的弹性体法适用于均匀型网格，该应用策略可为工程问题提供指导.改进方法效果明显，在大变形时相对基于距离的弹性体法的各项网格质量增幅均在10%以上，相对基于面积的弹性体法的最小网格质量增幅几乎达到4倍，并且在数值模拟中具有较高的精度，可显著增强弹性体法的鲁棒性.      

Solid rocket motor propellant grain burnback simulation based on fast minimum distance function calculation and improved marching tetrahedron method

To efficiently compute arbitrary propellant grain evolution of the burning surface with uniform and non-uniform burning rate for solid rocket motor, a unified framework of burning surface regression simulation has been developed based on minimum distance function. In order to speed up the computation of the mini-mum distance between grid nodes of grain and the triangular mesh of burning surface, a fast distance querying method based on the equal size cube voxel structure was employed. An improved marching tetrahedron method based on piecewise linear approximation was carried out on second-order tetrahedral elements, achieved high-efficiency and adequate accuracy of burning surface extraction simultaneously. The cases of star grain, finocyl grain, and non-uniform tube grain were studied to verify the proposed method. The observed result indicates that the grain burnback computation method could realize the accurate simulation on unstructured tetrahedral mesh with a desirable performance on computational time.     

有限差分法的坐标变换诱导误差

有限差分法应用于具有复杂外形的网格时需要进行坐标变换,在此过程中经常会引入坐标变换诱导误差。在柱坐标系下使用均匀网格进行均匀流场计算,计算结果表明,即使物理坐标对计算坐标的变换函数连续可导、计算过程中坐标变换系数直接采用准确的解析式、网格完全正交并且充分光滑,也无法避免坐标变换诱导误差。理论分析表明,产生坐标变换诱导误差的机理是笛卡尔坐标系下的守恒型欧拉方程变换至贴体坐标系下后增加了源项。针对该问题,目前国内外学者通常采用几何守恒律,构建与差分格式相匹配的坐标变换系数计算方法来消除源项。本文介绍了从包含源项的离散等价方程基础上直接进行离散的新算法,在此基础上针对非等距网格条件下MUSCL类格式重构过程进行误差分析,理论推导表明重构中需要考虑非等距插值公式的影响系数,将变量转换至计算空间内进行MUSCL重构才能保证该过程具有均匀网格下的插值精度。通过理论分析及数值实验证明新算法对于均匀流场完全不会引入坐标变换误差。     

航天器在轨全过程表面辐射热计算数值仿真研究

对在轨航天器表面辐射热计算进行了全过程数值仿真研究。航天器结构较复杂,针对不同结构进行区域分解,对几何模型进行相应的规则化,同时采用结构化网格和非结构化网格建立通用的计算网格自生成技术。仿真过程重点考虑了任意曲面的网格自动划分和任意形状交界面的数据传递,兼顾几何结构、物理过程、计算精度和计算速度。将有限元法和能束均匀分布法相结合计算角系数和辐射传递系数。将积分法和能束均匀分布法相结合计算外热流。由于在轨航天器表面多用多层隔热组件包裹,针对这部分结构采用节点网络法和控制容积法计算其表面温度,而未被包裹的结构采用有限元法计算其表面温度。对具有辐射换热关系的非连通区域温度场的有限元计算进行了分析和公式推导。最后,用Microsoft VC++6.0编程设计开发了近地轨道航天器表面辐射热计算仿真软件。     

一种基于求解椭圆型方程的结构动网格生成方法

基于椭圆型网格生成法,实现了一种简单高效的贴体结构动网格生成方法,可用于具有移动边界问题的非定常流动数值模拟。该方法提出,在网格变形过程中,Poisson方程需要的控制网格间距和正交性的源项可以通过提取已知的静态网格源项直接得到,并在整个动网格生成过程中保持不变。因此,在椭圆型网格生成中需要通过外迭代确定源项的过程可以得到省略,而且该方法不需要人工指定参数。这使得方法具有高效和易于嵌入到已有程序中的特点。数值模拟结果证明,采用这种方法获得的网格能够较好地保持静态网格原有的正交性和光滑性,在相同迭代步数约束下,网格求解效率低于传统弹簧模拟法,但鲁棒性优于弹簧模拟法。     

气动热数值模拟中壁面热流后处理方法

针对气动热数值模拟中壁面热流后处理方法进行了系统的数值研究,提出了一种有别于常规微分法和积分法的壁面热流后处理方法.通过3种算例进行了对比分析,结果表明:即使在流场参数相同情况下,采用不同的壁面热流后处理方法计算的热流结果也会有较大差异,且表现出不同的网格依赖性.微分法对网格的依赖性较大,单纯提高差分精度不一定能够减弱网格依赖性和提高计算精度,还与近壁处网格分布密切相关;积分法从能量平衡积分方程出发,极大地减弱了网格依赖性,具有较高的计算精度.所提出方法具有与积分法相当的网格依赖性和计算精度,尤其在驻点区域,其网格依赖性更弱于积分法,且计算简单,提高了计算效率.      

基于误差估计的伴随网格自适应方法

网格自适应是一种有效减少解与网格相关性,改善计算结果的重要途径。给出了一种基于误差估计的伴随网格自适应计算技术,定义了基于误差估计的自适应准则,给出了非结构四面体网格单元细化剖分规则,以及一般的网格自适应流程,介绍了一种基于Rechardson外推的网格收敛技术,运用本文提供的方法对NLR7301翼型进行了数值计算研究,数值试验表明给出的网格自适应具有较高的可有效改善计算结果。     

用动态非结构重叠网格法模拟三维多体相对运动绕流

 发展了动态非结构重叠网格技术。本文提出了基于阵面推进的相邻单元搜索算法和适合于高阶空间离散格式的网格间插值关系定义方法，前者提高了非结构重叠网格在模拟多体间具有大幅相对运动的非定常问题时网格处理的效率，而后者保证了网格间边界附近流场求解精度；本文非定常Euler方程的求解采用有限体积法，通量计算采用Osher格式，时间积分采用双时间步长推进。为验证本文提出的方法，分别对定常M6机翼绕流和机翼下外挂物投放过程的非定常绕流进行了模拟,结果表明本文提出的方法提高了非结构重叠网格在模拟多体间具有大幅相对运动问题时的效率，保证了网格边界处的精度，具有较高的工程应用价值。