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     

Solitary waves and their bifurcations of KdV like equation with higher order nonlinearity

We investigate the KdV like equation with higher order nonlinearity ut + a(1 +bun)unux + uxxx = 0with n ≥ 1, a, b ∈ R and α≠ 0. The bifurcations and explicit expressions of solitary wave solutions for theequation are discussed by using the bifurcation method and qualitative theory of dynamical systems. Thebifurcation diagrams, existence and number of the solitary waves are given.     

Semi‐analytical integration of the 8‐node plane element stiffness matrix using symbolic computation

The semi‐analytical integration of an 8‐node plane strain finite element stiffness matrix is presented in this work. The element is assumed to be super‐parametric, having straight sides. Before carrying out the integration, the integral expressions are classified into several groups, thus avoiding duplication of calculations. Symbolic manipulation and integration is used to obtain the basic formulae to evaluate the stiffness matrix. Then, the resulting expressions are postprocessed, optimized, and simplified in order to reduce the computation time. Maple symbolic‐manipulation software was used to generate the closed expressions and to develop the corresponding Fortran code. Comparisons between semi‐analytical integration and numerical integration were made. It was demonstrated that semi‐analytical integration required less CPU time than conventional numerical integration (using Gaussian‐Legendre quadrature) to obtain the stiffness matrix. © 2005 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2006     

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.     

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.     

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

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

Coupled Thermo-Hydro-Mechanical and Chemical Analysis of Expansive Clay Subjected to Heating and Hydration

A fully coupled formulation combining reactive transport and an existing thermo-hydro-mechanical (THM) code is presented. Special attention has been given to phenomena likely to be encountered in clay barriers used as part of containment systems of nuclear waste. The types of processes considered include hydrolysis, complex formation, oxidation/reduction reactions, acid/base reactions, precipitation/dissolution of minerals and cation exchange. Both kinetically-controlled and equilibrium-controlled reactions have been incorporated. The total analytical concentrations (including precipitated minerals) are adopted as basic transport variables and chemical equilibrium is achieved by minimizing Gibbs Free Energy. The formulation has been incorporated in a general purpose computer code capable of performing numerical analysis of engineering problems. A validation exercise concerning a laboratory experiment involving the heating and hydration of an expansive compacted clay is described.     

A fourth‐order compact finite difference scheme for the steady stream function–vorticity formulation of the Navier–Stokes/Boussinesq equations

A fourth‐order compact finite difference scheme on the nine‐point 2D stencil is formulated for solving the steady‐state Navier–Stokes/Boussinesq equations for two‐dimensional, incompressible fluid flow and heat transfer using the stream function–vorticity formulation. The main feature of the new fourth‐order compact scheme is that it allows point‐successive overrelaxation (SOR) or point‐successive underrelaxation iteration for all Rayleigh numbers Ra of physical interest and all Prandtl numbers Pr attempted. Numerical solutions are obtained for the model problem of natural convection in a square cavity with benchmark solutions and compared with some of the accurate results available in the literature. Copyright © 2003 John Wiley & Sons, Ltd.     

Formulation and survey of ALE method in nonlinear solid mechanics

This paper investigates the applicability and accuracy of existing formulation methods in general purpose finite element programs to the finite strain deformation problems. The basic shortcomings in using such programs in these applications are then pointed out and the need for a different type of formulation is discussed. An arbitrary Lagrangian-Eulerian (ALE) method is proposed and a concise survey of ALE formulation is given. A consistent and complete ALE formulation is derived from the virtual work equation transformed to arbitrary computational reference configurations. Differences between the proposed formulations and similar ones in the literature are discussed. The proposed formulation presents a general approach to ALE method. It includes load correction terms and is suitable for rate-dependent and rate-independent material constitutive law. The proposed formulation reduces to both updated Lagrangian and Eulerian formulations as special cases.