基于Gappy POD的二维叶片压力反设计

采用离散型函数和B样条函数两种参数化方法生成叶片快照,研究比较基于单快照叶片和双快照叶片的Gappy POD拟合能力。引入流场求解,发展一种基于迭代偏差修正的Gappy POD反设计方法,并开展二维叶片表面压力分布的反设计研究。结果表明:双快照叶片方法能够有效提高Gappy POD的拟合空间和拟合精度,基于迭代偏差修正的反设计方法能够显著改善Gappy POD反设计方法的精度。     

基于POD的受激励后向台阶分离流低阶模型

在对后向台阶流场进行合成射流激励并研究不同激励频率对流场发展影响的过程中,发现流场在低频激励条件下与中高频条件下表现完全不同.为了详细分析这一现象,使用本征正交分解法（proper orthogonal decomposition,POD）将由PIV方法测得的流速分布数据进行分解,采用相位平均手段对含湍流动能较大的主要模态间的关系进行分析,并使用主要模态对流场的主要运动形式进行还原.结果表明,流场在各条件下的主要运动形式均可用少量低阶模态加以基本概括,低频激励下低阶模态相图近似于Lissajous图形,并描述了剪切层在激励作用下的摆动过程.      

槽道流POD重构及湍流动能耗散率分析

采用二维粒子图像测速仪（2DPIV）对槽道内涡波流场进行实验研究,用POD技术对2DPIV瞬态速度矢量场进行主导模态重构,得到槽道内的平均流速和湍流动能分布;采用大涡PIV方法对湍流动能耗散率分布进行计算.结果表明：重构流场表征了原始流场的主导结构,剔除了噪声等干扰信息;大涡PIV方法能有效地估算动能耗散率的分布;湍流动能在壁面附近较小,在接近槽道中心区域湍流动能越来越大,呈现出射流的特征;动能耗散率的峰值出现在壁面附近和槽道中心区域,动能耗散率随着远离壁面程度的增加先降低后逐渐增加直至达到峰值.     

对称槽道涡波流场流动特征的POD分析

为深入分析层流状态下对称槽道内涡波流场的流动特性及其变化规律,对流场进行了二维粒子图像测速(2DPIV)测量获取瞬态速度矢量数据,利用本征正交分解(POD)技术进行模态分解以及涡波流场的重构,然后根据重构的流场对对称槽道内涡波流场进行了平均速度剖面、流场脉动强度以及特征点的速度和频谱分布等方面的分析。结果表明:POD的前15阶模态能够表征涡波流场的主导结构,第1,3阶模态主要表现为一对旋向相反的涡对特征,第2阶模态具有涡旋和波状主流的特征;提取了5个涡旋涡核的位置作为流场流动特性的特征点;根据POD重构流场分析发现流向平均速度呈抛物线形状分布,法向平均速度呈对称分布特征;流向脉动强度受壁面的影响较大,法向脉动强度呈现抛物线形状分布;距离中心主流较近的1#,4#,5#特征点的速度脉动程度受主流的脉动强度影响较大,速度的脉动主频0.15 Hz与次频、流场的自然频率0.35 Hz共同影响特征点的速度分布;2#,3#特征点的流向速度呈衰减趋势,法向速度在初期幅度变化较大。     

飞行器大攻角复杂流动的POD和DMD对比分析

基于非结构/混合网格、耗散自适应2阶混合格式以及脱体涡模拟（detached eddy simulation,DES）方法开展了现代战斗机模型复杂分离流动的数值模拟,并与有限的平均气动力试验数据进行了对比,结果表明计算具有合理性,在此基础上进一步应用本征正交分解（proper orthogonal decomposition,POD）和动力学模态分解（dynamic mode decomposition,DMD）方法对数值模拟流场的非定常特性进行了对比分析.研究表明飞行器背风区流场由一对边条涡的螺旋运动主导,旋涡破裂前在横向空间截面上流场是中性稳定的,同时主涡核的运动是多频耦合的.POD和DMD的对比分析则表明:两者模态配对的方式不同,但主要模态之间具有一定相关性;POD模态中包含多种频率的运动,而且能量较集中于主模态,流场重构效率更高;DMD则将流场的主要特征运动提取为一些单频模态的组合,同时能够给出模态的稳定性.      

基于POD模态分解的液环泵瞬态气液两相流分析

针对液环泵内气液两相流动的复杂时空规律,采用本征正交分解(POD)方法对其瞬态气液两相流场进行特征分解,分析其相态场、速度场的空间基模态特征及模态系数的时域特征,建立非定常流场降阶模型,并对流场进行预测分析.结果表明POD方法可实现对液环泵内复杂流场的时空解耦分析,相态场及速度场的各阶模态系数在时域内的变化能够反映各阶...     

大雷诺数下方柱绕流PIV试验及数值模拟

基于风洞试验利用粒子图像测速技术(PIV)研究了大雷诺数(Re=3.42×10~4)下方柱绕流流场,并应用本征正交分解(POD)法分解PIV瞬态速度场,获得了方柱尾流的低阶模态和重构流场结果,研究方柱绕流的流动机理.采用γ-Re_θ和k-ωSST两种湍流模型研究了大雷诺数下方柱尾流的速度分布以及湍动能变化。风洞PIV试验可准确的观测方柱尾流脉动流场,基于POD方法得到的前6阶模态占总能量的78.7%,前6阶模态即可较准确重构流场并捕获流场的主要特征。通过对比,γ-Re_θ模型是较优的湍流模型,可较好地得到尾流的湍动能及其变化规律。     

基于自适应本征正交分解混合模型的跨音速流场分析

采用非线性模型替代线性回归模型响应本征正交分解(POD)基函数的系数,并采用自适应抽样方法确定快照集合,实现了基于自适应POD混合模型的跨音速叶片复杂流动分析及流场拟合.首先通过比较基于线性回归模型和非线性回归模型的基函数系数响应精度,验证非线性回归模型的收敛性和精确性;之后通过与静态抽样方法进行对比,研究分析自适应抽样技术的优越性;最后开展基于自适应POD混合模型的全三维跨音速流场分析及流动拟合.结果表明,采用自适应POD混合模型,不仅能够清晰地识别三维跨音速流场中的敏感流动特征,还能精确地拟合设计空间内任意状态的流场及出口气动参数.     

风力机翼型在失速工况下非定常流场的本征正交分解分析

风力机气动性能受静态失速与动态失速影响很大,对风力机翼型的失速问题研究具有重要意义。本文通过计算流体力学方法得到的风力机翼型在固定大攻角工况,以及大攻角震荡工况下的非定常流场,来研究翼型静态失速与动态失速。采用本征正交分解方法(POD),对非定常流场降阶,得到流场的POD模态以及对应的系数。POD模态结果表明在静态失速下,主要非定常流动结构是尾迹区域交替脱落的涡结构;在动态失速下,除了尾迹区域,前缘和整个吸力面都存在流动分离结构。     

一种基于模态域波束形成的水平阵被动目标深度估计

水面水下目标分辨与识别一直是被动声呐探测领域的难题.利用一种水平阵模态域波束形成算法获得己知方位目标声源的各阶模态强度,将其与不同深度的各阶参考模态强度进行匹配,最终实现了对声源的深度估计.仿真结果表明,该算法可以在信噪比为-10 dB的情况下,用300Hz带宽的信号样本,实现对声源深度的有效估计.系统分析了不同参数和不同波导条件对该方法目标深度估计性能的影响.其中,阵元数越多,模态样本数越多,计算频段越宽,方位估计精度越高,有效阵长越长,深度估计的性能越好.阵元间距和波导深度的变化不会影响该方法的深度估计性能,并且该方法的深度估计性能在声速剖面、海底参数等波导条件存在扰动时具有鲁棒性.     

风力机翼型气动计算的降阶模型研究

将基于POD的降阶模型应用于风力机翼型的气动研究。首先应用CFD数值模拟得到一系列快照结果;应用基于本征正交分解(POD)的方法得到流场的一组基模态,认为对于所研究的问题,任一流场可以由这些基模态通过线性表达得到;对控制方程进行Galerkin投影,得到降阶模型,将离散求解N-S方程的问题转化为一组只有十几个自由度的常微分方程,从而减少计算时间,提高计算效率。并对二维翼型的绕流的定常和非定常问题进行了分析,计算结果表明,降阶模型可以较好地捕捉流动的特征,与直接CFD模拟相比计算精度相当,但大幅有效地提高了计算效率。     

Natural pixel decomposition for interferometric tomographic reconstruction from ill-posed data

Reconstruction of three-dimensional refractive-index fields from multidirectional holographic interferograms can be a powerful flow diagnostic tool in many engineering applications. It is noninvasive and can capture gross fields instantaneously. Interferometric tomography for reconstructing flow fields, however, confronts ill-posed problems, that is, last of sufficient data. As the flow fields are sampled in the presence of opaque models and enclosures, the data are sparse, nonuniform, and sometimes incomplete in projection and scanning. Here, a new method, termed the variable grid method, to incorporate the distinct features of interferometric data for flow field reconstruction is developed and tested. The method is based on natural pixel decomposition. It utilizes variable rectangular grid elements in representing a field with a finite series. The variable grid method allows better reconstruction resolution with greater accuracy in the region where sampled rays are more closely crossed. This region is usually of interest in engineering.     

DMD和POD对超燃冲压发动机凹腔流动的稳定性分析

开式凹腔作为超燃冲压发动机中增加掺混和稳焰的装置, 其流动稳定性的研究对深入理解凹腔增加掺混和稳焰机理以及凹腔的设计有着重要的学术意义和工程应用价值.基于大涡模拟方法对超燃冲压发动机开式凹腔流动进行数值模拟, 分别采用动力学模态分解(dynamic mode decomposition, DMD)和本征正交分解方法(proper orthogonal decomposition, POD)对自激振荡流动进行稳定性分析. DMD方法可准确提取凹腔的振荡频率, 与Rossiter模型以及压力脉动FFT分析得到的频率吻合较好, 且DMD中对应Rossiter前3阶频率的模态在流动中的主导作用顺序也与FFT分析结果一致, 自激振荡中RossiterⅢ模态占据主导作用, 同时DMD方法对Rossiter 3阶以上模态频率的预测能力明显强于FFT分析方法.在对低频的提取方面, DMD方法比Rossiter模型更具有优势.与前6阶Rossiter模态对应DMD模态均缓慢收敛, 主要表现为剪切层中的分离涡结构和中部及下游区域中的涡结构.前3阶不稳定模态中的分离涡结构主要集中在中部剪切层以及后缘附近区域. POD方法中较少的模态包含流场绝大部分的能量.但是, 通过POD方法提取的模态频率在分辨率上效果不佳, 提取到最低频率为Rossiter 3阶模态对应的频率, 且模态中均存在次频, 次频与主频之间的耦合导致模态的形态相差较大.另外, 与DMD方法相比POD方法无法判断所提取的模态的稳定性.      

3-D Phase Object Reconstruction Using Beam-Deflection Iterative Tomography

３－ＤＰｈａｓｅＯｂｊｅｃｔＲｅｃｏｎｓｔｒｕｃｔｉｏｎＵｓｉｎｇＢｅａｍ－ＤｅｆｌｅｃｔｉｏｎＩｔｅｒａｔｉｖｅＴｏｍｏｇｒａｐｈｙ￥ＧＡＯＹｉｑｉｎｇ（ＮａｎｃｈａｎｇＩｎｓｔｉｔｕｔｅｏｆＡｅｒｏｎａｕｔｉｃａｌＴｅｃｈｎｏｌｏｇｙ，Ｎａｎｃｈ...     

一种基于单水听器的浅海水下声源被动测距方法

针对浅海波导中宽带脉冲声源的被动测距问题,本文在群延迟理论的基础上,与warping变换处理相结合,提出了一种适用于浅海波导中宽带声源的单水听器被动测距方法.利用warping变换可以实现对脉冲声源接收信号各阶简正波的分离提取,对分离后的简正波进行时频分析处理可以得到各阶简正波到达时刻和频率之间的关系,即各阶简正波的频散曲线,从而得到任意两阶简正波到达接收水听器的时延差.海底相移参数P是描述海底地声参数的一个重要参量,包含了海底地声参数信息,在海底环境参数未知而P已知的情况下,利用P和简正波水平波数之间的关系可以求得任意两阶简正波的?S_(g,mn)(群慢差).根据群延迟理论,利用得到的任意两阶简正波的时延和?S_(g,mn)可实现利用单水听器对水下声源进行被动测距.本文提出的测距方法测量简单、计算方便,具有较强的实用意义.数值仿真和海上实验数据处理结果的测距误差都在10%以内,证明了该方法的有效性.     

Waveguide mode solver based on Neumann-to-Dirichlet operators and boundary integral equations

For optical waveguides with high index-contrast and sharp corners, existing full-vectorial mode solvers including those based on boundary integral equations typically have only second or third order of accuracy. In this paper, a new full-vectorial waveguide mode solver is developed based on a new formulation of boundary integral equations and the so-called Neumann-to-Dirichlet operators for sub-domains of constant refractive index. The method uses the normal derivatives of the two transverse magnetic field components as the basic unknown functions, and it offers higher order of accuracy where the order depends on a parameter used in a graded mesh for handling the corners. The method relies on a standard Nyström method for discretizing integral operators and it does not require analytic properties of the electromagnetic field (which are singular) at the corners.     

基于本征正交分解的气动优化设计外形数据挖掘

气动外形的全局优化设计会产生大量的过程数据,其中隐含的设计知识具有较高的挖掘价值.数据挖掘有助于获取直观、可定性描述的设计知识.本文采用基于本征正交分解的数据挖掘方法从气动优化设计的过程数据中获取设计知识,数据挖掘对象为跨音速压气机转子叶片NASA Rotor 37的优化过程数据,该数据由基于粒子群方法的绝热效率最大化优化设计产生.结果表明:基于本文数据挖掘方法获取的设计知识能够直接反映气动外形的变化规律,为叶片的气动外形设计提供参考;数据挖掘的设计知识成功地验证了优化设计结果的有效性.     

基于闪烁体3D剂量测量系统中光场重聚焦位置的标定方法

放疗中基于闪烁体光场成像的3D剂量测量方法是,利用光场相机对射线在闪烁体沉积能量所产生的闪烁光进行拍照,给出每一个聚焦层的光分布后进行三维重建,进而得出剂量的三维分布。光场相机所获取图像的重聚焦位置影响闪烁光三维重建的准确性,作者提出利用棋盘格标定板进行重聚焦位置标定的方法。使用Lytro Illum光场相机拍摄置于已知位置的棋盘格标定板,经光场数据处理得到焦点堆栈序列,利用刃边法得到所有重聚焦像的高斯离焦参数σ,经比较后获得最小σ对应的重聚焦参数α,最终给出若干α和重聚焦位置的对应关系,最终完成重聚焦位置标定。与常用的刻度尺标定方法相比,高斯离焦参数减小2%～10%,有利于放疗3D剂量测量精度的提高。