层流平板摩擦阻力的数值计算

为提高层流摩擦阻力数值计算的精度,采用数值求解N-S方程的方法,对层流平板摩擦阻力进行了计算.计算中采用4种先进上风格式,同时配合以14种不同的壁面法向网格分布和5种常用限制器,并将计算结果与布拉修斯解进行对比.结果表明,van Leer的FVS(Flux Vector Splitting)格式数值耗散大,对网格分布和限制器敏感,不适合用于摩擦阻力计算.限制器和计算格式的粘性分辨率越高,摩擦阻力的计算精度越高.在计算量允许的情况下,可以减小壁面法向第1层网格高度,以提高计算精度,但应同时兼顾网格长宽比.对于摩擦阻力计算,边界层内应至少布置4个网格点.      

大长宽比单元有限元误差分析及高精度计算

针对大长宽比单元将传统的基于单元整体尺度的有限元误差估计扩展为基于单元不同方向尺度的误差估计;根据不同方向误差应该同量级的思想,得到了误差匹配准则.该准则可以作为网格划分的判据.根据这一匹配准则,提出了一种提高计算精度的方法——单向高次插值法.数值实验表明该误差估计是正确的,误差匹配准则作为网格划分的判据是有效的,采用合理的有限元插值逼近能够有效地减小计算误差,提高计算的精度和效率.      

高超声速热流高精度数值模拟方法

通过高超声速下钝锥热流的计算,对影响热流计算精度的因素进行了综合研究.获得了壁面法向网格Re数对热流计算的影响规律.采用边界高阶插值,提高了热流计算精度.用Roe和AUSM(Advection Upstream Splitting Method)+格式搭配minmod和混合limiter两种限制器,探讨了空间离散方法及限制器的耗散性对热流结果的影响.应用五阶WENO(Weighted Essentially Non-Oscillatory)方法,对高阶格式计算热流的性能进行了细致研究,最终确立了网格依赖小,热流精度高的计算方法.研究认为:高阶格式不但可以放宽热流计算对网格Re数的限制,同时使得格式对本身计算方法的依赖性也变小.      

高超声速再入钝头体表面热流计算

基于跟踪流线的轴对称比拟法,采用纯工程算法、Euler数值计算与边界层内工程算法相结合的方法,对高超声速再入钝头体的表面热流进行了计算,并将计算结果与风洞实验数据进行了对比,两者吻合较好,验证了两种工程算法在计算高超声速飞行器热环境方面的正确性.将两种工程算法与数值求解N-S方程进行对比,表明工程算法在迎风面的热流计算方面有较高的精度,节约了计算时间,很好地满足高超声速飞行器概念研究和设计的需要.      

喷流干扰气动热数值模拟的若干影响因素

工程上对喷流气动热进行近似计算时,常常对真实的全化学反应流近似处理,这时需要了解喷流气体和喷管外形的近似方法对计算结果的影响规律。采用数值试验对比分析的方法,探讨了近似计算的2个影响因素:其一是使用量热完全的空气喷流近似燃气喷流进行数值模拟时,喷管几何形状对热流分布的影响;其二是使用空气喷流近似燃气喷流进行数值模拟的方法,分析不同的喷流热力学参数匹配方法对喷流形态及热流分布的影响。通过对不同方法进行对比计算,揭示了各个方法对喷流干扰气动热计算的影响规律,可以用于指导工程分析。      

高超声速二方程湍流模型的数值模拟对比

对基于Reynolds和Favré(密度加权)混合平均的二方程湍流模型进行了修正,同时根据数值模拟高超声速流动时必须具有高分辨率捕捉间断面与在边界层内抑制数值粘性的能力的要求,提出了新的总变差减小(TVD)格式熵修正函数.在此基础上,通过对压缩拐角的高超声速湍流的数值模拟,对基于Reynolds和Favré混合平均的二方程湍流模型,以及其它不可压缩模型及可压缩性修正模型进行了对比,显示了不同湍流模型及可压缩性修正在计算壁面压力分布和热流分布上的特点,说明了对高超声速压缩拐角型流动,湍流模型可压缩修正的必要性,得到了基于Reynolds和Favré混合平均的二方程湍流模型的计算结果最接近实验结果.      

再入飞行器表面磁流体发电装置数值模拟

提出了钝锥型再入飞行器上的表面磁流体发电装置方案,采用低磁雷诺数近似下的磁流体力学模型对其进行了流动与静电场耦合的数值模拟,由此建立了表面磁流体发电装置中的物理图像.数值模拟结果表明,在典型再入条件(飞行高度46 km,速度7 km/s)下,所提出表面磁流体发电装置方案能够实现兆瓦级(电功率1.28MW)能量输出,电能提取导致再入飞行器阻力增大13.7%,对飞行器壁面总热流值的影响并不显著,但发电区域及其下游壁面的热流密度分布发生明显变化,该区域内热流密度峰值发生于电极的前、后缘.     

TVD格式与高超音速气动加热数值模拟

为了在化学非平衡流动中获得准确的流场解以及表面热流分布,将总变差减小TVD(Total Variation Diminishing)格式中的熵修正函数,由各向同性分布改为各向异性分布,同时让熵修正函数中的参数与流场中的压力梯度分布相关.将改进后的熵修正函数运用到高超音速化学非平衡绕流流场的数值模拟中,获得了较使用原有熵修正函数更为准确的流场参数和表面热流分布.采用改进的熵修正函数,可以提高壁面附近的粘性分辨率,降低热流计算结果对壁面附近法向网格尺度的敏感性.      

太阳电池阵极月轨道在轨热分析

以极月轨道上太阳电池阵作为研究对象,详细分析了该轨道上电池阵所得各外热流的确定方法.以此为基础,对太阳电池阵进行在轨温度数值模拟,得到月球轨道上电池阵所得外热流以及温度的周期分布.计算结果与地球卫星的相应数据进行比较,明确了月球轨道上电池阵所得外热流以及温度的周期分布特点,为月球卫星上太阳电池阵以及整个卫星的设计提供了有利的数据参考.      

矩形通道中亚尺度肋片的流动换热数值分析

对装有不同结构亚尺度肋片矩形通道的流动和换热进行了数值模拟,获得了通道流场、温度场分布以及平均努塞尔数的基本特征,对各种肋片通道的换热特性进行了对比分析.计算结果表明:亚尺度肋片距主肋片根部越近越有利于散热;在远离热面区域,亚尺度肋片的比表面积越大传热效果越好.当扩展表面积相同时,亚尺度肋片的长宽比越大换热效果越好;计算结果及分析揭示了控制肋片最优几何形状的统一原则——广义温度梯度均匀化原则,并依照此原则确定了可以强化换热的较优的肋片结构.计算分析表明广义温度梯度均匀化原则可以控制对流换热过程的传热强化.      

混合重叠网格插值方法的改进及应用

混合重叠网格间合理插值是保证流场计算正确的基础之一。本文针对重叠区网格尺寸匹配度较差时插值误差较大问题,发展了一种新型混合重叠网格插值方法,通过使用二阶精度插值、按单元尺度区分和扩充模板,改善了插值区网格尺寸匹配度较差时的插值精度。本文方法适用于任意单元类型的混合网格重叠,对各类单元处理透明,实现简单。计算结果表明,采用本文方法在网格重叠区流场变量传递正确有效,插值区网格尺寸匹配度较差时,相比原始方法,等值线过渡更为光滑,变量经过插值区耗散更小,计算与试验值符合更好。     

风扇/增压级带间隙三维粘性流场数值模拟

使用三维数值模拟软件NUMECA(Numerical Mechanics Application)对双涵道风扇/增压级带转子叶尖间隙的粘性流场进行了三维定常流动数值模拟,获得了其设计转速下的流场特性.结果显示,随着外涵背压的升高,风扇顶部激波后阻塞区周向和轴向尺寸增大,激波变形加剧,继续提高外涵压比,则外涵效率和进出口流量下降.内涵背压的升高,使内涵各转子叶片顶部的泄漏流动影响范围向上游扩散.内涵背压过高,效率和进出口流量同样会有下降的趋势.同时通过和实验结果的比较表明,数值模拟软件模拟的风扇/增压级的三维定常粘性流场能够较真实的反映出转子叶尖间隙泄漏流场的结构特征,以及间隙对流场参数的影响.      

涡轮流量计转子内完全三元流场的变域变分有限元解

三、涡轮流量计仪表特性曲线的计算由动量矩原理可知,当转子处于平衡时,可得到力矩平衡方程 T_d=∑T_i=T_h T_b T_t T_m T_P T_f下面分别讨论各力矩的大小。     

低速流预处理Roe格式中的数值粘性

使用Roe格式计算低马赫数流动,需采用预处理技术.在预处理技术中,需要选择适合的预处理参数,但有时很难找到,因而得不到合理的计算结果.通过数值实验,发现使用预处理方法时,格式的数值粘性对低速流动的模拟有显著的影响.提出了在预处理Roe格式中加入可调参数,用以控制数值粘性对计算结果的影响.带有可调节数值粘性参数的预处理Roe格式的数值粘性更接近低速流动的真实耗散,反映了低速流动的真实物理现象.对方腔流动和圆柱绕流的低雷诺数定常流动所做的数值实验,得到了比较满意的结果,表明这一改进是有效的.对同一低雷诺数不同低马赫数的定常圆柱绕流所做的数值实验,均得到了满意的结果,表明了方法的适用性.     

Taylor-Couette流的DSMC数值模拟

使用圆柱坐标系网格的三维DSMC(Direct Simulation Monte Carlo)模型,对同轴圆柱间Taylor-Couette流中Taylor涡的形成进行了数值模拟,并分析了不同计算域和边界条件下稳定流场中Taylor涡的轴向排列结构.在网格设置和流场参数不变的情况下,使用轴对称DSMC模型对Taylor涡进行数值模拟,所得Taylor涡的稳定过程与三维结果一致,验证了使用三维DSMC方法来解决微尺度低速的稀薄气流问题的可行性.三维模拟结果表明Taylor涡以较大的圆周速度绕轴旋转,二维模拟则无法体现.对不同的内圆柱旋转速度进行数值模拟,确定了能够产生Taylor涡的临界速度值.     

高超声速翼面气动热与静气动弹性综合分析

研究了考虑热效应影响的高超声速飞行器的静气动弹性问题.建立了考虑热效应影响的高超声速飞行器从气动热计算到静气动弹性分析的快速分析方法,该方法采用分层求解思路对热气动弹性问题进行解耦,利用边界层外的无粘数值求解和边界层内的工程算法相结合的方法,计算高超声速飞行器的表面热流,获得模型的温度场,并基于该温度场计算结构热刚度矩阵.以高超声速飞行器小展弦比翼面为研究对象,计算了该模型在3种飞行状态下的温度场、热应力、热刚度和静气动弹性特性,从而验证了热静气动弹性快速分析方法的可行性和适用性.     

用GAO-YONG湍流模式对翼型跨音粘流的数值模拟

应用GAO-YONG可压缩湍流模式,数值模拟了NACA0012,RAE2822翼型的定常跨音速粘流算例.对流项采用三阶ROE格式,扩散项采用二阶中心格式,用多步Runge-Kutta显式时间推进法求解空间离散后的控制方程.计算结果很好地预测了翼型表面的压力系数的分布、激波的位置、马赫数等值线的分布等情况,并且对翼型表面激波与边界层相互干扰以及层流向湍流的转捩问题进行了分析计算.计算结果与实验值符合很好,表明GAO-YONG可压缩湍流模式应用合适的计算方法能够高精度模拟翼型跨音粘性流动问题,并且基于GAO-YONG可压缩湍流模式各向异性湍流粘性的机理,提供了一种预测转捩起始位置的判别方法.