Wavelet Multi-Resolution Interpolation Galerkin Method for Linear Singularly Perturbed Boundary Value Problems

In this study,a wavelet multi-resolution interpolation Galerkin method(WMIGM)is proposed to solve linear singularly perturbed boundary value problems.Unlike conventional wavelet schemes,the proposed algorithm can be readily extended to special node generation techniques,such as the Shishkin node.Such a wavelet method allows a high degree of local refinement of the nodal distribution to efficiently capture localized steep gradients.All the shape functions possess the Kronecker delta property,making the imposition of boundary conditions as easy as that in the finite element method.Four numerical examples are studied to demonstrate the validity and accuracy of the proposedwavelet method.The results showthat the use ofmodified Shishkin nodes can significantly reduce numerical oscillation near the boundary layer.Compared with many other methods,the proposed method possesses satisfactory accuracy and efficiency.The theoretical and numerical results demonstrate that the order of theε-uniform convergence of this wavelet method can reach 5.

Efficient simulation of spatially correlated non-stationary ground motions by wavelet-packet algorithm and spectral representation method

Although the classical spectral representation method(SRM)has been widely used in the generation of spatially varying ground motions,there are still challenges in efficient simulation of the non-stationary stochastic vector process in practice.The first problem is the inherent limitation and inflexibility of the deterministic time/frequency modulation function.Another difficulty is the estimation of evolutionary power spectral density(EPSD)with quite a few samples.To tackle these problems,the wavelet packet transform(WPT)algorithm is utilized to build a time-varying spectrum of seed recording which describes the energy distribution in the time-frequency domain.The time-varying spectrum is proven to preserve the time and frequency marginal property as theoretical EPSD will do for the stationary process.For the simulation of spatially varying ground motions,the auto-EPSD for all locations is directly estimated using the time-varying spectrum of seed recording rather than matching predefined EPSD models.Then the constructed spectral matrix is incorporated in SRM to simulate spatially varying non-stationary ground motions using efficient Cholesky decomposition techniques.In addition to a good match with the target coherency model,two numerical examples indicate that the generated time histories retain the physical properties of the prescribed seed recording,including waveform,temporal/spectral non-stationarity,normalized energy buildup,and significant duration.

Least square method based on Haar wavelet to solve multi-dimensional stochastic Ito-Volterra integral equations

This paper proposes a method combining blue the Haar wavelet and the least square to solve the multi-dimensional stochastic Ito-Volterra integral equation.This approach is to transform stochastic integral equations into a system of algebraic equations.Meanwhile,the error analysis is proven.Finally,the effectiveness of the approach is verified by two numerical examples.

Research on Sleeve Grouting Density Detection Based on the Impact Echo Method

Grouting defects are an inherent challenge in construction practices,exerting a considerable impact on the operational structural integrity of connections.This investigation employed the impact-echo technique for the detection of grouting anomalies within connections,enhancing its precision through the integration of wavelet packet energy principles for damage identification purposes.A series of grouting completeness assessments were meticulously conducted,taking into account variables such as the divergent material properties of the sleeves and the configuration of adjacent reinforcement.The findings revealed that:(i)the energy distribution for the highstrength concrete cohort predominantly occupied the frequency bands 42,44,45,and 47,whereas for other groups,it was concentrated within the 37 to 40 frequency band;(ii)the delineation of empty sleeves was effectively discernible by examining the wavelet packet energy ratios across the spectrum of frequencies,albeit distinguishing between sleeves with 50%and full grouting density proved challenging;and(iii)the wavelet packet energy analysis yielded variable detection outcomes contingent on the material attributes of the sleeves,demonstrating heightened sensitivity when applied to ultrahigh-performance concrete matrices and GFRP-reinforced steel bars.

A Two-step Inverse Method for Quantitative Damage Evaluation of Plate-like Structures Using Vibration Approach

This study presents a novel two-step approach to assess plate-like structural laminar damages,particularly for delamination damage detection of composite structures.Firstly,a 2-D continuous wavelet transform is employed to identify the damage location and sizes from vibration curvature data.An inverse method is subsequently then used to determine the bending stiffness reduction ratio along a specified direction,enabling the quantification of the delamination severity.The method employed in this study is an extension of the one-dimensional inverse method developed in a previous work of the authors.The applicability of the two-step inverse approach is demonstrated in a simulation analysis and by an experimental study on a cantilever composite plate containing a single delamination.The inverse method is shown to have the capacity to reveal the detailed damage information of delamination within a constrained searching space and can be used to determine the effective flexural stiffness of composite plate structures,even in cases of complex delamination damage.

Integrated interpretation of dual frequency induced polarization measurement based on wavelet analysis and metal factor methods

In mineral exploration, the apparent resistivity and apparent frequency （or apparent polarizability） parameters of induced polarization method are commonly utilized to describe the induced polarization anomaly. When the target geology structure is significantly complicated, these parameters would fail to reflect the nature of the anomaly source, and wrong conclusions may be obtained. A wavelet approach and a metal factor method were used to comprehensively interpret the induced polarization anomaly of complex geologic bodies in the Adi Bladia mine. Db5 wavelet basis was used to conduct two-scale decomposition and reconstruction, which effectively suppress the noise interference of greenschist facies regional metamorphism and magma intrusion, making energy concentrated and boundary problem unobservable. On the basis of that, the ore-induced anomaly was effectively extracted by the metal factor method.

B-Spline Wavelet on Interval Finite Element Method for Static and Vibration Analysis of Stiffened Flexible Thin Plate

A new wavelet finite element method(WFEM)is constructed in this paper and two elements for bending and free vibration problems of a stiffened plate are analyzed.By means of generalized potential energy function and virtual work principle,the formulations of the bending and free vibration problems of the stiffened plate are derived separately.Then,the scaling functions of the B-spline wavelet on the interval(BSWI)are introduced to discrete the solving field variables instead of conventional polynomial interpolation.Finally,the corresponding two problems can be resolved following the traditional finite element frame.There are some advantages of the constructed elements in structural analysis.Due to the excellent features of the wavelet,such as multi-scale and localization characteristics,and the excellent numerical approximation property of the BSWI,the precise and efficient analysis can be achieved.Besides,transformation matrix is used to translate the meaningless wavelet coefficients into physical space,thus the resolving process is simplified.In order to verify the superiority of the constructed method in stiffened plate analysis,several numerical examples are given in the end.

A high-order accurate wavelet method for solving Schrdinger equations with general nonlinearity

A sampling approximation for a function defined on a bounded interval is proposed by combining the Coiflet-type wavelet expansion and the boundary extension technique. Based on such a wavelet approximation scheme, a Galerkin procedure is developed for the spatial discretization of the generalized nonlinear Schr6dinger （NLS） equa- tions, and a system of ordinary differential equations for the time dependent unknowns is obtained. Then, the classical fourth-order explicit Runge-Kutta method is used to solve this semi-discretization system. To justify the present method, several widely considered problems are solved as the test examples, and the results demonstrate that the proposed wavelet algorithm has much better accuracy and a faster convergence rate in space than many existing numerical methods.

ADAPTIVE INTERVAL WAVELET PRECISE INTEGRATION METHOD FOR PARTIAL DIFFERENTIAL EQUATIONS

The quasi-Shannon interval wavelet is constructed based on the interpolation wavelet theory, and an adaptive precise integration method, which is based on extrapolation method is presented for nonlinear ordinary differential equations ( ODEs). And then, an adaptive interval wavelet precise integration method (AIWPIM) for nonlinear partial differential equations(PDEs) is proposed. The numerical results show that the computational precision of AIWPIM is higher than that of the method constructed by combining the wavelet and the 4th Runge-Kutta method, and the computational amounts of these two methods are almost equal. For convenience, the Burgers equation is taken as an example in introducing this method, which is also valid for more general cases.

A Method for Solving Fredholm Integral Equations of the First Kind Based on Chebyshev Wavelets

In this paper, we suggest a method for solving Fredholm integral equation of the first kind based on wavelet basis. The continuous Legendre and Chebyshev wavelets of the first, second, third and fourth kind on [0,1] are used and are utilized as a basis in Galerkin method to approximate the solution of integral equations. Then, in some examples the mentioned wavelets are compared with each other.

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.

A wavelet finite-difference method for numerical simulation of wave propagation in fluid-saturated porous media

In this paper, we consider numerical simulation of wave propagation in fluidsaturated porous media. A wavelet finite-difference method is proposed to solve the 2-D elastic wave equation. The algorithm combines flexibility and computational efficiency of wavelet multi-resolution method with easy implementation of the finite-difference method. The orthogonal wavelet basis provides a natural framework, which adapt spatial grids to local wavefield properties. Numerical results show usefulness of the approach as an accurate and stable tool for simulation of wave propagation in fluid-saturated porous media.

APPLICATION OF WAVELET THEORY IN RESEARCHON WEIGHT FUNCTION OF MESHLESS METHOD

Multiresolution analysis of wavelet theory can give an effective way to describe the information at various levels of approximations or different resolutions, based on spline wavelet analysis,so weight function is orthonormally projected onto a sequence of closed spline subspaces, and is viewed at various levels of approximations or different resolutions. Now, the useful new way to research weight function is found, and the numerical result is given.

Updating method of a high-speed maglev guideway model based on wavelet transform

The structure of a high-speed maglev guideway is taken as the research object.With the aim of identifying the inconsistency of modal parameters between the simulation model and the actual model,and based on the 600 km/h high-speed maglev vehicle and the high-speed maglev test line,the arrangement of sensors and the vibration acceleration data collection of the 12.384 m concrete guideway were conducted.The modal parameters were identified from the guideway response signal using wavelet transform,after which the wavelet ridge was extracted by using the maximum slope method.Next,the vibration modes and frequency parameters of the interaction vibration characteristics of the high-speed maglev guideway and 600 km/h maglev vehicle were analyzed.The updating objective function for the finite element model of the guideway was established,and the initial guideway finite element model was modified and updated by repeatedly iterating the parameters.In doing so,the model structure of the high-speed maglev guideway was obtained,which is consistent with the actual structure.The accuracy of the updated guideway model in the calculation of the dynamic response was verified by combining this with the vehicle-guideway coupling dynamic model of the high-speed maglev system with 18 degrees of freedom.The research results reveal that the model update method based on the wavelet transform and the maximum slope method has the characteristics of high accuracy and fast recognition speed.This can effectively obtain an accurate guideway model that ensures the correctness of the vehicle-guideway coupling dynamic analysis and calculation while meeting the parameters of the measured structure model.This method is also suitable for updating other structural models of high-speed maglev systems.

General Modified Split-Step Balanced Methods for Stiff Stochastic Differential Equations

A class of general modified split-step balanced methods proposed in the paper can be applied to solve stiff stochastic differential systems with m-dimensional multiplicative noise. Compared to some other already reported split-step balanced methods, the drift increment function of the methods can be taken from any chosen ane-step ordinary differential equations （ODEs） solver. The schemes is proved to be strong convergent with order one. For the mean-square stability analysis, the investigation is confined to two cases. Some numerical experiments are reported to testify the performance and the effectiveness of the methods.

Jacobian-free Newton-Krylov subspace method with wavelet-based preconditioner for analysis of transient elastohydrodynamic lubrication problems with surface asperities

This paper presents an investigation into the effect of surface asperities on the over-rolling of bearing surfaces in transient elastohydrodynamic lubrication(EHL) line contact. The governing equations are discretized by the finite difference method. The resulting nonlinear system of algebraic equations is solved by the Jacobian-free Newtongeneralized minimal residual(GMRES) from the Krylov subspace method(KSM). The acceleration of the GMRES iteration is accomplished by a wavelet-based preconditioner.The profiles of the lubricant pressure and film thickness are obtained at each time step when the indented surface moves through the contact region. The prediction of pressure as a function of time provides an insight into the understanding of fatigue life of bearings.The analysis confirms the need for the time-dependent approach of EHL problems with surface asperities. This method requires less storage and yields an accurate solution with much coarser grids. It is stable, efficient, allows a larger time step, and covers a wide range of parameters of interest.

A novel method for gamma spectrum analysis of low-level and intermediate-level radioactive waste

The uncertainty of nuclide libraries in the analysis of the gamma spectra of low-and intermediate-level radioactive waste(LILW)using existing methods produces unstable results.To address this problem,a novel spectral analysis method is proposed in this study.In this method,overlapping peaks are located using a continuous wavelet transform.An improved quadratic convolution method is proposed to calculate the widths of the peaks and establish a fourth-order filter model to estimate the Compton edge baseline with the overlapping peaks.Combined with the adaptive sensitive nonlinear iterative peak,this method can effectively subtracts the background.Finally,a function describing the peak shape as a filter is used to deconvolve the energy spectrum to achieve accurate qualitative and quantitative analyses of the nuclide without the aid of a nuclide library.Gamma spectrum acquisition experiments for standard point sources of Cs-137 and Eu-152,a segmented gamma scanning experiment for a 200 L standard drum,and a Monte Carlo simulation experiment for triple overlapping peaks using the closest energy of three typical LILW nuclides(Sb-125,Sb-124,and Cs-134)are conducted.The results of the experiments indicate that(1)the novel method and gamma vision(GV)with an accurate nuclide library have the same spectral analysis capability,and the peak area calculation error is less than 4%;(2)compared with the GV,the analysis results of the novel method are more stable;(3)the novel method can be applied to the activity measurement of LILW,and the error of the activity reconstruction at the equivalent radius is 2.4%;and(4)The proposed novel method can quantitatively analyze all nuclides in LILW without a nuclide library.This novel method can improve the accuracy and precision of LILW measurements,provide key technical support for the reasonable disposal of LILW,and ensure the safety of humans and the environment.

Application of Wavelet Finite Element Method to Simulation of the Temperature Field of Copier Paper

Simulation of the temperature field of copier paper in copier fusing is very important for improving the fusing property of reprography. The temperature field of copier paper varies with a high gradient when the copier paper is moving through the fusing rollers. By means of conventional shaft elements, the high gradient temperature variety causes the oscillation of the numerical solution. Based on the Daubechies scaling functions, a kind of wavelet based element is constructed for the above problem. The temperature field of the copier paper moving through the fusing rollers is simulated using the two methods. Comparison of the results shows the advantages of the wavelet finite element method, which provides a new method for improving the copier properties.

A wavelet multiscale method for inversion of Maxwell equations

This paper is concerned with estimation of electrical conductivity in Maxwell equations. The primary difficulty lies in the presence of numerous local minima in the objective functional. A wavelet multiscale method is introduced and applied to the inversion of Maxwell equations. The inverse problem is decomposed into multiple scales with wavelet transform, and hence the original problem is reformulated to a set of sub-inverse problems corresponding to different scales, which can be solved successively according to the size of scale from the shortest to the longest. The stable and fast regularized Gauss-Newton method is applied to each scale. Numerical results show that the proposed method is effective, especially in terms of wide convergence, computational efficiency and precision.

The Adaptive Wavelet Collocation Method and Its Application in Front Simulation

The adaptive wavelet collocation method （AWCM） is a variable grid technology for solving partial differential equations （PDEs） with high singularities. Based on interpolating wavelets, the AWCM adapts the grid so that a higher resolution is automatically attributed to domain regions with high singularities. Accuracy problems with the AWCM have been reported in the literature, and in this paper problems of efficiency with the AWCM are discussed in detail through a simple one-dimensional （1D） nonlinear advection equation whose analytic solution is easily obtained. A simple and efficient implementation of the AWCM is investigated. Through studying the maximum errors at the moment of frontogenesis of the 1D nonlinear advection equation with different initial values and a comparison with the finite difference method （FDM） on a uniform grid, the AWCM shows good potential for modeling the front efficiently. The AWCM is also applied to a two-dimensional （2D） unbalanced frontogenesis model in its first attempt at numerical simulation of a meteorological front. Some important characteristics about the model are revealed by the new scheme.