A robust triaxial localization method of AE source using refraction path

Acoustic emission(AE)localization algorithms based on homogeneous media or single-velocity are less accurate when applied to the triaxial localization experiments.To the end,a robust triaxial localization method of AE source using refraction path is proposed.Firstly,the control equation of the refraction path is established according to the sensor coordinates and arrival times.Secondly,considering the influence of time-difference-of-arrival(TDOA)errors,the residual of the governing equation is calculated to estimate the equation weight.Thirdly,the refraction points in different directions are solved using Snell’s law and orthogonal constraints.Finally,the source coordinates are iteratively solved by weighted correction terms.The feasibility and accuracy of the proposed method are verified by pencil-lead breaking experiments.The simulation results show that the new method is almost unaffected by the refraction ratio,and always holds more stable and accurate positioning performance than the traditional method under different ratios and scales of TDOA outliers.

High-Order Decoupled and Bound Preserving Local Discontinuous Galerkin Methods for a Class of Chemotaxis Models

In this paper,we explore bound preserving and high-order accurate local discontinuous Galerkin(LDG)schemes to solve a class of chemotaxis models,including the classical Keller-Segel(KS)model and two other density-dependent problems.We use the convex splitting method,the variant energy quadratization method,and the scalar auxiliary variable method coupled with the LDG method to construct first-order temporal accurate schemes based on the gradient flow structure of the models.These semi-implicit schemes are decoupled,energy stable,and can be extended to high accuracy schemes using the semi-implicit spectral deferred correction method.Many bound preserving DG discretizations are only worked on explicit time integration methods and are difficult to get high-order accuracy.To overcome these difficulties,we use the Lagrange multipliers to enforce the implicit or semi-implicit LDG schemes to satisfy the bound constraints at each time step.This bound preserving limiter results in the Karush-Kuhn-Tucker condition,which can be solved by an efficient active set semi-smooth Newton method.Various numerical experiments illustrate the high-order accuracy and the effect of bound preserving.

Stress-corrosion coupled damage localization induced by secondary phases in bio-degradable Mg alloys:phase-field modeling

In this study,a phase-field scheme that rigorously obeys conservation laws and irreversible thermodynamics is developed for modeling stress-corrosion coupled damage(SCCD).The coupling constitutive relationships of the deformation,phase-field damage,mass transfer,and electrostatic field are derived from the entropy inequality.The SCCD localization induced by secondary phases in Mg is numerically simulated using the implicit iterative algorithm of the self-defined finite elements.The quantitative evaluation of the SCCD of a C-ring is in good agreement with the experimental results.To capture the damage localization,a micro-galvanic corrosion domain is defined,and the buffering effect on charge migration is explored.Three cases are investigated to reveal the effect of localization on corrosion acceleration and provide guidance for the design for resistance to SCCD at the crystal scale.

Local singularity and S–A methods for analyzing ore-producing anomalies in the Jianbiannongchang area of Heilongjiang,China

The Heilongjiang Jianbiannongchang area is located at the confluence of the Great and Lesser Xing’an Ranges.This area has a complex magmatic and tectonic evolutionary history that has resulted in a complex and diverse geological background for mineralization.In this study,isometric logarithmic ratio(ILR)transformations of Au,Cu,Pb,Zn,and Sb contents were performed in the1:50,000 soil geochemical data of the Jianbiannongchang area.Robust principal component analysis(RPCA)was conducted based on ILR transformation.The local singularity and spectrum-area(S-A)methods were used to extract information on mineralogic anomalies.The results showed that:(1)the transformed data eliminated the influence of the original data closure effect,and the PC1and PC2 information obtained by applying RPCA reflected ore-producing element anomalies dominated by Au and Cu.(2)The local singularity method can enhance the information of the local strong and weak slow anomalies.After performing local singularity analysis on PC1 and PC2,the obtained local anomalies reflected the local singularity spatial anomaly patterns related to Cu and Au mineralization in this area,which is an effective method for trapping ore-producing anomalies.(3)Furthermore,the composite anomaly decomposition of PC1 and PC2 was performed using the S-A method,and the screened anomalous and background fields reflect the ore-producing anomalies related to Cu and Au mineralization.This information is in agreement with known Cu and Au mineralization.(4)The geochemical anomalies with mineralization potential were obtained outside the known mineralization sites by integrating the information of oreproducing anomalies extracted by the local singularity and S-A methods,providing the theoretical basis and exploration direction for future exploration in the study area.

A Novel Localized Meshless Method for Solving Transient Heat Conduction Problems in Complicated Domains

This paper first attempts to solve the transient heat conduction problem by combining the recently proposed local knot method(LKM)with the dual reciprocity method(DRM).Firstly,the temporal derivative is discretized by a finite difference scheme,and thus the governing equation of transient heat transfer is transformed into a non-homogeneous modified Helmholtz equation.Secondly,the solution of the non-homogeneous modified Helmholtz equation is decomposed into a particular solution and a homogeneous solution.And then,the DRM and LKM are used to solve the particular solution of the non-homogeneous equation and the homogeneous solution of the modified Helmholtz equation,respectively.The LKM is a recently proposed local radial basis function collocationmethod with themerits of being simple,accurate,and free ofmesh and integration.Compared with the traditional domain-type and boundary-type schemes,the present coupling algorithm could be treated as a really good alternative for the analysis of transient heat conduction on high-dimensional and complicated domains.Numerical experiments,including two-and three-dimensional heat transfer models,demonstrated the effectiveness and accuracy of the new methodology.

A LOCAL DISCONTINUOUS GALERKIN METHOD FOR TIME-FRACTIONAL DIFFUSION EQUATIONS

In this paper,a local discontinuous Galerkin(LDG)scheme for the time-fractional diffusion equation is proposed and analyzed.The Caputo time-fractional derivative(of orderα,with 0<α<1)is approximated by a finite difference method with an accuracy of order3-α,and the space discretization is based on the LDG method.For the finite difference method,we summarize and supplement some previous work by others,and apply it to the analysis of the convergence and stability of the proposed scheme.The optimal error estimate is obtained in the L2norm,indicating that the scheme has temporal(3-α)th-order accuracy and spatial(k+1)th-order accuracy,where k denotes the highest degree of a piecewise polynomial in discontinuous finite element space.The numerical results are also provided to verify the accuracy and efficiency of the considered scheme.

Comparison of two Bayesian-point-estimation methods in multiple-source localization

Environmental uncertainty represents the limiting factor in matched-field localization. Within a Bayesian framework, both the environmental parameters, and the source parameters are considered to be unknown variables. However, including environmental parameters in multiple-source localization greatly increases the complexity and computational demands of the inverse problem. In the paper, the closed-form maximumlikelihood expressions for source strengths and noise variance at each frequency allow these parameters to be sampled implicitly, substantially reducing the dimensionality and difficulty of the inversion. This paper compares two Bayesian-point-estimation methods： the maximum a posteriori（MAP） approach and the marginal posterior probability density（PPD） approach to source localization. The MAP approach determines the sources locations by maximizing the PPD over all source and environmental parameters. The marginal PPD approach integrates the PPD over the unknowns to obtain a sequence of marginal probability distribution over source range or depth.Monte Carlo analysis of the two approaches for a test case involving both geoacoustic and water-column uncertainties indicates that：（1） For sensitive parameters such as source range, water depth and water sound speed, the MAP solution is better than the marginal PPD solution.（2） For the less sensitive parameters, such as,bottom sound speed, bottom density, bottom attenuation and water sound speed, when the SNR is low, the marginal PPD solution can better smooth the noise, which leads to better performance than the MAP solution.Since the source range and depth are sensitive parameters, the research shows that the MAP approach provides a slightly more reliable method to locate multiple sources in an unknown environment.

A FAST SEARCH METHOD OF STEERED RESPONSE POWER WITH SMALL-APERTURE MICROPHONE ARRAY FOR SOUND SOURCE LOCALIZATION

The Steered Response Power(SRP)method works well for sound source localization in noisy and reverberant environment.However,the large computation complexity limits its practical application.In this paper,a fast SRP search method is proposed to reduce the computational complexity using small-aperture microphone array.The proposed method inspired by the SRP spatial spectrum includes two steps:first,the proposed method estimates the azimuth of the sound source roughly and determines whether the sound source is in far field or near field;then,different fine searching operations are performed according to the sound source being in far field or near field.Experiments both in simulation environments and real environments have been performed to compare the localization accuracy and computation complexity of the proposed method with those of the conventional SRP-PHAT algorithm.The results show that,the proposed method has a comparative accuracy with the conventional SRP algorithm,and achieves a reduction of 93.62%in computation complexity compared to the conventional SRP algorithm.

Robust SLAM localization method based on improved variational Bayesian filtering

Aimed at the problem that the state estimation in the measurement update of the simultaneous localization and mapping(SLAM)method is incorrect or even not convergent because of the non-Gaussian measurement noise,outliers,or unknown and time-varying noise statistical characteristics,a robust SLAM method based on the improved variational Bayesian adaptive Kalman filtering(IVBAKF)is proposed.First,the measurement noise covariance is estimated using the variable Bayesian adaptive filtering algorithm.Then,the estimated covariance matrix is robustly processed through the weight function constructed in the form of a reweighted average.Finally,the system updates are iterated multiple times to further gradually correct the state estimation error.Furthermore,to observe features at different depths,a feature measurement model containing depth parameters is constructed.Experimental results show that when the measurement noise does not obey the Gaussian distribution and there are outliers in the measurement information,compared with the variational Bayesian adaptive SLAM method,the positioning accuracy of the proposed method is improved by 17.23%,20.46%,and 17.76%,which has better applicability and robustness to environmental disturbance.

Restraint on the Effect of Friction by Localization Method

The characterization of friction in control systems and its restraint by localization method are analyzed. The minimum velocity and position precision of control systems with friction are obtained analytically. The condition of global stability is given.

Integrated multi-scale approach combining global homogenization and local refinement for multi-field analysis of high-temperature superconducting composite magnets

Second-generation high-temperature superconducting(HTS)conductors,specifically rare earth-barium-copper-oxide(REBCO)coated conductor(CC)tapes,are promising candidates for high-energy and high-field superconducting applications.With respect to epoxy-impregnated REBCO composite magnets that comprise multilayer components,the thermomechanical characteristics of each component differ considerably under extremely low temperatures and strong electromagnetic fields.Traditional numerical models include homogenized orthotropic models,which simplify overall field calculation but miss detailed multi-physics aspects,and full refinement(FR)ones that are thorough but computationally demanding.Herein,we propose an extended multi-scale approach for analyzing the multi-field characteristics of an epoxy-impregnated composite magnet assembled by HTS pancake coils.This approach combines a global homogenization(GH)scheme based on the homogenized electromagnetic T-A model,a method for solving Maxwell's equations for superconducting materials based on the current vector potential T and the magnetic field vector potential A,and a homogenized orthotropic thermoelastic model to assess the electromagnetic and thermoelastic properties at the macroscopic scale.We then identify“dangerous regions”at the macroscopic scale and obtain finer details using a local refinement(LR)scheme to capture the responses of each component material in the HTS composite tapes at the mesoscopic scale.The results of the present GH-LR multi-scale approach agree well with those of the FR scheme and the experimental data in the literature,indicating that the present approach is accurate and efficient.The proposed GH-LR multi-scale approach can serve as a valuable tool for evaluating the risk of failure in large-scale HTS composite magnets.

Stability and Time-Step Constraints of Implicit-Explicit Runge-Kutta Methods for the Linearized Korteweg-de Vries Equation

This paper provides a study on the stability and time-step constraints of solving the linearized Korteweg-de Vries(KdV)equation,using implicit-explicit(IMEX)Runge-Kutta(RK)time integration methods combined with either finite difference(FD)or local discontinuous Galerkin(DG)spatial discretization.We analyze the stability of the fully discrete scheme,on a uniform mesh with periodic boundary conditions,using the Fourier method.For the linearized KdV equation,the IMEX schemes are stable under the standard Courant-Friedrichs-Lewy(CFL)conditionτ≤λh.Here,λis the CFL number,τis the time-step size,and h is the spatial mesh size.We study several IMEX schemes and characterize their CFL number as a function ofθ=d/h^(2)with d being the dispersion coefficient,which leads to several interesting observations.We also investigate the asymptotic behaviors of the CFL number for sufficiently refined meshes and derive the necessary conditions for the asymptotic stability of the IMEX-RK methods.Some numerical experiments are provided in the paper to illustrate the performance of IMEX methods under different time-step constraints.

ROBUST ACOUSTIC SOURCE LOCALIZATION FOR DIGITAL HEARING AIDS IN NOISE AND REVERBERANT ENVIRONMENT

A new method in digital hearing aids to adaptively localize the speech source in noise and reverberant environment is proposed. Based on the room reverberant model and the multichannel adaptive eigenvalue decomposition （MCAED） algorithm, the proposed method can iteratively estimate impulse response coefficients between the speech source and microphones by the adaptive subgradient projection method. Then, it acquires the time delays of microphone pairs, and calculates the source position by the geometric method. Compared with the traditional normal least mean square （NLMS） algorithm, the adaptive subgradient projection method achieves faster and more accurate convergence in a low signal-to-noise ratio （SNR） environment. Simulations for glasses digital hearing aids with four-component square array demonstrate the robust performance of the proposed method.

A local pseudo arc-length method for hyperbolic conservation laws

A local pseudo arc-length method（LPALM）for solving hyperbolic conservation laws is presented in this paper.The key idea of this method comes from the original arc-length method,through which the critical points are bypassed by transforming the computational space.The method is based on local changes of physical variables to choose the discontinuous stencil and introduce the pseudo arc-length parameter,and then transform the governing equations from physical space to arc-length space.In order to solve these equations in arc-length coordinate,it is necessary to combine the velocity of mesh points in the moving mesh method,and then convert the physical variable in arclength space back to physical space.Numerical examples have proved the effectiveness and generality of the new approach for linear equation,nonlinear equation and system of equations with discontinuous initial values.Non-oscillation solution can be obtained by adjusting the parameter and the mesh refinement number for problems containing both shock and rarefaction waves.

A MESHLESS LOCAL PETROV-GALERKIN METHOD FOR GEOMETRICALLY NONLINEAR PROBLEMS

Nonlinear formulations of the meshless local Petrov-Galerkin （MLPG） method are presented for geometrically nonlinear problems. The method requires no mesh in computation and therefore avoids mesh distortion difficulties in the large deformation analysis. The essential boundary conditions in the present formulation axe imposed by a penalty method. An incremental and iterative solution procedure is used to solve geometrically nonlinear problems. Several examples are presented to demonstrate the effectiveness of the method in geometrically nonlinear problems analysis. Numerical results show that the MLPG method is an effective one and that the values of the unknown variable are quite accurate.

ELASTIC WAVE LOCALIZATION IN TWO-DIMENSIONAL PHONONIC CRYSTALS WITH ONE-DIMENSIONAL QUASI-PERIODICITY AND RANDOM DISORDER

The band structures of both in-plane and anti-plane elastic waves propagating in two-dimensional ordered and disordered （in one direction） phononic crystals are studied in this paper. The localization of wave propagation due to random disorder is discussed by introducing the concept of the localization factor that is calculated by the plane-wave-based transfer-matrix method. By treating the quasi-periodicity as the deviation from the periodicity in a special way, two kinds of quasi phononic crystal that has quasi-periodicity （Fibonacci sequence） in one direction and translational symmetry in the other direction are considered and the band structures are characterized by using localization factors. The results show that the localization factor is an effective parameter in characterizing the band gaps of two-dimensional perfect, randomly disordered and quasi-periodic phononic crystals. Band structures of the phononic crystals can be tuned by different random disorder or changing quasi-periodic parameters. The quasi phononic crystals exhibit more band gaps with narrower width than the ordered and randomly disordered systems.

NUMERICAL ANALYSIS OF MINDLIN SHELL BY MESHLESS LOCAL PETROV-GALERKIN METHOD

The objectives of this study are to employ the meshless local Petrov-Galerkin method （MLPGM） to solve three-dimensional shell problems. The computational accuracy of MLPGM for shell problems is affected by many factors, including the dimension of compact support domain, the dimension of quadrture domain, the number of integral cells and the number of Gauss points. These factors＇ sensitivity analysis is to adopt the Taguchi experimental design technology and point out the dimension of the quadrature domain with the largest influence on the computational accuracy of the present MLPGM for shells and give out the optimum combination of these factors. A few examples are given to verify the reliability and good convergence of MLPGM for shell problems compared to the theoretical or the finite element results.

Shift control method for the local time at descendi ngnode based on sun-synchronous orbit satellite

This article analyzes the shift factors of the descending node local time for sun-synchronous satellites and proposes a shift control method to keep the local time shift within an allowance range. It is found that the satellite orbit design and the orbit injection deviation are the causes for the initial shift velocity, whereas the atmospheric drag and the sun gravitational perturbation produce the shift acceleration. To deal with these shift factors, a shift control method is put forward, through such methods as orbit variation design, orbit altitude, and inclination keeping control. The simulation experiment and practical application have proved the effectiveness of this control method.

A STUDY ON THE WEIGHT FUNCTION OF THE MOVING LEAST SQUARE APPROXIMATION IN THE LOCAL BOUNDARY INTEGRAL EQUATION METHOD

The meshless method is a new numerical technique presented in recent years.It uses the moving least square(MLS)approximation as a shape function.The smoothness of the MLS approximation is determined by that of the basic function and of the weight function,and is mainly determined by that of the weight function.Therefore,the weight function greatly affects the accuracy of results obtained.Different kinds of weight functions,such as the spline function, the Gauss function and so on,are proposed recently by many researchers.In the present work,the features of various weight functions are illustrated through solving elasto-static problems using the local boundary integral equation method.The effect of various weight functions on the accuracy, convergence and stability of results obtained is also discussed.Examples show that the weight function proposed by Zhou Weiyuan and Gauss and the quartic spline weight function are better than the others if parameters c and α in Gauss and exponential weight functions are in the range of reasonable values,respectively,and the higher the smoothness of the weight function,the better the features of the solutions.