Micromixing enhanced by pulsating flows

We consider the micromixing enhancement by pulsating flows. Dimensionless governing equations and boundary conditions were developed for T-type micromixers with two inlet pulsating flows. The problem involves a set of parameters. Three key dimensionless parameters are identified: the Reynolds number, the Strouhal number, and the disturbance amplitude. Suitable Strouhal number or disturbance amplitude causes symmetrical meniscus-shape mixing interfaces, separating the whole mixing channel into a set of segments. Thus uniform exit species concentration can be reached. Too large or too small Strouhal number or disturbance amplitude yields the meniscus-shape mixing interfaces deviating from the centerline of the mixing channel, deteriorating the mixing performance. The optimized disturbance amplitude is increased with increases in Strouhal numbers. Low Reynolds number needs larger disturbance amplitude.

Study on methane hydration process in a semi-continuous stirred tank reactor

The methane hydration process is investigated in a semi-continuous stirred tank reactor. Liquid temperatures and reaction rates without stirrer are compared with those occurring with stirrer, while at the same time better stirring conditions of the methane hydration process are given by the experiments. Some basic data of fluid mechanics, for example, stirring Reynolds number, Froucle number and stirrer power, are calculated during the methane hydration process, which can be applied to evaluate stirrer capacity and provide some basic data for a scaled up reactor. Based on experiment and calculations in this work, some conclusions are drawn. First, the stirrer has great influence on the methane hydration process. Batch stirring is helpful to improve the mass transfer and heat transfer performances of the methane hydration process. Second, induction time can be shortened effectively by use of the stirrer. Third, in this paper, the appropriate stirring velocity and stirring time were 320 rpm and 30 min, respectively, at 5.0 MPa, for which the storage capacity and reaction time were 159.1 V/V and 370 min, respectively. Under the condition of the on-flow state, the initial stirring Reynolds number of the fluid and the stirring power were 12,150 and 0.54 W, respectively. Fourth, some suggestions, for example, the use of another type of stirrer or some baffles, are proposed to accelerate the methane hydration process. Comparing with literature data, higher storage capacity and hydration rate are achieved in this work. Moreover, some fluid mechanics parameters are calculated, which can provide some references to engineering application.

草本植物对土质路面径流水动力学特征及水沙过程的影响

通过室内人工模拟降雨,研究了草本植物不同覆盖度土质路面的径流水动力学特征以及路面的产流产沙过程。结果表明,路面径流的弗劳德数和雷诺数属层流的急流范畴,并且二者与径流流速及过水断面单位能量相同,均随草本植物覆盖度的增大而减小。而Darcy-weisbach阻力系数以及曼宁糙率系数则均随覆盖度的增大而增大。在不同覆盖度土质路面上,径流量随径流时间急剧增大后趋于稳定,且二者符合移轴双曲线关系,而输沙率则经历了先增大后减小的过程。随着植被覆盖度的增大,土质道路的水分入渗率逐渐增大,而径流量、输沙率及含沙率均逐渐减小。该研究表明土质路面种植的草本植物具有较好的水土保持功能,为阐明草本植物与土壤侵蚀的关系提供了理论依据。

薄层水流水力特性试验研究

对粘砂床面的薄层水流进行定床试验,研究了不同流量、坡度下薄层水流的水力学参数(流态、水深、流速及阻力系数)随坡度和流量的变化规律。结果表明:水流雷诺数Re随着坡度和流量的增大而增大,但Re主要是由流量决定的,坡度对Re的影响不大,水流流态基本上在过渡流区和紊流区;平均水深-单宽流量,流速-单宽流量,水流阻力系数-单宽流量均成很好的幂函数形式,平均水深和流速均随单宽流量的增大而增大,阻力系数随着单宽流量的增大而减小。

黄土坡面水力学特征参数与土壤侵蚀量间关系研究

通过室内人工降雨模拟试验,研究黄土坡面水流的水力学特征参数与土壤侵蚀量间的关系。结果表明:随着放水流量的增加,出口断面平均流速、雷诺数在逐渐增大,弗劳德数在逐渐减小,坡面流流态也从层流状态逐渐过渡到紊流状态,当放水流量在10.1~19.9 L/min范围内时,坡面水流的流态为层流与紊流之间的混和状态;出口断面的平均流速随着放水流量的增加而增大,二者之间呈良好的幂函数关系;黄土坡面的累计泥沙侵蚀总量随着雷诺数的增加而增大,二者呈幂函数关系;单宽土壤侵蚀产沙量与雷诺数之间拟合关系的好坏受雷诺数取值的影响,当雷诺数的值在临界雷诺数(取500)附近时,单宽土壤侵蚀产沙量与雷诺数之间的负相关性良好。研究结果不仅可以揭示土壤侵蚀水动力过程的内在机制,而且对黄土高原地区生态环境建设具有重要的现实意义。

Thrust generation and wake structure of wiggling hydrofoil

Marine animals and micro-machines often use wiggling motion to generate thrust. The wiggling motion can be modeled by a progressive wave where its wavelength describes the flexibility of wiggling animals. In the present study, an immersed boundary method is used to simulate the flows around the wiggling hydrofoil NACA 65-010 at low Reynolds numbers. One can find from the numerical simulations that the thrust generation is largely determined by the wavelength. The thrust coefficients decrease with the increasing wavelength while the propulsive efficiency reaches a maximum at a certain wavelength due to the viscous effects. The thrust generation is associated with two different flow patterns in the wake: the well-known reversed Karman vortex streets and the vortex dipoles. Both are jet-type flows where the thrust coefficients associated with the reversed Karman vortex streets are larger than the ones associated with the vortex diploes.

Mechanism Responsible for the Complete Suppression of Karman Vortex in Flows Past a Wavy Square-Section Cylinder

The Karman vortex shedding is totally suppressed in flows past a wavy square-section cylinder at a Reynolds number of 100 and the wave steepness of 0.025. Such a phenomenon is illuminated by the numerical simulations. In the present study, the mechanism responsible for it is mainly attributed to the vertical vorticity. The geometric disturbance on the rear surface leads to the appearance of spanwise flow near the base. The specific vertical vorticity is generated on the rear surface and convecting into the near wake. The wake flow is recirculated with the appearance of the pair of recirculating cells. The interaction between the upper and lower shear layers is weakened by such cells, so that the vortex rolls could not be formed and the near wake flow becomes stable.

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 developed 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.

数值摄动算法及其CFD格式

数值摄动算法将流体动力学效应耦合进NS方程组和对流扩散(CD)方程离散的数学基本格式(MBS),特别是耦合进最简单的一阶迎风和二阶中心格式之中,由此构建成一系列新的摄动格式(PS).构建PS的主要步骤是将MBS中的通量重构为步长的幂级数,利用空间分裂和导出的高阶流体动力学线性关系式,并引入下游不影响上游的对流运动规律,通过消除重构格式修正微分方程的截断误差诸项求出幂级数的待定系数,由此获得非线性PS.PS的项是MBS中对应项与R△x(及λR△x)之简单多项式的乘积,R△x和λ分别是网格Reynolds数和网格CFL数.PS和MBS使用相同结点,简单性彼此相当,但PS精度高,稳定范围大,例如PS包含了许多绝对稳定高阶迎风和中心有限差分(FD)格式和绝对正型有限体积(FV)格式,这些格式对网格Reynolds数的任意值均为不振荡格式.数值摄动算法因此是构建高精度不振荡CFD格式的新方法.PS用于计算不可压缩流,可压缩流,液滴萃取传质,微通道两相流等,均获得良好数值结果或与已有Benchmark解一致的数值结果.已有文献称数值摄动算法为新型高精度方法和高算法,文中也讨论了一些值得进一步研究的课题

对流扩散方程的三阶迎风格式的数值摄动高精度重构

利用高智提出的数值摄动算法,把求解对流扩散方程常用三阶迎风格式(3-UDS)(粘性项和对流项分别用二阶中心格式和3-UDS离散)进行了高精度重构,包括使用离散单元内所有节点的全域重构和分别使用上下游节点的上下游重构,得到两类新的更高阶精度迎风差分格式,称为高的迎风差分格式(记作GUDS)。讨论了GUDS的数学性质,GUDS比原来的3-UDS精度显著提高;全域重构的GUDS和3-UDS均为条件稳定,一些上下游重构GUDS为绝对稳定。本文通过稳定性分析和四个算例(一维常系数、变系数、非线性及二维变系数对流扩散方程)的计算证实了GUDS的优良性质。上下游重构GUDS为避免在3-UDS中使用人工粘性提供了一条有效途径,适合于求解高Reynolds数线性和非线性问题。

对流扩散方程的绝对稳定高阶中心差分格式

将作者提出的数值摄动算法改进为区分离散单元内上游和下游并分别对通量进行高精度重构的双重数值摄动算法,与原(单重)摄动算法相比,双重摄动算法既提高了格式精度又明显扩大了格式的稳定域范围,利用双重摄动算法,即分别利用上游和下游基点变量的摄动重构将高阶流体力学关系及迎风机制耦合进二阶中心格式之中,由此构建了对流扩散方程的对网格Reynolds数的任意值均稳定(绝对稳定)高精度(四阶和八阶精度)三基点中心TVD差分格式,通过解析分析以及3个算例计算证实了构建格式的优良性能;3个算例包括一维线性、非线性(Burgers方程)和二维变系数对流扩散方程,数值计算表明:构建的格式在粗网格下不振荡,构建格式在粗网格时的最大误差L∞和均方误差L2与二阶中心格式在细网格时的相应误差一致,对线性方程,构建格式在细网格下可达到L2精度阶

基于直接数值模拟的超/高超声速湍流模型评估与改进

对来流Mach数2.25及6的平板边界层湍流进行了直接数值模拟，并通过与理论、实验及他人计算结果的对比对数值结果进行了验证。基于直接数值模拟得到的湍流数据库，对常用的湍流模型进行了先验评估。评估的湍流模型有k-e模型（包括标准k-e模型、可实现的k-e模型及低Reynolds数k-模型）、SA模型及BL模型。结果显示，可实现的k-e模型的具有较好的预测能力，而标准k-e模型预测的湍流粘性系数偏高。SA模型在边界层内层预测准确度较高，而在外层预测值偏高。对于高Mach数情况，原始的BL模型严重低估了内-外层交界位置，造成湍流粘性系数预测值偏低。作者通过修改模型系数及内-外层交界位置对BL模型进行了修改，修改后模型预测的湍流粘性系数与DNS给出的值吻合较好

Thermocapillary Migration of Deformable Bubbles at Moderate to Large Marangoni Number in Microgravity

Using the level-set method and the continuum interface model, the axisymmetric thermocapillary migration of gas bubbles in an immiscible bulk liquid with a temperature gradient at moderate to large Marangoni number is simulated numerically. Constant material properties of the two phases are assumed. Steady state of the motion can always be reached. The terminal migration velocity decreases monotonously with the increase of the Marangoni number due to the wrapping of isotherms around the front surface of the bubble. Good agreements with space experimental data and previous theoretical and numerical studies in the literature are evident. Slight deformation of bubble is observed, but no distinct influence on the motion occurs. It is also found that the influence of the convective transport of heat inside bubbles cannot be neglected at finite Marangoni number, while the influence of the convective transport of momentum inside bubbles may be actually negligible.

Linear spatio-temporal instability analysis of ice growth under a falling water film

A linear spatio-temporal stability analysis is conducted for the ice growth under a falling water film along an inclined ice plane. The full system of linear stability equations is solved by using the Chebyshev collocation method. By plotting the boundary curve between the linear absolute and convective instabilities (AI/CI) of the ice mode in the parameter plane of the Reynolds number and incline angle, it is found that the linear absolute instability exists and occurs above a minimum Reynolds number and below a maximum inclined angle. Furthermore, by plotting the critical Reynolds number curves with respect to the inclined angle for the downstream and upstream branches, the convectively unstable region is determined and divided into three parts, one of which has both downstream and upstream convectively unstable wavepackets and the other two have only downstream or upstream convectively unstable wavepacket. Finally, the effect of the Stefan number and the thickness of the ice layer on the AI/CI boundary curve is investigated.