基于流固耦合的刷式密封泄漏特性理论与实验

理论与实验两方面研究了刷式密封泄漏流动特性及影响因素。建立了考虑刷丝变形的刷式密封泄漏流动特性流固耦合理论模型,实验研究了刷式密封泄漏量随进出口压比和转速等的变化规律,并与理论求解模型对比验证。在此基础上,运用理论模型分析了不同结构参数对密封泄漏量的影响规律,验证了构造的刷式密封泄漏量理论计算公式的准确性。研究结果表明:径向间隙一定时,刷式密封泄漏量随压比的增大而近似线性增加,转速对刷式密封的泄漏量影响不大;刷丝直径和刷丝间隙通过影响气体在刷丝束区域有效流通面积来影响泄漏量,有效流通面积越小,泄漏量越小;反之,则越大;随着刷丝直径、刷丝轴向排数的增大,刷式密封泄漏量逐渐降低;随着刷丝间隙的增大,泄漏量逐渐增大;随着后挡板保护高度、刷丝与后挡板轴向间隙的增大,泄漏量先显著增大后缓慢增加,最后趋于稳定;所构造的泄漏量理论公式可以准确预测刷式密封的泄漏量,为刷式密封结构优化设计提供理论依据。      

基于欧拉旋转的在轨服务航天器姿态跟踪算法

针对在轨服务航天器对目标器交会捕获的任务要求,提出了一种基于欧拉旋转的姿态     

基于基尔霍夫面的四维声学频域公式

流动诱发噪声问题是实际工程领域极为普遍的问题之一,经典的声比拟模型仅以声学压力为参考来评估声场特征分布还远远不够。从声压和声学速度矢量为变量的四维线性波动方程出发,选择包围非线性声源的基尔霍夫面为积分面,并结合对流格林函数,给出均匀运动介质的四维声学频域积分公式;针对静止、旋转单极子源和偶极子源开展数值预测研究。结果表明：本文获得的声压和声学速度分布与解析解吻合,均匀来流情形下静止点源的声场分布表现出典型的对流效应;受均匀来流、点源的自激频率、谐波阶次和旋转频率的共同影响,旋转点源的声场分布则表现出明显的多普勒效应和对流效应。本文对均匀流矢量气动噪声开展的精细化研究能够为声能量的评估和传输路径预测提供技术支撑,进一步为降噪研究提供理论依据。     

环形燃烧室两相喷雾燃烧的大涡模拟

在三维任意曲线坐标系下采用欧拉-拉格朗日两相大涡模拟方法对环形燃烧室火焰筒气液两相紊流瞬态反应流进行数值模拟.采用椭圆偏微分方程生成三维贴体网格,计算中所采用的数学度模型有:k方程亚网格尺度模型估算亚网格紊流粘性;亚网格EBU燃烧模型估算化学反应速率;热通量辐射模型估算辐射换热.并在非交错网格体系下气相采用SIMPLE算法和混合差分格式求解,液相采用随机离散模型(Stachasttc Separated Flow,简称SSF),在拉格朗日坐标系下追踪各油珠群沿各自轨道运动、质量损失及能量变化.通过计算结果与实验数据相比较,表明在三维贴体坐标系下对燃烧室火焰筒两相紊流油雾燃烧流场进行大涡模拟,采用欧拉-拉格朗日两相大涡模拟方法能反映两相紊流化学反应流流动及实际燃烧过程.     

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

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

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

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

基于准定常的飞行器动导数的高效计算方法

结合工程设计经验和工程计算方法,针对飞行器动导数计算方法展开研究。基于求解准定常欧拉方程发展了一种适合于亚跨超声速状态下,任意气动外形气动导数的高效计算方法,并将计算结果与一般非定常方法的求解结果以及国际标准实验数据进行了对比,研究结果表明该方法结果可靠、效率较高。     

蜂窝芯层非线性等效弹性参数

考虑蜂窝壁板面内荷重对壁板弯曲的影响,对Gibson公式进行了修正,提出了一种蜂窝芯层的非线性等效拉伸弹性模量的拟合方法.讨论了蜂窝层合板的非线性响应,给出了x向拉伸和压缩、y向拉伸时的参数计算公式,以及某结构的应力一应变关系拟合结果,建立了蜂窝芯层非线性等效泊松比的计算模型.为蜂窝层合板的非线性动力响应分析提供了参考.     

三维机翼非定常数值模拟

本文研究了三维非定常Euler方程的显式有限体积解法,并编制了相应的计算程序。在程序中,网格生成部分采用的是代数方法生成贴体网格,再用求解椭圆型方程的方法对网格进行优化;在计算部分采用二步Lax-Wendroff格式的一种变体进行迭代,它具有二阶时间精度,并在预估步和修正步采用同样的处理过程,避免了对中间点的定义,使程序简洁有效率。全部计算结果均与有关实验和其它资料进行了比较。     

欧拉—庞卡莱方程在“和平号”太空站姿态恢复分析中的应用

1997年因"进步号"飞船与空间站对接时发生意外碰撞而使得姿态失去控制,从而导致太阳翼无法朝向太阳而获得能量。使用欧拉方法和四元数方法建立运动学方程,通过欧拉—庞卡莱方程建立动力学方程来分析其姿态恢复的运动过程,研究得出太阳翼能够重新得到光照的可行性方案。