A new approach of a limiting process for multi-dimensional flows

An enhanced Multi-dimensional Limiting Process (e-MLP) is developed for the accurate and efficient computation of multi-dimensional flows based on the Multi-dimensional Limiting Process (MLP). The new limiting process includes a distinguishing step and an enhanced multi-dimensional limiting process. First, the distinguishing step, which is independent of high order interpolation and flux evaluation, is newly introduced. It performs a multi-dimensional search of a discontinuity. The entire computational domain is then divided into continuous, linear discontinuous and nonlinear discontinuous regions. Second, limiting functions are appropriately switched according to the type of each region; in a continuous region, there are no limiting processes and only higher order accurate interpolation is performed. In linear discontinuous and nonlinear discontinuous regions, TVD criterion and MLP limiter are respectively used to remove oscillation. Hence, e-MLP has a number of advantages, as it incorporates useful features of MLP limiter such as multi-dimensional monotonicity and straightforward extensionality to higher order interpolation. It is applicable to local extrema and prevents excessive damping in a linear discontinuous region through application of appropriate limiting criteria. It is efficient because a limiting function is applied only to a discontinuous region. In addition, it is robust against shock instability due to the strict distinction of the computational domain and the use of regional information in a flux scheme as well as a high order interpolation scheme. This new limiting process was applied to numerous test cases including one-dimensional shock/sine wave interaction problem, oblique stationary contact discontinuity, isentropic vortex flow, high speed flow in a blunt body, planar shock/density bubble interaction, shock wave/vortex interaction and, particularly, magnetohydrodynamic (MHD) cloud-shock interaction problems. Through these tests, it was verified that e-MLP substantially enhances the accuracy and efficiency with both continuous and discontinuous multi-dimensional flows.     

A stable and high-order accurate conjugate heat transfer problem

This paper analyzes well-posedness and stability of a conjugate heat transfer problem in one space dimension. We study a model problem for heat transfer between a fluid and a solid. The energy method is used to derive boundary and interface conditions that make the continuous problem well-posed and the semi-discrete problem stable. The numerical scheme is implemented using 2nd-, 3rd- and 4th-order finite difference operators on Summation-By-Parts (SBP) form. The boundary and interface conditions are implemented weakly. We investigate the spectrum of the spatial discretization to determine which type of coupling that gives attractive convergence properties. The rate of convergence is verified using the method of manufactured solutions.     

海表白冠覆盖率的散射计资料反演研究

针对波浪破碎海况,采用考虑波浪破碎的KHCC03散射模式和散射计卫星遥感资料,对白冠覆盖率进行了反演研究,并得到了白冠覆盖率的反演公式.经检验,在10 m高度处风速取值为7 m/s—16 m/s时,利用文中给出的公式得到的白冠覆盖率反演结果与实测数据相比,优于以往其他计算公式的对应结果.这表明,利用散射计资料和波浪破碎的散射模式,可以反演白冠覆盖率.而且,这种方法不同于人们通常采用的被动遥感方法. 关键词： 白冠覆盖率 KHCC03散射模式 波浪破碎     

Effect of media-induced modification of travel rates on disease transmission in a multiple patch setting

A general SIS epidemic model is formulated that incorporates media-induced modification of travel rates. Basic local properties of solutions to the model are established. In particular, it is shown that the basic reproduction number does not involve parameters related to the effect of media on travel. The general model is subsequently specialised to two-patch models, with two different scenarios regarding patch population size. Qualitative analyses show that the basic reproduction number acts as a sharp threshold between disease persistence and extinction. The concept of uniform weak persistence is used to prove the existence of an endemic equilibrium and disease uniform strong persistence under a certain condition. Numerical investigations are carried out to gain insight into the analytically tractable and intractable cases, highlighting the importance of considering not only the basic reproduction number but also other measures of disease severity.     

High-order compact solution of the one-dimensional heat and advection–diffusion equations

In this work, we propose a high-order accurate method for solving the one-dimensional heat and advection–diffusion equations. We apply a compact finite difference approximation of fourth-order for discretizing spatial derivatives of these equations and the cubic C¹$C^1$-spline collocation method for the resulting linear system of ordinary differential equations. The cubic C¹$C^1$-spline collocation method is an A-stable method for time integration of parabolic equations. The proposed method has fourth-order accuracy in both space and time variables, i.e. this method is of order O(h⁴,k⁴)$O(h^4\text{,}{k}^4)$. Additional to high-order of accuracy, the proposed method is unconditionally stable which will be proved in this paper. Numerical results show that the compact finite difference approximation of fourth-order and the cubic C¹$C^1$-spline collocation method give an efficient method for solving the one-dimensional heat and advection–diffusion equations.     

A new high accuracy finite difference discretization for the solution of 2D nonlinear biharmonic equations using coupled approach

In this article, using coupled approach, we discuss fourth order finite difference approximation for the solution of two dimensional nonlinear biharmonic partial differential equations on a 9‐point compact stencil. The solutions of unknown variable and its Laplacian are obtained at each internal grid points. This discretization allows us to use the Dirichlet boundary conditions only and there is no need to discretize the derivative boundary conditions. We require only system of two equations to obtain the solution and its Laplacian. The proposed fourth order method is used to solve a set of test problems and produce high accuracy numerical solutions. © 2009 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 2010     

High order conservative Lagrangian schemes with Lax–Wendroff type time discretization for the compressible Euler equations

In this paper, we explore the Lax–Wendroff (LW) type time discretization as an alternative procedure to the high order Runge–Kutta time discretization adopted for the high order essentially non-oscillatory (ENO) Lagrangian schemes developed in 3 and 5. The LW time discretization is based on a Taylor expansion in time, coupled with a local Cauchy–Kowalewski procedure to utilize the partial differential equation (PDE) repeatedly to convert all time derivatives to spatial derivatives, and then to discretize these spatial derivatives based on high order ENO reconstruction. Extensive numerical examples are presented, for both the second-order spatial discretization using quadrilateral meshes [3] and third-order spatial discretization using curvilinear meshes [5]. Comparing with the Runge–Kutta time discretization procedure, an advantage of the LW time discretization is the apparent saving in computational cost and memory requirement, at least for the two-dimensional Euler equations that we have used in the numerical tests.     

基底显微结构对薄膜生长影响的Monte Carlo模拟研究

利用Monte Carlo方法研究了基底显微结构对薄膜生长的影响. 对不同显微结构基底上薄膜生长的初始阶段岛的形貌和尺寸与薄膜覆盖度和入射粒子沉积速率之间的关系进行了模拟和分析. 模型中考虑了粒子沉积、吸附粒子扩散和蒸发等过程. 结果表明,基底显微结构对薄膜生长具有明显影响. 当沉积温度为300K、沉积速率为0.005ML/s（Monolayer/second,简称ML/s）、覆盖度为0.05ML时,四方基底上薄膜生长呈现凝聚生长. 随着覆盖度增加,岛的尺寸变大,岛的数目减少. 而对于六方基底,当覆盖度从0.05ML变化到0.25ML时,薄膜生长经历了一个从分散生长过渡到分形生长的过程. 无论是四方还是六方基底,随着沉积速率的增加,岛的形貌由少数聚集型岛核分布状态向众多各自独立的离散型岛核分布状态过渡.     

基于模板匹配法的枪械射击准确度验收系统

枪械射击准确度是枪械验收的一项重要指标。为提高验收精度、降低验收成本,针对现有验收方法的不足,提出了光电成像并结合图像处理的验收方法。根据成像的图像特点和验收的实际要求,采用模板匹配算法识别代表瞄准线的十字分划,拟合枪管轴线并用十字分划表示,通过检验两十字之间的角度关系来进行射击准确度的验收。研究结果表明,采用模板匹配算法,十字分划的识别精度小于1个像素,射击准确度验收系统的精度约为0.26 mil。