Method of lines with boundary elements for 1‐D transient diffusion‐reaction problems

Time‐dependent differential equations can be solved using the concept of method of lines (MOL) together with the boundary element (BE) representation for the spatial linear part of the equation. The BE method alleviates the need for spatial discretization and casts the problem in an integral format. Hence errors associated with the numerical approximation of the spatial derivatives are totally eliminated. An element level local cubic approximation is used for the variable at each time step to facilitate the time marching and the nonlinear terms are represented in a semi‐implicit manner by a local linearization at each time step. The accuracy of the method has been illustrated on a number of test problems of engineering significance. © 2005 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 2006     

Finite Element Solution of Conical Diffraction Problems

This paper is devoted to the numerical study of diffraction by periodic structures of plane waves under oblique incidence. For this situation Maxwell's equations can be reduced to a system of two Helmholtz equations in R ² coupled via quasiperiodic transmission conditions on the piecewise smooth interfaces between different materials. The numerical analysis is based on a strongly elliptic variational formulation of the differential problem in a bounded periodic cell involving nonlocal boundary operators. We obtain existence and uniqueness results for discrete solutions and provide the corresponding error analysis.     

Direct numerical simulation of settling of a large solid particle during bioconvection

Settling of a large solid particle in bioconvection flow caused by gyrotactic microorganisms is investigated. The particle is released from the top of the bioconvection chamber; its settling pattern depends on whether it is released in the centre of the bioconvection plume or at its periphery. The Chimera method is utilized; a subgrid is generated around a moving particle. The method suggested by Liu and Wang (Comput. Fluid 2004; 33 :223–255) is further developed to account for the presence of a moving boundary in the streamfunction‐vorticity formulation using the finite‐difference method. A number of cases for different release positions of the particle are computed. It is demonstrated that bioconvection can either accelerate or decelerate settling of the particle depending on the initial position of the particle relative to the plume centre. It is also shown that the particle impacts bioconvection plume by changing its shape and location in the chamber. Copyright © 2005 John Wiley & Sons, Ltd.     

On a new calcium vanadate: synthesis, structure and Li insertion behavior

A synthetic form of the mineral hewettite was prepared via a new route in aqueous medium, starting either from the crystalline compound Li_1.1V₃O₈, or from its amorphous precursor. The anhydrous, crystalline derivative Ca_0.5V₃O₈ was obtained by heating the synthetic hewettite at 250°C under dynamic vacuum. The diffraction studies show that the 2D structure of Ca_0.5V₃O₈ involves the same V₃O₈ layers as in the hewettite or in Li_1+αV₃O₈. The stacking of the layers is similar to that in the metahewettite. A structural model is proposed, where the Ca²⁺ ions occupy octahedral sites in the interlayer space. The electrochemical behavior of Ca_0.5V₃O₈ vs. lithium insertion is presented. It is original and reveals particularly good performances in terms of stability during cycling at C/5 rate. The homologues obtained with Mg or Ba, instead of Ca, are briefly presented.     

A subdomain boundary element method for high‐Reynolds laminar flow using stream function‐vorticity formulation

The paper presents a new formulation of the integral boundary element method (BEM) using subdomain technique. A continuous approximation of the function and the function derivative in the direction normal to the boundary element (further ‘normal flux’) is introduced for solving the general form of a parabolic diffusion‐convective equation. Double nodes for normal flux approximation are used. The gradient continuity is required at the interior subdomain corners where compatibility and equilibrium interface conditions are prescribed. The obtained system matrix with more equations than unknowns is solved using the fast iterative linear least squares based solver. The robustness and stability of the developed formulation is shown on the cases of a backward‐facing step flow and a square‐driven cavity flow up to the Reynolds number value 50 000. Copyright © 2004 John Wiley & Sons, Ltd.     

A complete boundary integral formulation for compressible Navier–Stokes equations

A complete boundary integral formulation for compressible Navier–Stokes equations with time discretization by operator splitting is developed using the fundamental solutions of the Helmholtz operator equation with different order. The numerical results for wall pressure and wall skin friction of two‐dimensional compressible laminar viscous flow around airfoils are in good agreement with field numerical methods. Copyright © 2004 John Wiley & Sons, Ltd.     

Long-term optimal operation of series parallel reservoirs for critical period with specified monthly generation and average monthly storage

We present in this paper an efficient approach for solving the problem of planning the long-term (multiyear) operation of a multireservoir hydroelectric power system for the critical period with a monthly variable load. This load is equal to a certain percentage of the total generation at the end of the year, subject to satisfying a number of constraints on the hydrosystem, using the minimum norm formulation.The proposed method is efficient in computing time and in calculating the total expected benefits from the system during the critical period. Numerical results are reported for a real system in operation consisting of two rivers. Each river has two series reservoirs.This work was supported by the Natural Science and Engineering Research Council of Canada, Grant No. A4146.     

无K1,r图中的哈密顿圈（英文）

本文借助对图的本质独立集和图的部分平方图的独立集的研究，对于K1,r图中哈密顿圈的存在性给出了八个充分条件。我们将利用T－插点技术对这八个充分条件给出统一的证明，本文的结果从本质上改进了C－Q.Zhang于1988年利用次形条件给出的k-连通无爪图是哈密顿图的次型充分条件，同时。G.Chen和R.H.Schelp在1995年利用次型条件给出的关于k-连通无K1，4图是哈密顿图的充分条件也被我们的结果改进并推广到无K1,r图。     

电子储存环中插入元件的非线性束流动力学的Lie代数方法

介绍了用能保证哈密顿系统“辛”性质的Lie映射方法来研究电子储存环中的插入元件对束流动力学的影响.并以此方法对合肥光源中将要安装的波荡器UD?–?1对储存环的影响作了初步的研究     

Gk^d图的平面性和平面图的色数

利用欧拉公式研究了Gdk图的平面性,获得了一个重要定理,并由此得到了关于平面图色数的一个结论.