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.

Exponential Time Differencing Method for a Reaction-Diffusion System with Free Boundary

For reaction-diffusion equations in irregular domains with moving boundaries,the numerical stability constraints from the reaction and diffusion terms often require very restricted time step sizes,while complex geometries may lead to difficulties in the accuracy when discretizing the high-order derivatives on grid points near the boundary.It is very challenging to design numerical methods that can efficiently and accurately handle both difficulties.Applying an implicit scheme may be able to remove the stability constraints on the time step,however,it usually requires solving a large global system of nonlinear equations for each time step,and the computational cost could be significant.Integration factor(IF)or exponential time differencing(ETD)methods are one of the popular methods for temporal partial differential equations(PDEs)among many other methods.In our paper,we couple ETD methods with an embedded boundary method to solve a system of reaction-diffusion equations with complex geometries.In particular,we rewrite all ETD schemes into a linear combination of specificФ-functions and apply one state-of-the-art algorithm to compute the matrix-vector multiplications,which offers significant computational advantages with adaptive Krylov subspaces.In addition,we extend this method by incorporating the level set method to solve the free boundary problem.The accuracy,stability,and efficiency of the developed method are demonstrated by numerical examples.

Integrating object-oriented methods and formal methods for requirement engineering

High quality software requirement specification is crucial for a software development. Although much efforts and research works have been done to address the problem, the errors in user requirement are still prevent us from developing high quality software. To address the problem, this paper proposes integrating graphical specification technique UML with formal specification technique to construct user requirement specification. We also present a prototype tool to perform the automatic translation from UML specification into Object-Z specification.

Integrating Geochemical and Geophysical Method in Coexistence of Oil and Potassium to Identify K-rich Brine:Research and Application in Southwestern Sichuan

Lithology of Triassic in southwestern Sichuan is consistent with the whole basin,and there is no discussion about stratum division,the difference is stratum denudation which is made by the uplifting of Luzhou uplift,especially

Multiresolution method for bending of plates with complex shapes

A high-accuracy multiresolution method is proposed to solve mechanics problems subject to complex shapes or irregular domains.To realize this method,we design a new wavelet basis function,by which we construct a fifth-order numerical scheme for the approximation of multi-dimensional functions and their multiple integrals defined in complex domains.In the solution of differential equations,various derivatives of the unknown function are denoted as new functions.Then,the integral relations between these functions are applied in terms of wavelet approximation of multiple integrals.Therefore,the original equation with derivatives of various orders can be converted to a system of algebraic equations with discrete nodal values of the highest-order derivative.During the application of the proposed method,boundary conditions can be automatically included in the integration operations,and relevant matrices can be assured to exhibit perfect sparse patterns.As examples,we consider several second-order mathematics problems defined on regular and irregular domains and the fourth-order bending problems of plates with various shapes.By comparing the solutions obtained by the proposed method with the exact solutions,the new multiresolution method is found to have a convergence rate of fifth order.The solution accuracy of this method with only a few hundreds of nodes can be much higher than that of the finite element method(FEM)with tens of thousands of elements.In addition,because the accuracy order for direct approximation of a function using the proposed basis function is also fifth order,we may conclude that the accuracy of the proposed method is almost independent of the equation order and domain complexity.

A model for extracting large deformation mining subsidence using D-InSAR technique and probability integral method

Due to the difficulties in obtaining large deformation mining subsidence using differential Interferometric Synthetic Aperture Radar （D-InSAR） alone, a new algorithm was proposed to extract large deformation mining subsidence using D-InSAR technique and probability integral method. The details of the algorithm are as follows：the control points set, containing correct phase unwrapping points on the subsidence basin edge generated by D-InSAR and several observation points （near the maximum subsidence and inflection points）, was established at first; genetic algorithm （GA） was then used to optimize the parameters of probability integral method; at last, the surface subsidence was deduced according to the optimum parameters. The results of the experiment in Huaibei mining area, China, show that the presented method can generate the correct mining subsidence basin with a few surface observations, and the relative error of maximum subsidence point is about 8.3%, which is much better than that of conventional D-InSAR （relative error is 68.0%）.

A new method for integration of a Birkhoffian system

The idea of the gradient method for integrating the dynamical equations of a nonconservative system presented by Vujanovic is transplanted to a Birkhoffian system, which is a new method for the integration of Birkhoff＇s equations. First, the differential equations of motion of the Birkhoffian system are written out. Secondly, 2n Birkhoff＇s variables are divided into two parts, and assume that a part of the variables is the functions of the remaining part of the variables and time. Thereby, the basic quasi-linear partial differential equations are established. If a complete solution of the basic partial differential equations is sought out, the solution of the problem is given by a set of algebraic equations. Since one can choose n arbitrary Birkhoff＇s variables as the functions of n remains of variables and time in a specific problem, the method has flexibility. The major difficulty of this method lies in finding a complete solution of the basic partial differential equation. Once one finds the complete solution, the motion of the systems can be obtained without doing further integration. Finally, two examples are given to illustrate the application of the results.

A SOLVING METHOD FOR A SYSTEM OF PARTIAL DIFFERENTIAL EQUATIONS WITH AN APPLICATION TO THE BENDING PROBLEM OF A THICK PLATE

A theorem of solving a system of linear non-homogeneous differential equations through integrating and adding its basic solutions is put forward and proved, the mathematical role and physical nature of the theorem is interpreted briefly. As an example, the theorem is applied to solve the problem of thermo-force bending of a thick plate.

Method of variation of parameters for solving a constrained Birkhoffian system

For an in-depth study on the integration problem of the constrained mechanical systems the method of integration for the Birkhoffian system with constraints is discussed and the method of variation of parameters for solving the dynamical equations of the constrained Birkhoffian system is provided.First the differential equations of motion for the constrained Birkhoffian system as well as for the corresponding free Birkhoffian system are established.Secondly a system of auxiliary equations is constructed and the general solution of the equations is found.Finally by varying the parameters and utilizing the properties of the generalized canonical transformation of the Birkhoffian system the solution of the problem can be obtained.The proposed method reveals the inherent relationship between the solution of a free Birkhoffian system and that of a constrained Birkhoffian system. The research results are of universal significance which can be further used in a variety of constrained mechanical systems such as non-conservative systems and nonholonomic systems etc.

Fuzzy Inventory Model under Selling Price Dependent Demand and Variable Deterioration with Fully Backlogged Shortages

The objective is to develop a model considering demand dependent on selling price and deterioration occurs after a certain period of time, which follows two-parameter Weibull distribution. Shortages are allowed and fully backlogged. Fuzzy optimal solution is obtained by considering hexagonal fuzzy numbers and for defuzzification Graded Mean Integration Representation Method. A numerical example is provided for the illustration of crisp and fuzzy, both models. To observe the effect of changes in parameters, sensitivity analysis is carried out.

3D elastic wave equation forward modeling based on the precise integration method

The Finite Difference （FD） method is an important method for seismic numerical simulations. It helps us understand regular patterns in seismic wave propagation, analyze seismic attributes, and interpret seismic data. However, because of its discretization, the FD method is only stable under certain conditions. The Arbitrary Difference Precise Integration （ADPI） method is based on the FD method and adopts an integration scheme in the time domain and an arbitrary difference scheme in the space domain. Therefore, the ADPI method is a semi-analytical method. In this paper, we deduce the formula for the ADPI method based on the 3D elastic equation and improve its stability. In forward modeling cases, the ADPI method was implemented in 2D and 3D elastic wave equation forward modeling. Results show that the travel time of the reflected seismic wave is accurate. Compared with the acoustic wave field, the elastic wave field contains more wave types, including PS- and PP- reflected waves, transmitted waves, and diffracted waves, which is important to interpretation of seismic data. The method can be easily applied to elastic wave equation numerical simulations for eoloical models.

Study on Characteristics of 3-D Translating-Pulsating Source Green Function of Deep-Water Havelock Form and Its Fast Integration Method

The singularities, oscillatory performances and the contributing factors to the 3-＇D translating-pulsating source Green function of deep-water Havelock form which consists of a local disturbance part and a far-field wave-like part, are analyzed systematically. Relative numerical integral methods about the two parts are presented in this paper. An improved method based on LOBATTO rule is used to eliminate singularities caused respectively by infinite discontinuity and jump discontinuous node from the local disturbance part function, which makes the improvement of calculation efficiency and accuracy possible. And variable substitution is applied to remove the singularity existing at the end of the integral interval of the far-field wave-like part function. Two auxiliary techniques such as valid interval calculation and local refinement of integral steps technique in narrow zones near false singularities are applied so as to avoid unnecessary integration of invalid interval and improve integral accordance. Numerical test results have proved the efficiency and accuracy in these integral methods that thus can be applied to calculate hydrodynamic performance of floating structures moving in waves.

PARAMETRIC EQUATIONS OF NONHOLONOMIC NONCONSERVATIVE SYSTEMS IN THE EVENT SPACE AND THE METHOD OF THEIR INTEGRATION

In this paper,the parametric equations with multipliers of nonholonomic nonconservative sys- tems in the event space are established,their properties are studied,and their explicit formulation is obtained. And then the field method for integrating these equations is given.Finally,an example illustrating the appli- cation of the integration method is given.

Measurement of mouse liver glutathione S-transferase activity by the integrated method

Objective: The integrated method was investigated to measure Vm/Km of mouse liver glutathione S-transfer-ase (GST) activity on GSH and 7-Cl-4-nitrobenzofurazozan. Methods: Presetting concentration of one substrate twenty-fold above the other's and taking maximum product absorbance Am as parameter while Km as constant, Vm/Km was obtained by nonlinear fitting of GST reaction curve to the integrated Michaelis-Menten equation In [Am/(Am -Ai)] + Ai/ ( ξ× Km ) = ( Vm/Km )×ti (1). Results: Vm/Km for GST showed slight dependence on initial substrate concentration and data range, but it was resistant to background absorbance, error in reaction origin and small deviation in presetting Km. Vm/Km was proportional to the amount of GST with upper limit higher than that by initial rate. There was close correlation between Vm/Km and initial rate of the same GST. Consistent results were obtained by this integrated method and classical initial rate method for the measurement of mouse liver GST. Conclusion: With the concentration of one substrate twenty-fold above the other's, this integrated method was reliable to measure the activity of enzyme on two substrates , and substrate concentration of the lower one close to its apparent Km was able to be used.

Calibration of Soil Electromagnetic Conductivity in Inverted Salinity Profiles with an Integration Method

Various calibration methods have been propounded to determine profiles of apparent bulk soil electrical conductivity (ECa) and soil electrical conductivity of a saturated soil paste extract (ECe) or a 1:5 soil water extract (EC1:5) using an electromagnetic induction instrument (EM38). The modeled coefficients, one of the successful and classical methods hitherto, were chosen to calibrate the EM38 measurements of the inverted salinity profiles of characteristic coastal saline soils at selected sites of Xincao Farm, Jiangsu Province, China. However, this method required three parameters for each depth layer. An integration approach, based on an exponential decay profile model, was proposed and the model was fitted to all the calibration sites. The obtained model can then be used to predict EC1:5 at a certain depth from electromagnetic measurements made using the EM38 device positioned in horizontal and vertical positions at the soil surface. This exponential decay model predicted the EC1:5 well according to the results of a one-way analysis of variance, and the further comparison indicated that the modeled coefficients appeared to be slightly superior to, but not statistically different from, this exponential decay model. Nevertheless, this exponential decay model was more significant and practical because it depended on less empirical parameters and could be used to perform point predictions of EC1:5 continuously with depth.

Fast precise integration method for hyperbolic heat conduction problems

A fast precise integration method is developed for the time integral of the hyperbolic heat conduction problem. The wave nature of heat transfer is used to analyze the structure of the matrix exponential, leading to the fact that the matrix exponential is sparse. The presented method employs the sparsity of the matrix exponential to improve the original precise integration method. The merits are that the proposed method is suitable for large hyperbolic heat equations and inherits the accuracy of the original version and the good computational efficiency, which are verified by two numerical examples.

Evaluation of integration methods for hybrid simulation of complex structural systems through collapse

This study examines the performance of integration methods for hybrid simulation of large and complex structural systems in the context of structural collapse due to seismic excitations. The target application is not necessarily for real-time testing, but rather for models that involve large-scale physical sub-structures and highly nonlinear numerical models. Four case studies are presented and discussed. In the first case study, the accuracy of integration schemes including two widely used methods, namely, modified version of the implicit Newmark with fixed-number of iteration （iterative） and the operator-splitting （non-iterative） is examined through pure numerical simulations. The second case study presents the results of 10 hybrid simulations repeated with the two aforementioned integration methods considering various time steps and fixed-number of iterations for the iterative integration method. The physical sub-structure in these tests consists of a single-degree-of-freedom （SDOF） cantilever column with replaceable steel coupons that provides repeatable highly- nonlinear behavior including fracture-type strength and stiffness degradations. In case study three, the implicit Newmark with fixed-number of iterations is applied for hybrid simulations of a 1：2 scale steel moment frame that includes a relatively complex nonlinear numerical substructure. Lastly, a more complex numerical substructure is considered by constructing a nonlinear computational model of a moment frame coupled to a hybrid model ofa 1：2 scale steel gravity frame. The last two case studies are conducted on the same porotype structure and the selection of time steps and fixed number of iterations are closely examined in pre-test simulations. The generated unbalance forces is used as an index to track the equilibrium error and predict the accuracy and stability of the simulations.

Modified precise time step integration method of structural dynamic analysis

The precise time step integration method proposed for linear time-invariant homogeneous dynamic systems can provide precise numerical results that approach an exact solution at the integration points. However, difficulty arises when the algorithm is used for non-homogeneous dynamic systems, due to the inverse matrix calculation and the simulation accuracy of the applied loading. By combining the Gaussian quadrature method and state space theory with the calculation technique of matrix exponential function in the precise time step integration method, a new modified precise time step integration method （e.g., an algorithm with an arbitrary order of accuracy） is proposed. In the new method, no inverse matrix calculation or simulation of the applied loading is needed, and the computing efficiency is improved. In particular, the proposed method is independent of the quality of the matrix H. If the matrix H is singular or nearly singular, the advantage of the method is remarkable. The numerical stability of the proposed algorithm is discussed and a numerical example is given to demonstrate the validity and efficiency of the algorithm.

Response Sensitivity Analysis of the Dynamic Milling Process Based on the Numerical Integration Method

As one of the bases of gradient-based optimization algorithms, sensitivity analysis is usually required to calculate the derivatives of the system response with respect to the machining parameters. The most widely used approaches for sensitivity analysis are based on time-consuming numerical methods, such as finite difference methods. This paper presents a semi-analytical method for calculation of the sensitivity of the stability boundary in milling. After transforming the delay-differential equation with time-periodic coefficients governing the dynamic milling process into the integral form, the Floquet transition matrix is constructed by using the numerical integration method. Then, the analytical expressions of derivatives of the Floquet transition matrix with respect to the machining parameters are obtained. Thereafter, the classical analytical expression of the sensitivity of matrix eigenvalues is employed to calculate the sensitivity of the stability lobe diagram. The two-degree-of-freedom milling example illustrates the accuracy and efficiency of the proposed method. Compared with the existing methods, the unique merit of the proposed method is that it can be used for analytically computing the sensitivity of the stability boundary in milling, without employing any finite difference methods. Therefore, the high accuracy and high efficiency are both achieved. The proposed method can serve as an effective tool for machining parameter optimization and uncertainty analysis in high-speed milling.

Three-dimensional numerical modeling of surface-to-borehole electromagnetic method for monitoring reservoir

Dynamic exploration for oil and gas requires careful monitoring of reservoir contents for safety and efficiency of oil extraction. This paper proposes a multi-source and multi-azimuth walk-around vertical electromagnetic profiling （MM-VEP） technique for surface-to-borehole electromagnetic surveying. Based on the difference in conductivities between reservoirs with different concentrations of oil and water, MM-VEP can be used to monitor reservoirs as they are injected with water. The MM-VEP response in five azimuth planes is modeled with three-dimensional （3D） integral equation calculations. The progress of waterflooding in four stages for enhanced oil recovery is shown to be indicated by field anomalies MM-VEP caused by variations in the reservoir resistivity. Numerical modeling demonstrates that MM-VEP measurements provides enough quantitative information from an underground reservoir to accurately detect oil deposits and monitor the progress of waterflooding.