湍流研究最近半世纪的一些发展

<正>内容一、绪言(Ⅰ)研究湍流的意义(Ⅱ)国外开展湍流研究的情况二、湍流理论的进展(Ⅰ)先把流体动力学方程组平均的理论(a)1950年以前发展情况简述(1)Reynolds方程和混合长度理论(2)各向同性湍流的统计理论(3)具有剪应力的普通湍流理论(b)最近的剪切湍流的半经验理论(c)最近的湍流统计理论(1)E.Hopf理论(2)R.H.Kraichnan直接相互作用理论(8)Lewis等人的分子运动理论(4)Meecham理论。(5)S.Grossmann重正化群法(8)陈善谟统计动力学重复级串法

剪切水-气界面下的湍流特征

利用氢气泡技术测量了剪切水-气界面下二维瞬时流速的分布,根据从流动显示图像获得的流速剖面样本,考察水面附近湍流场的平均特性,重点分析猝发事件的信号特征。氢气泡方法较为准确地测定了水流中平均流速、湍流脉动强度以及Reynolds切应力的分布,更重要的是,它能将流动显示图像中观察到猝发事件与定量的信号特征直接联系起来.在猝发过程中,水面近区的流向和垂向瞬时速度以及Reynolds切应力出现较大幅度的脉动,它们与低速流体条带的喷射行为对应;大幅度脉动在时间和空间垂直方向上表现出高度的相干性,这是猝发事件的一个显著特征。

高速剪切层线性稳定性空间模式参数分析

用一类高精度对称紧致差分格式数值离散一阶改型三维可压粘性扰动方程,对导出的非线性离散特征值问题采用二阶修正Newton-Raphson边值迭代局部解法,实现了超声速剪切流的线性空间稳定性分析。基于对流Mach数、Reynolds数、速度比、密度比等参数研究,讨论了压缩性效应、粘性效应、超声速扰动快/慢模态等,结果显示超声速粘性剪切流的弱不稳定性同多种影响因素密切相关。

微重力流体力学

目录

前言
第一章 微重力流体科学概论
一、微重力科学与微重力流体科学
1、微重力环境
2、重力的影响
3、微重力流体科学的发展
二、微重力流体力学概述
1、对流
2、扩散及输运现象
3、液滴和气泡动力学
4、多相流过程
5、残余重力效应
6、其他流体力学问题
三、微重力物理化学概述
1、临界现象
2、燃烧
3、分散体悬浮系统
4、晶体生长的物理化学问题
四、微重力流体科学的研究途径
1、微重力研究的一般途径
2、微重力实验手段
参考文献
第二章 基本方程组和流体运动特性
一、引言
二、连续性方程和迁移方程
三、动量方程
1、流体的粘性——Reynolds应力
2、动量守恒定律
3、Navier-Stokes方程
四、能量方程
1、总能量方程
2、动能方程
3、内能方程
4、粘性耗散函数
5、Fourier定律及另外二种形式的能量方程
6、不可压流体的导热方程
五、Newton流体的运动方程组及定解条件
1、基本方程组和适定性
2、定解条件
六、Boussinesq近似及适用范围
七、相似律和无量纲参数
1、利用Buckinghan〓定理导出相似参数
2、微重力流体力学的有关物理量和无量纲参数
参考文献
第三章 毛细现象以及界面的平衡和稳定
一、引言
二、表面张力的物理描述
三、液体射流的表面不稳定
1、基本方程组和基态
2、小扰动的线性化方程
3、本征值方程及其解
四、等温条件下液桥的平衡位型和稳定
1、表面张力作用下的平衡条件
2、毛细稳定性
3、旋转稳定性
第五章 液桥的流体动力学稳定理论
1、基本假设和液桥的平衡条件
2、稳定问题的数学提法
3、液桥的Liapunov稳定理论
4、特殊情形（Ω〓=μ=0）以及纯半波不稳定（n=1，m=1）
5、小扰动方程的变分方程
6、小Weber数和大Reynolds数情形的不稳定发展率
7、液桥微重力实验的结果的分析
8、讨论和结论
参考文献
第四章 对流和扩散
一、Pearson对流
1、自由面不变形时的小扰动分析
2、自由面可变形情形
3、非线性理论
4、多层不混溶液体系统
二、热毛细对流
1、矩形容器中的热毛细对流
2、柱形液桥的热毛细对流
3、半浮区液桥热毛细对流的数值模拟
4、薄层液体的热毛细对流
三、热毛细振荡对流的实验研究
1、液桥内部的温度振荡
2、热毛细对流的表面振荡
3、综合测量
四、热毛细对流的振荡机理
1、热流体波不稳定性
2、表面波不稳定性
3、有限高度液桥的线性不稳定性
4、三维不定常数值模拟
5、重力的影响
6、一种非稳定性理论
7、关于振荡的激发机制
参考文献
第五章 液滴动力学
一、等温液滴动力学
1、球形液滴的振荡
2、不混溶液体中球形液滴的振荡
3、弱非线性理论
4、实验模拟
二、非等温液滴的Marangoni迁移
1、定常线性化理论（小Reynolds数，小Marangoni数）
2、非线性理论
3、实验结果
三、液滴和气泡的相互作用
1、双气泡的轴对称理论
2、多液滴的轴对称理论
四、旋转液滴的演化序列和分叉理论
1、旋转液滴的演化
2、旋转液滴的Thomson-Tait稳定准则
3、长期稳定性和动力稳定性
4、长期稳定性真实性的实验证明
5、结论
参考文献
第六章 微重力材料流体力学
一、晶体生长过程
二、纯扩散过程
1、一维扩散过程
2、二维扩散过程
3、固-液界面弯曲对径向分凝的影响
三、浮区晶体生长
1、浮区的热毛细对流
2、浮区的熔质毛细对流
3、浮区对流的振荡特征（小Prandtl数对流）
4、耦合过程
四、溶液晶体生长
1、溶液晶体生长的相变界面过程
2、一维纯扩散过程
3、准定常溶液晶体生长过程
4、不定常溶液生长过程
五、气相晶体生长
1、气相晶体生长过程
2、一维模型
3、物理气相输运中的对流效应
4、化学气相沉积（CVD）过程
参考文献

郭永怀文集

可压平面混合层稳定性分析及数值模拟

可压平面混合层是包含复杂多时空尺度运动的非定常流体力学部问题，具有深刻的理论意义和广泛的应用背景。针对该问题所涉及内容的多面性，本文的目的是，基于高精度、高分辨率数值算法的构造、发展和数值行为分析，采用线性稳定性分析和直接数值模拟方法。从理论和计算两方面集中研究压缩性效应、粘性效应、初值效应以及燃烧反应放热效应等对可压平面混合层早期稳定性行为和大尺度拟序涡结构非线性演化的影响。以混合层已有研究成果的分析和综述为开端，论文主体共包括四部分：第一部分是可压平面混合层时间/空间模式数值线性稳定性分析。实现了高精度对称紧致差分格式（SCD）对可压粘性扰动线性稳定性边值问题的求解，对导出的线性和非线性离散特征值问题，提出了两个高效局部解法。研究涉及二维/三维扰动波、无粘/粘性扰动波、特征函数和特征值谱、第一/第二模态、超声速快/慢模态、速度比和密度比等。验证了对流Mach数Mc为一个合理的压缩性参数。指出压缩性效应和粘性效应对最不稳定扰动波的波数（频率）和增长率呈相拟的抑制作用，且时间模式稳定性分析结果在许多方面是可信的。从随机和线性扰动场出发，采用高精度五阶迎风紧致和六阶对称紧致混合差分算法（UCD5/SCD6）对可压平面混合层的稳定性特征进行了直接数值模拟，揭示了初始主导线性扰动与一些实际涡结构非线性作用形态间的内在关联，印证了线性稳定性分析方法的合理性和有效性。第二部分是高精度迎风紧致差分格式（UCD）时空全离散数值行为分析。导出了其一维/二维一般色散表达式。研究表明，UCD格式在高波数区具有内在的全离散耗散和色散特性；其数值群速度的快/慢特征可因CFL数不同而改变；在稳定CFL数下简单附加人工粘性可强化UCD格式在高波数区的耗散量；提高时间精度可放宽稳定CFL数限制；UCD格式的二维全离散色散介质中存在三个不同性质的数值波，其全离散稳定性由数值声波主控。第三部分实现了高精度UCD5/SCD6差分算法对空间发展可压平面混合层的直接数值模拟。通过亚谐扰动波的个数和扰动频率的控制，捕捉到了基频涡的饱和、一次和二次对并等现象，显示了大尺度涡结构与入中初始扰动方式之间的内在联系。利用参数Mc观察了压缩性效应对大尺度涡空间演化及其相互作用的影响。第四部分实现了高精度UCD5/SCD6差分算法对非预混扩散火焰化学反应平面混合层的直接数值模拟。研究指出，放热效应可抑制和延迟涡的形成，使基频涡卷拉伸甚至丧失，混合层Reynolds 应力ρu'v'和流向速度波动关联项u'v'下降，以致涡结构与外流动量交换和标量输运减少，脉动输运能力被削弱，从而混合效率、产物生成率和混合层增长率下降，放热主要通过膨胀效应和斜压效应来抑制大尺度涡的演化。

Effects of geometric shape on the hydrodynamics of a self-propelled flapping foil

The hydrodynamics of a free flapping foil is studied numerically. The foil undergoes a forced vertical oscillation and is free to move horizontally. The effect of chord-thickness ratio is investigated by varying this parameter while fixing other ones such as the Reynolds number, the density ratio, and the flapping amplitude. Three different flow regimes have been identified when we increase the chord-thickness ratio, i.e., left-right symmetry, back-and-forth chaotic motion, and unidirectional motion with staggered vortex street. It is observed that the chord-thickness ratio can affect the symmetry-breaking bifurcation, the arrangement of vortices in the wake, and the terminal velocity of the foil. The similarity in the symmetry-breaking bifurcation of the present problem to that of a flapping body under constraint is discussed. A comparison between the dynamic behaviors of an elliptic foil and a rectangular foil at various chord-thickness ratios is also presented.

Numerical study on heavy rigid particle motion in a plane wake flow by spectral element method

Experimental particle dispersion patterns in a plane wake flow at a high Reynolds number have been predicted numerically by discrete vortex method (Phys. Fluids A 1992; 4:2244-2251; Int. J. Multiphase Flow 2000; 26:1583-1607). To address the particle motion at a moderate Reynolds number, spectral element method is employed to provide an instantaneous wake flow field for particle dynamics equations, which are solved to make a detail classification of the patterns in relation to the Stokes and Froude numbers. It is found that particle motion features only depend on the Stokes number at a high Froude number and depend on both numbers at a low Froude number. A ratio of the Stokes number to squared Froude number is introduced and threshold values of this parameter are evaluated that delineate the different regions of particle behavior. The parameter describes approximately the gravitational settling velocity divided by the characteristic velocity of wake flow. In order to present effects of particle density but preserve rigid sphere, hollow sphere particle dynamics in the plane wake flow is investigated. The evolution of hollow particle motion patterns for the increase of equivalent particle density corresponds to that of solid particle motion patterns for the decrease of particle size. Although the thresholds change a little, the parameter can still make a good qualitative classification of particle motion patterns as the inner diameter changes.

Generalization of the JTZ model to open plane wakes

The JTZ model [C. Jung, T. T¶el and E. Ziemniak, Chaos 3, (1993) 555], as a theoretical model of a plane wake behind a circular cylinder in a narrow channel at a moderate Reynolds number, has previously been employed to analyze phenomena of chaotic scattering. It is ex- tended here to describe an open plane wake without the con¯ned nar- row channel by incorporating a double row of shedding vortices into the intermediate and far wake. The extended JTZ model is found in qualitative agreement with both direct numerical simulations and ex- perimental results in describing streamlines and vorticity contours. To further validate its applications to particle transport processes, the in- teraction between small spherical particles and vortices in an extended JTZ model °ow is studied. It is shown that the particle size has signif- icant in°uences on the features of particle trajectories, which have two characteristic patterns: one is rotating around the vortex centers and the other accumulating in the exterior of vortices. Numerical results based on the extended JTZ model are found in qualitative agreement with experimental ones in the normal range of particle sizes.

Large-eddy simulation of flows past a flapping airfoil using immersed boundary method

The numerical simulation of flows past flapping foils at moderate Reynolds numbers presents two challenges to computational fluid dynamics: turbulent flows and moving boundaries. The direct forcing immersed boundary (IB) method has been devel- oped to simulate laminar flows. However, its performance in simulating turbulent flows and transitional flows with moving boundaries has not been fully evaluated. In the present work, we use the IB method to simulate fully developed turbulent channel flows and transitional flows past a stationary/plunging SD7003 airfoil. To suppress the non-physical force oscillations in the plunging case, we use the smoothed discrete delta function for interpolation in the IB method. The results of the present work demonstrate that the IB method can be used to simulate turbulent flows and transitional flows with moving boundaries.

AUTOROTATION OF TWO TANDEM TRIANGULAR CYLINDERS

The autorotation of two tandem triangular cylinders at different gap distances is investigated by numerical simulations. At the Reynolds number of 200， three distinct regimes are observed with the increase of gap distance: namely, angular oscillation, quasi-periodic autorotation and ‘chaotic’ autorotation. For various gap distances, the characteristic of vortex shedding and vortex interaction are discussed. The phase graphs (angular acceleration vs. angular velocity) and the power spectra of moment are analyzed to characterize the motion of the cylinder. The Lyapunov exponent is also calculated to identify the existence of chaos.

聚二甲基硅氧烷（PDMS）表面液滴动力学的实验研究

液滴是自然界中普遍存在的一种物质形态。非连续微流体（液滴）是近年来微流体技术重要发展方向之一。对液滴的产生、启动、移动、合并、分离和碰撞过程的研究对于航天、微纳系统、电子显示、计算机冷却、喷墨、生物医学等学科领域有着重要的应用价值。液滴属于软物质，其力学性质介于流体和固体之间，其类固体（solid-like）行为来自于曲率产生的Laplace压力和表面张力的约束。对液滴动力学行为的研究有着重要的学术价值。 本文的主要工作是针对生物微电子机械系统（Bio-MEMS）以及柔性微纳电子加工中常用的材料聚二甲基硅氧烷（Polydimethylsiloxane，PDMS）为基底的液滴动力学实验研究。 液滴是一个理想的微反应器，许多实验可以集成在一个液滴或多个液滴内完成。液滴本身的动力学特性对于实验的完成效率和质量有着重要的影响。液滴的微操控技术包括多相流法、电润湿法、热毛细管法、介电泳法等。液滴的动力学特性受到基底的影响非常大，包括基频、振动模态、运动过程等均随基底的润湿性、弹性模量的变化而有所变化。 在Bio-MEMS以及柔性微纳电子加工当中，PDMS扮演着越来越重要的角色，尤其是PDMS的润湿性和电润湿特性。目前的PDMS在Bio-MEMS当中主要是用于制备各种微流道。常见的问题主要是一方面PDMS是疏水材料，影响流体的输运。另一方面是液体在这种低Reynolds数情况下不易混合，反应效率低。本文提出了在PDMS表面溅射纳米厚度的金来减小PDMS表观接触角的方法。这种方式在特定喷金量的情况下可以在PDMS表面产生多层次的压应力波纹。这种压应力波纹对于柔性微纳电子加工，以及微流道中加速流体混合有着非常重要的作用。 电润湿是另一种可以使PDMS亲水化的方法。实验证明，PDMS具有较好的电润湿性质。此外电润湿也是目前操纵液滴的主要方式。目前一个常见的问题是电击穿现象阻碍了驱动电压的低压化，且低Reynolds数情况下液滴的混合效率偏低。此外电极还会由于少量电解的发生导致腐蚀及对液体样品的污染。本文提出了接触式的电润湿，在电极逐渐触碰液滴的过程中，液滴发生百Hz的失稳振动，稳定后接触角减小。这种电润湿模式可以有效的提高临界击穿电压，避免液滴被腐蚀后的电极污染，同时可以加快液滴的混合效率。其失稳特征时间在10 ms量级，这恰是所用液滴特征尺度在1 mm左右的电润湿器件的最快响应时间。并采用液滴振动的理论估算了液滴的失稳时间，同时还考虑了基底润湿性对液滴振动过程的影响。 液滴的启动是电润湿操控液滴过程中的重要环节。通常的液滴启动都是在非连续基底上依靠逻辑电路产生的电势变化来驱动液滴。无论是逻辑电路的设计还是驱动装置的加工都非常复杂。本文首次实现了在超疏水生物样品荷叶上的液滴启动，启动速度为数十毫米/秒，启动时间为10 ms量级。并利用PDMS成功的仿制了荷叶结构实现了超疏水的PDMS表面，荷叶同仿荷叶的PDMS超疏水表面具有相近的润湿性。 在数字微流体操控液滴的过程中，液滴的合并涉及液滴的碰撞，而且MEMS系统当中利用液滴撞击进行冷却的实验已经有所开展。同时理解液滴碰撞还对许多领域包括生物、化学、喷墨、大气物理等有着非常重要的作用。本文实验研究了Weber数和毛细数对液滴碰撞过程的影响，通过改变Weber数和毛细数得到了四种不同的响应模式。

Numerical Investigation On The Flowfield Of "Swallowtail" Cavity For Supersonic Mixing Enhancement

A "swallowtail" cavity for the supersonic combustor was proposed to serve as an efficient flame holder for scramjets by enhancing the mass exchange between the cavity and the main flow. A numerical study on the "swallowtail" cavity was conducted by solving the three-dimensional Reynolds-averaged Navier-Stokes equations implemented with a k-epsilon turbulence model in a multi-block mesh. Turbulence model and numerical algorithms were validated first, and then test cases were calculated to investigate into the mechanism of cavity flows. Numerical results demonstrated that the certain mass in the supersonic main flow was sucked into the cavity and moved spirally toward the combustor walls. After that, the flow went out of the cavity at its lateral end, and finally was efficiently mixed with the main flow. The comparison between the "swallowtail" cavity and the conventional one showed that the mass exchanged between the cavity and the main flow was enhanced by the lateral flow that was induced due to the pressure gradient inside the cavity and was driven by the three-dimensional vortex ring generated from the "swallowtail" cavity structure.

Experimental Study Of Multi-Hole Cooling For Integrally-Woven, Ceramic Matrix Composite Walls For Gas Turbine Applications

In this paper, multi-hole cooling is studied for an oxide/oxide ceramic specimen with normal injection holes and for a SiC/SiC ceramic specimen with oblique injection holes. A special purpose heat transfer tunnel was designed and built, which can provide a wide range of Reynolds numbers (10(5)similar to 10(7)) and a large temperature ratio of the primary flow to the coolant (up to 2.5). Cooling effectiveness determined by the measured surface temperature for the two types of ceramic specimens is investigated. It is found that the multi-hole cooling system for both specimens has a high cooling efficiency and it is higher for the SiC/SiC specimen than for the oxide/oxide specimen. Effects on the cooling effectiveness of parameters including blowing ratio, Reynolds number and temperature ratio, are studied. In addition, profiles of the mean velocity and temperature above the cooling surface are measured to provide further understanding of the cooling process. Duplication of the key parameters for multi-hole cooling, for a representative combustor flow condition (without radiation effects), is achieved with parameter scaling and the results show the high efficiency of multi-hole cooling for the oblique hole, SiC/SiC specimen. (C) 2008 Elsevier Ltd. All rights reserved.

A CE/SE Scheme for Flows in Porous Media and Its Applications

In this paper, a new computational scheme for solving flows in porous media was proposed. The scheme was based on an improved CE/SE method (the space-time Conservation Element and Solution Element method). We described porous flows by adopting DFB (Brinkman-Forchheimer extended Darcy) equation. The comparison between our computational results and Ghia's confirmed the high accuracy, resolution, and efficiency of our CE/SE scheme. The proposed first-order CE/SE scheme is a new reliable way for numerical simulations of flows in porous media. After investigation of effects of Darcy number on porous flow, it shows that Darcy number has dominant influence on porous flow for the Reynolds number and porosity considered.