A High Order Accurate Bound-Preserving Compact Finite Difference Scheme for Two-Dimensional Incompressible Flow

For solving two-dimensional incompressible flow in the vorticity form by the fourth-order compact finite difference scheme and explicit strong stability preserving temporal discretizations,we show that the simple bound-preserving limiter in Li et al.(SIAM J Numer Anal 56:3308–3345,2018)can enforce the strict bounds of the vorticity,if the velocity field satisfies a discrete divergence free constraint.For reducing oscillations,a modified TVB limiter adapted from Cockburn and Shu(SIAM J Numer Anal 31:607–627,1994)is constructed without affecting the bound-preserving property.This bound-preserving finite difference method can be used for any passive convection equation with a divergence free velocity field.

SOME DISCRETE NONLINEAR INEQUALITIES AND APPLICATIONS TO DIFFERENCE EQUATIONS

In this article, the authors establish some new nonlinear difference inequalities in two independent variables, which generalize some existing results and can be used as handy tools in the study of qualitative as well as quantitative properties of solutions of certain classes of difference equations.

A Stable FE-FD Method for Anisotropic Parabolic PDEs with Moving Interfaces

In this paper,a new finite element and finite difference(FE-FD)method has been developed for anisotropic parabolic interface problems with a known moving interface using Cartesian meshes.In the spatial discretization,the standard P,FE discretization is applied so that the part of the coefficient matrix is symmetric positive definite,while near the interface,the maximum principle preserving immersed interface discretization is applied.In the time discretization,a modified Crank-Nicolson discretization is employed so that the hybrid FE-FD is stable and second order accurate.Correction terms are needed when the interface crosses grid lines.The moving interface is represented by the zero level set of a Lipschitz continuous function.Numerical experiments presented in this paper confirm second orderconvergence.

DIFFERENCE DISCRETE SYSTEM OF EULER-BEAM WITH ARBITRARY SUPPORTS AND SIGN-OSCILLATORY PROPERTY OF STIFFNESS MATRICES

The difference discrete system of Euler-beam with arbitrary supports was constructed by using the two order central difference formulas. This system is equivalent to the spring-mass-rigidrod model. By using the theory of oscillatory matrix, the signoscillatory property of stiffness matrices of this system was proved, and the necessary and sufficient condition for the system to be positive was obtained completely.

ELEMENT FUNCTIONS OF DISCRETE OPERATOR DIFFERENCE METHOD

The discrete scheme called discrete operator difference for differential equations was given. Several difference elements for plate bending problems and plane problems were given. By investigating these elements, the ability of the discrete forms expressing to the element functions was talked about. In discrete operator difference method, the displacements of the elements can be reproduced exactly in the discrete forms whether the displacements are conforming or not. According to this point, discrete operator difference method is a method with good performance.

NUMERICAL SOLUTION OF A SINGULARLY PERTURBED ELLIPTIC-HYPERBOLIC PARTIAL DIFFERENTIAL EQUATION ON A NONUNIFORM DISCRETIZATION MESH

In this paper, we consider the upwind difference scheme for singular perturbation problem (1.1). On a special discretization mesh, it is proved that the solution of the upwind difference scheme is first order convergent, uniformly in the small parameter e , to the solution of problem (1.1). Numerical results are finally provided.

New exact solutions to some difference differential equations

In this paper, we use our method to solve the extended Lotka-Volterra equation and discrete KdV equation. With the help of Maple, we obtain a number of exact solutions to the two equations including soliton solutions presented by hyperbolic functions of sinh and cosh, periodic solutions presented by trigonometric functions of sin and cos, and rational solutions. This method can be used to solve some other nonlinear difference-differential equations.

Equivalent low-order angular flux nonlinear finite difference equation of MOC transport calculation

The key issue in accelerating method of characteristics(MOC)transport calculations is in obtaining a completely equivalent low-order neutron transport or diffusion equation.Herein,an equivalent low-order angular flux nonlinear finite difference equation is proposed for MOC transport calculations.This method comprises three essential features:(1)the even parity discrete ordinates method is used to build a low-order angular flux nonlinear finite difference equation,and different boundary condition treatments are proposed;(2)two new defined factors,i.e.,the even parity discontinuity factor and odd parity discontinuity factor,are strictly defined to achieve equivalence between the low-order angular flux nonlinear finite difference method and MOC transport calculation;(3)the energy group and angle are decoupled to construct a symmetric linear system that is much easier to solve.The equivalence of this low-order angular flux nonlinear finite difference equation is analyzed for two-dimensional(2D)pin,2D assembly,and 2D C5G7 benchmark problems.Numerical results demonstrate that a low-order angular flux nonlinear finite difference equation that is completely equivalent to the pin-resolved transport equation is established.

An Algebraic Method for Constructing Exact Solutions to Difference-Differential Equations

In this paper, we present a method to solve difference differential equation（s）. As an example, we apply this method to discrete KdV equation and Ablowitz-Ladik lattice equation. As a result, many exact solutions are obtained with the help of Maple including soliton solutions presented by hyperbolic functions sinh and cosh, periodic solutions presented by sin and cos and rational solutions. This method can also be used to other nonlinear difference-differential equation（s）.

Traveling Waves for 2-1 Dimension Lattice Difference Equations

A definition is introduced about traveling waves of 2-1 dimension lattice difference equations. Discrete heat equation is introduced and a discussion is given for the existence of traveling waves. The theory of traveling waves is extended on 2-1 dimension lattice difference equations. As an application, an example is presented to illustrate the main results.

Finite difference streamline diffusion method using nonconforming space for incompressible time-dependent Navier-Stokes equations

This paper proposes a new nonconforming finite difference streamline diffusion method to solve incompressible time-dependent Navier-Stokes equations with a high Reynolds number. The backwards difference in time and the Crouzeix-Raviart （CR） element combined with the P0 element in space are used. The result shows that this scheme has good stabilities and error estimates independent of the viscosity coefficient.

The dynamics of a memristor-based Rulkov neuron with fractional-order difference

The exploration of the memristor model in the discrete domain is a fascinating hotspot.The electromagnetic induction on neurons has also begun to be simulated by some discrete memristors.However,most of the current investigations are based on the integer-order discrete memristor,and there are relatively few studies on the form of fractional order.In this paper,a new fractional-order discrete memristor model with prominent nonlinearity is constructed based on the Caputo fractional-order difference operator.Furthermore,the dynamical behaviors of the Rulkov neuron under electromagnetic radiation are simulated by introducing the proposed discrete memristor.The integer-order and fractional-order peculiarities of the system are analyzed through the bifurcation graph,the Lyapunov exponential spectrum,and the iterative graph.The results demonstrate that the fractional-order system has more abundant dynamics than the integer one,such as hyper-chaos,multi-stable and transient chaos.In addition,the complexity of the system in the fractional form is evaluated by the means of the spectral entropy complexity algorithm and consequences show that it is affected by the order of the fractional system.The feature of fractional difference lays the foundation for further research and application of the discrete memristor and the neuron map in the future.

Finite Difference Method of Modelling Groundwater Flow

In this study, finite difference method is used to solve the equations that govern groundwater flow to obtain flow rates, flow direction and hydraulic heads through an aquifer. The aim therefore is to discuss the principles of Finite Difference Method and its applications in groundwater modelling. To achieve this, a rectangular grid is overlain an aquifer in order to obtain an exact solution. Initial and boundary conditions are then determined. By discretizing the system into grids and cells that are small compared to the entire aquifer, exact solutions are obtained. A flow chart of the computational algorithm for particle tracking is also developed. Results show that under a steady-state flow with no recharge, pathlines coincide with streamlines. It is also found that the accuracy of the numerical solution by Finite Difference Method is largely dependent on initial particle distribution and number of particles assigned to a cell. It is therefore concluded that Finite Difference Method can be used to predict the future direction of flow and particle location within a simulation domain.

A New Fourth Order Difference Approximation for the Solution of Three-dimensional Non-linear Biharmonic Equations Using Coupled Approach

This paper deals with a new higher order compact difference scheme, which is, O(h4) using coupled approach on the 19-point 3D stencil for the solution of three dimensional nonlinear biharmonic equations. At each internal grid point, the solution u(x,y,z) and its Laplacian Δ4u are obtained. The resulting stencil algo-rithm is presented and hence this new algorithm can be easily incorporated to solve many problems. The present discretization allows us to use the Dirichlet boundary conditions only and there is no need to discretize the derivative boundary conditions near the boundary. We also show that special treatment is required to handle the boundary conditions. Convergence analysis for a model problem is briefly discussed. The method is tested on three problems and compares very favourably with the corresponding second order approximation which we also discuss using coupled approach.

Numerical Solution of Advection Diffusion Equation Using Semi-Discretization Scheme

Numerical diffusion and oscillatory behavior characteristics are averted applying numerical solutions of advection-diffusion equation are themselves immensely sophisticated. In this paper, two numerical methods have been used to solve the advection diffusion equation. We use an explicit finite difference scheme for the advection diffusion equation and semi-discretization on the spatial variable for advection-diffusion equation yields a system of ordinary differential equations solved by Euler’s method. Numerical assessment has been executed with specified initial and boundary conditions, for which the exact solution is known. We compare the solutions of the advection diffusion equation as well as error analysis for both schemes.

An Implicit-Explicit Computational Method Based on Time Semi-Discretization for Pricing Financial Derivatives with Jumps

This paper considers pricing European options under the well-known of SVJ model of Bates and related computational methods. According to the no-arbitrage principle, we first derive a partial differential equation that the value of any European contingent claim should satisfy, where the asset price obeys the SVJ model. This equation is numerically solved by using the implicit- explicit backward difference method and time semi-discretization. In order to explain the validity of our method, the stability of time semi-discretization scheme is also proved. Finally, we use a simulation example to illustrate the efficiency of the method.

基于DEM-FDM耦合的过渡段膨胀诱发钢轨上拱研究

研究目的:为分析涵洞过渡段地基膨胀引起的钢轨上拱响应,基于现场测试、室内膨胀试验数据,开展DEM-FDM耦合数值模拟,分析某涵洞附近路基土在膨胀范围为16 m,膨胀中心距离涵洞中心分别为0 m、5 m、10 m这三种工况下,不同膨胀率时基床填料的运动规律及钢轨的上拱响应。研究结论:(1)涵洞对于钢轨上拱位移的传递存在阻断作用,但会增大钢轨上拱的峰值,原位膨胀率下工况二的钢轨上拱峰值达到46 mm,当路基膨胀率为0.3%时钢轨上拱位移量达到无砟轨道钢轨可调节临界值(4mm);(2)过渡段钢轨上拱处同时产生轴向应力集中,其中原位膨胀率下工况二轴向应力峰值达到14.4 MPa;(3)对于膨胀区域位于涵洞下方的工况,钢轨轴向应力呈现出来的分布规律为钢轨上拱拱顶处为主拉应力状态,拱脚处为主压应力状态,因此一共包括三个压应力峰值点以及两个拉应力峰值点;(4)本文研究可为高铁涵洞过渡段路基膨胀病害解决方案的确定提供理论依据。

基于包装材料数字图像的RGB颜色空间色差评价方法研究

目的解决包装材料复杂颜色区域色差评价,目测检验方法存在判定标准复杂、结果一致性差等问题。方法提出一种基于包装材料高清数字图像的RGB颜色空间色差评价方法。通过搭建高清数字图像采集装置,获取包装材料标准样(上限、中限、下限)的全幅数字图像。通过采用RGB颜色空间降维的方法,对标准样测试区域的颜色空间进行离散边界拟合,构建标准样品测试区域的三限样RGB色差模型,与测试样对应测试区域的颜色空间进行对比;根据颜色空间模型的一致性判定测试样色差是否合格。结果对比多种离散边界拟合模型,采用滚球法边缘拟合模型的包含率和多出率分别达到100%、78.9%,模型拟合效果最佳。通过对包装材料色差缺陷样验证测试,采用上述色差评价方法准确率达到100%。结论本文提出的滚球法颜色空间边缘拟合模型和色差评价方法,可实现对包装材料复杂颜色区域的色差合格性判定,有利于提高企业生产的智能化检测水平。

一维平流方程迎风格式的最优时空步长

在考虑舍入误差影响的情况下,研究一维平流方程迎风格式最优时空步长的选取.首先,分析每一时间层产生的离散误差和舍入误差,以及2种误差向高时间层传播的累积,得到数值解总误差的理论上界;然后推导出最优时间步长和最优空间步长的理论公式,进而得到2种不同机器精度下最优时间步长之比满足的一个仅与机器精度有关的普适关系;最后通过数值算例验证了结论的可靠性.

基于L_(2)灵敏度最小化的数字滤波器状态空间实现稀疏化方法

为解决传统数字滤波器在有限精度实现时因有限字长(Finite Word Length,FWL)效应导致滤波器性能下降的问题,提出一种L_(2)灵敏度最小化的数字滤波器状态空间实现稀疏化方法.推导前向差分算子数字滤波器结构传输函数及其等效状态空间实现,根据可控及可观格莱姆矩阵得到基于相似变换矩阵的L_(2)灵敏度表达式,并进行稀疏化校准,将L_(2)灵敏度最小化问题转换为凸函数求最值问题,求导得到L_(2)灵敏度最小化表达式,代回即得前向差分算子数字滤波器的稀疏化状态空间实现.仿真结果表明,所提方法设计的数字滤波器具有更好的抗FWL效应.