结构振动对湍流近尾迹的影响

研究了圆柱绕流中流体与结构的相互作用,侧重结构振动对湍流尾迹的影响,用激光测振仪测量圆柱在升力方向的位移;用热线和LDA（二维）测量湍流的近尾迹,通过变化自由流的速度和圆柱体直径（特征尺寸）来变化雷诺数,用两个振动特性不同的（一个相对刚性,一个相对弹）圆柱来产生尾迹,研究固体结构振动对湍流近尾迹的平均速度场和湍流场的影响,结果表明,结构自由振动对湍流近尾迹场影响明显,该影响随雷诺数的变化不明显。     

湍流边界层相干结构剩余脉动雷诺应力项的实验测量

用IFA300恒温热线风速仪和X形二分量热线探针以采样间隔小于湍流耗散时间尺度的分辨率精细测量了风洞中平板湍流边界层不同法向位置的瞬时流向、法向速度分量的时间序列信号。用条件采样和相位平均技术提取了相干结构猝发过程中相干结构剩余脉动雷诺应力和随机脉动对相干结构贡献的雷诺应力的条件相位平均波形。基于理论上对湍流相干结构动量方程中随机脉动对相干结构贡献的雷诺应力和相干结构剩余脉动雷诺应力项的分析,对两种雷诺应力项进行了对比研究。研究发现,相干结构剩余脉动雷诺应力项在数值上具有和随机脉动对相干结构贡献的雷诺应力相同的数量级,表明在相干结构动力学模型方程中,相干结构剩余脉动雷诺应力项并不像以前估计的那样可以忽略不计。     

壁湍流猝发过程中速度分量的相位差对雷诺应力影响的实验研究

用IFA300恒温热线风速仪和×形二分量热线探针,以采样间隔小于最小湍流时间尺度的分辨率,精细测量了风洞中平板湍流边界层不同法向位置的瞬时流向、展向速度分量的时间序列信号.用子波分析辨识壁湍流相干结构猝发事件的能量最大准则,确定壁湍流相干结构猝发事件的时间尺度;用条件相位平均技术提取了相干结构猝发过程中流向、展向脉动速度分量条件相位平均波形,用互相关方法研究了相干结构猝发过程中流向、展向脉动速度分量条件相位平均波形的相位差关系及其对雷诺应力的影响,发现在缓冲层和对数律区,展向脉动速度与流向脉动速度的条件相位平均波形具有不同的相位;当两者相位基本一致时,雷诺应力达到正的最大值,此时湍流相干结构的产生非常活跃;当两者相位差分别集中在90°和270°附近时,雷诺应力的幅值减小并接近于零,此时湍流相干结构的产生和猝发都得到了抑制.     

壁面加热湍流相干结构的实验研究

用热线风速仪测量了风洞中壁面加热平板湍流边界层不同法向位置的流向脉动速度,用条件采样方法中的m-u-level法和VITA法研究了壁面加热对壁湍流相干结构猝发条件平均波形的影响,发现壁面另热能够提高猝发强度,缩短喷射时间,使一次猝发中出现多次喷射事件的现象明显增多,增强猝发的间歇性。     

TRPIV MEASUREMENT OF DRAG-REDUCTION IN THE TURBULENT BOUNDARY LAYER OVER RIBLETS PLATE

采用高时间分辨率粒子图像测速技术对沟槽壁面平板湍流边界层速度矢量场的时间序列及其统计量进行了实验测量,讨论了在同一来流速度下沟槽壁面对平均速度剖面﹑雷诺切应力及湍流强度的影响. 用流向速度分量的多尺度空间局部平均结构函数辨识壁湍流多尺度相干结构,用条件采样和相位平均技术提取壁湍流多尺度相干结构喷射和扫掠事件的脉动速度、展向涡量的二维空间拓扑形态. 结果表明,与同材料光滑壁面对比,沟槽壁面实现了10.73%的摩阻减小量;沟槽壁面湍流边界层湍流强度及雷诺切应力皆比光滑平板湍流边界层对应统计量小,说明沟槽壁面有效降低了湍流边界层内流体的脉动. 通过比较壁湍流相干结构猝发事件各脉动速度分量与展向涡量的空间分布特征,肯定了沟槽壁面的减阻效果,发现沟槽壁面通过抑制相干结构猝发事件实现减阻.     

流向磁场作用下圆柱绕流的直接数值模拟

绕流是托卡马克装置中液态包层内常见的流动形态,对流场与热量分布有着重要的影响.本文通过直接数值模拟(DNS),研究了不同磁场强度下$Re=3900$的圆柱绕流,分析了磁场强度对于湍流尾迹的影响.无磁场情况下,直接数值模拟的结果与前人的实验及模拟结果吻合很好.圆柱下游的尾迹中,随着流向距离的增大, 流向速度剖面逐渐从U型进化呈V型, 并慢慢趋于平缓,这表明尾迹中的流动结构受圆柱影响逐渐减小.圆柱后方两侧的剪切层中,由于Kelvin-Helmholtz不稳定性的影响,可以清晰地看到小尺度剪切层涡的脱落.通过对无磁场的计算结果施加流向磁场,本文计算了哈特曼数($Ha$)分别为20, 40和80的工况,以研究磁场效应对于湍流的影响.结果表明磁场较弱时,流动依然呈三维湍流状态.随着磁场增强, 近圆柱尾流区受磁场抑制明显,回流区被拉长,剪切层失稳位置向下游转移.圆柱后方的涡结构由于受到竖直方向洛伦兹力的挤压作用,随着哈特曼数的增加尾迹区域逐渐变窄.相比于无磁场情况的涡结构,由于磁场的耗散作用,相应的涡结构尺度变小.该研究不仅扩展了现有磁场下湍流运动的参数范围,对于液态包层的设计及安全运行同样具有重要的理论指导意义和工程应用价值.      

DIRECT NUMERICAL SIMULATIONS ON THE TURBULENT FLOW PAST A CONFINED CIRCULAR CYLINDER WITH THE INFLUENCE OF THE STREAMWISE MAGNETIC FIELDS1)

绕流是托卡马克装置中液态包层内常见的流动形态,对流场与热量分布有着重要的影响.本文通过直接数值模拟(DNS),研究了不同磁场强度下$Re=3900$的圆柱绕流,分析了磁场强度对于湍流尾迹的影响.无磁场情况下,直接数值模拟的结果与前人的实验及模拟结果吻合很好.圆柱下游的尾迹中,随着流向距离的增大, 流向速度剖面逐渐从U型进化呈V型, 并慢慢趋于平缓,这表明尾迹中的流动结构受圆柱影响逐渐减小.圆柱后方两侧的剪切层中,由于Kelvin-Helmholtz不稳定性的影响,可以清晰地看到小尺度剪切层涡的脱落.通过对无磁场的计算结果施加流向磁场,本文计算了哈特曼数($Ha$)分别为20, 40和80的工况,以研究磁场效应对于湍流的影响.结果表明磁场较弱时,流动依然呈三维湍流状态.随着磁场增强, 近圆柱尾流区受磁场抑制明显,回流区被拉长,剪切层失稳位置向下游转移.圆柱后方的涡结构由于受到竖直方向洛伦兹力的挤压作用,随着哈特曼数的增加尾迹区域逐渐变窄.相比于无磁场情况的涡结构,由于磁场的耗散作用,相应的涡结构尺度变小.该研究不仅扩展了现有磁场下湍流运动的参数范围,对于液态包层的设计及安全运行同样具有重要的理论指导意义和工程应用价值.     

雷诺应力各向异性涡黏模型的层析TRPIV测量

利用层析TRPIV测量水洞中平板湍流边界层3D-3C速度场的高分辨率时间序列数据库. 提出了空间局部平均多尺度速度结构函数的新概念, 描述湍流多尺度涡结构的空间拉伸、压缩、剪切变形和旋转. 用空间局部平均多尺度速度结构函数对湍流脉动速度进行了空间多尺度分解. 用空间流向局部平均多尺度速度结构函数, 根据湍流多尺度涡结构在流向的拉伸和压缩物理特征, 提出了新的湍流相干结构条件采样方法, 检测并提取了层析TRPIV数据中相干结构“喷射”和“扫掠”事件中的脉动速度、平均速度变形率、雷诺应力等物理量的空间拓扑形态. 通过研究平均速度变形率各分量与雷诺应力各分量之间的空间相位差异,肯定了壁湍流相干结构雷诺应力各向异性复涡黏模型的合理性.      

输气管道壁面涂料减阻机理的实验研究

用IFA-300热线风速仪以高于对应最小湍流时间尺度的分辨率精细测量了风洞中不同壁面涂料的管道湍流边界层不同法向位置流向速度分量的时间序列信号,利用湍流边界层近壁区域对数律平均速度剖面与壁面摩擦速度、流体黏性系数等内尺度物理量的关系和壁面摩擦速度与壁面摩擦切应力的关系,在准确测量湍流边界层近壁区域对数律平均速度剖面的基础上,间接测量湍流边界层的壁面摩擦阻力．对不同壁面涂料的壁湍流脉动速度信号用子波分析进行多尺度分解,用子波系数的瞬时强度因子和平坦因子检测管道湍流边界层中的多尺度相干结构,提取不同尺度相干结构的条件相位平均波形,对比研究输气管道壁面涂料的减阻机理．     

基于单个压电振子的湍流边界层主动控制

利用安装在壁面上的单个压电振子周期振荡,采用开环主动控制方案,实现了对平板湍流边界层相干结构猝发的主动控制和壁湍流减阻.根据不同的输入电压幅值和频率,完成了10种工况的实验.在压电振子下游2mm处,用热线风速仪和迷你热线单丝探针,精细测量湍流边界层不同法向位置瞬时流向速度信号的时间序列,分析了在Re_?=2183压电振子振动对湍流边界层平均速度剖面、减阻率和相干结构猝发过程的影响.实验结果表明,施加控制的工况使平均速度剖面对数律层上移,产生减阻效果;压电振子振幅越大,减阻率越高,减阻效果越明显;通过对施加控制前后流向瞬时速度的多尺度湍涡结构脉动动能的尺度分析,当压电振子振动频率与壁湍流能量最大尺度的猝发频率相近时,减阻率达到最大,为25%,说明控制壁湍流能量最大尺度相干结构的猝发是实现壁湍流减阻的关键;通过对比相干结构猝发的流向速度分量条件相位平均波形,发现施加控制的工况中相干结构猝发流向速度分量的波形幅值明显降低,且流向速度在扫掠后期高速阶段迅速衰减,缩短了高速流体的下扫过程,说明压电振子的振动能抑制相干结构的高速流体下扫过程,减弱高速流体与壁面的强烈剪切过程,并使近壁区域相干结构的振幅显著减弱,迁移速度加快,从而减小壁面摩擦阻力.      

The Role of Biofilm Growth in Bacteria-Facilitated Contaminant Transport in Porous Media

Groundwater contaminants adhered to colloid surfaces may migrate to greater distances than predicted by using the conventional advective-dispersive transport equation. Introduction of exogenous bacteria in a bioremediation operation or mobilization of indigenous bacteria in groundwater aquifers can enhance the transport of contaminants in groundwater by reducing the retardation effects. Because of their colloidal size and favorable surface conditions, bacteria can be efficient contaminant carriers. In cases where contaminants have low mobility because of their high partition with aquifer solids, facilitated contaminant transport by mobile bacteria can create high contaminant fluxes. In this paper, we developed a methodology to describe the bacteria-facilitated contaminant transport in a subsurface environment using the biofilm theory. The model is based on mass balance equations for bacteria and contaminant. The contaminant is utilized as a substrate for bacterial growth. Bacteria are attached to solid surfaces as a biofilm. We investigate the role of the contaminant adsorption on both biofilm and mobile bacteria on groundwater contaminant transport. Also, the effect of bacterial injection on the contaminant transport is evaluated in the presence of indigenous bacteria in porous media. The model was solved numerically and validated by experimental data reported in the literature. Sensitivity analyses were conducted to deduce the effect of critical model parameters. Results show that biofilm grows rapidly near the top of the column where the bacteria and contaminant are injected, and is detached by increasing fluid shear stress and re-attach downstream. The adsorption of contaminant on bacterial surfaces reduces contaminant mobility remarkably in the presence of a biofilm. The contaminant concentration decreases significantly along the biofilm when contaminant partition into bacteria. Bacterial injection and migration in subsurface environments can be important in bioremediation operations regardless of the presence of indigenous bacteria.     

Analysis of saltation characteristics of bed-load transport in flowing water

The study of bed-load transport is of great significance both in theory and in practice. This paper discusses the saltation of bed-load solid grains in flowing water. Experiments and theoretic analysis have been made by means of high-speed photographing and advanced data processing technique with a combined method based on mechanical and statistical theories. It indicates that the saltation is the main form of the bed-load transport under ordinary flowing conditions. In the meantime, taking successive saltation as the model of bed-load transport, systematic analysis has been made with regard to the kinematic properties and mechanism of saltation. The statistical analysis shows that the probability density functions of the relative height and length of saltation are in conformity with Γ-type distribution, while the probability density functions of the relative velocities of saltation are in conformity with the Gaussian distribution. The project supported by National Natural Science Foundation of China     

Modeling the Transport of Contaminants Originating from the Dissolution of DNAPL Pools in Aquifers in the Presence of Dissolved Humic Substances

A twodimensional finite difference numerical model was developed to describe the transport of dissolved organics originating from nonaqueous phase liquid (NAPL) pool dissolution in saturated porous media in the presence of dissolved humic substances. A rectangular NAPL pool was considered in a homogeneous porous medium with unidirectional interstitial groundwater velocity. It was assumed that dissolved humic substances and aqueous phase contaminants may sorb onto the solid matrix under local equilibrium conditions. The contaminant in the aqueous phase may undergo firstorder decay. Also, the dissolved contaminant may sorb onto humic substances. The transport properties of dissolved humic substances are assumed to be unaffected by sorbing contaminants, because dissolved humic macromolecules are much larger than dissolved contaminants and sorption of nonpolar contaminants onto humic substances do not affect the overall surface charge of humic substances. The sorption characteristics of dissolved humic substances onto clean sand were determined from column experiments. An effective local mass transfer rate coefficient accounting for the presence of dissolved humic substances was developed. Model simulations indicate that dissolved humic substances substantially increase NAPL pool dissolution, and consequently reduce the required pumpandtreat aquifer remediation time.     

On the control problem for a compound wheeled vehicle

A compound transport vehicle with one steerable wheel is considered as a controlled mechanical system with nonholonomic constraints. A control system for this system is synthesized within the framework of the kinematic approximation. A procedure is proposed to synthesize a robust controller for the system. The efficiency of this procedure is demonstrated with an example __________ Translated from Prikladnaya Mekhanika, Vol. 43, No. 11, pp. 105–112, November 2007.     

APPLICATION OF COMPUTATIONAL FLUID DYNAMICS IN THE PLANNED DITCHING OF A TRANSPORT AIRPLANE

论述了边界元法、有限元法、光滑粒子水动力学法、有限体积法、工程算法与混合分析法在飞行器水上迫降中的应用和发展,总结了不同方法的优缺点,指明了未来的研究重点,为进一步开展水上迫降研究提供参考.     

脑Willis环的一维血流动力学及氧输运特性的数值研究

Willis环是大脑侧枝循环的重要组成部分, 研究其血流动力学特性以及氧输运规律对脑缺血疾病的认知和预防有着非常重要的作用. 该文旨在利用一维血流动力学模型模拟整个Willis环的流量变化和压力分布, 并建立动脉内氧输运的一维模型以模拟Willis 环内氧分压的变化规律, 为研究脑组织内血液流动和氧输运打下基础. 首先, 基于弹性圆管内的一维非线性流动方程和状态方程建立血流动力学模型, 在一维对流扩散方程的基础上, 考虑由管腔向壁面的扩散和壁面细胞的新陈代谢消耗推导出氧输运特性方程. 通过 Lax-Wendroff两步法对血流动力学方程进行离散, 而在进行对流扩散方程的离散时, 则运用迎风格式. 通过数值计算得到了正常情况下Willis环各个血管任意位置的流量、压力和氧分压的变化曲线, 正常情况下各个位置的氧分压处于稳定的平衡状态. 最后, 还通过此模型进一步模拟了右侧颈内动脉狭窄对各个血管内流动的影响. 当狭窄程度达到80%时, 中脑动脉的流量和压力会明显下降, 造成其供应区域的血流减少. 同时, Willis环右侧血管内的氧分压会大大降低, 而左侧血管的氧分压会出现上升趋势, 但幅度要小于右侧血管降低的幅度.     

Induced flow and ascent of heavy particles when an air intake operates near a surface

The structure and the properties of the air flow developing when an intake system (in particular, a turbojet engine) operates near a surface are studied. On the basis of numerical calculations, the conditions of onset, the intensity, and the structure of the vortex flow ascending from the surface toward the intake are determined. From an analysis of the particle transport toward the intake by the vortex flow, the main factors responsible for the rise of large particles and objects are determined and the conditions of minimization of the probability of large particles entering the intake are found.     

Computing brine transport in porous media with an adaptive-grid method

An adaptive-grid finite-difference method is applied to a model for non-isothermal, coupled flow and transport of brine in porous media. In the vicinity of rock salt formations the salt concentration in the fluid becomes large, giving rise to disparate scales in the salt concentrations profiles. A typical situation one encounters is that of a sharp freshwater-saltwater interface that moves in time. In such situations adaptive-grid methods are more effective than standard fixed-grid methods, since they refine the space grid locally and, hence, provide for substantial reduction in the number of grid points, memory use and CPU time. The adaptive-grid method of this paper is a static, local uniform grid refinement method. Its main feature is that it integrates on nested sequences of locally uniformly refined Cartesian space grids, which are automatically adjusted in time to follow rapid spatial transitions. Variable time steps are used to cope with rapid temporal transitions, including a fast march to possible steady-state solutions. For time stepping, the implicit, second-order BDF scheme is used. Two specific example problems are numerically illustrated. The main physical properties involved here are advection and dispersion and in case of dominant advection sharp freshwater-saltwater interfaces arise.     

Non-passive Transport of Volatile Organic Compounds Infiltrated from a Surface Disk Source

An unsaturated flow and non-passive transport model for water-soluble organic compounds has been implemented in cylindrical coordinates, with a top boundary condition that accounts for different zones of the surface that can be under infiltration or volatilization independently. We simulated two-dimensional infiltration of aqueous mixtures of methanol from a disk source, its redistribution and volatilization in both homogeneous and heterogeneous soils. Simulations showed that the incoming composition significantly affects volumetric liquid content and concentration profiles, as well as the fraction of infiltrated mass of methanol that is released to the atmosphere. Concentration-dependent viscosity had the major impact on the liquid flow. The differences in volumetric liquid content and normalized concentration of methanol became more pronounced during transport through a soil composed of a clay lens embedded within a main matrix of sandy clay loam texture. Dispersion in the liquid-phase was the predominant transport mechanism when dispersivity at saturation was set to 7.8 cm. However, for dispersivity of 1.0 cm, changes in composition led to changes in surface tension inducing significantly higher liquid flow. In this case, liquid-phase advection was the most active transport mechanism for homogeneous soils and highly concentrated infiltrating mixtures.