带插值的自适应网格重构算法求解带移动热源的反应扩散方程的理论分析

本文研究一个带插值的网格重构算法求解一类带移动热源的反应扩散方程. 算法包括两步: 第一步是用旧时间网层上的计算解计算新时间层上的空间网格; 第二步是使用有限差分方法在新时间层 空间网格上离散方程, 并且将旧时间层上计算解的插值作为初始值. 对于时间, 我们获得了一阶收敛结果. 对于空间, 我们证明了使用线性插值算法的一阶收敛性和使用二次插值算法的二阶收敛性. 数值例子肯定了本文的理论结果.     

A high order moving boundary treatment for compressible inviscid flows

We develop a high order numerical boundary condition for compressible inviscid flows involving complex moving geometries. It is based on finite difference methods on fixed Cartesian meshes which pose a challenge that the moving boundaries intersect the grid lines in an arbitrary fashion. Our method is an extension of the so-called inverse Lax–Wendroff procedure proposed in [17] for conservation laws in static geometries. This procedure helps us obtain normal spatial derivatives at inflow boundaries from Lagrangian time derivatives and tangential derivatives by repeated use of the Euler equations. Together with high order extrapolation at outflow boundaries, we can impose accurate values of ghost points near the boundaries by a Taylor expansion. To maintain high order accuracy in time, we need some special time matching technique at the two intermediate Runge–Kutta stages. Numerical examples in one and two dimensions show that our boundary treatment is high order accurate for problems with smooth solutions. Our method also performs well for problems involving interactions between shocks and moving rigid bodies.     

分层流体中移动动量源生成准二维偶极子涡街特性实验

采用水平移动射流方法模拟移动动量源,实验研究了该移动动量源在密度分层流体中生成准二维偶极子涡街的机理,分析了偶极子涡街的演化特性. 在系列实验基础上,获得了移动动量源在密度分层流体中能够演化为偶极子涡街的(Re,Fr)组合条件. 对不同的Re取值,获得了偶极子涡街无因次形成时间及其无因次涡街平均波长倒数与Fr之间的相关关系,表明它们不依赖于Re数,而与Fr近似为幂指数关系. 关键词： 分层流体 移动动量源 偶极子涡街 涡街形成时间     

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

从平面几何中曲线之间的相切关系不依赖于坐标轴的选取这一基本事实去看数学分析中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模型相比较，新模型模拟的车流量较高，说明新模型减少了车流的总体停滞状态. 关键词： 交通流 元胞自动机 行车状态 反馈规则     

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

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