基于多孔介质模型的膜式氧合器内部流场分析

膜式氧合器内部流体运动特性对其性能有重要影响,利用计算流体力学（CFD）对氧合器模型进行流体动力学分析是预测其性能的重要方法之一。本文基于压降实验计算氧合器纤维束的黏性阻力系数,建立了各向同性多孔介质模型。采用RNGk-ε湍流模型对不同流量下氧合器内部流场进行计算,得到了血液速度、压力和壁面剪切应力分布云图。发现随着流量的增加,氧合器内部速度梯度分布形式基本保持不变,压力分布呈倾斜状态且逐渐减小,大部分压力损失位于纤维束内,其中53.3%位于氧合室,42.6%位于变温室。氧合器血液的入口及出口位置为血液损伤的高发区域。采用溶血数值预估模型计算得到了氧合器的标准溶血指数NIH。结果表明：在低流量1.65～3.00L/min下,各向同性多孔介质模型的模拟结果与实验结果基本一致,模拟数值与实验数值的偏差会随着液体流量的增加而变大;流量为1.65～6.00L/min时,标准溶血指数NIH为0.0084～0.0835g/100L,满足人体生理允许的使用范围。     

Measurement of ultralow permeability

A permeability measurement method based on perturbed pressure decay between a source and a sink communicating through a porous medium is formulated. Nonideality of gas as well as corrections due to adsorption are considered. The method allows us to infer permeability at a given rock and fluid state. The accuracy of the method may be evaluated through comparison with theoretical decay characteristics and its modal amplitudes. A quantitative evaluation of the experimental inference is therefore possible. Based on the pressure decay theory, an apparatus built in our laboratory for measuring permeability is presented. As a secondary output, data analysis also provides porosity. The lowest permeability that the present apparatus is capable of characterizing is about 0.3 nm², and with additional improvements, capability down to 0.05 nm² is anticipated. Where a steady‐state permeameter is pragmatic, excellent agreement between the two measurements is showed, validating the method. Measurements in synthetic samples also show that the method is accurate. Porosity data from a standard helium pycnometer confirm the porosity obtained from the decay method. © 2015 American Institute of Chemical Engineers AIChE J, 62: 1278–1293, 2016     

Absorption of picoliter droplets by thin porous substrates

Absorption of picoliter (pL) droplets into porous substrates is studied experimentally and numerically. In the case of pL droplets, major phenomena involved in the interaction between droplet and porous media develop at different time scales: spreading and wetting at microseconds, absorption and wicking at milliseconds, and evaporation at seconds. Therefore, one can decouple these processes to minimize the complexity of the study. A high‐speed imaging system capable of 1 million frames per second is used to visualize individual droplets impacting, spreading, and imbibing on substrates. To simulate droplet dynamics, the governing equations for flow outside and inside porous media are proposed and solved using an in‐house developed computational fluid dynamics solver. The simulation results are in good agreement with the experimental data. The effect of drop impact velocity and fluid properties on final dot shape in the porous substrates is investigated through a series of parametric numerical studies. © 2016 American Institute of Chemical Engineers AIChE J, 63: 1690–1703, 2017     

非对称陶瓷膜管渗透性能的CFD模拟研究

陶瓷膜因其化学稳定性好、机械强度大等优点得到广泛应用。计算流体力学(CFD)的快速发展使得计算模拟成为研究和优化陶瓷膜管结构性能的有效手段。为了优化非对称结构陶瓷膜管的结构和操作参数,对其渗透性能进行了CFD计算模拟。针对非对称结构陶瓷膜管的膜层和过渡层的厚度在10 μm级的特点,采用Navier-Stokes方程和Darcy定律来分别描述膜管内和膜多孔介质内的纯水流动,利用多孔介质模型描述膜管的主体支撑层,用多孔跳跃边界简化膜管的膜层和过渡层,利用Konzey-Carmen方程对膜元件各层的渗透率进行估算。计算结果与实验值吻合较好,为优化陶瓷膜管的通道结构提供了便捷的工具。     

Contact dynamic modeling of a liquid droplet between two approaching porous materials

A computational fluid dynamics model based on a finite difference solution to mass and momentum conservation equations (Navier–Stokes equations) for a liquid droplet transport between two porous or nonporous contacting surfaces (CSs) is developed. The CS dynamic (equation of motion) and the spread of the incompressible liquid available on the primary surface for transfer are coupled with the Navier–Stokes equations. The topologies of the spread dynamic between and inside both surfaces (primary and CSs) are compared with experimental data. The amount of mass being transferred into the CS, predicted by the model, is also compared to the experimental measurements. The impact of the initial velocity on the spread topology and mass transfer into the pores is addressed. © 2014 American Institute of Chemical Engineers AIChE J, 60: 2346–2353, 2014     

Characterization of Taylor vortex flow in a short liquid column

We present a study on Taylor vortex flow in the annulus between a rotating inner cylinder and a stationary outer cylinder, featured with a wide gap (radius ratio is 0.613) and a short column (aspect ratio is 5.17). A particle image velocimetry (PIV) system was used to determine the position, shape, and velocity distribution of the vortices, by which the flow was also confirmed to lie in the nonwavy Taylor vortex regime for all operating conditions explored in this study. Our results suggest that end boundary effects are important, in which the vortex number decreases with decreasing column length. For a system with an aspect ratio of 5.17, six vortices appear in the gap with their position, size, and shape varying at different Reynolds numbers. The fluid velocities show an asymmetric feature with respect to the vortex centers, while the maximum axial and radial velocities increase almost linearly with the increasing reduced Reynolds number (Re ? Re_c). In addition, computational fluid dynamics study was employed under the same conditions, and its results agree well with the PIV measurements. Overall, this study provides a quantitative understanding of the formation of Taylor vortices in a constrained space. © 2009 American Institute of Chemical Engineers AIChE J, 2009     

Direct numerical simulation of the turbulent flow in a baffled tank driven by a Rushton turbine

We present a direct numerical simulation (DNS) of the turbulent flow in a baffled tank driven by by a Rushton turbine. The DNS is compared to a Large Eddy Simulation (LES), a Reynolds Averaged Navier‐Stokes (RANS) simulation, Laser Doppler Velocimetry data, and Particle Image Velocimetry data from the literature. By Reynolds averaging the DNS‐data, we validate the turbulent viscosity hypothesis by demonstrating strong alignment between the Reynolds stress and the mean strain rate. Although the turbulent viscosity ν_T in the DNS is larger than in the RANS simulation, the turbulent viscosity parameter C_μ = ν_T?/k², is an order of magnitude smaller than the standard 0.09 value of the k‐? model. By filtering the DNS‐data, we show that the Smagorinsky constant C_S is uniformly distributed over the tank with C_S ≈ 0.1. Consequently, the dynamic Smagorisnky model does not improve the accuracy of the LES. © 2012 American Institute of Chemical Engineers AIChE J, 2012     

Transient response of two‐layer slot coating flows to periodic disturbances

Coating uniformity requirement is becoming more severe as new products come into the market. Coating processes have to be designed not only based on the steady‐state operation but also taking into account how the flow responds to ongoing disturbances on process conditions. These disturbances may lead to thickness variation on the deposited liquid layers that may be unacceptable for product performance. This study extends available transient analysis of single‐layer slot coating to determine the amplitude of the oscillation of each individual coated layer in two‐layer slot coating process in response to small periodic perturbation on different operating parameters. The predictions were obtained by solving the complete transient Navier–Stokes equations for free surface flows. The results show the most dangerous perturbations and how the deposited film thickness variations of each layer can be minimized by changing the geometry of the die lip and liquid viscosities. © 2015 American Institute of Chemical Engineers AIChE J, 61: 1699–1707, 2015     

基于详细反应机理的富燃多孔燃烧制氢的计算流体力学模拟

Computational fluid dynamics （CFD） combined with detailed chemical kinetics was employed to model the filtration combustion of a mixture of methane/air in a packed bed of uniform 3 mm diameter alumina spherical particles. The standard k-ε turbulence model and a methane oxidation mechanism with 23 species and 39 elemental reactions were used. Various equivalence ratios （1.47, 1.88, 2.12 and 2.35） were studied. The numerical results showed good agreement with the experimental data. For ultra-rich mixtures, the combustion temperature exceeds the adiabatic value by hundreds of centigrade degrees. Syngas （hydrogen and carbon monoxide） can be obtained up to a mole fraction of 23%. The numerical results also showed that the combination of CFD with detailed chemical kinetics gives good performance for modeling the pseudo-homogeneous flames of methane in porous media.     

Influence of the boundary conditions on capillary flow dynamics and liquid distribution in a porous medium

After depositing a wetting liquid onto a porous medium surface, and under the influence of the capillary pressure, the liquid is imbibed into the porous medium creating a wetted imprint. The flow within the porous medium does not cease once all the liquid is imbibed but continues as a secondary capillary flow, where the liquid flows from large pores into small pores along the liquid interface. The flow is solved using the capillary network model, and the influence of the boundary condition on the liquid distribution within the porous medium is investigated. The pores at the porous medium boundaries can be defined as open or closed pores, where an open pore is checked for the potential threshold condition for flow to take place. In contrast, the closed pore is defined as a static entity, in which the potential condition for flow to take place is never satisfied. By defining the pores at distinct porous medium boundaries as open or closed, one is able to obtain a very different liquid distribution within the porous medium. The liquid saturation profiles along the principal flow direction, ranging from constant to steadily decreasing, to the profile with a local maximum, are found numerically. It is shown that these saturation profiles are also related to the geometrical dimension that is perpendicular to the flow principal direction, and changing the boundary type from open to closed allows the liquid distribution within the porous medium to be controlled. In addition to the liquid distribution, the influence of the boundary conditions on capillary pressure and relative permeability is investigated, where both parameters are not influenced by variation of the boundary condition types. © 2011 American Institute of Chemical Engineers AIChE J, 2012     

Fate of a sessile droplet absorbed into a porous surface experiencing chemical degradation

A general‐purpose multiphase and multicomponent computer model was developed for simulation of the spread, evaporation, and chemical reaction of sessile droplet(s) in porous substrates. In the model, chemical reactions were allowed in or between any of the liquid, gas, or solid phases present. The species mass and momentum conservation equations were solved on a finite difference mesh representing the domain. These equations were marched in time using the Runge–Kutta fourth‐order method. The model's function was studied via simulation of experiments, both those performed by the authors and found in the literature. These simulations demonstrated a quantitative match to the time history of product evolution and a similar spread of liquid reactants. The model may be particularly beneficial for predicting the extent of contamination and the possible threat outcomes of those chemical agents that are harmful when introduced into the environment. © 2014 American Institute of Chemical Engineers AIChE J, 60: 2557–2565, 2014     

柱式膜组件结构的CFD优化设计

采用Euler模型与多孔介质模型对不同结构的柱式膜组件内的流体流动进行了计算模拟。研究了曝气孔数目（开孔率为1.92%保持不变）与膜组件高度对膜组件膜丝填充区域内的气液两相分布、壁面剪应力、湍流黏度以及液相速度场的影响。计算模拟数据与实验结果吻合良好。计算模拟表明：通过减小曝气孔直径,增加曝气孔数目的方式能够促进气液两相流场与液相速度场的均匀分布,以及壁面剪应力与湍流黏度的增强;增加膜组件的高度,有利于增加单支膜组件膜面积的同时充分利用曝气擦洗过程中气液两相流对膜丝壁面进行高效的气擦洗。综合考虑膜组件的安装运输、膜丝通量分布以及能耗等因素,对于直径250 mm的膜组件采用曝气孔的直径为6.32 mm,数目为30个,长度在2~2.5 m之间为最优。     

Supercritical fluid flow in porous media: modeling and simulation

The present work aims the modeling and simulation of supercritical fluid flow through porous media. This type of flow appears in several situations of interest in applied science and engineering, as the supercritical flow in porous materials employed in chromatography, supercritical extraction and petroleum reservoirs. The fluid is constituted of one pure substance, the flow is monophasic, highly compressible and isothermal. The porous media is isotropic, possibly heterogeneous, with rectangular format and the flow is two-dimensional. The heterogeneities of porous media are modeled by a simple power law, which describes the relationship between permeability and porosity. The modeling of the hydrodynamic phenomena incorporates the Darcy's law and the equation of mass conservation. Appropriated correlations are used to model, in a realistic form, the density and the viscosity of the fluid. A conservative finite-difference scheme is used in the discretization of the differential equations. The nonlinearity is treated by Newton method, together with the conjugate gradient method. The results of the simulation for pressure and mobility of supercritical and liquid propane flowing through porous media are presented, analyzed and graphically depicted.     

The effect of particle rotation on the motion and rejection of capsular particles in slit pores

The results from a two‐dimensional computational model describing the motion of capsule‐shaped particles in a slit pore under small Re conditions are reported. Average particle velocities and particle rejection coefficients were determined for capsules with aspect ratios of 2 and 4. Two different approaches were used to characterize particle rotation and hydrodynamic particle‐pore wall interactions. In one approach, all sterically allowed particle orientation angles had equal probability, i.e., infinite rotational diffusion was assumed. In the second approach, particles were allowed to freely rotate in the pore; particle orientations were dictated by hydrodynamic forces acting on the particle surface and rotational particle diffusion was neglected. Minimal lateral migration across the pore was observed for the freely rotating particles. Although particle alignment was observed for the freely rotating particles, rejections predicted from the two approaches were found to be in close agreement. © 2018 American Institute of Chemical Engineers AIChE J, 64: 2828–2836, 2018     

Formation of liquid bridges between porous matrix blocks

It is widely accepted that, in fluid flow and transport in fractured porous media, there exists some degree of block‐to‐block interaction that may lead to capillary continuity. The formation of liquid bridges causing interaction between blocks will affect oil recovery from naturally fractured reservoirs. However, the accurate modeling of the growth and detachment of liquid bridges that may cause capillary continuity between matrix blocks remains a controversial topic. In an attempt to improve our understanding of the problem, a mechanistic model is developed in this work for the formation of liquid bridges between porous blocks. The proposed model considers growth and detachment of pendant liquid droplets perpendicular to the horizontal and smooth fracture between porous matrix blocks. The liquid bridge model is then coupled with various upscaled fracture capillary pressure models to study the liquid bridge formation process. An expression is obtained that relates the commonly used fracture capillary pressure to the critical length of the liquid element. Results based on various fracture capillary pressure models reveal that the threshold Bond number is an important parameter in the formation of liquid bridges. We introduce a simple mechanistic model for the formation of liquid bridges in a horizontal fracture between two porous blocks, advancing our understanding of the two‐phase flow in fractured porous media. © 2010 American Institute of Chemical Engineers AIChE J, 2011     

Simulation of deformable preformed particle gel propagation in porous media

Preformed particle gel (PPG) treatment is a proven cost‐effective method for improving oil recovery. Although PPG system has a suspension‐like property, it has different propagation rules from the rigid particle suspension in porous media because of its good deformation property. In this study, an advanced phenomenological model of PPG propagation in porous media is presented. The model includes both PPG plugging and restarting behaviors. Log‐normal and normal distribution functions have been introduced in this model to calculate the PPG plugging probability. Power‐law equation is used to calculate the PPG restarting rate. This method can represent the commensurate relation between PPG and throat size. Then, the equations are solved numerically, using an explicit finite‐difference formulation in conjunction with a fourth‐order Runge‐Kutta method. The results match favorably with several laboratory experiments. Finally, the propagation rules and sensitivity analysis of PPG size, permeability and injection rate to propagation rules, and permeability reduction are performed. © 2017 American Institute of Chemical Engineers AIChE J, 63: 4628–4641, 2017     

When does a vessel become a pipe?

The conditions under which the transient outflow from a punctured pipeline may be approximated as that emanating from a vessel using a simplified analytically based vessel blowdown model (VBM) is investigated in this article. The above addresses the fundamental drawback of long computational run times associated with the numerically based techniques used for simulating pipeline puncture failures. The efficacy of the VBM is tested by comparison of its predictions against simulation data obtained using a validated rigorous but computationally demanding numerical technique based on the method of characteristics. The results show that the accuracy of the VBM increases with decreasing puncture/pipe diameter ratio, line pressure, and increasing pipeline length. Surprisingly, the VBM produces more accurate predictions for two‐phase mixtures when compared with permanent gases. This is found to be a consequence of the better applicability of the isothermal bulk fluid decompression assumption within the pipeline in the case of two‐phase mixtures. © 2011 American Institute of Chemical Engineers AIChE J, 2011     

The flow and displacement in porous media of fluids with yield stress

We study the mobilization and subsequent flow in a porous medium of a fluid with a yield stress, modeled as a Bingham plastic. We use single-capillary expressions for the mobilization and flow in a pore-throat, and a pore-network model that accounts for distributed yield-stress thresholds. First, we extend the statistical physics method of invasion percolation with memory, which models lattice problems with thresholds, to incorporate dynamic effects due to the viscous friction following the onset of mobilization. Macroscopic relations between the applied pressure gradient and the flow rate for single-phase flow are proposed as a function of the pore-network microstructure and the configuration of the flowing phase. Then, the algorithm is applied to model the displacement of a Bingham plastic by a Newtonian fluid in a porous medium. The results find application to a number of industrial processes including the recovery of oil from oil reservoirs and the flow of foam in porous media.     

湿法提钒浸出段搅拌反应器结构的优化

清洁提钒工艺中,熟料的湿法浸出是重要的操作单元.浸出搅拌反应器的合理设计与优化,可缩短浸矿时间以及提高浸矿效率.本文通过改变搅拌桨桨叶间的层间距、搅拌桨的安装层数以及安装导流筒等方法,对攀钢集团公司的浸出搅拌反应器进行优化和改进;并结合计算流体动力学(CFD)Fluent商业软件,分别模拟了原浸出搅拌反应器和改进后浸出搅拌反应器的宏观流场结构.结果表明:在层间距C₂为1100mm的原双层搅拌桨浸出搅拌反应器内,流体轴向速度较小,“死区”现象较严重;与原反应器相比,调整双层搅拌桨桨叶之间的层间距C₂为1800mm以及安装3层搅拌桨或导流筒,都可加强反应器内流体的轴向流动,减小“死区”范围,进而改善流场结构的均匀分布,有助于强化流体混合.