Development and validation of a SUPG finite element scheme for the compressible Navier–Stokes equations using a modified inviscid flux discretization

This paper considers the streamline‐upwind Petrov–Galerkin (SUPG) method applied to the unsteady compressible Navier–Stokes equations in conservation‐variable form. The spatial discretization, including a modified approach for interpolating the inviscid flux terms in the SUPG finite element formulation, and the second‐order accurate time discretization are presented. The numerical method is discussed in detail. The performance of the algorithm is then investigated by considering inviscid flow past a circular cylinder. Validation of the finite element formulation via comparisons with experimental data for high‐Mach number perfect gas laminar flows is presented, with a specific focus on comparisons with experimentally measured skin friction and convective heat transfer on a 15° compression ramp. Copyright © 2007 John Wiley & Sons, Ltd.     

Mortar element methods for parabolic problems

In this article a standard mortar finite element method and a mortar element method with Lagrange multiplier are used for spatial discretization of a class of parabolic initial‐boundary value problems. Optimal error estimates in L^∞(L²) and L^∞(H¹)‐norms for semidiscrete methods for both the cases are established. The key feature that we have adopted here is to introduce a modified elliptic projection. In the standard mortar element method, a completely discrete scheme using backward Euler scheme is discussed and optimal error estimates are derived. The results of numerical experiments support the theoretical results obtained in this article. © 2008 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 2008     

Heavy tailed durations of regional rainfall

Durations of rain events and drought events over a given region provide important information about the water resources of the region. Of particular interest is the shape of upper tails of the probability distributions of such durations. Recent research suggests that the underlying probability distributions of such durations have heavy tails of hyperbolic type, across a wide range of spatial scales from 2 km to 120 km. These findings are based on radar measurements of spatially averaged rain rate (SARR) over a tropical oceanic region. The present work performs a nonparametric inference on the Pareto tail-index of wet and dry durations at each of those spatial scales, based on the same data, and compares it with conclusions based on the classical Hill estimator. The results are compared and discussed. The authors express sincere thanks to the Mathematisches Forschungsinstitut Oberwolfach (MFO) for facilitating their collaboration under a “Research in Pairs” project hosted at MFO during March 5–25, 2006. The research of the second and third authors was supported by the project LC06024.     

最高阶元素个数为40的非可解群的分类

设$\varphi$为群${\rm Aut}(N)$的同态,记$H_\varphi\times N$为群$N$借助于群$H$的半直积.设$G$为有限不可解群,本文证明: 若$G$中最高阶元素个数为40, 则$G$同构于下列群之一:(1)~$Z_{4\varphi}\times A_5$,\,${\rm ker}\varphi=Z_2$; (2)~$D_{8\varphi}\times A_5,\,{\rm ker}\varphi=Z_2\times Z_2$; (3)~$G/N=S_5$, $N=Z(G)=Z_2$; (4)~$G/N=S_5$, $N=Z_2\times Z_2,\,N\cap Z(G)=Z_2$.     

Hochschild and ordinary cohomology rings of small categories

Let C be a small category and k a field. There are two interesting mathematical subjects: the category algebra kC and the classifying space |C|=BC. We study the ring homomorphism HH^∗(kC)→H^∗(|C|,k) and prove it is split surjective, using the factorization category of Quillen [D. Quillen, Higher algebraic K-theory I, in: Lecture Notes in Math., vol. 341, Springer-Verlag, Berlin, 1973, pp. 85-147] and certain techniques from functor cohomology theory. This generalizes the well-known theorems for groups and posets. Based on this result, we construct a seven-dimensional category algebra whose Hochschild cohomology ring modulo nilpotents is not finitely generated, disproving a conjecture of Snashall and Solberg [N. Snashall, Ø. Solberg, Support varieties and Hochschild cohomology rings, Proc. London Math. Soc. 88 (3) (2004) 705-732].     

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     

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.     

负离子源引出系统束光学的数值模拟

根据在表面-等离子体型负离子源中,负离子主要由转换电极表面产生这一特点,参考正离子源引出系统数值计算程序,建立了负离子源引出系统的数值计算模型和计算程序;对表面-等离子体型桶式负离子源引出系统束光学的性质进行了数值模拟。对计算结果的检验和分析表明,这个计算模型和程序反映了负离子源引出系统束光学的基本特性。