关于定积分近似计算中矩形法的误差估计

利用拉格朗日中值定理给出复合矩形法的误差估计,并指出复合矩形法只有在等分的次数很大的时候才能比较精确的估计所求定积分的值.最后,给出复合矩形法、复合梯形法及复合抛物线法的误差估计实例.     

五维热传导方程的一族两层显格式

构造了五维热传导方程的一族两层显格式,证明了当截断误差阶为O（τ＋h2）时,其稳定性条件为网比r=hτ2≤21,优于同类的其它显格式,当截断误差阶为O（τ2＋h2）时,可以得到一个简洁而实用的二阶精度的两层显格式.     

Rounded data analysis based on multi-layer ranked set sampling

Observations of sampling are often subject to rounding, but are modeled as though they were unrounded. This paper examines the impact of rounding errors on parameter estimation with multi-layer ranked set sampling. It shows that the rounding errors seriously distort the behavior of covariance matrix estimate, and lead to inconsistent estimation. Taking this into account, we present a new approach to implement the estimation for this model, and further establish the strong consistency and asymptotic normality of the proposed estimators. Simulation experiments show that our estimates based on rounded multi-layer ranked set sampling are always more efficient than those based on rounded simple random sampling.     

A numerical approach for solving weakly singular partial integro‐differential equations via two‐dimensional‐orthonormal Bernstein polynomials with the convergence analysis

In this paper, we develop an efficient matrix method based on two‐dimensional orthonormal Bernstein polynomials (2D‐OBPs) to provide approximate solution of linear and nonlinear weakly singular partial integro‐differential equations (PIDEs). First, we approximate all functions involved in the considerable problem via 2D‐OBPs. Then, by using the operational matrices of integration, differentiation, and product, the solution of Volterra singular PIDEs is transformed to the solution of a linear or nonlinear system of algebraic equations which can be solved via some suitable numerical methods. With a small number of bases, we can find a reasonable approximate solution. Moreover, we establish some useful theorems for discussing convergence analysis and obtaining an error estimate associated with the proposed method. Finally, we solve some illustrative examples by employing the presented method to show the validity, efficiency, high accuracy, and applicability of the proposed technique.     

Maximum‐norm a posteriori error bounds for a collocation method applied to a singularly perturbed reaction–diffusion problem in three dimensions

Collocation with triquadratic C¹‐splines for a singularly perturbed reaction–diffusion problem in three dimension is studied. A posteriori error bound in the maximum norm is derived for the collocation method on arbitrary tensor‐product meshes which is robust in the perturbation parameter. Numerical results are presented that support our theoretical estimate.     

A fitted numerical method for parabolic turning point singularly perturbed problems with an interior layer

The objective of this paper is to construct and analyze a fitted operator finite difference method (FOFDM) for the family of time‐dependent singularly perturbed parabolic convection–diffusion problems. The solution to the problems we consider exhibits an interior layer due to the presence of a turning point. We first establish sharp bounds on the solution and its derivatives. Then, we discretize the time variable using the classical Euler method. This results in a system of singularly perturbed interior layer two‐point boundary value problems. We propose a FOFDM to solve the system above. Through a rigorous error analysis, we show that the scheme is uniformly convergent of order one with respect to both time and space variables. Moreover, we apply Richardson extrapolation to enhance the accuracy and the order of convergence of the proposed scheme. Numerical investigations are carried out to demonstrate the efficacy and robustness of the scheme.     

Convergence analysis of virtual element methods for semilinear parabolic problems on polygonal meshes

In this article, we discuss and analyze new conforming virtual element methods (VEMs) for the approximation of semilinear parabolic problems on convex polygonal meshes in two spatial dimension. The spatial discretization is based on polynomial and suitable nonpolynomial functions, and a Euler backward scheme is employed for time discretization. The discrete formulation of both the proposed schemes—semidiscrete and fully discrete (with time discretization) is discussed in detail, and the unique solvability of the resulted schemes is discussed. A priori error estimates for the proposed schemes (semidiscrete and fully discrete) in H¹‐ and L²‐norms are derived under the assumption that the source term f is Lipschitz continuous. Some numerical experiments are conducted to illustrate the performance of the proposed scheme and to confirm the theoretical convergence rates.     

A posteriori error analysis of an augmented fully mixed formulation for the nonisothermal Oldroyd–Stokes problem

In this article, we consider an augmented fully mixed variational formulation that has been recently proposed for the nonisothermal Oldroyd–Stokes problem, and develop an a posteriori error analysis for the 2‐D and 3‐D versions of the associated mixed finite element scheme. More precisely, we derive two reliable and efficient residual‐based a posteriori error estimators for this problem on arbitrary (convex or nonconvex) polygonal and polyhedral regions. The reliability of the proposed estimators draws mainly upon the uniform ellipticity of the bilinear forms of the continuous formulation, suitable assumptions on the domain and the data, stable Helmholtz decompositions, and the local approximation properties of the Clément and Raviart–Thomas operators. On the other hand, inverse inequalities, the localization technique based on bubble functions, and known results from previous works are the main tools yielding the efficiency estimate. Finally, several numerical results confirming the properties of the a posteriori error estimators and illustrating the performance of the associated adaptive algorithms are reported.     

Stabilized mixed finite element methods for the Navier‐Stokes equations with damping

In this paper, the stabilized mixed finite element methods are presented for the Navier‐Stokes equations with damping. The existence and uniqueness of the weak solutions are proven by use of the Brouwer fixed‐point theorem. Then, optimal error estimates for the H¹‐norm and L²‐norm of the velocity and the L²‐norm of the pressure are derived. Moreover, on the basis of the optimal L²‐norm error estimate of the velocity, a stabilized two‐step method is proposed, which is more efficient than the usual stabilized methods. Finally, two numerical examples are implemented to confirm the theoretical analysis.     

On the finite element approximation of elliptic QVIs with noncoercive operators

In this paper, we extend the approach developed by the author for the standard finite element method in the L^∞‐norm of the noncoercive variational inequalities (VI) (Numer Funct Anal Optim.2015;36:1107‐1121.) to impulse control quasi‐variational inequality (QVI). We derive the optimal error estimate, combining the so‐called Bensoussan‐Lions Algorithm and the concept of subsolutions for VIs.