On the absolute instability of semi-implicit schemes for hydrostatic models

The dependence of the linear stability of two-time-level finite-difference semi-implicit schemes on the choice of reference temperature profile is studied. Particular vertical profiles of the temperature are considered to derive analytical conditions of stability. Analysis is made for general form of different model parameters such as the number of vertical levels and their distribution, the time step size, and the values of the viscosity coefficients. The derived conditions of stability are more restrictive than those for three-time-level schemes, but obtained necessary and sufficient condition for constant vertical lapse rates of the temperature has the form frequently applied to three-time-level schemes: the basic temperature profile should be warmer than the actual one. Performed numerical experiments show that the last restriction is neither necessary nor sufficient condition of stability for general temperature profiles.     

Penalty methods for the numerical solution of American multi-asset option problems

We derive and analyze a penalty method for solving American multi-asset option problems. A small, non-linear penalty term is added to the Black–Scholes equation. This approach gives a fixed solution domain, removing the free and moving boundary imposed by the early exercise feature of the contract. Explicit, implicit and semi-implicit finite difference schemes are derived, and in the case of independent assets, we prove that the approximate option prices satisfy some basic properties of the American option problem. Several numerical experiments are carried out in order to investigate the performance of the schemes. We give examples indicating that our results are sharp. Finally, the experiments indicate that in the case of correlated underlying assets, the same properties are valid as in the independent case.     

Weak second order S-ROCK methods for Stratonovich stochastic differential equations

It is well known that the numerical solution of stiff stochastic ordinary differential equations leads to a step size reduction when explicit methods are used. This has led to a plethora of implicit or semi-implicit methods with a wide variety of stability properties. However, for stiff stochastic problems in which the eigenvalues of a drift term lie near the negative real axis, such as those arising from stochastic partial differential equations, explicit methods with extended stability regions can be very effective. In the present paper our aim is to derive explicit Runge–Kutta schemes for non-commutative Stratonovich stochastic differential equations, which are of weak order two and which have large stability regions. This will be achieved by the use of a technique in Chebyshev methods for ordinary differential equations.     

Revisit of Semi-Implicit Schemes for Phase-Field Equations

It is a very common practice to use semi-implicit schemes in various computations, which treat selected linear terms implicitly and the nonlinear terms explicitly. For phase-field equations, the principal elliptic operator is treated implicitly to reduce the associated stability constraints while the nonlinear terms are still treated explicitly to avoid the expensive process of solving nonlinear equations at each time step. However, very few recent numerical analysis is relevant to semi-implicit schemes, while "stabilized" schemes have become very popular. In this work, we will consider semi-implicit schemes for the Allen-Cahn equation with $general$ $potential$ function. It will be demonstrated that the maximum principle is valid and the energy stability also holds for the numerical solutions. This paper extends the result of Tang & Yang (J. Comput. Math., 34(5) (2016), pp. 471-481), which studies the semi-implicit scheme for the Allen-Cahn equation with $polynomial$ $potentials$.     

On an Unconditionally Stable Scheme for the Unsteady Navier-Stokes Equations

Theoretical time step constraints of semi-implicit schemes are known to be more restrictive than should be in practice. We intend to alleviate the constraints with more smoothness assumptions on the solutions. By introducing a new scheme with modification on the treatment of the nonlinear term, we are able to prove that the scheme is unconditionally stable and convergent. Furthermore, we show that the modified scheme and the original semi-implicit one are equivalent under a weak condition on the time step and the number of space discretization points.     

Improved stability techniques for the solution of Navier-Stokes equations

When solving the Navier-Stokes equations for transient incompressible viscous flow problems, one normally encounters a decrease in numerical stability of the time integration algorithm with an increase in Reynolds number. This instability cannot be easily overcome due to the non-linearities present. The present paper, using the finite element method to integrate the equations in the spacial dimension, incorporates a time-staggered semi-implicit fractional step technique to improve stability at the higher Reynolds numbers. Unlike the upwind or directional differencing schemes normally employed to increase numerical stability, the present scheme does not introduce numerical damping or artificial viscocity, and becomes more stable as the Reynolds number increases. Results for this scheme are compared with various explicit integration schemes for the case of flow around a circular cylinder at Reynolds numbers of 100 to 400. For comparable accuracy the time-staggered semi-implicit fractional step technique was found to be up to 25 times more efficient than the other explicit integration schemes.     

对流扩散方程的高效稳定差分格式

基于二阶修正Dennis格式 ,提出了采用时间相关法求解定常对流扩散方程的一种具有节省内存空间和提高定常解收敛速度的有理式型优化半隐和松驰半隐紧致格式 .本文建立的差分格式具有运算量小、无网格雷诺数限制的优点 ,是无条件稳定和无条件单调的。通过对非线性Burgers方程进行的数值计算结果表明 ,文中构造的有理式型优化半隐和松驰半隐紧致格式适合于非线性问题计算 ,且保持了无条件稳定和无条件单调的特性 ,尤其能使定常解收敛速度加快 ,精度提高 .     

Higher-order semi-implicit Taylor schemes for Itô stochastic differential equations

The paper considers the derivation of families of semi-implicit schemes of weak order N=3.0 (general case) and N=4.0 (additive noise case) for the numerical solution of Itô stochastic differential equations. The degree of implicitness of the schemes depends on the selection of N parameters which vary between 0 and 1 and the families contain as particular cases the 3.0 and 4.0 weak order explicit Taylor schemes. Since the implementation of the multiple integrals that appear in these theoretical schemes is difficult, for the applications they are replaced by simpler random variables, obtaining simplified schemes. In this way, for the multidimensional case with one-dimensional noise, we present an infinite family of semi-implicit simplified schemes of weak order 3.0 and for the multidimensional case with additive one-dimensional noise, we give an infinite family of semi-implicit simplified schemes of weak order 4.0. The mean-square stability of the 3.0 family is analyzed, concluding that, as in the deterministic case, the stability behavior improves when the degree of implicitness grows. Numerical experiments confirming the theoretical results are shown.     

A stabilized semi-implicit Galerkin scheme for Navier-Stokes equations

By introducing a time relaxation term for the time derivative of higher frequency components, we proposed a stabilized semi-implicit Galerkin scheme for evolutionary Navier-Stokes equations in this paper. Analysis shows that such a scheme has weaker stability conditions than that of a classical semi-implicit Galerkin scheme and, when a suitable relaxation parameter σ is chosen, it generates an approximate solution with the same accuracy as the classical one. That means the proposed scheme might use a larger time step to generate a bounded approximate solution. Thus it is more suitable for long time simulations.     

High order accurate semi-implicit WENO schemes for hyperbolic balance laws

In this paper we propose a family of well-balanced semi-implicit numerical schemes for hyperbolic conservation and balance laws. The basic idea of the proposed schemes lies in the combination of the finite volume WENO discretization with Roe’s solver and the strong stability preserving (SSP) time integration methods, which ensure the stability properties of the considered schemes [S. Gottlieb, C.-W. Shu, E. Tadmor, Strong stability-preserving high-order time discretization methods, SIAM Rev. 43 (2001) 89-112]. While standard WENO schemes typically use explicit time integration methods, in this paper we are combining WENO spatial discretization with optimal SSP singly diagonally implicit (SDIRK) methods developed in [L. Ferracina, M.N. Spijker, Strong stability of singly diagonally implicit Runge-Kutta methods, Appl. Numer. Math. 58 (2008) 1675-1686]. In this way the implicit WENO numerical schemes are obtained. In order to reduce the computational effort, the implicit part of the numerical scheme is linearized in time by taking into account the complete WENO reconstruction procedure. With the proposed linearization the new semi-implicit finite volume WENO schemes are designed.A detailed numerical investigation of the proposed numerical schemes is presented in the paper. More precisely, schemes are tested on one-dimensional linear scalar equation and on non-linear conservation law systems. Furthermore, well-balanced semi-implicit WENO schemes for balance laws with geometrical source terms are defined. Such schemes are then applied to the open channel flow equations. We prove that the defined numerical schemes maintain steady state solution of still water. The application of the new schemes to different open channel flow examples is shown.     

The composite Milstein methods for the numerical solution of Ito stochastic differential equations

In this paper, we present the composite Milstein methods for the strong solution of Ito stochastic differential equations. These methods are a combination of semi-implicit and implicit Milstein methods. We give a criterion for choosing either the implicit or the semi-implicit scheme at each step of our numerical solution. The stability and convergence properties are investigated and discussed for the linear test equation. The convergence properties for the nonlinear case are shown numerically to be the same as the linear case. The stability properties of the composite Milstein methods are found to be more superior compared to those of the Milstein, the Euler and even better than the composite Euler method. This superiority in stability makes the methods a better candidate for the solution of stiff SDEs.     

Finite volume schemes for nonhomogeneous scalar conservation laws: error estimate

Summary. In this paper, we study finite volume schemes for the nonhomogeneous scalar conservation law with initial condition . The source term may be either stiff or nonstiff. In both cases, we prove error estimates between the approximate solution given by a finite volume scheme (the scheme is totally explicit in the nonstiff case, semi-implicit in the stiff case) and the entropy solution. The order of these estimates is in space-time -norm (h denotes the size of the mesh). Furthermore, the error estimate does not depend on the stiffness of the source term in the stiff case. Received October 21, 1999 / Published online February 5, 2001     

Three-stage stochastic Runge–Kutta methods for stochastic differential equations

In this paper we discuss three-stage stochastic Runge–Kutta (SRK) methods with strong order 1.0 for a strong solution of Stratonovich stochastic differential equations (SDEs). Higher deterministic order is considered. Two methods, a three-stage explicit (E3) method and a three-stage semi-implicit (SI3) method, are constructed in this paper. The stability properties and numerical results show the effectiveness of these methods in the pathwise approximation of several standard test problems.     

Balanced semi-implicit Roe scheme for open channel flow equations

In this paper, we present an extension of balanced upwind explicit finite volume schemes designed in Vuković and Sopta (SIAM J Sci Comput 24(5):1630–1649, 2003) to semi-implicit ones. Particularly, those schemes are applied to open-channel flows with general geometries and we verify the exact conservation property (C-property). We present the algorithm, the proof of exact C-property and results for several test cases. Also, we test and compare balanced semi-implicit and explicit schemes on standard test cases and on cases involving friction, non-uniform bed slopes and strong channel width variations.      

Global error estimation and control in linearly-implicit parallel two-step peer W-methods

The class of linearly-implicit parallel two-step peer W-methods has been designed recently for efficient numerical solutions of stiff ordinary differential equations. Those schemes allow for parallelism across the method, that is an important feature for implementation on modern computational devices. Most importantly, all stage values of those methods possess the same properties in terms of stability and accuracy of numerical integration. This property results in the fact that no order reduction occurs when they are applied to very stiff problems. In this paper, we develop parallel local and global error estimation schemes that allow the numerical solution to be computed for a user-supplied accuracy requirement in automatic mode. An algorithm of such global error control and other technical particulars are also discussed here. Numerical examples confirm efficiency of the presented error estimation and stepsize control algorithm on a number of test problems with known exact solutions, including nonstiff, stiff, very stiff and large-scale differential equations. A comparison with the well-known stiff solver RODAS is also shown.     

The application of Rosenbrock-Wanner type methods with stepsize control in differential-algebraic equations

Summary Two Rosenbrock-Wanner type methods for the numerical treatment of differential-algebraic equations are presented. Both methods possess a stepsize control and an index-1 monitor. The first method DAE34 is of order (3)4 and uses a full semi-implicit Rosenbrock-Wanner scheme. The second method RKF4DA is derived from the Runge-Kutta-Fehlberg 4(5)-pair, where a semi-implicit Rosenbrock-Wanner method is embedded, in order to solve the nonlinear equations. The performance of both methods is discussed in artificial test problems and in technical applications.     

Stability of Galerkin and Inertial Algorithms with variable time step size

This paper provides Galerkin and Inertial Algorithms for solving a class of nonlinear evolution equations. Spatial discretization can be performed by either spectral or finite element methods; time discretization is done by Euler explicit or Euler semi-implicit difference schemes with variable time step size. Moreover, the boundedness and stability of these algorithms are studied. By comparison, we find that the boundedness and stability of Inertial Algorithm are superior to the ones of Galerkin Algorithm in the case of explicit scheme and the boundedness and stability of two algorithms are same in the case of semi-implicit scheme.     

A semi-implicit tidal model of the North European Continental Shelf

A semi-implicit scheme for the numerical solution of the shallow water equations is proposed. An example of successful application is the simulation of the major tidal constituent (M₂) on the European Shelf. Good agreement is found both qualitatively and quantitatively. The basic outlines of the method are presented and some practical aspects of computation are discussed. Apart from accuracy, compared to explicit or to existing semi-implicit schemes, this approach has another important feature: its economy and simplicity.     

Comparing the contractivity properties of semi-implicit methods

Contractivity is a desirable property of numerical integration methods for stiff systems of ordinary differential equations. In this paper, numerical parameters are used to allow a direct and quantitative comparison of the contractivity properties of various methods for non-linear stiff problems. Results are provided for popular Rosenbrock methods and some more recently developed semi-implicit methods.