Rational approximation to trigonometric operators

We consider the approximation of trigonometric operator functions that arise in the numerical solution of wave equations by trigonometric integrators. It is well known that Krylov subspace methods for matrix functions without exponential decay show superlinear convergence behavior if the number of steps is larger than the norm of the operator. Thus, Krylov approximations may fail to converge for unbounded operators. In this paper, we propose and analyze a rational Krylov subspace method which converges not only for finite element or finite difference approximations to differential operators but even for abstract, unbounded operators. In contrast to standard Krylov methods, the convergence will be independent of the norm of the operator and thus of its spatial discretization. We will discuss efficient implementations for finite element discretizations and illustrate our analysis with numerical experiments. AMS subject classification (2000)  65F10, 65L60, 65M60, 65N22     

Bottlebrush‐shaped copolymers with cellulose diacetate backbone by a combination of ring opening polymerization and ATRP

Graft copolymers with cellulose diacetate (CDA) backbone and both the poly(ε‐caprolactone) and polystyrene, or poly(butyl acrylate) or PMMA grafts were prepared by two‐step process. First, ε‐caprolactone (CL) was polymerized by ring‐opening polymerization (ROP) initiated with CDA, partly funcionalized with 2‐bromo‐isobutyryl groups (degree of functionalization was 0.5). The p(CDA‐g‐CL) copolymers were used in the second step as polyfunctional macroinitiators of ATRP of the vinyl monomer, giving densely grafted copolymers with polyester and PSt, or PBuA, or PMMA grafts. The prepared copolymers were characterized by SEC, some of them also by FTIR spectroscopy and atomic force microscopy (AFM). © 2007 Wiley Periodicals, Inc. J Polym Sci Part A: Polym Chem 46: 564–573, 2008     

High‐order 𝒞1 finite‐element interpolating schemes—Part I: Semi‐Lagrangian linear advection

This paper is devoted to the development of accurate high‐order interpolating schemes for semi‐Lagrangian advection. The characteristic‐Galerkin formulation is obtained by using a semi‐Lagrangian temporal discretization of the total derivative. The semi‐Lagrangian method requires high‐order interpolators for accuracy. A class of ??¹ finite‐element interpolating schemes is developed and two semi‐Lagrangian methods are considered by tracking the feet of the characteristic lines either from the interpolation or from the integration nodes. Numerical stability and analytical results quantifying the amount of artificial viscosity induced by the two methods are presented in the case of the one‐dimensional linear advection equation, based on the modified equation approach. Results of test problems to simulate the linear advection of a cosine hill illustrate the performance of the proposed approach. Copyright © 2007 John Wiley & Sons, Ltd.     

A stabilized incremental projection scheme for the incompressible Navier–Stokes equations

It is well known that any spatial discretization of the saddle‐point Stokes problem should satisfy the Ladyzhenskaya–Brezzi–Babuska (LBB) stability condition in order to prevent the appearance of spurious pressure modes. Particularly, if an equal‐order approximation is applied, the Schur complement (or, as called some times, the Uzawa matrix) of the pressure system has a non‐trivial null space that gives rise to such modes. An idea in the past was that all the schemes that solve a Poisson equation for the pressure rather than the Uzawa pressure equation (splitting/projection methods) should overcome this difficulty; this idea was wrong. There is numerical evidence that at least the so‐called incremental projection scheme still suffers from spurious pressure oscillations if an equal‐order approximation is applied. The present paper tries to distinguish which projection requires LBB‐compliant approximation and which does not. Moreover, a stabilized version of the incremental projection scheme is derived. Proper bounds for the stabilization parameter are also given. The numerical results show that the stabilized scheme does indeed achieve second‐order accuracy and does not produce spurious (node to node) pressure oscillations. Copyright © 2001 John Wiley & Sons, Ltd.     

A hybrid vortex method for the simulation of three‐dimensional flows

This paper presents an integral vorticity method for solving three‐dimensional Navier–Stokes equations. A finite volume scheme is implemented to solve the vorticity transport equation, which is discretized on a structured hexahedral mesh. A vortex sheet algorithm is used to enforce the no‐slip boundary condition through a vorticity flux at the boundary. The Biot–Savart integral is evaluated to compute the velocity field, in conjunction with a fast algorithm based on multipole expansion. This method is applied to the simulation of uniform flow past a sphere. Copyright © 2007 John Wiley & Sons, Ltd.     

关于3阶Carmichael数

朱文余和孙琦（见《数学进展》，2004，33（4）：505-507）提出了关于3阶Carmichael数的三个问题，我们（见《四川大学学报（自然科学版）》，2006，43（6）：1197-1201）肯定地回答了问题1．本文模仿Howe的寻找严格2阶Carmichael数（见Mathematics of Computation，2000，69（232）：1711—1719）的方法，提出寻找满足某种条件的3阶Carmichael数的方法，并用这种方法确实找到了几百个这样的数，因而完全肯定地回答了问题2．     

A partial differential equation connected to option pricing with stochastic volatility: Regularity results and discretization

This paper completes a previous work on a Black and Scholes equation with stochastic volatility. This is a degenerate parabolic equation, which gives the price of a European option as a function of the time, of the price of the underlying asset, and of the volatility, when the volatility is a function of a mean reverting Orstein-Uhlenbeck process, possibly correlated with the underlying asset. The analysis involves weighted Sobolev spaces. We give a characterization of the domain of the operator, which permits us to use results from the theory of semigroups. We then study a related model elliptic problem and propose a finite element method with a regular mesh with respect to the intrinsic metric associated with the degenerate operator. For the error estimate, we need to prove an approximation result.

     

热传导型半导体瞬态问题的特征有限体积方法及分析

在这篇文章中,将特征线方法与有限体积方法相结合,构造了热传导型半导体瞬态问题的全离散特征有限体积方法(CFVM).在一般的条件下进行理论分析,得到了次优阶L2误差估计结果.     

Wave reflection in slightly compressible, finitely deformed elastic media

Summary  In this paper, the reflection of a plane wave at an incrementally traction-free boundary of a half-space composed of nearly incompressible elastic material is considered. It is shown that two distinct cases exist, these being dependent on the underlying primary deformation. In the first case, the appropriate slowness sections are each approximately elliptical, and the corresponding reflection phenomena closely mirrors that associated with the corresponding linear isotropic theory. Specifically, an angular range of direction of incident wave exists, for which both a quasi-longitudinal and quasi-shear wave are reflected, the former being replaced by a surface wave outside this angular range. In the second case, the outer slowness section is re-enrant and, in addition to the scenarios previously mentioned, it is possible for two quasi-shear waves to be reflected. Numerical illustrations of reflection coefficients are presented in respect of a modified Varga material and the case of increasing bulk modulus is investigated. Received 17 January 2000; accepted for publication 22 February 2000