The band-structure of a one-dimensional,periodic system in a scaling limit

The band-structure of one-dimensional Schrödinger operators is calculated when there are (a) high potential barriers (or deep wells), or (b) a wide lattice-spacing (i.e., the distance between minima of the potential). Explicit power-series formulae and error estimates are rigorously proved. The procedure used, a rigorous semi-classical method, is actually convergent for nonzero scaling parameters. Some general facts about the spectra are also discussed.     

势问题的无单元Galerkin方法的误差估计

在高维情况下,首先研究了无单元Galerkin方法的形函数构造方法——移动最小二乘法在Sobolev空间W^k,p(Ω)中的误差估计.然后,在势问题的无单元Galerkin方法的基础上,研究了势问题的通过罚函数法施加本质边界条件的无单元Galerkin方法在Sobolev空间中的误差估计.当节点和形函数满足一定条件时,证明了该误差估计是最优阶的.从误差分析中可以看出,数值解的误差与权函数的影响半径密切相关.最后,通过算例验证了结论的正确性. 关键词： 无网格方法 无单元Galerkin方法 势问题 误差估计     

Upper bound of errors in solving the inverse problem of identifying a voice source

The paper considers the inverse problem of finding the shape of a voice-source pulse from a specified segment of a speech signal using a special mathematical model that relates these quantities. A variational method for solving the formulated inverse problem for two new parametric classes of sources is proposed: a piecewise-linear source and an A-source. The error in the obtained approximate solutions of the inverse problem is considered, and a technique to numerically estimate this error is proposed, which is based on the theory of a posteriori estimates of the accuracy in solving ill-posed problems. A computer study of the adequacy of the proposed models of sources, and a study of the a posteriori estimates of the accuracy in solving inverse problems for such sources were performed using various types of voice signals. Numerical experiments for speech signals showed satisfactory properties of such a posteriori estimates, which represent the upper bounds of possible errors in solving the inverse problem. The estimate of the most probable error in determining the source-pulse shapes for the investigated speech material is on average ~7%. It is noted that the a posteriori accuracy estimates can be used as a criterion for the quality of determining the voice-source pulse shape in the speaker-identification problem.     

A new anisotropic mesh adaptation method based upon hierarchical a posteriori error estimates

A new anisotropic mesh adaptation strategy for finite element solution of elliptic differential equations is presented. It generates anisotropic adaptive meshes as quasi-uniform ones in some metric space, with the metric tensor being computed based on hierarchical a posteriori error estimates. A global hierarchical error estimate is employed in this study to obtain reliable directional information of the solution. Instead of solving the global error problem exactly, which is costly in general, we solve it iteratively using the symmetric Gauß–Seidel method. Numerical results show that a few GS iterations are sufficient for obtaining a reasonably good approximation to the error for use in anisotropic mesh adaptation. The new method is compared with several strategies using local error estimators or recovered Hessians. Numerical results are presented for a selection of test examples and a mathematical model for heat conduction in a thermal battery with large orthotropic jumps in the material coefficients.     

Kinetic Description of a Rayleigh Gas with Annihilation

In this paper, we consider the dynamics of a tagged point particle in a gas of moving hard-spheres that are non-interacting among each other. This model is known as the ideal Rayleigh gas. We add to this model the possibility of annihilation (ideal Rayleigh gas with annihilation), requiring that each obstacle is either annihilating or elastic, which determines whether the tagged particle is elastically reflected or removed from the system. We provide a rigorous derivation of a linear Boltzmann equation with annihilation from this particle model in the Boltzmann–Grad limit. Moreover, we give explicit estimates for the error in the kinetic limit by estimating the contributions of the configurations which prevent the Markovianity. The estimates show that the system can be approximated by the Boltzmann equation on an algebraically long time scale in the scaling parameter.

     

Error estimates for the fast multipole method

Error estimates for algorithms based on truncations for evaluating electrostatic interactions in molecular dynamics applications are very important for several reasons. For example, the estimates are necessary to establish the validity of the simulations and can be used to estimate various simulation parameters. Very precise estimates have been found for the Ewald method and the related particle mesh Ewald method. However, for the very popular fast multipole method such a precise estimate is not available. In this paper, we illustrate the rather complicated error behavior of the fast multipole method and we use statistical methods to derive an estimate for the root mean square error on the forces. Furthermore, the expected maximum error on the force acting on a single particle is studied. The estimates are tested against errors obtained from simulations and are found to be very precise.     

Output error estimation for summation-by-parts finite-difference schemes

The paper develops a posteriori error estimates of integral output functionals for summation-by-parts finite-difference methods. The error estimates are based on the adjoint-weighted residual method and take advantage of a variational interpretation of summation-by-parts discretizations. The estimates are computed on a fixed grid and do not require an embedded grid or explicit interpolation operators. For smooth boundary-value problems containing first and second derivatives the error estimates converge to the exact error as the mesh is refined. The theory is verified using linear boundary-value problems and the Euler equations.     

Spectral accuracy in fast Ewald-based methods for particle simulations

A spectrally accurate fast method for electrostatic calculations under periodic boundary conditions is presented. We follow the established framework of FFT-based Ewald summation, but obtain a method with an important decoupling of errors: it is shown, for the proposed method, that the error due to frequency domain truncation can be separated from the approximation error added by the fast method. This has the significance that the truncation of the underlying Ewald sum prescribes the size of the grid used in the FFT-based fast method, which clearly is the minimal grid. Both errors are of exponential-squared order, and the latter can be controlled independently of the grid size. We compare numerically to the established SPME method by Essmann et al. and see that the memory required can be reduced by orders of magnitude. We also benchmark efficiency (i.e. error as a function of computing time) against the SPME method, which indicates that our method is competitive. Analytical error estimates are proven and used to select parameters with a great degree of reliability and ease.     

Analysis of the Element-Free Galerkin Method with Penalty for Stokes Problems

The element-free Galerkin (EFG) method with penalty for Stokes problems is proposed and analyzed in this work. A priori error estimates of the penalty method, which is used to deal with Dirichlet boundary conditions, are derived to illustrate its validity in a continuous sense. Based on a feasible assumption, it is proved that there is a unique weak solution in the modified weak form of penalized Stokes problems. Then, the error bounds with the penalty factor for the EFG discretization are derived, which provide a rationale for choosing an efficient penalty factor. Numerical examples are given to confirm the theoretical results.     

一类带扰动的Sine-Gordon方程的谱方法

研究一类具扰动的Sine-Gordon方程u_tt-u_xx+αsinu-βu_xxt=g(u),t>0,-∞< x< ∞的周期初值问题,提出了谱方法,并用先验估计方法作了误差估计,证明了近似方法的收敛性,并得到了该问题广义解的存在、唯一性。     

A posteriori error estimation and adaptive mesh refinement for a multiscale operator decomposition approach to fluid–solid heat transfer

We analyze a multiscale operator decomposition finite element method for a conjugate heat transfer problem consisting of a fluid and a solid coupled through a common boundary. We derive accurate a posteriori error estimates that account for all sources of error, and in particular the transfer of error between fluid and solid domains. We use these estimates to guide adaptive mesh refinement. In addition, we provide compelling numerical evidence that the order of convergence of the operator decomposition method is limited by the accuracy of the transferred gradient information, and adapt a so-called boundary flux recovery method developed for elliptic problems in order to regain the optimal order of accuracy in an efficient manner. In an appendix, we provide an argument that explains the numerical results provided sufficient smoothness is assumed.     

Sequential reconstruction of driving-forces from nonlinear nonstationary dynamics

This paper describes a functional analysis-based method for the estimation of driving-forces from nonlinear dynamic systems. The driving-forces account for the perturbation inputs induced by the external environment or the secular variations in the internal variables of the system. The proposed algorithm is applicable to the problems for which there is too little or no prior knowledge to build a rigorous mathematical model of the unknown dynamics. We derive the estimator conditioned on the differentiability of the unknown system’s mapping, and smoothness of the driving-force. The proposed algorithm is an adaptive sequential realization of the blind prediction error method, where the basic idea is to predict the observables, and retrieve the driving-force from the prediction error. Our realization of this idea is embodied by predicting the observables one-step into the future using a bank of echo state networks (ESN) in an online fashion, and then extracting the raw estimates from the prediction error and smoothing these estimates in two adaptive filtering stages. The adaptive nature of the algorithm enables to retrieve both slowly and rapidly varying driving-forces accurately, which are illustrated by simulations. Logistic and Moran-Ricker maps are studied in controlled experiments, exemplifying chaotic state and stochastic measurement models. The algorithm is also applied to the estimation of a driving-force from another nonlinear dynamic system that is stochastic in both state and measurement equations. The results are judged by the posterior Cramer-Rao lower bounds. The method is finally put into test on a real-world application; extracting sun’s magnetic flux from the sunspot time series.     

Numerical solution of certain classes of transport equations in any dimension by Shannon sampling

A method is developed for computing solutions to some class of linear and nonlinear transport equations (hyperbolic partial differential equations with smooth solutions), in any dimension, which exploits Shannon sampling, widely used in information theory and signal processing. The method can be considered a spectral or a wavelet method, strictly related to the existence of characteristics, but allows, in addition, for some precise error estimates in the reconstruction of continuous profiles from discrete data. Non-dissipativity and (in some case) parallelizability are other features of this approach. Monotonicity-preserving cubic splines are used to handle nonuniform sampling. Several numerical examples, in dimension one or two, pertaining to single linear and nonlinear (integro-differential) equations, as well as to certain systems, are given.     

一种基于序贯估计的直达声区水面舰船 被动测距方法 *

该文利用基于射线的盲解卷积方法,从直达声区的水面舰船噪声中提取出船和锚系于海底的垂直接收阵之间的时域信道响应,并利用直达波在不同阵元相对于参考阵元的到达时间差,通过序贯方法,利用射线模型和声速剖面信息,对水面舰船距接收阵的距离进行了估计。通过处理海深约为580 m的2016年美国圣巴巴拉海峡的实验数据,对1.6~3.5 km直达声区范围内Anna Maersk商船与垂直阵之间的距离进行了估计,验证了测距方法的有效性,并将结果与系统测量值和几何方法的估计值进行了比较。由于该方法不需要对海底参数进行估计,所以在海底参数未知时要优于传统匹配场方法;在声速剖面存在跃层且海底为多层分布的复杂信道条件下,该方法仍能对距离进行有效估计,且与测量值的相对误差在6%以内,小于几何方法的估计误差,测距结果精度较高。     

Bayesian ensemble approach to error estimation of interatomic potentials

Using a Bayesian approach a general method is developed to assess error bars on predictions made by models fitted to data. The error bars are estimated from fluctuations in ensembles of models sampling the model-parameter space with a probability density set by the minimum cost. The method is applied to the development of interatomic potentials for molybdenum using various potential forms and databases based on atomic forces. The calculated error bars on elastic constants, gamma-surface energies, structural energies, and dislocation properties are shown to provide realistic estimates of the actual errors for the potentials.     

Combined phase-time encoding method in quantum cryptography

A new combined phase-time encoding method is proposed for fiber optic quantum cryptography systems. Preliminary estimates are given for the probability of the critical error to which cryptographic key distribution is possible.     

Maximum a posteriori deconvolution of ultrasonic signals using multiple transducers

A new method for deconvolution of ultrasonic pulse-echo measurements employing multiple-transducer setup is proposed in the paper. An optimal way of estimating the material reflection sequence for a linear signal generation model using maximum a posteriori estimation is proposed. The method combines the measurements from a number of transducers covering different frequency bands yielding an optimal estimate of the reflection sequence. The main idea of this approach is to complement the information unavailable from one transducer in some frequency bands with the information from the other transducers. The method is based on the assumption that the measurements are performed using transducers with identical apertures and apodization, which are located exactly at the same position relative to the test object during the measurement. An error analysis presented in the paper proves that when the above assumptions are fulfilled, the proposed method, by utilizing more data for estimation, consistently yields more accurate reflection sequence estimates than the classical Wiener filter. Experimental evidence is presented using both simulated and real ultrasonic data as a verification of the correctness of the multiple-transducer model and the estimation scheme. An illustration of the advantages of the method is also given using real ultrasonic data.     

Gaps and eigenfrequencies in the spectrum of a periodic acoustic waveguide

Sufficient conditions are found for the existence of eigenfrequencies within the spectral gaps which occur as a result of local regular perturbation of the hard walls of a periodic acoustic waveguide. Asymptotic formulas are derived and error estimates shown.     

Microwave Polarimeter with Synchronous Recording of Stokes Parameters

We describe a new method for designing a polarimeter which estimates simultaneously the four polarization parameters of the measured microwave signal, using complex phase and amplitude modulation. It is shown that, by choosing an appropriate modulation law, one can reduce significantly the dependence of the accuracy of measurement of polarization parameters on the systematic error of the meter. A block diagram of the polarimeter is presented and its main characteristics are described.     

Uncertainty quantification in DIC with Kriging regression

A Kriging regression model is developed as a post-processing technique for the treatment of measurement uncertainty in classical subset-based Digital Image Correlation (DIC). Regression is achieved by regularising the sample-point correlation matrix using a local, subset-based, assessment of the measurement error with assumed statistical normality and based on the Sum of Squared Differences (SSD) criterion. This leads to a Kriging-regression model in the form of a Gaussian process representing uncertainty on the Kriging estimate of the measured displacement field. The method is demonstrated using numerical and experimental examples. Kriging estimates of displacement fields are shown to be in excellent agreement with ‘true’ values for the numerical cases and in the experimental example uncertainty quantification is carried out using the Gaussian random process that forms part of the Kriging model. The root mean square error (RMSE) on the estimated displacements is produced and standard deviations on local strain estimates are determined.