Parallel arithmetic numerical simulation and application of secondary migration-accumulation of oil resources

From such actual conditions as the effects of characteristics of miltilayer petroleum geology and permeation fluid mechanics, a new numerical model is put forward and coupling splitting-up implicit interactive scheme is formed. For the actual situation of Dongying hollow ( four-layer ) and Tanhai region ( three-layer ) of Shengli Petroleum Field, the numerical simulation test results and the actual conditions are coincident.     

Numerical simulation analysis for migration-accumulation of oil and water

ＩｎｔｒｏｄｕｃｔｉｏｎＴｈｅｎｕｍｅｒｉｃａｌｓｉｍｕｌａｔｉｏｎｏｆｔｈｅｅｖｏｌｕｔｉｏｎａｒｙｈｉｓｔｏｒｙｏｆｔｈｅｂａｓｉｎｉｓｃｏｎｄｕｃｔｅｄａｃｏｒｄｉｎｇｔｏｔｈｅｍｅｃｈａｎｉｓｍｏｆｐｅｔｒｏｌｅｕｍｇｅｏｌｏｇｙａｎｄｍｅｃ...     

NUMERICAL SIMULATION OF OIL MIGRATION-ACCUMULATION OF MULTILAYER AND ITS APPLICATION

ＩｎｔｒｏｄｕｃｔｉｏｎＴｈｅｏｉｌｆｏｒｍａｔｉｏｎｉｎｓｅｄｉｍｅｎｔｂａｓｉｎｓ ,ｉｔｓｄｉｓｐｌａｃｅｍｅｎｔ,ｔｒａｎｓｐｏｒｔａｎｄａｃｃｕｍｕｌａｔｉｏｎ ,ａｎｄｔｈｅｆｉｎａｌｆｏｒｍａｔｉｏｎｏｆｏｉｌｄｅｐｏｓｉｔｓｈａｖｅｂｅｅｎｏｎｅｏｆｔｈｅｋｅｙｐｒｏｂｌｅｍｓｉｎｔｈｅｅｘｐｌｏｒａｔｉｏｎｏｆｏｉｌ＿ｇａｓｒｅｓｏｕｒｃｅｓ.Ｈｏｗｈａｓｏｉｌｂｅｅｎａｃｃｕｍｕｌａｔｅｄｉｎｔｈｅｐｒｅｓｅｎｔｌｏｏｐａｃｃｏｒｄｉｎｇｔｏｔｈｅｍｅｃｈａｎｉｃｓｏｆｉｍｍｉｓｃ…     

泥水平衡盾构机掘进并行数值模拟方法研究

以某大型越江隧道工程为背景,建立该工程场地土层的大规模三维非线性数值模型。同时基于共享内存式并行计算平台(SMP),利用刚度迁移法,完成了泥水平衡盾构机逐步掘进百米长距离,并穿越大堤这一复杂力学过程的大规模数值模拟。为模拟泥水平衡盾构机掘进各施工阶段,不同施工因素之间的复杂力学作用关系,提出了对各施工因素模拟的简化方法。具体包括:开挖面泥水压力对前方土体的作用、盾构机机头的超挖、盾构机机身对土体的作用、盾尾注浆材料的时空变化特性、盾构机后方百米长拖车对已铺成隧道的作用等施工因素。通过对计算结果与工程试验段实测数据对比验证的基础上,预测了地表沉降在盾构机掘进过程中的变化规律。通过对SMP并行计算平台下盾构法施工数值模型之并行计算效率的对比研究,给出了各并行算法的加速比、计算效率以及接触模型对并行计算效率的影响关系。     

A model for simulation of friction phenomenon between dies and workpiece

Based on the interaction of asperities and upperbound approach a mathematical model for simulation of friction phenomenon between dies and workpiece is proposed. Optimizing the mathematical model with respect to several variables it is found that in addition to adhering, tearing, ploughing, etc., asperities workpiece can move wave-like along the surface layer and under certain circumstances they may disappear. If the asperities wavily move along the surface layer the friction coefficient depends on the geometry of asperities. However, the bonding strength of asperities has no significant influence on friction coefficient. The depth of the plastic deformation layer is related to the geometry of asperities, too. The soundness of the prerequisite of the proposed model and some analytical results were verified by experiments.Department of Materials Engineering, Northwestern Polytechnical UniversityThis is a project of the cooperative program between the Department of Materials Engineering of Northwestern Polytechnical University of the People's Republic of China and the Institute for Metal Forming at Lehigh University of the United States of America.     

建筑风压数值模拟的几种并行化策略

针对四个处理机的Ｔｒａｎｓｐｕｔｅｒ并行计算机系统，建立了建筑风压数值模拟问题基于ＳＩＭ－ＰＬＥＣ算法的几种并行化策略：分区并行策略、方程并行策略和双重并行策略。对各种策略的计算流程、数据通讯及并行效率等进行了分析和比较，并通过实例计算作了验证。     

无钴合金钢的层裂断裂及数值模拟研究

利用一级轻气炮作为加载手段,对无钴合金钢的层裂特性进行了研究,获得了Hugoniot弹性极限、层裂强度、层裂片厚度以及塑性应变率等动力学参数。对回收的样品进行了断口分析和金相分析,从宏微观角度分析了无钴合金在不同应变率下的断裂机制。利用LS-DYNA有限元程序对层裂现象进行了数值模拟,计算结果与实验结果基本吻合。     

岩石弹塑性破裂过程的数值模拟研究

分析了对于非均质、非连续且材料参数分布具有分散性和随机性的岩体介质,采用传统有限元法解决其工程问题的局限性．阐述了用岩石力学性质的分布特性来描述岩石材料不均匀性,并用蒙特卡洛模拟方法,将符合Weibull分布规律的材料参数赋予岩体结构中的各个单元体,进而用计算机进行数值模拟的研究方法．应用了在Windows95环境下利用FORTRAN Power Station4．0平台开发的具有可视性的“岩石弹塑性破裂过程分析”软件——REPFPA软件,对典型的岩石试样和巷道模型的弹塑性破裂过程进行数值模拟研究,得到了有价值的实验结果．证明把材料参数随机赋值方法引入到岩土类材料的弹塑性破裂过程分析中去,并用计算机进行数值模拟是研究岩石弹塑性破裂过程的一种有效方法。     

船舶密闭舱室火灾过程及其模拟计算

密闭舱室火灾是一个三维非定常过程，压力升高是这一过程的重要特征。本文采用质量守恒的方法进行压力校正，建立了计算密闭舱室火灾过程的数学物理模型，数值计算采用非均匀交错网格，结果与实验和理论分析相符较好。     

Study of numerical errors in direct numerical simulation and large eddy simulation

By comparing the energy spectrum and total kinetic energy, the effects of numerical errors (which arise from aliasing and discretization errors), subgrid-scale (SGS) models, and their interactions on direct numerical simulation (DNS) and large eddy simulation (LES) are investigated. The decaying isotropic turbulence is chosen as the test case. To simulate complex geometries, both the spectral method and Pade compact difference schemes are studied. The truncated Navier-Stokes (TNS) equation model with Pade discrete filter is adopted as the SGS model. It is found that the discretization error plays a key role in DNS. Low order difference schemes may be unsuitable. However, for LES, it is found that the SGS model can represent the effect of small scales to large scales and dump the numerical errors. Therefore, reasonable results can also be obtained with a low order discretization scheme.     

Direct numerical simulation of turbulent pipe flow

Fully developed turbulent pipe flow at low Re-number is studied by means of direct numerical simulation (DNS). In contrast to many previous DNS's of turbulent flows in rectangular geometries, the present DNS code, developed for a cylindrical geometry, is based on the finite volume technique rather than being based on a spectral method. The statistical results are compared with experimental data obtained with two different experimental techniques. The agreement between numerical and experimental results is found to be good which indicates that the present DNS code is suited for this kind of numerical simulations.     

薄膜结构流固耦合的CFD数值模拟研究

基于弱耦合分区求解策略,在CompaqVisualFortran6．5环境下搭建了薄膜结构三维流固耦合效应的CFD数值模拟平台。程序采用模块化编程思想,主要包含几何建模、流体分析、结构分析和数据交换四个模块。其中几何建模模块采用自行编制的膜结构找形分析程序,流体分析模块采用经过二次开发的计算流体力学软件FLUENT6．0,结构分析模块采用自行编制的膜结构动力分析程序MDLFX;在数据交换模块中,编制了基于薄板样条法的插值计算程序,以实现流固交界面上不同区域网格间的数据传递问题,编制了基于代数法和迭代法的动网格变形程序,以实现流固耦合运算中的动网格更新。基于该软件平台,对单向柔性屋盖和鞍形膜结构屋盖进行了流固耦合数值模拟,验证了方法的有效性。     

入水问题的Euler-Lagrange流固耦合数值模以技术研究

针对结构入水过程的特点,建立了一种把Eulerian计算过程和Lagrangian计算过程耦合起来的数值技术,编写了计算程序,把结构对水作用的过程和结构动力学响应过程的计算结合起来,实现了结构与水之间流固耦合作用过程的模拟.其中水、空气的运动利用Eulerian过程的流体弹塑性有限差分计算程序进行模拟,结构的运动变形利用Lagrangian过程的结构动力学有限元计算程序模拟结构响应.通过板入水问题的计算结果与商业软件计算结果进行比较,验证了计算方法对流固间相互作用处理的正确性.本文计算程序可用于一般入水问题的结构响应-流场运动分析中.     

聚能装药的欧拉数值模拟

用二维有限差分欧拉程序MEPH2Y模拟了聚能装药的作用过程,包括爆轰波的形成、传播及与其他介质的相互作用,高温高压下射流（或射弹）的形成、延展、减压、断裂,射流（或射弹）对靶的侵彻及靶的成坑和动态响应等过程。介绍了程序所用的数学模型、数值方法,以及模拟的部分问题与实验结果的比较。结果表明,数值模拟结果与实验结果符合较好。     

Crack tip superblunting: experiment,theory and numerical simulation

Experimental revelation ofsuperblunting along plane strain predominant crack tip segment is reported here for modified polyproplene. As elucidated by a heuristic model of progressive circumferential cold-draw, the formation of superblunting crack tip depends critically on the ratio between the cold-draw propagation speed and the loading speed, and contributes significantly to the material toughening, especially for the improvement in impact toughness. Detailed numerical calculations are conducted based on a hyperelastic-viscoplastic and anisotropic damage constitutive model at finite deformation. The simulated results recapitulate the essential features of crack tip superblunting. The project supported by the National Natural Science Foundation of China     

Direct numerical simulation of confined swirling jets

This study investigates the Lagrangian acceleration and velocity of fluid particles in swirling flows via direct numerical simulation. The intermittency characteristics of acceleration and velocity of fluid particles are investigated at different swirl numbers and Reynolds numbers. The flatness factor and trajectory curvature are used to analyse the effect of Lagrangian intermittency. The joint probability density function of Lagrangian acceleration and turbulence intensity is shown to explain the augmentation effect of Lagrangian intermittency by the strongly swirling levels under the relatively low intensity of turbulence. In addition, the correlation between the Lagrangian acceleration and the turbulence intensity is enhanced as the swirl level increases. It shows the important effect of swirl on the motion behaviour of fluid particles in the strongly swirling flows.     

Direct numerical simulation of compressible turbulent flows

This paper reviews the authors＇ recent studies on compressible turbulence by using direct numerical simulation （DNS）,including DNS of isotropic（decaying） turbulence, turbulent mixing-layer,turbulent boundary-layer and shock/boundary-layer interaction.Turbulence statistics, compressibility effects,turbulent kinetic energy budget and coherent structures are studied based on the DNS data.The mechanism of sound source in turbulent flows is also analyzed. It shows that DNS is a powerful tool for the mechanistic study of compressible turbulence.     

Direct numerical simulation of particle interaction with ejections in turbulent channel flows

Direct numerical simulations (DNS) of incompressible turbulent channel flows coupled with Lagrangian particle tracking are performed to study the characteristics of ejections that surround solid particles. The behavior of particles in dilute turbulent channel flows, without particle collisions and without feedback of particles on the carrier fluid, is studied using high Reynolds number DNS (Re = 12,500). The results show that particles moving away from the wall are surrounded by ejections, confirming previous studies on this issue. A threshold value separating ejections with only upward moving particles is established. When normalized by the square root of the Stokes number and the square of the friction velocity, the threshold profiles follow the same qualitative trends, for all the parameters tested in this study, in the range of the experiments. When compared to suspension thresholds proposed by other studies in the Shields diagram, our simulations predict a much larger value because of the measure used to characterize the fluid and the criterion chosen to decide whether particles are influenced by the surrounding fluid. However, for intermediate particle Reynolds numbers, the threshold proposed here is in fair agreement with the theoretical criterion proposed by Bagnold (1966) [Bagnold, R., 1966. Geological Survey Professional Paper, vol. 422-1]. Nevertheless, further studies will be conducted to understand the normalization of the threshold.     

Simulation and application of three-dimensional migration-accumulation of oil resources

ＩｎｔｒｏｄｕｃｔｉｏｎＴｈｅｏｉｌｆｏｒｍａｔｉｏｎｉｎｓｅｄｉｍｅｎｔｂａｓｉｎｓ,ｉｔｓｄｉｓｐｌａｃｅｍｅｎｔ,ｔｒａｎｓｐｏｒｔａｎｄａｃｃｕｍｕｌａｔｉｏｎ,ａｎｄｔｈｅｆｉｎａｌｆｏｒｍａｔｉｏｎｏｆｏｉｌｄｅｐｏｓｉｔｓｈａｖｅｂｅｅｎｔｈｅｋｅｙｐｒｏｂｌｅｍｓｉｎｔｈｅｅｘｐｌｏｒａｔｉｏｎｏｆｏｉｌ＿ｇａｓｒｅｓｏｕｒｃｅｓ．Ｈｏｗｈａｓｏｉｌｂｅｅｎａｃｃｕｍｕｌａｔｅｄｉｎｔｈｅｐｒｅｓｅｎｔｌｏｏｐａｃｃｏｒｄｉｎｇｔｏｔｈｅｍｅｃｈａｎｉｃｓｏｆｉｍｍｉｓｃｉｂｌｅ…     

A numerical method for simulation of attached cavitation flows

A new numerical algorithm for attached cavitation flows is developed. A cavitation model is implemented in a viscous Navier–Stokes solver. The liquid–vapour interface is assumed as a free surface boundary of the computation domain. Its shape is determined with an iterative procedure to match the cavity surface to a constant pressure boundary. The pressure distribution, as well as its gradient along the wall, is taken into account in updating the cavity shape iteratively. A series of computations are performed for the cavitating flows across three kinds of headform/cylinder bodies: conic, ogival and hemispheric heads. A range of cavitation numbers is investigated for each headform/cylinder body. The obtained results are reasonable and the iterative procedure of cavity shape updating is quite stable. The superiority of the developed cavitation model and algorithm is demonstrated. Copyright © 2006 John Wiley & Sons, Ltd.