壁面展向匀速运动的槽道湍流的直接数值模拟

为了解展向存在匀速流动时湍流槽道流动的特性,应用直接数值模拟方法研究了壁面的展向运动对湍流槽道流动的影响.采用标准谱方法对固定质量流率条件下,上壁面展向运动速度分别等于零、流向流量平均速度以及两倍于流向流量平均速度等3种情况进行了直接数值模拟.结果表明,壁面的展向运动会导致表面摩擦力的增大,二阶湍流统计量诸如均方根速度脉动、雷诺剪切应力和均方根涡量脉动等都会随着壁面展向运动速度的增大而显著增长,同时,槽道截面上的均方根速度脉动也变得更均匀.      

基于涡方法生成大涡模拟进口条件的数值计算

探索改进涡方法来生成大涡模拟的非定常进口条件.改进方法中,为了避免局部漩涡数量过多或者过少,采用密度分布方式放置漩涡场.并用Langevin方程控制漩涡运动,模拟实验方法中的蜂窝器,使改进后的涡方法生成的脉动速度场更加符合湍流的特征.在已知雷诺平均的流场结果下,利用涡方法产生漩涡场,进而生成能满足大涡模拟所需要的非定常进口流场.为了检验改进后的涡方法在生成脉动速度过程中的效果,在槽道中进行了对比数值试验,并借助于直接数值模拟数据做对比,对比分析槽道进出口的平均速度、涡量以及雷诺应力统计,证明改进后的涡方法在生成大涡模拟进口条件下是非常有效的.      

槽道内涡波流场展向涡的分布特征

通过2DPIV测量了槽道内涡波流场的瞬态速度矢量场,利用雷诺分解和旋转强度的λci准则进行了展向涡的识别;并从识别的逆时针和顺时针展向涡的尺度属性、力学属性和运动属性方面分析了展向涡的分布特征。结果表明,流场中的顺时针展向涡的数量大于逆时针展向涡的数量;在法向位置0.05~0.45范围内,逆时针展向涡的旋转强度呈抛物线形状变化,顺时针展向涡整体呈下降趋势;逆时针展向涡的平均直径和椭圆平均长轴长度随着法向位置的增加而缓慢增加,而顺时针展向涡趋于下降;逆时针展向涡的椭圆长轴倾角最大值为67.28°,顺时针展向涡的等效椭圆的离心率呈下降趋势,似圆性较好;逆时针展向涡对全局平均的〈vωz〉分量在法向位置0.15处出现最大值61%;展向涡对流场雷诺应力的影响较小;逆时针展向涡的数量密度在0~0.08法向位置范围内呈急剧上升趋势,并达到最大值0.36,顺时针展向涡数量密度整体呈倒U形分布;不同旋向的展向涡占总展向涡数量密度的比值都大于0.3;展向涡的流向对流速度整体上小于流场的平均速度,且在槽道中心区域差值较大;展向涡的法向对流速度随着法向距离增加,从负值上升到正值,在中心主流两侧分别发生喷射和扫掠2种现象。     

高压高剪切率下密封环缝隙流体动态润滑特性

为减少高压高剪切率下热效应对密封环缝隙流体润滑性能的影响,以密封环缝隙流体剪切速度梯度、壁面剪切应力和油膜温升表征摩擦副润滑性能,采用RNG k-ε湍流模型,通过Workbench建立流-固-热多场模型,计算得到不同主轴转角下缝隙流动形态变化规律和温升前后油膜厚度变化量、不同L型槽内外径比及长径比下流速流型分布规律、壁面剪切应力、温度分布和热变形量。研究结果表明:密封环L型槽黏性底层随剪切速度增大而增厚,油膜厚度随运动周期增加而变薄;当L型槽内外径比小于1.07时,随比值减小密封环壁面剪切应力不断增大,最大变化率为7%;密封环L型槽长径比在0.19时平均壁面剪切应力达到最大,当长径比继续增大时,油膜区剪切速度梯度逐渐减小,油膜温升和热变形亦随之减小。研究结果可为马达优化以减少能量损失和改善密封环润滑条件提供理论指导和依据。      

同轴旋转圆台环隙流动机制及实验研究

同轴旋转圆台环隙内流体的流动是基于经典的两同轴旋转圆柱环隙内泰勒库特流的一种扩展研究。通过流动显示实验和PIV粒子图像测速技术对圆台环隙的内部流场进行可视化和定量化的研究,分析涡运动的周期性规律,探究圆台环隙内雷诺应力分布和水位高度对流场的影响。研究表明,随着时间的发展,涡列周期性明显且整体在下移,当脉冲数为100脉冲/s时,环隙内形成均匀分布的正反交替的涡;当脉冲数为200~500脉冲/s时,存在快慢交替的分裂周期;在3种水位高度下也都存在明显周期性分裂,只是周期时间和涡的个数不同;平均流场存在上凸型外向流和下凹型内向流2种流态,流态的差异是离心力与静压力双重作用的强弱变化所致;雷诺应力分布中,径向正应力占主导,并主要集中在环隙中部。      

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

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

基于DFT的水射流红外热像频域时空分析

基于二维离散傅里叶变换及空间频谱分析,对水射流湍流脉动的空间尺度进行了研究,得到了由红外辐射温度表征的被动标量湍流场在对流区、耗散区、惯性子区的特征空间尺度及其时间演化规律.对射流不同区域的关心点重新采样,得到湍流场中关心点的时间序列,利用一维离散傅里叶变换,分析了对流区大尺度涡中心、惯性子区小尺度涡中心、耗散区及射流轴心线上各关心点的湍流波动特征.计算了时间序列频谱的分形维数,研究了自由湍流不同尺度区间上述各关心点湍流脉动的分形特征.      

基于共振三波的大涡模拟壁面模型

大涡模拟的一个重要发展方向是处理工程流动问题,但高雷诺数下近壁湍流边界层对计算网格的过高要求成为制约其应用的主要瓶颈;要减少计算量,构建合理且能够正确反映近壁湍流动力学特性的壁面模型是一条有效途径.在此背景下,构造了一种计算中便于实现的基于共振三波的壁面模型.通过对充分发展槽道湍流算例的研究,完成了对该模型构造方法、方案的论证和验证计算.结果表明使用壁面模型后,在减少计算量的同时,定性和定量的结果均较为合理,能够得到重要的瞬态特征结构和正确的一阶、二阶统计量,从而部分验证了该方法的有效性.     

向阳面磁顶区可压缩涡旋诱发重联模拟研究:通量传输事件的形成和结构

本文用二维可压缩MHD模型,模拟研究了向阳面磁顶区的涡旋诱发重联过程。重联过程的结果形成了同心的磁流体涡旋。同时还对可压缩模型与不可压缩模型的异同点作了对比研究。结果表明,在同样条件下,可压缩情况的涡旋发展的增长率比不可压缩情况慢;结构特性无显著变化。本文指出,涡旋诱发重联是产生通量传输事件的一种重要机制。通量传输管可看作为磁流体涡旋管。对通量传输事件诸物理量(磁场、速度、密度、压力和温度)的分布特性作了模拟研究,可解释FTEs的主要观测特性。      

耀斑后环珥系速度场的定量研究

本文研究了在太阳重力、偶极磁压强梯度和气体压力梯度联合作用下环珥系内物质的下落运动,并在等温、定常假设下用数值方法计算了该环系的二维视向速度场。通过将速度场的理论计算跟文[1]中给出的观测速度场拟合,导出了该环系的有关参数及分布。计算表明环系内物质密度和磁场强度对物质下落运动影响较显著。     

Numerical study of magnetic reconnection possibly occurring near heliospheric current sheet

The third-order accurate upwind compact difference scheme has been applied for the numerical study of the magnetic reconnection driven by a plasma blob impacting the heliospheric current sheet, under the framework of the two-dimensional compressible magnetohydrodynamics. The results show that the driven reconnection near the current sheet could occur in about 10–30 min for the interplanetary high magnetic Reynolds number, R_M = 2000–10,000, a stable magnetic reconnection structure can be formed in hour order of magnitude, and there appear some basic properties such as the multiple X-line reconnections, vortex structures, filament current systems, splitting and collapse of the high-density plasma blob. These results are helpful in understanding and identifying the magnetic reconnection phenomena possibly occurring near the heliospheric current sheets.     

应用修正湍流模型对旋转通道内换热的模拟

针对旋转光滑矩形通道分别应用针对旋转状态修正的 k-ε、标准k-ω以及提出的针对旋转状态修正的k-ω湍流模型进行流动和换热的数值模拟,通过与实验结论的对比,讨论了采用不同湍流模型对计算结果的影响.计算工况为旋转数Ro=0.24,流体进口雷诺数Re=25000.计算结果表明:采用所提出的针对旋转状态修正的k-ω湍流模型的计算结果要比采用针对旋转状态修正的k-ε以及标准k-ω湍流模型的计算结果更接近实验结论.     

超临界压力RP-3在竖直细圆管内混合对流研究

研究了超临界压力下碳氢燃料航空煤油RP-3在竖直细圆管内混合对流,分析了浮升力及热物性对碳氢燃料在垂直管中对流换热的影响。实验中控制热流密度从200~500 kW/m2变化,进口压力变化范围为3~5 MPa,进口雷诺数从5 000~10 500范围内变化。研究表明:在向上流动情况中进口段存在较为明显的入口效应,换热出现恶化现象,而在向下流动中未出现;对于向上和向下流动,由于热物性的综合影响,换热系数沿流动方向增大;在较低进口雷诺数(Re=5 700)时,对于向下流动,随着浮升力影响的增大,浮升力改变了流体径向速度分布,出现了换热强化;在较高进口雷诺数(Re=10 500)时,浮升力对换热的影响依然显著;判别式Bo*数小于5.6×10-7未能预测浮升力对碳氢燃料换热影响。      

Wobble of a spin stabilized satellite with cross products of inertia and magnetic attitude control

Asymmetrical spin stabilized satellite dynamics in the vicinity of the required motion is considered. The principal axis of the maximum moment of inertia slightly deviates from its assumed direction in the satellite reference frame. This is formalized in the cross products of inertia. This inertial uncertainty results in a wobble, that is undesired angular velocity components perpendicular to the rotation axis, and oscillations of this axis near the required direction. The torque-free motion is investigated first. Expressions that explicitly relate satellite inertia parameters to wobble are provided. Wobble evolution under the action of magnetic damping control is analyzed next. Its amplitude approximate exponential decay behavior and residual unavoidable wobble level are derived. These expressions are compared with numerical simulation results of nonlinear equations of motion including various disturbance sources.     

PANS方法在双圆柱绕流数值模拟中的性能分析

PANS方法以传统的雷诺平均Navier-Stokes(RANS)方程为模板,引入控制参数对原方程进行修改,是一类能够有效模拟分离湍流流动的混合方法.采用基于Menter 剪切应力输运(SST)湍流模型的PANS方法,选取双圆柱绕流作为研究对象, 分析双圆柱流动结构,验证方法的准确性与可靠性,考察PANS方法中模型参数选取对于计算结果的影响,评估不同类型PANS方法的性能.研究表明:模化湍动能比例这一模型参数对计算结果有显著影响,且计算域内统一的模型参数取值难以处理复杂流动问题.在此基础上使方法中模化湍动能比例取值随流场物理信息和网格尺度变化,所得结果与SST DES计算结果及实验结果吻合较好,表明该方法适用于复杂湍流数值模拟.