半规则三角网格模型细分曲面重构

针对逆向工程中的三角网格重构问题,提出了一种保持尖锐特征的半规则三角网格模型细分曲面重构算法,以充分利用细分曲面的多分辨特性.首先经提取尖锐特征和删除最大独立点集得到基网格,然后利用插值Loop细分和最近点法向投影法不断调整半规则网格得到重采样网格,最后运用渐进插值(PI,Progressive Interpolation)算法生成半规则细分曲面.实例表明,重构后获得的细分曲面连接性好,可以直接进行细分小波分析.     

基于两相流理论滑动轴承动力特性求解

将计算流体动力学(CFD)两相流与动网格技术应用于滑动轴承动力特性数值求解,建立了基于CFD两相流滑动轴承动力特性求解模型,该模型无需设定油膜破裂边界条件且更能准确模拟滑动轴承流场特性.比较了单相流与两相流滑动轴承压力分布特性,计算分析了滑动轴承气穴分布特征及其影响因素,研究了两相流模型对滑动轴承动力特性的影响.计算结果表明:气化比例随着转速、偏心率和气化压力的增加而迅速增大,随进口压力的增加而缓慢减小.考虑两相流后,直接刚度系数增加,交叉刚度系数减小,直接与交叉阻尼系数均减小.随着偏心率的增加,单相流与两相流动力特性系数求解结果偏差增大.      

非结构混合网格高雷诺数流动计算中一种改进的LU-SGS隐式格式(英文)

The lower-upper symmetric Gauss-Seidel (LU-SGS) implicit relaxation has been widely used because it has the merits of less dependency on grid topology,low numerical complexity and modest memory require...     

二维平板水漂运动数值模拟

探索平板水漂运动的流体力学现象机理对飞行器着水问题的研究有重要的参考价值。基于有限体积法和k-ε RNG湍流模型求解非定常雷诺平均Navier-Stokes（URANS）方程,采用流体体积分数（VOF）模型与速度入口造波法构造数值波浪水池,并结合整体动网格方法,完成二维平板在静水面与波浪水面的水漂运动数值模拟。在与试验值和理论值对比的基础上,讨论初始姿态角、投掷角度与投掷速度对水漂运动的影响。进一步研究不同波浪参数和波浪位置对平板水漂的影响,并从能量角度展开分析。结果表明：20°姿态角平板能以最小的投掷速度实现水漂运动;运动中平板能量的相对损失率受初始投掷速度的影响较小,主要受投掷角和姿态角作用,随投掷角或姿态角的增加而增大。波浪情况下,在上行波位置触水的平板能够获得更大的接触面积,而更长的触水时间发生在波谷处;因此,这两个位置触水的平板在水漂运动过程中的相对能量损失大,其数值变化比波峰位置大5%左右;在上行波位置触水时,平板的速度衰减与相对能量损失随着波高的增加而增大,而波峰位置则有相反的变化趋势。     

倾转涵道倾转过渡阶段的非定常气动力

为研究倾转涵道动力装置在倾转过渡阶段的非定常气动力,使用基于滑移网格技术的非定常计算方法,利用基于内场信息和叶素理论的压力盘模型模拟风扇螺旋桨,通过求解Navier-Stokes(N-S)方程,对涵道风扇俯仰拉起过程进行数值模拟.结果表明:在倾转过程中涵道风扇非定常气动力的迟滞特性明显,俯仰角速度延迟了流场分离,增加了涵道风扇的升力和阻力;低速飞行时,涵道风扇在整个倾转过渡阶段气动性能优良;在高速大迎角飞行时,涵道风扇气动性能恶化,俯仰力矩曲线紊乱,不利于进行倾转过渡飞行.      

螺旋桨/机翼相互干扰的非定常数值模拟

针对太阳能无人机分布式螺旋桨滑流问题,采用基于结构/非结构混合网格的CFD方法进行了非定常数值模拟.通过截取翼段消除有限翼展三维效应的影响并简化研究对象,将滑移网格和转捩模型应用于数值模拟;分别研究了螺旋桨位于翼段前方和后方的情况,并与干净螺旋桨和干净翼段进行对比.结果表明:螺旋桨位于翼段前方或后方均使得翼段升力、阻力、低头力矩增加并呈周期性波动,同时螺旋桨位于翼段后方时对翼段的影响较小;而翼段的存在使得螺旋桨的拉力、吸收功率和效率增加并呈周期性波动,流场的非均匀性导致了螺旋桨振动.      

采用近似边界求解壁板颤振的CBS有限元算法

为了简化Euler方程求解壁板颤振问题,将近似边界条件成功地推广到CBS(characteristic-based split)有限元方法求解壁板颤振领域,把运动壁面的影响转化为Euler方程的边界条件,避免了对动网格的处理.通过对近似边界方法、动网格方法以及3阶活塞理论的数值模拟结果进行对比发现:首先,加载近似边界的CBS有限元方法操作非常简便,且模拟结果与动网格方法完全吻合;其次,Euler方程可以准确捕捉到壁板前后两端不光滑尖点所引起的压力突跳,而3阶活塞理论则无法做到;最后,采用近似边界的CBS有限元方法在求解无黏可压缩流动中细长体小变形的流-固耦合问题上具有潜力和优势.      

Delaunay graph-based moving mesh method with damping functions

The fluid–structure interaction and aerodynamic shape optimization usually involve the moving or deforming boundaries, thus the dynamic mesh techniques are the key techniques to cope with such deformation. A novel dynamic mesh method was developed based on the Delaunay graph in this paper. According to the Delaunay graph, the mesh points were divided into groups. In each group, a factor ranging from 0 to 1 was calculated based on the area/volume ratio. By introducing a proper function for this factor, this method can control the mesh quality with high efficiency. Several test cases were compared with other dynamic mesh methods regarding mesh quality and CPU time, such as radial basis function method and Delaunay graph mapping method.     

RBF整体网格变形技术与多体轨迹仿真

在多体分离的CFD动态轨迹仿真中,计算网格需要准确刻画固壁边界之间幅度很大的相对运动。理想的状况是采用整体网格,在不增删网格节点和不改变原网格拓扑结构的条件下,通过网格变形来刻画多体边界间的相对运动。为此,以先前发展的径向基函数(RBF)网格变形技术为基础,通过引入一种改进的拉普拉斯光顺算法,大幅改善了单独RBF网格变形技术在处理大幅度边界位移时的网格质量下降和网格拓扑破坏问题,成功实现了多体分离中的超大幅度网格变形。在改进的拉普拉斯算法中,针对四面体、三棱柱、六面体、金字塔四种网格单元类型提出了一套统一的网格质量评估标准,并将其用于非结构混合网格下的拉普拉斯光顺。通过耦合求解非定常Navier-Stokes方程和刚体六自由度运动方程,利用典型跨声速下机翼与外挂物分离问题对以上网格变形和光顺技术进行验证,将计算获得的外挂物投放过程的空间位置、姿态角和气动力与风洞实验结果进行了对比分析,结果显示:在不采用网格嵌套和局部重生的前提下,与改进的拉普拉斯光顺相结合的RBF网格变形技术可以有效地解决多体分离过程中的网格大幅度变形问题,且网格质量能够满足精确刻画外挂物投放飞行轨迹的要求。     

多激励下某直升机传动系统动载特性

针对由弧齿锥齿轮和行星轮系构成的直升机传动系统,构建了纯扭振动模型,采用集中参数法建立了齿侧间隙非线性动力学方程.通过有限元方法求得了时变啮合刚度,采用4-5阶变步长Runge-Kutta法对动力学方程进行了数值求解,借助动载系数、相图、Poincaré截面图、快速傅里叶变换频谱图等分析手段,研究了传动系统在时变啮合刚度、齿侧间隙、综合传动误差、外载荷等多种激励作用下系统的动载特性.结果表明啮合刚度对传动系统的影响最大,动载系数最大值为1.5;齿侧间隙对系统响应特性的影响是有限的;啮合误差在一定程度上抑制了齿轮系统的振动;外载荷波动对不同速级的影响不同,动载系数最大值发生在并车传动.