求解P0函数互补问题的一种微分方程方法

本针对P0函数互补问题，给出了一种微方程方法，并且证明了P0函数互补问题的解是微分方程系统的渐进稳定平衡点。在适当的假设条件下，证明了所给出的算法具有二次收敛速度。几个数值例子表明了该算法的有效性。     

Numerical analysis of time-dependent Boussinesq models

In this paper we analyse numerical models for time-dependent Boussinesq equations. These equations arise when so-called Boussinesq terms are introduced into the shallow water equations. We use the Boussinesq terms proposed by Katapodes and Dingemans. These terms generalize the constant depth terms given by Broer. The shallow water equations are discretized by using fourth-order finite difference formulae for the space derivatives and a fourth-order explicit time integrator. The effect on the stability and accuracy of various discrete Boussinesq terms is investigated. Numerical experiments are presented in the case of a fourth-order Runge-Kutta time integrator.     

Electromagnetic properties of strange baryons in a relativistic quark model

We present some of our results for the electromagnetic properties of excited Σ hyperons, computed within the framework of the Bonn constituent-quark model, which is based on the Bethe-Salpeter approach. The seven parameters entering the model are fitted against the best-known baryon masses. Accordingly, the results for the form factors and helicity amplitudes are genuine predictions. We compare with the scarce experimental data available and discuss the processes in which Σ ^*'s may play an important role.     

Superconvergence of discontinuous Galerkin finite element method for the stationary Navier‐Stokes equations

This article focuses on discontinuous Galerkin method for the two‐ or three‐dimensional stationary incompressible Navier‐Stokes equations. The velocity field is approximated by discontinuous locally solenoidal finite element, and the pressure is approximated by the standard conforming finite element. Then, superconvergence of nonconforming finite element approximations is applied by using least‐squares surface fitting for the stationary Navier‐Stokes equations. The method ameliorates the two noticeable disadvantages about the given finite element pair. Finally, the superconvergence result is provided under some regular assumptions. © 2006 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 23: 421–436, 2007     

On the Use of the Restitution Condition in Flexible Body Dynamics

The difference between the classical treatment offlexible body impact and the treatment of impact in flexiblemultibody dynamics is due to several fundamental reasons. Inthe classical impact theory, simple structures such as beamsand plates are used. Infinite dimensional models can bedeveloped for these simple structural elements to study theimpact dynamics and the wave propagation problem. Flexiblemultibody impact problems, on the other hand, involve bodieswith complex geometry that cannot be modeled using infinitenumber of degrees of freedom. Furthermore, the classicalimpact theory has been mainly concerned with the impactbetween a rigid mass that moves without constraints beforeit impacts a simple flexible structure. This is not amultibody simulation scenario in which the impact occursbetween kinematically constrained bodies that are subjectedto impulsive constraint forces in addition to the impactforces. These constraint forces can influence the motion ofthe two bodies immediately after impact, and as aconsequence, the simple classical theory scenario of impactdoes not apply. It is the objective of this paper to discussthe use of the restitution condition in flexible multibodyimpact problems and demonstrate that the use of thisapproach does not exclude the classical formulation.Nonetheless, the impulse momentum balance approach can serveas an effective and efficient procedure for solving theimpact problem in finite dimensional models that do not obeythe classical wave theory. Energy results of simplestructural elements are presented in order to demonstratethe consistency of using the impulse momentum balanceapproach in solving impact problems in finite dimensionalflexible body applications.     

四阶方程两点边值问题Hermite有限元解的渐近展式与外推

1引言有限元解的渐近展式是提高微分方程数值解精度的重要工具,比如亏量校正和外推就是建立在有限元解的渐近展式的基础之上．许多作者对此进行了大量的研究(见[1]-[4]),特别是文[1],提出了在研究有限元解的渐近展式中十分有用的能量嵌入技巧．本文利用能量嵌入定理得到了四阶方程两点边值问题Hermite有限元解及其二阶平均导数的渐近展式,进一步我们还讨论了它们的Richardson外推公式．考虑四阶方程两点边值问题     

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.     

General synchronization dynamics of coupled Van der Pol–Duffing oscillators

This paper considers the general synchronization dynamics of coupled Van der Pol–Duffing oscillators. The linear and nonlinear stability analysis on the synchronization process is derived through the Whittaker method and the Floquet theory in addition to the multiple time scales method. A stability map displaying different dynamical states of the system is performed. Numerical simulation is carried out to support and to complement the accuracy of the analytical treatment.     

On regularity of a weak solution to the Navier–Stokes equation with generalized impermeability boundary conditions

We consider the 3D Navier–Stokes equation with generalized impermeability boundary conditions. As auxiliary results, we prove the local in time existence of a strong solution (‘strong’ in a limited sense) and a theorem on structure. Then, taking advantage of the boundary conditions, we formulate sufficient conditions for regularity up to the boundary of a weak solution by means of requirements on one of the eigenvalues of the rate of deformation tensor. Finally, we apply these general results to the case of an axially symmetric flow with zero angular velocity.