透平级S_2流面内的反问题湍流计算

本文发展了应用湍流流函数方程计算湍流流场的方法。计算了扭曲叶片Ｓ２流面的湍流流场．在非正交曲线坐标可压湍流方程组中包含了以无量纲形式表示的能量方程，熵方程和Ｋ－ε湍流模型方程。固壁边界上采用了壁函数方程。通过计算得到了透平级内部的湍流流动合理的收敛解。此方法可以通过调节子午通道的形状及环量的大小分布来防止流动中负速度的出现，为高效透平级的设计提供了一种新的途径。     

三维可压缩圆孔单相射流近场DNS并行数值模拟研究

本文对圆孔射流近场的湍流统计量进行了DNS模拟研究.模拟的Re=5000,计算网格3000万,计算CPU个数为60.采用MPI并行求解了可压缩的射流守恒方程,包括连续性方程、动量方程、能量方程和理想气体的状态方程.对射流近场(30d以内,d为喷口的直径)的湍流统计量的迁移规律进行了系统的分析.研究了入口湍流条件对射流统...     

k－ε－PDF两相湍流模型和台阶后方气粒两相流动的模拟

本文提出了两相湍流的ｋ－ε－ＰＤＦ模型。ＰＤＦ模型所得的湍流两相流动的统计平均方程，与雷诺时均方程有相似的形式，但ＰＤＦ模型可精确计算出颗粒相各方程的脉动关联项。本文将该模型用于预报台阶后方湍流两相流动，与ｋ－ε－ｋｐ模型相比，它可以更合理地预报出颗粒湍流的各向异性。     

用视密度加权平均二阶矩模型模拟旋流两相流动

本文用视密度加权平均代替时平均,建立了视密度加权平均的统一二阶矩两相湍流模型方程组（ＭＵＳＭ）,其中用体积分数代替了数密度,用颗粒驰豫时间作为封闭两相脉动速度关联方程耗散项的时间尺度,并引入了颗粒视在的气体速度脉动的输运方程。用ＭＵＳＭ模型模拟了旋流数为０．４７的气粒两相流动。并和实验结果及时间平均的ＵＳＭ模型的模拟结果进行了对照,两种模型均能较好地预报的两相的轴向和切向速度,轴向和切向脉动速度。此外,ＭＵＳＭ模型可以减少所用方程数,节省计算量。因此视密度加权平均的统一二阶矩两相湍流模型是一种对时间平均的统一二阶矩模型的改进,今后可以进一步扩大应用。     

开敞式蜗壳轴流式水轮机湍流计算及性能预估

本文在进行模型水轮机试验的基础上,完成了开敞式蜗壳轴流式模型水轮机的全流道三维定常湍流计算。计算程序以雷诺时均连续方程和N-S方程作为控制方程,选用标准κ-ε双方程湍流模型使方程组封闭。数值求解通过CIMPLEC算法实现速度、压力的分离求解,离散差分格式具有二阶精度。通过计算获得了各过流部件的流动细节,预估了水轮机的能量性能,并与试验结果进行了比较。     

DSM-LPDF两相湍流模型及旋流两相流动的模拟

本文由流体－颗粒速度的拉氏联合概率密度函数（ＰＤＦ）输运方程出发,用Ｓｉｍｏｎｉｎ建议的Ｌａｎｇｅｖｉｎ模型封闭颗粒所遇到流体瞬时速度的条件期望项,并用Ｍｏｎｔｅ Ｃａｒｌｏ方法直接求解 ＰＤＦ输运方程,将其和求解流体雷诺应力方程模型的有限差分方法结合,建立了雷诺应力－拉氏ＰＤＦ（ＤＳＭ－ＬＰＤＦ,简称ＤＬ）两相湍流模型．用此模型模拟了旋流数为０．４７的突扩旋流气粒两相流动,并与文献中ＰＤＰＡ实验和用类似于单相流动湍流模型封闭方法的时平均统一二阶矩（ＵＳＭ）模型的预报进行了对比．     

基于Lagrange方法封闭的两相湍流场方程模型

两相湍流场方程模型采用基于Euler方法的一阶矩方程,而二阶矩方程由Lagrange方法得到,新模型的封闭不需要附加其它假设.首先基于概率密度函数给出颗粒运动的连续方程和动量方程,其次由基于平均Langevin方程的Lagrange模型推导得到颗粒二阶矩方程,最终获得封闭的二阶矩模型.将新模型用于气固两相壁面射流的数值模拟,结果表明新模型合理有效.     

一种新的k方程转捩湍流模型

本文在分析了针对充分发展湍流模型建立的k方程湍流模型模拟转捩流动的能力之后，构筑了一种新的Fμ函数及相应的k方程湍流模型．将该模型应用于数值求解NS方程，模拟了中等来流湍流度下零压力梯度平板捩报流动，并将部分计算结果与实验结果及经验关系进行了比较．     

湍流模型在复杂流场数值模拟中的应用

采用4种湍流模型:代数Baldwin Lomax(B-L)模型、半方程Johnson King(J-K)模型的两个版本(J-K90A和J K92)以及两方程k-g模型,分别数值模拟了导弹超音速流动、NASATND-712标模和民机翼身组合体(两区C-O网格)跨音速流动.采用中心有限体积和多步Runge-Kutta方法数值积分三维可压缩雷诺平均Navier-Stokes(N-S)方程组.k-g湍流模型方程的求解采用类似于N-S方程组的方法进行.所有湍流模型均能很好地模拟附体及小分离流动;对于大攻角、分离剧烈的导弹流动,k-g和J-K92模型与实验吻合更好;B-L模型在模拟民机跨音速流动时,它所捕捉的激波位置较其余3种模型靠后.利用多块网格模拟民机翼身组合体流场时,k-g模型的模拟能力强于其余3种模型.     

用求解概率密度函数输运方程的方法模拟湍流自由射流

用求解概率密度函数输运方程的方法模拟湍流自由射流周向阳，郑楚光，马毓义（华中理工大学煤燃烧国家重点实验室武汉４３００７４）关键词概率密度函数；湍流；自由射流１前言和通常的湍流统计矩模型（如ｋ－ε模型或Ｒｅｙｎｏｌｄｓ应力模型等）相比，用求解速度和标量...     

Arterial stenosis murmurs: an analysis of flow and pressure fields

The flow field distal to an arterial stenosis is simulated by a confined turbulent jet with moderate Reynolds numbers. The wall pressure fluctuations are related to the momentum fluctuations of the jet by the Poisson equation. A Green's function was derived to satisfy the boundary conditions on a cylindrical surface. This allows the solution of the Poisson's equation to include only a volume integral of the fluctuating momentum, weighed by the relative distance between the source and the sensor. The velocity fluctuations on the jet centerline and at the middle of the shear layer were measured using a laser Doppler anemometer. The wall pressure fluctuations were detected simultaneously by an array of nine wall-mounted pressure transducers along the axial direction. Cross correlation performed between the velocity and pressure fluctuations reveals that the pressure fluctuations were mostly imposed by the passage of turbulent eddies with a convective velocity that is a function of the jet exit velocity. The Strouhal number, defined by the frequency of the passing large-scale structure, is a function of the initial conditions only very close to the jet exit. Further downstream, where the effect of the initial conditions is lost, the Strouhal number approaches a constant irrespect of the jet Reynolds number. The contribution of a source near the jet exit to wall pressure fluctuation near the reattachment is rather weak due to the rapidly decaying weighting function in the axial direction. However, for sources located within one nozzle diameter from the sensor, the cross-spectral density function has a high magnitude with maximum coherence where the pressure spectral changes its slope.     

湍流理论的近代发展

这篇论文的内容是关於不可压缩流体的湍流理论近代发展的综合性介绍舆分析。我们首先评述了根据Reynolds的平均运动方程所建立的混合长度理论。其次,分析关於均匀各向同性湍流的主要理论工作。第三,讨论了运用Reynolds的平均运动方程和根据速度涨落方程求得的速度关联函数的动力学方程来处理具有Reynolds剪应力的普通湍流运动问题。同时说明这个方法虽然能够给出比混合长度理论舆实验较为接近的理论结果并能提出速度涨落平方平均值的理论分布,但是由於在求出的速度关联的动力学方程中出现高次元的速度关联,它继续地导致不封闭的微分方程组因而遇到不易克服的困难。因此,从以上湍流理论发展的回顾和最近关於均匀各向同性湍流在后期衰变运动的涡性结构工作,我们在最后提出了对今后湍流理论研究工作的新看法:湍流运动的基本组成部分是流体粘性作用所引起的涡旋运动;这个涡旋运动的动力学根据是用平均的方法后Navier-Stokes方程所导出的Reyonlds的平均运动方程典带度涨落方程。我们并着重说明Reynolds认识到湍流运动可分作平均运动与涨落运动的重要性。今后的理论工作则在於求这两组动力学方程的涡旋运动解,而这种类型的解并须满足像Колмогоров在高Reynolds数运动的局部各向同性湍流理论中所提出的统计条件,方能使解满足惟一性并可舆实验结果相比较。     

Equations for describing a family of self-similar solutions of the far field of a round submerged jet

The flow of a turbulent round submerged jet in far field has been considered. The modified Navier?Stokes equations have been derived with the aid of known self-similar properties of the time-mean flow as well as the spatial and temporal characteristics of the entire spectrum of the scales of turbulent fluctuations. A technique has been proposed for the numerical computation of these equations to describe a family of self-similar solutions in a cylindrical region with periodic boundary conditions in the longitudinal direction, and the velocity of the jet expansion enters the obtained equations explicitly as a parameter.     

Jets in tokamaks: A theoretical study

A new class of bifurcation of the momentum balance equations for a tokamak plasma is presented. The solution exhibits a monopolar localized jet of ExB flow. The jet is generated by the reduction of turbulent viscosity due to ExB velocity shear. Strong jets of localized plasma flow have been observed in tokamaks as a precursor to the development of a transport barrier region with reduced turbulent transport. The jet solution is shown to fit well with the experimental observations.     

Flame-acoustic coupling of combustion instability in a non-premixed backward-facing step combustor: the role of acoustic-Reynolds stress

Combustion instability in a laboratory scale backward-facing step combustor is numerically investigated by carrying out an acoustically coupled incompressible large eddy simulation of turbulent reacting flow for various Reynolds numbers with fuel injection at the step. The problem is mathematically formulated as a decomposition of the full compressible Navier–Stokes equations using multi-scale analysis by recognising the small length scale and large time scale of the flow field relative to a longitudinal mode acoustic field for low mean Mach numbers. The equations are decomposed into those for an incompressible flow with temperature-dependent density to zeroth order and linearised Euler equations for acoustics as a first order compressibility correction. Explicit coupling terms between the two equation sets are identified to be the flow dilatation as a source of acoustic energy and the acoustic Reynolds stress (ARS) as a source of flow momentum. The numerical simulations are able to capture the experimentally observed flow–acoustic lock-on that signifies the onset of combustion instability, marked by a shift in the dominant frequency from an acoustic to a hydrodynamic mode and accompanied by a nonlinear variation of pressure amplitude. Attention is devoted to flow conditions at two Reynolds numbers before and after lock-on to show that, after lock-on, the ARS causes large-scale vortical rollup resulting in the evolution of a compact flame. As compared to acoustically uncoupled simulations at these Reynolds numbers that show an elongated flame with no significant roll up and disturbance in the upstream flow field, the ARS is seen to alter the shear layer dynamics by affecting the flow field upstream of the step as well, when acoustically coupled.     

Investigation of a turbulent wall jet in forced convection issuing into a directed coflow stream

In the present paper, we have studied numerically the directed coflow stream effects on mean and turbulent flow properties of a turbulent plane wall jet in forced convection emerging into a directed coflow stream. The system of equations governing the studied configuration is solved with a finite difference scheme using a staggered grid for numerical stability, not uniform in the two directions of the flow. The modified version of the first-order low Reynolds number k–? turbulence model is used and compared to existing experimental findings. It is found that predicted results are in satisfactory agreement with the experimental data and that the wall jet fluid decays faster in presence of a directed coflow stream. Results show also that the increase of coflow deviation angles causes an increase of the growth rates of the dynamic and thermal half-width of the jet and enhances the turbulent mixing. It is found that the longitudinal development of normalised forms of the jet characteristics parameters at different directed coflow velocity ratios can be reasonably well collapsed onto universal trends through the use of momentum length scale.     

托卡马克边缘湍流中带状流的负粘滞阻尼

基于Hasegawa-Wakatani湍流模型,数值模拟了托卡马克边缘等离子体中漂移波湍流和相关的反常粒子输运。从等离子体动量守恒方程出发导出了不采用常规的布辛涅斯克近似的带状流方程,论证了大振幅密度扰动和湍性粒子流对激发带状流的贡献可等效地对应于低阶负粘滞阻尼效果。数值模拟表明,大振幅密度扰动的非线性大大增强了带状流饱和振幅,从而有效抑制了湍性粒子输运。研究结果阐明了托卡马克边缘等离子体大振幅密度扰动的非线性对驱动等离子体旋转、动量输运和带状流的重要性。     

�п���˱�Ե�����д�״���ĸ�ճ������

基于Hasegawa-Wakatani湍流模型,数值模拟了托卡马克边缘等离子体中漂移波湍流和相关的反常粒子输运.从等离子体动量守恒方程出发导出了不采用常规的布辛涅斯克近似的带状流方程,论证了大振幅密度扰动和湍性粒子流对激发带状流的贡献可等效地对应于低阶负粘滞阻尼效果.数值模拟表明,大振幅密度扰动的非线性大大增强了带状流饱和振幅,从而有效抑制了湍性粒子输运.研究结果阐明了托卡马克边缘等离子体大振幅密度扰动的非线性对驱动等离子体旋转、动量输运和带状流的重要性.     

Zone conditional analysis of a freely propagating one-dimensional turbulent premixed flame

The zone conditional conservation equations are derived and validated against the DNS data of a freely propagating one-dimensional turbulent premixed flame. Conditional flow velocities are calculated by the conditional continuity and momentum equations, and a modeled transport equation for the Reynolds average reaction progress variable. An asymptotic formula for turbulent burning velocity is obtained with the effects of a finite Damköhler number accounted for as an additional factor. It is shown that flame generated turbulence is primarily due to correlations between fluctuating gas velocities and fluctuating unit normal vector on a flame surface. More investigation is required to validate general predictive capability of the derived conditional conservation equations and the relationships modeled for closure.