Multilevel preconditioners for a quadrature Galerkin solution of a biharmonic problem

Efficient multilevel preconditioners are developed and analyzed for the quadrature finite element Galerkin approximation of the biharmonic Dirichlet problem. The quadrature scheme is formulated using the Bogner–Fox–Schmit rectangular element and the product two‐point Gaussian quadrature. The proposed additive and multiplicative preconditioners are uniformly spectrally equivalent to the operator of the quadrature scheme. The preconditioners are implemented by optimal algorithms, and they are used to accelerate convergence of the preconditioned conjugate gradient method. Numerical results are presented demonstrating efficiency of the preconditioners. © 2005 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 2006     

逆（反）对策问题与“奇门遁甲”预测理论（术）

本文首先提出逆（反）对策这一新问题，给出了数学模型；探讨了“奇门遁甲”预测理论（术）中的数学问题；通过系统分析“专门遁甲”预测过程，可知它的预测过程隐含着一个特殊的逆（反）对策问题；最后指出逆（反）对策问题的广泛存在并给出案例分析．     

On weak spreadability of distributed-parameter systems and its achievement via linear--quadratic control techniques

In this paper, the concept of weak spreadability is introduced.It constitutes an extension of the idea developed by El Jai& Kassara. In the case of linear distributed systems weconsider quadratic control techniques with a conveniently penalizedcriterion which makes the system weakly spreadable. The approachis outlined for a convection—diffusion system, and theresults of a numerical study are also included to illustratethe main features of the considered problem.     

On a Question of Phillips

In [5] Phillips proved that one can obtain the additive group of any nonstandard model *? of the ring ? of integers by using a linear mod 1 function h : F ?, where F is the α-dimensional vector space over ? when α is the cardinality of *?. In this connection it arises the question whether there are linear mod 1 functions which are neither addition nor quasi-linear. We prove that this is the case.     

Some geometrical aspects of semidefinite linear complementarity problems

This article studies some geometrical aspects of the semidefinite linear complementarity problem (SDLCP), which can be viewed as a generalization of the well-known linear complementarity problem (LCP). SDLCP is a special case of a complementarity problem over a closed convex cone, where the cone considered is the closed convex cone of positive semidefinite matrices. It arises naturally in the unified formulation of a pair of primal-dual semidefinite programming problems. In this article, we introduce the notion of complementary cones in the semidefinite setting using the faces of the cone of positive semidefinite matrices and show that unlike complementary cones induced by an LCP, semidefinite complementary cones need not be closed. However, under R₀-property of the linear transformation, closedness of all the semidefinite complementary cones induced by L is ensured. We also introduce the notion of a principal subtransformation with respect to a face of the cone of positive semidefinite matrices and show that for a self-adjoint linear transformation, strict copositivity is equivalent to strict semimonotonicity of each principal subtransformation. Besides the above, various other solution properties of SDLCP will be interpreted and studied geometrically.     

An experiment with arithmetic precision in linear algebra computations

Several techniques for experimental determination of floating point precision in practical computations are examined, and applied to linear algebra algorithms. These techniques are simple enough to be directly applicable to existing production codes, requiring a very limited amount of software on many machines, and yet they yield interesting information on the numerical precision of a computation. Our choice of linear algebra algorithms includes a direct solver (namely the MA32 program from the Harwell Library) and several variants of preconditioned conjugate gradients (the methods DIAG, INV, MINV and POL of Reference 1). The results may be of interest as method selection criteria, and thus complement Mflop performance data available from several sources.     

Numerical analysis of finite amplitude motion of waves and a moored floating body under severe storm conditions

The motion of a moored floating body under the action of wave forces, which is influenced by fluid forces, shape of the floating body and mooring forces, should be analysed as a complex coupled motion system. Especially under severe storm conditions or resonant motion of the floating body it is necessary to consider finite amplitude motions of the waves, the floating body and the mooring lines as well as non-linear interactions of these finite amplitude motions. The problem of a floating body has been studied on the basis of linear wave theory by many researchers. However, the finite amplitude motion under a correlated motion system has rarely been taken into account. This paper presents a numerical method for calculating the finite amplitude motion when a floating body is moored by non-linear mooring lines such as chains and cables under severe storm conditions.     

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.     

A posteriori H1 error estimation for a singularly perturbed reaction diffusion problem on anisotropic meshes

The paper deals with a singularly perturbed reaction diffusionmodel problem. The focus is on reliable a posteriori error estimatorsfor the H¹ seminorm that can be applied to anisotropic finiteelement meshes. A residual error estimator and a local problemerror estimator are proposed and rigorously analysed. They arelocally equivalent, and both bound the error reliably. Threemodifications of these estimators are introduced and discussed. Much attention is given to the performance of the error estimatorin numerical experiments. This helps to identify those estimatorsthat are suitable for practical applications.     

A new VOF‐based numerical scheme for the simulation of fluid flow with free surface. Part II: application to the cavity filling and sloshing problems

Finite element analysis of fluid flow with moving free surface has been performed in 2‐D and 3‐D. The new VOF‐based numerical algorithm that has been proposed by the present authors (Int. J. Numer. Meth. Fluids, submitted) was applied to several 2‐D and 3‐D free surface flow problems. The proposed free surface tracking scheme is based on two numerical tools; the orientation vector to represent the free surface orientation in each cell and the baby‐cell to determine the fluid volume flux at each cell boundary. The proposed numerical algorithm has been applied to 2‐D and 3‐D cavity filling and sloshing problems in order to demonstrate the versatility and effectiveness of the scheme. The proposed numerical algorithm resolved successfully the free surfaces interacting with each other. The simulated results demonstrated applicability of the proposed numerical algorithm to the practical problems of large free surface motion. It has been also demonstrated that the proposed free surface tracking scheme can be easily implemented in any irregular non‐uniform grid systems and can be extended to 3‐D free surface flow problems without additional efforts. Copyright © 2003 John Wiley & Sons, Ltd.