一种改进的浸入运动边界算法

浸入运动边界是数值处理流固耦合作用的主要方法之一，其原始权重函数根据经验确定且形式固定，难以准确地定量描述复杂的流固现象。本文提出了一种新的权重函数对格子Boltzmann方法框架下的浸入运动边界进行修正，通过假定零固含率处的权重因子多阶导数为0，引入可调参数b增强权重函数对固含率的关联，实现了对固体边界渗透性能的改进。对静止圆柱绕流、Taylor-Couette流和振动圆柱绕流3个算例开展数值实验，结果表明：随着雷诺数的增加，原始权重函数预测的升阻力系数均高于文献结果，严重低估了边界处流体组分的作用；当新权重函数中的参数b为3时，改进的浸入运动边界效果最佳；新权重函数增强了固体边界的渗透率，提高了流场的预测精度，使平均误差降低了38.5%；最后，运用新权重函数成功地复现了振动圆柱绕流的锁频现象，表明改进的浸入运动边界在处理复杂的流固耦合问题方面具有广阔的应用前景。     

一种包含运动边界的高精度流场数值计算方法

为了准确快速地模拟运动边界的流场,提出一种反馈力源形式的包含运动边界的非定常流场数值计算方法。该方法采用完全正交的网格,以反馈力源作用点的运动来模拟边界的运动。采用物理量及其各阶导数在边界两侧的突跃修正中心差分格式,使之达到二阶精度,以此离散求解二维不可压Navier-Stokes方程。并且提出了与运动边界相适应的反馈力源构造方法及对边界上速度进行插值的方法。基于此数值计算方法,对低雷诺数的圆柱绕流、静止流体中的振荡圆柱以及昆虫振翅运动的二维非定常流场进行了数值计算,计算结果与以往的数值及实验结果非常吻合,表明本文方法与Peskin的浸入式边界方法在处理运动边界问题时具有同样的高效率,且精度高于浸入式边界方法。     

任意运动变形边界流场测量PIV算法的研究：薄膜涡激振动

提出了一种精确测量任意运动变形边界流场的粒子图像速度场PIV算法,能智能识别强流固耦合问题流场测量中任意运动变形的边界信息,然后生成贴体自适应的图像互相关计算窗口,计算获取运动变形边界附近的速度场数据。这一算法可大大提高任意运动边界附近流场的测量精度,为流固耦合问题的理论分析和数值计算验证提供可靠的高质量实验数据。针对数字合成给定流场的粒子图像序列,采用所发展的PIV算法对粒子图像区域中的运动变形边界进行了精确识别,高质量地复现了原始流场信息。最后,对低速闭式循环水洞中的钝体尾流柔性薄膜涡激振动现象进行了PIV实验测量,获得了柔性薄膜大变形运动状态下的瞬态流场特征。     

一种求解低速跨流域流固耦合问题的离散统一气体动理学格式

运动边界及流固耦合问题是低速及高速、连续流域及稀薄流域流动中常面临的多场耦合问题。本文在原始离散统一气体动理学格式基础上，采用任意拉格朗日-欧拉方法及动网格技术，并耦合计算结构动力学，实现适用于低速流动运动边界及流固耦合问题的计算方法。在方法验证方面，进行连续流域圆柱绕流强迫振动和自由振动数值模拟，计算结果与参考文献结果吻合良好，准确捕捉圆柱在锁频区和非锁频区的流动特性。在稀薄流域范围，将方法中的连续流D2Q9离散速度模型替换成跨流域离散速度模型，对比低速圆柱绕流强迫振动算例在两种离散速度模型下的数值结果。虽主要关注低速流动问题，但发展的算法仍具备模拟高超声速运动边界问题的能力，对此进行了Ma=5.0稀薄流振动圆柱研究。低速和高速运动边界模拟能力显示了方法在考虑稀薄气体效应流固耦合问题方面的潜在应用价值。     

曲壁面两相流的简单扩散界面浸没边界法

针对在欧拉-拉格朗日网格框架下曲边界上的两相流润湿边界条件实施问题，本文提出了一种简单扩散界面浸没边界法。在此方法的框架下，静止和移动曲边界的湿润边界条件都可以采用平板边界方式处理。在该方法中，我们在实际边界内外一个拉格朗日网格间距处添加了两层虚拟贴体拉格朗日网格，并通过插值得到三层拉格朗日点上的序参数，从而得到实际的曲边界润湿条件。该方案实施思路简单，不仅适用于静止曲边界，还适用于复杂的运动曲边界。为了验证该方法的有效性，本文在广义守恒相场格子玻尔兹曼方法（LBM）框架下，模拟了静态液滴在圆柱上的铺展、动态液滴撞击圆柱、圆柱在水池中沉降三个问题。结果表明该方法能够有效地模拟曲面上的两相流动，具有良好的精度、鲁棒性和质量守恒性。     

基于深度神经网络的含运动边界非定常流场预测方法研究

针对飞行器设计过程中对流固耦合系统快速预测的需求,探索基于数据驱动的非定常流场建模策略,缩短流场演化求解耗时,从而加快流固耦合系统模拟速度。流固耦合系统中流场演化求解部分等价于含运动边界的非定常流场演化。本文提出了一种基于神经网络的流场预测模型,来学习并预测运动边界附近的非定常流场演化过程。此神经网络可以基于当前流场状态及边界运动信息预测出下一个时刻的流场状态。通过不同振动频率及振幅下的运动圆柱绕流问题,测试了本文提出的神经网络模型的预测精度及泛化能力。神经网络的预测结果和计算流体力学仿真结果中流场结构吻合度较高,且通过对预测的流场数据中边界上的压力积分得到的气动力也具有较高的精度。测试结果展示了该神经网络模型具有良好的预测性能,因此该方法可以用于快速、准确获得运动边界周围非定常流场状态。     

基于DDB的二维/轴对称三波交汇流场数值模拟

采用运动间断边界拟合算法,在运动网格下求解二维与轴对称欧拉方程组,对脱体激波、锥面波、三波交汇等多种稳态波系进行数值模拟,精确得到了多种典型激波形状、波后流场及波系结构。本文提出的运动间断边界拟合算法,系基于经典计算流体力学方法、R-H运动间断条件、运动边界条件和Delaunay背景网格技术,由指定的网格块边界拟合待求激波,并以网格分块方式解决多道波相交问题。通过引入轴对称形式欧拉方程组,直接将此算法推广至轴对称流动领域,并精确拟合了锥面波相交问题。算法精度高,易于在CFD程序基础上实现,可用于超声速飞行器进气道压缩面、唇口设计以及乘波体设计。     

单个中性球形颗粒在三维方腔中的运动

采用多松弛格子Boltzmann方法模拟了单个中性球形颗粒在三维顶盖驱动方腔流中的运动。考虑了展向弱受限的对称边界和强受限的固壁边界两种情况下初始位置、颗粒大小以及雷诺数的影响。对于展向弱受限的情形，发现颗粒的初始位置显著影响着最终的运动轨迹。根据相图被划分为三个区域：外层稳定区、内层稳定区以及涡中心区域。通过对颗粒受力的分解解释了其在极限环上运动的机理。此外，还详细介绍了球形颗粒在极限环上的顺时针旋转运动。随着雷诺数的增加，颗粒逐渐向外围靠近，不断旋转达到相应的极限环轨迹。对于选定的初始位置，观察到在高雷诺数时大颗粒向外迁移，而在低雷诺数时大颗粒的极限环靠近涡中心。对于展向强受限的情形，极限环与颗粒的初始位置无关。随着雷诺数的增加，除方腔左上角外，极限环轨迹有向外迁移的趋势。最后，随着颗粒尺寸的增大，极限环向方腔内部收缩。     

数值模拟中的边界处理方法研究

为适应于模拟超声速流场中的一些特殊物理现象,研究了不同物理现象对边界处理的要求。针对激波、旋涡和声波等流场中典型的流动现象,采用常用的外推、特征和缓冲三类边界处理方法分别进行数值模拟。计算结果表明超声速流场边界采用外推技术,只适于计算带有激波的边界。特征方法还能适于旋涡运动的计算。声波的计算只能采用缓冲层技术。     

欧拉方程计算非定常流动的一种数值方法

提出了利用欧拉方程计算非定常流动的一种数值方法。在利用中心差分的有限体积法和双时间推进法求解非定常欧拉方程过程中，通过在固定物面边界上满足完全的非定常边界条件，计算出非定常流场。算例结果表明，本文方法具有高效、简便的特点，而且能描述激波的非定常运动和变化情况。     

差值法与比值法的误差修正效果比较研究及应用

差值法与比值法对不同测量系统的误差修正效果不一样。本文将两种方法的比较用一个比值来描述,并用推得的公式进行量化表示。在进一步的应用研究中,理论分析与实验数据相互印证,同时得到许多有益的结论,有助于全面掌握测量系统的误差来源,从而灵活选择各种修正方法以获得更加准确的测量结果。     

卫星通信设备抗直升机旋翼遮挡技术测试验证方法研究

抗旋翼遮挡技术是直升机载卫星通信设备的核心技术.根据直升机载卫星通信系统旋翼遮挡的特点,分析研究了直升机无遮挡通信时间比率与直升机航向和卫星波束入射角的关系,以此分析设计了试飞方案.发现了抗旋翼遮挡技术的重大设计缺陷,并通过试飞中的测试数据对抗旋翼遮挡算法进行了优化改进,有效实现了抗旋翼遮挡技术的有效性验证,达到直升机载卫星通信设备定型试飞的目的,为后续该类型设备的试飞奠定了扎实的理论与实践基础,具有较高的工程应用价值.     

基于SWT模型的概率寿命分析体方法

为了研究大应力体积对航空发动机轮盘低周循环疲劳破坏的影响,提高轮盘概率寿命预估精度,针对应变体方法对中、长寿命轮盘预测精度不准的缺陷,提出基于SWT模型的概率寿命分析体方法。该方法以最大主应变平面为临界面的SWT模型为基础,结合考虑尺寸效应的体方法,且不需要进行平均应力修正。通过对某风扇盘进行概率寿命分析,得到的轮盘中值寿命与试验结果吻合良好,预测精度高于应变点方法、应力修正系数法得到的结果。表明基于SWT模型的概率寿命分析体方法可以用来预估航空发动机轮盘概率寿命。     

提高测量RCS的扫频方法精度研究

为提高测量RCS值的扫频方法的精度,从环境要求入手,提出了改进环境的基本措施,给出了一些解决的办法;最后通过测试技术手段消除了隔离误差.从结果可以得出测量的精度达到要求.     

Modelling of the electromagnetic field in the interplanetary space and in the Earth's magnetosphere

A magnetohydrodynamic model of the solar wind flow is constructed using a kinematic approach. It is shown that a phenomenological conductivity of the solar wind plasma plays a key role in the forming of the interplanetary magnetic field (IMF) component normal to the ecliptic plane. This component is mostly important for the magnetospheric dynamics which is controlled by the solar wind electric field. A simple analytical solution for the problem of the solar wind flow past the magnetosphere is presented. In this approach the magnetopause and the Earth's bow shock are approximated by the paraboloids of revolution. Superposition of the effects of the bulk solar wind plasma motion and the magnetic field diffusion results in an incomplete screening of the IMF by the magnetopause. It is shown that the normal to the magnetopause component of the solar wind magnetic field and the tangential component of the electric field penetrated into the magnetosphere are determined by the quarter square of the magnetic Reynolds number. In final, a dynamic model of the magnetospheric magnetic field is constructed. This model can describe the magnetosphere in the course of the severe magnetic storm. The conditions under which the magnetospheric magnetic flux structure is unstable and can drive the magnetospheric substorm are discussed. The model calculations are compared with the observational data for September 24–26, 1998 magnetic storm (Dst _min=−205 nT) and substorm occurred at 02:30 UT on January 10, 1997. This revised version was published online in August 2006 with corrections to the Cover Date.     

数控机床铣螺纹用盘形铣刀的近似造型及误差计算

在数控铣床上加工螺纹，需要设计盘形铣刀，本文求导了铣刀宽度计算公式，分析了产生误差的原因，并提出了改进措施。     

准入审批对基础产业供给波动的影响--一个理论模式的探讨

对基础产业的准入审批进行了经济分析，指出准入审批是主管部门在其目标、激励约束、信息等条件下的必然选择。基于对主管部门危机管理的阐述，提出了一个准入审批对基础产业供给波动影响的模型和两个拓展模型。其结论是：准入审批是导致某些基础产业供给波动的重要原因，并进而创造了自身存在的理由。由此，提出了四点改善基础产业供给的建议。     

航空发动机燃油泵试验台入口油温控制

航空发动机燃油泵试验台入口油温控制系统是一种非线性、时滞系统,难以建立数学模型,常规的PID控制方法不能实现精确的控制。提出采用改进的单神经元自适应PID控制方法实现对入口油温的控制。该控制方法结合了神经网络和PID控制的优点。仿真结果表明,所采用的单神经元自适应PID控制与常规PID控制相比,具有更好的控制品质。     

地火转移轨道误差分析与第1次中途修正计算

针对地火转移过程中出现的各种误差,基于地火转移轨道的误差传递矩阵分析误差发散的性质,在此基础上讨论如何选取轨道中途修正的时机,并基于该矩阵对地火转移轨道第一次中途轨道修正的速度增量进行估算。与微分修正方法的严格计算结果的比较表明,基于该方法定性研究地火转移轨道第1次中途修正速度增量变化和选取合适的轨道机动时机是可行的。使用蒙特卡洛数值模拟对上述方法和微分修正方法进行计算和比较,结果表明,第1次中途修正速度增量大小差异不超过1.2m/s,相对误差不超过6%。在轨道控制精度大约为1m/s的情况下使用该方法代替微分修正方法进行计算,可以节省大量的计算时间。     

基于PSO_LS-SVM的装备研制费用预测研究

在简要介绍了粒子群优化算法基本原理的基础上,建立了基于POS_LS-SVM的装备研制费用预测模型,并应用实例进行了验证。结果表明,采用PSO_LS-SVM方法进行装备研制费用预测,得到的预测结果精度更高、计算速度更快。