一种渐变无反射边界条件的建立与验证

建立并验证了一种适用于非线性流动的新无反射边界条件,以简单的迁移方程模型为基础,提出了一种通量渐变消失无反射边界条件。通过乘以合适的函数,在渐变层内偏微分方程中的空间偏导项逐渐消失,从而在边界点上得到只含有时间变量的常微分方程,来解决自由区域处边界造成的定解问题。将渐变边界条件推广到了一维和二维非线性Euler方程,其方程形式简单,便于实现。在计算域边界处采用空间滤波和二维情形时角区内引入牛顿阻尼区提高了渐变边界条件的数值稳定性。在非线性情况下对二维Euler方程渐变边界条件进行了一系列的数值验证,得到了令人满意的数值结果。     

轴向受载的Bernoulli-Euler薄壁梁的弯扭耦合动力响应

建立了一种普遍的解析理论用于研究确定性载荷作用下轴向受载的单对称Bernoulli—Euler薄壁梁的弯扭耦合动力响应。首先通过直接求解轴向受载的单对称均匀Bernoulli-Euler薄壁梁单元弯扭耦合振动的运动偏微分方程，给出了计算其自由振动的精确方法，并导出了轴向受载的Bernoulli—Euler薄壁梁自由振动主模态的正交条件。然后利用简正模态法研究了确定性载荷作用下轴向受载的Bernoulli—Euler薄壁梁的弯扭耦合动力响应，该梁所受到的载荷可以是集中载荷或分布载荷。最后假定确定性载荷是谐波变化的，得到了各种激励下封闭形式的解，并针对具体算例讨论了动力弯曲位移和扭转位移的数值结果。     

无碰撞激波的理想MHD模拟

采用了WENO格式数值求解一维理想MHD方程组,模拟了行星际无碰撞激波,研究了垂直无碰撞激波与行星际反向磁场结构和高密度等离子体团的相互作用过程,并与粒子模拟的结果进行比对,两者的结果非常类似.模拟结果表明,对涉及无碰撞激波的大部分现象中理想MHD模拟是准确且可行的,同时相对于粒子模拟又有很好的计算效率,便于扩展至二维或三维的情形.     

Rankine涡作用下激波变形和分叉的数值研究

数值模拟了激波和Rankine涡相互作用,重点研究了波涡相互作用引起的激波结构的变化过程.应用非结构化适应网格下的二阶精度Godunov型的PLM格式来求解Euler方程,对一平面激波与不同涡强Rankine涡相互作用下瞬态激波的运动、发展进行了研究.漩涡与激波间相对强度的大小对激波结构的影响明显,弱相互作用诱导激波变形,强相互作用诱导激波变形并产生分叉.计算结果表明,应用非结构化适应网格的PLM格式所得到的数值解较好地反映了瞬态激波结构的有关特征和信息.     

径向基函数在三维Euler方程数值计算中的应用

基于三维Euler方程,对非结构四面体网格给出了一种基于紧支径向基函数重构的ENO型有限体积格式,方法的主要思想是先对每一个四面体单元构造插值径向基函数,而在计算交界面的流通量采用高斯积分公式以保证格式的整体精度,时间离散采用三阶TVD Runge-Kutta方法。最后用该格式对一些典型算例进行了数值模拟,结果表明该方法计算速度快,对间断有很好的分辨能力。     

运输机短舱挂架纵向位置优选数值研究

针对运输机的翼/身/架/舱干扰问题,本文基于Euler方程,采用笛卡尔直角网格,应用有限体积法对全机绕流流场进行数值计算,分析和比较了短舱挂架不同前伸量、下沉量的影响,提出了基于最佳前伸量和下沉量的优选配合位置,并与初始位置的全机气动特性进行对比。结果表明,获得的短舱挂架优选配合位置使全机气动特性得到了明显改善和提高。     

Radical recombination in a hydrocarbon-fueled scramjet nozzle

To reveal the radical recombination process in the scramjet nozzle flow and study the effects of various factors of the recombination, weighted essentially non-oscillatory(WENO)schemes are applied to solve the decoupled two-dimensional Euler equations with chemical reactions to simulate the hydrocarbon-fueled scramjet nozzle flow. The accuracy of the numerical method is verified with the measurements obtained by a shock tunnel experiment. The overall model length is nearly 0.5 m, with inlet static temperatures ranging from 2000 K to 3000 K, inlet static pressures ranging from 75 k Pa to 175 k Pa, and inlet Mach numbers of 2.0 ± 0.4 are involved.The fraction Damkohler number is defined as functions of static temperature and pressure to analyze the radical recombination progresses. Preliminary results indicate that the energy releasing process depends on different chemical reaction processes and species group contributions. In hydrocarbon-fueled scramjet nozzle flow, reactions with H have the greatest contribution during the chemical equilibrium shift. The contrast and analysis of the simulation results show that the radical recombination processes influenced by inflow conditions and nozzle scales are consistent with Damkohler numbers and potential dissociation energy release. The increase of inlet static temperature improves both of them, thus making the chemical non-equilibrium effects on the nozzle performance more significant. While the increase of inlet static pressure improves the former one and reduces the latter, it exerts little influence on the chemical non-equilibrium effects.     

预估-校正LU-SGS的隐式算法

为了进一步提高Euler方程求解效率,在传统LU-SGS格式基础上提出了预估-校正LU-SGS。结合高分辨率迎风通量格式,发展了一套高效的Euler方程求解器。以NACA0012翼型和RAE2822翼型跨音速无粘流动作为算例,对比分析了预估-校正LU-SGS、传统LU-SGS和GMRES格式的计算效率。结果表明:预估-校正LU-SGS算法能显著提高求解效率,其效率接近GMRES算法,比传统LU-SGS方法提高了近三倍。     

A New Single-blade Based Hybrid CFD Method for Hovering and Forward-flight Rotor Computation

A hybrid Euler/full potential/Lagrangian wake method,based on single-blade simulation,for predicting unsteady aerodynamic flow around helicopter rotors in hover and forward flight has been developed.In this method,an Euler solver is used to model the near wake evolution and transonic flow phenomena in the vicinity of the blade,and a full potential equation(FPE) is used to model the isentropic potential flow region far away from the rotor,while the wake effects of other blades and the far wake are incorporated into the flow solution as an induced inflow distribution using a Lagrangian based wake analysis.To further reduce the execution time,the computational fluid dynamics(CFD) solution and rotor wake analysis(including induced velocity up-date) are conducted parallelly,and a load balancing strategy is employed to account for the information exchange between two solvers.By the developed method,several hover and forward-flight cases on Caradonna-Tung and Helishape 7A rotors are per-formed.Good agreements of the loadings on blade surface with available measured data demonstrate the validation of the method.Also,the CPU time required for different computation runs is compared in the paper,and the results show that the pre-sent hybrid method is superior to conventional CFD method in time cost,and will be more efficient with the number of blades increasing.     

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

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