离心泵并联运行工况点的数值求解

建立了求解多台不同型号离心泵并联运行工况点的数学模型,给出了用数值计算中的“二分法”求解该模型的方法,并附有用Exce l计算的实例。     

HEGF中CO井筒流动及扩散规律研究

目前,高能气体压裂广泛的应用于低渗透油气田,可以有效清除近井地带由于钻井、射孔和各种措施造成的污染和堵塞,达到油气井增产、注水井增注的目的。但最近一些低渗透油气田却出现了频繁的CO气体中毒事件,给油田及员工带来巨大的损失。针对这一现象,运用现场数据以及物理模拟和数学模拟方法首次通过对CO气体流动、扩散等的综合研究,建立了CO井筒流动及大气扩散模型。同时,也考察了井筒压力、气油比、风速、气体泄放速率、大气稳定度等主要因素对CO气体扩散的影响。通过现场实例计算,证实模型具有较高的准确性。此项工作的完成为建立监控系统提供基础技术依据,同时对合理、有效的开发油气田以及煤层气的开发将具有重要的实际意义。     

A comparison of transformation methods for evaluating two‐dimensional weakly singular integrals

Accurate numerical evaluation of integrals arising in the boundary element method is fundamental to achieving useful results via this solution technique. In this paper, a number of techniques are considered to evaluate the weakly singular integrals which arise in the solution of Laplace's equation in three dimensions and Poisson's equation in two dimensions. Both are two‐dimensional weakly singular integrals and are evaluated using (in a product fashion) methods which have recently been used for evaluating one‐dimensional weakly singular integrals arising in the boundary element method. The methods used are based on various polynomial transformations of conventional Gaussian quadrature points where the transformation polynomial has zero Jacobian at the singular point. Methods which split the region of integration into sub‐regions are considered as well as non‐splitting methods. In particular, the newly introduced and highly accurate generalized composite subtraction of singularity and non‐linear transformation approach (GSSNT) is applied to various two‐dimensional weakly singular integrals. A study of the different methods reveals complex relationships between transformation orders, position of the singular point, integration kernel and basis function. It is concluded that the GSSNT method gives the best overall results for the two‐dimensional weakly singular integrals studied. Copyright © 2002 John Wiley & Sons, Ltd.     

储层泥浆侵入深度预测方法研究

泥浆侵入半径的确定尚缺乏严格的验证标准.以油水两相渗流理论和离子扩散方程为基础,结合储集层特点,研究了不同储层参数下泥浆滤液对地层的侵入特性.数值模拟侵入时间选取10 d和20 d.数值模拟结果表明,侵入半径在渗透率不变的情况下随孔隙度的增大而减小,在孔隙度不变的情况下随渗透率的增大而增大;当渗透率和孔隙度都发生变化时,泥浆侵入半径一般随孔隙度的增加呈幂函数增加.依据这种关系对测井资料约束处理,得到的泥浆侵入半径较客观地反映了地层的真实情况.     

Prediction of brittle-to-ductile transitions in polystyrene

In this study it is attempted to predict brittle-to-ductile transitions (BDTs) in polystyrene blends, induced either by an increase in temperature or by a decrease in inter-particle distance. A representative, two-dimensional volume element (RVE) of a polystyrene matrix with 20% circular voids, is deformed in tension. During deformation a hydrostatic-stress based craze-nucleation criterion [1] is evaluated. The simulations demonstrate that crazes initiate at low temperatures while a transition from crazing to shear yielding (BDT) is found around 75 °C. The numerical results correlate well with tensile tests on similar heterogeneous polystyrene. The presence of an absolute length, as experimentally found, is more difficult to explain. Near a free surface a T_g-depression is measured for polystyrene and also the resistance to indentation in polystyrene is lower than expected from bulk properties. Both observations are rationalised by an enhanced segmental mobility of chains near a free surface. As a consequence of these findings, an absolute length-scale could be incorporated in the numerical simulations. For simplicity, the length-scale is modelled by taking a temperature gradient over a thin layer near the internal free surfaces of the RVE. Deformation of the RVE with different absolute length-scales shows that indeed also the experimentally found brittle-to-ductile transition can be predicted if the ligament thickness between the inclusions (‘voids’) in polystyrene is below a critical value of ca. 15 nm.     

数值试井在白6-2井测试资料分析评价中的应用

数值试井是在近来发展的一项新的试井解释技术，具有处理非均质，复杂边界油气藏问题的优点，弥补了常规解析试井技术的不足。运用数值试井解释技术，结合动、静态资料，通过对白6-1断块白6-2井测试资料的综合分析评价，确定了储层渗流参数，落实了该气藏的边界情况，为该断块下一步开发提供可靠依据。     

伶仃洋三维潮流输沙的数值模拟

作者使用垂向坐标变换技术，结合有限差分方法建立了三维潮流输沙的数学模型。为了增加模型的稳定性和减少数值耗散，对悬沙对流扩散方程中的水平对流项采用三阶迎风格式离散。应用本模型对伶仃洋的枯季和洪季的潮流输沙的悬沙场进行了数值模拟，将模拟结果与实测资料进行比较，二者吻合较好。模拟结果分析表明：伶仃洋悬沙运动主要集中在西滩，对伶仃航道影响不大；洪季水体含沙浓度高于枯季；水体含沙浓度沿水深存在浓度梯度，且洪季的浓度梯度明显大于枯季。     

Dynamic crack propagation with cohesive elements: a methodology to address mesh dependency

In this paper, two brittle fracture problems are numerically simulated: the failure of a ceramic ring under centrifugal loading and crack branching in a PMMA strip. A three‐dimensional finite element package in which cohesive elements are dynamically inserted has been developed. The cohesive elements' strength is chosen to follow a modified weakest link Weibull distribution. The probability of introducing a weak cohesive element is set to increase with the cohesive element size. This reflects the physically based effect according to which larger elements are more likely to contain defects. The calculations illustrate how the area dependence of the Weibull model can be used to effectively address mesh dependency. On the other hand, regular Weibull distributions have failed to reduce mesh dependency for the examples shown in this paper. The ceramic ring calculations revealed that two distinct phenomena appear depending on the magnitude of the Weibull modulus. For low Weibull modulus, the fragmentation of the ring is dominated by heterogeneities. Whereas many cracks were generated, few of them could propagate to the outer surface. Monte Carlo simulations revealed that for highly heterogeneous rings, the number of small fragments was large and that few large fragments were generated. For high Weibull modulus, signifying that the ring is close to being homogeneous, the fragmentation process was very different. Monte Carlo simulations highlighted that a larger number of large fragments are generated due to crack branching. Copyright © 2003 John Wiley & Sons, Ltd.     

Hydrothermal Reaction Mechanism and Pathway for the Formation of K2Ti6O13 Nanowires

Calculations and detailed first principle and thermodynamic analyses have been performed to understand the formation mechanism of K₂Ti₆O₁₃ nanowires (NWs) by a hydrothermal reaction between bulk Na₂Ti₃O₇ crystals and a KOH solution. It is found that direct ion exchange between K⁺ and Na⁺ plus H⁺ interactions with [TiO₆] octahedra in Na₂Ti₃O₇ promote the formation of an intermediate H₂K₂Ti₆O₁₄ phase. The large lattice mismatch between this intermediate phase and the bulk Na₂Ti₃O₇ structure, and the large energy reduction associated with the formation of this intermediate phase, drive the splitting of the bulk crystal into H₂K₂Ti₆O₁₄ NWs. However, these NWs are not stable because of large [TiO₆] octahedra distortion and are subject to a dehydration process, which results in uniform K₂Ti₆O₁₃ NWs with narrowly distributed diameters of around 10 nm.     

Analysis of S-SEED''s characteristics in optical switch

This article introduces the basic structure of a symmetric self-electrooptic effect device （S-SEED）, and applies the Kirchoff＇ s current law and a purely equivalent capacitive model, to analyze S-SEED＇s switch characteristics. Linear approximation and N-segment approximation are utilized to obtain S-SEED＇s voltage-time （V-T） and characteristics. Theoretical analysis is verified by simulations, and the results demonstrate that the precision of S-SEED＇s switch time can satisfy the requirement in applications with linear approximation. Moreover, the simulations compare S-SEED＇s switch characteristics with different input powers and input contrast ratios, which reveal that increasing input contrast ratio is an effective way to improve S-SEED＇s switch characteristics.