A modified predictor–corrector scheme for the two-dimensional sine-Gordon equation

A three-time level finite-difference scheme based on a fourth order in time and second order in space approximation has been proposed for the numerical solution of the nonlinear two-dimensional sine-Gordon equation. The method, which is analysed for local truncation error and stability, leads to the solution of a nonlinear system. To avoid solving it, a predictor–corrector scheme using as predictor a second-order explicit scheme is proposed. The procedure of the corrector has been modified by considering as known the already evaluated corrected values instead of the predictor ones. This modified scheme has been tested on the line and circular ring soliton and the numerical experiments have proved that there is an improvement in the accuracy over the standard predictor–corrector implementation. This research was co-funded by E.U. (75%) and by the Greek Government (25%).     

A second order numerical scheme for the solution of the one-dimensional Boussinesq equation

A predictor–corrector (P-C) scheme is applied successfully to a nonlinear method arising from the use of rational approximants to the matrix-exponential term in a three-time level recurrence relation. The resulting nonlinear finite-difference scheme, which is analyzed for local truncation error and stability, is solved using a P-C scheme, in which the predictor and the corrector are explicit schemes of order 2. This scheme is accelerated by using a modification (MPC) in which the already evaluated values are used for the corrector. The behaviour of the P-C/MPC schemes is tested numerically on the Boussinesq equation already known from the bibliography free of boundary conditions. The numerical results are derived for both the bad and the good Boussinesq equation and conclusions from the relevant known results are derived.      

Sixth algebraic order trigonometrically fitted predictor–corrector methods for the numerical solution of the radial Schrödinger equation

Our new trigonometrically fitted predictor–corrector (P–C) schemes presented here are based on the well known Adams–Bashforth–Moulton methods: the predictor is based on the fifth order Adams–Bashforth scheme and the corrector on the sixth order Adams–Moulton scheme. We tested the efficiency of our newly developed schemes against well known methods, with excellent results. The numerical experiments showed that at least one of our schemes is noticeably more efficient compared to other methods, some of which are specially designed for this type of problem. It is also worth mentioning that this is the first time that sixth algebraic order trigonometrically fitted Adams–Bashforth–Moulton P–C schemes are used to efficiently solve the radial Schrödinger equation.Active Member of the European Academy of Sciences and Arts     

A finite volume cell‐centered Lagrangian hydrodynamics approach for solids in general unstructured grids

A finite volume cell‐centered Lagrangian hydrodynamics approach, formulated in Cartesian frame, is presented for solving elasto‐plastic response of solids in general unstructured grids. Because solid materials can sustain significant shear deformation, evolution equations for stress and strain fields are solved in addition to mass, momentum, and energy conservation laws. The total stress is split into deviatoric shear stress and dilatational components. The dilatational response of the material is modeled using the Mie‐Grüneisen equation of state. A predicted trial elastic deviatoric stress state is evolved assuming a pure elastic deformation in accordance with the hypo‐elastic stress‐strain relation. The evolution equations are advanced in time by constructing vertex velocity and corner traction force vectors using multi‐dimensional Riemann solutions erected at mesh vertices. Conservation of momentum and total energy along with the increase in entropy principle are invoked for computing these quantities at the vertices. Final state of deviatoric stress is effected via radial return algorithm based on the J‐2 von Mises yield condition. The scheme presented in this work is second‐order accurate both in space and time. The suitability of the scheme is evinced by solving one‐ and two‐dimensional benchmark problems both in structured grids and in unstructured grids with polygonal cells. Copyright © 2013 John Wiley & Sons, Ltd.     

A numerical analysis of a class of generalized Boussinesq‐type equations using continuous/discontinuous FEM

An improved class of Boussinesq systems of an arbitrary order using a wave surface elevation and velocity potential formulation is derived. Dissipative effects and wave generation due to a time‐dependent varying seabed are included. Thus, high‐order source functions are considered. For the reduction of the system order and maintenance of some dispersive characteristics of the higher‐order models, an extra O(μ²ⁿ⁺²) term (n ∈ ?) is included in the velocity potential expansion. We introduce a nonlocal continuous/discontinuous Galerkin FEM with inner penalty terms to calculate the numerical solutions of the improved fourth‐order models. The discretization of the spatial variables is made using continuous P₂ Lagrange elements. A predictor‐corrector scheme with an initialization given by an explicit Runge–Kutta method is also used for the time‐variable integration. Moreover, a CFL‐type condition is deduced for the linear problem with a constant bathymetry. To demonstrate the applicability of the model, we considered several test cases. Improved stability is achieved. Copyright © 2011 John Wiley & Sons, Ltd.     

双贴合无光焦度校正板在双球面镜次镜后的系统的像差特性研究

讨论了双贴合无光焦度校正板放在次镜后的双反射球面系统.根据像差理论分析了消像差的条件,研究了主镜的相对孔径变化、校正板在次镜后的位置变化以及无光焦度校正板的正负透镜的不同组合形式对系统高级残余像差的影响.结果表明：当主镜的相对孔径增加或校正板靠近次镜时,系统的残余像差增加;无论校正板的第一个透镜是正透镜还是负透镜,系统在某一光焦度处残余像差最小,并且靠近次镜的透镜最好采用负透镜.     

大型Schmidt望远镜光学系统调整方法

安装于紫金山天文台盱眙观测站用以观测近地天体的大型Schmidt望远镜,采用了平场Schmidt光学系统,改正镜通光口径为1m,球面反射镜口径为1.2m,焦距为1.8m,接收器用了4K×4K的高灵敏度CCD。本文叙述了该仪器光轴调整的方法及调整后观测获得的初步结果。叙述了望远镜光轴校正方法、校正结果、CCD靶面的调整及望远镜极轴高度、方位的调整及相应的照片。经实际观测的结果是:露光1s可拍到18等星;露光4s可拍到19.3等星;露光20s可拍摄到21.2等星。     

合肥光源直线加速器BPM系统的研制和测量

合肥光源直线加速器新改造的BPM系统已经安装, 并做了在线的束流位置测量实验. 本文介绍该BPM系统的结构、作为系统核心部分的信号处理模块以及离线标定的实验结果和在线应用的实验结果, 表明该BPM系统具有不小于40dBm的动态范围, 不大于100μm的束流位置分辨率, 水平方向和垂直方向分别具有1.61dB/mm和1.51dB/mm的灵敏度. 校正铁强度与BPM读数关系的实验, 表明该BPM系统能正确地反映束流的位置变化情况, 并且确定了校正铁在水平方向和垂直方向的校正系数分别为0.46mm/A和0.58mm/A, 校正铁设置为默认值时该BPM处的束流中心水平位置和垂直位置分别为0.83mm和－0.57mm.     

Solving the logit-based stochastic user equilibrium problem with elastic demand based on the extended traffic network model

This paper proposes a novel extended traffic network model to solve the logit-based stochastic user equilibrium (SUE) problem with elastic demand. In this model, an extended traffic network is established by properly adding dummy nodes and links to the original traffic network. Based on the extended traffic network, the logit-based SUE problem with elastic demand is transformed to the SUE problem with fixed demand. Such problem is then further converted to a linearly constrained convex programming and addressed by a predictor–corrector interior point algorithm with polynomial complexity. A numerical example is provided to compare the proposed model with the method of successive averages (MSA). The numerical results indicate that the proposed model is more efficient and has a better convergence than the MSA.     

Robust spectral method for numerical valuation of european options under Merton's jump‐diffusion model

We propose a novel numerical method based on rational spectral collocation and Clenshaw–Curtis quadrature methods together with the “” transformation for pricing European vanilla and butterfly spread options under Merton's jump‐diffusion model. Under certain assumptions, such model leads to a partial integro‐differential equation (PIDE). The differential and integral parts of the PIDE are approximated by the rational spectral collocation and the Clenshaw–Curtis quadrature methods, respectively. The application of spectral collocation method to the PIDE leads to a system of ordinary differential equations, which is solved using the implicit–explicit predictor–corrector (IMEX‐PC) schemes in which the diffusion term is integrated implicitly, whereas the convolution integral, reaction, advection terms are integrated explicitly. Numerical experiments illustrate that our approach is highly accurate and efficient for pricing financial options.Copyright © 2014 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 30: 1169–1188, 2014