In this paper, a new extrapolation economy cascadic multigrid method is proposed to solve the image restoration model. The new method combines the new extrapolation formula and quadratic interpolation to design a nonl...In this paper, a new extrapolation economy cascadic multigrid method is proposed to solve the image restoration model. The new method combines the new extrapolation formula and quadratic interpolation to design a nonlinear prolongation operator, which provides more accurate initial values for the fine grid level. An edge preserving denoising operator is constructed to remove noise and preserve image edges. The local smoothing operator reduces the influence of staircase effect. The experiment results show that the new method not only improves the computational efficiency but also ensures good recovery quality.展开更多
The quasi-Shannon interval wavelet is constructed based on the interpolation wavelet theory, and an adaptive precise integration method, which is based on extrapolation method is presented for nonlinear ordinary diffe...The quasi-Shannon interval wavelet is constructed based on the interpolation wavelet theory, and an adaptive precise integration method, which is based on extrapolation method is presented for nonlinear ordinary differential equations ( ODEs). And then, an adaptive interval wavelet precise integration method (AIWPIM) for nonlinear partial differential equations(PDEs) is proposed. The numerical results show that the computational precision of AIWPIM is higher than that of the method constructed by combining the wavelet and the 4th Runge-Kutta method, and the computational amounts of these two methods are almost equal. For convenience, the Burgers equation is taken as an example in introducing this method, which is also valid for more general cases.展开更多
An efficient high-order immersed interface method (IIM) is proposed to solve two-dimensional (2D) heat problems with fixed interfaces on Cartesian grids, which has the fourth-order accuracy in the maximum norm in ...An efficient high-order immersed interface method (IIM) is proposed to solve two-dimensional (2D) heat problems with fixed interfaces on Cartesian grids, which has the fourth-order accuracy in the maximum norm in both time and space directions. The space variable is discretized by a high-order compact (HOC) difference scheme with correction terms added at the irregular points. The time derivative is integrated by a Crank-Nicolson and alternative direction implicit (ADI) scheme. In this case, the time accuracy is just second-order. The Richardson extrapolation method is used to improve the time accuracy to fourth-order. The numerical results confirm the convergence order and the efficiency of the method.展开更多
Assume that a convergent matrix sequence{A<sub>n</sub>}:A<sub>n</sub>→A(n→∞), A<sub>n</sub>,A∈C<sup>3×3</sup>.We want to form a new matrix sequence {H<sub&...Assume that a convergent matrix sequence{A<sub>n</sub>}:A<sub>n</sub>→A(n→∞), A<sub>n</sub>,A∈C<sup>3×3</sup>.We want to form a new matrix sequence {H<sub>n</sub>}, derived from {A<sub>n</sub>}, which has also A aslimit and whose convergence is faster than the of {A<sub>n</sub>}. Three rational extrapolation meth-ods for accelerating the convergence of matrix sequences {A<sub>n</sub>} are presented in this paper.The underlying methods are based on the generalized inverse for matrices which is展开更多
1. Introduction It is known that the following Cauchy problem for a parabolic partial differential equation (where the values at the right boundary, u.(1, t)=v(t) are unknown and sought for) is ill-posed: the solution...1. Introduction It is known that the following Cauchy problem for a parabolic partial differential equation (where the values at the right boundary, u.(1, t)=v(t) are unknown and sought for) is ill-posed: the solution (v) does not depend continuously on the data (g). In order to treat the ill-posedness and develop the numerical method, one reformulates the problem as a Volterra integral equation of the first kind wish a convolution type kernel (see Sneddon [1], Carslaw and Jaeger [2])展开更多
A mathematical model of two-phase fluid nonlinear flow in the direction of normal of ellipse through low-permeability porous media was established according to a nonlinear flow law expressed in a continuous function w...A mathematical model of two-phase fluid nonlinear flow in the direction of normal of ellipse through low-permeability porous media was established according to a nonlinear flow law expressed in a continuous function with three parameters, a mass conservation law and a concept of turbulent ellipses. A solution to the model was obtained by using a finite difference method and an extrapolation method. Formulas of calculating development index not only before but also after water breaks through an oil well in the condition of two-phase fluid nonlinear flow in the media were derived. An example was discussed. Water saturation distribution was presented. The moving law of drainage front was found. Laws of change of pressure difference with time were recognized. Results show that there is much difference of water saturation distribution between nonlinear flow and linear flow; that drainage front by water moves faster, water breaks through sooner and the index gets worse because of the nonlinear flow; and that dimensionless pressure difference gets larger at the same dimensionless time and difficulty of oil development becomes bigger by the nonlinear flow. Thus, it is necessary that influence of nonlinear flow on development indexes of the oil fields be taken into account. The results provide water-flooding development of the oilfields with scientific basis.展开更多
The frequency domain analysis of systems is an important topic in control theory. Powerful graphical tools exist in classic control, such as the Nyquist plot, Bode plots, and Nichols chart. These methods have been wid...The frequency domain analysis of systems is an important topic in control theory. Powerful graphical tools exist in classic control, such as the Nyquist plot, Bode plots, and Nichols chart. These methods have been widely used to evaluate the frequency domain behavior of system. A literature survey shows that various approaches are available for the computation of the frequency response of control systems under different types of parametric dependencies, such as affine, multi-linear, polynomial, etc. However, there is a lack of tools in the literature to construct the Bode envelopes for the general nonlinear type of parametric dependencies. In this paper, we address the problem of computation of the envelope of Bode frequency response of a non-rational transfer function with nonlinear parametric uncertainties varying over a box. We propose two techniques to compute the Bode envelopes:first, based on the natural interval extensions (NIE) combined with uniform subdivision and second, based on the existing Taylor model combined with subdivision strategy. We also propose the algorithms to further speed up both methods through extrapolation techniques.展开更多
This paper is concerlled with the investigation of a twrvparametric linear stationary iterative method, called Modified Extrapolated Jacobi (MEJ) method, for solving linear systems Ax = b, where A is a nonsingular con...This paper is concerlled with the investigation of a twrvparametric linear stationary iterative method, called Modified Extrapolated Jacobi (MEJ) method, for solving linear systems Ax = b, where A is a nonsingular consistently ordered 2-cyclic matrix. We give sufficient and necessary conditions for strong convergence of the MEJ method and we determine the optimum extrapolation parameters and the optimum spectral radius of it, in the case where all the efornvalues of the block Jacobi iteration matrir associated with A are real. In the last section, we compare the MEJ with other known methods.展开更多
The classical composite rectangle (constant) rule for the computation of Cauchy principle value integral with the singular kernel is discussed. We show that the superconvergence rate of the composite midpoint ru...The classical composite rectangle (constant) rule for the computation of Cauchy principle value integral with the singular kernel is discussed. We show that the superconvergence rate of the composite midpoint rule occurs at certain local coodinate of each subinterval and obtain the corresponding superconvergence error estimate. Then collation methods are presented to solve certain kind of Hilbert singular integral equation. At last, some numerical examples are provided to validate the theoretical analysis.展开更多
The steady, asymmetric and two-dimensional flow of viscous, incompressible and Newtonian fluid through a rectangular channel with splitter plate parallel to walls is investigated numerically. Earlier, the position of ...The steady, asymmetric and two-dimensional flow of viscous, incompressible and Newtonian fluid through a rectangular channel with splitter plate parallel to walls is investigated numerically. Earlier, the position of the splitter plate was taken as a centreline of channel but here it is considered its different positions which cause the asymmetric behaviour of the flow field. The geometric parameter that controls the position of splitter is defined as splitter position parameter a. The plane Poiseuille flow is considered far from upstream and downstream of the splitter. This flow-problem is solved numerically by a numerical scheme comprising a fourth order method, followed by a special finite-method. This numerical scheme transforms the governing equations to system of finite-difference equations, which are solved by point S.O.R. iterative method. In addition, the results obtained are further refined and upgraded by Richardson Extrapolation method. The calculations are carried out for the ranges -1 α R < 10<sup>5</sup>. The results are compared with existing literature regarding the symmetric case (when a = 0) for velocity, vorticity and skin friction distributions. The comparison is very favourable. Moreover, the notable thing is that the decay of vorticity to its downstream value takes place over an increasingly longer scale of x as R increases for symmetric case but it is not so for asymmetric one.展开更多
We discuss semiconvergence of the extrapolated iterative methods for solving singular linear systems. We obtain the upper bounds and the optimum convergence factor of the extrapolation method as well as its associated...We discuss semiconvergence of the extrapolated iterative methods for solving singular linear systems. We obtain the upper bounds and the optimum convergence factor of the extrapolation method as well as its associated optimum extrapolation parameter. Numerical examples are given to illustrate the theoretical results.展开更多
In this paper, we develop a new technique called multiplicative extrapolation method which is used to construct higher order schemes for ordinary differential equations. We call it a new method because we only see add...In this paper, we develop a new technique called multiplicative extrapolation method which is used to construct higher order schemes for ordinary differential equations. We call it a new method because we only see additive extrapolation method before. This new method has a great advantage over additive extrapolation method because it keeps group property. If this method is used to construct higher order schemes from lower symplectic schemes, the higher order ones are also symplectic. First we introduce the concept of adjoint methods and some of their properties. We show that there is a self-adjoint scheme corresponding to every method. With this self-adjoint scheme of lower order, we can construct higher order schemes by multiplicative extrapolation method, which can be used to construct even much higher order schemes. Obviously this constructing process can be continued to get methods of arbitrary even order.展开更多
Finite element approximation of the elliptic operator on a non-convexdomain composed of rectangles is considered using a graded mesh.Some errorestimates and error expansion are presented.
文摘In this paper, a new extrapolation economy cascadic multigrid method is proposed to solve the image restoration model. The new method combines the new extrapolation formula and quadratic interpolation to design a nonlinear prolongation operator, which provides more accurate initial values for the fine grid level. An edge preserving denoising operator is constructed to remove noise and preserve image edges. The local smoothing operator reduces the influence of staircase effect. The experiment results show that the new method not only improves the computational efficiency but also ensures good recovery quality.
文摘The quasi-Shannon interval wavelet is constructed based on the interpolation wavelet theory, and an adaptive precise integration method, which is based on extrapolation method is presented for nonlinear ordinary differential equations ( ODEs). And then, an adaptive interval wavelet precise integration method (AIWPIM) for nonlinear partial differential equations(PDEs) is proposed. The numerical results show that the computational precision of AIWPIM is higher than that of the method constructed by combining the wavelet and the 4th Runge-Kutta method, and the computational amounts of these two methods are almost equal. For convenience, the Burgers equation is taken as an example in introducing this method, which is also valid for more general cases.
基金supported by the National Natural Science Foundation of China(No.51174236)the National Basic Research Program of China(973 Program)(No.2011CB606306)the Opening Project of State Key Laboratory of Porous Metal Materials(No.PMM-SKL-4-2012)
文摘An efficient high-order immersed interface method (IIM) is proposed to solve two-dimensional (2D) heat problems with fixed interfaces on Cartesian grids, which has the fourth-order accuracy in the maximum norm in both time and space directions. The space variable is discretized by a high-order compact (HOC) difference scheme with correction terms added at the irregular points. The time derivative is integrated by a Crank-Nicolson and alternative direction implicit (ADI) scheme. In this case, the time accuracy is just second-order. The Richardson extrapolation method is used to improve the time accuracy to fourth-order. The numerical results confirm the convergence order and the efficiency of the method.
基金The works is supported by the National Natural Science Foundation of China(19871054)
文摘Assume that a convergent matrix sequence{A<sub>n</sub>}:A<sub>n</sub>→A(n→∞), A<sub>n</sub>,A∈C<sup>3×3</sup>.We want to form a new matrix sequence {H<sub>n</sub>}, derived from {A<sub>n</sub>}, which has also A aslimit and whose convergence is faster than the of {A<sub>n</sub>}. Three rational extrapolation meth-ods for accelerating the convergence of matrix sequences {A<sub>n</sub>} are presented in this paper.The underlying methods are based on the generalized inverse for matrices which is
文摘1. Introduction It is known that the following Cauchy problem for a parabolic partial differential equation (where the values at the right boundary, u.(1, t)=v(t) are unknown and sought for) is ill-posed: the solution (v) does not depend continuously on the data (g). In order to treat the ill-posedness and develop the numerical method, one reformulates the problem as a Volterra integral equation of the first kind wish a convolution type kernel (see Sneddon [1], Carslaw and Jaeger [2])
文摘A mathematical model of two-phase fluid nonlinear flow in the direction of normal of ellipse through low-permeability porous media was established according to a nonlinear flow law expressed in a continuous function with three parameters, a mass conservation law and a concept of turbulent ellipses. A solution to the model was obtained by using a finite difference method and an extrapolation method. Formulas of calculating development index not only before but also after water breaks through an oil well in the condition of two-phase fluid nonlinear flow in the media were derived. An example was discussed. Water saturation distribution was presented. The moving law of drainage front was found. Laws of change of pressure difference with time were recognized. Results show that there is much difference of water saturation distribution between nonlinear flow and linear flow; that drainage front by water moves faster, water breaks through sooner and the index gets worse because of the nonlinear flow; and that dimensionless pressure difference gets larger at the same dimensionless time and difficulty of oil development becomes bigger by the nonlinear flow. Thus, it is necessary that influence of nonlinear flow on development indexes of the oil fields be taken into account. The results provide water-flooding development of the oilfields with scientific basis.
文摘The frequency domain analysis of systems is an important topic in control theory. Powerful graphical tools exist in classic control, such as the Nyquist plot, Bode plots, and Nichols chart. These methods have been widely used to evaluate the frequency domain behavior of system. A literature survey shows that various approaches are available for the computation of the frequency response of control systems under different types of parametric dependencies, such as affine, multi-linear, polynomial, etc. However, there is a lack of tools in the literature to construct the Bode envelopes for the general nonlinear type of parametric dependencies. In this paper, we address the problem of computation of the envelope of Bode frequency response of a non-rational transfer function with nonlinear parametric uncertainties varying over a box. We propose two techniques to compute the Bode envelopes:first, based on the natural interval extensions (NIE) combined with uniform subdivision and second, based on the existing Taylor model combined with subdivision strategy. We also propose the algorithms to further speed up both methods through extrapolation techniques.
文摘This paper is concerlled with the investigation of a twrvparametric linear stationary iterative method, called Modified Extrapolated Jacobi (MEJ) method, for solving linear systems Ax = b, where A is a nonsingular consistently ordered 2-cyclic matrix. We give sufficient and necessary conditions for strong convergence of the MEJ method and we determine the optimum extrapolation parameters and the optimum spectral radius of it, in the case where all the efornvalues of the block Jacobi iteration matrir associated with A are real. In the last section, we compare the MEJ with other known methods.
文摘The classical composite rectangle (constant) rule for the computation of Cauchy principle value integral with the singular kernel is discussed. We show that the superconvergence rate of the composite midpoint rule occurs at certain local coodinate of each subinterval and obtain the corresponding superconvergence error estimate. Then collation methods are presented to solve certain kind of Hilbert singular integral equation. At last, some numerical examples are provided to validate the theoretical analysis.
文摘The steady, asymmetric and two-dimensional flow of viscous, incompressible and Newtonian fluid through a rectangular channel with splitter plate parallel to walls is investigated numerically. Earlier, the position of the splitter plate was taken as a centreline of channel but here it is considered its different positions which cause the asymmetric behaviour of the flow field. The geometric parameter that controls the position of splitter is defined as splitter position parameter a. The plane Poiseuille flow is considered far from upstream and downstream of the splitter. This flow-problem is solved numerically by a numerical scheme comprising a fourth order method, followed by a special finite-method. This numerical scheme transforms the governing equations to system of finite-difference equations, which are solved by point S.O.R. iterative method. In addition, the results obtained are further refined and upgraded by Richardson Extrapolation method. The calculations are carried out for the ranges -1 α R < 10<sup>5</sup>. The results are compared with existing literature regarding the symmetric case (when a = 0) for velocity, vorticity and skin friction distributions. The comparison is very favourable. Moreover, the notable thing is that the decay of vorticity to its downstream value takes place over an increasingly longer scale of x as R increases for symmetric case but it is not so for asymmetric one.
基金supported by the National Natural Science Foundation of China under grant 10371056the Foundation for the Authors of the National Excellent Doctoral Thesis Award of China under grant 200720+1 种基金the Natural Science Foundation of Jiangsu Province under grant BK2006725the College Natural Science Foundation of Jiangsu Province under grant 05KJB110062
文摘We discuss semiconvergence of the extrapolated iterative methods for solving singular linear systems. We obtain the upper bounds and the optimum convergence factor of the extrapolation method as well as its associated optimum extrapolation parameter. Numerical examples are given to illustrate the theoretical results.
文摘In this paper, we develop a new technique called multiplicative extrapolation method which is used to construct higher order schemes for ordinary differential equations. We call it a new method because we only see additive extrapolation method before. This new method has a great advantage over additive extrapolation method because it keeps group property. If this method is used to construct higher order schemes from lower symplectic schemes, the higher order ones are also symplectic. First we introduce the concept of adjoint methods and some of their properties. We show that there is a self-adjoint scheme corresponding to every method. With this self-adjoint scheme of lower order, we can construct higher order schemes by multiplicative extrapolation method, which can be used to construct even much higher order schemes. Obviously this constructing process can be continued to get methods of arbitrary even order.
基金supported by the Deutsche Forschungsgemeinschaft(DFG),SFB 123."Stochatistische Mathematische Modelle",Universitat Heidelberg
文摘Finite element approximation of the elliptic operator on a non-convexdomain composed of rectangles is considered using a graded mesh.Some errorestimates and error expansion are presented.