基于节点分离Lagrange有限元方法的超高速碰撞碎片云数值模拟

提出了一种基于节点分离概念的Lagrange有限元冲击破碎分析的算法,采用节点分离技术实现 了网格断裂,并用畸变侵蚀技术处理严重畸变单元。利用C++编程实现了节点分离计算模型的创建和畸变 侵蚀,结合LS-dyna的重启动分析功能,编程调用LS-dyna求解器以及畸变侵蚀程序实现了破碎分析。分别 用节点分离的Lagrange有限元算法、SPH 无网格算法对超高速碰撞问题进行了数值模拟,并与文献中的实 验结果和Euler有限元数值模拟结果进行了对比。结果表明,节点分离算法具有计算速度快、数值稳定、边界 明确等优点,能准确有效模拟超高速碰撞问题。     

热喷干扰气体模型对飞行器气动特性影响分析

针对不同气体模型对高超声速飞行器喷流反作用控制系统(RCS)热喷干扰流场模拟的计算效率和准确性问题, 基于喷流燃气物理化学模型, 通过数值求解含化学反应源项的三维N-S方程, 建立了飞行器RCS热喷干扰流场数值模拟方法, 分别采用化学反应流、反应冻结流、二元异质流以及空气喷流四种气体模型开展了典型外形热喷干扰流场的数值模拟, 研究了不同气体模型对热喷干扰流场结构、飞行器气动力热特性的影响, 分析了不同马赫数、飞行高度下的变化规律. 研究表明: 化学反应流模型计算精度较高, 计算与风洞试验数据的吻合程度优于其他三种简化模型; 在本文的低空条件下, 采用简化模型进行热喷干扰流场数值模拟, 会低估分离区大小, 使飞行器气动力特性预测出现偏差, 同时也会低估表面热环境, 对防热系统设计不利, 随着马赫数增加, 简化模型对气动力热特性预估的误差进一步增大, 同时不同简化模型之间的差异也进一步增大; 飞行高度较高时, 模型之间的差异减小, 此时可采用简化模型进行计算以提高计算效率. 本文的研究结果可为飞行器热喷干扰流场数值模拟及喷流反作用控制系统设计提供参考.      

一种冲击波作用下结构毁伤算法研究

针对爆炸冲击波与建筑物结构相互作用过程,分析了冲击波与结构碎块作用机理,发展了一种能够模拟建筑物结构破坏及冲击波传播过程的计算模型和方法。采用建筑物结构工程毁伤载荷作为判据,处理结构在冲击波作用下的破坏问题;利用流固耦合界面算法处理结构运动引起的泄压效应,利用“虚拟网格通气技术”处理结构碎块对冲击波的阻碍作用,模拟了冲击波作用下典型建筑物的毁伤过程及冲击波传播过程。结果表明,该模型在模拟冲击波与结构的作用过程中,压力计算结果与非结构动网格模拟结果符合较好;在典型建筑物毁伤过程的数值模拟中,计算得到的建筑物毁伤效果和冲击波超压分布与建筑物物理毁伤特点符合。     

高超声速尾迹流场稳定性数值研究

通过数值模拟, 对高超声速尾迹流场进行了研究, 对其尾迹流动的失稳过程进行了分析.选取计算模型为圆球,Ma= 6.0, Re = 1.71\times 10^6(Re以球头半径为参考长度). 通过数值模拟,首先得到的流动是稳定解,在底部发展出一个主分离区和一个二次分离区,流动是轴对称状态. 不添加任何扰动继续进行计算,发现底部流场缓慢发展出微弱的非定常流动. 随后,该现象继续发展,出现明显的结构失稳,得到了无量纲周期为12.0的周期解. 给出了高超声速圆球绕流尾迹结构的周期性演化过程,对其涡系结构的演化及奇点特征进行了分析. 研究表明该数值模拟方法可用于底部流动稳定性问题的研究,同时证实了高超声速底部流动也存在流动不稳定性.      

基于滑移速度壁模型的复杂边界湍流大涡模拟

采用滑移速度壁模型实现了浸入边界方法与壁模型相结合的大涡模拟.本文首先分别采用平衡层模型和非平衡壁模型对周期山状流进行数值模拟,以考查在壁模型中考虑切向压力梯度的作用.数值结果表明,流场的压力对本文所采用的壁模型形式并不敏感,但是考虑切向压力梯度可以显著改进壁面摩擦力的计算结果,并且能够准确的预测强压力梯度区以及分离区内的流动平均统计特性.不考虑压力梯度效应的平衡层模型显著低估了壁面摩擦力的分布,同时无法准确预测分离区内的平均速度剖面.非平衡模型的修正项正比于切向压力梯度和壁面法向距离,因此在强压力梯度区或者网格较粗时,计算得到的平均压力和摩擦力分布以及流动的低阶统计量均与参考的实验和计算结果吻合.在此基础上,通过回转体绕流的大涡模拟考查了该方法用于模拟高雷诺数壁湍流的适用性,非平衡壁模型可以准确地捕捉流动的物理结构并较准确地预测其水动力学特性.结果表明,将浸入边界方法与非平衡滑移速度壁模型相结合的大涡模拟,有望成为数值模拟复杂边界高雷诺数壁湍流的工具.      

关于交通流中扰动传播和发展的数值模拟

针对不同拥挤程度交通流中扰动的传播和发展，进行了系统的数值模拟研究. 讨论了PW模型、一维管流模型、单侧传播模型等几种典型交通流模型的不同特征；指出了 不同差分格式对数值结果的影响；也分析了数学模型中不同形式的平衡函数的作用. 结果发 现在适当的数学模型、平衡函数和离散格式下，能够对交通流中扰动的传播和发展，特别是 扰动波的传播速度，取得与实际测量数据相当接近的模拟效果. 建议了分别适用于不同密度 交通流的计算格式.     

基于滑移速度壁模型的复杂边界湍流大涡模拟

采用滑移速度壁模型实现了浸入边界方法与壁模型相结合的大涡模拟.本文首先分别采用平衡层模型和非平衡壁模型对周期山状流进行数值模拟,以考查在壁模型中考虑切向压力梯度的作用.数值结果表明,流场的压力对本文所采用的壁模型形式并不敏感,但是考虑切向压力梯度可以显著改进壁面摩擦力的计算结果,并且能够准确的预测强压力梯度区以及分离区内的流动平均统计特性.不考虑压力梯度效应的平衡层模型显著低估了壁面摩擦力的分布,同时无法准确预测分离区内的平均速度剖面.非平衡模型的修正项正比于切向压力梯度和壁面法向距离,因此在强压力梯度区或者网格较粗时,计算得到的平均压力和摩擦力分布以及流动的低阶统计量均与参考的实验和计算结果吻合.在此基础上,通过回转体绕流的大涡模拟考查了该方法用于模拟高雷诺数壁湍流的适用性,非平衡壁模型可以准确地捕捉流动的物理结构并较准确地预测其水动力学特性.结果表明,将浸入边界方法与非平衡滑移速度壁模型相结合的大涡模拟,有望成为数值模拟复杂边界高雷诺数壁湍流的工具.     

冲击载荷作用下岩体破坏规律的数值流形方法模拟研究

为了更好地模拟岩体这种典型的连续与非连续介质共存的复杂结构体在冲击载荷作用下的破坏规律,在近年来新出现的数值流形方法基本理论的基础上,利用其算法程序对岩体中存在的不连续面如节理、裂隙等对岩体破裂效果的影响进行了数值模拟分析,计算结果充分显示了数值流形方法在模拟这类问题中的准确性和有效性。     

带尾翼的翻转型爆炸成形弹丸的三维数值模拟

利用LS-DYNA三维动力有限元程序对三点起爆的翻转型EFP形成过程及三个尾翼进行了数值模拟。计算结果与实验结果的比较表明:数值计算结果和实验结果吻合较好,可为弹丸优化设计提供重要参考依据。对药型罩材料分别选择Johnson-Cook和Steinberg两种本构模型进行了数值模拟。计算结果表明:药型罩材料的本构模型选取对形成的尾翼效果有一定影响,Steinberg本构模型与实验结果符合更好。     

满足几何守恒律的WENO格式及其应用

对几何守恒律的来源进行了分析,发展了一种满足几何守恒律的WENO格式,并应用于翼型层流分离现象的数值模拟中。为消除网格质量影响,采用守恒型方法计算网格导数,并将标准的WENO格式分解为中心差分部分和数值耗散部分。算例计算结果表明,几何守恒律对高精度有限差分WENO格式计算结果具有重要影响,本文方法能够消除网格导数计算误差,保证来流保持性。将本文方法应用于SD7003翼型层流分离现象的数值模拟中,计算结果与文献中计算及试验数据吻合较好,同时能够精细捕捉小尺度流场结构,准确模拟翼型层流分离现象中的复杂流动过程。     

油水重力分离实验模拟系统的研制

所研制的油水重力分离实验模拟系统，由介质循环、乳化液制备、激光片光源、图象采集及设备模型等部分组成，可对系统流量、介质物性、分散相含量、乳化状态、设备结构等多种因素进行实验模拟，从而能有效地将实际介质在工程设备中的流动与分离过程再现到实验室中．该系统与ＰＩＶ、停留时间分布等有关检测技术相配合，可对多种油田地面分离设备进行多工况、多参数的流动诊断分析，为揭示介质在设备中的流动与分离规律、研究设备工作理论、优化设备结构等提供了重要手段．     

NUMERICAL SIMULATION OF FINE PARTICLE SEPARATION IN A ROTATIONAL TUBE SEPARATOR

This paper presents a numerical analysis of gas-solid separation in a rotational tube separator. This separator which collects fine particles from gas in laminar flow is effective for fine particle separation. The separation efficiency and critical particle diameter of the separator were simulated using CFD package （FLUENT 6.0）. The simulation showed that separation efficiency can be significantly decreased due to the presence of turbulence. The simulation also showed that the Saffman lift force has little effect on the efficiency of this separator. The critical particle diameter of this tube separator was also calculated theoretically, Some experimental data were provided to validate the simulation results. Comparison between experimental results and simulation predictions on separation efficiency showed satisfactory agreement.     

On the effect of the nozzle design on the performances of gas–liquid cylindrical cyclone separators

This paper constitutes an experimental study of the separation performances of a gas–liquid cylindrical cyclone (GLCC) separator that interests the oil industry. The global hydrodynamics behavior in the GLCC is characterized by flow visualization under various inflow operating conditions. The effect of the inlet nozzle design on the performances of the separator is studied by using three different nozzles, and it proves to be a key parameter. With an insufficient nozzle restriction, low swirl intensity is imparted to the flow. Due to inadequate centrifugal effects, liquid is prematurely carried over by the gas as flooding occurs in the separator upper part. High amounts of gas are also carried under by the liquid stream. On the other hand, with a too severe nozzle convergence, the important drag applied by the gas leads to liquid “short circuiting” the cyclone toward the gas outlet. In addition to the nozzle design, the separator performances are influenced by phenomena such as liquid bridging or the occurrence of the slug flow regime at the cyclone inlet. This paper leads to a better understanding of the links between the hydrodynamics in the GLCC and its operational limits, which is necessary to enable reliable scaling up tools.     

Upwind finite difference method for miscible oil and water displacement problem with moving boundary values

The research of the miscible oil and water displacement problem with moving boundary values is of great value to the history of oil-gas transport and accumulation in the basin evolution as well as to the rational evaluation in prospecting and exploiting oil-gas resources. The mathematical model can be described as a coupled system of nonlinear partial differential equations with moving boundary values. For the twodimensional bounded region, the upwind finite difference schemes are proposed. Some techniques, such as the calculus of variations, the change of variables, and the theory of a priori estimates, are used. The optimal orderl2-norm estimates are derived for the errors in the approximate solutions. The research is important both theoretically and practically for the model analysis in the field, the model numerical method, and the software development.     

Oil recovery in the presence of microbubbles in the filtration flow

This paper presents mathematical models for oil-gas flow taking into account the various processes due to the formation of gas micronuclei (microbubbles) in oil: slip of oil relative to the walls of the pore channels (gas lubrication), changes in oil viscosity, and motion of microbubles with respect to oil. We consider examples of oil flow in the near-wellbore zone for the case where a reduction in pressure to the saturation pressure leads to the formation of gas microbubbles and micronuclei and examples of the action of a water-gas mixture in the case where oil foams in the contact area with the injected gas, i.e., a finely dispersed mixture of oil and microbubbles is formed. The behavior of indicator curves for an oil well with the formation of microbubbles is simulated, and the effect of microbubbles on the oil recovery factor in a water-alternating-gas injection process is studied.     

静电旋风分离器气相流场的数值模拟

对一内部安装电晕极的切向进口旋风分离器,以三维贴体坐标为基础,应用Bradshaw的修正k-ε湍流模型,用非交错的SIMPLE算法对静电旋风分离器的气相流场进行求解,计算结果与实验结果进行了对比。分析了电晕极的不同安装位置对旋风分离器流场的影响。从流场的角度来看,电晕极安装在筒体与排气管之间并靠近排气管的位置有利于提高静电旋风分离器的分离效率。     

可降解性表面活性剂体系下水平管内气液两相螺旋流实验研究

通过气液两相螺旋流实验仪器,研究具有可降解性的天然椰子油新型添加剂对于气液两相螺旋流流型影响以及流型的转变规律,并与表面活性剂十二烷基苯磺酸钠(SDBS)进行对比研究。实验工况设定为:实验介质为空气和水,含气率10%~90%,气相折算速度0.01~4.0m/s,液相折算速度0.01~4.0m/s,表面活性剂采用从植物提取的可降解性椰子油和SDBS,起旋装置为叶轮。实验观察到天然椰子油对于螺旋轴状流、螺旋团状流、螺旋弥散流转换特性的影响与SDBS的效果相类似,该三种流型发生条件相比于以往都有所提前,且存在范围被拓宽。浓度为500ppm时椰子油体系下的主要流型为螺旋弥散流,而SDBS体系下则以螺旋团状流为主。     

Simulation of immiscible multiphase flow in porous media: A focus on the capillary fringe of oil-contaminated aquifers

This paper deals with the analysis of some aspects of the vertical and lateral migration of oil spills in the unsaturated and the capillary zone of a phreatic aquifer. Our motivation stems from the fact that such contamination represents a severe danger for ground-water resources all over the world and from the present acute problem of jet-fuel contamination in some location of Israel. In the present study, we shall focus our efforts on the analysis of the upper layers of the aquifer which are often subjected to the most significant oil contamination. Neglecting coupled processes effects such as dilution, adsorption and volatilization, also adopting Richard's assumption, a three-phase flow model is introduced with capillary heads of the water and the oil as variables. The resulting model which is coupled and strongly non-linear is solved using a vertical two-dimensional Finite-Element procedure together with a quasi-Newton optimization algorithm. Applying that scheme, various scenarios of oil migration in the unsaturated and the capillary zone were simulated. Some migration characteristics prediced by the numerical simulations are discussed. In particular, the dynamics of the water and oil phases during the migration process is discussed.     

水平井变质量流研究进展

综述了水平井筒沿管壁具有流体流入的变质量流动的最新研究进展.内容包括:考虑油藏流体流入时井筒流动压力及产量的计算,具有不同流入时管壁摩擦系数的经验计算公式,以及流体径向流入井筒后流动阻力损失产生的机理分析等.重点介绍了水平井筒变质量流动阻力系数的计算,同时文中还介绍了井筒两相变质量流动损失产生的机理和数值模拟研究进展.      

Characterizing the Role of Shale Geometry and Connate Water Saturation on Performance of Polymer Flooding in Heavy Oil Reservoirs: Experimental Observations and Numerical Simulations

Many heavy oil reservoirs contain discontinuous shales which act as barriers or baffles to flow. However, there is a lack of fundamental understanding about how the shale geometrical characteristics affect the reservoir performance, especially during polymer flooding of heavy oils. In this study, a series of polymer injection processes have been performed on five-spot glass micromodels with different shale geometrical characteristics that are initially saturated with the heavy oil. The available geological characteristics from one of the Iranian oilfields were considered for the construction of the flow patterns by using a controlled-laser technology. Oil recoveries as a function of pore volumes of injected fluid were determined from analysis of continuously recorded images during the experiments. We observed a clear bypassing of displacing fluid which results in premature breakthrough of injected fluid due to the shale streaks. Moreover, the results showed a decrease of oil recovery when shales’ orientation, length, spacing, distance of the shale from production well, and density of shales increased. In contrast, an increase of shale discontinuity or distance of the shale streak from the injection well increased oil recovery. The obtained experimental data have also been used for developing and validating a numerical model where good matching performance has been observed between our experimental observations and simulation results. Finally, the role of connate water saturation during polymer flooding in systems containing flow barriers has been illustrated using pore level visualizations. The microscopic observations confirmed that besides the effect of shale streaks as heterogeneity in porous medium, when connate water is present, the trapped water demonstrates another source of disturbance and causes additional perturbations to the displacement interface leading to more irregular fingering patterns especially behind the shale streaks and also causes a reduction of ultimate oil recovery. This study reveals the application of glass micromodel experiments for studying the effects of barriers on oil recovery and flow patterns during EOR processes and also may provide a set of benchmark data for recovery of oil by immiscible polymer flood around discontinuous shales.