Highly Accurate Golden Section Search Algorithms and Fictitious Time Integration Method for Solving Nonlinear Eigenvalue Problems

This study sets up two new merit functions,which are minimized for the detection of real eigenvalue and complex eigenvalue to address nonlinear eigenvalue problems.For each eigen-parameter the vector variable is solved from a nonhomogeneous linear system obtained by reducing the number of eigen-equation one less,where one of the nonzero components of the eigenvector is normalized to the unit and moves the column containing that component to the right-hand side as a nonzero input vector.1D and 2D golden section search algorithms are employed to minimize the merit functions to locate real and complex eigenvalues.Simultaneously,the real and complex eigenvectors can be computed very accurately.A simpler approach to the nonlinear eigenvalue problems is proposed,which implements a normalization condition for the uniqueness of the eigenvector into the eigenequation directly.The real eigenvalues can be computed by the fictitious time integration method(FTIM),which saves computational costs compared to the one-dimensional golden section search algorithm(1D GSSA).The simpler method is also combined with the Newton iterationmethod,which is convergent very fast.All the proposed methods are easily programmed to compute the eigenvalue and eigenvector with high accuracy and efficiency.

CAUCHY TYPE INTEGRALS AND A BOUNDARY VALUE PROBLEM IN A COMPLEX CLIFFORD ANALYSIS

Clifford analysis is an important branch of modern analysis;it has a very important theoretical significance and application value,and its conclusions can be applied to the Maxwell equation,Yang-Mill field theory,quantum mechanics and value problems.In this paper,we first give the definition of a quasi-Cauchy type integral in complex Clifford analysis,and get the Plemelj formula for it.Second,we discuss the H?lder continuity for the Cauchy-type integral operators with values in a complex Clifford algebra.Finally,we prove the existence of solutions for a class of linear boundary value problems and give the integral representation for the solution.

Element nodal computation-based radial integration BEM for non-homogeneous problems

This paper presents a new strategy of using the radial integration boundary element method （RIBEM） to solve non-homogeneous heat conduction and thermoelasticity problems. In the method, the evaluation of the radial in-tegral which is used to transform domain integrals to equivalent boundary integrals is carried out on the basis of elemental nodes. As a result, the computational time spent in evaluating domain integrals can be saved considerably in comparison with the conventional RIBEM. Three numerical examples are given to demonstrate the correctness and computational efficiency of the proposed approach.

Precise integration method for solving singular perturbation problems

This paper presents a precise method for solving singularly perturbed boundary-value problems with the boundary layer at one end. The method divides the interval evenly and gives a set of algebraic equations in a matrix form by the precise integration relationship of each segment. Substituting the boundary conditions into the algebraic equations, the coefficient matrix can be transformed to the block tridiagonal matrix. Considering the nature of the problem, an efficient reduction method is given for solving singular perturbation problems. Since the precise integration relationship introduces no discrete error in the discrete process, the present method has high precision. Numerical examples show the validity of the present method.

A new simple method of implicit time integration for dynamic problems of engineering structures

This paper presents a new simple method of implicit time integration with two control parameters for solving initial-value problems of dynamics such that its accuracy is at least of order two along with the conditional and unconditional stability regions of the parameters. When the control parameters in the method are optimally taken in their regions, the accuracy may be improved to reach of order three. It is found that the new scheme can achieve lower numerical amplitude dissipation and period dispersion than some of the existing methods, e.g. the Newmark method and Zhai＇s approach, when the same time step size is used. The region of time step dependent on the parameters in the new scheme is explicitly obtained. Finally, some examples of dynamic problems are given to show the accuracy and efficiency of the proposed scheme applied in dynamic systems.

Precise integration method for a class of singular two-point boundary value problems

In this paper we present a precise integration method based on high order multiple perturbation method and reduction method for solving a class of singular twopoint boundary value problems.Firstly,by employing the method of variable coefficient dimensional expanding,the non-homogeneous ordinary differential equations（ODEs） are transformed into homogeneous ODEs.Then the interval is divided evenly,and the transfer matrix in every subinterval is worked out using the high order multiple perturbation method,and a set of algebraic equations is given in the form of matrix by the precise integration relation for each segment,which is worked out by the reduction method.Finally numerical examples are elaboratedd to validate the present method.

Consistency and Stability Issues in the Numerical Integration of the First and Second Order Initial Value Problem

In this note we consider some basic, yet unusual, issues pertaining to the accuracy and stability of numerical integration methods to follow the solution of first order and second order initial value problems (IVP). Included are remarks on multiple solutions, multi-step methods, effect of initial value perturbations, as well as slowing and advancing the computed motion in second order problems.

Integration of axiomatic design and theory of inventive problem solving for conceptual design

Axiomatic design (AD) and theory of inventive problem solving (TRIZ) are widely used in conceptual design. Both of them have limitations, however. We presented an integrated model of these two methods to increase the efficiency and quality of the problem solving process for conceptual design. AD is used for systematically defining and structuring a problem into a hierarchy. Sometimes, the design matrix is coupled in AD which indicates the functional requirements are coupled. TRIZ separation principles can be used to separate non-independent design parameters, which provide innovative solutions at each hierarchical level. We applied the integrated model to the heating and drying equipment of bitumen reproduction device. The result verifies that the integrated model can work very well in conceptual design.

Volterra Integral Equations and Some Nonlinear Integral Equations with Variable Limit of Integration as Generalized Moment Problems

In this paper we will see that, under certain conditions, the techniques of generalized moment problem will apply to numerically solve an Volterra integral equation of first kind or second kind. Volterra integral equation is transformed into a one-dimensional generalized moment problem, and shall apply the moment problem techniques to find a numerical approximation of the solution. Specifically you will see that solving the Volterra integral equation of first kind f（t） = {a^t K（t, s）x（s）ds a ≤ t ≤ b or solve the Volterra integral equation of the second kind x（t） =f（t）＋{a^t K（t,s）x（s）ds a ≤ t ≤ b is equivalent to solving a generalized moment problem of the form un = {a^b gn（s）x（s）ds n = 0,1,2… This shall apply for to find the solution of an integrodifferential equation of the form x＇（t） = f（t） ＋ {a^t K（t,s）x（s）ds for a ≤ t ≤ b and x（a） = a0 Also considering the nonlinear integral equation： f（x）= {fa^x y（x-t）y（t）dt This integral equation is transformed a two-dimensional generalized moment problem. In all cases, we will find an approximated solution and bounds for the error of the estimated solution using the techniques ofgeneralized moment problem.

An Eight Component Integrable Hamiltonian Hierarchy from a Reduced Seventh-Order Matrix Spectral Problem

We present an eight component integrable Hamiltonian hierarchy, based on a reduced seventh order matrix spectral problem, with the aim of aiding the study and classification of multicomponent integrable models and their underlying mathematical structures. The zero-curvature formulation is the tool to construct a recursion operator from the spatial matrix problem. The second and third set of integrable equations present integrable nonlinear Schrödinger and modified Korteweg-de Vries type equations, respectively. The trace identity is used to construct Hamiltonian structures, and the first three Hamiltonian functionals so generated are computed.

Solving Large Scale Unconstrained Minimization Problems by a New ODE Numerical Integration Method

In reference [1], for large scale nonlinear equations , a new ODE solving method was given. This paper is a continuous work. Here has gradient structure i.e. , is a scalar function. The eigenvalues of the Jacobian of;or the Hessian of , are all real number. So the new method is very suitable for this structure. For quadratic function the convergence was proved and the spectral radius of iteration matrix was given and compared with traditional method. Examples show for large scale problems (dimension ) the new method is very efficient.

FRACTIONAL ORDER BOUNDARY VALUE PROBLEM WITH INTEGRAL BOUNDARY CONDITIONS INVOLVING PETTIS INTEGRAL

In this article, we investigate the existence of Pseudo solutions for some fractional order boundary value problem with integral boundary conditions in the Banach space of continuous function equipped with its weak topology. The class of such problems constitute a very interesting and important class of problems. They include two, three, multi-point and nonlocal boundary-value problems as special cases. In our investigation, the right hand side of the above problem is assumed to be Pettis integrable function. To encompass the full scope of this article, we give an example illustrating the main result.

Location and allocation problem for spare parts depots on integrated logistics support

In equipment integrated logistics support(ILS), the supply capability of spare parts is a significant factor. There are lots of depots in the traditional support system, which makes too many redundant spare parts and causes high cost of support. Meanwhile,the inconsistency among depots makes it difficult to manage spare parts. With the development of information technology and transportation, the supply network has become more efficient. In order to further improve the efficiency of supply-support work and the availability of the equipment system, building a system of one centralized depot with multiple depots becomes an appropriate way.In this case, location selection of the depots including centralized depots and multiple depots becomes a top priority in the support system. This paper will focus on the location selection problem of centralized depots considering ILS factors. Unlike the common location selection problem, depots in ILS require a higher service level. Therefore, it becomes desperately necessary to take the high requirement of the mission into account while determining location of depots. Based on this, we raise an optimal depot location model. First, the expected transportation cost is calculated.Next, factors in ILS such as response time, availability and fill rate are analyzed for evaluating positions of open depots. Then, an optimization model of depot location is developed with the minimum expected cost of transportation as objective and ILS factors as constraints. Finally, a numerical case is studied to prove the validity of the model by using the genetic algorithm. Results show that depot location obtained by this model can guarantee the effectiveness and capability of ILS well.

NEW NUMERICAL METHOD FOR VOLTERRA INTEGRAL EQUATION OF THE SECOND KIND IN PIEZOELASTIC DYNAMIC PROBLEMS

The elastodynamic problems of piezoelectric hollow cylinders and spheres under radial deformation can be transformed into a second kind Volterra integral equation about a function with respect to time, which greatly simplifies the solving procedure for such elastodynamic problems. Meanwhile, it becomes very important to find a way to solve the second kind Volterra integral equation effectively and quickly. By using an interpolation function to approximate the unknown function, two new recursive formulae were derived, based on which numerical solution can be obtained step by step. The present method can provide accurate numerical results efficiently. It is also very stable for long time calculating.

Second Order Periodic Boundary Value Problems Involving the Distributional Henstock-Kurzweil Integral

We apply the distributional derivative to study the existence of solutions of the second order periodic boundary value problems involving the distributional Henstock-Kurzweil integral. The distributional Henstock-Kurzweil integral is a general intergral, which contains the Lebesgue and Henstock-Kurzweil integrals. And the distributional derivative includes ordinary derivatives and approximate derivatives. By using the method of upper and lower solutions and a fixed point theorem, we achieve some results which are the generalizations of some previous results in the literatures.

ON THE UNIQUENESS OF BOUNDARY INTEGRAL EQUATION FOR THE EXTERIOR HELMHOLTZ PROBLEM

From the point of view of energy analysis, the cause that the uniqueness of the boundary integral equation induced from the exterior Helmholtz problem does not hold is investigated in this paper. It is proved that the Sommerfeld's condition at the infinity is changed so that it is suitable not only for the radiative wave but also for the absorptive wave when we use the boundary integral equation to describe the exterior Helmholtz problem. There fore, the total energy of the system is conservative. The mathematical dealings to guarantee the uniqueness are discussed based upon this explanation.

NOVEL REGULARIZED BOUNDARY INTEGRAL EQUATIONS FOR POTENTIAL PLANE PROBLEMS

The universal practices have been centralizing on the research of regularization to the direct boundary integal equations （DBIEs）. The character is elimination of singularities by using the simple solutions. However, up to now the research of regularization to the first kind integral equations for plane potential problems has never been found in previous literatures. The presentation is mainly devoted to the research on the regularization of the singular boundary integral equations with indirect unknowns. A novel view and idea is presented herein, in which the regularized boundary integral equations with indirect unknowns without including the Cauchy principal value （CPV） and Hadamard-finite-part （HFP） integrals are established for the plane potential problems. With some numerical results, it is shown that the better accuracy and higher efficiency, especially on the boundary, can be achieved by the present system.

Integral Global Minimization of Constrained Problems with Discontinuous Penalty Functions

A class of discontinuous penalty functions was proposed to solve constrained minimization problems with the integral approach to global optimization, m-mean value and v-variance optimality conditions of a constrained and penalized minimization problem were investigated. A nonsequential algorithm was proposed. Numerical examples were given to illustrate the effectiveness of the algorithm.

Numerical Solution of Mean-Square Approximation Problem of Real Nonnegative Function by the Modulus of Double Fourier Integral

A nonlinear problem of mean-square approximation of a real nonnegative continuous function with respect to two variables by the modulus of double Fourier integral dependent on two real parameters with use of the smoothing functional is studied. Finding the optimal solutions of this problem is reduced to solution of the Hammerstein type two-dimensional nonlinear integral equation. The numerical algorithms to find the branching lines and branching-off solutions of this equation are constructed and justified. Numerical examples are presented.

THE EXACT INTEGRAL EQUATION OF HERTZ’S CONTACT PROBLEM

This paper presents the exact integral equation of Hertz's contact problem, which is obtained by taking into account the horizontal displacement of points in the contacted surfaces due to pressure.