基于欧拉方程的二维振荡机翼非定常气动设计反命题方法

建立了基于欧拉方程组的二维振荡机翼非定常气动设计反命题方法的数学模型.通过一系列变换,将物理时空中的求解域转换成映象坐标系中的规范区域;并导出了映象坐标下的欧拉方程,结合非定常反命题的边界条件,便可用有限差分方法求解.本方法引入映象坐标系解决反命题边界形状的不确定性,并能利用欧拉方程现有的各种差分格式于反命题求解.     

螺旋桨滑流与机翼气动干扰数值模拟研究

分别采用非定常方法和激励盘理论研究螺旋桨滑流对机翼的干扰作用,并对计算结果进行对比分析。非定常方法基于动态面搭接结构化网格,运用有限体积法求解非定常RANS方程。对于激励盘理论,首先采用点对点结构化网格,随单独螺旋桨旋转的方式进行非定常模拟得到螺旋桨桨叶压力分布,然后将桨叶压强分布转化为激励盘载荷分布,最后将载荷分布作为边界条件对干扰流场进行定常模拟。数值结果表明:除了在较大来流攻角时有所差异,两种方式得到的机翼气动力吻合良好;非定常方法可以得到螺旋桨与机翼的瞬时干扰流场,而激励盘理论则可以大大缩短数值模拟的时间并降低计算要求。     

变形飞机动态气动特性数值模拟研究

根据近年来国内外提出的折叠翼变形飞机方案,设计了一种折叠翼变形飞机的3D数模,利用作者发展的动态混合网格技术和非定常计算方法,对该折叠翼变形飞机机翼折叠和展开过程中的动态气动特性进行了数值模拟,并与准定常计算结果进行了对比分析.计算结果表明:机翼在展开/折叠的过程中,飞机的气动特性会发生剧烈变化,在变形飞行控制率设计时需要引起高度的重视.     

超声速机翼－尾翼组合体的非定常气动力边界元法

本文以Ｇｒｅｅｎ函数为基础的边界元，应用时间历程法，计算机翼－尾翼组合体的非定常气动性能。这为飞行器的设计及其非定常气动特性分析提供了一种有效和完善的计算方法。本文介绍的方法，特别适用于研究、分析各类飞行器非定常动态响应过程的数值计算，例如飞机对突风动态响应等。本文提供的三维机翼的超声速启动过程非定常气动力数值结果，与Ｌｏｍａｘ理论解作了比较，可以看到吻合得很好，说明本文所使用的方法是可行的。对于机翼－尾翼组合体的数值结果，则进一步提供了尾翼受到机翼的洗流影响过程。     

基于非定常面元/黏性涡粒子法的低雷诺数滑流气动干扰

针对太阳能无人机螺旋桨滑流与机翼的气动干扰,考虑了低雷诺数流动下气体黏性和压缩性影响,并根据黎曼边界条件和涡量等效原则建立了能够快速计算分析螺旋桨-机翼气动干扰的非定常面元/黏性涡粒子的混合方法。首先使用有试验数据的风洞模型以及数值模拟技术对混合方法进行验证,在此基础上研究了不同安装位置与工况下螺旋桨与机翼的气动干扰。结果表明:螺旋桨对轴向气流的加速以及滑流诱导的上洗和下洗效应使机翼气动力呈现出增升增阻的现象,机翼升阻比有所下降。较大的弦向间距以及较高的垂直安装位置在减缓机翼升阻比下降的同时也使得螺旋桨拉力有所减小。对于多个螺旋桨的气动干扰,不同的桨叶旋转方向导致机翼气动力不同的变化规律,当旋转方向与机翼翼尖涡反向时,螺旋桨滑流能够抑制翼尖涡的强度,提高机翼气动效率。     

螺旋桨滑流与机翼气动干扰的非定常数值模拟

基于动态面搭接网格,求解非定常雷诺平均Navier-Stokes方程模拟螺旋桨滑流与机翼的干扰流动.分别对3种构型,即:单独螺旋桨十短舱构型、螺旋桨十短舱十机翼构型以及短舱十机翼构型进行流场数值模拟,对比分析了机翼对螺旋桨滑流流动结构的干扰,在滑流作用下机翼气动特性的改变以及机翼对螺旋桨桨叶气动力的影响.模拟结果表明:...     

基于黏性伴随方法的跨声速机翼气动优化设计

 进行了基于黏性伴随方法和Navier Stokes方程的跨声速机翼气动优化设计研究。分别推导了适用于三维跨声速机翼气动反设计和减阻设计的黏性伴随方程、边界条件和梯度求解表达式，并研究了伴随方程的数值求解方法。通过将网格生成、流场计算、黏性伴随方程数值求解、梯度求解和拟牛顿优化算法等几方面的有效结合，发展了一种跨声速机翼气动优化设计方法。为了提高计算效率，将多重网格方法应用到方程的数值求解中来加速收敛。跨声速机翼反设计和减阻设计算例验证了本文所发展的方法的正确性。采用本文的方法进行优化设计，一般通过20~30次迭代就能得到满意的结果。     

基于直接力控制的阵风响应及阵风减缓研究

考察某无人机模型的阵风气动响应特性和基于直接力操纵的阵风减缓效果.在求解非定常Euler方程时引入“网格速度”方法模拟阵风边界条件,通过动态嵌套网格技术实现舵面运动.首先对NACA0006翼型迎角阶跃型阵风的气动力响应进行了验证计算,计算结果与理论结果、文献计算结果吻合良好.在此基础上对某无人机模型在迎角阶跃型、One-minus-cosine型阵风作用下的气动力响应过程进行了数值模拟,并分析评估了不同舵面运动方式对飞机气动性能的影响.最终研究比较了不同舵面运动方式下阵风气动响应的减缓效果.结果表明:通过合理设计舵面运动,可以有效达到抑制阵风引起的非定常气动干扰的目的.     

弹性飞机跨声速机动载荷计算方法

基于跨声速非定常气动力求解的重叠场源法,开展了弹性飞机跨声速机动载荷计算方法研究,为现代飞机结构强度设计提供更加可靠和精确的临界载荷计算方法。首先通过采用重叠场源法求解关于62%根弦俯仰振荡的LANN机翼在马赫数为0.822的非定常气动力并与实验结果进行对比,验证了重叠场源法对跨声速激波效应预测和非定常气动力计算的能力;其次应用频域气动力有理近似技术拟合场源法计算得到的广义气动力系数矩阵,建立了机动载荷分析的状态空间模型;然后完成了某型民用飞机俯仰机动载荷分析,研究了俯仰机动飞行情况下飞机机体状态量及飞机部件载荷响应规律。计算结果表明：考虑机体弹性变形后,机翼和平尾气动载荷响应最大值分别减小了5.1%和10.6%,升降舵气动载荷响应最大值增大了16.2%,在飞机结构强度设计中必须考虑机体弹性效应对飞机部件载荷的影响。     

低雷诺数螺旋桨滑流对机翼气动特性的影响研究

深入研究低雷诺数滑流对机翼的影响,能够推进临近空间低速流动机理性研究,提供可靠的气动参数。参考某太阳能无人机,建立单螺旋桨计算模型,采用两叶螺旋桨,通过ICEM网格软件生成具有两个计算域的高质量结构网格,应用滑移网格边界条件,对模型进行数值模拟;分析低雷诺数螺旋桨滑流的发展和机翼在滑流作用下的非定常气动特性,研究不同螺旋桨位置对机翼气动特性的影响,计算结果表明螺旋桨滑流会很大程度地改变机翼表面压力分布和沿翼展的升力分布,对机翼升阻特性有显著影响,同时螺旋桨滑流可以抑制机翼表面层流分离泡的产生。     

基于非定常面元/黏性涡粒子混合法的旋翼/平尾非定常气动干扰

旋翼/平尾非定常气动干扰是导致直升机纵向“抬头（Pitch-up）”现象的主要原因。为在直升机设计阶段准确分析旋翼/平尾非定常气动干扰以及由此引起的低速纵向操纵特性变化,通过涡量等效原则和Neumann物面边界条件建立了适用于旋翼/平尾气动干扰分析的非定常面元/黏性涡粒子混合法。该方法耦合了考虑尾迹时变效应的非定常面元法、黏性涡粒子法及涡量镜面法,以准确模拟旋翼和平尾的非定常气动载荷、旋翼尾迹的非定常特性以及旋翼尾迹对平尾的气动干扰效应。首先通过计算NASA ROBIN（Rotor Body Interaction）旋翼尾迹几何和诱导速度分布,并与实验测量值、时间精确自由尾迹及CFD计算结果对比验证方法的准确性。相比于时间精确自由尾迹,本文方法计算精度更高。随后分析了旋翼/平尾非定常气动干扰对平尾向下气动载荷和气动导数的影响,并分析了平尾构型对旋翼/平尾非定常气动干扰的影响规律。分析表明：旋翼尾迹与平尾干扰导致低速状态的平尾载荷突增,气动导数反号;低平尾气动载荷突增较大,高平尾较小,但高速气动导数反号;前置平尾载荷突增量减小,但对应速度范围较宽;右旋直升机右平尾载荷突增量较小,但气动导数特性基本不变。     

基于CFD方法的机翼阵风响应研究

通过引入＂网格速度＂方法模拟阵风条件,求解非定常欧拉方程实现了不同展弦比0012平直机翼阵风响应的数值模拟。首先采用该方法对NACA0006翼型迎角阶跃型阵风的气动力响应进行了计算,计算结果与文献结果、理论值吻合良好。进一步对展弦比分别为5、10的0012平直机翼在迎角阶跃型、One-Minus-Cosine型阵风作用下的气动力响应过程进行了模拟分析。研究结果表明,当展弦比增大时,阵风作用下机翼的升力系数响应会增加,机翼翼根部位气动响应幅值大于翼尖部位响应特性。     

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

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

微型飞行器低雷诺数矩形扑翼非定常气动特性的数值模拟

采用双时间步方法求解三维可压缩非定常N-S方程,数值模拟了微型飞行器低雷诺数矩形扑翼的非定常绕流,首先将得到的结果与文献进行了对比,数据间具有较好的一致性.然后针对不同的展弦比、减缩频率及初始攻角,计算了矩形扑翼的非定常气动特性及表面流态和动态压力分布,并分析了翼尖涡对扑翼非定常气动特性的影响.     

螺旋桨滑流对带后缘襟翼机翼气动特性影响的数值分析

采用动量理论结合有限体积法求解N-S方程的方法,进行了螺旋桨滑流对带后缘襟翼机翼气动特性干扰的数值模拟分析.螺旋桨模型简化为一个激励盘,即一个无厚度可穿透的圆盘,以给定的圆盘拉力载荷,即压强差来模拟螺旋桨滑流的效果.采用嵌套网格技术处理桨盘与机翼间复杂的几何关系.算例计算表明,该计算方法能够较好的模拟滑流对带后缘襟翼机翼气动特性的影响.     

边条翼作俯仰运动时翼面吹气的作用

在３ｍ低速风洞中通过测力与微丝显示方法，研究了在边条翼的边条和主翼上吹气对机翼作大攻角快速俯仰运动时空气动力特性的影响，同时还研究了在不同迎角下开始吹气的作用。结果表明，翼面吹气能有效提高机翼的非定常气动特性和缩小机翼的非定常气动特性迟滞回环，特别是在小迎角下开始吹气效果较好。     

大尖削机翼非定常跨音速绕流的数值解法(英文)

本文给出一种计算飞机机翼上定常和非定常跨音速气动力的数值解法,使用了一种特殊设计的坐标变换,是用时间精确交替方向隐式(ADI)有限差分算法来求解非定常跨音速修正三元小扰动位势方程。给出了F5战斗机机翼的数值结果并与XTRAN3S,ATRAN3S及试验结果进行了比较,表明本方法是有效的和经济的。     

机翼和叶栅非定常流的Kutta条件

本文在现有的实验研究和数值计算的基础上，对经典Ｋｕｔｔａ条件在非定常流动中的适用性作了进一步探讨，并对文献〔１，２〕中提出的一种新的Ｋｕｔｔａ条件－Ｇｉｅｓｉｎｇ－Ｍａｓｋｅｌｌ（Ｇ－Ｍ）模型进行了分析与比较，本文给出了二维振荡机翼和叶栅绕流的变域变分有限元解法中Ｋｕｔｔａ条件的处理方法。本文旨在为机翼和叶栅非定常绕流的数值计算提供更合理可靠的尾缘边界条件。     

共性半环翼跨介质航行器变体气动特性研究

基于一种特殊的双层半环形闭合翼构型,利用流体仿真软件对该构型跨越水空介质后变体过程中的不同状态进行了仿真,并对其非定常气动特性进行了分析;比较了特定条件下非定常状态和定常状态的气动参数随折叠角的变化曲线,得到了机翼展开和回收至不同角度的压力和速度分布云图以及对应位置的压力分布曲线;研究了不同变体速率对机翼变体的非定常气动特性的影响,并提出了附加动力学影响、流场迟滞影响等物理效应来解释该构型的非定常气动特性的形成原理。仿真结果表明,机翼回收过程的气动参数大于定常状态,展开过程气动参数小于定常状态,变体速率越大,非定常效应越明显,产生这种差异的原因是来自于流场迟滞的影响。     

基于伴随方法的机翼多设计点气动反设计方法

针对三维机翼多点多约束气动反设计问题,基于伴随理论方法和粘性流雷诺平均N-S方程,通过粘性流数值模拟、伴随方程与梯度精确数值求解、计算网格高效算法及梯度类优化算法等有效结合,并采用考虑多设计点梯度权重系数的并行计算近似模式,开展了一种三维机翼多设计点多约束气动反设计方法研究,进行了典型算例验证。研究表明:所发展的机翼多设计点气动反设计方法具有较好的鲁棒性及优化效率。