Positive numerical solution for a nonarbitrage liquidity model using nonstandard finite difference schemes

In this article, we construct a numerical method based on a nonstandard finite difference scheme to solve numerically a nonarbitrage liquidity model with observable parameters for derivatives. This nonlinear model considers that the parameters involved are observable from order book data. The proposed numerical method use a exact difference scheme in the linear convection‐reaction term, and the spatial derivative is approximated using a nonstandard finite difference scheme. It is shown that the proposed numerical scheme preserves the positivity as well as stability and consistence. To illustrate the accuracy of the method, the numerical results are compared with those produced by other methods. © 2013 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 30: 210‐221, 2014     

Discrete artificial boundary conditions for nonlinear Schrödinger equations

In this work we construct and analyze discrete artificial boundary conditions (ABCs) for different finite difference schemes to solve nonlinear Schrödinger equations. These new discrete boundary conditions are motivated by the continuous ABCs recently obtained by the potential strategy of Szeftel. Since these new nonlinear ABCs are based on the discrete ABCs for the linear problem we first review the well-known results for the linear Schrödinger equation. We present our approach for a couple of finite difference schemes, including the Crank–Nicholson scheme, the Dùran–Sanz-Serna scheme, the DuFort–Frankel method and several split-step (fractional-step) methods such as the Lie splitting, the Strang splitting and the relaxation scheme of Besse. Finally, several numerical tests illustrate the accuracy and stability of our new discrete approach for the considered finite difference schemes.     

The Conservation and Convergence of Two Finite Difference Schemes for KdV Equations with Initial and Boundary Value Conditions

Korteweg-de Vries equation is a nonlinear evolutionary partial differential equation that is of third order in space. For the approximation to this equation with the initial and boundary value conditions using the finite difference method, the difficulty is how to construct matched finite difference schemes at all the inner grid points. In this paper, two finite difference schemes are constructed for the problem. The accuracy is second-order in time and first-order in space. The first scheme is a two-level nonlinear implicit finite difference scheme and the second one is a three-level linearized finite difference scheme. The Browder fixed point theorem is used to prove the existence of the nonlinear implicit finite difference scheme. The conservation, boundedness, stability, convergence of these schemes are discussed and analyzed by the energy method together with other techniques. The two-level nonlinear finite difference scheme is proved to be unconditionally convergent and the three-level linearized one is proved to be conditionally convergent. Some numerical examples illustrate the efficiency of the proposed finite difference schemes.     

An Unconditionally Stable Difference Approximation for a Class of Nonlinear Dispersive Equations

An unconditionally stable leapfrog finite difference scheme for a class of nonlinear dispersive equations is presented and analyzed. The solvability of the difference equation which is a tridiagonal circular linear system is discussed. Moreover, the convergence and stability of the difference scheme are also investigated by a standard argument so that more difficult priori estimations are avoided. Finally, numerical examples are given.     

Numerical investigation into nonlinear dynamic behavior of electrically-actuated clamped–clamped micro-beam with squeeze-film damping effect

A numerical investigation is performed into the nonlinear dynamic behavior of a clamped–clamped micro-beam actuated by a combined DC/AC voltage and subject to a squeeze-film damping effect. An analytical model based on a nonlinear deflection equation and a linearized Reynolds equation is proposed to describe the deflection of the micro-beam under the effects of the electrostatic actuating force. The deflection of the micro-beam is investigated under various actuating conditions by solving the analytical model using a hybrid numerical scheme comprising the differential transformation method and the finite difference approximation method. It is shown that the numerical results for the dynamic pull-in voltage of the clamped–clamped micro-beam deviate by no more than 2.04% from those presented in the literature based on the conventional finite difference scheme. The effects of the AC voltage amplitude, excitation frequency, residual stress, and ambient pressure on the center-point displacement of the micro-beam are systematically explored. Moreover, the actuation conditions which ensure the stability of the micro-beam are identified by means of phase portraits. Overall, the results presented in this study confirm that the hybrid numerical method provides an accurate means of analyzing the complex nonlinear behavior of common electrostatically-actuated microstructures.     

Compact difference scheme for two-dimensional fourth-order nonlinear hyperbolic equation

High-order compact finite difference method for solving the two-dimensional fourth-order nonlinear hyperbolic equation is considered in this article. In order to design an implicit compact finite difference scheme, the fourth-order equation is written as a system of two second-order equations by introducing the second-order spatial derivative as a new variable. The second-order spatial derivatives are approximated by the compact finite difference operators to obtain a fourth-order convergence. As well as, the second-order time derivative is approximated by the central difference method. Then, existence and uniqueness of numerical solution is given. The stability and convergence of the compact finite difference scheme are proved by the energy method. Numerical results are provided to verify the accuracy and efficiency of this scheme.     

微分本构粘弹性轴向运动弦线横向振动分析的差分法

给出了微分本构粘弹性轴向运动弦线横向振动数值仿真的一种差分法．文中建立了具有微分本构的粘弹性运动弦线的横向振动模型;通过对系统的控制方程和本构方程在不同的分数节点离散,得到一种新的差分方法．利用这一方法,弦线振动方程的数值计算过程可以交替地显式进行,且有较小的截断误差和好的数值稳定性．与通用的方法比较,新的方法计算简单、方便．文中利用方程的不变量检验了数值结果的可靠性,并利用这一方法给出了一类弦线模型的参数振动分析．     

New conservative difference schemes for a coupled nonlinear Schrödinger system

In this paper, two conservative difference schemes for solving a coupled nonlinear Schrödinger (CNLS) system are numerically analyzed. Firstly, a nonlinear implicit two-level finite difference scheme for CNLS system is studied, then a linear three-level difference scheme for CNLS system is presented. An induction argument and the discrete energy method are used to prove the second-order convergence and unconditional stability of the linear scheme. Numerical examples show the efficiency of the new scheme.     

Evaluating American put options on zero-coupon bonds by a penalty method

In this paper, American put options on zero-coupon bonds are priced under a single factor model of short-term rate. The linear complementarity problem of the option value is solved numerically by a penalty method, by which the problem is transformed into a nonlinear PDE by adding a power penalty term. The solution of the penalized problem converges to that of the original problem. A numerical scheme is established by using the finite volume method and the corresponding stability and convergence are discussed. Numerical results are presented to show the usefulness of the method.     

铁磁链Landau-Lifshitz方程的显式差分法

正如在研究流体动力学时,Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程起着十分重要的作用一样,在对于非平衡态磁学的研究中,描述连续铁磁体自旋场发展过程的 Ｌａｎｄａｎ－Ｌｉｆｓｈｉｔｚ方程［１］起着十分重要的作用［２］．一九九三年,美国和印度签署了一个大约 ２８０万美元的合作研究计划,在三年的时间里,对Ｌａｎｄａｕ－Ｌｉｆｓｈｉｔｚ铁磁链方程进行研究．在无阻尼的情况下,它为一完全可积的孤立子系统［３,４,５］。很多物理学家研究了它的孤立子解的存在性、逆散射方法以及相互碰撞［３,４,５］．关于解的存在性, Ａｌｏｎ…     

Conservative finite difference methods for fractional Schrödinger–Boussinesq equations and convergence analysis

In this paper, two conservative finite difference schemes for fractional Schrödinger–Boussinesq equations are formulated and investigated. The convergence of the nonlinear fully implicit scheme is established via discrete energy method, while the linear semi‐implicit scheme is analyzed by means of mathematical induction method. Our schemes are proved to preserve the total mass and energy in discrete level. The numerical results are given to confirm the theoretical analysis.     

Nonlinear scheme with high accuracy for nonlinear coupled parabolic-hyperbolic system

A nonlinear finite difference scheme with high accuracy is studied for a class of two-dimensional nonlinear coupled parabolic-hyperbolic system. Rigorous theoretical analysis is made for the stability and convergence properties of the scheme, which shows it is unconditionally stable and convergent with second order rate for both spatial and temporal variables. In the argument of theoretical results, difficulties arising from the nonlinearity and coupling between parabolic and hyperbolic equations are overcome, by an ingenious use of the method of energy estimation and inductive hypothesis reasoning. The reasoning method here differs from those used for linear implicit schemes, and can be widely applied to the studies of stability and convergence for a variety of nonlinear schemes for nonlinear PDE problems. Numerical tests verify the results of the theoretical analysis. Particularly it is shown that the scheme is more accurate and faster than a previous two-level nonlinear scheme with first order temporal accuracy.     

A numerical method for European Option Pricing with transaction costs nonlinear equation

This paper deals with the construction of a finite difference scheme and the numerical analysis of its solution for a nonlinear Black–Scholes partial differential equation modelling stock option pricing in the realistic case when transaction costs arising in the hedging of portfolios are taken into account. The analysed model is the Barles–Soner one for which an appropriate fully nonlinear numerical method has not still applied. After construction of the numerical solution, consistency and stability are studied and some illustrative examples are included.     

非线性Schrdinger方程的高精度守恒差分格式

本文首先分析线性Schrodinger方程一种高阶差分格式的构造方法,得到方程的耗散项.在此基础上对三次非线性Schrodinger方程,提出了一种精度为O(r2 h2)的差分格式,证明了该格式保持了连续方程的两个守恒量,且是收敛的与稳定的.并通过数值例子与已有隐格式进行了比较,结果表明,本文格式在计算量类似的情况下,提高了数值精度.     

A FULLY DISCRETE IMPLICIT-EXPLICIT FINITE ELEMENT METHOD FOR SOLVING THE FITZHUGH-NAGUMO MODEL

This work develops a fully discrete implicit-explicit finite element scheme for a parabolic-ordinary system with a nonlinear reaction term which is known as the FitzHugh-Nagumo model from physiology. The first-order backward Euler discretization for the time derivative, and an implicit-explicit discretization for the nonlinear reaction term are employed for the model, with a simple linearization technique used to make the process of solving equations more efficient. The stability and convergence of the fully discrete implicit-explicit finite element method are proved, which shows that the FitzHugh-Nagumo model is accurately solved and the trajectory of potential transmission is obtained. The numerical results are also reported to verify the convergence results and the stability of the proposed method.     

A nonlinear iteration method for solving a two-dimensional nonlinear coupled system of parabolic and hyperbolic equations

A nonlinear iteration method for solving a class of two-dimensional nonlinear coupled systems of parabolic and hyperbolic equations is studied. A simple iterative finite difference scheme is designed; the calculation complexity is reduced by decoupling the nonlinear system, and the precision is assured by timely evaluation updating. A strict theoretical analysis is carried out as regards the convergence and approximation properties of the iterative scheme, and the related stability and approximation properties of the nonlinear fully implicit finite difference (FIFD) scheme. The iterative algorithm has a linear constringent ratio; its solution gives a second-order spatial approximation and first-order temporal approximation to the real solution. The corresponding nonlinear FIFD scheme is stable and gives the same order of approximation. Numerical tests verify the results of the theoretical analysis. The discrete functional analysis and inductive hypothesis reasoning techniques used in this paper are helpful for overcoming difficulties arising from the nonlinearity and coupling and lead to a related theoretical analysis for nonlinear FI schemes.     

一类多维非线性反应扩散方程有限差分格式的稳定性问题(英文)

本文研究了一类多维线性反应扩散方程差分格式的稳定性.利用量未知元方法,建立了具有增量未知元的有限差分格式;然后利用非线性Galerkin方法,得到该差分格式的稳定性条件.通过对该格式的稳定性分析,说明和经典的差分格式的稳定性相比较,带有增量未知元的有限差分格式的稳定性得到了提高.     

IIM-based ADI finite difference scheme for nonlinear convection–diffusion equations with interfaces

Based on Li’s immersed interface method (IIM), an ADI-type finite difference scheme is proposed for solving two-dimensional nonlinear convection–diffusion interface problems on a fixed cartesian grid, which is unconditionally stable and converges with two-order accuracy in both time and space in maximum norm. Correction terms are added to the right-hand side of standard ADI scheme at irregular points. The nonlinear convection terms are treated by Adams–Bashforth method, without affecting the stability of difference schemes. A new method for computing the correction terms is developed, in which the Adams–Bashforth method is employed. Thus we can get an explicit approximation for the computation of corrections, when the jump condition is solution-dependent. Three numerical experiments are displayed and analyzed. The numerical results show good agreement with the exact solutions and confirm the convergence order.     

Second Order Unconditionally Stable and Convergent Linearized Scheme for a Fluid-Fluid Interaction Model

In this paper, a fully discrete finite element scheme with second-order temporal accuracy is proposed for a fluid-fluid interaction model, which consists of two Navier-Stokes equations coupled by a linear interface condition. The proposed fully discrete scheme is a combination of a mixed finite element approximation for spatial discretization, the second-order backward differentiation formula for temporal discretization, the second-order Gear's extrapolation approach for the interface terms and extrapolated treatments in linearization for the nonlinear terms. Moreover, the unconditional stability is established by rigorous analysis and error estimate for the fully discrete scheme is also derived. Finally, some numerical experiments are carried out to verify the theoretical results and illustrate the accuracy and efficiency of the proposed scheme.     

Strong convergence of a stochastic Rosenbrock-type scheme for the finite element discretization of semilinear SPDEs driven by multiplicative and additive noise

This paper aims to investigate the numerical approximation of a general second order parabolic stochastic partial differential equation(SPDE) driven by multiplicative and additive noise. Our main interest is on such SPDEs where the nonlinear part is stronger than the linear part, usually called stochastic dominated transport equations. Most standard numerical schemes lose their good stability properties on such equations, including the current linear implicit Euler method. We discretize the SPDE in space by the finite element method and propose a novel scheme called stochastic Rosenbrock-type scheme for temporal discretization. Our scheme is based on the local linearization of the semi-discrete problem obtained after space discretization and is more appropriate for such equations. We provide a strong convergence of the new fully discrete scheme toward the exact solution for multiplicative and additive noise and obtain optimal rates of convergence. Numerical experiments to sustain our theoretical results are provided.