A modified predictor–corrector scheme for the two-dimensional sine-Gordon equation

A three-time level finite-difference scheme based on a fourth order in time and second order in space approximation has been proposed for the numerical solution of the nonlinear two-dimensional sine-Gordon equation. The method, which is analysed for local truncation error and stability, leads to the solution of a nonlinear system. To avoid solving it, a predictor–corrector scheme using as predictor a second-order explicit scheme is proposed. The procedure of the corrector has been modified by considering as known the already evaluated corrected values instead of the predictor ones. This modified scheme has been tested on the line and circular ring soliton and the numerical experiments have proved that there is an improvement in the accuracy over the standard predictor–corrector implementation. This research was co-funded by E.U. (75%) and by the Greek Government (25%).     

Clifford分析中无界域上正则函数的边值问题

在引入修正Cauchy核的基础上,讨论了无界域上正则函数的带共轭值的边值问题：a(t)Φ (t) b(t)Φ (t) c(t)Φ-(t) d(t)Φ=g(t)．首先给出了无界域上正则函数的Plemelj公式,然后利用积分方程方法和压缩不动点原理证明了问题解的存在唯一性．     

Instabilities in the Chapman-Enskog Expansion and Hyperbolic Burnett Equations

It is well-known that the classical Chapman-Enskog procedure does not work at the level of Burnett equations (the next step after the Navier-Stokes equations). Roughly speaking, the reason is that the solutions of higher equations of hydrodynamics (Burnett's, etc.) become unstable with respect to short-wave perturbations. This problem was recently attacked by several authors who proposed different ways to deal with it. We present in this paper one of possible alternatives. First we deduce a criterion for hyperbolicity of Burnett equations for the general molecular model and show that this criterion is not fulfilled in most typical cases. Then we discuss in more detail the problem of truncation of the Chapman-Enskog expansion and show that the way of truncation is not unique. The general idea of changes of coordinates (based on analogy with the theory of dynamical systems) leads finally to nonlinear Hyperbolic Burnett Equations (HBEs) without using any information beyond the classical Burnett equations. It is proved that HBEs satisfy the linearized H-theorem. The linear version of the problem is studied in more detail, the complete Chapman-Enskog expansion is given for the linear case. A simplified proof of the Slemrod identity for Burnett coefficients is also given.     

两种拔尿管方法对膀胱功能损伤的临床观察

为观察两种不同的拔尿管方法对膀胱功能损伤的差异,将80例患者随机分成两组,分别采用改良法(40例)和常规法(40例)拔除留置尿管,分析比较了两组患者排尿顺利率、逼尿肌功能恢复率及尿道刺激症、尿失禁、尿潴留。结果表明,两组患者拔管后上述指标与情况有显著性差异,P<0.05。提示改良组膀胱功能损伤明显低于对照组。改良的拔管方法是一种值得护理推广应用的好方法。     

TWO-POINT BOUNDARY VALUE PROBLEM FOR A CLASS OF SECOND-ORDER DIFFERENCE EQUATIONS

In this paper, we apply a critical point theorem and obtain the condition for the existence of three solutions to two-point boundary value problem of a second order nonlinear difference equation.     

Irregularities of distribution with respect to strips

Summary A lower bound is established for the strip discrepancy of a broad class of point distributions. The bound implies unbounded strip discrepancy for equally weighted point distributions under favorable conditions. The methods of proof use notions from integral geometry.     

Mathematical modeling of flows through inelastic porous solids

We propose a framework, based on classical mixture theory, to describe the isothermal flow of an incompressible fluid through a deformable inelastic porous solid. The modeling of the behavior of the inelastic solid takes into account changes in the elastic response due to evolution in the microstructure of the material. We apply the model to a compression layer problem. The mathematical problem generated by the model is a free boundary problem.     

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     

Smart Geodesic Tracing in GRworkbench

We have developed a new tool for numerical work in General Relativity: GRworkbench. We discuss how GRworkbench's implementation of a numerically-amenable analogue to Differential Geometry facilitates the development of robust and chart-independent numerical algorithms. We consider, as an example, geodesic tracing on two charts covering the exterior Schwarzschild space-time.     

Adaptive Algorithm for Constrained Least-Squares Problems

This paper is concerned with the implementation and testing of an algorithm for solving constrained least-squares problems. The algorithm is an adaptation to the least-squares case of sequential quadratic programming (SQP) trust-region methods for solving general constrained optimization problems. At each iteration, our local quadratic subproblem includes the use of the Gauss–Newton approximation but also encompasses a structured secant approximation along with tests of when to use this approximation. This method has been tested on a selection of standard problems. The results indicate that, for least-squares problems, the approach taken here is a viable alternative to standard general optimization methods such as the Byrd–Omojokun trust-region method and the Powell damped BFGS line search method.