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.     

Nondensely defined evolution impulsive differential inclusions with nonlocal conditions

In this paper, we shall establish sufficient conditions for the existence of integral solutions for some nondensely defined evolution impulsive differential inclusions in Banach spaces with nonlocal conditions.     

Runge–Kutta Methods for Itô Stochastic Differential Equations with Scalar Noise

A general class of stochastic Runge–Kutta methods for Itô stochastic differential equation systems w.r.t. a one-dimensional Wiener process is introduced. The colored rooted tree analysis is applied to derive conditions for the coefficients of the stochastic Runge–Kutta method assuring convergence in the weak sense with a prescribed order. Some coefficients for new stochastic Runge–Kutta schemes of order two are calculated explicitly and a simulation study reveals their good performance.     

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

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

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.     

The driving torque of a rotating levitated cylindrical permanent magnet above a superconductor

Due to the Meissner effect, a permanent magnet is levitated, when released above a high temperature superconductor. When there is an inhomogeneous temperature field around the magnet, the magnet might start to oscillate with increasing amplitude until it remains in a continuous rotation. A mathematical model for the described effect is presented which couples heat transfer and electromagnetic forces with the equation of motion, yielding to a multiphysics task. In a detailed analysis it is found, that the torque which drives the rotation of the magnet, is explicitly given in terms of Bessel functions and the Fourier coefficients of order zero and one of the temperature field of the surrounding air.     

A Relaxation Scheme for Solving the Boltzmann Equation Based on the Chapman-Enskog Expansion

Abstract In [16] a visco-elastic relaxation system, called the relaxed Burnett system, was proposed by Jinand Slemrod as a moment approximation to the Boltzmann equation. The relaxed Burnett system is weaklyparabolic, has a linearly hyperbolic convection part, and is endowed with a generalized eotropy inequality. Itagrees with the solution of the Boltzmann equation up to the Burnett order via the Chapman-Enskog expansion. We develop a one-dimensional non-oscillatory numerical scheme based on the relaxed Burnett system forthe Boltzmann equation. We compare numerical results for stationary shocks based on this relaxation scheme,and those obtained by the DSMC (Direct Simulation Monte Carlo), by the Navier-Stokes equations and bythe extended thermodynamics with thirteen moments (the Grad equations). Our numerical experiments showthat the relaxed Burnett gives more accurate approximations to the shock profiles of the Boltzmann equationobtained by the DSMC, for a range of Mach numbers for hypersonic flows, th     

Linear Relations, Differential Inclusions, and Degenerate Semigroups

A generalized linear differential equation in a Banach space is studied. The construction of a phase space and solutions with the help of the spectral theory of linear operators, ergodic theorems, and degenerate semigroups of linear operators is carried out.     

掺铒光纤放大器中上转换过程的影响

掺杂浓度的提高影响了掺铒光纤放大器(EDFA)的性能,对经典Giles模型提出了新的要求。首先研究了高掺杂下铒离子距离较小引起的上转换过程的机理,通过优化速率方程模型分析了上转换过程对EDFA的影响,得出了在泵浦光和信号光远远大于ASE光的假设下泵浦光和信号光沿光纤长度变化的近似表达式,并将该近似分析解与解析解进行了比较。仿真结果表明上转换过程会降低随光纤传输的信号光和泵浦光,影响EDFA的放大性能,并得出了当掺铒浓度大于2.95×1025m-3时,上转换过程的影响不能忽略的结论。