基于梯度下降的水滴收集率计算方法

拉格朗日方法在计算三维复杂外形的水滴撞击特性时，一般需要通过插值获得物面网格上的水滴收集率，而这一过程容易造成计算结果失真。建立了一种基于梯度下降的水滴收集率计算方法，通过确定撞击物面网格角点的水滴轨迹实现物面网格中心处的水滴收集率计算，避免插值运算。该方法的计算思路如下：根据撞击点到目标点的距离，采用梯度下降法自适应地调控水滴的初始释放位置，使撞击点落在壁面网格角点附近的给定阈值范围内；采用邻点检索算法，快速确定迎风面内所有撞击网格角点的水滴轨迹；基于流管定义，计算出网格中心处的水滴收集率。采用该方法对典型三维外形的水滴收集率进行了计算，并与相关文献结果进行对比。结果显示计算结果与文献数据吻合良好，有效克服了现有技术中因插值导致的结果失真问题，能够为飞机结冰及防除冰技术研究提供参考。     

结冰面水滴收集率欧拉计算方法研究及应用

发展了一套用欧拉法计算飞机结冰表面水滴收集率的算法。先用CFD方法计算空气场,在获得空气流场分布的基础上,求解水滴相控制方程,进而得到水滴收集率。空气相控制方程和水滴相控制方程均写成典型输运方程的形式,采用一致的有限体积法离散求解,方便了计算程序的编制。计算了典型条件下圆柱和NACA0012翼型的水滴收集率,得到了与实验或文献一致的结果,初步表明本文发展的有限体积算法以及相应的计算代码是有效的。最后给出了本文方法在三维复杂外形水滴收集率计算中的应用,对三段翼、T型尾翼的水滴收集率进行了计算,获得了相应的分布规律。     

大型商用运输机机翼增升构型水滴撞击特性计算

基于欧拉-拉格朗日数值模拟框架,发展了一套适合求解三维多体模型水滴撞击特性的计算方法.采用粒子统计与面积比率计算相结合的方式,实现了水滴收集率的高效精确预测.在常规粒径条件下,该方法计算结果与国外标模试验数据吻合良好.应用上述算法,开展了某大型商用运输机增升构型的水滴撞击特性数值模拟研究,计算结果显示:水滴在缝翼背面缝...     

三维复杂表面水滴撞击特性计算

为预测多段翼和发动机进气道等三维复杂形状表面在结冰气象条件下的水滴撞击特性,提出了一种基于欧拉两相流模型求解水滴收集系数的方法.空气相和水滴相认为是单相耦合的,空气流场由Euler或Navier-Stokes(N-S)方程独立求得.给出了一种水滴相撞击壁面边界条件的处理方法.为避免因局部水滴容积分数异常而导致计算发散,提出了一种自适应的数值扩散模型以增加高阶计算的稳定性.验证计算了圆柱、MS(1)-0317翼型、球面等典型二三维物面的收集系数,并计算了水滴直径呈某种分布时的情形.与公开发表的试验数据对比表明,计算方法准确有效,能很好的应用于三维复杂表面的水滴撞击特性预测.      

一种高效的壁面距离计算方法

采用CFD(Computational Fluid Dynamic)技术进行湍流效应数值模拟时,经常需要计算流场网格点到最近壁面的距离。当网格规模很大,特别是针对一些动网格问题,壁面距离的计算量很大且费时,影响整个流场求解效率。本文对壁面距离计算的直接算法和基于二叉树技术的ADT(Alternating Digital Tree)快速算法进行了对比分析和研究,提出了一种高效的、基于方盒切割技术的快速计算方法,并采用不同外形的CFD计算网格对该方法进行了验证。计算结果表明,新方法的壁面距离计算效率明显高于直接算法和ADT算法,并具有较好的鲁棒性和通用性,可以便捷地移植到现有的CFD计算程序中,从而提高整个流场的数值计算效率。     

考虑遮蔽区影响的旋翼三维水滴撞击特性计算新方法

针对直升机旋翼三维黏性流场特有的复杂环境,建立了一种基于欧拉法的旋翼三维水滴撞击特性计算的新方法。首先,在旋翼桨叶嵌套网格的基础上,发展了一套用于预测旋翼绕流流场的计算流体力学(CFD)模拟方法。然后,为克服传统直升机旋翼二维水滴撞击特性计算方法的不足,充分考虑旋翼流场的三维效应,在嵌套网格中基于欧拉法求解旋翼三维水滴撞击流场。其中,为解决尾流等区域的密度脉冲现象所引起的稳定性和收敛性问题,提出并建立了遮蔽区扩散模型。该模型通过判断遮蔽区变量,在计算域中动态生成遮蔽区域,并随迭代步数逐渐扩散。最后,通过与NACA0012翼型及国外UH-1H桨叶的试验和计算结果的对比,验证了旋翼三维水滴撞击特性计算新方法的可靠性,并进行了温度和水滴当量直径(MVD)对旋翼三维水滴撞击特性的影响分析。结果表明:遮蔽区扩散模型的加入,使二维情况的计算时间减少了22%,并增加了三维情况的计算稳定性,显著提高了旋翼三维水滴撞击特性的计算效率;沿着旋翼桨叶展向位置增大的方向,旋翼桨叶剖面水滴撞击范围有所增大,最大水滴局部收集系数呈先增加后减少再增加的变化趋势,其变化幅度接近50%;旋翼桨叶表面的水滴撞击区域和水滴局部收集系数随水滴当量直径的增加而增加。     

过冷水滴撞击三维机翼的数值模拟

基于欧拉两相流理论,提出一种数值模拟过冷水滴撞击三维机翼的方法.根据水滴拟流体模型建立三维水滴控制方程;提出一种可穿透型壁面边界模拟水滴对机翼表面的撞击;由水滴流场的求解结果得到机翼表面的水滴收集特性.对ONERA M6机翼表面的水滴收集特性进行了研究,计算结果表明所提出的数值模拟方法是有效、可行的,可以为三维积冰的数值模拟研究奠定良好的基础.      

动力对全机水滴收集率的影响计算

针对飞机在飞行中遭遇过冷水滴撞击并结冰现象,建立了适合于发动机带动力情况下结冰过程水滴收集率计算的三维数值方法和计算程序.其基本思路为:采用多块技术与SIMPLE方法计算空气流场,以流场分布的计算结果为基础,求解水滴相的控制方程,进而获得物体表面的水滴收集率.空气相控制方程和水滴相控制方程均写成典型输运方程的形式,采用一致的有限体积法离散求解,方便了计算程序的编制.对某型运输机巡航构型有/无动力条件的水滴收集率进行了比较计算,获得了不同直径水滴在飞机表面的撞击特征以及水滴收集率在飞机机翼、平尾、垂尾和发动机进气道唇口上的分布规律.研究表明:(1)发动机是否带动力对机翼、平尾、垂尾的水滴收集率基本无影响;(2)飞机带动力主要影响发动机进气道唇口处的水滴收集率,带动力后唇口的收集率比无动力情况高,水滴撞击范围增大,在进行防除冰研究和设计时需引起重视.     

霜状冰结冰数值模拟研究

提出了一种三维霜冰结冰的数值模拟方法。基于欧拉-拉格朗日法计算空气-水滴两相流流场,获得了表面水滴撞击特性以及撞击量,求解质量守恒和能量守恒方程获得网格单元内的结冰量,三维霜冰的结冰计算程序采用用户自定义（UDF）编写,FLUENT提供的动网格功能实现了结冰后的网格重构。计算了4°攻角,液态水含量为1 g/m3,水滴直径20μm,结冰时间6 min,温度分别为-19.4°和-28.3°时NACA0012翼型的结冰,并与相关文献结果进行了对比,计算结果与文献实验结果吻合较好,验证了方法的合理性和可靠性。分析了不同的粒径、液态水含量、来流速度等对撞击特性以及结冰的影响,随着粒径、来流速度的增加,水滴收集系数、撞击区和结冰厚度增大,随着液态水含量的增加,水滴收集系数和撞击区是不变的,但是结冰厚度增加。     

非结构网格在积冰数值模拟中的应用

基于非结构网格对翼型表面的积冰进行了数值模拟研究.通过求解气体欧拉方程得到翼型的绕流流场,使用穿透型壁面边界模拟水滴对翼型表面的撞击,采用迭代法对水滴控制方程进行数值求解,得到翼型表面的水滴收集特性;采用不同的积冰冻结模型模拟霜冰和明冰的冻结过程;利用冰层时间推进法模拟积冰过程,基于积冰法向生长假设生成积冰外形.对NA...     

Mesh deformation on 3D complex configurations using multistep radial basis functions interpolation

The Radial Basis Function (RBF) method with data reduction is an effective way to perform mesh deformation. However, for large deformations on meshes of complex aerodynamic configurations, the efficiency of the RBF mesh deformation method still needs to be further improved to fulfill the demand of practical application. To achieve this goal, a multistep RBF method based on a multilevel subspace RBF algorithm is presented to further improve the efficiency of the mesh deformation method in this research. A whole deformation is divided into a series of steps, and the supporting radius is adjusted in accordance with the maximal displacement error. Furthermore, parallel computing is applied to the interpolation to enhance the efficiency. Typical deformation problems of the NASA Common Research Model (CRM) configuration, the DLR-F6 wing-body-nacelle-pylon configuration, and the DLR-F11 high-lift configuration are tested to verify the feasibility of this method. Test results show that the presented multistep RBF mesh deformation method is efficient and robust in dealing with large deformation problems over complex geometries.     

基于CFD动网格技术的飞艇动导数计算方法

由于飞艇和飞机这两种飞行器在飞行原理、气动外形等方面都存在较大的差别,因此现有的飞机动导数工程估算公式和经验数据并不适用于飞艇动导数求解。为了解决上述问题,提出了一种新的基于CFD动网格技术的飞艇动导数计算方法。通过非结构化动网格技术,建立飞艇在流场中的强迫振动模型,采用有限体积法对N—S方程求解,研究飞艇在强迫振动中的气动特性变化规律,并基于数值计算结果求解飞艇的动导数。研究发现,得到的计算结果与风洞试验结果基本一致,表明基于CFD动网格技术进行飞艇的动导数求解是有效可行的。     

Aerodynamic optimization design of general parameters for cycloidal propeller in hover based on surrogate model

A surrogate-model-based aerodynamic optimization design method for cycloidal propeller in hover was proposed, in order to improve its aerodynamic efficiency, and analyze the basic criteria for its aerodynamic optimization design. The reliability and applicability of overset mesh method were verified. An optimization method based on Kriging surrogate model was proposed to optimize the geometric parameters for cycloidal propeller in hover with the use of genetic algorithm. The optimization results showed that the thrust coefficient was increased by 3.56%, the torque coefficient reduced by 12.05%, and the figure of merit (FM) increased by 19.93%. The optimization results verified the feasibility of this design idea. Although the optimization was only carried out at a single rotation speed, the aerodynamic efficiency was also significantly improved over a wide range of rotation speeds. The optimal configuration characteristics for micro and small-sized cycloidal propeller were: solidity of 0.2-0.22, maximum pitch angle of 25°-35°, pitch axis locating at 35%-45% of the blade chord length.      

一类二维非结构网格DSMC方法的实现策略及其应用

研究了一类二维非结构网格DSMC方法的实现策略。在数据结构方面,设计了局部化的数据组织方式,节约了内存与计算时间。发展了一种跟踪模拟分子迁移的算法,该算法仅需少量的逻辑运算与代数运算,不仅可以快速跟踪模拟分子在网格之间的迁移,而且可以准确判别分子与物面是否相互作用,搜索过程中的附带信息给出了分子与物面碰撞的精确时间与位置,避免了重新计算。引入碰撞距离的思想,既减少存储又保证正确的模拟结果。在程序编制过程中,我们充分展现了Fortran90高级语言的主要特性,引入动态数组、指针、链表以及派生类型数据,编制了计算程序。最后对过渡流域高超声绕流进行了数值试验。     

A mesh optimization method using machine learning technique and variational mesh adaptation

Computational mesh is an important ingredient that affects the accuracy and efficiency of CFD numerical simulation. In light of the introduced large amount of computational costs for many adaptive mesh methods, moving mesh methods keep the number of nodes and topology of a mesh unchanged and do not increase CFD computational expense. As the state-of-the-art moving mesh method, the variational mesh adaptation approach has been introduced to CFD calculation. However, quickly estimating the flow fi...     

在MIMD型机上求解板的临界载荷的一种并行算法

提出了一种在ＭＩＭＤ型并行机上求解板的临界载荷的子结构方法的并行算法。此算法实际上就是解决Ａｘ＝λＢｘ广义特征值问题（其中Ａ、Ｂ为正定的箭头型矩阵）。主要通过矩阵变换将此问题转化成适合并行算法的一般特征值问题。给出了此并行算法的并行加速及效率的分析和算例,证明了此并行算法的优越性。     

机翼和翼身融合体跨音速绕流有限元算法的研究

从全速位方程出发研究了机翼和翼身融合体跨音速绕流的有限元算法。给出适用于有限元法和复杂外形的块结构网格生成技术和使单元系数矩阵积分数目减少2/3的方法;并发展了一种内存量少的总体系数矩阵存贮方式。为了提高有限元离散方程组解的精度和加速收敛,提出α—GMRES算法代替传统全速位方程有限元解法中所习用收敛慢的线松弛迭代解法。通过机翼和翼身融合体实例考核,效果是好的,为飞行器先进气动布局的数值计算提供了一种新的有效算法。     

一种物联网的时钟服务器配置方法

为了进一步提高物联网时钟同步的准确性和有效性,文中提出了一种时钟服务器的配置方法.此方法是在建立一个时钟服务器配置优化模型的基础上,利用粒子群优化算法求解出时钟服务器的配置结果,并对此方法进行了仿真实验验证.实验结果表明,此时钟服务器配置方法能够满足物联网时钟同步的要求,并且降低了配置代价.     

翼上螺旋桨构型耦合气动特性及翼型优化设计

针对半环形式翼上螺旋桨构型,研究了螺旋桨-机翼耦合流场特性,并以短距起降（STOL）状态最优升阻特性为目标对机翼翼型进行全局优化。首先,针对螺旋桨-气动面耦合构型,通过动量源法与真实桨叶模型CFD的计算对比,分析动量源法用于该构型设计分析的可行性。其次,为得到有利于桨-翼耦合特征的新翼型,建立了翼上螺旋桨构型自由型面变形（FFD）参数化模型,采用遗传算法对翼上螺旋桨构型机翼翼型进行全局寻优设计,分析了优化翼型参数及流场变化规律。最后,将优化翼型用于三维半环形机翼,分析其流场特性与二维计算结果的异同,验证二维翼型优化的有效性。结果表明：真实桨叶多重参考系（MRF）方法不能准确计算翼上螺旋桨构型下的流场结构,而动量源法计算结果与真实桨叶滑移网格非定常方法较为吻合;采用二维动量源CFD方法进行翼型的遗传算法优化是有效的,受半涵道的保护,二维优化翼型的优势在三维构型中得到了有效继承;翼上螺旋桨构型的翼型优化应当着重关注翼面曲率变化,在本文计算状态下,通过增加桨盘附近翼面曲率、保持附着流动来加强Coanda效应,有效实现了气动增升,优化后机翼升力提高了22.51%,显著减弱桨盘后高压区并产生二次吸力峰值,同时保持了机翼负阻力特性。