可压缩流动的高分辨率低耗散混合格式研究(英文)

针对可压缩流场的计算,提出了一种混合型激波捕捉格式来获得对于光滑区域和流场间断都具有高精度的解。该混合方法采用有限体积形式的五阶WENO格式来计算流场的光滑区域,而采用MUSCL格式来捕捉流场中的间断。该混合格式结合WENO格式的低耗散特性和MUSCL格式的高分辨率特性。两种格式之间通过一个间断感受因子进行切换。文中选择了三种典型的间断感受因子进行对比研究。WENO格式和MUSCL格式均为激波捕捉格式,因此混合开关中的常数的选择范围相对较宽。选取了若干典型算例对于混合格式的性能进行了研究,并对不同格式的计算精度和效率进行了对比。数值结果表明,与WENO和MUSCL相比,本文提出的混合方法在较广的范围内均表现出良好的特性。     

悬停旋翼无粘流场数值模拟中的多重网格方法

发展了一种基于不同空间离散格式的多重网格算法,并应用于悬停旋翼无粘绕流的Euler方程数值模拟。由于悬停旋翼流场中存在不可压区域,同时旋翼尾涡系统的发展需要较长的时间,使得旋翼流场的收敛速度远低于固定翼流场,因此研究旋翼流场的多重网格算法具有重要意义。空间离散采用了Roe s FDS格式和Jameson中心有限体积格式,时间推进应用了五步Runge-Kutta方法。采用多重网格V循环方式,对一跨声速悬停旋翼无粘流场进行了数值计算。计算结果表明:多重网格算法可以显著加速悬停旋翼无粘流场的数值计算收敛速度;无论在激波分辨率还是在计算精度上,Roe s FDS格式都优于JST格式。     

三角翼大攻角绕流场数值模拟研究

本文基于N-S方程采用Jameson的显式中心有限体积格式和Roc的上风格式对三角翼大攻角绕流场进行初步数值模拟，空间网格为H-0型网格。计算结果表明两种计算格式均可较好的模拟大攻角绕流场。本文计算结果为今后开展全机大攻角绕流场数值模拟提供了计算经验。     

塞式喷管流场计算

采用NND差分格式，对塞式喷管轴对称流场进行了数值模拟。通过流场分析，可以对塞式喷管火箭发动机有清楚的认识。计算表明，NND格式有很好的分辨率；无粘和有粘的情况流场结构和性能相差悬珠；塞式喷管的自适应能力很好。     

ENO-AUSMPW格式对激波涡干扰的数值模拟

高阶ENO格式在复杂流场计算中具有重要的地位,它理论上可以达到任何阶的精度。本文采用二维高阶ENO插值方法,构造了高阶AUSMPW格式,并对包含复杂激波结构的激波涡干扰流场进行了数值模拟。计算结果表明了高阶ENO插值下的AUSMPW格式具有较高的激波分辨率和较低的数值耗散。     

高精度ENO格式在射流数值模拟中的应用

本文分析了ENO格式的特点，并应用于全NS方程的迁移项和压力顶，模拟了轴对称和三维射流。首先对轴对称自由射流为算例比较了高阶TVD格式和ENO格式，证实了ENO格式捕获复杂激波波系的能力；其次计算了四喷管湍射流干扰流场，给出了相应的计算结果；最后计算了四股射流流入圆管内的流场，给出了物理量沿壁面和截面的分布。     

跳点有限体积法及其在跨音速流场计算中的应用

[1]和[2]曾用几种差分格式求解一维N-S方程和欧拉方程,计算表明,Hopscotch格式(跳点格式)计算速度最快。为探讨这种格式在多维流场特别是跨音速流场计算中应用的效果,本文将此差分格式应用于有限体上,成为跳点有限体法,并用该方法求解二维欧拉方程组。对若干个喷管和叶栅所做的计算得到了满意的结果,计算时间较目前常用的差分格式少得多。     

介绍一种适合于计算跨音速流动的混合差分格式

本文介绍MacCormack格式和迎风格式的一种混合差分格式。它适用于跨音速流计算。本文用它计算了绕圆弧翼型的跨音速流场,并详述了处理方法和结果。     

用隐式WENO格式计算悬停旋翼跨声速流场

发展了一种基于高阶迎风格式的计算悬停旋翼跨声速流场的隐式有限体积法。对流项采用Roe通量差分分裂格式,使用五阶WENO格式进行左右状态重构,并与MUSCL插值进行比较;粘性项采用中心有限体积法。为提高收敛到定常解的效率,时间推进采用LU-SGS隐式方法。数值模拟采用了一种能够有效传递网格间流场信息的重叠网格,其中使用了三层内边界和贡献边界的方法以便插值的直接进行。用该方法对一跨声速悬停旋翼粘性流场进行了数值计算,数值结果表明:与MUSCL格式相比,WENO格式对激波位置捕捉得更准确,具有更强的涡捕捉能力,同时还表明了WENO格式在很大程度上能够克服涡耗散问题。     

MLP高阶重构格式应用

在有限体积框架下,利用MLP(multi-dimensional limiting process)系列重构格式结合HLLHLLC(Harten-Lax-van Leer with contact)近似黎曼求解器,同时引入激波探测函数进一步降低MLP在光滑流动区域的数值耗散,数值模拟了超声速前台阶流动、结冰翼型的非定常流动、高超声速双楔流动和DLR F6-WB跨声速流场,研究了MLP系列格式在可压缩复杂流场中的表现.结果表明:在多维空间中,MLP格式能够在如强斜激波与网格线不重合等复杂流场数值模拟中保持严格的流场单调性;具有和传统MUSCL(monotone upstream-centered schemes for conservation laws)格式类似的计算效率,可以实现5阶,甚至更高阶重构;数值耗散更低,捕获更准确的激波位置,对航空工程数值模拟具有重要意义.     

含不等长双裂纹无限大板应力强度因子

采用交替法求解了含任意相对位置不等长双裂纹的无限大极受均布拉应力作用的应力强度因子。此法以平面弹性力学的一个问题为基础,即含单裂纹无限大板,裂纹面受集中载荷,此问题的解已知,然后交替地满足边界条件。各种计算结果已绘成曲线。     

含裂纹半无限大板受集中力作用的应力强度因子

本文釆用交替法求解了含边缘裂纹或内部裂纹或半无限长裂纹的半无限大板在板面内受集中力作用的应力强度因子。此法以平面弹性力学的两个问题为基础,而此两问题的解是已知的。第一个是裂纹表面受任意载荷的无限大板,第二个是半平面在边界上受任意载荷作用。然后交替地满足边界条件。计算结果已绘成曲线。     

无纬碳布增强针刺毡C/C复合材料性能的研究

比较了几种不同工艺制备的针刺毡C/C复合材料。对针刺碳毡织物首先进行预增密处理,得到初始密度和碳纤维含量较高的坯料,然后用树脂浸渍法进一步致密化。研究表明,用该方法制备的C/C复合材料比未经预处理的试样,拉伸强度提高39%,压缩强度提高14%,层间剪切强度提高36%。通过SEM观察和常温力学性能的测试,分析表明工艺的改进是强度提高的主要原因。     

多重网格法求二元跨声速流在一种半自适应网格上的欧拉解

本文采用MacCormack两步显式格式,用有限体积法求解了二元跨声速欧拉流。推导了物面边界条件,采用了特征远场边界条件及远场环量修正。利用保角变换方法生成O型贴体网格,并修改得到了一种在激波处局部加密的半自适应贴体网格。采用多重网格及焓阻尼加速收敛技术计算了NACA0012翼型的跨声速气动特性,得到了十分满意的收敛过程和计算结果。     

Motion equations of hemispherical resonator and analysis of frequency split caused by slight mass non-uniformity

The mass non-uniformity of hemispherical resonator is one of reasons for frequency split, and frequency split can cause gyroscope to drift. Therefore, it is of great significance to analyze the relationship between mass non-uniformity and frequency split, which can provide a theoretical basis for mass balance of imperfect resonator. The starting point of error mechanism analysis for gyroscope is the motion equations of resonator. Firstly, based on the Kirchhoff-Love hypothesis in the elastic thin shell theory, the geometric deformation equations of resonator are deduced. Secondly, the deformation energy equation of resonator is derived according to the vibration mode and relationship between the stress and strain of hemispherical thin shell. Thirdly, the kinetic energy equation of resonator is deduced by the Coriolis theorem. Finally, the motion equations of resonator are established by the Lagrange mechanics principle. The theoretical values of precession factor and natural frequency are calculated by the motion equations, which are substantially consistent with the ones by the finite element method and practical measurement, the errors are within a reasonable range. Simultaneously, the varying trend of natural frequency with respect to the geometrical and physical parameters of resonator by the motion equations is consistent with that by the finite element analysis. The above conclusions prove the correctness and rationality of motion equations. Similarly, the motion equations of resonator with mass non-uniformity are established by the same modeling method in case of ignoring the input angular rate and damping, and the state equations with respect to the velocity and displacement of vibration system are derived, then two natural frequencies are solved by the characteristic equation. It is concluded that one of reasons for frequency split is the 4th harmonic of mass non-uniformity, and thus much attention should be paid to minimizing the 4th harmonic of mass non-uniformity in the course of mass balancing for imperfect resonator.     

二维机翼二自由度跨音速颤振分析

 <正> On the assumption that airfoils take small amplitude ha-rmonic oscillating and middle reduced frequency, the two-dimensional, inviscid, small-disturbance nonlinear unsteady transonic flow equation can be split into two parts. One is the nonlinear steady small-disturbance equation, the other is the time-linearized unsteady small-disturbance equation. The steady equation is solved by Carlson code. The unsteady equation is solved by integral method. The standard v-g method is used to solve the flutter eigenvalue equations. A transonic flutter analysis is performed for a NACA 64A006 airfoil with pitching and plunging degrees of freedom.The Mach numbers considered are 0.7, 0.8, and 0.85. The aerodynamic coefficients are obtained by the above method. For each Mach number, the flutter speed and the corresponding flutter reduced frequencies are achieved by varying the airfoil-airmass ratio, plunge-pitch frequency ratio. The results are compared with other computational method and good agreement has been observe     

高速飞行器翼面结构热振动试验的TARMA模型方法

高速飞行器翼面结构的热振动试验研究对这类飞行器的设计和安全飞行具有重要的意义。采用时变自回归滑动平均(TARMA)模型方法建立了受热时变结构系统模态频率辨识的数学模型,并用一个数值算例进行了验证。将地面振动测试系统与瞬态热环境模拟系统相结合,设计了翼面结构热振动试验系统并模拟结构的瞬态温度场,同时对纯随机激振力激励下受热时变结构系统的振动位移信号进行测量,并用TARMA模型对时变固有频率进行了辨识,获得了前4阶固有频率随加热时间的变化规律,并将辨识结果与数值计算结果进行了比较,两者误差在5%以内。另外,在稳态均匀热环境下辨识得到的结构系统固有频率变化与数值计算结果也吻合得很好。通过将均匀温度场与瞬态温度场下的结果进行对比分析,指出了瞬态热环境下时变结构的固有频率随加热时间变化的趋势主要由结构材料属性的退化和结构内部不均匀热应力的影响共同决定。     

复合材料波纹板剪切载荷作用下的屈曲试验与分析

复合材料波纹板稳定性是飞机机翼结构设计需要考虑的重要方面,是目前飞机结构设计中的一大难点.为了建立波纹板屈曲的数值分析方法,对复合材料波纹板在剪切载荷作用下的屈曲进行试验、工程算法和有限元方法研究.通过试验,探索了3种不同波长的复合材料波纹板的屈曲载荷和破坏载荷;采用工程算法和有限元方法计算了8种不同波长的复合材料波纹...     

辉光放电等离子体对边界层流动控制的机理研究

本文求解电位势方程得到电场分布,假定电场强度大于击穿阈值的区域为等离子体区,给定电荷密度,得到了作用于流体上的电场力。通过求解带源项的NS方程,研究了一个大气压下的均匀辉光放电等离子体对边界层流动的影响,考察了电场力做功对流动的影响。本文研究结果同文献[6]一致,即电场力总体上使边界层流动加速。另外,电场力做功对流动参数的影响可以忽略。文献[6]给出的线化电场模型和本文得到的电场相比,具有较大的差别,该差别引起了流动参数显著差别。     

劈绕射系数的一种通用表达式及应用

 在讨论矩形波导辐射场、飞行器翼边电磁散射场时,要用到二维金属尖劈的绕射系数表达式。许多学者对尖劈绕射场进行了详细讨论,给出了明确结果,其中级数形式的严格解收敛速度较慢,简化公式是当源置于无限远处取得的,对于尖劈近区场源问题没有触及,文献[2]给出的简化解析式在应用上有一定的困难。因此我们重新进行了推导,获得一种新的表达式,适合于近、远场源问题的数值计算。矩形波导棱边的绕射是普通劈尖的一个特例。