A corrected particle method with high-order Taylor expansion for solving the viscoelastic fluid flow

In this paper, a corrected particle method based on the smoothed particle hydrodynamics (SPH) method with high-order Taylor expansion (CSPH-HT) for solving the viscoelastic flow is proposed and investigated. The validity and merits of the CSPH-HT method are first tested by solving the nonlinear high order Kuramoto-Sivishinsky equation and simulating the drop stretching, respectively. Then the flow behaviors behind two stationary tangential cylinders of polymer melt, which have been received little attention, are investigated by the CSPH-HT method. Finally, the CSPH-HT method is extended to the simulation of the filling process of the viscoelastic fluid. The numerical results show that the CSPH-HT method possesses higher accuracy and stability than other corrected SPH methods and is more reliable than other corrected SPH methods.

TAYLOR EXPANSION METHOD FOR NONLINEAR EVOLUTION EQUATIONS

A new numerical method of integrating the nonlinear evolution equations, namely the Taylor expansion method, was presented. The standard Galerkin method can be viewed as the 0_th order Taylor expansion method; while the nonlinear Galerkin method can be viewed as the 1_st order modified Taylor expansion method. Moreover, the existence of the numerical solution and its convergence rate were proven. Finally, a concrete example, namely, the two_dimensional Navier_Stokes equations with a non slip boundary condition,was provided. The result is that the higher order Taylor expansion method is of the higher convergence rate under some assumptions about the regularity of the solution.

Taylor expansion MUSIC method for joint DOD and DOA estimation in a bistatic MIMO array

We propose a Taylor expansion multiple signal classification(TE MUSIC) method for joint direction of departure(DOD) and direction of arrival(DOA) estimation in a bistatic multiple-input multiple-output(MIMO)array. First, using a Taylor expansion of the steering vector, a two-dimensional(2D) search in the conventional MUSIC method for MIMO arrays is reduced to a two-step one-dimensional(1D) search in the proposed TE MUSIC method. Second, DOAs of the targets can be achieved via Lagrange multiplier by a 1D search. Finally, substituting the DOA estimates into the 2D MUSIC spectrum function, DODs of the targets are obtained by another 1D search.Thus, the DOD and DOA estimates can be automatically paired. The performance of the proposed method is better than that of the MIMO ESPRIT method, and is similar to that of the 2D MUSIC method. Furthermore, due to the 1D search, the TE MUSIC method avoids the high computational complexity of the 2D MUSIC method. Simulation results are presented to show the effectiveness of the proposed method.

An Approximate Approach for Systems of Singular Volterra Integral Equations Based on Taylor Expansion

In this article, an extended Taylor expansion method is proposed to estimate the solution of linear singular Volterra integral equations systems. The method is based on combining the m-th order Taylor polynomial of unknown functions at an arbitrary point and integration method, such that the given system of singular integral equations is converted into a system of linear equations with respect to unknown functions and their derivatives. The required solutions are obtained by solving the resulting linear system. The proposed method gives a very satisfactory solution,which can be performed by any symbolic mathematical packages such as Maple, Mathematica, etc. Our proposed approach provides a significant advantage that the m-th order approximate solutions are equal to exact solutions if the exact solutions are polynomial functions of degree less than or equal to m. We present an error analysis for the proposed method to emphasize its reliability. Six numerical examples are provided to show the accuracy and the efficiency of the suggested scheme for which the exact solutions are known in advance.

An Efficient Calculation of Photonic Crystal Band Structures Using Taylor Expansions

In this paper we present an efficient algorithm for the calculation of photonic crystal band structures and band structures of photonic crystal waveguides.Our method relies on the fact that the dispersion curves of the band structure are smooth functions of the quasi-momentum in the one-dimensional Brillouin zone.We show the derivation and computation of the group velocity,the group velocity dispersion,and any higher derivative of the dispersion curves.These derivatives are then employed in a Taylor expansion of the dispersion curves.We control the error of the Taylor expansion with the help of a residual estimate and introduce an adaptive scheme for the selection of nodes in the one-dimensional Brillouin zone at which we solve the underlying eigenvalue problem and compute the derivatives of the dispersion curves.The proposed algorithm is not only advantageous as it decreases the computational effort to compute the band structure but also because it allows for the identification of crossings and anti-crossings of dispersion curves,respectively.This identification is not possible with the standard approach of solving the underlying eigenvalue problem at a discrete set of values of the quasi-momentum without taking the mode parity into account.

Q estimation using multifrequency point average method based on the Taylor series expansion with a different order

The quality factor Q is an important parameter because it can refl ect the reservoir attenuated features and can be used for inverse-Q filtering to compensate for the seismic wave energy.The accuracy of the Q estimation is greatly significant for improving the precision of the reservoir prediction and the resolution of seismic data.In this paper,the Q estimation formulas of the single-frequency point are derived on the basis of a diff erent-order Taylor series expansion of the amplitude attenuated factor.Moreover,the multifrequency point average(MFPA)method is introduced to obtain a stable Q estimation.The model tests demonstrate that the MFPA method is less aff ected by the frequency band,travel time diff erence,time window width,and noise interference than the logical spectrum ratio(LSR)method and the energy ratio(ER)method and has a higher Q estimation accuracy.In addition,the proposed method can be applied to post-stack seismic data and obtain eff ective Q values of complex models.When the MFPA method was applied to real marine seismic data,the Q values estimated by the MFPA method with the 1st–4th order showed good consistency with each other.In contrast,the Q values obtained by the ER method were larger than those of the proposed method,while those estimated by the LSR method signifi cantly deviated from the average values.In conclusion,the MFPA method has superior stability and practicability for the Q estimation.

High-Precision Time Delay Estimation Based on Closed-Form Offset Compensation

To improve the estimation accuracy,a novel time delay estimation(TDE)method based on the closed-form offset compensation is proposed.Firstly,we use the generalized cross-correlation with phase transform(GCC-PHAT)method to obtain the initial TDE.Secondly,a signal model using normalized cross spectrum is established,and the noise subspace is extracted by eigenvalue decomposition(EVD)of covariance matrix.Using the orthogonal relation between the steering vector and the noise subspace,the first-order Taylor expansion is carried out on the steering vector reconstructed by the initial TDE.Finally,the offsets are compensated via simple least squares(LS).Compared to other state-of-the-art methods,the proposed method significantly reduces the computational complexity and achieves better estimation performance.Experiments on both simulation and real-world data verify the efficiency of the proposed approach.

A TSE based design for MMSE and QRD of MIMO systems based on ASIP

A Taylor series expansion(TSE) based design for minimum mean-square error(MMSE) and QR decomposition(QRD) of multi-input and multi-output(MIMO) systems is proposed based on application specific instruction set processor(ASIP), which uses TSE algorithm instead of resource-consuming reciprocal and reciprocal square root(RSR) operations.The aim is to give a high performance implementation for MMSE and QRD in one programmable platform simultaneously.Furthermore, instruction set architecture(ISA) and the allocation of data paths in single instruction multiple data-very long instruction word(SIMD-VLIW) architecture are provided, offering more data parallelism and instruction parallelism for different dimension matrices and operation types.Meanwhile, multiple level numerical precision can be achieved with flexible table size and expansion order in TSE ISA.The ASIP has been implemented to a 28 nm CMOS process and frequency reaches 800 MHz.Experimental results show that the proposed design provides perfect numerical precision within the fixed bit-width of the ASIP, higher matrix processing rate better than the requirements of 5G system and more rate-area efficiency comparable with ASIC implementations.

Accelerating spectral digital image correlation computation with Taylor series image reconstruction

In this paper,we introduce an accelerating algorithm based on the Taylor series for reconstructing target images in the spectral digital image correlation method(SDIC).The Taylor series image reconstruction method is employed instead of the previous direct Fourier transform(DFT)image reconstruction method,which consumes the majority of the computational time for target image reconstruction.The partial derivatives in the Taylor series are computed using the fast Fourier transform(FFT)of the entire image,following the principles of Fourier transform theory.To examine the impact of different orders of Taylor series expansion on accuracy and efficiency,we employ third-and fourth-order Taylor series image reconstruction methods and compare them with the DFT image reconstruction method through simulated experiments.As a result of these enhancements,the computational efficiency using the third-and fourth-order Taylor series improves by factors of 57 and 46,respectively,compared to the previous method.In terms of analysis accuracy,within a strain range of 0–0.1 and without the addition of image noise,the accuracy of the proposed method increases with higher expansion orders,surpassing that of the DFT image reconstruction method when the fourth order is utilized.However,when different levels of Gaussian noise are applied to simulated images individually,the accuracy of the third-or fourth-order Taylor series expansion method is superior to that of the DFT reconstruction method.Finally,we present the analyzed experimental results of a silicone rubber plate specimen with bilateral cracks under uniaxial tension.

Thin-bed thickness calculation formula and its approximation using peak frequency

Quantitative thickness estimation below tuning thickness is a great challenge in seismic exploration. Most studies focus on the thin-beds whose top and bottom reflection coefficients are of equal magnitude and opposite polarity. There is no systematic research on the other thin-bed types. In this article, all of the thin-beds are classified into four types： thin-beds with equal magnitude and opposite polarity, thin-beds with unequal magnitude and opposite polarity, thin-beds with equal magnitude and identical polarity, and thin-beds with unequal magnitude and identical polarity. By analytical study, an equation describing the general relationship between seismic peak frequency and thin-bed thickness was derived which shows there is a Complex implicit non-linear relationship between them and which is difficult to use in practice. In order to solve this problem, we simplify the relationship by Taylor expansion and discuss the precision of the approximation formulae with different orders for the four types of thin-beds. Compared with the traditional amplitude method for thin-bed thickness calculation, the method we present has a higher precision and isn＇t influenced by the absolute value of top or bottom reflection coefficient, so it is convenient for use in practice.

On band gaps of nonlocal acoustic lattice metamaterials:a robust strain gradient model

We have proposed an"exact"strain gradient(SG)continuum model to properly predict the dispersive characteristics of diatomic lattice metamaterials with local and nonlocal interactions.The key enhancement is proposing a wavelength-dependent Taylor expansion to obtain a satisfactory accuracy when the wavelength gets close to the lattice spacing.Such a wavelength-dependent Taylor expansion is applied to the displacement field of the diatomic lattice,resulting in a novel SG model.For various kinds of diatomic lattices,the dispersion diagrams given by the proposed SG model always agree well with those given by the discrete model throughout the first Brillouin zone,manifesting the robustness of the present model.Based on this SG model,we have conducted the following discussions.(Ⅰ)Both mass and stiffness ratios affect the band gap structures of diatomic lattice metamaterials,which is very helpful for the design of metamaterials.(Ⅱ)The increase in the SG order can enhance the model performance if the modified Taylor expansion is adopted.Without doing so,the higher-order continuum model can suffer from a stronger instability issue and does not necessarily have a better accuracy.The proposed SG continuum model with the eighth-order truncation is found to be enough to capture the dispersion behaviors all over the first Brillouin zone.(Ⅲ)The effects of the nonlocal interactions are analyzed.The nonlocal interactions reduce the workable range of the well-known long-wave approximation,causing more local extrema in the dispersive diagrams.The present model can serve as a satisfactory continuum theory when the wavelength gets close to the lattice spacing,i.e.,when the long-wave approximation is no longer valid.For the convenience of band gap designs,we have also provided the design space from which one can easily obtain the proper mass and stiffness ratios corresponding to a requested band gap width.

Side Wall Effects on the Hydrodynamics of a Floating Body by ImageGreen Function Based on TEBEM

A novel numerical model based on the image Green function and first-order Taylor expansion boundary element method(TEBEM), which can improve the accuracy of the hydrodynamic simulation for the non-smooth body, was developed to calculate the side wall effects on first-order motion responses and second-order drift loads upon offshore structures in the wave tank. This model was confirmed by comparing it to the results from experiments on hydrodynamic coefficients, namely the first-order motion response and second-order drift load upon a hemisphere, prolate spheroid, and box-shaped barge in the wave tank. Then,the hydrodynamics of the KVLCC2 model were also calculated in two wave tanks with different widths. It was concluded that this model can predict the hydrodynamics for offshore structures effectively, and the side wall has a significant impact on the firstorder quantities and second-order drift loads, which satisfied the resonant rule.

SOME PROPERTIES OF HOLOMORPHIC CLIFFORDIAN FUNCTIONS IN COMPLEX CLIFFORD ANALYSIS

In this article, we mainly develop the foundation of a new function theory of several complex variables with values in a complex Clifford algebra defined on some subdomains of C^n＋l, so-called complex holomorphic Cliffordian functions. We define the complex holomorphic Cliffordian functions, study polynomial and singular solutions of the equation D△^mf= 0, obtain the integral representation formula for the complex holo-morphic Cliffordian functions with values in a complex Clifford algebra defined on some submanifolds of C^n＋1, deduce the Taylor expansion and the Laurent expansion for them and prove an invariance under an action of Lie group for them.

Low complexity DOA estimation for massive UCA with single snapshot

In this paper, a low complexity direction of arrival(DOA) estimation method for massive uniform circular array(UCA) with single snapshot is proposed.Firstly, the coarse DOAs are estimated by finding the peaks from the circular convolution between a fixed coefficient vector and the received data vector.Thereafter, in order to refine coarse DOA estimates, we reconstruct the direction matrix based on the coarse DOA estimations and take the first order Taylor expansion with DOA estimation offsets into account.Finally, the refined estimations are obtained by compensating the offsets, which are obtained via least squares(LS) without any complex searches.In addition, the refinement can be iteratively implemented to enhance the estimation results.Compared to the offset search method, the proposed method achieves a better estimation performance while requiring lower complexity.Numerical simulations are presented to demonstrate the effectiveness of the proposed method.

A lattice Boltzmann model with an amending function for simulating nonlinear partial differential equations

This paper proposes a lattice Boltzmann model with an amending function for one-dimensional nonlinear partial differential equations （NPDEs） in the form ut ＋αuux ＋βu^nuz ＋γuxx ＋δuzxx ＋ζxxxx = 0. This model is different from existing models because it lets the time step be equivalent to the square of the space step and derives higher accuracy and nonlinear terms in NPDEs. With the Chapman-Enskog expansion, the governing evolution equation is recovered correctly from the continuous Boltzmann equation. The numerical results agree well with the analytical solutions.

DOA Estimation Based on Subspace Compensation for Unfolded Coprime Array

Direction of arrival(DOA)estimation for unfolded coprime array(UFCA)is discussed,and a method based on subspace compensation is proposed.Conventional DOA estimation meth-ods partition the UFCA into two subarrays for separate estimations,which are then combined for unique DOA determination.However,the DOA estimation performance loss is caused as only the partial array aperture is exploited.We use the estimations from one subarray as initial estimations,and then enhance the estimation results via a compensation based on the whole array,which is im-plemented via a simple least squares(LS)operation constructed from the initial estimation and first-order Taylor expansion.Compared to conventional methods,the DOA estimation performance is improved while the computational complexity is in the same level.Multiple simulations are con-ducted to verify the efficiency of the proposed approach.

Numerical Solution of System of Fractional Delay Differential Equations Using Polynomial Spline Functions

The aim of this paper is to approximate the solution of system of fractional delay differential equations. Our technique relies on the use of suitable spline functions of polynomial form. We introduce the description of the proposed approximation method. The error analysis and stability of the method are theoretically investigated. Numerical example is given to illustrate the applicability, accuracy and stability of the proposed method.

Higher-Order Corrections to Algebraic Derivation of Electric Dipole-Dipole Interaction

Nucleons are fermions with intrinsic spins exhibiting dipole character. Dipole-dipole interaction via their dipole moments is the key feature quantifying the short-range nucleonics interaction in two-body physics. For a pair of interacting dipoles, the energy of a pair is the quantity of interest. The same is true for chemical polar molecules. For both cases, derivation of energy almost exclusively is carried out vectorially [1]. Although uncommon the interacting energy can be derived algebraically too. For the latter Taylor, expansion is applied [2]. The given expression although appears to be correct it is incomplete. In our report, by applying Taylor’s expansion up to the 4th order and utilizing a Computer Algebra System we formulate the missing terms. Our report highlights the impact of correcting missing terms by giving two explicit examples.

LAI scale effect research based on compact airborne spectrographic imager data in the Heihe Oasis

As one of the key parameters for characterizing crop canopy structure, Leaf Area Index(LAI) has great significance in monitoring the crop growth and estimating the yield. However, due to the nonlinearity and spatial heterogeneity of LAI inversion model, there exists scale error in LAI inversion result, which limits the application of LAI product from different remote sensing data. Therefore, it is necessary to conduct studies on scale effect. This study was based on the Heihe Oasis, Zhangye city, Gansu province, China and the following works were carried out: Airborne hyperspectral CASI(Compact Airborne Spectrographic Imager) image and LAI statistic models were adopted in muti-scale LAI inversion. The overall difference of muti-scale LAI inversion was analyzed in an all-round way. This was based on two aspects, "first inversion and then integration" and "first integration and then inversion", and on scale difference characteristics of three scale transformation methods. The generation mechanism of scale effect was refined, and the optimal LAI inversion model was expanded by Taylor expansion. By doing so, it quantitatively analyzed the contribution of various inversion processes to scale effect. It was found that the cubic polynomial regression model based on NDVI(940.7 nm, 712 nm) was the optimal model, where its coefficient of determination R2 and the correlation coefficient of test samples R reached 0.72 and 0.936, respectively. Combined with Taylor expansion, it analyzed the scale error generated by LAI inversion model. After the scale effect correction of one-dimensional and twodimensional variables, the correlation coefficient of CCD-LAI(China Environment Satellite HJ/CCD images) and CASI-LAI products(Compact Airborne Spectro graphic Imager products) increased from 0.793 to 0.875 and 0.901, respectively. The mean value, standard deviation, and relative true value of the two went consistent. Compared with onedimensional variable correction method, the twodimensional method had a better correction result. This research used the effective information in hyperspectral data as sub-pixels and adopted Taylor expansion to correct the scale error in large-scale and low-resolution LAI product, achieving large-scale and high-precision LAI monitoring.

Cosmological applications in Kaluza-Klein theory

The field equations of Kaluz-Klein （KK） theory have been applied in the domain of cosmology. These equations are solved for a fiat universe by taking the gravitational and the cosmological constants as a function of time t. We use Taylor＇s expansion of cosmological function, A（t）, up to the first order of the time t. The cosmological parameters are calculated and some cosmological problems are discussed.