逆压梯度湍流边界层的新涡粘模式

文献[1]提出了着重体现湍流应力各向异性的湍流新涡粘模式,并提出应用新涡粘模式预报湍流平均流动的途径:由每类剪切流比较基本流动已知的实验结果来寻找涡粘张量分量分布的规律,再去预报较复杂的剪切流动.本文对二维湍流边界层,按此途径找到了涡粘张量分量的分布规律,并准确、简便地预报了逆压梯度二维湍流边界层的平均流动,与实验符合得相当好.1方程和边界条件     

湍流运动的非随机场封闭数学模型

本文把湍流场定义为存在离散涡的流场。并定义涡旋特征尺度λ→∞时的流动状态为层流或湍流场的基态。由于湍流场存在拟序结     

激波与湍流相互作用的实验研究

在超声速湍流混合层风洞中, 以NPLS技术为基础研究了激波与湍流的相互作用, 观察到了斜激波与湍流边界层的相互作用以及湍流混合层中的大涡对斜激波的影响. 所得到的NPLS实验图像不仅再现了复杂的流场结构, 而且不同图像之间还具有时间相关性. 在NPLS图像基础上探讨了激波与湍流相互作用的规律.     

不可压湍流燃烧大涡模拟的分步投影方法

通过对Ma<<1的湍流燃烧过程物理现象的分析, 从数学和物理上讨论了低Ma数湍流燃烧的“不可压”性质, 在此基础上采用改进的不可压湍流燃烧数值模拟的分步投影方法对甲烷-空气平面射流扩散火焰进行了大涡模拟. 该方法不需对密度进行预估而是直接显式求得, 仅需对动量方程进行一次投影即可实现对N-S方程的显式求解, 计算量小, 求解效率高, 在Ma<<1的变密度有反应或无反应湍流流动的细观数值模拟中具有明显优势. 甲烷-空气的化学反应采用简化的四步反应机理模拟, 应用动态亚格子模式及动态相似模式分别模拟湍流流动的亚格子项及化学反应速率的亚格子项. 计算结果得到了射流扩散火焰中瞬态发展变化的涡运动以及涡-火焰相互作用的细节结构.     

关于地表湍流热通量计算中函数φ_p(x)的一个实验研究

一、问题的提出地表湍流热通量对于物理气象和应用气象都是重要的。实际工作中一般都根据贴地气层的温度及风速梯度资料简接计算。以热量湍流系数K_T和动量湍流系数K相等为条件的“扩散法”曾被广泛应用。随着贴地层物理概念的发展和实验资料的积累,人们发现这两个湍流系     

近壁圆柱绕流尾流的特性

近壁圆柱绕流是许多工程问题分析简化的物理模型。圆柱尾流与壁面边界层的相互作用产生了复杂而丰富的多尺度湍流结构。文章对工程中由这种湍流结构引起的圆柱涡脱落频率、圆柱所受阻力和升力的变化，以及这种湍流结构的产生和演化的规律进行了综述和分析，并对进一步开展的研究提出了建议。     

受损伴流湍流氢气射流火焰的直接数值模拟

采用直接数值模拟DNS的方法对受损伴流湍流氢气射流火焰进行了数值模拟, 采用16步的氢气燃烧详细化学反应机理, 冷的高速H₂/N₂燃料射流喷入热的低速伴随流, 伴随流由贫燃氢气预混火焰燃烧产生, 温度1045 K, 氧量较低. 化学反应源项由主程序在每一时间步长内动态调用CHEMKIN库函数获得. 计算采用消息传递MPI的并行计算方法, 采用12颗CPU在并行计算机上完成. 作为与实验对比的Faver平均结果由DNS瞬态结果做长时间的统计平均后获得. 火焰中涡结构的卷起以及发展过程均能很好地被捕捉, 可以观察到同旋向涡结构之间的相互吸引和反旋向涡结构之间的相互排斥过程, 伴随射流两侧涡结构彼此复杂的吸引、合并、挤压和撕裂过程, 湍流拟序结构由最初的轴对称模式开始向非对称模式演化. 流场中5.67 ms时刻瞬态的H, OH和H₂O分布, 表征了燃料射流自点燃过程中的详细火焰结构. 计算中获取的火焰抬升高度为9d ~ 11d, 与实验结果相吻合. 计算发现由OH和H粒子表征的火焰锋面中, 在火焰锋面转角位置, 燃烧过程得到强化, 可能与火焰面的拉伸以及较长的停留时间有关. 从湍流强度的分布曲线来看, 火焰的传播应该是从两侧向中心发展的. 这里的DNS结果可以作为今后发展更准确通用湍流燃烧模型的参考.     

基于NPLS的超声速层流/湍流后台阶流动精细结构研究

在Ma=3.0低噪声、吸气式超声速风洞中,对台阶高度h=5mm的超声速后台阶流场进行了精细结构测量.通过改变台阶上游壁面的表面粗糙度,实现了超声速层流、湍流两种后台阶流动.采用NPLS技术对流场整体结构的时空演化特性以及4个局部典型区域的细节结构等方面进行了实验研究.瞬态流场揭示了扇形膨胀波系、再附激波、超声速边界层及其分离、再附和恢复等结构的空间特征.通过比较时间平均的结果,可知超声速湍流后台阶流动分离后的膨胀角较大、回流区的长度相对更短,而再附后重新发展的边界层厚度以较小的倾角增长,但两种流动的再附激波角度大致相同.在时间演化上,超声速层流后台阶流动主要表现为K-H涡结构的变形受剪切、膨胀、再附以及三维效应等影响;而湍流后台阶流动则主要表现为大尺度结构在再附点前后受膨胀、黏性以及再附后逆压梯度的作用而倾斜和变形.对局部区域的研究表明,在超声速层流后台阶流场中微弱压缩波与当地对流马赫数和K-H涡结构的诱导作用有关,并且在下游汇聚成再附激波的现象明显;而湍流后台阶流场中则未有明显的压缩波和K-H涡结构,其再附激波的形成主要与壁面的压缩效应有关.     

竖直平板间湍流自然对流中的螺旋羽流结构

通过直接数值模拟研究竖直平板间温差湍流自然对流的相干结构。研究表明，流场中同时存在大尺度展向涡结构和螺旋羽流结构。螺旋羽流结构对应于流场中相对螺度、法向涡量及脉动温度较大的位置，定性分析了 螺旋羽流结构的形状、特点、成因以及与Rayleigh-Benard对流的异同，并应用条件采样分析说明了螺旋状热羽流结构的主要特性。     

中性大气边界层中风力机的湍流演化及叶根载荷分析

针对一台外场试验风电机组,采用大涡模拟(large eddy simulation, LES)耦合致动线模型的方法,构建了中性大气边界层和风力机风轮的气动耦合求解模型,模拟风力机在中性大气边界层中的流场.通过连续小波分析、频谱分析和相关性分析,研究了中性大气边界层中风力机前、后的湍流演化过程及其与叶根载荷的相关性.研究发现,自然来流从风轮前1D(D为风轮直径)处运动到后1D处时,大气中的湍流强度逐渐增大;在风轮平面处出现了较强的小尺度湍流结构,这些小尺度的湍流结构在向下游运动过程中不断耗散,并在风轮后1D处能量基本耗散殆尽;叶尖位置处的高频湍流出现频率约为1.82 Hz,此频率正好与叶片通过频率相对应.风力机的叶根挥舞载荷对大气中的湍流结构响应明显,低频湍流结构对叶根挥舞载荷的低频段影响显著,高频湍流结构对叶根挥舞载荷的高频段影响明显;叶尖高频湍流结构相对于叶根高频湍流结构,频率更高,能量更大,其对叶根挥舞载荷高频段的影响更为明显;同时,叶尖高频湍流与叶根挥舞载荷的高频部分表现出了一致的周期性变化规律.     

湍流提升火焰条件矩封闭模拟

采用k-ε-g湍流模型、23步反应机理和径向平均条件矩封闭(Conditional Moment Closure)CMC技术, 模拟了氢气湍流提升火焰结构, 并探讨了火焰稳定机理. 作为针对湍流提升火焰进行耦合CMC模型的多维大涡模拟LES的第一步, 计算考虑了湍流对燃烧的作用, 数值格式采用高精度Padé格式, 推导并采用了基于特征波分析的多组分条件矩Navier-Stokes特征边界条件NSCBC, 本研究暂不考虑热释放对湍流的影响, 模拟集中于火焰提升区域. 与国际上公布的相关实验数据的比较表明, 模型精确地预测了提升高度、Favre径向平均温度和组分浓度. 研究结果表明, 径向平均CMC模型能够较好重构目前仍争议较大的提升火焰稳定机理.     

气固两相流动中颗粒扩散的转捩现象

不引入任何湍流模型, 采用有限容积方法对三维气固两相湍流射流进行了直接数值模拟, 着重考察湍流结构的转捩行为对不同Stokes数颗粒扩散的影响. 为了得到高精度的结果, 同时降低计算量和存储量, 对气相流体控制方程组的求解采用分步投影算法, 对时间积分采用低存储、三阶精度的Runge-Kutta积分格式; 对颗粒控制方程的求解在拉格朗日框架中进行. 流场的统计结果与相关的实验数据吻合良好, 证实了数值算法的可靠性. 对颗粒扩散的研究发现, 展向涡结构对颗粒扩散的影响比较突出. 而在流场的演化过程中, 观察到颗粒扩散的新行为, 即中、小Stokes数的颗粒在流场中的分布出现了 “转捩”现象.     

近地面稳定大气层结中通量-廓线模式的适用性

在近代湍流理论研究和实际计算中,苏联学者和提出的相似理论方法获得了相当广泛的应用。正如yaglom在文献[1]中所评述那样:在处理大气表面层实验资料时,Monin-Obukhov相似理论方法现在起着主导作用。随着大气表面层观测技术的发展,特别是大气湍流通量直接测定方法的出现(如涡度相关法),不同大气稳定度条件     

多场协同原理在管内对流强化传热性能评价中的应用

针对管内层流和湍流换热,分析了流体质点矢量物理量的协同性,揭示了对流换热多场协同规律与强化传热机理之间的关系,提出了评价强化传热综合性能的效能评价系数EEC,同时,基于协同角α,β,θ,γ,η,初建了判断层流和湍流强化传热性能的统一评价体系以及相应的评价指标.此外,针对内插三角杆强化传热管建立了数学模型,通过数值模拟对所建立的评价体系及其指标进行了验证.计算结果表明,基于多场协同原理的协同角与强化传热效果的评价准则之间相互对应,在Re数为300~1800的范围内,性能评价系数EEC值在1.3~2.3之间,而效能评价系数EEC值则在0.33~0.45之间.     

卫星通讯中信号幅度闪烁现象的统计规律

近年来的研究表明,对流层中出现的湍流现象将造成空气折射系数在局部范围内的不均匀,而由这种不均匀所产生的电磁波前向散射将是Ku波段卫星通讯和广播中信号幅度闪烁现象的主要原因。闪烁的强度与卫星地面接收站的仰角、天线大小、工作频率及气象条件有关。由闪烁现象形成的信号起伏实际上是很多持续时间很短的“衰减过程”。在闪烁足够强     

风应力系数与波高熵关系的研究初探

1 风应力关系式 海面风应力系数C_D定义为 其中U_z为海平面上高度z处的风速,U__*是摩擦风速,k是Karman常数,z_0为海面粗糙度。风应力系数C_D一般认为随海面风速增加,而且与海面状况有关,不过CD增加的速率因不同的研究而有很大的变化。Charnock通过量纲分析给出一个z_0的经验公式     

格子玻尔兹曼方法模拟大气边界层中的阵风结构

东亚北方冬春季常会出现寒潮冷锋天气,冷锋过后有明显的下沉气流,且伴随强阵风.观测数据分析表明,阵风结构具有相干性,是春季沙尘暴能够被输送到边界层上部,从而远距离传输的重要条件.格子玻尔兹曼方法(lattice Boltzmann method,LBM)是一种基于Boltzmann分子输运方程的数值流体计算技术,自从出现以来,由于其在微观水平描述运动的特点,一直被用来研究湍流等复杂运动,本文在格子玻尔兹曼方法中引入大涡模式,使其能够应用到大气边界层湍流场,模拟阵风的产生及发展,得到了具有波动和涡旋相结合的相干结构,从而解释了阵风相干结构起沙扬尘的机理.     

天然橡胶的粘-弹性质——Ⅱ.结晶对流动的影响

高分子化合物的本体流动方式有两种:粘流及塑流,在给定条件下可以互相转变。例如我们曾发现,天然橡胶在中温时的流动方式系非牛顿粘流,于中等应力下且服从 Eisenschitz 公式;而当温度降低时即转变为塑流,且系非宾哈塑流。此外,对于这种从粘流转变为塑流的原因则迄未探明。例如李泰圭及 Eyring 在重新处理 Saunder 及 Treloar 对天然橡胶的流动数据时,曾认为这是一种不明原因的“二级转变”。     

湍流的间歇性和多分维

本文利用湍流间歇性的图象介绍单分维(各向同性湍流)和多分维(非各向同性湍流)的概念。使我们将分维概念从自相似和点集的几何性质,扩展到物理场的强弱区域有不同的标度,并将物理场的功率谱和分维、标度律联系起来。将分维和物理图象联系起来,必将丰富其概念,扩大其应用价值。     

能否发展关于湍流动力学和颗粒材料运动学的综合理论?

自然界和工程应用中遇到的流动绝大部分是湍流问题.湍流流动中包含了不同尺度的旋涡,如何准确地定义湍流仍十分困难,湍流动力学的理论体系仍尚未完成.颗粒流是固相颗粒以气相或液相为载体的流动,缺乏普适的本构方程,处理颗粒运动时无法将微观的颗粒尺度与宏观的流场尺度分开.本文介绍了湍流和颗粒流的基本特性、研究的理论和方法,分析了湍流和颗粒流研究所面临的挑战.