从几何的角度看微分中值定理

从平面几何中曲线之间的相切关系不依赖于坐标轴的选取这一基本事实去看数学分析中Rolle中值定理、Lagrange中值定理以及Cauchy中值定理.     

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.     

Numerical models for vehicle exhaust dispersion in complex urban areas

Two numerical models are presented for predicting vehicle exhaust dispersion in complex urban areas with or without the wind field. The models not only reflect the effect of building and street canyon configuration on the pollutant propagation, but also are able to predict the turbulent energy produced by moving vehicles on the road. In particular, in the discrete model, turbulent energy and pollutant concentration produced by each vehicle are dynamically described in the Lagrangian method. The pollutant propagation is calculated with the advection–diffusion equation. The Reynolds averaged Navier–Stokes equations are numerically solved for the wind flows. The movement and heat release rate of the vehicles are treated as sources of the turbulent energy equation for the computation of turbulent energy produced by the moving vehicles. This paper reports the detailed implementation of the models. Four typical numerical tests were carried out to represent the performance of the proposed numerical models. Copyright © 2010 John Wiley & Sons, Ltd.     

An evaluation of a 3D free‐energy‐based lattice Boltzmann model for multiphase flows with large density ratio

In this paper, the 3D Navier–Stokes (N–S) equation and Cahn–Hilliard (C–H) equations were solved using a free‐energy‐based lattice Boltzmann (LB) model. In this model, a LB equation with a D3Q19 velocity model is used to recover continuity and N–S equations while another LB equation with D3Q7 velocity model for solving C–H equation (Int. J. Numer. Meth. Fluids, 2008; 56 :1653–1671) is applied to solve the 3D C–H equation. To avoid the excessive use of computational resources, a moving reference frame is adopted to allow long‐time simulation of a bubble rising. How to handle the inlet/outlet and moving‐wall boundary conditions are suggested. These boundary conditions are simple and easy for implementation. This model's performance on two‐phase flows was investigated and the mass conservation of this model was evaluated. The model is validated by its application to simulate the 3D air bubble rising in viscous liquid (density ratio is 1000). Good agreement was obtained between the present numerical results and experimental results when Re is small. However, for high‐Re cases, the mass conservation seems not so good as the low‐Re case. Copyright © 2009 John Wiley & Sons, Ltd.     

Simulation of two‐phase flow with moving immersed boundaries

A two‐dimensional multi‐phase model for immiscible binary fluid flow including moving immersed objects is presented. The fluid motion is described by the incompressible Navier–Stokes equation coupled with a phase‐field model based on van der Waals' free energy density and the Cahn–Hilliard equation. A new phase‐field boundary condition was implemented with minimization of the free energy in a direct way, to specifically improve the physical behavior of the contact line dynamics for moving immersed objects. Numerical stability and execution time were significantly improved by the use of the new boundary condition. Convergence toward the analytical solution was demonstrated for equilibrium contact angle, the Lucas–Washburn theory and Stefan's problem. The proposed model may be used for multi‐phase flow problems with moving boundaries of complex geometry, such as the penetration of fluid into a deformable, porous medium. Copyright © 2010 John Wiley & Sons, Ltd.     

考虑行车状态的一维元胞自动机交通流模型

在Nagel Schrekenberg单车道元胞自动机交通流模型（简称NS模型）的基础上，考虑车辆之间的相对运动薛郁等提出了一种改进的单车道元胞自动机交通流模型（简称改进的NS模型）.通过两种情况列出了改进的NS模型存在不尽周严的地方，随之在新模型中引入了行车状态 变量和反馈规则，从而控制车辆出现倒车和刹车过急等现象.通过计算机对新模型进行模拟 ，发现减速概率和车流密度对车流状态的演化影响很大，当减速概率高（如道路条件差）时 ，即使车流密度低，车流也会出现局部堵塞状态；而当减速概率一定时，随着车流密度增加 ，车流的运动相与堵塞相发生了全局性的交替出现，此时类似于波的波峰和波谷的传播.与 改进的NS模型相比较，新模型模拟的车流量较高，说明新模型减少了车流的总体停滞状态. 关键词： 交通流 元胞自动机 行车状态 反馈规则     

与运动原子相互作用的场的压缩效应

采用量子光学中时间演化算符方法，研究了与运动原子相互作用的场的压缩效应，揭示了光场的初始平均光子数和场模结构参数对光场压缩特性的影响，结果表明：适当选择系统参数，可获得在时间上持续压缩的压缩光。     

Bootstrapping Autoregressive and Moving Average Parameter Estimates of Infinite Order Vector Autoregressive Processes

We consider anr-dimensional multivariate time series {y_t, tZ} which is generated by an infinite order vector autoregressive process. We show that a bootstrap procedure which works by generating time series replicates via an estimated finitek-order vector autoregressive process (k→∞ at an appropriate rate with the sample size) gives asymptotically valid approximations to the joint distribution of the growing set of estimated autoregressive coefficients and to the corresponding set of estimated moving average coefficients (impuls responses).     

Boundary element-free method for elastodynamics

The moving least-square approximation is discussed first. Sometimes the method can form an ill-conditioned equation system, and thus the solution cannot be obtained correctly. A Hilbert space is presented on which an orthogonal function system mixed a weight function is defined. Next the improved moving least-square approximation is discussed in detail. The improved method has higher computational efficiency and precision than the old method, and cannot form an ill-conditioned equation system. A boundary element-free method (BEFM) for elastodynamics problems is presented by combining the boundary integral equation method for elastodynamics and the improved moving least-square approximation. The boundary element-free method is a meshless method of boundary integral equation and is a direct numerical method compared with others, in which the basic unknowns are the real solutions of the nodal variables and the boundary conditions can be applied easily. The boundary element-free method has a higher computational efficiency and precision. In addition, the numerical procedure of the boundary element-free method for elastodynamics problems is presented in this paper. Finally, some numerical examples are given.     

层流中固体颗粒运动方程的初值问题和稳定性分析

层流中固体颗粒运动方程的初值问题解的存在唯一性表明悬浮颗粒在圆管层流人口段剪切区中的运动，由进入该区的颗粒初始状态（初始位置和初始速度）确定，其运动轨迹也由该初始状态唯一确定。层流中固体颗粒运动方程稳定性分析表明达到稳定状态的颗粒必定满足条件：up＝uf（∞，c），up＝0，Re（α／R）2（c／R）＜2．8728．