The Stability of Runge－Kutta Methods for Systems of Delay Differential Equations

ＴｈｅＳｔａｂｉｌｉｔｙｏｆＲｕｎｇｅ－ＫｕｔｔａＭｅｔｈｏｄｓｆｏｒＳｙｓｔｅｍｓｏｆＤｅｌａｙＤｉｆｆｅｒｅｎｔｉａｌＥｑｕａｔｉｏｎｓ￥ＷＡＮＧＸｉａｏｂｉａｏ；ＬＩＵＭｉｎｇｚｈｕ(王晓彪)(刘明珠)（Ｄｅｐｔ．ｏｆＭａｔｈｅｍａｔｉｃｓ，Ｈａｒ...     

Stress analysis with arbitrary body force by boundary element method

Linear stress analysis without body force can be easily carried out by means of the boundary element method. Some cases of linear stress analysis with body force can also be solved without the domain integral. However domain integrals are generally necessary to solve the linear stress problems with complicated body forces. This paper shows that the linear stress problems with complicated body forces can be solved approximately without the domain integral. In order to solve these problems, the domain is divided into small areas using contour lines of body force. In these areas, the distributions of body force are assumed approximately to satisfy the Laplace equation.     

NUMERICAL SIMULATION OF PARTICLE SEPARATION IN AN OIL-SAND SEPARATOR

ＮＯＴＡＴＩＯＮＳＡ　Ａｒｅａｏｆｉｎｌｅｔ ,ｍ2ａ　Ｃｅｎｔｒｉｆｕｇａｌａｃｃｅｌｅｒａｔｉｏｎ ,ｍ2 /ｓｉ　Ｇｒｉｄｏｆａｘｉａｌｄｉｒｅｃｔｉｏｎｍｌ 　ＭａｘｉｍｕｍｇｒｉｄｎｕｍｂｅｒｏｆｒａｄｉａｌｄｉｒｅｃｔｉｏｎＱ　Ｆｌｕｉｄｍａｓｓｒａｔｅ ,ｍ3/ｓＲｃ 　Ｍａｘｉｍｕｍｒａｄｉｕｓｏｆｃｙｃｌｏｎｅ ,ｍＲｅ　Ｒｅｙｎｏｌｄｓｎｕｍｂｅｒｒ　Ｄｉｍｅｎｓｉｏｎｌｅｓｓｒａｄｉａｌｄｉｓｔａｎｃｅｔ　Ｄｉｍｅｎｓｉｏｎｌｅｓｓｔｉｍｅｕ　Ｄｉｍｅｎｓｉｏｎｌｅｓｓａｘｉａｌｖｅｌｏｃ…     

井眼剪切带分叉研究

井眼基本破坏形式是由于不稳定性造成的井眼岩石崩落坍塌和剪切带。井眼破坏的非均匀性采用传统的理论和方法不能得到较好的解释和预测。采用材料的分叉分析方法和变形局部化的数值模拟研究 ,可望对井眼失效的机理分析、预测和控制提供新的思路     

挤压六通管及其连接结构的安全性研究

六通管是石油化工设备及压力管道中的重要元件,传统的制作及连接方法是在集流管上开孔,然后将支管插入并焊接而成或采用锻制的方法制成。采用挤压的方法制作六通管,将其与支管对接焊接。并通过试验及数值模拟对该结构的安全性进行了研究,实践证明,这种挤压六通管的强度性能优于普通的焊接六通管,可以广泛应用于工业生产。     

孤东油田的海油陆采进展

孤东油田极浅海区域新发现几个边际油藏。其中垦东401块水深-2m,距陆地2.7km,为稀油、疏松砂岩油藏,试油获14～25t/d工业油流。在三维地震、测井约束地震反演、数值模拟和油藏经营管理评价基础上编制了产能建设方案,通过修筑海公路、人工岛、丛式井组、大位移水平井、采用数据采集与控制的远程测控和紧急关断技术,能够无人职守、遥测遥控,整体实现了海油陆采。有效动用地质储量233×104t,钻井12口,建成年产油4×104t能力,累产油8.6×104t。目前在垦东12人工岛上钻打了26口井,推广了海油陆采经验。     

A stability classification method and its application to pipelined solution of linear recurrences

A new method of classification for numerical stability of parallel algorithms is proposed based on the theoretical foundation of forward error analysis. It partitions the algorithms according to their asymptotic stability—a measure introduced to relate the limiting behavior of the stability to the size of the problem. Using this method, the stability aspect of the pipelined solution technique for first-order and second-order linear recurrences—the core of a tridiagonal linear equation solver—is studied. In particular, it shows that the pipelined solution method of the first-order linear recurrences has the same degree of stability as the commonly used sequential evaluation algorithms. The stability problems of sequential and pipelined solution methods of the second-order linear recurrences are also studied.     

关于无限棒束间流动和传热特性的数值研究

应用非线性ｋε湍流模式，采用非正交曲线坐标系下求解三维ＮＳ方程的非交错网格有限体积方法，数值模拟了充分发展条件下，三角形排列无限棒束间通道内，不同的几何参数（Ｐ／Ｄ），雷诺数（Ｒｅ）下的流动和传热问题。给出了不同参数下的速度和温度分布以及湍流二次流动，分析了几何参数、雷诺数及二次流对棒束内流动和传热特性的影响，得到了不同参数下通道的摩擦系数和Ｎｕｓｅｌｔ数，并与经验关系式作了比较     

Isometric Piecewise Polynomial Curves

The main preoccupations of research in computer-aided geometric design have been on shape-specification techniques for polynomial curves and surfaces, and on the continuity between segments or patches. When modelling with such techniques, curves and surfaces can be compressed or expanded arbitrarily. There has been relatively little work on interacting with direct spatial properties of curves and surfaces, such as their arc length or surface area. As a first step, we derive families of parametric piecewise polynomial curves that satisfy various positional and tangential constraints together with arc-length constraints. We call these curves isometric curves. A space curve is defined as a sequence of polynomial curve segments, each of which is defined by the familiar Hermite or Bézier constraints for cubic polynomials; as well, each segment is constrained to have a specified arc length. We demonstrate that this class of curves is attractive and stable. We also describe the numerical techniques used that are sufficient for achieving real time interaction with these curves on low-end workstations.     

Vertical bridgman techniques to homogenize zinc composition of CdZnTe substrates

In order to improve the Zn homogeneity along the axial direction of CdZnTe boule, we have employed a modified Bridgman technique using a (Cd, Zn) alloy source in communication with the melt, whose temperature has been gradually changed from 800 to 840°C during growth. Electron probe microanalysis (EPMA) measurements of Zn composition in the boule shows an excellent homogeneity of Zn along the axis of the CdZnTe boule compared with results in a boule grown by using a fixed source temperature. We have performed a numerical simulation to obtain the approximate temperatures of additional heating and cooling needed to improve the radial Zn homogeneity. CdZnTe boule has been grown by seeded vertical Bridgman furnace with two zones of heater and cooler. Ultraviolet/visible spectroscopic measurements of Zn composition over the length of the boule indicate that the radial distribution of Zn composition is very homogeneous in the body region of the boule, where the radial variation of Zn composition is ±0.0005.