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.     

求解非光滑方程组的Newton-GMRES方法的全局收敛性分析

1引言对于非光滑方程组F(y)=0(1.1)的求解,这里F:D(?)R~n→R~n是一个非光滑映射,目前主要有两种求解方法.一种方法是由Pang提出的.     

Maximum Principles for Infinite Dimensional Diffusion Equations

In this paper we prove the validity of the Maximum Principle for some class of elliptic and parabolic equations of diffusion type in infinite dimension. The main tools are Asplund’s theorem and Preiss’ theorem on differentiability of Lipschitz functions in Banach space.      

Differential equations uniquely determined by algebras of point symmetries

We continue to investigate strongly and weakly Lie remarkable equations, which we defined in a recent paper. We consider some relevant algebras of vector fields on ℝ^k (such as the isometric, affine, projective, or conformal algebras) and characterize strongly Lie remarkable equations admitted by the considered Lie algebras. __________ Translated from Teoreticheskaya i Matematicheskaya Fizika, Vol. 151, No. 3, pp. 486–494, June, 2007.     

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

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

Identification of a dissipation coefficient by a variational method

A generalized inverse problem for the identification of the absorption coefficient for a hyperbolic system is considered. The well-posedness of the problem is examined. It is proved that the regular part of the solution is an L ₂ function, which reduces the inverse problem to minimizing the error functional. The gradient of the functional is determined in explicit form from the adjoint problem, and approximate formulas for its calculation are derived. A regularization algorithm for the solution of the inverse problem is considered. Numerical results obtained for various excitation sources are displayed.     

Linearisation and potential symmetries of certain systems of diffusion equations

We consider systems of two pure one-dimensional diffusion equations that have considerable interest in Soil Science and Mathematical Biology. We construct non-local symmetries for these systems. These are determined by expressing the equations in a partially and wholly conserved form, and then by performing a potential symmetry analysis on those systems that can be linearised. We give several examples of such systems, and in a specific case we show how linearising and hodograph-type mappings can lead to new solutions of the diffusion system.