运动激波与气泡串相互作用的多介质数值模拟

通过对激波和流体界面相互作用而诱导的大变形界面演化的数值模拟,验证了Level set 方法精确模拟多个流体界面的有效性.采用间断有限元Galerkin方法求解欧拉方程得到流场解,采用5阶WENO格式求解Level set方程追踪多流体界面,界面附近的边界条件由虚拟流体方法处理.对运动激波和两个气泡相互作用过程进行了数值模拟,得到了不同时刻的压力和密度等值线分布,并分析了计算域中两个气泡同是氦气泡,以及一个是氦气泡,一个是R22气泡情况下的计算结果.计算结果表明:利用多界面Level set方程可高质量地捕捉多个流体界面,处理3种多介质流场数值模拟问题.     

NUMERICAL SIMULATION ON MULTI-FLUID INTERFACES

INTRODUCTIONThe study on phenomenon of multi- fluid in-terface driven by shock waves include:( 1 ) totrack the motion of multimaterial interfaces;( 2 )to simulate the chemical reaction on the inter-faces. As the pioneering work,chemical reactionis neglected in this paper. The objective is totrack the interfaces. Of course,the fluids aretreated as compressible medium under such ahigh impact.Although Eulerian schemes workwell for most gas flows,they have been shownto be too dissipative nea…     

基于反馈定理的IB-LB方法(英文)

格子波尔兹曼方法(Lattice Boltzmann method,LBM)和浸入边界方法(Immersed boundary method,IBM)皆为近年来发展的可替代Navier-Stokes(N-S)方程求解复杂流体力学问题的数值模拟方法。本文采用浸入边界-格子波尔兹曼方法(IB-LBM),将IBM作为一种边界处理格式应用到LBM中,模拟了静止单圆柱及并排双圆柱绕流,与文献对比结果吻合良好,证明该方法具有模拟复杂边界物体绕流的能力。与采用虎克定理、直接力法计算体积力项的传统方法不同的是本方法运用反馈定理,即采用流场和流场内物体之间的速度、位移反馈来计算体积力项,使得该方法易于实施并具有优良的并行性。     

GHOST FLUID METHOD FOR RICHTMYER-MESHKOV INSTABILITY SIMULATION

INTRODUCTIONThe energy of the Richtmyer- Meshkov in-stability is excited by the interaction between ashock wave and a contact interface across thatfluids are different.Itis very important thatcap-ture the interface in multi- material flows exact-ly.Here we will apply the ghost fluid method tosimulate two dimensional Richtmyer- Meshkovinstability.The mathematical equation used to simulatethe Ritchmyer- Meshkov instability is the com-pressible Euler equation for two ideal gases.ρρuρvE…     

Stiffened状态方程下γ-model方法跟踪二维多介质流动

在Stiffened状态方程下，运用Level—set方程跟踪界面运动变化，把界面捕捉的等效方程、Level-set函数和欧拉方程组耦合，求解包含比热比和材料参数的耦合形式的流体力学方程组。计算方法采用二阶精度Wave-propagation算法，通过对气体-液体两种流体的一维Riemann问题、激波和气泡相互作用以及气液两种流体Richtmyer—Meshkov界面不稳定性问题进行数值模拟，抑制了不同介质界面两侧的非物理振荡。     

二维多介质流动问题的界面处理方法

就二维可压缩多介质流动问题的数值模拟,给出了一种新的界面处理方法。通过在界面处构造Riemann问题,利用Riemann问题的解分别定义界面两边流体的边界条件,由于Riemann问题的解准确地描述了界面处流体的流动状态,因此得到了更加准确的界面边界条件。本文将由Riemann问题的解得到的界面速度外推到整个流场,重新定义速度场,避免了由于速度的大梯度变化而导致的Level—Set等值线相互交错,因而得到了更加精确的界面位置。利用该方法对水下激波与柱型气泡相互作用问题进行数值模拟,结果表明该方法能准确地捕捉各种物理现象。     

超机动飞行的神经网络动态逆控制

根据反馈线性化理论，讨论了神经网络自适应非线性动态逆控制设计。首先根据时标分离的原则，采用动态逆方法设计了快回路和慢回路控制器；其次提出了模型的神经网络非线性直接自适应控制方案，其中设计一种在线神经网络用于补偿模型逆误差。仿真表明，该控制方案具有较好的自适应能力的鲁棒性。     

基于自适应非结构嵌套网格的旋翼流场模拟

针对直升机旋翼CFD仿真的复杂性,提出了改进的适合于格心格式求解器的非结构嵌套网格算法。采用自适应网格技术在旋翼流场仿真的整个过程中进行网格的自适应加密和疏化操作,以更好地捕捉桨尖涡等流动细节。对于频繁的自适应过程中产生的大量重复点和无用点,采用了高效的交替数字树算法(Alternating digital tree,ADT)和标记-删除-移动算法(Mark,delete,move,MDM)进行删除,节约了不必要的存储。针对格心格式的求解器,采用了基于梯度的网格间插值方式,简化了网格间数值传递的复杂性,同时不降低求解器的精度。对CaradonnaTung旋翼悬停算例和HLISHAPE 7A旋翼悬停算例进行了模拟验证,计算值与实验值吻合,表明本文建立的方法具有良好的鲁棒性和有效性。最后,与未采用自适应时求解器对桨尖涡的捕捉效果进行了对比,结果表明本文所采用的方法可以明显地提高求解器对桨尖涡的捕捉。     

改进的CFD/CSD耦合系统设计方法及其计算

论述了界面上的数据交换方法和时间推进格式在流体动力学（CFD）和结构动力学（CSD）耦合求解非线性气动弹性问题中的关键作用;改进了CFD／CSD耦合系统．包括基于边界元方法设计了CFD和CSD耦合界面的数据转换方法,该方法可在同一映射矩阵处理结构响应和非定常气动载荷转换,保证了耦合边界上的能量守恒;改进了一种松耦合方法流程,在保持模块化基础上提高了计算精度和效率。最后将本文方法应用于二维位移外插算例和A—GARD445．6机翼的动响应分析中,结果表明本文方法具有较高精度。     

一类二阶守恒单调重映算法

在大变形流体力学问题的数值模拟中，经常会涉及到计算网格的重分。基于不同网格的物理量传递便是所谓的重映技术。重映算法是任意拉格朗日一欧拉方法的重要组成部分，本文描述了一类适用于任意网格的二阶守恒单调重映算法。该算法分为网格内物理量的多项式重构、近似积分计算、物理量的单调修正三个部分。本文采用了最小二乘法的思想构造了网格内物理量的梯度，并且通过对物理量的单调修正，保证了算法的精度和单调性。该算法计算公共表面上的网格间质量的改变，所以可以在任意网格上使用。通过一系列数值算例验证了该算法，并说明本文给出的算法是有效可行的。     

Roe格式在多组元燃烧流场数值模拟中的应用

对Roe通量分裂格式进行了详细讨论，并将其推广应用到基于非结构网格的高超声速多组元燃烧流场的数值模拟中。控制方程采用三维多组元Euler方程，重点针对Roe格式开展研究；考虑到燃烧流场求解过程中出现的刚性问题，对化学反应源项采用点隐式方法处理。对超声速条件下钝头体激波诱导燃烧流场进行了数值模拟，比较分析了3种经典化学动力学模型对流场结构的影响，获得与数值模拟、试验数据相符的结果。结果表明发展的方法可以用于多组元燃烧流场的数值分析，且化学动力学模型的选取直接影响激波诱导燃烧流场中诱导区厚度。     

应用确定性优化算法求解多目标气动设计问题的Pareto阵面

将确定性优化算法和Pareto阵面概念结合起来处理了多目标优化设计问题；给出了结合算法及数值过程的细节，并将其应用到了气动优化设计中；描述了如何用确定性优化算法快速抓获多目标优化问题的Pareto阵面以及能够抓获哪些类型的Pareto阵面。数值实验结果表明，确定性优化算法可以准确高效地抓获任意凸的和某些凹的Pareto阵面，故对于此类多目标气动优化问题，可用确定性算法代替进化算法。     

Lattice Boltzmann Flux Solver： An Efficient Approach for Numerical Simulation of Fluid Flows

A lattice Boltzmann flux solver （LBFS） is presented for simulation of fluid flows. Like the conventional computational fluid dynamics （CFD） solvers, the new solver also applies the finite volume method to discretize the governing differential equations, but the numerical flux at the cell interface is not evaluated by the smooth function approximation or Riemann solvers. Instead, it is evaluated from local solution of lattice Boltzmann equation （LBE） at cell interface. Two versions of LBFS are presented in this paper. One is to locally apply one-dimensional compressible lattice Boltzmann （LB） model along the normal direction to the cell interface for simulation of compressible inviscid flows with shock waves. The other is to locally apply multi-dimensional LB model at cell interface for simulation of incompressible viscous and inviscid flows. The present solver removes the drawbacks of conventional lattice Boltzmann method （LBM） such as limitation to uniform mesh, tie-up of mesh spacing and time interval, limitation to viscous flows. Numerical examples show that the present solver can be well applied to simulate fluid flows with non-uniform mesh and curved boundary.     

圆柱-球体三自由度超声电机及其控制

研制了一种一阶纵振和二阶弯振模态组合的圆柱-球体三自由度超声电机。研究了其中两个自由度的运动轨迹跟踪控制策略,并将其应用于电机样机的控制。其中,运转速度的控制采用激励电压幅值调节法,针对位置控制提出了两种运动轨迹控制策略一种是逐点比较控制,另一种是矢量分解控制。本文提出了调节激励电压幅值的3种脉冲宽度调制(PWM)方法,对这些方法作了比较,并通过一系列实验得到验证。结果表明当改变激励电压幅值时,保持3个振动模态的相位差不变和激励电压信号不失真,是简化控制过程、提高控制性能的关键;矢量控制法效果较好,并有轨迹跟踪误差小、运动平稳和噪声小等优点。     

现代空战单机多目标攻击火控计算与决策实现

研究了现代空战中近距格斗与超视距攻击的策略性；介绍了多目标攻击火控系统原理，并用快速模拟法对导弹攻击区进行了计算仿真；阐述了单机多目标攻击的战术决策思想，以及层次分析法单机多目标攻击战术决策方法，并用仿真实现了典型单机多目标攻击空战的攻击决策。仿真证明，本文论述的战术决策思想和决策方法是正确的，有效的。它对现代歼击和火力控制系统具有重要的参考价值。     

基于VXI总线的飞机供电系统参数测试处理系统的实现

介绍了一种基于VXIbus的飞机供电参数的自动测试处理系统的设计和实现，系统设计中采用符合VXI部总线规范的仪器模块，实现多通道高精度及快速的数据采集，采集方式包括并行和巡回检测两种，采用即插即用结构功能扩展非常方便。在HPVEE5．0软件平台下，利用软前面板技术可用交互方式控制VXI仪器模块，用户界面直观明了，大大简化了VXI仪器的操作使用，根据相关标准采用VC^ 和HPVEE相结合的方法设计了数据处理软件，建立了动态链接库，通过单相交流电源系统的测试试验，表明系统性能先进，稳定可靠，测试精度高，具有广泛的应用前景。     

基于基因算法与博弈论的气动高升力优化(英文)

讨论了基于基因算法与博弈论的组合优化算法在高升力气动优化问题中的应用。引入基因算法与博弈论相结合的分布式组合算法 ,可以将复杂的优化问题分解为几个简单的局部优化问题。文中论述了组合优化算法的构造方法 ,并应用于高升力多段翼型气动优化。与传统基因算法的数值计算结果进行了比较 ,表明本文构造的方法具有高效收敛性及强的鲁棒性 ,可广泛应用于先进气动设计问题。     

共轴刚性旋翼气动干扰数值计算方法

建立了一个适用于共轴刚性旋翼气动特性分析的数值模拟方法。该方法采用任意拉格朗日欧拉方法(Arbitrary Lagrange Euler,ALE)描述的可压缩Navier-Stokes(N-S)方程求解流场,采用低数值耗散的Roe格式进行空间离散;使用多重嵌套网格方法以模拟双旋翼的运动。针对共轴刚性旋翼配平,引入"差量修正"策略解决了传统配平中雅克比矩阵计算复杂的问题。首先,对Harrington-2共轴双旋翼的悬停气动性能进行了计算,然后,对某2 m直径共轴双旋翼的悬停及前飞状态进行了计算,并与试验值进行了对比。结果表明:在典型状态下拉力系数的计算结果与试验值误差在3%以内,扭矩系数的计算结果与试验值误差基本在5%以内;所采用的数值计算方法对旋翼涡尾迹特征具有较高的捕捉精度,可以有效模拟共轴刚性旋翼悬停和小速度前飞下的复杂流场及其细节特征。     

毛细流动聚焦的实验方法及过程控制

作为一种基于毛细流动制备微纳尺度液滴的技术，毛细流动聚焦（Capillary flow focusing）在工程领域具有重要应用。回顾了基于吹气式和吸气式核心装置的毛细流动聚焦实验方法，展示了完整的测试平台，介绍了施加外部激励控制流动聚焦射流破碎的实验方法。同时，给出了同轴界面的拍摄方法，可获得清晰的内外层流体界面图像。在对流体锥形收缩阶段的研究中，探讨了几何参数与流动控制参数对流体锥形的形态与稳定性的影响。在稳定锥形下，研究了流动控制参数对液体射流直径、扰动波长及复合射流界面耦合的影响，并基于光的折射定律，对复合射流外层界面透镜效应所导致的内界面失真进行了修正。在对激励作用下射流破碎的研究中，考察了射流长度随振幅的变化，建立了尺度率关系，探讨了频率对生成液滴的单分散性及粒径的影响规律，为在实际应用中可控制备单分散性微液滴提供了理论与技术支持。     

一种充分局部化的子结构并行分析方法

机群系统是一种低成本的松散耦合型分布式并行平台,它要求并行算法设计时遵循"分而治之"的原则,尽量降低节点任务之间的相关性。本文在FETI方法和A-FETI方法的基础上,直接从力学概念出发,提出了一种充分局部化的FETI方法。该方法在进行子域界面处理时引入三重变量:界面节点位移、界面节点力、分区框架上耦合节点位移,由此得到一组近似解耦的界面方程,使得各个子域的计算相对于经典的FETI算法更加独立。对得到的界面方程采用预处理共轭投影梯度法(PCPG)并行求解,所采用的预处理算子为局部化的集中型Dirichlet算子。分别在自建的两套4节点PC机群上进行了两组算例的验算,结果表明,本文方法具有很好的计算精度和收敛速度。并行加速比达到3.76。