分形粗糙面的电磁散射－扩展边界条件

用扩展边界条件方法对分形粗糙良导体面及介质面的电磁散射问题进行了分析。用推广的Floquet模式，在分界面处将场分量用Fourier级数展开，根据边界条件及扩展边界条件得到了水平极化和垂直极化散射场的幅度分量的表达式。用其它近似方法（Kirchhoff和Rayleigh方法）及能量守恒准则验证了此方法的有效性。     

河口水流、波浪、潮流、泥沙、河床变形二维数学模型

根据河口“波浪掀沙，潮流输沙”的泥沙运动机制，引入“波浪幅射应力”反应波浪对水体和泥沙的驱动力．考虑波浪对底部泥沙作用力、波浪对水流挟沙能力的影响及风场对水流和泥沙运动的影响．同时考虑泥沙粒径大小、水体含盐度、含沙量及水流紊动程度等因素对絮凝沉速的影响．用ＡＤＩ算法求解沿水深平均方程，破开算子法解悬移质不平衡输沙方程及河床变形方程，引入通度概念处理岸滩不规则边界，建立河口水流、波浪、潮流、泥沙、河床变形二维数学模型．通过实测资料及卫星摇感图像对所建数学模型进行检验，表明它可以较好地预测河口泥沙的运动     

A semi-analytic boundary element method for parabolic problems

A new semi-analytic solution method is proposed for solving linear parabolic problems using the boundary element method. This method constructs a solution as an eigenfunction expansion using separation of variables. The eigenfunctions are determined using the dual reciprocity boundary element method. This separation of variables-dual reciprocity method (SOV-DRM) allows a solution to be determined without requiring either time-stepping or domain discretisation. The accuracy and computational efficiency of the SOV-DRM is found to improve as time increases. These properties make the SOV-DRM an attractive technique for solving parabolic problems.     

Symmetric BEM formulations for 2-D steady-state and transit potential problems

This paper presents a new concept for symmetric boundary element method (SBEM) applicable to 2-D steady-state and transit potential problems. Two kinds of SBEM formulations are derived. Symmetry is obtained simply through matrix manipulation, and no hypersingularity appears. Therefore, SBEM is much easier than the traditional symmetric Galerkin BEM. Compared with the traditional asymmetric BEM, the present SBEM can reduce the computational cost for time domain problems only. However, when applied to BEM/FEM coupling procedure, SBEM can reduce the computational cost for both steady-state and time domain problems. Three numerical examples are included to illustrate the effectiveness and accuracy of the present formulations.     

Exact integration in the boundary element method for two-dimensional elastostatic problems

This paper derives the exact integrations for the integrals in the boundary element analysis of two-dimensional elastostatics. For facilitation, the derivation is based on the simple forms of the fundamental functions by taking constant, discontinuous linear and discontinuous quadratic elements as examples. The efficiency and accuracy of the derived exact integrations are verified against five benchmark problems; the results indicate that the derived exact integrations significantly reduces the CPU time for forming the matrices of the boundary element analysis and solving the internal displacements.     

An inverse solution methodology for estimating the diffusion coefficient of gases in Athabasca bitumen from pressure-decay data

An inverse solution methodology is developed for the estimation of diffusion coefficient of gases in highly viscous, oil-sands bitumens from isothermal, pressure-decay measurements. The approach involves modeling the rate of change in pressure using the diffusion equation for the liquid phase coupled with a mass balance equation for the gas phase. The inverse solution framework is utilized to arrive at two graphical techniques for estimating the diffusion coefficient. Both techniques involve the determination of the slope of a straight line resulting from plotting the experimental data in accordance with the developed model. An advantage of the proposed techniques is that the diffusion coefficient is estimated directly, i.e. without making it an adjustable parameter. The novelty of the proposed method is in its simplicity as well as its ability to isolate portions of the pressure-decay data affected by experimental fluctuations. The effect of the initial pressure on the predicted diffusion coefficient and pressure-decay profile was also investigated. The diffusion coefficients of CO₂, CH₄, C₂H₆ and N₂ in Athabasca bitumen at 50–90 °C and about 8 MPa were estimated and compared with literature values.     

地下水运动数值模拟过程中边界条件问题探讨

本文对地下水运动数值模拟过程中边界条件的涵义 和处理方法进行了分析和讨论。阐述了边界条件所包含的双重意义。指出随着人类活动影响 强度的日益增大，边界条件的处理要面临一些新的更为复杂的问题。在模型预报之前必须首 先对边界条件做出预报。边界条件的预报既要考虑自然因素的作用，同时也要考虑人类活动 (人工开采和人工补给)的影响及由于邻区水流条件变化而产生的耦合效应。之后，给出了两 个应用实例。     

建筑制图标准中图形符号适宜性、区分性的研究

应用线段长短、间距大小视觉反应的量化规律，提出设计图形符号大小的原则——既更使其大小满足“长”“短”“大”“小”的心理要求，又要互相之间在大小上有区分性。根据图形组织的相似、接近、连续三原则提出设计图形符号形状的原则——图形符号之间尽量减少相似、接近、连续，以保证较好的区分性。通过分析发现，现行建筑制图标准中的图形符号绝大部分在大小和形状具有适宜性、区分性。     

Boundary knot method for 2D and 3D Helmholtz and convection–diffusion problems under complicated geometry

The boundary knot method (BKM) of very recent origin is an inherently meshless, integration‐free, boundary‐type, radial basis function collocation technique for the numerical discretization of general partial differential equation systems. Unlike the method of fundamental solutions, the use of non‐singular general solution in the BKM avoids the unnecessary requirement of constructing a controversial artificial boundary outside the physical domain. The purpose of this paper is to extend the BKM to solve 2D Helmholtz and convection–diffusion problems under rather complicated irregular geometry. The method is also first applied to 3D problems. Numerical experiments validate that the BKM can produce highly accurate solutions using a relatively small number of knots. For inhomogeneous cases, some inner knots are found necessary to guarantee accuracy and stability. The stability and convergence of the BKM are numerically illustrated and the completeness issue is also discussed. Copyright © 2003 John Wiley & Sons, Ltd.     

Fundamental solutions for half plane with an oblique edge crack

An elastic half plane with an oblique edge crack is considered in this paper. A pair of concentrated forces or point dislocations is assumed to act at an arbitrary point in the half plane. The half plane with an edge crack is first mapped into a unit circle by a rational mapping function so that the following analysis can be carried out on the mapped plane analytically. Then the complex stress functions are derived by separating the whole problem into two parts; one is the principal part corresponding to the infinite plane acted on by concentrated forces or dislocations, the other is the holomorphic part, which can be determined by making use of the property of regularity of complex stress functions. The stress intensity factors of the crack can be calculated with different inclined angles of the crack, and the displacement and stress components at an arbitrary position in the half plane can be expressed explicitly.