Dissipative particle dynamics simulation of wettability alternation phenomena in the chemical flooding process

Wettability alternation phenomena is considered one of the most important enhanced oil recovery （EOR） mechanisms in the chemical flooding process and induced by the adsorption of surfactant on the rock surface. These phenomena are studied by a mesoscopic method named as dissipative particle dynamics （DPD）. Both the alteration phenomena of water-wet to oil-wet and that of oil-wet to water-wet are simulated based on reasonable definition of interaction parameters between beads. The wetting hysteresis phenomenon and the process of oil-drops detachment from rock surfaces with different wettability are simulated by adding long-range external forces on the fluid particles. The simulation results show that, the oil drop is liable to spread on the oil-wetting surface and move in the form of liquid film flow, whereas it is likely to move as a whole on the waterwetting surface. There are the same phenomena occuring in wettability-alternated cases. The results also show that DPD method provides a feasible approach to the problems of seepage flow with physicochemical phenomena and can be used to study the mechanism of EOR of chemical flooding.     

高温激波管化学非平衡流动数值模拟研究

为了预测氢氧定容燃烧驱动的高温激波管性能,需要准确分析激波管非定常化学非平衡流动过程.本文在破膜前的驱动段定容燃烧以及破膜后的化学非平衡流动数值模拟中,引入双时间步长方法,发展高温激波管化学非平衡流动数值模拟方法,该方法在时间上具有二阶精度.计算结果与目前存在的激波管流动解析解以及零维化学反应系统的数值解进行了比较,吻合较好.对于典型高温激波管状态,采用有限体积方法离散准一维流动Euler控制方程,并通过将流动过程和化学反应动力学过程耦合求解,获得了激波管内部的化学非平衡流动特征.     

含化学反应膛口流场的无网格数值模拟

基于无网格方法,对包含大位移运动边界和非平衡化学反应的膛口流场进行了数值模拟。所发展算法是基于线性基函数最小二乘显式无网格方法,忽略黏性及湍流的影响,对流场采用ALE(arbitrary Lagrangian-Eulerian)形式的Euler方程描述,对流通量和化学反应源项采用多组分HLLC(Harten-Lax-van Leer-Contact)格式和有限速率反应模型计算,对于运动边界造成的点云畸形采用局部点云重构方法处理,重构过程中采用虚拟边阵面推进。对圆柱绕流和激波诱导燃烧流场进行了数值模拟,验证了重构方法和化学反应计算的有效性。最后对12.7 mm口径机枪膛口流场进行了模拟,结果同实验照片、非结构网格方法结果吻合较好,数值结果清晰地再现了膛口初始冲击波、膛口冲击波、欠膨胀射流波系结构的动力学发展过程,以及膛口焰的时间、空间分布特征。

     

Steady-state two-phase flow through planar and nonplanar model porous media

Transport in Porous Media - A comparative experimental study of ‘steady-state’ two-phase flow in two types of model porous media is made to determine the effects of nonplanarity on the...     

弹塑性变形油藏中多相渗流的数值模拟

基于流固耦合力学理论,建立了弹性变形油藏中多相渗流的数学模型,该模型考虑了渗流与变形的耦合作用,以及注采交变载荷作下油藏多孔介质的弹性变形特征,给出了耦合数值模拟方法和算例。     

Numerical simulation of silencers

The problem of attenuating the noise from weapons firing is studied experimentally and numerically. As a possible method of attenuating the noise significantly, a silencer with no internal baffles is attached to the M242 cannon. The internal pressures inside the muffler are measured. The near-field overpressures outside the muffler at various polar angles are also measured. A numerical simulation of the flow through the muffler is performed, using Harten's shock-capturing method to solve the Euler equations of ideal compressible flow. The numerical simulation yields a detailed picture of the flow field as displayed by the pressure and Mach contours. Pressure–time curves at selected locations are obtained and compared with experimental data. There is good agreement, except that the numerical simulation generates more vigorous oscillations.     

NUMERICAL SIMULATION OF COMPLICATED PIPE FLOW WITH k-εMODEL

In this paper, turbulence in a complicated pipe is simulated by using the k-ε model. The ladder-like mesh approximation is used to solve the problem of complicated boundary with the result of numerical simulation favorable. Two computational examples are given to validate the strong adaptability and stability of k-ε model.     

垂直上升气液环状流流动参数的数值计算

提出了一个新的气核-液膜耦合模型来求解垂直上升气液环状流在充分发展段的流动参数.本模型考虑了液膜、气核以及它们之间的相互影响和作用.模型中基本的气核区域和液膜区域的质量和动量方程由Fluent6.3.26进行求解,而液滴方程以及相界面上的夹带和沉积作用通过用户自定义接口函数UDF(User Defined Functi...     

NUMERICAL SIMULATION OF THREE DIMENSIONAL GAS-PARTICLE FLOW IN A SPIRAL CYCLONE

The three-dimension gas-particle flow in a spiral cyclone is simulated numerically in this paper. The gas flow field was obtained by solving the three-dimension Navier-Stokes equations with Reynolds Stress Model (RSM). It is shown that there are two regions in the cyclone, the steadily tangential flow in the spiral channel and the combined vortex flow in the centre. Numerical results for particles trajectories show that the initial position of the particle at the inlet plane substantially affects its trajectory in the cyclone. The particle collection efficiency curves at different inlet velocities were obtained and the effects of inlet flow rate on the performance of the spiral cyclone were presented. Numerical results also show that the increase of flow rate leads to the increase of particles collection efficiency, but the pressure drop increases sharply.     

Numerical investigations on dynamic process of muzzle flow

The integrative process of a quiescent projectile accelerated by high-pressure gas to shoot out at a supersonic speed and beyond the range of a precursor flow field was simulated numerically. The calculation was based on ALE equations and a second-order precision Roe method that adopted chimera grids and a dynamic mesh. From the predicted results, the coupling and interaction among the precursor flow field, propellant gas flow field and high-speed projectile were discussed in detail. The shock-vortex interaction, shockwave reflection, shock-projectile interaction with shock diffraction, and shock focus were clearly demonstrated to explain the effect on the acceleration of the projectile.     

Numerical simulation of the flow field in a thermal plasma reactor

A numerical simulation is presented for a thermal plasma reactor with particle-trajectory model in this paper.Turbulance is considered by using simple SGS model.Thegoverning equations are solved by means of the algorithm of SIMPLER.The calculatedresults give the velocity and the temperature fields within plasma reactor,and thetrajectories of the injected particles.     

A dissipative particle dynamics study of flow in periodically grooved nanochannels

The effect of periodic rectangular wall roughness on planar nanochannel flow is investigated by dissipative particle dynamics simulation. The wall protrusion length is varied, and its effect on the flow is examined. Analysis of particle trajectories and average residence time reveals temporary trapping of fluid particles inside the rectangular cavities for a considerable amount of time. This trapping affects the density, velocity, pressure, and temperature distribution inside and close to the cavities. Inside the cavities, low‐velocity regions and regions of high density related to high pressure and high temperature are observed. When compared with that of the channel with flat walls case, lower flow velocities, temperatures, and pressures are observed for grooved channels. The reduction of the above quantities is more pronounced as the protrusion length, that is, the roughness characteristic length, decreases. Finally, the relation of friction factor, f, with the flow Reynolds number is discussed. The model predicts = constant in the range . The results of this work are of direct relevance to the design of nanofluidic devices. Copyright © 2011 John Wiley & Sons, Ltd.     

二次燃烧对底排装置尾部流场影响的数值模拟

为研究二次燃烧对底排尾部流场的影响,建立了底排装置尾部流场的化学非平衡流数学物理模型。其中二次燃烧模型采用10组分25步反应的H₂-CO燃烧模型,运用统一算法的思路编程求解二维轴对称方程组。数值模拟结果与实验结果较吻合。在此基础上,对尾部流场以及燃烧特性进行了数值预测, 结果表明:二次燃烧所释放的热能远大于排气本身的热能,对增压减阻的贡献可达78%。二次燃烧改变了尾部的温度分布规律,使温度峰值分布在两个回流区内。排气沿着两回流区间的狭缝流入剪切层发生燃烧。一部分混气回流入底部附近,其中氧气不充足,存在大量CO和少量H₂未直接反应。一部分混气沿着剪切层流入下游以及主回流区内,氧含量逐渐增多,H₂和CO被反应殆尽。结果可为进一步研究底排增压减阻提供参考。     

三维溃坝流的光滑粒子动力学方法模拟

采用光滑粒子动力学SPH(Smoothed Particle Hydrodynamics)方法对三维溃坝流问题进行了数值模拟。为了逼真地模拟出坝内水体与壁面间相互作用而产生的水花飞溅、融合以及近壁面流动等现象,加入了混合长度形式的湍流模型。为了有效地防止粒子穿透固壁,提出了一种新型的适合三维数值模拟的固壁边界处理方法。应用SPH方法对三维溃坝流进行了数值模拟,并分别考虑了未添加障碍物和添加圆柱障碍物两种情形。计算结果表明,改进SPH方法能够精细地捕捉溃坝流在不同时刻的自由液面,并获得稳定而精确的数值结果。     

空气炮碰撞实验台发射过程的数值模拟

利用计算流体动力学方法对单级空气炮碰撞实验台的发射过程进行了数值仿真,所得碰撞车发射速度与实验结果吻合。在此基础上,对碰撞车发射过程的流场变化、车体前后压力以及储气罐压力变化进行分析。结果表明,泄漏气体先于碰撞车充满整个发射空间,形成初始流场,使得车前压力出现正负交替现象,但其数值较小,对车体加速过程的影响可以忽略。当碰撞车进入泄压段后,受冲击射流作用,碰撞车仍处于加速状态,且速度增量约为2 m/s。

     

NUMERICAL SIMULATION OF FLOW OF HIGHLY VISCOELASTIC FLOW IN THREE－DIMENSIONAL VARYING THICK SLIT CHANNEL

（杨伯源）（李勇）ＮＵＭＥＲＩＣＡＬＳＩＭＵＬＡＴＩＯＮＯＦＦＬＯＷＯＦＨＩＧＨＬＹＶＩＳＣＯＥＬＡＳＴＩＣＦＬＯＷＩＮＴＨＲＥＥ－ＤＩＭＥＮＳＩＯＮＡＬＶＡＲＹＩＮＧＴＨＩＣＫＳＬＩＴＣＨＡＮＮＥＬ￥ＹａｎｇＢｏｙｕａｎ；ＬｉＹｏｎｇ（Ｄｅｐａｒｔ...     

Numerical simulation and homogenization of two-phase flow in heterogeneous porous media

A mathematically rigorous method of homogenization is presented and used to analyze the equivalent behavior of transient flow of two incompressible fluids through heterogeneous media. Asymptotic expansions and H-convergence lead to the definition of a global or effective model of an equivalent homogeneous reservoir. Numerical computations to obtain the homogenized coefficients of the entire reservoir have been carried out via a finite element method. Numerical experiments involving the simulation of incompressible two-phase flow have been performed for each heterogeneous medium and for the homogenized medium as well as for other averaging methods. The results of the simulations are compared in terms of the transient saturation contours, production curves, and pressure distributions. Results obtained from the simulations with the homogenization method presented show good agreement with the heterogeneous simulations.     

非结构混合网格化学非平衡粘性绕流数值模拟

在非结构混合网格上对化学非平衡粘性绕流进行了数值模拟。控制方程为考虑了化学非平衡效应的二维Navier-Stokes方程,化学动力学模型为7组元、7反应模型。控制方程中的对流项采用VanLeer逆风分裂格式处理,并应用MUSCL方法及Minmod限制器扩展到二阶精度,粘性项用中心差分格式处理。时间推进采用显式5步龙格-库塔方法。为了适应高超声速流场计算,对VanLeer通量分裂方法进行了改进,并引入了化学反应时间步长。对RAMC-II模型的飞行试验流场进行了数值模拟,计算结果与试验测量数据符合较好,并与参考文献中的数值模拟结果吻合。     

高超声速非定常流动的数值模拟与气动热计算

高超声速飞行器研究中的一个重点问题是飞行器表面的气动加热,它对飞行器的气动、热特性及安全性有重要的影响.受到当前实验技术的限制,地面实验无法准确模拟真实飞行条件,所以采用数值模拟研究气动加热问题成为目前重要的研究手段.本文采用数值方法求解三维N-S方程,得到钝头体再入模型绕流的瞬态流场,驻点温度及表面热流沿轨道变化规律.计算中采用变边界条件模拟沿轨道飞行的非定常性.