空腔流动的大涡模拟及气动噪声控制

大涡模拟(LES)和三维的Ffowcs Williams-Hawkings声学比拟方法相结合,研究空腔过流的一种噪声控制措施.空腔的底板/后墙使用多孔壁板,因此流体可以穿透空腔壁面,多孔效果使用Darcy压力-速度关系模拟.声源流场由LES计算,声辐射和远声场由声学比拟获得.结果表明,这种措施有效地减弱了空腔内的压力脉动和远场声辐射,低频脉动Rossiter模数对应的波动幅值被有效抑制,声源中偶极子占优项大幅度减小,从而抑制了声辐射.     

圆柱近尾迹湍流涡结构的直接数值模拟

本文运用谱元方法对圆柱尾流进行直接数值模拟。通过对中等Re数（1000）时的圆柱近尾迹湍流涡结构的演变系统细致的研究,经过长时间流动计算,用速度矢量及相应涡结构示意图显示主涡,二次涡结构的相互作用,显示整个尾迹涡结构在一个周期内的演变过程,获得了在实验中难以清楚观察到的近壁流动特征。本文得到圆柱近壁处二次涡形成有两种不同的机理,并总结出涡的合并所需要的条件。     

椭圆涡环相互作用的离散涡数值模拟

本文发展了三维离散涡模型研究了单个椭圆涡环以及两个沿同一轴线运动涡环的相互作用.首先用两个相向运动圆涡环的碰撞过程与实验的对比验证了模型的有效性,随后研究了单个椭圆涡环的运动,得到了其长短轴交替互换周期的表达式.最后以此为基础分别研究了沿同一轴线两个椭圆涡环的相向碰撞和交替穿越过程,发现碰撞后生成的新涡环变大而形状与涡环的初始形状相似,而在两个椭圆涡环的相随运动中其长短轴交替互换则始终贯穿于交替穿越过程.     

涡发生器对高超声速轴对称进气道外部流动的影响

采用三维CFD黏性模拟考察涡发生器对高超声速轴对称进气道外部流动的影响.针对前缘钝化半径0.8 mm和3.2 mm的轴对称进气道外部流场,以涡发生器高度与当地位移边界层厚度比值为影响参数,考察流场结构与性能参数的影响规律.结果表明,涡发生器产生的干扰波系使得前缘激波向外偏移,下游近壁面流动与主流区出现明显的交换,下游流动出现明显的展向非均匀性.涡发生器对流动的影响沿流向逐渐减弱.在气流压缩性能方面,涡发生器下游压比、动压比沿流向开始增大,随后逐渐恢复到无涡发生器工况;Mach数、总压恢复系数开始降低,随后逐渐向无涡发生器工况趋近.涡发生器高度与当地位移边界层厚度的比值h可作为衡量其影响的重要参数.当h≤1.5时,进气道流场结构、性能参数的变化几乎可忽略,h≤3.0时进气道入口处性能参数几乎能够恢复到无涡发生器工况.      

水轮机导叶分离流动的数值模拟

1引言分离流在本质上是粘性流动和非粘性流动互相干扰的一类复杂流动,绕流中的分离现象往往使水力机械叶片水力损失急剧增加、效率下降,形成的压力脉动还将造成水力机械的振动山。尾水管涡带的形成和分离密切相关,是造成水力机械在部分负荷时机组振动和出力摆动的主要根源［‘j,在偏离最优工况时,叶片内的分离流动将十分严重。本文结合某电站机组出现的振动问题,求解Navier－Stokes的二维方程组,对导叶分离流场进行了分析。2控制微分方程组设导叶内流动为二维不可压粘性流动,控制微分方程组的守恒形式为：。、。分别为。,y两方向…     

内破裂不稳定性及其数值模拟

用磁重联模型与Ｊａｈｎｓ等人提出的锯齿触发条件模拟托卡马克的锯齿振荡现象,对Ｔ－１０托卡马克上的一特定放电,将数值结果与实验作了比较。结果发现,在反转半径外的振幅太大,并且在破裂之前温度剖面太平坦,而中心温度较低。这一结果与重联模型中要求中心安全因子接近１有关。     

混合层流动的格子涡方法数值模拟

采用格子涡方法对二维平板混合层流动进行了数值模拟,重点考察了格子涡方法定量模拟混合层流动流场下游不同流向位置各湍流统计平均量的能力。模拟结果表明:除横向脉动速度均方根比实验值大外,混合层流动在下游不同流向位置的流向平均速度,流向脉动速度均方根和雷诺应力的模拟结果都与实验结果定量符合得很好。此外,模拟结果还很好地反应出混合层流动的自相似性质。从而表明,格子涡方法具有对空间发展混合层流动进行较精确的定量模拟的能力。     

圆柱和直桨叶突然启动瞬态流动的数值研究

采用基于动网格方法的有限体积法对圆柱和直桨叶突然启动引起的二维非定常不可压粘性流进行了数值模拟,给出了计算方法.通过计算得到了Re=5000和9500时圆柱突然启动后的流场随时间演化的非定常过程.将计算结果与实验结果进行了对比,两者吻合较好;采用大涡模拟对高雷诺数下的直桨叶突然旋转启动的流场进行了计算,得到了完成启动后叶片内部瞬态流场的分布结果.     

旋流燃烧室内同向和反向旋转射流湍流流动的数值模拟

针对发展高效低污染旋流燃烧技术的需要,对旋流燃烧室内两股同轴旋转射流相互作用的湍流旋流流动进行了数值模拟.计算中采用了一种新的代数Reynolds应力模型和QUICK离散格式.在两股射流同向旋转和反向旋转两种条件下,将模拟得到的燃烧室内湍流旋流流动的时均气体轴向速度场、切向速度场和静压场与实验数据进行了比较.     

三角翼涡破裂的高精度数值模拟

采用5阶精度的加权紧致非线性格式(WCNS-E-5)数值模拟65°后掠角尖前缘三角翼的大攻角跨声速绕流流场,考察低耗散、高分辨率的WCNS-E-5格式对于三角翼涡破裂模拟的适用性,及激波旋涡干扰对涡破裂点位置的影响,重点研究三角翼大攻角旋涡破裂点的突然前移.通过求解任意坐标系下的非定常雷诺平均N-S方程,采用WCNS-E-5和SST两方程湍流模型,与试验结果和文献计算结果对比,表明既有高阶精度又能光滑捕捉激波的WCNS格式在模拟三角翼旋涡破裂方面具有一定优势,其数值结果与试验结果吻合较好,三角翼大攻角旋涡破裂点的突然前移是由于跨声速流场的激波旋涡干扰.     

介质表面高功率微波击穿的数值模拟

研究了用于模拟高功率微波条件下介质表面击穿的静电PIC-MCC模型,并通过自行编写的数值模拟程序模拟了真空及不同气压条件下介质表面击穿过程中的次级电子倍增和气体电离等过程.模拟结果发现,在真空及低气压条件下,电子的主要来源是次级电子倍增,电子数量以两倍于入射场的频率振荡;在高气压情况下,电子的主要来源是气体电离. 关键词： 介质表面击穿 高功率微波 数值模拟 次级电子倍增     

环形磁场金属等离子体源冷却流场的数值模拟与优化

环形磁场金属等离子体源作为一种全新的等离子体源结构,可用于产生高度离化、无大颗粒、高密度的离子束流,但传统流道结构不能保证其高效、均匀散热,大功率工作时可能引起密封胶圈的烧蚀失效,需对其冷却流场进行优化设计.利用Solidworks Flow Simulation软件对等离子体源冷却流道进行模拟,分析出入水孔分布角度、...     

介质表面高功率微波击穿中释气现象的数值模拟研究

建立了一个简单的高功率微波(HPM)介质表面击穿释气模型,并采用PIC(partiele-in-cell)-MCC(Monte Carlo collisions)方法,通过自行编写的介质表面击穿数值模拟程序对不同释气条件下的介质表面HPM击穿过程进行了数值模拟研究,得到了击穿过程中电子数量等的时间图像和不同释气速度下的击穿延迟时间.模拟结果表明,对于具有一定时间宽度的HPM脉冲,当介质表面气体脱附速度较小时,由于介质表面气体层形成太慢而不会发生击穿;只有当脱附速度大于一定值时,击穿才会发生且击穿延迟时间在一定范围内随着脱附速度的增加而缩短.最后,将数值模拟得到的介质表面HPM击穿数据,与单极性表面击穿的实验诊断图像进行了对比,两者的发展趋势符合很好. 关键词： 释气现象 介质表面击穿 高功率微波 数值模拟     

Counter rotating vortex pair structure in a reacting jet in crossflow

This paper analyzes the time averaged flow structure of a reacting jet in cross flow (RJICF), emphasizing the structure of the counter-rotating vortex pair (CVP) by using simultaneous tomographic particle image velocimetry (TPIV) and hydroxyl radical planar laser induced fluorescence (OH-PLIF). It was performed to determine the extent to which heat release, and the associated effects of gas expansion and baroclinic vorticity production, impact the structure of the CVP. These results show the clear presence of a CVP in the time averaged flow field, whose trajectory lies below the jet centerline on either side of the centerplane. Consistent with other measurements of high momentum flux ratio JICF in nonreacting flows, there is significant asymmetry in strength of the two vortex cores. The strength and structure of the CVP was quantified with vorticity and swirling strength (λ_ci), showing that some regions of the flow with high shear are not necessarily accompanied by large scale bulk flow rotation and vice-versa. The OH PLIF measurement allows for correlation of the flame position with the dominant vortical structures, showing that the leeward flame branch lies slightly above, as well as, in the region between the CVP cores.     

激光辐照下旋转柱壳温度场的数值模拟

采用有限元方法数值模拟在连续激光辐照下旋转柱壳温度场的变化和分布情况,并分析了热性能参数对温度场造成的影响,同时还比较分析了不同旋转频率对柱壳温度场分布的影响。结果表明,激光作用下旋转柱壳的温升大大低于静止柱壳的温升, 外表面温度呈现出与旋转频率相符的周期性上升过程,而内表面温升由于热传导的原因在较小频率下才表现出这种周期性,当频率增大到一定值时,内表面温升不出现周期性的台阶而是曲线上升。     

Evolution of a single vortex in a uniform nonequilibrium medium

The problem of the evolution of a single vortex in a nonequilibrium medium with a temperature-or density-dependent energy release was solved. In a medium that is stable against small perturbations, the vortex passes to a new stationary state, and no vortex breakdown occurs. For a temperature-dependent energy release, analytical solutions that define changes in the characteristics of the vortex were obtained for high and low relaxation times. The solutions were compared with the results of numerical simulation by the Godunov method.     

用晶格玻尔兹曼方法研究颗粒在涡流中的运动

用晶格玻尔兹曼方法研究小颗粒在涡流中的运动.涡流由流经空腔的流体产生.用动量交换法和压力张量积分法计算颗粒在涡流中的运动轨迹、速度和角速度.最后用张量积分法计算两个不同半径的颗粒在涡流中的运动. 关键词： 晶格玻尔兹曼方法 涡流 颗粒     

激光辐照下充压柱壳的破坏能量阈值数值模拟

采用有限元方法模拟了激光辐照下充压柱壳的热力学响应,计算了不同工况下结构的瞬态温度场和应力场,根据材料强度准则判断了柱壳的破坏时刻,并提出了一种预测激光辐照下充压柱壳破坏能量阈值的数值方法。研究了壳体厚度和内压大小对柱壳破坏能量阈值的影响,并给出了典型工况下柱壳破坏能量阈值同壳体厚度以及充压大小的关系。数值计算结果表明:破坏能量密度阈值与壳体厚度、内压大小近似成线性关系,壳体厚度比内压大小对柱壳的激光破坏能量阈值影响更大。     

亚纳秒气体开关中气体击穿的数值计算

对亚纳秒气体开关的放电过程,提出了气体击穿阶段的物理与数学模型以及数值计算方法,并对氢气和氮气的这一击穿过程进行了数值仿真。对于1mm氮气间隙,计算了充电电压波形为纳秒级上升的斜角平顶波时,开关放电时延以及击穿电压随气压和充电电压上升时间的变化。对于氢气,利用美国Sandia实验室的一个实验结果对数值仿真方法进行了验证,所得到的计算结果与实验结果很好地吻合,初步表明所建立的物理模型与计算方法的正确性。     

激光辐照下充压柱壳的破坏能量阈值数值模拟

采用有限元方法模拟了激光辐照下充压柱壳的热力学响应,计算了不同工况下结构的瞬态温度场和应力场,根据材料强度准则判断了柱壳的破坏时刻,并提出了一种预测激光辐照下充压柱壳破坏能量阈值的数值方法。研究了壳体厚度和内压大小对柱壳破坏能量阈值的影响,并给出了典型工况下柱壳破坏能量阈值同壳体厚度以及充压大小的关系。数值计算结果表明：破坏能量密度阈值与壳体厚度、内压大小近似成线性关系,壳体厚度比内压大小对柱壳的激光破坏能量阈值影响更大。