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.     

Numerical study on flanging and spring-back of stamped thick metal plate

A quasi-three-dimensional shell element model, which can be effectively used to simulate the flanging and spring-back deformation, is introduced into the independently developed CAE software, KMAS. In this model, a double surface contact algorithm, which allows the gap between punch and die to change, and a spring-back treatment scheme based on finite element meshing are described. And then the flanging and spring-back deformations of the retractor's kickstand of a railcar made of stamped thick metal plate are numerically simulated. The simulation results of flanging deformation are compared with those of international commercial software, PAM-STAMP, and experimental ones. Finally, a predicting scheme of spring-back quantily for this problem is given. Project supported by the National Natural Science Foundation of China (No. 19832020) and the Ministry of Education of China.     

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