THE NUMERICAL STABILITY OF THE BLOCK θ-METHODS FOR DELAY DIFFERENTIAL EQUATIONS

This paper focuses on the numerical stability of the block θ methods adapted to differential equations with a delay argument. For the block θ methods, an interpolation procedure is introduced which leads to the numerical processes that satisfy an important asymptotic stability condition related to the class of test problems y′(t)=ay(t)+by(t-τ) with a,b∈C, Re(a)<-|b| and τ>0. We prove that the block θ method is GP stable if and only if the method is A stable for ordinary differential equations. Furthermore, it is proved that the P and GP stability are equivalent for the block θ method.

THE NUMERICAL STABILITIES OF MULTIDERIVATIVE BLOCK METHOD

In [1], a class of multiderivative block methods (MDBM) was studied for the numerical solutions of stiff ordinary differential equations. This paper is aimed at solving the problem proposed in [1] that what conditions should be fulfilled for MDBMs in order to guarantee the A-stabilities. The explicit expressions of the polynomialsP(h) and Q(h) in the stability functions h(h)=P(h)/Q(h)are given. Furthermore, we prove P(-h)-Q(h). With the aid of symbolic computations and the expressions of diagonal Fade approximations, we obtained the biggest block size k of the A-stable MDBM for any given l (the order of the highest derivatives used in MDBM,l>1)

High Order Block Method for Third Order ODEs

Many initial value problems are difficult to be solved using ordinary,explicit step-by-step methods because most of these problems are considered stiff.Certain implicit methods,however,are capable of solving stiff ordinary differential equations(ODEs)usually found in most applied problems.This study aims to develop a new numerical method,namely the high order variable step variable order block backward differentiation formula(VSVOHOBBDF)for the main purpose of approximating the solutions of third order ODEs.The computational work of the VSVO-HOBBDF method was carried out using the strategy of varying the step size and order in a single code.The order of the proposed method was then discussed in detail.The advancement of this strategy is intended to enhance the efficiency of the proposed method to approximate solutions effectively.In order to confirm the efficiency of the VSVO-HOBBDF method over the two ODE solvers in MATLAB,particularly ode15s and ode23s,a numerical experiment was conducted on a set of stiff problems.The numerical results prove that for this particular set of problem,the use of the proposed method is more efficient than the comparable methods.VSVO-HOBBDF method is thus recommended as a reliable alternative solver for the third order ODEs.

BLOCK PIVOT METHODS FOR SOLVING FRICTIONAL CONTACT PROBLEMS

Based on elementary group theory, the block pivot methods for solving two-dimensional elastic frictional contact problems are presented in this paper. It is proved that the algorithms converge within a finite number of steps when the friction coefficient is ''relative small''. Unlike most mathematical programming methods for contact problems, the block pivot methods permit multiple exchanges of basic and nonbasic variables.

KRYLOV SUBSPACE PROJECTION METHOD AND ITS APPLICATION ON OIL RESERVOIR SIMULATION

Krylov subspace projection methods are known to be highly efficient for solving large linear systems. Many different versions arise from different choices to the left and right subspaces. These methods were classified into two groups in terms of the different forms of matrix H-m, the main properties in applications and the new versions of these two types of methods were briefly reviewed, then one of the most efficient versions, GMRES method was applied to oil reservoir simulation. The block Pseudo-Elimination method was used to generate the preconditioned matrix. Numerical results show much better performance of this preconditioned techniques and the GMRES method than that of preconditioned ORTHMIN method, which is now in use in oil reservoir simulation. Finally, some limitations of Krylov subspace methods and some potential improvements to this type of methods are further presented.

Runge-Kutta Method and Bolck by Block Method to Solve Nonlinear Fredholm-Volterra Integral Equation with Continuous Kernel

In this paper, the existence and uniqueness of the solution of Fredholm-Volterra integral equation is considered (NF-VIE) with continuous kernel;then we used a numerical method to reduce this type of equations to a system of nonlinear Volterra integral equations. Runge-Kutta method (RKM) and Bolck by block method (BBM) are used to solve the system of nonlinear Volterra integral equations of the second kind (SNVIEs) with continuous kernel. The error in each case is calculated.

Applications of Bootstrap in Analyzing General Extreme Value Distributions

The bootstrap method is one of the new ways of studying statistical math which this article uses but is a major tool for studying and evaluating the values of parameters in probability distribution.Our research is concerned overview of the theory of infinite distribution functions.The tool to deal with the problems raised in the paper is the mathematical methods of random analysis(theory of random process and multivariate statistics).In this article,we introduce the new function to find out the bias and standard error with jackknife method for Generalized Extreme Value distributions.

SOLVING NONLINEAR DELAY-DIFFERENTIAL-ALGEBRAIC EQUATIONS WITH SINGULAR PERTURBATION VIA BLOCK BOUNDARY VALUE METHODS

Block boundary value methods(BBVMs)are extended in this paper to obtain the numerical solutions of nonlinear delay-differential-algebraic equations with singular perturbation(DDAESP).It is proved that the extended BBVMs in some suitable conditions are globally stable and can obtain a unique exact solution of the DDAESP.Besides,whenever the classic Lipschitz conditions are satisfied,the extended BBVMs are preconsistent and pth order consistent.Moreover,through some numerical examples,the correctness of the theoretical results and computational validity of the extended BBVMs is further confirmed.

On Stochastic Error and Computational Efficiency of the Markov Chain Monte Carlo Method

InMarkov ChainMonte Carlo(MCMC)simulations,thermal equilibria quantities are estimated by ensemble average over a sample set containing a large number of correlated samples.These samples are selected in accordance with the probability distribution function,known from the partition function of equilibrium state.As the stochastic error of the simulation results is significant,it is desirable to understand the variance of the estimation by ensemble average,which depends on the sample size(i.e.,the total number of samples in the set)and the sampling interval(i.e.,cycle number between two consecutive samples).Although large sample sizes reduce the variance,they increase the computational cost of the simulation.For a given CPU time,the sample size can be reduced greatly by increasing the sampling interval,while having the corresponding increase in variance be negligible if the original sampling interval is very small.In this work,we report a few general rules that relate the variance with the sample size and the sampling interval.These results are observed and confirmed numerically.These variance rules are derived for theMCMCmethod but are also valid for the correlated samples obtained using other Monte Carlo methods.The main contribution of this work includes the theoretical proof of these numerical observations and the set of assumptions that lead to them.

DIRECT ITERATIVE METHODS FOR RANK DEFICIENT GENERALIZED LEAST SQUARES PROBLEMS

The generalized least squares (LS) problem ... (Ax - b)[sup TW[sup -1](Ax - b) appears in, many application areas. Here W is an m × m symmetric positive definite matrix and A is an m × n matrix with m ≥ n. Since the problem has many solutions in rank deficient case, some special preconditioned techniques are adapted to obtain the minimum 2-norm solution. A block SOR method and the preconditioned conjugate gradient (PCG) method are proposed here. Convergence and optimal relaxation parameter for the block SOR method are studied. An error bound for the PCG method is given. The comparison of these methods is investigated. Some remarks on the implementation of the methods and the operation cost are given as well. [ABSTRACT FROM AUTHOR]

Text extraction method for historical Tibetan document images based on block projections

Text extraction is an important initial step in digitizing the historical documents. In this paper, we present a text extraction method for historical Tibetan document images based on block projections. The task of text extraction is considered as text area detection and location problem. The images are divided equally into blocks and the blocks are filtered by the information of the categories of connected components and corner point density. By analyzing the filtered blocks' projections, the approximate text areas can be located, and the text regions are extracted. Experiments on the dataset of historical Tibetan documents demonstrate the effectiveness of the proposed method.

Numerical Solution for Neutral Delay Differential Equation of Constant or Proportional Type using Hybrid Block Method

Neutral Delay Differential Equation(NDDE)is a differential problem that has regularly existed in numerous occurrences and has represented a significant role in dealing with real-life phenomena,especially on their application in biological and physiological processes.A fifth-order two-point hybrid implicit multistep block method(2PIH5)has been formulated in this research for the numerical solution of Neutral Delay Differential Equation(NDDE).A Taylor series interpolation polynomial has been implemented in the formulation of the proposed 2PIH5.The order,consistency,and zero-stability for 2PIH5 have been illustrated.The analyses of convergence and stability test have been performed and discussed.The initial value problems for the first-order NDDE with constant or proportional delay have been solved using the proposed block method.Some numerical results for the proposed method have been presented to prove the adaptability and applicability of the proposed method in solving NDDE.The proposed method is proved to be comparable with the other existing methods.It is assumed to be reliable and efficient for solving the first-order NDDE with constant or proportional delay.

A Parallel Line Search Subspace Correction Method for Composite Convex Optimization

In this paper,we investigate a parallel subspace correction framework for composite convex optimization.The variables are first divided into a few blocks based on certain rules.At each iteration,the algorithms solve a suitable subproblem on each block simultaneously,construct a search direction by combining their solutions on all blocks,then identify a new point along this direction using a step size satisfying the Armijo line search condition.They are called PSCLN and PSCLO,respectively,depending on whether there are overlapping regions between two imme-diately adjacent blocks of variables.Their convergence is established under mild assumptions.We compare PSCLN and PSCLO with the parallel version of the fast iterative thresholding algorithm and the fixed-point continuation method using the Barzilai-Borwein step size and the greedy coordinate block descent method for solving the l1-regularized minimization problems.Our numerical results showthatPSCLN andPSCLOcan run fast and return solutions notworse than those from the state-of-theart algorithms on most test problems.It is also observed that the overlapping domain decomposition scheme is helpful when the data of the problem has certain special structures.

A Spectral Time-Domain Method for Computational Electrodynamics

Ever since its introduction by Kane Yee over forty years ago,the finitedifference time-domain(FDTD)method has been a widely-used technique for solving the time-dependent Maxwell’s equations that has also inspired many other methods.This paper presents an alternative approach to these equations in the case of spatially-varying electric permittivity and/or magnetic permeability,based on Krylov subspace spectral(KSS)methods.These methods have previously been applied to the variable-coefficient heat equation and wave equation,and have demonstrated high-order accuracy,as well as stability characteristic of implicit timestepping schemes,even though KSS methods are explicit.KSS methods for scalar equations compute each Fourier coefficient of the solution using techniques developed by Golub and Meurant for approximating elements of functions of matrices by Gaussian quadrature in the spectral,rather than physical,domain.We show how they can be generalized to coupled systems of equations,such as Maxwell’s equations,by choosing appropriate basis functions that,while induced by this coupling,still allow efficient and robust computation of the Fourier coefficients of each spatial component of the electric and magnetic fields.We also discuss the application of block KSS methods to problems involving non-self-adjoint spatial differential operators,which requires a generalization of the block Lanczos algorithm of Golub and Underwood to unsymmetric matrices.

NonlinearAsynchronousBlockIterativeMethod

This paper proposes a class of asynchronous block iterative methods for solving large scale nonlinear equations F(x)=0 and proves local convergence. This method splits F into p blocks, then does the asynchronous parallel iteration on the p multiprocessor with shared memory. Because each processor need only solve equations with a low dimension and there is no synchronous waiting time, the parallel efficiency can be increased. Finally, we give the results of the numerical test of three kinds of Newton like asynchronous block iteration methods which run well on a multiprocessor system. These results show that the parallel efficiency is very high.

On the Linear Convergence of the Approximate Proximal Splitting Method for Non-smooth Convex Optimization

Consider the problem of minimizing the sum of two convex functions,one being smooth and the other non-smooth.In this paper,we introduce a general class of approximate proximal splitting(APS)methods for solving such minimization problems.Methods in the APS class include many well-known algorithms such as the proximal splitting method,the block coordinate descent method(BCD),and the approximate gradient projection methods for smooth convex optimization.We establish the linear convergence of APS methods under a local error bound assumption.Since the latter is known to hold for compressive sensing and sparse group LASSO problems,our analysis implies the linear convergence of the BCD method for these problems without strong convexity assumption.

A VARIATION ON THE BLOCK ARNOLDIMETHOD FOR LARGE UNSYMMETRIC MATRIX EIGENPROBLEMS

The approximate eigenvectors or Ritz vectors obtained by the block Arnoldi method may converge very slowly and even fail to converge even if the approximate eigenvalues do. In order to improve the quality of the Ritz vectors, a modified strategy is proposed such that new approximate eigenvectors are certain combinations of the Ritz vectors and the waSted (m+1) th block basis vector and their corresponding residual norms are minimized in a certain sense. They can be cheaply computed by solving a few small 'dimensional minimization problems. The resulting modified m-step block Arnoldi method is better than the standard m-step one in theory and cheaper than the standard (m+1)-step one. Based on this strategy, a modified m-step iterative block Arnoldi algorithm is presented. Numerical experiments are reported to show that the modified m-step algorithm is often considerably more efficient than the standard (m+1)-step iterative one.

Plume aerodynamic effects of cushion engine in lunar landing

During the second period of China ＂Tanyue＂ Project, the explorer will softland on the moon. The cushion engines are used to decelerate the explorer and reduce the impact on the lunar ground. It is necessary to study its plume effects on the explorer component. The self-developed PWS （Plume WorkStation） software based on direct simulation Monte Carlo （DSMC） method is used to simulate the plume effects of two 150 N engines. Due to the complex structure of the explorer, PWS uses a decoupling method to treat the boundary mesh, which mainly interacts with simulation particles, and has no relation with the computational grids. After the analytical expressions of plane surfaces and curved surfaces of each boundary block are given, the particle position within or without the boundary blocks can be easily determined. Finally the 3D plume field of two 150 N engines is simulated. The pressure, temperature and velocity distributions of plume field are clearly presented by three characteristic slices. The aerodynamic effects on the explorer bottom, the landfall legs and antenna are separately shown. The compression influence on the plume flow of four landfall legs can be observed.

A Two-stage Adaptive Robust Model for Residential Micro-CHP Expansion Planning

This paper addresses the planning problem of residential micro combined heat and power (micro-CHP) systems (including micro-generation units, auxiliary boilers, and thermal storage tanks) considering the associated technical and economic factors. Since the accurate values of the thermal and electrical loads of this system cannot be exactly predicted for the planning horizon, the thermal and electrical load uncertainties are modeled using a two-stage adaptive robust optimization method based on a polyhedral uncertainty set. A solution method, which is composed of column-and-constraint generation (C&CG) algorithm and block coordinate descent (BCD) method, is proposed to efficiently solve this adaptive robust optimization model. Numerical results from a practical case study show the effective performance of the proposed adaptive robust model for residential micro-CHP planning and its solution method.

ON PRECONDITIONING OF INCOMPRESSIBLE NON-NEWTONIAN FLOW PROBLEMS

This paper deals with fast and reliable numerical solution methods for the incompress- ible non-Newtonian Navier-Stokes equations. To handle the nonlinearity of the governing equations, the Picard and Newton methods are used to linearize these coupled partial dif- ferential equations. For space discretization we use the finite element method and utilize the two-by-two block structure of the matrices in the arising algebraic systems of equa- tions. The Krylov subspace iterative methods are chosen to solve the linearized discrete systems and the development of computationally and numerically efficient preconditioners for the two-by-two block matrices is the main concern in this paper. In non-Newtonian flows, the viscosity is not constant and its variation is an important factor that effects the performance of some already known preconditioning techniques. In this paper we examine the performance of several preconditioners for variable viscosity applications, and improve them further to be robust with respect to variations in viscosity.