用瑞利激光雷达16年观测数据分析合肥地区大气密度分布特征

介绍了L625瑞利激光雷达系统结构以及基于瑞利散射理论探测大气分子数密度的原理.提出了反复迭代修正大气透过率的计算方法,并通过模拟仿真验证了该算法的可靠性.通过误差分析得到影响大气分子密度不确定度的主要因素为回波信号信噪比以及参考点处大气分子数密度值,给出了回波信号误差产生的主要来源以及参考点选取方法.最后,分析了激光雷达16年观测数据反演的结果,得到合肥地区大气分子数密度的月份以及年份分布状况,结果表明:中层大气分子数密度分布呈现明显的季节性分布特征,冬季分布稀疏,夏季分布密集,随年份分布则较为平稳.通过将统计平均得到的密度廓线与1976美国标准大气模式比对分析,发现由激光雷达观测反演得到的结果较模式值大,二者的比值在1.05～1.13之间.     

飞秒激光冲击AZ31B镁合金过程的数值模拟

采用有限元分析法对飞秒激光冲击AZ31B镁合金进行数值模拟,研究了激光冲击处理对镁合金变形过程的影响,分析了单脉冲激光冲击下材料内部的位移、动能、应力和应变的分布情况,得到了材料的瞬态速度和应变率变化过程.仿真结果表明,单脉冲飞秒激光冲击镁合金产生的塑性变形,可在材料表面形成微米级凹坑,中心点处最大位移为34μm,最大变形速度390m/s;在冲击初期,材料表面的应力和应变主要分布在冲击区域中心节点和边缘附近,并且得到镁合金的最大应力和最大应变率分别为955 MPa和1.8×106 s-1.研究结果能够为深入分析飞秒激光与镁合金作用时材料变形参量的变化规律提供数值理论依据.     

反对称分子动力学模型中的冻结概念（英文）

给出了反对称分子动力学模型(AMD)计算的50 Me V/nucleon112Sn+112Sn反应的分析结果。该研究是反对称分子动力学模型中统计冻结概念的部分研究结果。利用自洽法结合修正的Fisher模型,提取了发射源的温度和密度分别为T=(6.1±0.2)Me V,ρ/ρ0=0.69±0.03。通过与AMD模型计算的系统在时间演化过程中的最大密度比较,得出碎片发射源的密度远小于系统的最大密度。利用自洽法提取的温度和密度与35 Me V/nucleon的40Ca+40Ca反应系统及40 Me V/nucleon的64Zn+112Sn反应系统所提取的温度和密度非常接近。该结果表明反对称分子动力学模型中,系统在中等质量碎片形成时刻处于统计冻结体积。     

初始状态对自组织演化过程的影响

从实验上系统研究了在强电场下处于介质中的大量金属小球从不同初始状态通过自组织演化为分形的过程;用Sandbox方法定量分析了最终形成的稳定树枝状分形的维数。结果表明:不同初始状态形成分形的过程完全不同,但最终的树枝状分形的维数基本相同。分析认为,小球之间的相互作用存在临界作用距离,不同初始状态相邻小球之间的平均间距不同,因而演化过程不同;而维数相同正是系统耗散相同的体现。     

Decoupled two‐grid finite element method for the time‐dependent natural convection problem I: Spatial discretization

In this article, a decoupled two grid finite element method (FEM) is proposed and analyzed for the nonsteady natural convection problem using the coarse grid numerical solutions to decouple the nonlinear coupled terms, and the corresponding optimal error estimates are derived. Compared with the standard Galerkin FEM and the usual two‐grid FEM, our algorithm not only keeps good accuracy but also saves a lot of computational cost. Some numerical examples are provided to verify the performances of the decoupled two‐grid FEM. Both theoretical analysis and numerical experiments show the efficiency and effectiveness of the decoupled two‐grid FEM for the nonsteady natural convection problem. © 2015 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 31: 2135–2168, 2015     

A new method for numerical differentiation based on direct and inverse problems of partial differential equations

In this paper, we mainly study a numerical differentiation problem which aims to approximate the second order derivative of a single variable function from its noise data. By transforming the problem into a combination of direct and inverse problems of partial differential equations (heat conduction equations), a new method that we call the PDEs-based numerical differentiation method is proposed. By means of the finite element method and the Tikhonov regularization, implementations of the proposed PDEs-based method are presented with a posterior strategy for choosing regularization parameters. Numerical results show that the PDEs-based numerical differentiation method is highly feasible and stable with respect to data noise.     

Positive numerical splitting method for the Hull and White 2D Black–Scholes equation

We consider the locally one‐dimensional backward Euler splitting method to solve numerically the Hull and White problem for pricing European options with stochastic volatility in the presence of a mixed derivative term. We prove the first‐order convergence of the time‐splitting. The parabolic equation degenerates on the boundary x = 0 and we apply a fitted finite volume scheme to the equation to resolve the degeneracy and derive the fully discrete problem as we also investigate the discrete maximum principle. Numerical experiments illustrate the efficiency of our difference scheme. © 2014 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 31: 822–846, 2015     

A local parallel superconvergence method for the incompressible flow by coarsening projection

In this article, we investigate a local parallel superconvergence method by coarsening projection for the incompressible Stokes flow. The method is a combination of the local superconvergence technique and the given framework of local parallel method. For the smooth subdomains, the local superconvergence method is applied in a higher order finite dimensional space corresponding to an appropriate coarse mesh on interior domain. Moreover, a useful and flexible local parallel method is designed to obtain the local parallel superconvergence results of presented method, which offset theoretical limitation of the model without the smoothness of the exact solution and a priori regularity of the underlying problem over the whole domain. © 2014 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 31: 1209–1223, 2015