非易燃易爆溶剂在喷涂烧蚀隔热涂层中的应用

为消除易燃易爆的风险,使用非易燃易爆溶剂SK-410B代替汽油作为烧蚀隔热涂料的稀释剂。测试了涂层的剪切强度、低温性能、耐湿热性能。电弧风洞试验结果表明:涂层抗冲刷性能良好;试片的最高背面温度为55℃,可满足使用要求。     

存在烧蚀的高超声速压缩拐角流动数值研究

耦合求解热化学非平衡流控制方程和烧蚀壁面边界条件,进行存在石墨烧蚀的压缩拐角流场数值模拟。流场化学反应采用16组元(N2,O 2,NO,N,O,NO +,N+2,O +,N+,CO,CO 2,C,C2,C3,CN,e-)29个反应的非平衡模型,热力非平衡的双温度模型下,不同反应采用不同控制温度。石墨材料表面反应包括碳的氧化反应、碳催化的O 原子复合反应和碳的升华反应。对15°、18°、24°压缩拐角模型,在自由流 Ma =10~30,总焓值6~55 MJ/kg 范围,分别进行无烧蚀的壁面催化与非催化条件和石墨烧蚀条件下的流场计算,分析各类条件下的流场结构、流动分离特性以及流场热化学参数分布特点,研究壁面条件对流动特性的影响。结果表明：流动分离可能性和分离区范围随着压缩拐角斜面倾角增大而增大,随来流马赫数增大而减小;相对于低壁温条件,无烧蚀的辐射平衡壁温和壁面烧蚀条件下流动分离区增大,斜面上压力、摩阻和热流峰值点也有所后移。     

含碱金属杂质的高超声速烧蚀流场特性数值模拟

针对高超声速飞行器的烧蚀现象,采用Newton迭代的空间推进方法,数值求解化学非平衡的抛物化Navier-Stokes(PNS)方程.以黏胶基碳布-氨酚醛树脂为烧蚀材料,建立了含碱金属杂质的18组分33反应的化学反应模型,计算了小钝锥体的壁面温度,并和现有文献结果进行对比验证.在此基础上,对含碱金属杂质的高超声速烧蚀流场进行了模拟,深入分析了含碱金属杂质的烧蚀反应对流场温度、压力和电子数密度的影响.结果表明:烧蚀反应使得流场温度较非烧蚀情况增加了10%~15%;考虑碱金属杂质Na后,流场电子数密度增大了1~2个数量级.      

线性GMRES算法在求解全速位方程中的应用

从全速位方程出发,利用压强极小积分得到有限元方程组。采用适用于有限元法的块结构网格生成技术,给出求解有限元总体方程组的线性GMRES算法。该算法比线松弛迭代法的收敛速度快得多。对M6机翼跨音速有升力情况进行了计算,其结果与实验值符合较好。     

非线性区间迭代法在结构后屈曲分析中的应用

使用区间迭代法跟踪结构后屈曲平衡路径。区间迭代法可以在接近临界载荷处自动解出屈曲前和屈曲后的两个可能平衡位形,从而避免了由于刚度矩阵奇异而使迭代停止的数值困难。一些典型例题的数值解,证明了非线性区间迭代法可以穿过极值点,并能追踪到极值点之后的解曲线,而且对所谓“位移回弹”（Ｓｎａｐ－ｂａｃｋ）问题,区间迭代法仍然是有效的。     

高超声速化学非平衡欧拉方程数值模拟

本文给出了求解高超声速化学非平衡无粘流场的数值模拟方法。计算中采用守恒变量型的ＮＮＤ差分格式和全流场捕捉技术;建立了（数学）奇性轴上的控制方程,克服了通常在其上插值方法所引起的扰动,改善了收敛性和计算精度;化学反应模型由五组元五反应组成;温度和压力的计算是非迭代的,有别于国内外惯用的迭代算法;此外,还发展了背风过低压区的处理方法,从而提高了欧拉方程的应用范围,使无粘非平衡绕流的计算攻角达到３０°     

共轭梯度方法在求解不可压N-S方程中的应用

采用非结构网格上的SIMPLE(Semi-implicit Method for Pressure-linked Equations)算法、k-ε湍流模型求解了三维不可压Navier-stokes方程,对低速不可压粘性流场进行数值模拟,在求解压强方程时采用共轭梯度方法,求解其它变量方程时,采用预处理的BICG(biconjugate gradients)算法.对NACA0012翼型绕流流场和飞艇绕流流场进行数值模拟并对结果进行分析,取得较好的结果.     

隐式Runge-Kutta方法在非定常化学非平衡流动中的应用

推导了求解NS方程组的隐式龙格-库塔方法,并从稳定性、精度和效率三方面展开讨论。通过定常和非定常算例对该方法进行验算,分析表明该方法计算结果正确、稳定性好、精度高、效率高,适合求解大型复杂非定常化学非平衡流动。     

预处理方法在化学非平衡流场数值模拟中的应用

将预处理方法发展应用于化学非平衡流场的数值模拟中。采用有限差分,LU-SGS隐式方法对二维可压化学非平衡控制方程进行了耦合求解,对化学源项和无粘项隐式处理,粘性项则采用显式。通过化学反应剪切层和激波点火的算例表明,所采用的数值方法能够有效求解化学非平衡流场,证明了本算法的可行性,为下一步的工程应用奠定了基础。     

相位延迟边界条件在叶轮机械颤振分析中的应用

基于带相位延迟的周期边界条件,建立了某跨声速转子的双通道高效气动阻尼计算模型.数值计算了该转子的气动性能、颤振边界和叶片模态,和实验数据吻合较好.通过传统的多通道能量法以及双通道方法计算了叶片在一弯模态,不同叶片间相位角条件下的气动阻尼,获得了基本一致的计算结果,而双通道方法相比于传统的多通道能量法计算效率提升约7.7倍,内存需求约为后者的0.45倍.不同叶片振幅对气动阻尼结果的影响研究表明,对于较小的叶片振幅,流动非线性对气动阻尼计算结果仍然有显著的影响.不同工况的计算结果表明:叶片间相位角对转子叶片的气动阻尼有显著的影响,对于该转子最小的气动阻尼均在叶片间相位角为-42.4°时得到;同时,在近颤振状态,不同叶片间相位角对应的气动阻尼均小于近设计状态.     

基于欧拉方程的机翼设计方法及应用

以Takanashi提出的三维机翼设计理论为基础，研究与发展了一个基于欧方程和“正－反迭代、余量修正原理”的机翼设计方法。用改进的无限插值方法生成绕机翼的O－O贴体网格，采用三维欧拉方程作为充动分析计算的基本方程，该设计方法已用于某无人机机翼和一个超临界机翼设计，设计结果达到了预期的目标。     

基于动态边界控制的气动优化设计方法研究

本文根据控制变形物体边界变化的思想,建立了一种基于动态边界控制的气动优化设计方法,此方法通过在求解非定常流动方程过程中计算设计目标函数随设计变量变化的敏感性导数,这样在优化过程中不再需要求解定常流动方程,或共轭方程,从而使得优化过程得到简化并且便于在不同条件下应用.文中以二维翼型为例,用非定常欧拉方程为基本流动控制方程,针对两种优化设计问题,开展了翼型的优化设计,取得了良好的优化设计结果.     

采用附面层抽吸(BLS)控制流动分离的数值模拟

采用数值模拟的方法,研究了附面层抽吸(BLS)对两种具有不同分离特征的叶栅流场的气动影响。通过比较抽吸前后极限流线谱和叶栅损失的变化,对抽吸前后的流场进行定性分析。结果表明,采用附面层抽吸技术可以明显改善大转角、大分离叶栅的气动性能。其中,对以闭式分离为主要特征的叶栅流场,其最佳抽吸位置为主分离区的起始位置,即鞍点附近;对以开式分离为主要特征的叶栅流场,其最佳抽吸位置位于主分离区的上游。     

基于高阶物面近似的自适应间断有限元法欧拉方程数值模拟

将高阶间断有限元与网格自适应相结合,于非结构网格上求解二维Euler方程。将数值解多项式的高阶项贡献用人工粘性项系数的形式进行量化,网格自适应过程中以人工粘性项系数作为网格自适应的指示器。在系数达到设定的上限阀值的区域进行网格加密,在系数达到设定的下限阀值的区域将迭代过程中加密过的网格稀疏以减少网格量。所有自适应均在高阶曲线逼近真实物面的基础上进行,以保证数值结果的精度。典型数值算例结果与实验结果进行了对比,表明采用该自适应间断有限元法可以保证以尽可能小的计算量得到高精度结果。     

扇翼飞行器翼型附面层控制数值模拟

基于扇翼飞行器翼型特殊的几何形状及流场特性,在原有翼型的弧形槽下方和后缘加装控制阀门,通过调节阀门开启及开启尺寸的大小,利用弧形槽低压涡所产生的吸力对翼型后缘的附面层进行一定的控制,达到增升减阻的效果。通过采用计算流体力学的方法对其机理及阀门开启尺寸的影响进行了详细计算和分析,研究表明当阀门开启的尺寸为10 mm时,修改翼型的最大升力系数、失速迎角及相同迎角下的升力系数和推力系数均大于基本翼型;随着阀门开启尺寸的增大,修改翼型的最大升力系数和失速迎角均减小,但是在失速前,修改翼型在相同迎角下的升力系数大于基本翼型。此方法可以改变先前通过增大横流风扇的转速来提高其气动性能的做法,减小了能量的消耗,增大了整个飞行器的航程,为扇翼飞行器能够早日投入实际运用奠定了一定的理论基础。     

基于控制理论和NS方程的气动设计方法研究

研究了基于控制理论和NS方程的气动设计方法,针对给定的目标函数表达形式,应用该设计理论在计算坐标下详细推导了相应的共轭方程及边界条件具体表达形式,以及梯度方程求解表达式,通过合理的数学变换,导出了共轭方程在笛卡尔坐标系下的直观表达形式,发展了有效的共轭方程数值求解方法,通过流动控制方程数值求解、共轭方程数值求解、目标函数对设计变量的梯度求解和优化算法等方面的有效结合,研究与发展了一种新的气动设计方法,以二维机翼气动设计为例,成功进行了亚、跨音速情形下的相关设计算例研究,研究结果表明应用控制理论和NS方程的气动设计方法在设计理论、适用性以及时间花费等方面都有着很好的特色和优点,且设计结果也更为可靠.     

Derivation of the Prandtl-Sirazetdinov equations for optimally controlled boundary layer based on the Navier-Stokes equations

The Euler-Lagrange equations are derived based on the Navier-Stokes equations. The limit transition is performed at high Reynolds numbers.     

On the application of the boundary element method in coronal magnetic field reconstruction

Solar magnetic field is believed to play a central role in solar activities and flares, filament eruptions as well as CMEs are due to the magnetic field re-organization and the interaction between the plasma and the field. At present the reliable magnetic field measurements are still confined to a few lower levels like in photosphere and chromosphere. Although IR technique may be applied to observe the coronal field but the technique is not well-established yet. Radio techniques may be applied to diagnose the coronal field but assumptions on radiation mechanisms and propagations are needed. Therefore extrapolation from photospheric data upwards is still the primary method to reconstruction coronal field. Potential field has minimum energy content and a force-free field can provide the required excess energy for energy release like flares, etc. Linear models have undesirable properties and it is expected to consider non-constant-alpha force-free field model. As the recent result indicates that the plasma beta is sandwich-ed distributed above the solar surface (Gary, 2001), care must be taken in modeling the coronal field correctly. As the reconstruction of solar coronal magnetic fields is an open boundary problem, it is desired to apply some technique that can incorporate this property. The boundary element method is a well-established numerical techniques that has been applied to many fields including open-space problems. It has also been applied to solar magnetic field problems for potential, linear force-free field and non-constant-alpha force-free field problems. It may also be extended to consider the non-force-free field problem. Here we introduce the procedure of the boundary element method and show its applications in reconstruction of solar magnetic field problems. This revised version was published online in August 2006 with corrections to the Cover Date.     

一种抑制激波-边界层相互作用的新型无源被动控制

为了探索一种抑制激波-边界层相互作用的新型无源、被动控制方法,采用有限体积法计算流体力学程序,分别对压缩拐角和激波入射平板两种典型流动的附面层分离进行了无源被动控制的数值研究.数值结果表明,新型无源被动控制方法消除了20°压缩拐角流动所产生的分离泡;对激波入射平板所诱导的附面层分离,采用的自适应无源控制方法可将分离区长度减小为无控制时的58％、总压恢复可比无控制时提高7.2％、x=0.85m截面最大回流速度较无控制时减小了69.3％.从而验证了无源被动抑制控制激波-边界层相互作用的可行性.     

矩阵函数在解线性方程组中的应用

本文利用矩阵函数及其微积分运算来求解一阶线性常系数非齐次微分方程组。