A Constitutive Equation of Thermoelasticity for Solid Materials

A new constitutive equation of thermoelasticity for crystals is presented based on the interatomic potential and solid mechanics at finite temperature. Using the new constitutive equation, the calculations for crystal copper and graphene are carried out under different loading paths at different temperatures. The calculated results are in good agreement with those of the previous thermoelasticity constitutive equation based on quantum mechanics, which clearly indicates that our new constitutive equation of thermoelasticity is correct. A lot of comparisons also show that the present theory is more concise and efficient than the previous thermal stress theory in the practical application.     

The Marginal Density of a TMA(1) Process

This note reconsiders the marginal density of a threshold moving average process and proposes a simple yet effective numerical algorithm to implement that by solving an associated integral equation. This algorithm can also be applied to calculate stationary probability density or distribution functions of a few other types of nonlinear stationary stochastic processes numerically.     

平面阵列波导光栅波导模式参量的分析

用周期波导的方法分析平面阵列波导光栅(AWG)的本征方程及波导模式参数,结果表明波导传播常数是带状分布的,与波导间隔有关。分析结果对器件的优化设计有一定的指导意义。     

大牛地气田异常修正等时试井资料求取产能方法

大牛地气田为低渗-致密气田,因存在启动压力梯度等因素,修正等时试井曲线按常规方法处理出现异常,无法准确获取气井的产能方程.采用LIT法和异常曲线校正法对气田异常试井资料进行了重新处理,最终获得气井产能方程,为同类气藏产能试井可能出现的异常资料处理解释探索了新的途径.     

An overlapping domain decomposition dual reciprocity method

We present a domain decomposition boundary integral equation method for two-dimensional partial differential equations. The overlapping Schwarz method is employed to improve the dual reciprocity boundary element method. The resulting algorithm turns out to be more accurate than a non-overlapping approach previously proposed. Some numerical results showing the improved accuracy and efficiency of the methods are given.     

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.     

大牛地气田气井无阻流量的确定方法研究

大牛地气田气井试井大多采用一点法，由于气层渗透性差，测试时间短，流压很难达到稳定，计算的无阻流量往往偏高，且无法获得正确的产能方程。通过气田多口井修正等时试井数据进行统计分析，获得了一点法试井在不同生产时间下计算无阻流量的经验公式；利用气田修正等时试井的二项式方程系数曰，计算惯性阻力系数并经过曲线回归处理获得了经验关系式，很好地解决了一点法试井计算产能方程的问题。     

波动方程数值模拟的隐式差分法

在有限差分波动方程数值模拟中,通常采用高阶差分方法来提高空间导数的数值逼近精度,以实现降低数值频散,提高数值模拟精度的目的。首先对差分频散进行了理论分析;然后讨论了估计一阶空间导数的隐式差分格式,并与通常采用的高阶精度显式差分格式进行了对比分析,结果表明,隐式差分格式能够在更宽的波数范围使差分频散控制在可接受的水平,如8阶精度的显式差分格式所适应的波数带宽约为O．55kmax,而隐式差分格式所适应的波数带宽约为0．7kmax;最后通过模型试算,对隐式差分格式的有效性进行了验证。模拟结果表明,用隐式差分格式在一定程度上降低了差分频散,提高了模拟精度。     

A multigrid accelerated hybrid unstructured mesh method for 3D compressible turbulent flow

 A cell vertex finite volume method for the solution of steady compressible turbulent flow problems on unstructured hybrid meshes of tetrahedra, prisms, pyramids and hexahedra is described. These hybrid meshes are constructed by firstly discretising the computational domain using tetrahedral elements and then by merging certain tetrahedra. A one equation turbulence model is employed and the solution of the steady flow equations is obtained by explicit relaxation. The solution process is accelerated by the addition of a multigrid method, in which the coarse meshes are generated by agglomeration, and by parallelisation. The approach is shown to be effective for the simulation of a number of 3D flows of current practical interest. Sponsored by The Research Council of Norway, project number 125676/410 Dedicated to the memory of Prof. Mike Crisfield, a respected colleague     

透水坡面降雨径流过程的运动波近似解析解

本文从坡面流运动波理论的基本方程出发，利用运动波特征线法和分级叠加法，推导出了实验条件下考虑净雨强随时间变化的坡面流运动波方程近似解析解。用实验实测资料对模拟计算结果进行了验证。结果表明，坡面漫流过程的模拟计算结果与实测资料符合良好。