Linear stability of a fluid channel with a porous layer in the center简

We perform a Poiseuille flow in a channel linear stability analysis of a inserted with one porous layer in the centre, and focus mainly on the effect of porous filling ratio. The spectral collocation technique is adopted to solve the coupled linear stability problem. We investigate the effect of permeability, σ, with fixed porous filling ratio ψ = 1/3 and then the effect of change in porous filling ratio. As shown in the paper, with increasing σ, almost each eigenvalue on the upper left branch has two subbranches at ψ = 1/3. The channel flow with one porous layer inserted at its middle （ψ = 1/3） is more stable than the structure of two porous layers at upper and bottom walls with the same parameters. By decreasing the filling ratio ψ, the modes on the upper left branch are almost in pairs and move in opposite directions, especially one of the two unstable modes moves back to a stable mode, while the other becomes more instable. It is concluded that there are at most two unstable modes with decreasing filling ratio ψ. By analyzing the relation between ψ and the maximum imaginary part of the streamwise phase speed, Cimax, we find that increasing Re has a destabilizing effect and the effect is more obvious for small Re, where ψ a remarkable drop in Cimax can be observed. The most unstable mode is more sensitive at small filling ratio ψ, and decreasing ψ can not always increase the linear stability. There is a maximum value of Cimax which appears at a small porous filling ratio when Re is larger than 2 000. And the value of filling ratio 0 corresponding to the maximum value of Cimax in the most unstable state is increased with in- creasing Re. There is a critical value of porous filling ratio （= 0.24） for Re = 500; the structure will become stable as ψ grows to surpass the threshold of 0.24; When porous filling ratio ψ increases from 0.4 to 0.6, there is hardly any changes in the values of Cimax. We have also observed that the critical Reynolds number is especially sensitive for small ψ where the fastest drop is observed, and there may be a wide range in which the porous filling ratio has less effect on the stability （ψ ranges from 0.2 to 0.6 at σ = 0.002）. At larger permeability, σ, the critical Reynolds number tends to converge no matter what the value of porous filling ratio is.     

一种投影浸入边界法的快速直接求解方法

根据投影浸入边界法分步投影求解的特点,同时针对压力泊松方程离散后的大型稀疏线性方程组是非奇异非对称的特点,结合开源函数库UMFPACK,在传递线性方程组的系数矩阵和右端向量时,采用函数库Eigen将系数矩阵的数据结构改写优化,大大降低了存储空间,实现对高维大型稀疏线性方程组的快速求解,同时求解保持良好的稳定性。本文首先利用一具有解析解的数值算例验证了求解泊松方程数值方法的准确性和网格依赖性,进而利用VC++编写投影浸入边界法的数值计算程序,以单圆柱绕流为基准数值算例,通过与其他文献和实验结果的对比,验证了投影浸入边界法数值计算结果的可靠性,并进一步分析了不同雷诺数下圆柱绕流的流场结构特征和尾涡结构的动态演化过程。     

A qualitative study on the pulsatile flow phenomenon in a dense fly ash pneumatic conveyor

Understanding of the dynamic particulate flow structures within a dense gas-fly ash pneumatic conveyor must be improved in order to better aid its design guidance.The complex pulsatile movement of the gas-fly ash mixture dominates the flow performance within the pipeline,and historically,non-invasive measurement devices such as the electrical capacitance tomography（ECT） were often used to sufficiently capture the flow dynamics.However,inadequate studies have been conducted on the pulsatile flow phenomenon,which directly relate to the gas-fly ash two-phase flow performance.This paper aims to investigate the pulsatile flows using an ECT device.Initially,pulsatile flow patterns under various experimental conditions were obtained through ECT.Pulses within a flow were then characterised into pulse growth and decay segments,which represent the superficial fluidisation and deaeration processes during conveying.Subsequently,structural and statistical analyses were performed on the pulse growth and decay segments.Results suggested that the increasing air mass flow rate led to the decrease of the superficial fluidisation/deaeration magnitude,however,the increase of the superficial fluidisation/deaeration durations.Also,the air mass flow rate was indicated as the dominant factor in determining the pulsing statistical parameters.This research provides fundamental insights for further modelling the dense fly ash pneumatic flows.     

壁面结构对三维可压缩气泡群影响的数值模拟研究

基于流体体积(VOF)法追踪自由液面,研究了壁面结构对三维可压缩气泡群流动的影响.通过在待测壁面上设置不同形状的壁面结构(长方体、椭球体和圆锥体)并改变它们各自的几何参数(位置和长度),来研究壁面结构对壁面附近的气泡群流动的影响,该影响表现为气泡群对壁面的空间平均压力.研究发现,壁面结构对气泡群的拓扑结构的影响会造成壁...     

基于转捩SST模型凸起圆柱绕流数值研究

为了研究局部凸起对边界层转捩的影响,采用转捩SST模型分别对亚临界、临界和超临界状态下带突起的圆柱绕流问题进行了数值模拟,分析了不同Reynolds数下带突起的圆柱绕流问题的近壁面流动特征以及表面时均压力与摩擦力系数的分布和凸起对圆柱表面流动分离以及转捩的影响,对比了有无凸起两侧圆柱表面时均压力、摩擦力系数的不同. 结果表明:当来流Reynolds数处于临界区时,气流在圆柱上表面凸起处形成了3个反向旋转的漩涡,之后随着θ的增大,发生了流动分离和流动转捩现象;对于不同Reynolds数下的来流,圆柱上表面的凸起可以使气流发生转捩的位置提前;圆柱上表面的凸起使流速增大、压强降低,从而导致圆柱产生升力,随着来流Reynolds数的增大,其升力逐渐变大.      

格子Boltzmann方法分析气泡的运动及其相互作用

利用格子Boltzmann方法对气泡在液体中的运动规律及相互作用的气液两相流问题进行了研究;采用了基于自由能模型且适用于大密度比多相流问题的Z-S-C模型,计算模拟了处于不可压缩流体中的气泡在浮力作用下的运动特性;通过算例验证了Z-S-C模型的可行性;研究了不同放置位置两气泡的运动形态及其相互作用情况,分别分析了沿y方向同轴放置的两气泡间距变化对合并过程和沿x方向并列放置的两气泡的间距变化对气泡相互作用的影响。结果表明:两气泡的初始间距在一定范围内时,初始间距对上升过程中气泡间的相互作用具有重要影响;当初始间距超过一定量值时,上升过程中两气泡相互独立,不存在相互作用。     

玉米根茬剪切力学性能试验与分析

为提高免耕播种机破茬刀的工作性能,需要分析剪切速度、切刀刃型及剪切方式对根茬剪切性能的影响规律,采用微机控制电子万能试验机对玉米根茬进行上述三因素随机区组剪切性能试验。结果表明:剪切方式对最大剪切力的影响最大,其次是剪切速度;切刀刃型对最大剪切力的影响较小;最大剪切力随着剪切速度的增大而减小;同一刃型切刀在相同剪切速度下,横切时最大剪切力最大,劈切时最大剪切力最小;同一剪切方式在相同剪切速度下,直刃切刀的最大剪切力最大,凸圆刃切刀的最大剪切力最小;不同剪切方式,切茬过程的负荷-变形曲线不同。     

斜齿轮弹流润滑下的接触疲劳寿命计算

经典齿轮接触疲劳强度理论是基于光滑表面赫兹干接触理论,而实际齿面具有粗糙度,且啮合轮齿多数处于混合润滑状态.本文基于齿轮润滑接触分析建立了渐开线斜齿轮的接触疲劳寿命计算模型.模型由齿轮润滑接触分析模型和基于次表面应力分布的疲劳寿命模型组成.首先将斜齿圆柱齿轮一对齿的瞬时啮合等效为两反向圆锥的接触问题,建立了齿轮的有限长弹流润滑计算模型,考虑了齿轮啮合周期内瞬时载荷、接触线长、卷吸速度等因素的影响,基于统一雷诺方程方法求得啮合齿对间的润滑压力和油膜厚度分布;在此基础上,计算轮齿接触区次表面的米歇斯应力分布,根据Zaretsky接触疲劳寿命计算模型,对齿轮组的接触疲劳寿命进行模拟预测.针对不同工况参数下接触疲劳寿命计算表明:润滑油黏度、轮齿表面粗糙度等因素对齿面接触疲劳寿命均有显著的影响.     

酰胺接枝ACF的制备及其在处理水体中重金属Cu2+的应用

利用聚丙烯腈基活性碳纤维(PAN-ACF)表面存在的羧基与氯化亚砜SOCl2发生酰氯化反应后与乙二胺110℃加热回流在其表面接枝胺基、酰胺基团对其进行改性,用傅立叶转换红外光谱仪(FTIR)、X射线光电子能谱仪(XPS)对改性前后活性碳纤维的表面结构和元素结合态及含量变化进行了表征,并将改性后的活性碳纤维(NH-ACF)应用于水体中Cu2+离子的去除处理,用电感耦合等离子体原子发射光谱仪(ICP-AES)检测被处理溶液中的Cu2+离子浓度变化,进而得出Cu2+离子在NH-ACF上的最大吸附量qmax.结果表明,NH-ACF比未改性处理的PAN-ACF对Cu2+吸附去除能力具有显著的优势,前者约为后者的1.5倍.