基于修正面元法的机翼静气动弹性计算

针对静气动弹性弱耦合解法多次迭代耗时较长的问题,提出一种斜率修正面元法计算气动力.该方法以机翼在不同迎角下的CFD数据作为修正基础,以直线段近似曲线段升力线,弥补面元法无法考虑非线性因素的缺陷.通过M6机翼气动力损失以及迎角补偿等静气动弹性的计算表明,相对于CFD技术,该方法保持了较高的准确度并且计算效率较高,适用于复杂结构优化.     

CFD在高速飞行器热气弹问题中的应用

采用气动力/热/结构耦合的方法对高速细长体飞行器结构热静气动弹性问题进行了研究.为保证耦合计算精度，达到准确预测热气动弹性特性的能力，气动力和气动热计算采用CFD数值模拟方法，热应力和热变形计算采用有限元方法并通过热考核试验验证.以该简单细长体飞行器模型为研究对象，对其热静气动弹性特性进行了计算与分析，计算结果表明:CFD/CSD耦合可准确模拟热气弹问题，且气动加热造成结构温升不均衡是结构变形的主导因素，力热耦合静气弹变形与单纯受力分析变形形式不同，对飞行器气动特性影响规律不同.准确预测飞行器热气动弹性特性对飞行器结构设计十分必要.      

大展弦比气弹机翼柔性非线性变形的气动力效应分析与风洞试验

从飞行器刚弹耦合动力学模型出发，引入柔性机翼准定常假设，建立大柔性飞行器非线性静气动弹性气动力方程，利用非线性迭代求解思路模拟了柔性飞行器的静气动弹性响应行为，开展了大展弦比飞机静气动弹性风洞试验验证，采用气动力有限基本解与机翼的耦合计算，发现了大柔性飞机大变形状态下载荷及结构变形形式随风速的变化规律.传统基于小变形假设的线性分析方法和刚体分析由于无法考虑气动面随结构变形的曲面气动力因素和结构变形后的非线性刚度特性，均与风洞试验存在一定的误差.对于大展弦比柔性飞机的非线性静气动弹性分析十分必要.      

带支撑装置的CRM-WBH模型流固耦合数值模拟

共同研究模型（common research model，CRM）机翼/机身/平尾组合体构型（CRM-WBH）是第4届AIAA阻力预测研讨会（Drag Prediction Workshop IV，DPW IV）选择的基准构型，DPW IV会议的统计分析结果显示计算结果与试验结果之间存在显著差异.采用CFD方法和流固耦合方法数值模拟了带支撑装置的CRM机翼/机身/平尾组合体模型（CRM-WBHS）的气动特性，主要目的是评估支撑装置和静气动弹性变形对CRM-WBH模型气动特性数值模拟结果的影响.通过与CRM-WBH模型CFD数值模拟结果和欧洲ETW风洞（European Transonic Wind Tunnel）测力、测压和模型变形测量结果的对比分析，表明模型支撑装置导致机翼上翼面激波位置前移，升力系数、阻力系数下降，俯仰力矩系数增加；静气动弹性变形主要影响机翼上表面激波位置和外翼处激波位置前负压，导致升力系数、阻力系数进一步下降，俯仰力矩系数进一步增加.CRM-WBHS模型的流固耦合数值模拟结果更加接近试验结果.      

三维机翼跨音速颤振的分布式并行计算

以多台微机组成分布式计算系统,利用无限插值方法(TFI)生成三维多块贴体运动网格,以Navier Stokes方程为控制方程,求解三维机翼的跨音速非定常气动力,并与颤振方程耦合迭代计算,求解机翼广义位移响应的时间历程,根据广义位移的时间历程的衰减、等幅和发散振荡等情况确定机翼跨音速颤振临界条件.通过算例验证,计算结果与实验结果和理论分析相吻合.     

面向飞机各设计阶段考虑静气动弹性效应的面元法飞行载荷分析方法

针对现代飞机设计对于高效率、高精度飞行载荷分析的迫切需求，建立了面向飞机各设计阶段考虑静气动弹性效应的面元法飞行载荷分析方法，并以此为基础构建了飞机各设计阶段气动弹性优化框架，以加快各阶段的迭代，提高设计效率.概念设计阶段的气动力分析采用低阶面元法，可适应该阶段为确定气动外形参数范围所需的大量计算.初步设计阶段采用高阶面元法开展气动力分析，以兼顾求解精度和效率的需求.详细设计阶段的飞行载荷分析方法在基本面元法的基础上引入了外部高精度CFD气动力或试验气动力，以进一步提高飞行载荷的精度并兼顾求解效率.研究结果表明，所建立考虑静气动弹性效应的面元法飞行载荷分析方法具有较强的工程实用性，可以满足现代飞机各个设计阶段的需要.      

CFD/CSD耦合飞行器多学科优化设计新方法

文章以飞行器巡航外形为设计对象, 构建了一种新的飞行器的气动和结构特性评估方法, 即结构模型反迭代方法.该方法较传统的松耦合静气动弹性方法效率提高了4倍以上.以此为基础建立了一种新的飞行器气动/结构耦合多学科优化设计框架, 将优化效率提高4倍以上.采用数值求解N-S方程和结构有限元方程方法作为气动和结构学科的分析工具, 保证了设计结果的可信性.算例表明以巡航外形作为设计对象能够获得与传统方法一致的飞行器气动与结构特性, 以此为基础开展的无人机气动外形优化设计也获得了良好的设计结果.      

吸气式高超声速飞行器机体/推进一体化内外流非定常耦合方法

基于吸气式高超声速飞行器机体/推进一体化的气动布局设计方式，文章提出了一种内外流一体化流场的耦合求解方法，其中燃烧室内流场采用考虑有限速率化学反应动力学模型的一维非稳态方法求解，进气道和尾喷管外流场采用二维CFD软件计算，进气道与燃烧室在耦合界面处通过一维平均方法实现静温、静压和Mach数等参数传递.并分别以日本国家航空与航天实验室（NAL）的氢燃料燃烧室模型作为内流场验证算例，以某典型高超声速飞行器一体化模型作为内外耦合流场验证算例.研究结果表明:有限速率化学反应准一维方法能较为准确地模拟燃烧室内燃烧流场，提出的内外流场耦合方法能够有效地计算出内外流耦合效应，计算后体压力分布与理论值较接近.该方法可为超燃冲压发动机的性能快速分析和吸气式高超声速飞行器机体/推进一体化的初步分析设计提供重要参考.      

激光辐照高超声速舵面热力响应研究

建立了基于流固耦合的高超声速飞行器舵面结构在激光辐照下静气动弹性模型,流体控制方程为三维雷诺平均N-S方程,分别采用了中心格式和AUSM+up格式对粘性项和对流项进行空间离散,时间推进采用了高斯-塞德尔隐式推进方法,湍流粘性系数求解使用Menter SST模型。利用冯·卡门研究所高超风洞实验结果对模型进行了校核,预估了激光辐照对高超舵面的热力影响。结果表明,气动力/热计算模型与实验数据符合较好,能够准确模拟高超飞行器的热力参数,根据模型外推结果,激光在较低功率下加热高超飞行器舵面可能导致舵面材料弹性模量大大降低,继而发生弯曲发散而折断,高超飞行器可能因此发生气动失稳而坠毁。     

三维拉氏理想磁流体数值模拟方法

在Z箍缩驱动ICF过程中，磁场流体耦合作用是其整个物理过程中非常重要的部分.针对Z箍缩过程中多介质、大变形的特点，发展了三维相容拉氏理想磁流体交错型以及单元中心型格式，两种格式均具有一阶时间与空间精度.通过数值算例，验证其精度和强壮性，并比较分析两种格式的特点.     

Aeroelastic analysis and nonlinear dynamics of an elastically mounted wing

The response of an elastically mounted wing that is free to plunge and pitch, supported by nonlinear translational and torsional springs, and interacting with an incoming stream is analyzed. A tightly coupled model of the wing flow interaction is developed. A three-dimensional code based on the unsteady vortex lattice method is used for the prediction of the unsteady aerodynamic loads. The response of the wing shows a sequence of static and dynamic bifurcations and chaotic motions when increasing the flow speed. Pairs of stable solutions are observed over the different response regimes. The effects of the gust and structural nonlinearity on the wing's response are also investigated. The results show that gust may lead to jumps between the pairs of solutions for static and dynamic equilibrium responses without impacting the boundaries of the different response regimes. As for the effect of the structural nonlinearity, increasing the nonlinear coefficient of the stiffness of the torsional spring yields lower static deflections and amplitudes of the limit cycle oscillations.     

Lyapunov-Schmidt reduction algorithm for three-dimensional discrete vortices

We address the persistence and stability of three-dimensional vortex configurations in the discrete nonlinear Schrödinger equation and develop a symbolic package based on Wolfram’s MATHEMATICA for computations of the Lyapunov-Schmidt reduction method. The Lyapunov-Schmidt reduction method is a theoretical tool which enables us to study continuations and terminations of the discrete vortices for small coupling between lattice nodes as well as the spectral stability of the persistent configurations. The method was developed earlier in the context of the two-dimensional lattice and applied to the onsite and offsite configurations (called the vortex cross and the vortex cell) by using semianalytical computations [D.E. Pelinovsky, P.G. Kevrekidis, D. Frantzeskakis, Physica D 212 (2005) 20-53; P.G. Kevrekidis, D.E. Pelinovsky, Proc. R. Soc. A 462 (2006) 2671-2694]. The present treatment develops a full symbolic computational package which takes a desired waveform at the anticontinuum limit of uncoupled sites, performs a required number of Lyapunov-Schmidt reductions and outputs the predictions on whether the configuration persists, for finite coupling, in the three-dimensional lattice and whether it is stable or unstable. It also provides approximations for the eigenvalues of the linearized stability problem. We report a number of applications of the algorithm to important multisite three-dimensional configurations, such as the simple cube, the double cross and the diamond. For each configuration, we identify exactly one solution, which is stable for small coupling between lattice nodes.     

前缘弯掠风扇叶尖涡的数值分析与PIV实验

本文采用三维Navier-stokes方程、κ-ε两方程湍流模型和SIMPLE算法,对空调用前缘弯掠开式斜流风扇转子叶尖涡的三维流动特性进行了数值分析与PIV实验测量。结果表明:沿子午面流动的主气流与从叶顶外侧吸入的气流之间的卷吸作用导致在吸力面叶顶区域产生了叶尖涡。叶尖涡从25％叶顶弦长位置到50％弦长位置不断发展,同时沿吸力面向下游移动;叶尖涡随后逐渐向压力面移动;叶尖涡从75％叶顶弦长位置至转子下游逐渐耗散,沿着一条与叶轮旋向相反的斜线向转子下游发展,大致持续到转子下游65％叶顶弦长位置。计算结果和PIV的实测结果吻合,表明用 CFD工具可以预测流场,同时为前缘弯掠开式斜流转子在大型中央空调室外机上的应用和风机系统的优化设计及其降噪提供了重要的内流数据。     

A Lattice Boltzmann Model for Fluid-Solid Coupling Heat Transfer in Fractal Porous Media

We report a lattice Boltzmann model that can be used to simulate fluid-solid coupling heat transfer in fractal porous media. A numerical simulation is conducted to investigate the temperature evolution under different ratios of thermal conductivity of solid matrix of porous media to that of fluid. The accordance of our simulation results with the solutions from the conventional CFD method indicates the feasibility and the reliability for the developed lattice Boltzmann model to reveal the phenomena and rules of fluid-solid coupling heat transfer in complex porous structures.     

Nonlinear analysis and enhancement of wing-based piezoaeroelastic energy harvesters

We investigate the level of harvested power from aeroelastic vibrations for an elastically mounted wing supported by nonlinear springs. The energy is harvested by attaching a piezoelectric transducer to the plunge degree of freedom. The considered wing has a low-aspect ratio and hence three dimensional aerodynamic effects cannot be neglected. To this end, the three dimensional unsteady vortex lattice method for the prediction of the unsteady aerodynamic loads is developed. A strong coupling scheme that is based on Hamming's fourth-order predictor–corrector method and accounts for the interaction between the aerodynamic loads and the motion of the wing is employed. The effects of the electrical load resistance, nonlinear torsional spring and eccentricity between the elastic axis and the gravity axis on the level of the harvested power, pitch and plunge amplitudes are investigated for a range of operating wind speeds. The results show that there is a specific wind speed beyond which the pitch motion does not pick any further energy from the incident flow. As such, the displacement in the plunge direction grows significantly and causes enhanced energy harvesting. The results also show that the nonlinear torsional spring plays an important role in enhancing the level of the harvested power. Furthermore, the harvested power can be increased by an order of magnitude by properly choosing the eccentricity and the load resistance. This analysis is helpful in designing piezoaeroelastic energy harvesters that can operate optimally at specific wind speeds.     

Phase transitions in a two-dimensional vortex system with defects: Monte Carlo simulation

AbstractThe phase states and phase transitions in a system consisting of a two-dimensional vortex lattice with defects are studied by the Monte Carlo method. It is shown that a “rotating lattice” phase, which is an intermediate phase between the vortex crystal and vortex liquid phases, is present. The dependence of the temperature of the transition from the rotating lattice phase into a vortex liquid on the strength of the defect potential is determined. The current-voltage characteristics of the system are calculated at various temperatures for point, square, and linear defects. It is shown that the phase state of the system strongly affects its transport properties.