An adaptive finite-difference method for seismic traveltime modeling based on 3D eikonal equation

3D eikonal equation is a partial differential equation for the calculation of first-arrival traveltimes and has been widely applied in many scopes such as ray tracing,source localization,reflection migration,seismic monitoring and tomographic imaging.In recent years,many advanced methods have been developed to solve the 3D eikonal equation in heterogeneous media.However,there are still challenges for the stable and accurate calculation of first-arrival traveltimes in 3D strongly inhomogeneous media.In this paper,we propose an adaptive finite-difference(AFD)method to numerically solve the 3D eikonal equation.The novel method makes full use of the advantages of different local operators characterizing different seismic wave types to calculate factors and traveltimes,and then the most accurate factor and traveltime are adaptively selected for the convergent updating based on the Fermat principle.Combined with global fast sweeping describing seismic waves propagating along eight directions in 3D media,our novel method can achieve the robust calculation of first-arrival traveltimes with high precision at grid points either near source point or far away from source point even in a velocity model with large and sharp contrasts.Several numerical examples show the good performance of the AFD method,which will be beneficial to many scientific applications.

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.

A novel adaptive harmonic balance method with an asymptotic harmonic selection

The harmonic balance method(HBM)is one of the most widely used methods in solving nonlinear vibration problems,and its accuracy and computational efficiency largely depend on the number of the harmonics selected.The adaptive harmonic balance(AHB)method is an improved HBM method.This paper presents a modified AHB method with the asymptotic harmonic selection(AHS)procedure.This new harmonic selection procedure selects harmonics from the frequency spectra of nonlinear terms instead of estimating the contribution of each harmonic to the whole nonlinear response,by which the additional calculation is avoided.A modified continuation method is proposed to deal with the variable size of nonlinear algebraic equations at different values of path parameters,and then all solution branches of the amplitude-frequency response are obtained.Numerical experiments are carried out to verify the performance of the AHB-AHS method.Five typical nonlinear dynamic equations with different types of nonlinearities and excitations are chosen as the illustrative examples.Compared with the classical HBM and Runge-Kutta methods,the proposed AHB-AHS method is of higher accuracy and better convergence.The AHB-AHS method proposed in this paper has the potential to investigate the nonlinear vibrations of complex high-dimensional nonlinear systems.

Coupled-generalized nonlinear Schr¨odinger equations solved by adaptive step-size methods in interaction picture

We extend two adaptive step-size methods for solving two-dimensional or multi-dimensional generalized nonlinear Schr ¨odinger equation(GNLSE): one is the conservation quantity error adaptive step-control method(RK4IP-CQE), and the other is the local error adaptive step-control method(RK4IP-LEM). The methods are developed in the vector form of fourthorder Runge–Kutta iterative scheme in the interaction picture by converting a vector equation in frequency domain. By simulating the supercontinuum generated from the high birefringence photonic crystal fiber, the calculation accuracies and the efficiencies of the two adaptive step-size methods are discussed. The simulation results show that the two methods have the same global average error, while RK4IP-LEM spends more time than RK4IP-CQE. The decrease of huge calculation time is due to the differences in the convergences of the relative photon number error and the approximated local error between these two adaptive step-size algorithms.

An Adaptive Sequential Replacement Method for Variable Selection in Linear Regression Analysis

With the rapid development of DNA technologies, high throughput genomic data have become a powerful leverage to locate desirable genetic loci associated with traits of importance in various crop species. However, current genetic association mapping analyses are focused on identifying individual QTLs. This study aimed to identify a set of QTLs or genetic markers, which can capture genetic variability for marker-assisted selection. Selecting a set with k loci that can maximize genetic variation out of high throughput genomic data is a challenging issue. In this study, we proposed an adaptive sequential replacement (ASR) method, which is considered a variant of the sequential replacement (SR) method. Through Monte Carlo simulation and comparing with four other selection methods: exhaustive, SR method, forward, and backward methods we found that the ASR method sustains consistent and repeatable results comparable to the exhaustive method with much reduced computational intensity.

基于Power Extrapolation和Adaptive Method的网页评估新算法

Google的PageRank算法通过对超链接结构的分析,有效地提高了搜索结果的排序质量。PowerExtrapolation算法通过特征值直接求解马尔可夫超链接矩阵的主特征向量,但该算法的迭代次数与参数d的选择密切相关,而参数d的确定目前无明显规律可寻。另一方面,AdaptiveMethod通过将马尔可夫超链接矩阵稀疏化以达到节省迭代时间的目的。文章在PowerExtrapolation算法的基础上引入AdaptiveMethod,实验结果初步证明了新算法可以减少迭代运算的时间。

An Enhanced Adaptive Differential Evolution Approach for Constrained Optimization Problems

Effective constrained optimization algorithms have been proposed for engineering problems recently.It is common to consider constraint violation and optimization algorithm as two separate parts.In this study,a pbest selection mechanism is proposed to integrate the current mutation strategy in constrained optimization problems.Based on the improved pbest selection method,an adaptive differential evolution approach is proposed,which helps the population jump out of the infeasible region.If all the individuals are infeasible,the top 5%of infeasible individuals are selected.In addition,a modified truncatedε-level method is proposed to avoid trapping in infeasible regions.The proposed adaptive differential evolution approach with an improvedεconstraint processmechanism(IεJADE)is examined on CEC 2006 and CEC 2010 constrained benchmark function series.Besides,a standard IEEE-30 bus test system is studied on the efficiency of the IεJADE.The numerical analysis verifies the IεJADE algorithm is effective in comparisonwith other effective algorithms.

Global Optimization Method Using SLE and Adaptive RBF Based on Fuzzy Clustering

High fidelity analysis models,which are beneficial to improving the design quality,have been more and more widely utilized in the modern engineering design optimization problems.However,the high fidelity analysis models are so computationally expensive that the time required in design optimization is usually unacceptable.In order to improve the efficiency of optimization involving high fidelity analysis models,the optimization efficiency can be upgraded through applying surrogates to approximate the computationally expensive models,which can greately reduce the computation time.An efficient heuristic global optimization method using adaptive radial basis function（RBF） based on fuzzy clustering（ARFC） is proposed.In this method,a novel algorithm of maximin Latin hypercube design using successive local enumeration（SLE） is employed to obtain sample points with good performance in both space-filling and projective uniformity properties,which does a great deal of good to metamodels accuracy.RBF method is adopted for constructing the metamodels,and with the increasing the number of sample points the approximation accuracy of RBF is gradually enhanced.The fuzzy c-means clustering method is applied to identify the reduced attractive regions in the original design space.The numerical benchmark examples are used for validating the performance of ARFC.The results demonstrates that for most application examples the global optima are effectively obtained and comparison with adaptive response surface method（ARSM） proves that the proposed method can intuitively capture promising design regions and can efficiently identify the global or near-global design optimum.This method improves the efficiency and global convergence of the optimization problems,and gives a new optimization strategy for engineering design optimization problems involving computationally expensive models.

COMBINED DELAUNAY TRIANGULATION AND ADAPTIVE FINITE ELEMENT METHOD FOR CRACK GROWTH ANALYSIS

The paper presents the utilization of the adaptive Delaunay triangulation in the finite element modeling of two dimensional crack propagation problems, including detailed description of the proposed procedure which consists of the Delaunay triangulation algorithm and an adaptive remeshing technique. The adaptive remeshing technique generates small elements around crack tips and large elements in the other regions. The resulting stress intensity factors and simulated crack propagation behavior are used to evaluate the effectiveness of the procedure. Three sample problems of a center cracked plate, a single edge cracked plate and a compact tension specimen, are simulated and their results assessed.

Adaptive element free Galerkin method applied to analysis of earthquake induced liquefaction

An automatically adaptive element free method is presented to analyze the seismic response of liquefiable soils. The method is based on the element free Galerkin method （EFGM） and the fission procedure that is part of h-refinement, indicated by error estimation. In the proposed method, a posteriori error estimate procedure that depends on the energy norm of stress and the T-Belytschko （TB） stress recovery scheme is incorporated. The effective cyclic elasto-plastic constitutive model is used to describe the nonlinear behavior of the saturated soil. The governing equations are established by u-p formulation. The proposed method can effectively avoid the volumetric locking due to large deformation that usually occurs in numerical computations using the finite element method （FEM）. The efficiency of the proposed method is demonstrated by evaluating the seismic response of an embankment and comparing it to results obtained through FEM. It is shown that the proposed method provides an accurate seismic analysis of saturated soil that includes the effects of liquefaction .

Adaptive Multiscale Method for Two-dimensional Nanoscale Adhesive Contacts

There are two separate traditional approaches to model contact problems: continuum and atomistic theory. Continuum theory is successfully used in many domains, but when the scale of the model comes to nanometer, continuum approximation meets challenges. Atomistic theory can catch the detailed behaviors of an individual atom by using molecular dynamics (MD) or quantum mechanics, although accurately, it is usually time-consuming. A multiscale method coupled MD and finite element (FE) is presented. To mesh the FE region automatically, an adaptive method based on the strain energy gradient is introduced to the multiscale method to constitute an adaptive multiscale method. Utilizing the proposed method, adhesive contacts between a rigid cylinder and an elastic substrate are studied, and the results are compared with full MD simulations. The process of FE meshes refinement shows that adaptive multiscale method can make FE mesh generation more flexible. Comparison of the displacements of boundary atoms in the overlap region with the results from full MD simulations indicates that adaptive multiscale method can transfer displacements effectively. Displacements of atoms and FE nodes on the center line of the multiscale model agree well with that of atoms in full MD simulations, which shows the continuity in the overlap region. Furthermore, the Von Mises stress contours and contact force distributions in the contact region are almost same as full MD simulations. The method presented combines multiscale method and adaptive technique, and can provide a more effective way to multiscale method and to the investigation on nanoscale contact problems.

Adaptive Backstepping Terminal Sliding Mode Control Method Based on Recurrent Neural Networks for Autonomous Underwater Vehicle

The trajectory tracking control problem is addressed for autonomous underwater vehicle(AUV) in marine environ?ment, with presence of the influence of the uncertain factors including ocean current disturbance, dynamic modeling uncertainty, and thrust model errors. To improve the trajectory tracking accuracy of AUV, an adaptive backstepping terminal sliding mode control based on recurrent neural networks(RNN) is proposed. Firstly, considering the inaccu?rate of thrust model of thruster, a Taylor’s polynomial is used to obtain the thrust model errors. And then, the dynamic modeling uncertainty and thrust model errors are combined into the system model uncertainty(SMU) of AUV; through the RNN, the SMU and ocean current disturbance are classified, approximated online. Finally, the weights of RNN and other control parameters are adjusted online based on the backstepping terminal sliding mode controller. In addition, a chattering?reduction method is proposed based on sigmoid function. In chattering?reduction method, the sigmoid function is used to realize the continuity of the sliding mode switching function, and the sliding mode switching gain is adjusted online based on the exponential form of the sliding mode function. Based on the Lyapu?nov theory and Barbalat’s lemma, it is theoretically proved that the AUV trajectory tracking error can quickly converge to zero in the finite time. This research proposes a trajectory tracking control method of AUV, which can e ectively achieve high?precision trajectory tracking control of AUV under the influence of the uncertain factors. The feasibility and e ectiveness of the proposed method is demonstrated with trajectory tracking simulations and pool?experi?ments of AUV.

Adaptive split-step Fourier method for simulating ultrashort laser pulse propagation in photonic crystal fibres

In this paper, the generalized nonlinear Schrodinger equation （GNLSE） is solved by an adaptive split-step Fourier method （ASSFM）. It is found that ASSFM must be used to solve GNLSE to ensure precision when the soliton selffrequency shift is remarkable and the photonic crystal fibre （PCF） parameters vary with the frequency considerably. The precision of numerical simulation by using ASSFM is higher than that by using split-step Fourier method in the process of laser pulse propagation in PCFs due to the fact that the variation of fibre parameters with the peak frequency in the pulse spectrum can be taken into account fully.

Key Techniques and Applications of Adaptive Growth Method for Stiffener Layout Design of Plates and Shells

The application of the adaptive growth method is limited because several key techniques during the design process need manual intervention of designers. Key techniques of the method including the ground structure construction and seed selection are studied, so as to make it possible to improve the effectiveness and applicability of the adaptive growth method in stiffener layout design optimization of plates and shells. Three schemes of ground structures, which are comprised by different shell elements and beam elements, are proposed. It is found that the main stiffener layouts resulted from different ground structures are almost the same, but the ground structure comprised by 8-nodes shell elements and both 3-nodes and 2-nodes beam elements can result in clearest stiffener layout, and has good adaptability and low computational cost. An automatic seed selection approach is proposed, which is based on such selection rules that the seeds should be positioned on where the structural strain energy is great for the minimum compliance problem, and satisfy the dispersancy requirement. The adaptive growth method with the suggested key techniques is integrated into an ANSYS-based program, which provides a design tool for the stiffener layout design optimization of plates and shells. Typical design examples, including plate and shell structures to achieve minimum compliance and maximum bulking stability are illustrated. In addition, as a practical mechanical structural design example, the stiffener layout of an inlet structure for a large-scale electrostatic precipitator is also demonstrated. The design results show that the adaptive growth method integrated with the suggested key techniques can effectively and flexibly deal with stiffener layout design problem for plates and shells with complex geometrical shape and loading conditions to achieve various design objectives, thus it provides a new solution method for engineering structural topology design optimization.

Adaptive energy-preserving algorithms for guiding center system

We develop two types of adaptive energy preserving algorithms based on the averaged vector field for the guiding center dynamics,which plays a key role in magnetized plasmas.The adaptive scheme is applied to the Gauss Legendre’s quadrature rules and time stepsize respectively to overcome the energy drift problem in traditional energy-preserving algorithms.These new adaptive algorithms are second order,and their algebraic order is carefully studied.Numerical results show that the global energy errors are bounded to the machine precision over long time using these adaptive algorithms without massive extra computation cost.

Adaptive H-infinity control of synchronous generators with steam valve via Hamiltonian function method

Based on Hamiltonian formulation, this paper proposes a design approach to nonlinear feedback excitation control of synchronous generators with steam valve control, disturbances and unknown parameters. It is shown that the dynamics of the synchronous generators can be expressed as a dissipative Hamiltonian system, based on which an adaptive H-infinity controller is then designed for the systems by using the structure properties of dissipative Hamiltonian systems. Simulations show that the controller obtained in this paper is very effective.

Combined indirect and direct method for adaptive fuzzy output feedback control of nonlinear system

A novel control method for a general class of nonlinear systems using fuzzy logic systems （FLSs） is presertted. Indirect and direct methods are combined to design the adaptive fuzzy output feedback controller and a high-gain observer is used to estimate the derivatives of the system output. The closed-loop system is proven to be semiglobally uniformly ultimately bounded. In addition, it is shown that if the approximation accuracy of the fuzzy logic system is high enough and the observer gain is chosen sufficiently large, an arbitrarily small tracking error can be achieved. Simulation results verify the effectiveness of the newly designed scheme and the theoretical discussion.

Adaptive hp finite element method for fluorescence molecular tomography with simplified spherical harmonics approximation

Recently,the simplified spherical harmonics equations(SP)model has at tracted much att entionin modeling the light propagation in small tissue ggeometriesat visible and near-infrared wave-leng ths.In this paper,we report an eficient numerical method for fluorescence moleeular tom-ography(FMT)that combines the advantage of SP model and adaptive hp finite elementmethod(hp-FEM).For purposes of comparison,hp-FEM and h-FEM are,respectively applied tothe reconstruction pro cess with diffusion approximation and SPs model.Simulation experiments on a 3D digital mouse atlas and physical experiments on a phantom are designed to evaluate thereconstruction methods in terms of the location and the reconstructed fluorescent yield.Theexperimental results demonstrate that hp-FEM with SPy model,yield more accurate results thanh-FEM with difusion approximation model does.The phantom experiments show the potentialand feasibility of the proposed approach in FMT applications.

h-ADAPTIVE ANALYSIS BASED ON MESHLESS LOCAL PETROV-G ALERKIN METHOD WITH B SPLINE WAVELET FOR PLATES AND SHELLS

Using the two-scale decomposition technique, the h-adaptive meshless local Petrov- Galerkin method for solving Mindlin plate and shell problems is presented. The scaling functions of B spline wavelet are employed as the basis of the moving least square method to construct the meshless interpolation function. Multi-resolution analysis is used to decompose the field variables into high and low scales and the high scale component can commonly represent the gradient of the solution according to inherent characteristics of wavelets. The high scale component in the present method can directly detect high gradient regions of the field variables. The developed adaptive refinement scheme has been applied to simulate actual examples, and the effectiveness of the present adaptive refinement scheme has been verified.

An adaptive method for high-resolution topology design

For the purpose of achieving high-resolution optimal solutions this paper proposes a nodal design variablebased adaptive method for topology optimization of continuum structures. The analysis mesh-independent density field, interpolated by the nodal design variables at a given set of density points, is adaptively refined/coarsened accord- ing to a criterion regarding the gray-scale measure of local regions. New density points are added into the gray regions and redundant ones are removed from the regions occupied by purely solid/void phases for decreasing the number of de- sign variables. A penalization factor adaptivity technique is employed-to prevent premature convergence of the optimiza- tion iterations. Such an adaptive scheme not only improves the structural boundary description quality, but also allows for sufficient further topological evolution of the structural layout in higher adaptivity levels and thus essentially enables high-resolution solutions. Moreover, compared with the case with uniformly and finely distributed density points, the proposed adaptive method can achieve a higher numerical efficiency of the optimization process.