基于人工鱼群算法求任意函数的数值积分

将不等距离分割方法与人工鱼群算法相结合,提出一种基于人工鱼群算法求任意函数数值积分的方法,该方法除能计算通常意义下任意函数的定积分外,还能计算奇异函数积分、振荡函数积分以及原函数不易求得的被积函数的积分.最后给出几个数值积分算例,并与传统数值积分方法作了比较,仿真结果分析表明,该算法十分有效,能够快速有效地获得任意函数的数值积分值.     

一类非线性延迟微分方程数值解的振动性分析

本文考虑一类非线性延迟微分方程-带有单调造血率的造血模型数值解的振动性.通过研究特征方程根的情况得到数值解振动的条件并且讨论了非振动的数值解的一些性质.为了更有力的说明我们的结果,最后给出了相应的算例.     

A finite volume element formulation and error analysis for the non-stationary conduction–convection problem

The non-stationary conduction–convection problem including the velocity vector field and the pressure field as well as the temperature field is studied with a finite volume element (FVE) method. A fully discrete FVE formulation and the error estimates between the fully discrete FVE solutions and the accuracy solution are provided. It is shown by numerical examples that the results of numerical computation are consistent with theoretical conclusions. Moreover, it is shown that the FVE method is feasible and efficient for finding the numerical solutions of the non-stationary conduction–convection problem and is one of the most effective numerical methods by comparing the results of the numerical simulations of the FVE formulation with those of the numerical simulations of the finite element method and the finite difference scheme for the non-stationary conduction–convection problem.     

Error analysis of the numerical solution of split differential equations

The operator splitting method is a widely used approach for solving partial differential equations describing physical processes. Its application usually requires the use of certain numerical methods in order to solve the different split sub-problems. The error analysis of such a numerical approach is a complex task. In the present paper we show that an interaction error appears in the numerical solution when an operator splitting procedure is applied together with a lower-order numerical method. The effect of the interaction error is investigated by an analytical study and by numerical experiments made for a test problem.     

A non-standard finite difference method to solve a model of HIV–Malaria co-infection

In this paper, we design and analyse a non-standard finite difference numerical scheme for the numerical solution of the HIV–malaria co-infection model with a distributed delay representing the incubation period of malaria parasite in the mosquitoes vector. To come up with the efficient numerical method for the full co-infection model, we study a number of qualitative properties of sub-models and then use the information while designing the numerical methods for these sub-models. One of the salient features of these methods is that they preserve positivity of the solution which is very essential while studying epidemiological models. We also present numerical simulations to confirm the theoretical findings.     

Numerical stability of higher-order derivative methods for the pantograph equation

Dealing with numerical stability of higher-order derivative methods with variable stepsize is the purpose of this paper for pantograph equations. A new way to compute this kind of equation is provided, and a sufficient condition for the numerical stability of high order derivative forms is given. Some numerical examples are presented to confirm our theoretical analysis.     

A numerical approach for solving an extended Fisher-Kolomogrov-Petrovskii-Piskunov equation

In the present paper a numerical method, based on finite differences and spline collocation, is presented for the numerical solution of a generalized Fisher integro-differential equation. A composite weighted trapezoidal rule is manipulated to handle the numerical integrations which results in a closed-form difference scheme. A number of test examples are solved to assess the accuracy of the method. The numerical solutions obtained, indicate that the approach is reliable and yields results compatible with the exact solutions and consistent with other existing numerical methods. Convergence and stability of the scheme have also been discussed.     

Numerical methods for scattering problems from multi-layers with different periodicities

In this paper, we consider a numerical method to solve scattering problems with multi-periodic layers with different periodicities. The main tool applied in this paper is the Bloch transform. With this method, the problem is written into an equivalent coupled family of quasi-periodic problems. As the Bloch transform is only defined for one fixed period, the inhomogeneous layer with another period is simply treated as a non-periodic one. First, we approximate the refractive index by a periodic one where its period is an integer multiple of the fixed period, and it is decomposed by finite number of quasi-periodic functions. Then the coupled system is reduced into a simplified formulation. A convergent finite element method is proposed for the numerical solution, and the numerical method has been applied to several numerical experiments. At the end of this paper, relative errors of the numerical solutions will be shown to illustrate the convergence of the numerical algorithm.     

用线方法求解带混合边值条件的抛物方程

研制了分别用显式Euler法、隐式Euler法、Crank-Nicolson格式(梯形方法)求解带第一、第二及混合边值条件的抛物问题的应用软件,通过求解若干抛物问题对该软件作了测试,获得了预期的数值结果,讨论了时间和空间步长的变化对格式计算结果的影响,得到了三种方法的稳定性、收敛精度和计算量.     

A problem on heat conduction in an insulated wire solved by numerical inverse Laplace transform

A problem of transient heat conduction in an insulated wire is solved by use of Laplace transform and numerical inversion. The problem is solved for the radiation boundary condition and also for the boundary condition of no heat flux through the outer surface of the insulation. The results are presented both numerically with four significant figures and graphically. Asymptotic expansions are derived for small and large values of the time variable. The numerical inversion of the Laplace transform is checked by comparison with the asymptotic expansions and with the numerical results obtained by a numerical inversion formula utilizing one more abscissa than the previous one.     

The numerical method of successive interpolations for Fredholm functional integral equations

A new numerical method for Fredholm functional integral equations is proposed. The method combines the fixed point technique with numerical integration and cubic spline interpolation. The convergence and the numerical stability of the method are proved and tested on some numerical examples.     

A Fourier method and an extrapolation technique for Stokes’ first problem for a heated generalized second grade fluid with fractional derivative

In this article, we consider Stokes’ first problem for a heated generalized second grade fluid with fractional derivative (SFP-HGSGF). Implicit and explicit numerical approximation schemes for the SFP-HGSGF are presented. The stability and convergence of the numerical schemes are discussed using a Fourier method. In addition, the solvability of the implicit numerical approximation scheme is also analyzed. A Richardson extrapolation technique for improving the order of convergence of the implicit scheme is proposed. Finally, a numerical test is given. The numerical results demonstrate the good performance of our theoretical analysis.     

The Numerical Simulation of Space-Time Variable Fractional Order Diffusion Equation

Many physical processes appear to exhibit fractional order behavior that may vary with time or space. The continuum of order in the fractional calculus allows the order of the fractional operator to be considered as a variable. Numerical methods and analysis of stability and convergence of numerical scheme for the variable fractional order partial differential equations are quite limited and difficult to derive. This motivates us to develop efficient numerical methods as well as stability and convergence of the implicit numerical methods for the space-time variable fractional order diffusion equation on a finite domain. It is worth mentioning that here we use the Coimbra-definition variable time fractional derivative which is more efficient from the numerical standpoint and is preferable for modeling dynamical systems. An implicit Euler approximation is proposed and then the stability and convergence of the numerical scheme are investigated. Finally, numerical examples are provided to show that the implicit Euler approximation is computationally efficient.     

Multiresolution algorithms for the numerical solution of hyperbolic conservation laws

Given any scheme in conservation form and an appropriate uniform grid for the numerical solution of the initial value problem for one-dimensional hyperbolic conservation laws we describe a multiresolution algorithm that approximates this numerical solution to a prescribed tolerance in an efficient manner. To do so we consider the grid-averages of the numerical solution for a hierarchy of nested diadic grids in which the given grid is the finest, and introduce an equivalent multiresolution representation. The multiresolution representation of the numerical solution consists of its grid-averages for the coarsest grid and the set of errors in predicting the grid-averages of each level of resolution in this hierarchy from those of the next coarser one. Once the numerical solution is resolved to our satisfaction in a certain locality of some grid, then the prediction errors there are small for this particular grid and all finer ones; this enables us to compress data by setting to zero small components of the representation which fall below a prescribed tolerance. Therefore instead of computing the time-evolution of the numerical solution on the given grid we compute the time-evolution of its compressed multiresolution representation. Algorithmically this amounts to computing the numerical fluxes of the given scheme at the points of the given grid by a hierarchical algorithm which starts with the computation of these numerical fluxes at the points of the coarsest grid and then proceeds through diadic refinements to the given grid. At each step of refinement we add the values of the numerical flux at the center of the coarser cells. The information in the multiresolution representation of the numerical solution is used to determine whether the solution is locally well-resolved. When this is the case we replace the costly exact value of the numerical flux with an accurate enough approximate value which is obtained by an inexpensive interpolation from the coarser grid. The computational efficiency of this multiresolution algorithm is proportional to the rate of data compression (for a prescribed level of tolerance) that can be achieved for the numerical solution of the given scheme.     

Numerical simulation for the variable-order Galilei invariant advection diffusion equation with a nonlinear source term

In this paper, we consider the variable-order Galilei advection diffusion equation with a nonlinear source term. A numerical scheme with first order temporal accuracy and second order spatial accuracy is developed to simulate the equation. The stability and convergence of the numerical scheme are analyzed. Besides, another numerical scheme for improving temporal accuracy is also developed. Finally, some numerical examples are given and the results demonstrate the effectiveness of theoretical analysis.     

High-resolution schemes for bubbling flow computations

This paper is concerned with the effects of numerical schemes on the simulation of dense gas-particle two-phase flows. The first-order upwind difference, the central difference, the second-order upwind difference, the central difference plus artificial dissipation, the deferred correction, the quadratic upstream interpolation (for convective kinematics) and the monotone upstream-centered schemes for conservation laws with different flux limiters were all accomplished to simulate the multi-phase flows. It was found that numerical schemes may significantly affect the solution accuracy and numerical convergence. The monotone upstream-centered schemes for conservation laws are the best choice of all, because they can effectively suppress the non-physical oscillations with the introduction of adaptive numerical dissipation into numerical solutions.     

带漂移的单侧正规化回火分数阶扩散方程的Crank-Nicolson拟紧格式

基于已有的针对单侧正规化回火分数阶扩散方程的三阶拟紧算法,将该算法的思想应用于带漂移的单侧正规化回火分数阶扩散方程的数值模拟,并结合Crank-Nicolson方法导出数值格式.证明了数值格式的稳定性与收敛性,且数值格式的时间收敛阶和空间收敛阶分别是二阶和三阶.通过数值试验验证了数值格式的有效性和理论结果.     

Numerical solution of an inverse initial boundary-value problem for the full time-dependent Maxwell's equations in the presence of imperfections of small volume

We consider the numerical solution, in a three-dimensional bounded domain, of the inverse problem for identifying the location of small electromagnetic imperfections in a medium with homogeneous background. Our numerical algorithm is based on the coupling of a discontinuous Galerkin method for the time-dependent Maxwell's equations, on the exact controllability method and on a Fourier inversion. Several numerical results are given with one and two imperfections and the robustness and accuracy of the numerical method used for the dynamic detection problem are shown.     

A meshless method for Burgers’ equation using MQ-RBF and high-order temporal approximation

In this paper, a new numerical method is proposed to solve one-dimensional Burgers’ equation using multiquadric (MQ) radial basis function (RBF) for spatial approximation and a second-order compact finite difference scheme for temporal approximation. The numerical results obtained by this way for different Reynolds number have been compared with the existing numerical schemes to show the accuracy and efficiency of the approach. To show the superiority of this meshless method, numerical experiments with non-uniform MQ interpolation node distribution are also performed.     

Numerical pseudodifferential operator and Fourier regularization

The concept of numerical pseudodifferential operator, which is an extension of numerical differentiation, is suggested. Numerical pseudodifferential operator just is calculating the value of the pseudodifferential operator with unbounded symbol. Many ill-posed problems can lead to numerical pseudodifferential operators. Fourier regularization is a very simple and effective method for recovering the stability of numerical pseudodifferential operators. A systematically theoretical analysis and some concrete examples are provided.