Extrapolated Tikhonov method and inversion of 3D density images of gravity data

Tikhonov regularization（TR） method has played a very important role in the gravity data and magnetic data process. In this paper, the Tikhonov regularization method with respect to the inversion of gravity data is discussed. and the extrapolated TR method（EXTR） is introduced to improve the fitting error. Furthermore, the effect of the parameters in the EXTR method on the fitting error, number of iterations, and inversion results are discussed in details. The computation results using a synthetic model with the same and different densities indicated that. compared with the TR method, the EXTR method not only achieves the a priori fitting error level set by the interpreter but also increases the fitting precision, although it increases the computation time and number of iterations. And the EXTR inversion results are more compact than the TR inversion results, which are more divergent. The range of the inversion data is closer to the default range of the model parameters, and the model features and default model density distribution agree well.

SnO_2-based solid solutions for CH_4 deep oxidation: Quantifying the lattice capacity of SnO_2 using an X-ray diffraction extrapolation method

A series of SnO2‐based catalysts modified by Mn, Zr, Ti and Pb oxides with a Sn/M （M=Mn, Zr, Ti and Pb） molar ratio of 9/1 were prepared by a co‐precipitation method and used for CH4 and CO oxidation. The Mn3＋, Zr4＋, Ti4＋and Pb4＋cations are incorporated into the lattice of tetragonal rutile SnO2 to form a solid solution structure. As a consequence, the surface area and thermal stability of the catalysts are improved. Moreover, the oxygen species of the modified catalysts become easier to be reduced. Therefore, the oxidation activity over the catalysts was improved, except for the one modified by Pb oxide. Manganese oxide demonstrates the best promotional effects for SnO2. Using an X‐ray diffraction extrapolation method, the lattice capacity of SnO2 for Mn2O3 was 0.135 g Mn2O3/g SnO2, which indicates that to form stable solid solution, only 21%Sn4＋cations in the lattice can be maximally replaced by Mn3＋. If the amount of Mn3＋cations is over the capacity, Mn2O3 will be formed, which is not favorable for the activity of the catalysts. The Sn rich samples with only Sn‐Mn solid solution phase show higher activity than the ones with excess Mn2O3 species.

ADAPTIVE INTERVAL WAVELET PRECISE INTEGRATION METHOD FOR PARTIAL DIFFERENTIAL EQUATIONS

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.

Fatigue life evaluation of girth butt weld within welded cast steel joints based on the extrapolation notch stress method

The fatigue life evaluation of the girth butt weld within the welded cast steel joint was studied based on the extrapolation notch stress method.Firstly,the mesh sensitivity of the finite element model of the welded cast steel joint was analyzed to determine the optimal mesh size.Based on the stress field analysis of the finite element model of the welded cast steel joint at the weld toe and weld root,the sharp model of the extrapolation notch stress method was applied to derive the effective notch stress of the rounded model belonging to the effective notch stress method,in which the key problem is to calculate the extrapolation point C,and the extrapolation point C has an exponential function relationship with the geometric parameters of the welded cast steel joint.By setting different values of geometric parameters,the corresponding value of parameter C is calculated,and then the functional relationship between the extrapolation point C and the geometric parameters can be obtained by the multiple linear regression analysis.Meanwhile,the fatigue life evaluation of the girth butt weld within welded cast steel joints based on the effective notch stress was performed according to the guideline recommended by the IIW(International Institute of Welding).The results indicate that the extrapolation notch stress method can effectively simplify the process of calculating the effective notch stress and accurately evaluate the fatigue life of the girth butt weld within welded cast steel joints.

GENERALIZED INVERSE RATIONAL EXTRAPOLATION METHODS FOR MATRIX SEQUENCES

Assume that a convergent matrix sequence{An}:An→A（n→∞）, An,A∈C3×3.We want to form a new matrix sequence {Hn}, derived from {An}, which has also A aslimit and whose convergence is faster than the of {An}. Three rational extrapolation meth-ods for accelerating the convergence of matrix sequences {An} are presented in this paper.The underlying methods are based on the generalized inverse for matrices which is

Efficient high-order immersed interface methods for heat equations with interfaces

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.

EXTRAPOLATION FOR COLLOCATION METHOD OF THE FIRST KIND VOLTERRA INTEGRAL EQUATIONS

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 New Extrapolation Economy Cascadic Multigrid Method for Image Restoration Problems

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.

MATHEMATICAL MODEL OF TWO-PHASE FLUID NONLINEAR FLOW IN LOW-PERMEABILITY POROUS MEDIA WITH APPLICATIONS

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.

Construction of Bode Envelopes Using REP Based Range Finding Algorithms

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.

Investigation of lattice capacity effect on Cu2+-doped SnO2 solid solution catalysts to promote reaction performance toward NOx-SCR with NH3

To understand the effect of the doping amount of Cu^2+ on the structure and reactivity of SnO2 in NOx-SCR with NH3, a series of Sn-Cu-O binary oxide catalysts with different Sn/Cu ratios have been prepared and thoroughly characterized. Using the XRD extrapolation method, the SnO2 lattice capacity for Cu^2+ cations is determined at 0.10 g Cu O per g of SnO2, equaling a Sn/Cu molar ratio of 84/16. Therefore, in a tetragonal rutile SnO2 lattice, only a maximum of 16% of the Sn4+ cations can be replaced by Cu^2+ to form a stable solid solution structure. If the Cu content is higher, Cu O will form on the catalyst surface, which has a negative effect on the reaction performance. For samples in a pure solid solution phase, the number of surface defects increase with increasing Cu content until it reaches the lattice capacity, as confirmed by Raman spectroscopy. As a result, the amounts of both active oxygen species and acidic sites on the surface, which critically determine the reaction performance, also increase and reach the maximum level for the catalyst with a Cu content close to the lattice capacity. A distinct lattice capacity threshold effect on the structure and reactivity of Sn-Cu binary oxide catalysts has been observed. A Sn-Cu catalyst with the best reaction performance can be obtained by doping the SnO2 matrix with the lattice capacity amount of Cu^2+.

The Rectangle Rule for Computing Cauchy Principal Value Integral on Circle

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.

Numerical Simulation of Asymmetric Merging Flow in a Rectangular Channel

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 < 105. 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.

Extrapolation cascadic multigrid method on piecewise uniform grid

The triangular linear finite elements on piecewise uniform grid for an elliptic problem in convex polygonal domain are discussed. Global superconvergence in discrete Hi-norm and global extrapolation in discrete L2-norm are proved. Based on these global estimates the conjugate gradient method （CG） is effective, which is applied to extrapolation cascadic multigrid method （EXCMG）. The numerical experiments show that EXCMG is of the global higher accuracy for both function and gradient.

ON THE OPTIMIZATION OF EXTRAPOLATION METHODS FOR SINGULAR LINEAR SYSTEMS

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.

OPTIMUM MODIFIED EXTRAPOLATED JACOBI METHOD FOR CONSISTENTLY ORDERED MATRICES

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.

MULTIPLICATIVE EXTRAPOLATION METHOD FOR CONSTRUCTING HIGHER ORDER SCHEMES FOR ORDINARY DIFFERENTIAL EQUATIONS

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 USING A GRADED MESH ON DOMAINS WITH REENTRANT CORNERS

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.