Accelerating Galerkin BEM for linear elasticity using adaptive cross approximation

The adaptive cross approximation (ACA) algorithm (Numer. Math. 2000; 86 :565–589; Computing 2003; 70 (1):1–24) provides a means to compute data‐sparse approximants of discrete integral formulations of elliptic boundary value problems with almost linear complexity. ACA uses only few of the original entries for the approximation of the whole matrix and is therefore well‐suited to speed up existing computer codes. In this article we extend the convergence proof of ACA to Galerkin discretizations. Additionally, we prove that ACA can be applied to integral formulations of systems of second‐order elliptic operators without adaptation to the respective problem. The results of applying ACA to boundary integral formulations of linear elasticity are reported. Furthermore, we comment on recent implementation issues of ACA for non‐smooth boundaries. Copyright © 2006 John Wiley & Sons, Ltd.     

关于求解常微方程组的离散波形松驰方法的加速收敛

提出一种关于求解常微线性系统的离散波形松驰方法的新的加速收敛技巧。通过对系统矩阵A的分裂，该技巧使迭代矩阵（(zI M)^-1N）具有理想的较小谱半径。在LU分解的基础上给出了一个迭代算法以及用该法与Gauss-Seidel方不垢敛速进行比较的数值例子。     

具常数特征拟线性双曲型方程组Cauchy问题

考虑具常数特征拟线性双曲型方程,提出一个新的可化约方程组的方法,证明了具常特征方程组Cauchy问题经典解的整体存在性定理.同时构造一些例子说明一些有趣的现象.     

Flow in a waterfall with large Froude number

An approximate method is developed to solve the full nonlinear equations governing two-dimensional irrotational flow in a free waterfall, falling under the influence of gravity, at high Froude number based on conditions far upstream. Schwarz—Christoffel transformation is used to map the region, in the complex potential-plane, onto the upper half-plane. The Hilbert transformation as well as the perturbation technique, for large Froude number, are used as a basis for the approximate solution of the problem. A complete solution, up to second-order approximation, for the downstream free-surfaces profiles, for different Froude number, is discussed and illustrated. The obtained approximate solutions are compared with those of other authors. Favourable agreement with other results suggests that this method is effective in dealing with flow problems strongly influenced by gravity and high Froude number. The results obtained by this method are sufficiently accurate for practical purposes.     

OPTIMAL INSTAR PARASITIZATION IN A STAGE STRUCTURED HOST-PARASITOID MODEL

A stage structured host-parasitoid model is derived and the equilibria studied. It is shown under what conditions the parasitoid controls an exponentially growing host in the sense that a coexistence equilibrium exists. Furthermore, for host populations whose inherent growth rate is not too large it is proved that in order to minimize the adult host equilibrium level it is necessary that the parasitoids attack only one of the larval stages. It is also proved in this case that the minimum adult host equilibrium level is attained when the parasitoids attack that larval stage which also maximizes the expected number of emerging adult parasitoid per larva at equilibrium. Numerical simulations tentatively indicate that the first conclusion remains in general valid for the model. However, numerical studies also show that it is not true in general that the optimal strategy will maximize the number of emerging adult parasitoid per larva at equilibrium.     

泊松类型方程边界元解法

本文采用高阶拉普拉斯算子基本解将泊松类型方程的区域积分全部变换成边界积分，使计算问题的维数减少一维．通过斯托克斯方程的算例，表明本文所用的方法是有效的方法。     

Dynamics of liquid membranes. II: Adaptive finite difference methods

Two domain-adaptive finite difference methods are presented and applied to study the dynamic response of incompressible, inviscid, axisymmetric liquid membranes subject to imposed sinusoidal pressure oscillations. Both finite difference methods map the time-dependent physical domain whose downstream boundary is unknown onto a fixed computational domain. The location of the unknown time-dependent downstream boundary of the physical domain is determined from the continuity equation and results in an integrodifferential equation which is non-linearly coupled with the partial differential equations which govern the conservation of mass and linear momentum and the radius of the liquid membrane. One of the finite difference methods solves the non-conservative form of the governing equations by means of a block implicit iterative method. This method possesses the property that the Jacobian matrix of the convection fluxes has an eigenvalue of algebraic multiplicity equal to four and of geometric multiplicity equal to one. The second finite difference procedure also uses a block implicit iterative method, but the governing equations are written in conservation law form and contain an axial velocity which is the difference between the physical axial velocity and the grid speed. It is shown that these methods yield almost identical results and are more accurate than the non-adaptive techniques presented in Part I. It is also shown that the actual value of the pressure coefficient determined from linear analyses can be exceeded without affecting the stability and convergence of liquid membranes if the liquid membranes are subjected to sinusoidal pressure variations of sufficiently high frequencies.     

两个非线性方程的准确解

借助于计算机代数系统Ｍａｔｈｅｍａｔｉｃａ利用直接代数方法获得了两个有数学物理意义的非线性耗散－色散方程的准确解，这种方法也适用于高维非线性方程。     

双包层椭圆光波导解析解

解析求解了双包层椭圆光纤中的波动方程，得到了模式精确解有模式特征方程。对基模的特征方程进行了数值计算，给出了不同椭圆比下的归一化双折射和模间色散随归一化频率的变化关系曲线，并与高斯近似解的结果进行了比较。     

ON DIRECT METHOD OF SOLUTION FOR A CLASS OF SINGULAR INTEGRAL EQUATIONS

In this article, by introducing characteristic singular integral operator and associate singular integral equations (SIEs), the authors discuss the direct method of solution for a class of singular integral equations with certain analytic inputs. They obtain both the conditions of solvability and the solutions in closed form. It is noteworthy that the method is different from the classical one that is due to Lu.