Fast algorithm for viscous Cahn-Hilliard equation

The main purpose of this paper is to solve the viscous Cahn-Hilliard equation via a fast algorithm based on the two time-mesh (TT-M) finite element (FE) method to ease the problem caused by strong nonlinearities. The TT-M FE algorithm includes the following main computing steps. First, a nonlinear FE method is applied on a coarse time-mesh τ_c. Here, the FE method is used for spatial discretization and the implicit second-order θ scheme (containing both implicit Crank-Nicolson and second-order backward difference) is used for temporal discretization. Second, based on the chosen initial iterative value, a linearized FE system on time fine mesh is solved, where some useful coarse numerical solutions are found by Lagrange’s interpolation formula. The analysis for both stability and a priori error estimates is made in detail. Numerical examples are given to demonstrate the validity of the proposed algorithm. Our algorithm is compared with the traditional Galerkin FE method and it is evident that our fast algorithm can save computational time.     

An alternative to the Bathe algorithm

This paper presents a new composite sub-steps algorithm for solving reliable numerical responses in structural dynamics. The newly developed algorithm is a two sub-steps, second-order accurate and unconditionally stable implicit algorithm with the same numerical properties as the Bathe algorithm. The detailed analysis of the stability and numerical accuracy is presented for the new algorithm, which shows that its numerical characteristics are identical to those of the Bathe algorithm. Hence, the new sub-steps scheme could be considered as an alternative to the Bathe algorithm. Meanwhile, the new algorithm possesses the following properties: (a) it produces the same accurate solutions as the Bathe algorithm for solving linear and nonlinear problems; (b) it does not involve any artificial parameters and additional variables, such as the Lagrange multipliers; (c) The identical effective stiffness matrices can be obtained inside two sub-steps; (d) it is a self-starting algorithm. Some numerical experiments are given to show the superiority of the new algorithm and the Bathe algorithm over the dissipative CH-α algorithm and the non-dissipative trapezoidal rule.     

Implicit algorithms for solving the Cauchy problem for the equations of the dynamics of mechanical systems

It is shown that in the numerical solution of the Cauchy problem for systems of second-order ordinary differential equations, when solved for the highest-order derivative, it is possible to construct simple and economical implicit computational algorithms for step-by-step integration without using laborious iterative procedures based on processes of the Newton-Raphson iterative type. The initial problem must first be transformed to a new argument — the length of its integral curve. Such a transformation is carried out using an equation relating the initial parameter of the problem to the length of the integral curve. The linear acceleration method is used as an example to demonstrate the procedure of constructing an implicit algorithm using simple iterations for the numerical solution of the transformed Cauchy problem. Propositions concerning the computational properties of the iterative process are formulated and proved. Explicit estimates are given for an integration stepsize that guarantees the convergence of the simple iterations. The efficacy of the proposed procedure is demonstrated by the numerical solution of three problems. A comparative analysis is carried out of the numerical solutions obtained with and without parametrization of the initial problems in these three settings. As a qualitative test the problem of the celestial mechanics of the “Pleiades” is considered. The second example is devoted to modelling the non-linear dynamics of an elastic flexible rod fixed at one end as a cantilever and coiled in its initial (static) state into a ring by a bending moment. The third example demonstrates the numerical solution of the problem of the “unfolding” of a mechanical system consisting of three flexible rods with given control input.     

An Implicit Solver for the Time-Dependent Kohn-Sham Equation

The implicit numerical methods have the advantages on preserving the physical properties of the quantum system when solving the time-dependent Kohn-Sham equation. However, the efficiency issue prevents the practical applications of those implicit methods. In this paper, an implicit solver based on a class of Runge-Kutta methods and the finite element method is proposed for the time-dependent Kohn-Sham equation. The efficiency issue is partially resolved by three approaches, i.e., an $h$-adaptive mesh method is proposed to effectively restrain the size of the discretized problem, a complex-valued algebraic multigrid solver is developed for efficiently solving the derived linear system from the implicit discretization, as well as the OpenMP based parallelization of the algorithm. The numerical convergence, the ability on preserving the physical properties, and the efficiency of the proposed numerical method are demonstrated by a number of numerical experiments.     

求解非线性时滞脉冲双曲微分方程的隐式差分格式

构造了一类求解非线性时滞脉冲双曲型偏微分方程的隐式差分格式.在一定条件下,获得了该差分格式的唯一可解性、收敛性和无条件稳定性,且空间和时间均二阶精度.最后,数值实验表明了所得格式的精度和有效性.     

NONLINEAR GALERKIN METHOD AND CRANK-NICOLSONMETHOD FOR VISCOUS INCOMPRESSIBLE FLOW

1.Introduction'NonlinearGalerkinmethodisnumericalmethodfordissipativeevolutionpartialdifferentialequationswherethespatialdiscretizationreliesonanonlinearmanifoldinsteadofalinearspaceasintheclassicalGalerkinmethod.Morepreciselygoneconsidersafinitedimension…     

Numerical discretization-based kernel type estimation methods for ordinary differential equation models

We consider the problem of parameter estimation in both linear and nonlinear ordinary differential equation(ODE) models. Nonlinear ODE models are widely used in applications. But their analytic solutions are usually not available. Thus regular methods usually depend on repetitive use of numerical solutions which bring huge computational cost. We proposed a new two-stage approach which includes a smoothing method(kernel smoothing or local polynomial fitting) in the first stage, and a numerical discretization method(Eulers discretization method, the trapezoidal discretization method,or the Runge–Kutta discretization method) in the second stage. Through numerical simulations, we find the proposed method gains a proper balance between estimation accuracy and computational cost.Asymptotic properties are also presented, which show the consistency and asymptotic normality of estimators under some mild conditions. The proposed method is compared to existing methods in term of accuracy and computational cost. The simulation results show that the estimators with local linear smoothing in the first stage and trapezoidal discretization in the second stage have the lowest average relative errors. We apply the proposed method to HIV dynamics data to illustrate the practicability of the estimator.     

Frequency-dependent vibration analysis of symmetric cross-ply laminated plate of Levy-type by spectral element and finite strip procedures

This research describes spectral finite element formulation for vibration analysis of rectangular symmetric cross-ply laminated composite plates of Levy-type based on classical lamination plate theory (CLPT). Formulation based on SFEM includes partial differential equations of motion, spectral displacement field, dynamic shape functions, and spectral element stiffness matrix (SESM). In this paper, vibration analysis of composite plate is investigated in two sections: free vibrations and forced vibrations. In free vibrations, natural frequencies are calculated for different Young’s moduli ratios and boundary conditions. In forced vibrations, plate vibrations are investigated under high-frequency concentrated impulsive loads. The resulting responses due to spectral element formulation are compared with those of (time-domain) finite element and analytical formulations, whenever available. The results demonstrate the superiority of SFEM with respect to FEM, in reducing computational burden, simultaneously increasing numerical accuracy, specifically for excitations of high-frequency content. By reducing the time duration of impulsive loads, and consequently increasing the modal contribution of higher modes in the transient response of plate, the accuracy of FEM responses decreases substantially accompanied with a high volume of computations, while the accuracy of the SFEM response results is very high and simultaneously, with a low computational burden. Practically, SFEM follows very closely exact analytical solutions.     

Multiscale computational method for thermoelastic problems of composite materials with orthogonal periodic configurations

This study develops a novel multiscale computational method for thermoelastic problems of composite materials with orthogonal periodic configurations. Firstly, the multiscale asymptotic analysis for these multiscale problems is given successfully, and the formal second-order two-scale approximate solutions for these multiscale problems are constructed based on the above-mentioned analysis. Then, the error estimates for the second-order two-scale (SOTS) solutions are obtained. Furthermore, the corresponding SOTS numerical algorithm based on finite element method (FEM) is brought forward in details. Finally, some numerical examples are presented to verify the feasibility and effectiveness of our multiscale computational method. Moreover, our multiscale computational method can accurately capture the local thermoelastic responses in composite block structure, plate, cylindrical and doubly-curved shallow shells.     

Numerical algorithm for solving diffusion-type equations on the basis of multigrid methods

A new numerical algorithm based on multigrid methods is proposed for solving equations of the parabolic type. Theoretical error estimates are obtained for the algorithm as applied to a two-dimensional initial-boundary value model problem for the heat equation. The good accuracy of the algorithm is demonstrated using model problems including ones with discontinuous coefficients. As applied to initial-boundary value problems for diffusion equations, the algorithm yields considerable savings in computational work compared to implicit schemes on fine grids or explicit schemes with a small time step on fine grids. A parallelization scheme is given for the algorithm.     

广义sine-Gordon方程高精度隐式紧致差分方法

本文研究一类二维非线性的广义sine-Gordon(简称SG)方程的有限差分格式.首先构造三层时间的紧致交替方向隐式差分格式,并用能量分析法证明格式具有二阶时间精度和四阶空间精度.然后应用改进的Richardson外推算法将时间精度提高到四阶.最后,数值算例证实改进后的算法在空间和时间上均达到四阶精度.     

Error Analysis of SAV Finite Element Method to Phase Field Crystal Model

In this paper, we construct and analyze an energy stable scheme by combining the latest developed scalar auxiliary variable (SAV) approach and linear finite element method (FEM) for phase field crystal (PFC) model, and show rigorously that the scheme is first-order in time and second-order in space for the $L^2$ and $H^{-1}$ gradient flow equations. To reduce efficiently computational cost and capture accurately the phase interface, we give a simple adaptive strategy, equipped with a posteriori gradient estimator, i.e., $L^2$ norm of the recovered gradient. Extensive numerical experiments are presented to verify our theoretical results and to demonstrate the effectiveness and accuracy of our proposed method.     

An unconditionally stable spatial sixth-order CCD-ADI method for the two-dimensional linear telegraph equation

The telegraph equation is one of the important models in many physics and engineering. In this work, we discuss the high-order compact finite difference method for solving the two-dimensional second-order linear hyperbolic equation. By using a combined compact finite difference method for the spatial discretization, a high-order alternating direction implicit method (ADI) is proposed. The method is O(τ² + h⁶) accurate, where τ, h are the temporal step size and spatial size, respectively. Von Neumann linear stability analysis shows that the method is unconditionally stable. Finally, numerical examples are used to illustrate the high accuracy of the new difference scheme.     

Fast Evaluation of Linear Combinations of Caputo Fractional Derivatives and Its Applications to Multi-Term Time-Fractional Sub-Diffusion Equations

In the present work, linear combinations of Caputo fractional derivatives are fast evaluated based on the efficient sum-of-exponentials (SOE) approximation for kernels in Caputo fractional derivatives with an absolute error $\epsilon,$ which is a further work of the existing results in [13] (Commun. Comput. Phys., 21 (2017), pp. 650-678) and [16] (Commun. Comput. Phys., 22 (2017), pp. 1028-1048). Both the storage needs and computational amount are significantly reduced compared with the direct algorithm. Applications of the proposed fast algorithm are illustrated by solving a second-order multi-term time-fractional sub-diffusion problem. The unconditional stability and convergence of the fast difference scheme are proved. The CPU time is largely reduced while the accuracy is kept, especially for the cases of large temporal level, which is displayed by numerical experiments.     

第二类Abel型积分方程的高精度组合算法与后验估计

本文提出解第二类Abel型积分方程的高精度组合算法与后验估计,理论与算例皆表明本文方法省计算时间、省存贮,精度高且近似解有后验估计。     

二阶锥线性互补问题的广义模系矩阵分裂迭代算法

通过将二阶锥线性互补问题转化为等价的不动点方程,介绍了一种广义模系矩阵分裂迭代算法,并研究了该算法的收敛性.进一步,数值结果表明广义模系矩阵分裂迭代算法能够有效地求解二阶锥线性互补问题.     

求解一维非定常对流扩散反应方程的高精度紧致差分格式

提出了数值求解一维非定常对流扩散反应方程的一种高精度紧致隐式差分格式,其截断误差为O(τ~4+τ~2h~2+h~4),即格式整体具有四阶精度.差分方程在每一时间层上只用到了三个网格节点,所形成的代数方程组为三对角型,可采用追赶法进行求解,最后通过数值算例验证了格式的精确性和可靠性.     

Novel numerical techniques for the finite moment log stable computational model for European call option

Option pricing models are often used to describe the dynamic characteristics of prices in financial markets. Unlike the classical Black–Scholes (BS) model, the finite moment log stable (FMLS) model can explain large movements of prices during small time steps. In the FMLS, the second-order spatial derivative of the BS model is replaced by a fractional operator of order α which generates an α-stable Lévy process. In this paper, we consider the finite difference method to approximate the FMLS model. We present two numerical schemes for this approximation: the implicit numerical scheme and the Crank–Nicolson scheme. We carry out convergence and stability analyses for the proposed schemes. Since the fractional operator routinely generates dense matrices which often require high computational cost and storage memory, we explore three methods for solving the approximation schemes: the Gaussian elimination method, the bi-conjugate gradient stabilized method (Bi-CGSTAB) and the fast Bi-CGSTAB (FBi-CGSTAB) in order to compare the cost of calculations. Finally, two numerical examples with exact solutions are presented where we also use extrapolation techniques to achieve higher-order convergence. The results suggest that the proposed schemes are unconditionally stable and convergent, and the FMLS model is useful for pricing options.     

Semi-Lagrangian semi-implicit time-splitting scheme for a regional model of the atmosphere

Semi-Lagrangian semi-implicit (SLSI) method is currently one of the most efficient approaches for numerical solution of the atmosphere dynamics equations. In this research we apply splitting techniques in the context of a two-time-level SLSI scheme in order to simplify the treatment of the slow physical modes and optimize the solution of the elliptic equations related to implicit part of the scheme. The performed numerical experiments show the accuracy and computational efficiency of the scheme.