非局部核与应变局部化问题解的相关性

在非局部理论框架下,研究了应变局部化现象模拟结果与所用非局部核函数的相关性。首先,在非局部损伤理论框架中,提出一个以应变与损伤参量为独立变量的新的自由能泛函形式,并由热力学定律,推导非局部弹性损伤本构方程。然后,通过二次线性化的方法,给出一维情况下的支配方程。最后,给出对应不同形式非局部核的损伤局部化现象模拟结果。通过与模拟结果的比较,分析非局部核与解的相关性,得到了一些有价值的结论。     

数值模拟有外挂物的空腔流动

数值模拟了有外挂物存在的开式空腔的超音速流动。数值计算采用全隐式无分裂方法,湍流模型采用修正的B-L代数模型。研究了有外挂物存在的空腔流动的剪切层自持性振荡,并与无外挂物存在时的现象进行类比。压强与时间的变化曲线、空腔底部流向的静压系数分布以及空腔后壁面法向的静压系数分布表明,空腔内的压力脉动与外挂物所处位置是有关的。对于Rossiter振荡模态,数值计算结果能够与实验结果相吻合。     

一类基于非等距单元平均重构的高效差分格式

构造了一维非线性双曲型守恒律的一类基于非等距单元平均值的点值重构的高精度高分辨率守恒型差分格式。其构造思想是:首先,将计算区间划分为若干个互不重叠的小区间,再根据格式精度的要求利用Gauss-Lobatto点和Gauss-Chebyshev点划分小区间,通过各非等距细小区间上的单元平均值,重构各细小区间交界面上的点值,并加以校正;其次,利用近似Riemann解计算细小区间交界面上的数值通量,并结合高阶Runge-Kutta TVD时间离散方法,得到了高精度的全离散方法。证明了该格式的无振荡特性。然后,将格式推广到一维双曲型守恒方程组情形。最后,给出了几个标准数值算例,验证了格式的高效性。     

A transition prediction method for flow over airfoils based on high-order dynamic mode decomposition

This article presents a novel approach for predicting transition locations over airfoils, which are used to activate turbulence model in a Reynolds-averaged Navier-Stokes flow solver. This approach combines Dynamic Mode Decomposition (DMD) with e^N criterion. The core idea is to use a spatial DMD analysis to extract the modes of unstable perturbations from a steady flowfield and substitute the local Linear Stability Theory (LST) analysis to quantify the spatial growth of Tollmien–Schlichting (TS) waves. Transition is assumed to take place at the stream-wise location where the most amplified mode’s N-factor reaches a prescribed threshold and a turbulence model is activated thereafter. To improve robustness, the high-order version of DMD technique (known as HODMD) is employed. A theoretical derivation is conducted to interpret how a spatial high-order DMD analysis can extract the growth rate of the unsteady perturbations. The new method is validated by transition predictions of flows over a low-speed Natural-Laminar-Flow (NLF) airfoil NLF0416 at various angles of attack and a transonic NLF airfoil NPU-LSC-72613. The transition locations predicted by our HODMD/e^N method agree well with experimental data and compare favorably to those obtained by some existing methods (LST/e^N or $\gamma -\stackrel{￣}{R{e}_{\theta \text{t}}}$). It is shown that the proposed method is able to predict transition locations for flows over different types of airfoils and offers the potential for application to 3D wings as well as more complex configurations.     

共轴刚性旋翼桨毂阻力特性及流动机理

共轴刚性旋翼直升机在高速飞行时,桨毂流动复杂、分离强、阻力大。为明晰其阻力特性和流动机理,采用CFD方法针对已完成风洞试验的共轴桨毂组合模型进行数值模拟研究,获得了桨毂组合模型各单独部件的阻力、表面流动和空间流场特征,阐明了产生阻力最大的部件和影响阻力的主要因素,揭示了中间轴整流罩和塔座设计参数的减阻机制。分析结果表明:上、下旋翼桨毂是产生阻力的主要部件;中间轴和塔座的分离尾流对桨毂表面流动产生较大的干扰作用,使桨毂整流罩表面受干扰区域产生气流分离;具有较缓和逆压梯度的中间轴整流罩和塔座能有效减小分离尾流对桨毂整流罩的干扰,从而降低整个共轴桨毂系统的阻力。     

变来流速度下旋翼翼型非定常气动特性分析

为分析变来流速度状态下的旋翼翼型气动特性,提出了利用翼型平移来模拟来流速度变化的数值方法.在此方法基础上,采用基于隐式LU-SGS(lower upper symmetric Gauss-Seidal)方法的非定常雷诺平均N-S(Navier-Stokes)(RANS)方程,模拟了SC1095旋翼翼型在定迎角 变来流速度及变迎角 变来流速度状态下的非定常气动特性.通过对比分析发现:翼型在变速度-定迎角状态下会表现出明显的非定常现象,产生了前缘分离涡,气动特性会出现明显的迟滞效应及波动现象,脉动速度越大,非定常效果越明显.并且基准速度越大,翼型气动特性的峰值越大;翼型迎角越大,非定常涡出现的也越早.考虑直升机旋翼翼型实际工作环境,在变速度-动态失速状态下,翼型最大迎角处的气动力会得到一定程度的削弱,在小迎角下的气动力得到一定程度的增强,且脉动速度越大,翼型的非定常特性也越强.      

倾转旋翼机多部件对机翼气动干扰的分析及优化

为研究倾转旋翼机各部件对机翼的气动干扰效应,建立了一套适用于倾转旋翼机流场CFD求解方法,提出并生成了一套由三棱柱/四面体组成的适用于倾转旋翼机多种飞行状态的非结构混合网格系统.以三维Navier-Stokes方程为主控方程,采用动量源方法进行倾转旋翼的模拟,并引入了高效的OpenMP并行加速技术等,提出了一套新型的动量源项添加及搜索方法.计算分析了倾转旋翼机旋翼/机身/短舱对于机翼气动特性的干扰影响,得出了一些对设计有指导意义的结论.采用基于径向基函数（RBF）的代理模型方法,根据机翼不同展向位置的干扰程度对各段翼型配置进行气动优化.优化结果表明:考虑气动干扰作用,在巡航速度下优化后的全机升阻比增长达到了36.78%.      

一种求解欧拉方程改进的通量分裂方法

本文以Steger-Warming提出求解Euler方程的通量分裂方法为基础,利用改进的特征值分裂方法,建立了一种连续的通量分裂方法。该方法克服了Steger-Warming通量分裂法存在的分裂后的通量在个别点导数不连续的问题,同时方法简单,便于推广。对于通量分裂后的Euler方程,采用迎风格式,建立了相应的二阶隐式有限体积方程。为了减少求解所得块矩阵方程组所需要的工作量,本文利用特征变换,在不改变原方程组离散化精度的前提下,将所得块矩阵方程组化为分离的代数方程组进行求解,从而大大减少了计算量。文中分别对一维管流,二维绕NACA-0012和NACA64A-10翼型在亚声速和超声速来流的跨声速流场进行了计算,并将所得结果与有关实验和计算结果进行了比较,均符合良好。 本文所建立的方法具有占机内存小,计算时间短的优点。在普通中小型计算机上亦可使用。     

管道中复杂马赫反射流场的高精度数值模拟

本文分析了一类高精度迎风TVD欧拉方程差分格式的精度,并用240×80个网格点进行马赫反射的计算,获得了目前实验和计算不易得到的两个激波三叉点的复杂波系结构。本文还研究了激波管上壁面对马赫反射的影响,并计算了轴对称扩张管道中激波的运动及激波与半球形障碍物的相互作用。计算结果分析指出,本文所采用的高精度迎风TVD格式与计算方法可以用来模拟具有复杂波系的工程问题研究。     

叠层环板轴对称自由振动的精确解

抛弃任何有关位移或应力模式的人为假设,从三维弹性力学基本方程出发,导出横观各向同性环板轴对称振动问题的状态方程。在多种边界条件下,给出单层、叠层环板轴对称振动问题的精确解。