三维海洋波导中多声散射的数值解

研究三维海洋波导中多障碍物的声散射问题.基于Galerkin变分原理和DtN映射方法,发展一种自然边界元与有限元耦合的求解方法,提出有效的数值实现方法,避免常规有限元中需要采用吸收边界的缺点.对多种形状多障碍物的数值计算结果表明,所建立的理论方法是有效的.     

两个相邻目标对平面波、高斯波束的光散射

基于等效原理和互易性定理研究了两个靠近目标对平面波、高斯波束的光散射问题,给出了这一复合光散射模型的二阶散射结果。通常一阶散射结果容易求解,但由于耦合效应的复杂性,很难给出二阶散射结果的解析形式。为了解决这一问题,应用互易性定理给出了求解任意相邻介质目标二阶散射场的公式,同时借助等效原理将求解散射场公式中的体积分简化为面积分的形式,从而降低了求解难度。求解了两相邻球形粒子的复合散射场,并将求解结果与应用时域积分方程法求得的结果进行了比较。同时,还讨论了束腰半径、目标位置对散射截面及偏振度的影响。     

Noniterative analytical formula for inverse scattering of multiply scattering point targets

This paper derives, in the exact framework of multiple scattering theory for point targets, a noniterative analytical formula for the nonlinear inversion of the target scattering strengths from the scattering or response matrix that can be applied after the target positions have been estimated in a previous step via, e.g., time-reversal multiple signal classification or another approach. The new formula provides a noniterative analytical alternative to the iterative numerical solution approach for the same problem presented in a recent paper [A. J. Devaney, E. A. Marengo, and F. K. Gruber, "Time-reversal-based imaging and inverse scattering of multiply scattering point targets," J. Acoust. Soc. Am. 118, 3129-3138 (2005)]. The two methods (noniterative versus iterative) are comparatively investigated with two numerical examples.     

Iterative Method for Solving a Problem with Mixed Boundary Conditions for Biharmonic Equation

The solution of boundary value problems (BVP) for fourth order differential equations by their reduction to BVP for second order equations, with the aim to use the available efficient algorithms for the latter ones, attracts attention from many researchers. In this paper, using the technique developed by the authors in recent works we construct iterative method for a problem with complicated mixed boundary conditions for biharmonic equation which is originated from nanofluidic physics. The convergence rate of the method is proved and some numerical experiments are performed for testing its dependence on a parameter appearing in boundary conditions and on the position of the point where a transmission of boundary conditions occurs.     

等效原理和互易性定理在两个相邻球形目标电磁散射中的应用

基于等效原理和互易性定理研究了两个相邻目标的电磁散射问题，给出了这一复合电磁散射模型的二阶散射结果. 通常平面入射波的一阶散射结果容易求解，但由于耦合效应的复杂性，很难给出二阶散射结果的解析形式. 应用互易性定理给出了求解任意相邻导体/介质目标二阶电磁散射场的公式，并利用等效原理将求解散射场公式中的体积分简化为面积分的形式，从而降低了求解难度. 同时还推导了两个目标的二次散射场之间的关系. 最后应用给出的公式，求解了两相邻球形目标的复合散射场，对双站散射结果进行了讨论，同时与应用时域积分方程法求得的结果进行了比较. 关键词： 等效原理 互易性定理 电磁散射 相邻圆球     

计算应力强度因子的离散分离变量法

提出一种用于数值求解带有一条边界裂纹的多角形区域上的Navier's方程组边值问题的半离散方法。做一个适当的坐标变换后,将原边值问题化为半无限长条上的不连续系数问题。将其半离散化以后,等价于一个常系数常微分方程组的边值问题。进一步,用直接法来求解这个边值问题,便得到原问题的半离散近似解。值得指出的是,这个用分离变量形式给出的半离散近似解自然地具有原问题的奇性。数值例子显示,用该方法可以很方便地计算出在裂纹顶端的应力强度因子的近似值。     

Index to volume 72

A numerical analysis of wave propagation in one dimensional inhomogeneous media is presented in this paper, the first (Part I) of two companion papers. The aim is to provide a natural and physical basis for the solution of the inverse scattering problem in one dimension. The paper includes (i) a discussion of some aspects of the numerical solution of the direct problem, (ii) typical results of calculation and (iii) a comparison between exact solutions and forward scattering approximation (FSA) solutions of the direct problem. It is shown that the FSA allows accurate calculations of the reflected wave and thus constitutes a sound basis for the solution of the inverse problem treated in Part II.     

A Second-Order Three-Level Difference Scheme for a Magneto-Thermo-Elasticity Model

This article deals with the numerical solution to the magneto-thermo-elasticity model, which is a system of the third order partial differential equations. By introducing a new function, the model is transformed into a system of the second order generalized hyperbolic equations. A priori estimate with the conservation for the problem is established. Then a three-level finite difference scheme is derived. The unique solvability, unconditional stability and second-order convergence in $L_{\infty}$-norm of the difference scheme are proved. One numerical example is presented to demonstrate the accuracy and efficiency of the proposed method.     

Phase reduction models for improving the accuracy of the finite element solution of time-harmonic scattering problems I: General approach and low-order models

This paper introduces a new formulation of high frequency time-harmonic scattering problems in view of a numerical finite element solution. It is well-known that pollution error causes inaccuracies in the finite element solution of short-wave problems. To partially avoid this precision problem, the strategy proposed here consists in firstly numerically computing at a low cost an approximate phase of the exact solution through asymptotic propagative models. Secondly, using this approximate phase, a slowly varying unknown envelope is introduced and is computed using coarser mesh grids. The global procedure is called Phase Reduction. In this first paper, the general theoretical procedure is developed and low-order propagative models are numerically investigated in detail. Improved solutions based on higher order models are discussed showing the potential of the method for further developments.     

Acoustical scattering by arrays of cylinders in waveguides

Multiple scattering of acoustic waves in a planar horizontal waveguide by finite-length cylinders is considered. Cylinder height equals the waveguide depth, and both are vertically constrained by the pressure-release boundaries. An analytically exact solution is obtained via normal mode expansion method in conjunction with the concept of the T matrix. The problem is decomposed into an infinite number of two-dimensional multiple scattering problems, modulated by waveguide mode shapes. Examples are presented for an isovelocity waveguide in which the medium is uniform and the waveguide depth is constant. It is found that, in numerical computations, including one or two evanescent modes captures the essence of the evanescent modes. Multiple scattering in the waveguide is compared with the corresponding two-dimensional case. It is concluded that, in low frequencies, the wave patterns in the two cases are very similar, with a shift in the frequency. The similarity diminishes when there are multiple propagating modes. Despite the mode mixing, some key features in the scattering as observed in the two-dimensional problem remain observable in the waveguide case.     

奇异退化扩散反应方程的高阶紧致差分及网格自适应方法

利用余项修正法建立奇异退化扩散反应方程非均匀网格上的高阶紧致差分式,其时间具有二阶精度,空间具有三阶至四阶精度. 利用等分布原理建立时间和空间的网格自适应方法.最后通过具有精确解的数值算例验证方法的可靠性和精确性,并研究一维爆破问题.     

Regularization of boundary integral equations in the problems of wave field diffraction by curved boundaries

A regularization of the exact Fredholm integral equations for the field or its derivative on a scattering surface is proposed. This approach allows one to calculate the scattering or diffraction of pulsed wave fields by curved surfaces of arbitrary geometry. Mathematically, the method is based on the replacement of the exact Fredholm integral equations by their truncated analogs, in which the contributions of the geometrically shadowed regions are cancelled. This approach has a clear physical meaning and provides stable solutions even when the direct numerical solution of mathematically exact initial integral equations leads to unstable results. The method is mathematically substantiated and tested using the problem of plane-wave scattering by a cylinder as an example.     

On the theory of X-ray coherent reflection from a rough surface

A solution of the problem of X-ray specular reflection from a statistically rough surface is presented. It is based on using the Green function formalism. The Kirchhoff formula is used to describe the transmission of the wave field through a rough interface. Generally, microscopically exact expressions for the coefficients of transmission through a rough surface and reflection from it are obtained by taking multiple scattering effects into account. Averaging of the obtained expressions over possible realizations of random roughness of the interface between media allows to obtain rigorous expressions for specular reflection and transmission coefficients. The behavior of exact solutions in the limiting case of infinite correlation lengths is studied. It is shown that, in this case, the obtained solution corresponds to the Debye-Waller normalization. Expressions for the reflection and transmission coefficients making it possible to carry out numerical calculations are obtained in the Bourret approximation of multiple-scattering theory.     

The factorization method in a rough-surface scattering problem: I

In the present paper the wave scattering problem on rough surface is considered for the Helmholtz equation with the Dirichlet boundary condition. An approximate solution is derived with using a factorization approach to the original Helmholtz equation. As a result, the system of two equations of parabolic type appears. The first system equation has an exact analytical solution whereas for the second one, an approximate solution, is considered in terms of perturbation series. It is shown that the obtained approximate solution is the modified classical small perturbation series with respect to small Rayleigh parameter. In Appendix A it is demonstrated that, when the derived perturbation series is converged, it is possible to summarize it and to represent the exact solution of original boundary problem in an analytical symbolical form.     

The factorization method in a rough-surface scattering problem: I

In the present paper the wave scattering problem on rough surface is considered for the Helmholtz equation with the Dirichlet boundary condition. An approximate solution is derived with using a factorization approach to the original Helmholtz equation. As a result, the system of two equations of parabolic type appears. The first system equation has an exact analytical solution whereas for the second one, an approximate solution, is considered in terms of perturbation series. It is shown that the obtained approximate solution is the modified classical small perturbation series with respect to small Rayleigh parameter. In Appendix A it is demonstrated that, when the derived perturbation series is converged, it is possible to summarize it and to represent the exact solution of original boundary problem in an analytical symbolical form.     

静态电磁场边值问题计算方法

介绍常用的电磁场分析方法——解析法、有限差分法、有限元法.以静态电磁场边值问题为例,分别用自己编写的MATLAB有限差分程序及大型有限元分析软件ANSYS求出了数值解,并与解析法得到的精确解进行比较,得出了用数值法求解电磁场问题基本满足工程需要的结论.采用数值计算时,对如何减小计算机内存,取消冗余数据,也做了进一步的讨论.     

Exact boundary conditions for the initial value problem of convex conservation laws

The initial value problem of convex conservation laws, which includes the famous Burgers’ (inviscid) equation, plays an important rule not only in theoretical analysis for conservation laws, but also in numerical computations for various numerical methods. For example, the initial value problem of the Burgers’ equation is one of the most popular benchmarks in testing various numerical methods. But in all the numerical tests the initial data have to be assumed that they are either periodic or having a compact support, so that periodic boundary conditions at the periodic boundaries or two constant boundary conditions at two far apart spatial artificial boundaries can be used in practical computations. In this paper for the initial value problem with any initial data we propose exact boundary conditions at two spatial artificial boundaries, which contain a finite computational domain, by using the Lax’s exact formulas for the convex conservation laws. The well-posedness of the initial-boundary problem is discussed and the finite difference schemes applied to the artificial boundary problems are described. Numerical tests with the proposed artificial boundary conditions are carried out by using the Lax–Friedrichs monotone difference schemes.     

Characteristics of exact retrieval of the structure of plane-stratified media using reflected and transmitted pulse signals

We consider the exact iterative algorithm for solution of the direct and inverse scattering problems in 1D inhomogeneous half-space using piecewise-constant approximation of the constitutive parameter profiles (e.g., refractive index). A novel modification of the Gelfand–Levitan equation is derived for solution of the inverse scattering problem. Numerical examples are given and the effects related to a finite frequency band of the sounding pulse are discussed.     

X-ray photoelectron spectroscopy: Exact solution to the problem with internal sources

A fundamentally new model of multiple photoelectron scattering in solids is described in this paper. It takes into account the nonuniform depth distribution of photoelectron sources and the anisotropic character of single scattering and yields an exact solution to the problem of multiple elastic scattering. Highly accurate and fast algorithms for calculating X-ray photoelectron spectroscopy signals can be developed on the basis of this model.     

用近场声学测量信息可视化多个三维障碍物的一种快速算法

提出了用点声源散射场的近场信息来可视化多个三维障碍物的一种简单快速算法,它是通过利用一个线性方程组的解的三维等值线图来实现这种可视化的.这个算法的一个特别吸引人的性质是,它避开了迭代技术的应用,而且不需要关于障碍物的任何几何和物理的先验信息.此外,它只需要散射场在测量表面的某个有限范围中若干个观测点处的近场测量信息,就能得到障碍物的一个很理想的可视化.数值算例保证了这个可视化算法是有效和实用的 关键词： 声散射 反问题 物形反演