Adaptive multi-step piecewise interpolation reproducing kernel method for solving the nonlinear time-fractional partial differential equation arising from financial economics

This paper is aimed at solving the nonlinear time-fractional partial differential equation with two small parameters arising from option pricing model in financial economics.The traditional reproducing kernel(RK)method which deals with this problem is very troublesome.This paper proposes a new method by adaptive multi-step piecewise interpolation reproducing kernel(AMPIRK)method for the first time.This method has three obvious advantages which are as follows.Firstly,the piecewise number is reduced.Secondly,the calculation accuracy is improved.Finally,the waste time caused by too many fragments is avoided.Then four numerical examples show that this new method has a higher precision and it is a more timesaving numerical method than the others.The research in this paper provides a powerful mathematical tool for solving time-fractional option pricing model which will play an important role in financial economics.

Analysis of Landscape Vitality of Historical and Cultural Blocks Based on AHPFuzzy Comprehensive Evaluation Method:A Case Study of Daopashi Street in Anqing City

Historical and cultural blocks are witnesses of history and inheritors of culture. As one of the main spaces for outdoor interaction in historical and cultural blocks, the improvement of its vitality is of great significance for the improvement of residential environment and the better inheritance of history and culture. Taking Daopashi Street in Anqing City as an example, an evaluation model of landscape spatial vitality of historical and cultural blocks was constructed from three aspects of viewing function, store status and service facilities, and analytic hierarchy process was used to determine the index weight and vaguely evaluate the landscape spatial vitality of historical and cultural blocks. The results show that through the comparison of weight, architectural style(0.317), the practicability of service facilities(0.168) and plant landscape(0.165) had a significant impact on the landscape spatial vitality of historical and cultural blocks,and the landscape spatial vitality of historical and cultural blocks in Daopashi Street in Anqing City was at a good level.

Investigation of the block toppling evolution of a layered model slope by centrifuge test and discrete element modeling

Primary toppling usually occurs in layered rock slopes with large anti-dip angles.In this paper,the block toppling evolution was explored using a large-scale centrifuge system.Each block column in the layered model slope was made of cement mortar.Some artificial cracks perpendicular to the block column were prefabricated.Strain gages,displacement gages,and high-speed camera measurements were employed to monitor the deformation and failure processes of the model slope.The centrifuge test results show that the block toppling evolution can be divided into seven stages,i.e.layer compression,formation of major tensile crack,reverse bending of the block column,closure of major tensile crack,strong bending of the block column,formation of failure zone,and complete failure.Block toppling is characterized by sudden large deformation and occurs in stages.The wedge-shaped cracks in the model incline towards the slope.Experimental observations show that block toppling is mainly caused by bending failure rather than by shear failure.The tensile strength also plays a key factor in the evolution of block toppling.The simulation results from discrete element method(DEM)is in line with the testing results.Tensile stress exists at the backside of rock column during toppling deformation.Stress concentration results in the fragmented rock column and its degree is the most significant at the slope toe.

Failure evaluation and control factor analysis of slope block instability along traffic corridor in Southeastern Tibet

The instability of slope blocks occurred frequently along traffic corridor in Southeastern Tibet(TCST),which was primarily controlled by the rock mass structures.A rapid method evaluating the control effects of rock mass structures was proposed through field statistics of the slopes and rock mass structures along TCST,which combined the stereographic projection method,modified M-JCS model,and limit equilibrium theory.The instabilities of slope blocks along TCST were then evaluated rapidly,and the different control factors of instability were analyzed.Results showed that the probabilities of toppling(5.31%),planar(16.15%),and wedge(35.37%)failure of slope blocks along TCST increased sequentially.These instability modes were respectively controlled by the anti-dip joint,the joint parallel to slope surface with a dip angle smaller than the slope angle(singlejoint),and two groups of joints inclined out of the slope(double-joints).Regarding the control effects on slope block instability,the stabilization ability of doublejoints(72.7%),anti-dip joint(67.4%),and single-joint(57.6%)decreased sequentially,resulting in different probabilities of slope block instability.Additionally,nearby regional faults significantly influenced the joints,leading to spatial heterogeneity and segmental clustering in the stabilization ability provided by joints to the slope blocks.Consequently,the stability of slope blocks gradually weakened as they approached the fault zones.This paper can provide guidance and assistance for investigating the development characteristics of rock mass structures and the stability of slope blocks.

BLOCK BIDIAGONALIZATION METHODS FOR MULTIPLE NONSYMMETRIC LINEAR SYSTEMS

The symmetric linear system gives us many simplifications and a possibility to adapt the computations to the computer at hand in order to achieve better performance. The aim of this paper is to consider the block bidiagonalization methods derived from a symmetric augmented multiple linear systems and make a comparison with the block GMRES and block biconjugate gradient methods.

Application of Multi-Step Differential Transform Method on Flow of a Second-Grade Fluid over a Stretching or Shrinking Sheet

In this study, a reliable algorithm to develop approximate solutions for the problem of fluid flow over a stretching or shrinking sheet is proposed. It is depicted that the differential transform method (DTM) solutions are only valid for small values of the independent variable. The DTM solutions diverge for some differential equations that extremely have nonlinear behaviors or have boundary-conditions at infinity. For this reason the governing boundary-layer equations are solved by the Multi-step Differential Transform Method (MDTM). The main advantage of this method is that it can be applied directly to nonlinear differential equations without requiring linearization, discretization, or perturbation. It is a semi analytical-numerical technique that formulizes Taylor series in a very different manner. By applying the MDTM the interval of convergence for the series solution is increased. The MDTM is treated as an algorithm in a sequence of intervals for finding accurate approximate solutions for systems of differential equations. It is predicted that the MDTM can be applied to a wide range of engineering applications.

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.

Two-Step Hybrid Block Method for Solving Nonlinear Jerk Equations

In this paper, a block method with one hybrid point for solving Jerk equations is presented. The hybrid point is chosen to optimize the local truncation errors of the main formulas for the solution and the derivative at the end of the block. Analysis of the method is discussed, and some numerical examples show that the proposed method is efficient and accurate.

Automatic identification of rock discontinuity and stability analysis of tunnel rock blocks using terrestrial laser scanning

Local geometric information and discontinuity features are key aspects of the analysis of the evolution and failure mechanisms of unstable rock blocks in rock tunnels.This study demonstrates the integration of terrestrial laser scanning(TLS)with distinct element method for rock mass characterization and stability analysis in tunnels.TLS records detailed geometric information of the surrounding rock mass by scanning and collecting the positions of millions of rock surface points without contact.By conducting a fuzzy K-means method,a discontinuity automatic identification algorithm was developed,and a method for obtaining the geometric parameters of discontinuities was proposed.This method permits the user to visually identify each discontinuity and acquire its spatial distribution features(e.g.occurrences,spac-ings,trace lengths)in great detail.Compared with hand mapping in conventional geotechnical surveys,the geometric information of discontinuities obtained by this approach is more accurate and the iden-tification is more efficient.Then,a discrete fracture network with the same statistical characteristics as the actual discontinuities was generated with the distinct element method,and a representative nu-merical model of the jointed surrounding rock mass was established.By means of numerical simulation,potential unstable rock blocks were assessed,and failure mechanisms were analyzed.This method was applied to detection and assessment of unstable rock blocks in the spillway and sand flushing tunnel of the Hongshiyan hydropower project after a collapse.The results show that the noncontact detection of blocks was more labor-saving with lower safety risks compared with manual surveys,and the stability assessment was more reliable since the numerical model built by this method was more consistent with the distribution characteristics of actual joints.This study can provide a reference for geological survey and unstable rock block hazard mitigation in tunnels subjected to complex geology and active rockfalls.

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.

Optimized Hybrid Block Adams Method for Solving First Order Ordinary Differential Equations

Multistep integration methods are being extensively used in the simulations of high dimensional systems due to their lower computational cost.The block methods were developed with the intent of obtaining numerical results on numerous points at a time and improving computational efficiency.Hybrid block methods for instance are specifically used in numerical integration of initial value problems.In this paper,an optimized hybrid block Adams block method is designed for the solutions of linear and nonlinear first-order initial value problems in ordinary differential equations(ODEs).In deriving themethod,the Lagrange interpolation polynomial was employed based on some data points to replace the differential equation function and it was integrated over a specified interval.Furthermore,the convergence properties along with the region of stability of the method were examined.It was concluded that the newly derived method is convergent,consistent,and zero-stable.The method was also found to be A-stable implying that it covers the whole of the left/negative half plane.From the numerical computations of absolute errors carried out using the newly derived method,it was found that the method performed better than the ones with which we compared our results with.Themethod also showed its superiority over the existing methods in terms of stability and convergence.

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)

OPTIMIZATION AND METHOD OF DRAWING CONTROL IN BLOCK CAVING AT TONGKUANGYU MINE

ＯＰＴＩＭＩＺＡＴＩＯＮＡＮＤＭＥＴＨＯＤＯＦＤＲＡＷＩＮＧＣＯＮＴＲＯＬＩＮＢＬＯＣＫＣＡＶＩＮＧＡＴＴＯＮＧＫＵＡＮＧＹＵＭＩＮＥ①ＺｈｏｕＡｉｍｉｎＣｈａｎｇｓｈａＩｎｓｔｉｔｕｔｅｏｆＭｉｎｉｎｇＲｅｓｅａｒｃｈ，Ｃｈａｎｇｓｈａ４１００１２...

A BLOCK VARIANT OF THE GMRES METHOD FOR UNSYMMETRIC LINEAR SYSTEMS

Iterative methods that take advantage of efficient block operations and block communications are popular research topics in parallel computation. These methods are especially important on Massively Parallel Processors (MPP). This paper presents a block variant of the GMRES method for solving general unsymmetric linear systems. It is shown that the new algorithm with block size s, denoted by BVGMRES(s,m), is theoretically equivalent to the GMRES(s. m) method. The numerical results show that this algorithm can be more efficient than the standard GMRES method on a cache based single CPU computer with optimized BLAS kernels. Furthermore, the gain in efficiency is more significant on MPPs due to both efficient block operations and efficient block data communications. Our numerical results also show that in comparison to the standard GMRES method, the more PEs that are used on an MPP, the more efficient the BVGMRES(s,m) algorithm is.

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.

Chaos control of gear system with elastomeric web based on multi-parameter multi-step method

This paper employs a multi-parameter multi-step chaos control method, which is built up on the OGY method, to stabilize desirable UPOs of a gear system with elastomeric web as a high-dimensional and non-hyperbolic chaotic system, and the analyses are carried out. Three types of relations between components of a certain control parameter combination are defined in a certain control process. Special emphasis is put on the comparison of control efficiencies of the multi-parameter multi-step method and single-parameter multi-step method. The numerical experiments show the ability to switch between different orbits and the method can be a good chaos control alternative since it provides a more effective UPOs stabilization of high-dimensional and non-hyperbolic chaotic systems than the single-parameter chaos control, and according to the relation between components of each parameter combination, the best combination for chaos control in a certain UPO stabilization process are obtained.

New approaches to quantify progressive damage and associated dynamic rock mass blockiness

In the past decade, numerical modelling has been increasingly used for simulating the mechanical behaviour of naturally fractured rock masses. In this paper, we introduce new algorithms for spatial and temporal analyses of newly generated fractures and blocks using an integrated discrete fracture network (DFN)-finite-discrete element method (FDEM) (DFN-FDEM) modelling approach. A fracture line calculator and analysis technique (i.e. discrete element method (DEM) fracture analysis, DEMFA) calculates the geometrical aspects of induced fractures using a dilation criterion. The resultant two-dimensional (2D) blocks are then identified and characterised using a graph structure. Block tracking trees allow track of newly generated blocks across timesteps and to analyse progressive breakage of these blocks into smaller blocks. Fracture statistics (number and total length of initial and induced fractures) are then related to the block forming processes to investigate damage evolution. The combination of various proposed methodologies together across various stages of modelling processes provides new insights to investigate the dependency of structure's resistance on the initial fracture configuration.

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.

A Note on the Adaptive Simpler Block GMRES Method

The adaptive simpler block GMRES method was investigated by Zhong et al.(J Comput Appl Math 282:139-156, 2015) where the condition number of the adaptively chosen basis for the Krylov subspace was evaluated. In this paper, the new upper bound for the condition number is investigated. Numerical tests show that the new upper bound is tighter.

Half-Step Continuous Block Method for the Solutions of Modeled Problems of Ordinary Differential Equations

In this paper, we developed a new continuous block method using the approach of collocation of the differential system and interpolation of the power series approximate solution. A constant step length within a half step interval of integration was adopted. We evaluated at grid and off grid points to get a continuous linear multistep method. The continuous linear multistep method is solved for the independent solution to yield a continuous block method which is evaluated at selected points to yield a discrete block method. The basic properties of the block method were investigated and found to be consistent and zero stable hence convergent. The new method was tested on real life problems namely: SIR model, Growth model and Mixture Model. The results were found to compete favorably with the existing methods in terms of accuracy and error bound.