合成双射流控制机翼水滴撞击特性

在低速来流条件下,针对前缘位置嵌有合成射流/合成双射流激励器的机翼的水滴撞击特性开展了数值模拟研究,基于Fluent软件,采用Euler气液两相模型和欧拉壁面液膜（Eulerian wall film,EWF）模型,得到的计算结果表明:在合成射流或合成双射流的主动控制下,阻挡了机翼前缘等积冰重点防护区域内的水滴撞击,从而大幅降低了该区域的结冰强度.其机理是:在高频合成射流的作用下,机翼前缘上游附近形成了一对稳定的闭合回流区,形成了水滴的"真空区域".由于回流区内部水滴速度和质量分数较低,改变了机翼前缘水滴运动轨迹和水滴收集率分布,能够减少机翼前缘结冰程度并改变冰形,起到了虚拟气动外形的作用.      

机翼热气防冰及冰脊形成数值模拟

应用N-S方程求解空气外流场和防冰腔内流场，用Euler法获得过冷水滴撞击特性，将内外流场进行耦合传热稳定后，开始结冰，来实现机翼热气防冰及冰脊形成的数值模拟.计算结果表明，热气防冰开启时，加热区蒙皮最低温度为286 K保证加热区没有结冰，但在加热区后的上下表面由于蒙皮温度降低至273.15 K以下，防冰区的溢流水到此处形成冰脊，这个计算结果表明了热气防冰数值模拟的可行性与合理性.但在将前缘积冰除去的情况下，防冰区外形成冰脊对气动特性影响很大，对带冰脊的机翼绕流流场进行数值模拟，对计算结果进行分析，得出环境温度越低、形成冰脊的位置越靠近防冰区，冰脊的高度越高；结冰时间越长，冰脊后的非定常特性越明显，流动细节的捕捉越困难，对气动特性的影响越大.      

发动机整流支板及帽罩积冰的数值模拟

航空发动机在高空结冰气象条件下工作时,其进口部件,如进气道、整流支板、帽罩以及压气机前几级叶片会因为冷空气中液态水滴的撞击和聚集而在表面产生积冰,严重影响发动机的工作性能和飞行安全。采用数值模拟方法,对整流支板和帽罩在相应工作状态下的积冰过程进行了三维数值计算。较小尺寸水滴会增大部件表面的水滴收集区和积冰覆盖区,而大尺寸水滴与固壁面的相互作用带来子水滴的返流和再附着,造成局部表面水滴收集率下降.模拟得到不同因素在相应条件下对部件表面的积冰覆盖范围和厚度分布影响,为发动机防/除冰系统的设计提供理论指导。     

低温光滑壁面上水滴撞击结冰行为

采用高速摄像技术测试低温光滑壁面上水滴撞击结冰过程, 分析了撞击速度、壁面温度和材料热导率对水滴撞击铺展、振荡及结冰行为的影响规律. 结果表明, 低温壁面造成水滴最大铺展直径缩小, 且结冰时间随温度降低而缩短; 当撞击We数提高时, 水滴最大铺展直径增大, 而振荡和结冰时间减小; 同时材料热导率越高, 最大铺展直径越小, 结冰越迅速. 另外, 从热力学角度推导出水滴撞击结冰时间的理论公式, 预测误差<5.3%.     

水滴撞击结冰过程的分子动力学模拟

利用三维分子动力学模拟方法,研究了纳米尺度水滴撞击冷壁面的结冰过程.数值模拟中,统计系统采用微正则系综,势能函数选用TIP4P/ice模型,温度校正使用速度定标法,牛顿运动方程的求解采用文莱特算法,水滴内部结冰过程则通过统计垂直方向水分子温度分布来判定.研究发现,当冷壁面温度降低时,水滴完全结冰的时间减小,但水滴降至壁面温度的时间却增大;同时随着壁面亲水性降低,水滴内部热传递速度减慢(尤其是冷壁面与水滴底端分子层间),水滴内部温度趋于均匀,但水滴完全结冰时间延长.     

过冷水滴碰撞结冰的实验与模拟研究

以宏观结冰/霜过程中过冷水滴的碰撞结冰现象为背景,实验对比了亲水和超疏水表面上常温水滴碰撞、常温水滴碰撞结冰和过冷水滴碰撞结冰的过程,建立了过冷水滴碰撞结冰过程的数值模型,研究了We数和接触角对碰撞结冰的影响。结果表明:相比于常温水滴的碰撞及其碰撞结冰过程,过冷水滴碰撞结冰过程的稳态铺展系数更大;随着过冷度和We数的增大以及接触角的减小,过冷水滴的碰撞结冰与常温水滴的碰撞在水滴形态和铺展系数上的差异逐渐增大。     

三维机翼跨音速颤振的分布式并行计算

以多台微机组成分布式计算系统,利用无限插值方法(TFI)生成三维多块贴体运动网格,以Navier Stokes方程为控制方程,求解三维机翼的跨音速非定常气动力,并与颤振方程耦合迭代计算,求解机翼广义位移响应的时间历程,根据广义位移的时间历程的衰减、等幅和发散振荡等情况确定机翼跨音速颤振临界条件.通过算例验证,计算结果与实验结果和理论分析相吻合.     

螺旋线行波管三维频域非线性注波互作用的计算

基于行波管中慢电磁行波的周期传输特性和能流坡印廷定理,考虑高频结构的衰减、切断、螺距的渐变与跳变对互作用的影响,建立了螺旋线行波管三维场论非线性自洽工作方程组.通过傅里叶展式将时域中的电流变换为与频率相关的交流电流分量,采用等离子体粒子模拟的方法,求解离散化的亥姆霍兹方程获得空间电荷场的三维数值解.计算三维电子轨迹,得到精确的互作用后的电子能谱结构,为多级降压收集极的设计提供关键参数.计算结果与电子所Ku波段的测试值比较具有较好的一致性,并分析了互作用后的电子能谱结构,与多级降压收集极的实验基本符合. 关键词： 行波管 非线性注波互作用 空间电荷场 等离子体粒子模拟     

三维非结构网格Euler方程的LU-SGS算法及其改进

将LU-SGS隐式时间推进格式运用到非结构网格Euler方程的求解中,并对传统LU-SGS格式进行改进,结合网格重排序,发展了一套效率更高的三维Euler方程求解器.以M6机翼及超临界LANN机翼的跨音速无粘流场为算例,将改进LU-SGS格式与四步龙格-库塔显式格式及传统LU-SGS格式进行了比较.计算结果表明:所有格式的计算结果与实验结果都符合很好;传统LU-SGS格式计算效率为显式格式的3倍多,而改进LU-SGS格式计算效率为显式格式的7倍多.     

垂直轴风力机叶片表面结冰的风洞试验

为研究垂直轴风力机叶片表面结冰的规律以及结冰对其性能的影响,对采用NACA0018翼型的风力机叶片进行了风洞结冰试验研究。在风洞试验段内安装了喷水装置,室外的寒冷空气被吸入风洞后与过冷水滴一起吹向叶片并碰撞结冰。测试了不同水滴流量和叶片攻角下的叶表结冰分布及叶片的升阻力系数变化。在一定攻角范围内,叶表结冰量随翼型迎风面积增加而增加;结冰后的阻力系数增大,升力系数减小,叶片的气动特性降低。     

基于大密度差多相流SPH方法的二维液滴碰撞数值模拟

该文结合了Ott提出的修正连续性方程和Adami改进的动量方程, 对空气中的液滴碰撞问题进行了二维数值模拟. 为有效提高计算精度, 推导了适用于大密度差多相流的人工黏性和人工应力方程. 通过表面张力作用下方形液滴自然变化和空气中两液滴互溶的算例, 验证了算法的有效性; 对不同韦伯数 (8.8, 19.8)、不同碰撞参数 (0, 0.5)下的液滴碰撞过程进行了数值模拟, 并与VOF方法对比,取得了较为一致的结果; 进一步计算多个韦伯数、多个碰撞参数下的液滴碰撞, 得到了空气中二维液滴碰撞结果分布图,与实验结果相符合. 结果表明, 该算法对于求解涉及大密度差多相流的液滴碰撞破碎问题十分有效,而且该方法容易拓展到三维, 从而为进一步模拟火箭发动机的二次雾化过程奠定了基础. 关键词： 光滑粒子流体动力学 大密度差 多相流 液滴碰撞     

A new approach for evaluating the current and charge density distributions in electron guns and beams

Current and space-charge density distribution calculation is of great significance for numerical analysis and design of high-current electron guns and beams. When the electrons’ thermal initial velocities are taken into account, though there have been some numerical methods published, the calculation is very complicated. By introducing equivalent meridional potential and projection trajectory theory, the curvilinear axis trajectory equation for electrons neighboring to a central curved trajectory in rotationally symmetric electro-magnetic fields is derived in first- and second-order approximations. The evolution equation of the current density distribution of toroidal electron sub-beams is derived and it can be used to calculate the current and charge density distribution in electron beams and guns in iteration calculation. A compact numerical algorithm for calculating round high-current electron guns and beams was developed and related program was written as well. As examples, the evolution of the current density distribution of a Pierce gun and a periodic magnetic focusing high-current electron beam is simulated. This proves that this method is effective and practical.     

Structure overset grid method and its applications to simulation of multi-body separation

This paper proposes an automatic structure overset grid method, which utilizes the hole-surface optimization with one-step searching, wall-surface grid oversetting, and dynamic overset grid approaches to achieve the high adaptability of overset grids for complex multi-body aircrafts. Specifically, based on the automatic structure overset grids, the method first solves the coupling of Navier-Stokes(N-S) unsteady flow equation and 6DOF motion equation, and establishes the multi-body collision model. Then, the numerical simulation of unsteady flow for complex aircrafts' multi-body separation, the simulation of multi-body separating trajectory and the separation safety analysis are accomplished. Thus, the method can properly handle practical engineering problems including the wing/drop tank separation, aircraft/mount separation, and cluster bomb projection. Experiments show that our numerical results match well with experimental results, which demonstrates the effectiveness of our methods in solving the multi-body separation problem for aircrafts with complex shapes.     

电荷运动轨迹方程解法的讨论

文章讨论了在库仑有心力作用下点电荷二维运动轨迹方程的不同解法：比耐方程法、Runge-Lenz矢量法以及速度积分法.比耐方程法是利用比耐公式求解微分方程,得出点电荷的运动轨迹;Runge-Lenz矢量法和速度积分法都是从点电荷的动力学方程出发,利用矢量积分得出一个常矢量,并应用该矢量分析得出点电荷的轨迹方程.3种方法得到的轨迹方程是一致的.计算表明,库仑有心力作用下的点电荷的二维运动轨迹为圆锥曲线,并分析了不同初始条件下圆锥曲线的类型.     

二甲基醚液滴超临界蒸发的数值模拟

建立了二甲基醚(DME)在N2介质中超临界蒸发的数学模型.采用状态方程法计算了DME-N2体系的高压气液相平衡,在相关文献的基础上得出一套流体热物性和输运参数的计算方法.完成了DME超临界蒸发过程的数值模拟.从数值模拟的结果讨论了超临界条件下环境压力、环境温度和液滴运动等因素对液滴升温、滴径变化以及生存时间的影响.     

一种模拟大密度比多相流的混合算法

在扩散界面法(diffuse interface method,DIM)的基础上提出一种能够处理大密度比(large density ratio)的多相流混合算法.流场信息通过格子波尔兹曼方法(lattice Boltzmann method,LBM)获得;相界面通过直接求解Cahn-Hilliard(C-H)方程确定.为保证在大密度比情况下求解界面方程的稳定性,采用二阶迎风格式来离散方程的对流项.通过对Rayleigh-Taylor(R-T)不稳定、液体中的气泡上升及液滴撞击干燥壁面的数值模拟,验证了方法的可行性.     

Droplet transient migration and dynamic force balance mechanism on vibration-controlled micro-texture surfaces

In this paper, forced vibration was used to regulate the droplet migration, fully recording the transient migration of droplets on a micro-textured substrate under the resonance frequency by a high-speed camera. The influence of resonance frequency and dynamic migration characteristics of droplets on the solid micro-texture surface under lateral vibration were researched. The experiment demonstrates that the driving force is caused by the difference between the left and right contact angles made the droplet oscillate and migrate, and as time t increases, the left and right contact points are periodically shifted and the amplitude of migration increases. Therefore, based on the droplet migration behavior and its force balance mechanism, a spring vibration model of migration behavior of the vibrating droplet micro unit was set up to predict the complete trajectory of its migration on a solid surface. The calculation results show that the theoretical displacement is less than the experimental displacement, and the longer the time, the larger the difference. Affected by the vibration, part of the droplet permeates through the micro-texture, resulting in the droplet losing height and the contact angle becoming smaller as well. While the other part of droplet overcomes the internal surface tension to migrate.     

超短脉冲激光与微小水滴相互作用中电子密度和光场的时空分布

为了研究超短激光脉冲和液滴相互作用过程中电子密度和光场的变化,基于非线性麦克斯韦方程组和电离速率方程,构建了激光等离子体非线性瞬态时域耦合模型,对飞秒激光脉冲击穿微米量级水滴时的电子密度和光场的时空分布进行了计算.结果显示水滴的击穿阈值最小可达2 TW/cm~2,为同等条件下无边界水介质击穿阈值的1/4.随着脉冲能量增强,水滴内自由电子密度峰值区域逆着激光入射方向移动,且入射光越强,水滴对光传播的屏蔽越明显.光束在水滴出射端外部汇聚,汇聚点的光功率密度可达入射光的5倍,且时域波形出现压缩和变形.另外,水滴对激光能量的吸收系数随光强增大而增大,并最终趋于饱和.     

Radiative transfer equation for graded index medium in cylindrical and spherical coordinate systems

In graded index medium, the ray goes along a curved path determined by Fermat principle, and the curved ray-tracing is very difficult and complex. To avoid the complicated and time-consuming computation of curved ray trajectory, the methods not based on ray-tracing technique need to be developed for the solution of radiative transfer in graded index medium. For this purpose, in this paper the streaming operator along a curved ray trajectory in original radiative transfer equation for graded index medium is transformed and expressed in spatial and angular ordinates and the radiative transfer equation for graded index medium in cylindrical and spherical coordinate systems are derived. The conservative and the non-conservative forms of radiative transfer equation for three-dimensional graded index medium are given, which can be used as base equations to develop the numerical simulation methods, such as finite volume method, discrete ordinates method, and finite element method, for radiative transfer in graded index medium in cylindrical and spherical coordinate systems.