A nodal coupling of finite and boundary elements

This article presents and analyzes a simple method for the exterior Laplace equation through the coupling of finite and boundary element methods. The main novelty is the use of a smooth parametric artificial boundary where boundary elements fit without effort together with a straight approximate triangulation in the bounded area, with the coupling done only in nodes. A numerically integrated version of the algorithm is also analyzed. Finally, an isoparametric variant with higher order is proposed. © 2003 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 19: 555–570, 2003     

Molecular dynamics simulation of entropy and surface tension for grain boundary of α-Fe

The grain boundary is an interface and the surface tension is one of its important thermodynamic properties. In this paper, the surface tension of the ∑9 grain boundary for α-Fe at various temperatures and pressures is calculated by means of Computer Molecular Dynamics (CMD). The results agree satisfactorily with the experimental data. It is shown that the contribution of entropy to surface tension of grain boundary can be ignored. The project supported by the National Natural Science Foundation of China and the Science Foundation of Chinese Academy of Sciences.     

On the properties of solutions to the Goursat-Darboux problem with boundary and distributed controls

We consider a control system described by the Goursat-Darboux equation. The system is controlled by distributed and boundary controls. The controls are subject to the constraints given as multivalued mappings with closed, possibly nonconvex, values depending on the phase variable. Alongside the initial constraints, we consider the convexified constraints and the constraints whose values are the extreme points of the convexified constraints. We study the questions of existence of solutions and establish connections between the solutions under various constraints.     

Front roughening in three-dimensional imbibition

We investigate the structure and dynamics of the interface between two immiscible liquids in a three-dimensional disordered porous medium. We apply a phase-field model that includes explicitly disorder and discuss both spontaneous and forced imbibition. The structure of the interface is dominated by a length scale ξ_× which arises from liquid conservation. We further show that disorder in the capillary and permeability act on different length scales and give rise to different scalings and structures of the interface properties. We conclude with a range of applications.     

Fixed-sign solutions for a system of singular focal boundary value problems

We consider a system of focal boundary value problems where the nonlinearities may be singular in the independent variable and may also be singular in the dependent arguments. Using Schauder fixed point theorem, we establish criteria such that the system of boundary value problems has at least one fixed-sign solution.     

Generalized 3G theorem and application to relativistic stable process on non-smooth open sets

Let G(x,y) and GD(x,y) be the Green functions of rotationally invariant symmetric α-stable process in Rd and in an open set D, respectively, where 0<α<2. The inequality GD(x,y)GD(y,z)/GD(x,z)?c(G(x,y)+G(y,z)) is a very useful tool in studying (local) Schrödinger operators. When the above inequality is true with c=c(D)∈(0,∞), then we say that the 3G theorem holds in D. In this paper, we establish a generalized version of 3G theorem when D is a bounded κ-fat open set, which includes a bounded John domain. The 3G we consider is of the form GD(x,y)GD(z,w)/GD(x,w), where y may be different from z. When y=z, we recover the usual 3G. The 3G form GD(x,y)GD(z,w)/GD(x,w) appears in non-local Schrödinger operator theory. Using our generalized 3G theorem, we give a concrete class of functions belonging to the non-local Kato class, introduced by Chen and Song, on κ-fat open sets. As an application, we discuss relativistic α-stable processes (relativistic Hamiltonian when α=1) in κ-fat open sets. We identify the Martin boundary and the minimal Martin boundary with the Euclidean boundary for relativistic α-stable processes in κ-fat open sets. Furthermore, we show that relative Fatou type theorem is true for relativistic stable processes in κ-fat open sets. The main results of this paper hold for a large class of symmetric Markov processes, as are illustrated in the last section of this paper. We also discuss the generalized 3G theorem for a large class of symmetric stable Lévy processes.     

Non-self-adjoint sturm-liouville operators with matrix potentials

We obtain asymptotic formulas for non-self-adjoint operators generated by the Sturm-Liouville system and quasiperiodic boundary conditions. Using these asymptotic formulas, we obtain conditions on the potential for which the system of root vectors of the operator under consideration forms a Riesz basis.     

On the convergence rate and optimization of a numerical method with splitting of boundary conditions for the stokes system in a spherical layer in the axisymmetric case: Modification for thick layers

The convergence rate of a fast-converging second-order accurate iterative method with splitting of boundary conditions constructed by the authors for solving an axisymmetric Dirichlet boundary value problem for the Stokes system in a spherical gap is studied numerically. For R/r exceeding about 30, where r and R are the radii of the inner and outer boundary spheres, it is established that the convergence rate of the method is lower (and considerably lower for large R/r) than the convergence rate of its differential version. For this reason, a really simpler, more slowly converging modification of the original method is constructed on the differential level and a finite-element implementation of this modification is built. Numerical experiments have revealed that this modification has the same convergence rate as its differential counterpart for R/r of up to 5 × 10³. When the multigrid method is used to solve the split and auxiliary boundary value problems arising at iterations, the modification is more efficient than the original method starting from R/r ～ 30 and is considerably more efficient for large values of R/r. It is also established that the convergence rates of both methods depend little on the stretching coefficient η of circularly rectangular mesh cells in a range of η that is well sufficient for effective use of the multigrid method for arbitrary values of R/r smaller than ～ 5 × 10³.     

Error estimate for the approximation of nonlinear conservation laws on bounded domains by the finite volume method

In this paper we derive a priori and a posteriori error estimates for cell centered finite volume approximations of nonlinear conservation laws on polygonal bounded domains. Numerical experiments show the applicability of the a posteriori result for the derivation of local adaptive solution strategies.

     

介绍一种二元阈值方法在股票指数上的应用

二元极值的阈值方法的一个发展是用来考虑两个变量的联合分布。这个方法是建立在二元极值的点过程表示法的基础上。本文用参数 (Logistic模型 )和非参数模型对 1992 1999年的上海、深圳日收盘指数对数收益进行分析并给出分析结果。