具有弱对称应力的线弹性方程的稳定化格式

利用稳定化方法讨论拉格朗日乘子法得到的具有弱对称应力的线弹性问题. 用线性元和分片常数分别逼近变分问题的应力和位移. 并通过添加稳定项$G_1(\cdot,\cdot)$, $G_2(\cdot,\cdot)$和$G_3(\cdot,\cdot)$ 使相应混合离散变分问题满足弱BB条件. 接着详细研究了变分问题的解与稳定混合有限元解之间的误差估计,最后用两个数值算例验证理论分析的有效性.     

A finite element formulation satisfying the discrete geometric conservation law based on averaged Jacobians

In this article, a new methodology for developing discrete geometric conservation law (DGCL) compliant formulations is presented. It is carried out in the context of the finite element method for general advective–diffusive systems on moving domains using an ALE scheme. There is an extensive literature about the impact of DGCL compliance on the stability and precision of time integration methods. In those articles, it has been proved that satisfying the DGCL is a necessary and sufficient condition for any ALE scheme to maintain on moving grids the nonlinear stability properties of its fixed‐grid counterpart. However, only a few works proposed a methodology for obtaining a compliant scheme. In this work, a DGCL compliant scheme based on an averaged ALE Jacobians formulation is obtained. This new formulation is applied to the θ family of time integration methods. In addition, an extension to the three‐point backward difference formula is given. With the aim to validate the averaged ALE Jacobians formulation, a set of numerical tests are performed. These tests include 2D and 3D diffusion problems with different mesh movements and the 2D compressible Navier–Stokes equations. Copyright © 2011 John Wiley & Sons, Ltd.     

Statistical stationary states for a two-layer quasi-geostrophic system

Existence of a family of locally invariant probability measures for large scale flows in enclosed temperate sea is proved. This model is extremely important for understanding the meso-scale phenomena in oceans. The techniques used are those developed by Albeverio and his collaborators.     

High‐order finite difference schemes for incompressible flows

This paper presents a new high‐order approach to the numerical solution of the incompressible Stokes and Navier–Stokes equations. The class of schemes developed is based upon a velocity–pressure–pressure gradient formulation, which allows: (i) high‐order finite difference stencils to be applied on non‐staggered grids; (ii) high‐order pressure gradient approximations to be made using standard Padé schemes, and (iii) a variety of boundary conditions to be incorporated in a natural manner. Results are presented in detail for a selection of two‐dimensional steady‐state test problems, using the fourth‐order scheme to demonstrate the accuracy and the robustness of the proposed methods. Furthermore, extensions to higher orders and time‐dependent problems are illustrated, whereas the extension to three‐dimensional problems is also discussed. Copyright © 2010 John Wiley & Sons, Ltd.     

Optimal transport with convex obstacle

We consider the Monge transportation problem when the cost is the squared geodesic distance around a convex obstacle. We show that there exists at least one—and in general infinitely many—optimal transport maps.     

The characteristic‐based split (CBS) meshfree method for free surface flow problems in ALE formulation

A semi‐implicit characteristic‐based split (CBS) meshfree algorithm in the arbitrary Lagrangian Eulerian (ALE) framework is proposed for the numerical solution of incompressible free surface flow problem in the paper. The algorithm is the extension of general CBS method which was initially introduced in finite element framework, this is due to the fact that CBS method not only can enhance the stability, but also avoid LBB condition when equal order basis function is used to approximate velocity and pressure variables. Meanwhile, a simple way for node update and node speed calculation is developed which is used to capture the free surface exactly. The numerical solutions are compared with available analytical and numerical solutions, which shows that the proposed method has better ability to simulate the free surface incompressible flow problem. Copyright © 2009 John Wiley & Sons, Ltd.     

Linear-programming extended formulations for the single-item lot-sizing problem with backlogging and constant capacity

Recently, several authors [8, 10] have argued for the use of extended formulations to tighten production planning models. In this work we present two linear-programming extended formulations of the constant-capacity lot-sizing problem with backlogging. The first one applies to the problem with a general cost function and has O(n³) variables and constraints. This improves on the more straightforward O(n⁴) Florian and Klein [2] type formulation. The second one applies when the costs satisfy the Wagner-Whitin property but it has O(n²) variables and O(n³) constraints. As a by-product, we positively answer an open question of Miller and Wolsey [4] about the tightness of an extended formulation for the continuous mixing set. This text presents research results of the Belgian Program on Interuniversity Poles of Attraction initiated by the Belgian State, Prime Minister's Office, Science Policy Programming. The research was carried out with financial support of the Growth Project G1RD-1999-00034 (LISCOS) of the European Community. The scientific responsibility is assumed by the author.     

Effect of numerical integration in the DGFEM for nonlinear convection‐diffusion problems

This paper is concerned with the effect of numerical integration applied to the discontinuous Galerkin finite element discretization of nonlinear convection‐diffusion problems in 2D. In the space semidiscretization the volume and line integrals are evaluated by numerical quadratures. Our goal is to estimate the error caused by the numerical integration and to show what numerical quadratures guarantee that the accuracy of the method with exact integration is preserved. © 2007 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2007     

Time-harmonic acoustic propagation in the presence of a shear flow

This work deals with the numerical simulation, by means of a finite element method, of the time-harmonic propagation of acoustic waves in a moving fluid, using the Galbrun equation instead of the classical linearized Euler equations. This work extends a previous study in the case of a uniform flow to the case of a shear flow. The additional difficulty comes from the interaction between the propagation of acoustic waves and the convection of vortices by the fluid. We have developed a numerical method based on the regularization of the equation which takes these two phenomena into account. Since it leads to a partially full matrix, we use an iterative algorithm to solve the linear system.     

面积,行列式,积分因子和Green公式

利用二阶行列式的几何意义是有向面积及积分因子的存在性给Green公式一个新的证明.尽管技术上走得远了些,但从概念上揭示了Green公式异常简明的几何意义,即Green公式只是面积的两种不同表达方式.同时这也蕴含了一个更深刻的哲学含义,即一般性隐含于特殊性(或特例)之中.