关于水气界面处的气体传输

水气界面处的气体传输问题由于其对于地球大气中温室效应气体收支的重要意义,近年来受到广泛的关注．综述国际上在水气界面处传质研究的新进展,在简评各种基本的传质模型基础上,着重介绍在剪切水气界面附近湍流相干结构及其与气体传输的关系方面的一些新发现、新观点．对于实际海气界面气体传输的若干机制进行了评述．为实际估算海气界面处气体传输速度的目的,建议了一种综合考虑剪切与波破碎两种机制的气体传输复合计算方法,得到与实测较满意的符合．      

低密度流体界面不稳定性大涡模拟

采用拉伸涡亚格子尺度应力模型对湍流输运中的亚格子作用项进行模式化处理,发展了适用于可压多介质黏性流动和湍流的大涡模拟方法和代码MVFT（multi-viscous flow and turbulence）。利用MVFT代码对低密度流体界面不稳定性及其诱发的湍流混合问题进行了数值模拟。详细分析了扰动界面的发展,流场中冲击波的传播、相互作用、湍流混合区边界的演化规律,以及流场瞬时密度和湍动能的分布和发展。数值模拟获得的界面演化图像和流场中波系结构与实验结果吻合较好。三维和二维模拟结果的比较显示,两者得到的扰动界面位置、波系及湍流混合区边界基本一致,只是后期的界面构型有所不同,这也正说明湍流具有强三维效应。     

可压缩多介质粘性流体和湍流的大涡模拟

在可压缩多介质粘性流体动力学高精度计算方法MVPPM(multi-viscous-fluid piecewise parabolicmethod)基础上,引入Smagorinsky和Vreman亚格子湍流模型,采用大涡数值模拟方法求解可压缩粘性流体NS(Navier-Stokes)方程,给出适用于可压缩多介质流体界面不稳定性发展演化至湍流阶段的计算方法和二维计算程序MVFT(multi-viscosity-fluid and turbulence)。在2种亚格子湍流模型下计算了LANL(Los Ala-mos National Laboratory)激波管单气柱RM不稳定性实验,分析了气柱的形状、流场速度以及涡的特征,通过与LANL实验和计算结果的比较可知,Vreman模型略优于Smagorinsky模型,MVFT方法和计算程序可用于对界面不稳定性发展演化至湍流阶段的数值模拟。     

在1.4米航空风洞中模拟大气边界层

在风洞中正确模拟大气边界层的流动特性是风工程风洞试验结果可信的必要条件,本次试验研究的目的是在短试验段的航空风洞中建立大比例的大气边界层模拟流场。通过适当的方式延长1.4m×1.4m航空风洞的试验段长度,并利用尖塔、粗糙元等边界层发生装置,在该风洞中建立了边界层流场,测量了流场的平均风速剖面、湍流强度剖面、脉动风速的自相关系数、风谱等参数,讨论了湍流积分尺度的处理和大气边界层几何模拟比例的确定,用谱拟合法和自相关系数积分法求出了湍流积分尺度。结果分析表明:试验所得流场是合理的大气边界层模拟流场,其平均风速剖面幂指数α=0.3,大气边界层模拟比例为1∶500,为后续的建筑物模型动态风荷载试验提供了前提条件。     

雷诺数对湍流特性的影响

主要对近年来关于雷诺数（Ｒｅ）对湍流特性的影响的研究进行了回顾．研究表明除径向偏斜系数和平坦系数分布外,确实存在一个不受Ｒｅ数影响的区域,但对于不同的物理量这一区域的范围是不同的,即使是同一物理量不同研究者的结果亦不尽相同,显示了近壁区流动的复杂性和测量的困难,激光测速技术（ＬＤＶ）由于其具有测量精度高、不干扰流场、空间分辨率细等优点,近年来广泛地应用于对湍流,特别是对湍流近壁区流动特性及Ｒｅ数对湍流统计特性影响等的精细测量,随着ＬＤＶ自身及其它相关测试技术的进一步改进,ＬＤＶ在了解湍流产生机理、能量转换等方面发挥着越来越重要的作用．     

FV/MC混合算法求解轴对称钝体后湍流流场

介绍一种有限容积／Monte Carlo结合求解湍流流场的相容的混合算法．有限容积法求解Reynolds平均的动量方程和能量方程，Monte Carlo方法求解模化的脉动速度—频率—标量联合的PDF方程．将该算法发展到无结构网格，探讨了在无结构网格中实现两种方法的耦合，包括颗粒定位，颗粒场和平均场之间数据交换等问题．并以二维轴对称钝体后湍流流场作为算例，比较了计算结果与实验结果．     

用Level Set方法追踪运动界面

首先介绍了近年来发展起来的界面追踪技术Level Set方法,然后采用五阶WENO格式和积分平均型TVD格式计算Level Set方程,用修正的Godunov方法求解重新初始化的Level Set方程,数值求解了同心圆在常数流场,圆和矩形界面在剪切流场,缺口圆在旋转流场中的界面变形和还原效果,比较了时间导数离散精度和Level Set函数有无重新初始化对界面追踪效果的影响．最后,通过和其它界面处理方法的比较可以看出,Level Set方法不仅能够比较准确地追踪运动界面,而且无须进行复杂繁琐的界面重构技术,容易编程,具有较大的通用性．     

床面上直立圆柱的三维湍流数值模拟

从数值预报桥墩等结构物床面局部冲刷的角度发展绕直立圆柱的三维湍流的数值模拟技术. 基于Wilcox的k-ω两方程湍流模式，采用基于有限体积法的压力修正SIMPLE算法, 计算了绕床面直立圆柱的三维湍流流场，分析了光滑和粗糙床面两种情况下的流动情况. 通过系列的验证计算，表明该计算模型能够比较准确地反映不同外来流条件下绕直立圆柱的流场. 计算结果揭示了床面粗糙度对绕圆柱的湍流流动的影响.     

一类平衡大气边界层边界条件研究

基于标准k-ε湍流模型,首先利用湍流粘度方程和剪切应力在整个边界层内恒定的假设,推导出一类耗散率表达式,并根据常用的湍动能入口剖面方程以及平均风速剖面方程,计算获得相应的耗散率方程;然后在输运方程中添加自定义源项,通过已经确定的平均速度方程、湍动能方程、耗散率方程计算得到相应输运方程的自定义源项表达式,并进行空风洞数值模拟,从而得到了一类满足平衡大气边界层的来流边界条件．通过将这种边界条件与由湍流平衡条件得到的边界条件进行比较,表明本方法获得的边界条件更适用．并且,本方法无需考虑修正壁面函数和修正湍流模型常数,因而计算更为简单,可为平衡大气边界层的研究提供一种新的思路．     

用平均速度剖面法测量壁湍流摩擦阻力

用IFA300恒温热线风速仪精细测量风洞中不同雷诺数流动条件下的平板湍流边界层近壁区域对数律平均速度剖面．利用平板湍流边界层近壁区域的对数律平均速度剖面与壁面摩擦速度、流体黏性系数等内尺度物理量的关系和壁面摩擦速度与壁面摩擦切应力的关系，在准确测量平板湍流边界层近壁区域对数律平均速度剖面的基础上，测量平板湍流边界层的壁面摩擦阻力．实现了平板湍流边界层壁面摩擦阻力的无干扰或微小干扰测量．该种方法操作简便，不需要在流场中安装测力天平、传感器等复杂的测量装置，不需要对湍流边界层的壁面进行破坏，不会影响湍流边界层壁面附近区域原有的流场条件，是一种切实可行的测量平板湍流边界层壁面摩擦阻力的简便方法．     

Numerical study on the performance of a twin-raft wave energy dissipator in a stilling basin

In this paper, a raft-typed wave energy dissipator is proposed, and a mathematical model for the hydrodynamics of such a dissipator is presented, based on Reynolds-averaged Navier–Stokes equations. The model is validated by a comparison of the numerical results with the results of other investigators. The validated model is then utilized to examine the effect of wave height, wave frequency, damping coefficient, flow velocity on wave energy dissipation ratio and wave transmission coefficient for a hinged twin-raft wave energy dissipator. Our results reveal that the differences in behaviour exhibited by an inviscid fluid and a viscous fluid can be large and vary considerably, depending on the flow velocity.     

Plastic deformation modelled as a self-excited wave process at the meso- and macro-level

A self-excited wave model is developed to describe plastic flow phenomena in crystalline solids. Experimental observations suggest that by plastic flow in single crystals and polycrystalline materials, different underlying mechanisms are responsible for key features of strain localisation corresponding to different stages of the deformation curve. The major autowave (self-excited wave) types manifest themselves in plastically deforming materials. The self-excited wave model could explain plastic flow pattern behaviour corresponding to different physical mechanisms.     

A model for pressure loss, film thickness, and entrained fraction for gas-liquid annular flow

A two-component (air-water) annular flow model is presented requiring only flow rates, absolute pressure, temperature, and tube diameter. Film thicknesses (base film and wave height) are calculated from a critical film thickness model. Modeled pressure gradient is weighted by wave intermittency to compute average pressure gradient. Film flow rate and wave velocity are estimated using the universal velocity profile in the waves and a piecewise linear profile in the base film. For vertical flow, mean absolute errors for film thickness, wave velocity, and pressure gradient are 9%, 9%, and 19%, respectively. In horizontal flow, mean absolute errors for pressure gradient, base film thickness, and disturbance wave velocity are 17%, 10%, and 14%, respectively, on par with those from single-behavior models that require additional film thickness or other data as inputs.     

颗粒材料中致密波结构研究

采用一维两相流模型与相应颗粒构形应力函数,研究了致密波的形成及其结构．用简化两相流模型系统地讨论致密波对有关因素的依赖关系．分析指出：小于基体材料音速的致密波仅能在非理想颗粒材料中存在,从波前到波后,所有状态物理量光滑过渡．大于基体材料音速的致密波,波头可能存在间断．应力函数与致密粘性确定后,致密波速度决定致密波结构、宽度、终态压实度．采用一维两相流模型模拟了活塞驱动颗粒床形成致密波这一动态过程．用线方法（MOL）对该方程组求数值解．计算表明,经过短暂的非稳态过程,颗粒床中形成一稳态致密波．分析了活塞速度与初始孔隙率对致密波结构的影响,并对简化两相流模型与两相流模型的计算结果进行了对比．     

基于势流理论的内孤立波与立管作用数值研究

内孤立波是一种发生在水面以下的在世界各个海域广泛存在的大幅波浪, 其剧烈的波面起伏所携带的巨大能量对以海洋立管为代表的海洋结构物产生严重威胁, 分析其传播演化过程的流场特征及立管在内孤立波作用下的动力响应规律对于海洋立管的设计具有重要意义. 本文基于分层流体的非线性势流理论, 采用高效率的多域边界单元法, 建立了内孤立波流场分析计算的数值模型, 可以实时获得内孤立波的流场特征. 根据获得的流场信息, 采用莫里森方程计算内孤立波对海洋立管作用的载荷分布. 将内孤立波流场非线性势流计算模型与动力学有限元模型结合来求解内孤立波作用下海洋立管的动力响应特征, 讨论了内孤立波参数、顶张力大小以及内部流体密度对立管动力响应的影响. 发现随着内孤立波波幅的增大, 海洋立管的流向位移和应力明显增大. 由于上层流体速度明显大于下层, 且在所研究问题中拖曳力远大于惯性力, 因此管道顺流向的最大位移发生在上层区域. 顶张力通过改变几何刚度阵的值进而对立管的响应产生明显影响. 对于弱约束立管, 内部流体的密度对管道的流向位移影响较小.      

KINEMATIC WAVE PROPERTIES OF ANISOTROPIC DYNAMICS MODEL FOR TRAFFIC FLOW

ＩｎｔｒｏｄｕｃｔｉｏｎＲｅｃｅｎｔｌｙｔｈｅａｕｔｈｏｒｓｐｒｅｓｅｎｔｅｄａｎｅｗｄｙｎａｍｉｃｓｍｏｄｅｌｆｏｒｔｒａｆｆｉｃｆｌｏｗ[1].Ｔｈｅｍｏｄｅｌｃｏｍｐｒｉｓｅｓａｃｏｎｔｉｎｕｕｍｅｑｕａｔｉｏｎａｎｄａｄｙｎａｍｉｃｓｅｑｕａｔｉｏｎ ｋ/ ｔ+ (ｋｕ) / ｘ=0 ,( 1 ) ｕ ｔ+ｕ ｕ ｘ=ｕｅ(ｋ) -ｕＴ +ａ ｕ ｘ,( 2 )ｗｈｅｒｅｋｉｓｔｒａｆｆｉｃｄｅｎｓｉｔｙ ;ｕｉｓｍｅａｎｓｐｅｅｄ ;ｘ ,ｔａｒｅｓｐａｃｅａｎｄｔｉｍｅｃｏｏｒｄｉｎａｔｅｓｒｅｓｐｅｃｔｉｖｅｌｙ .Ｔｉｓｒｅｌ…     

The mathematical model of the radial artery blood pressure waveform through monitoring of the age-related changes

Age-related changes in blood vessels affect the pulse wave propagation. These changes may cause an increase in wave reflection which leads to amplification of pulse pressure. The pulse pressure changes are associated with certain vascular diseases. Here we present a mathematical model of blood pressure for different age groups. The model is based on one-dimensional wave propagation theory and assumes that the pressure waveform is a superposition of forward propagating wave and backward waves from many reflection sites. The model is based on experimental data obtained by direct measurement of radial artery blood flow. The model clearly shows age-related changes in blood pressure waveform. The results of mathematical model correlate with the radial pressure waveform data. The model can be used in cases where it is not possible to measure the pressure due to movement of subject. Application of model to the direct blood flow measurements data allows the real-time pressure waveform monitoring. Furthermore, this approach enables monitoring of changes in pressure waveform due to the effects of medications.     

Generation of linear and nonlinear waves in numerical wave tank using clustering technique-volume of fluid method

A two-dimensional(2D) numerical model is developed for the wave simulation and propagation in a wave flume.The fluid flow is assumed to be viscous and incompressible,and the Navier-Stokes and continuity equations are used as the governing equations.The standard k-ε model is used to model the turbulent flow.The NavierStokes equations are discretized using the staggered grid finite difference method and solved by the simplified marker and cell(SMAC) method.Waves are generated and propagated using a piston type wave maker.An open boundary condition is used at the end of the numerical flume.Some standard tests,such as the lid-driven cavity,the constant unidirectional velocity field,the shearing flow,and the dam-break on the dry bed,are performed to valid the model.To demonstrate the capability and accuracy of the present method,the results of generated waves are compared with available wave theories.Finally,the clustering technique(CT) is used for the mesh generation,and the best condition is suggested.     

Numerical study of breaking waves by a two‐phase flow model

A two‐phase flow model, which solves the flow in the air and water simultaneously, is presented for modelling breaking waves in deep and shallow water, including wave pre‐breaking, overturning and post‐breaking processes. The model is based on the Reynolds‐averaged Navier–Stokes equations with the k ?ε turbulence model. The governing equations are solved by the finite volume method in a Cartesian staggered grid and the partial cell treatment is implemented to deal with complex geometries. The SIMPLE algorithm is utilised for the pressure‐velocity coupling and the air‐water interface is modelled by the interface capturing method via a high resolution volume of fluid scheme. The numerical model is validated by simulating overturning waves on a sloping beach and over a reef, and deep‐water breaking waves in a periodic domain, in which good agreement between numerical results and available experimental measurements for the water surface profiles during wave overturning is obtained. The overturning jet, air entrainment and splash‐up during wave breaking have been captured by the two‐phase flow model, which demonstrates the capability of the model to simulate free surface flow and wave breaking problems.Copyright © 2011 John Wiley & Sons, Ltd.     

A coupled boundary element–finite difference model of surface wave motion over a wall turbulent flow

An effective numerical technique is presented to model turbulent motion of a standing surface wave in a tank. The equations of motion for turbulent boundary layers at the solid surfaces are coupled with the potential flow in the bulk of the fluid, and a mixed BEM–finite difference technique is used to model the wave motion and the corresponding boundary layer flow. A mixing‐length theory is used for turbulence modelling. The model results are in good agreement with previous physical and numerical experiments. Although the technique is presented for a standing surface wave, it can be easily applied to other free surface problems. Copyright © 2005 John Wiley & Sons, Ltd.