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%).     

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。提示改良组膀胱功能损伤明显低于对照组。改良的拔管方法是一种值得护理推广应用的好方法。     

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     

线性分式规划最优解集的求法

本文使用多面集的表示定理，导出了线性分式规划最优解集的结构，并给出确定全部最优解的计算步骤。     

Truncated Newton method for sparse unconstrained optimization using automatic differentiation

When solving large complex optimization problems, the user is faced with three major problems. These are (i) the cost in human time in obtaining accurate expressions for the derivatives involved; (ii) the need to store second derivative information; and (iii), of lessening importance, the time taken to solve the problem on the computer. For many problems, a significant part of the latter can be attributed to solving Newton-like equations. In the algorithm described, the equations are solved using a conjugate direction method that only needs the Hessian at the current point when it is multiplied by a trial vector. In this paper, we present a method that finds this product using automatic differentiation while only requiring vector storage. The method takes advantage of any sparsity in the Hessian matrix and computes exact derivatives. It avoids the complexity of symbolic differentiation, the inaccuracy of numerical differentiation, the labor of finding analytic derivatives, and the need for matrix store. When far from a minimum, an accurate solution to the Newton equations is not justified, so an approximate solution is obtained by using a version of Dembo and Steihaug's truncated Newton algorithm (Ref. 1).This paper was presented at the SIAM National Meeting, Boston, Massachusetts, 1986.     

不完全信息群体决策专家权重的集结

本文对于属性权重信息和属性效用信息都不完全的群体多属性决策问题，通过构造属性值区间和运用系统聚类分析法，对群体决策中的专家进行分类，并确定每位专家的权重．     

Efficient computation of natural convection in a concentric annulus between an outer square cylinder and an inner circular cylinder

In this work, the natural convection in a concentric annulus between a cold outer square cylinder and a heated inner circular cylinder is simulated using the differential quadrature (DQ) method. The vorticity‐stream function formulation is used as the governing equation, and the coordinate transformation technique is introduced in the DQ computation. It is shown in this paper that the outer square boundary can be approximated by a super elliptic function. As a result, the coordinate transformation from the physical domain to the computational domain is set up by an analytical expression, and all the geometrical parameters can be computed exactly. Numerical results for Rayleigh numbers range from 104 to 106 and aspect ratios between 1.67 and 5.0 are presented, which are in a good agreement with available data in the literature. It is found that both the aspect ratio and the Rayleigh number are critical to the patterns of flow and thermal fields. The present study suggests that a critical aspect ratio may exist at high Rayleigh number to distinguish the flow and thermal patterns. Copyright © 2002 John Wiley & Sons, Ltd.     

A method for calculating invariant subspaces of symmetric hyperbolic equations

An algorithm is constructed for calculating invariant subspaces of symmetric hyperbolic systems arising in electromagnetic, acoustic, and elasticity problems. Discrete approximations are calculated for subspaces that correspond to minimal eigenvalues and smooth eigenfunctions. Difficulties related to the presence of an infinite-dimensional kernel in the differential operator are successfully handled. The efficiency of the algorithm is demonstrated using acoustics equations.     

有限时滞微分系统依照双测度的实际稳定性

本文对具有有限时滞的泛函微分方程建立了关于依照两种测度的实际稳定性的Razumikhin型判定定理，其中未采用通常的辅助函数，且可运用多个含有状态变量x的部分变元的Lyapunov函数，得出部分变远实际稳定性的判定定理，从而改进了已有的结果。