水平层状介质中电磁场并矢Green函数的一种快速新算法

本文利用提取直射波并结合自适应数字滤波等技术提出一种计算水平层状介质中电磁场并矢Green函数的快速算法. 首先将谱域Green函数中表征均匀介质作用的直射波提取出来并对其积分进行解析计算,这种处理降低了谱域Green函数的奇异性,可在很大程度上缩短其积分收敛区间. 然后在将谱域Green函数剩余部分对应积分转化为三个快速下降积分的基础上,引入一种自适应数字滤波算法对其进行快速求解. 最后通过具体算例验证了本文所述算法的有效性. 关键词： 并矢Green函数 快速算法 水平层状介质     

各向异性介质三维电磁响应模拟的Ho-GEBA算法

本文基于积分方程法研究并建立了一种模拟横向同性介质中任意各向异性异常 体三维电磁响应的高阶广义扩展Born近似(Ho-GEBA)算法. 首先利用逐次迭代技术给出积分方程的广义级数展开解, 为保证其收敛性, 引入一种各向异性条件下满足压缩映射的迭代算子. 然后利用异常体区域分解技术, 并结合扩展Born近似原理, 得到各向异性介质三维电磁响应的Ho-GEBA解. 为提高效率, 计算过程中采用并矢Green函数的解析表达式. 最后通过数值计算实例对比验证了本文算法的有效性. 关键词： 高阶广义扩展Born近似 积分方程 电磁模拟 解析Green函数     

用改进的传输线算法计算水平层状横向同性地层中海洋可控源电磁响应

利用二维Fourier变换与电磁场分解技术将层状横向同性地层中Maxwell方程转化成两个独立的关于横磁(TM)波和横电(TE)波的传输线方程; 借助传输线理论与叠加原理, 仅利用电流源传输线Green函数得到TM波和TE波的解, 改进传输线算法, 建立横向同性地层中频率-波数域电流源电场和磁场并矢Green函数的新算法与新的解析表达式, 提高海洋可控源电磁响应数值模拟效率. 在此基础上, 利用传输线Green函数的基本解以及边界条件, 推导出广义反射系数与振幅递推公式, 得到各个地层中传输线Green函数的解析解; 然后利用Fourier逆变换与Bessel公式将海洋可控源电磁响应表示为Sommerfeld形式的积分, 借助三次样条插值与Lommel积分公式快速计算其数值解. 通过数值模拟结果考察工作频率以及地层各向异性电阻率变化等对海洋电磁响应的影响. 关键词： 传输线法 横向同性地层 海洋可控源电磁 Sommerfeld积分     

水平层状介质并矢Green函数的递推矩阵方法

采用递推矩阵方法计算任意数目水平层状介质的并矢Green函数.根据层界面处电场和磁场的连续性条件得到3个确定Sommerfeld积分待定系数的矩阵方程组,分别对应于垂向单位电偶极子产生的TM波、水平方向单位电偶极子产生的TE波和TM波,这些方程组均可通过递推方法快速求解.只需改变3个方程组中源项元素的位置,就可以方便地得到当源点和场点在任意层时的并矢Green函数.本文给出的并矢Green函数表达式形式简洁且不含指数增加项,计算时不会出现溢出现象.     

水平层状横向同性地层中频率测深资料全参数快速迭代反演

提出水平层状横向同性地层中频率测深资料的全参数快速迭代反演算法,以便从测量资料中同时确定各个地层的横向、 纵向电导率及层界面深度.首先,利用水平层状介质中电磁场并矢Green函数在频率波数域中的解析解和Sommerfeld积分的 快速计算技术确定正演响应.然后,利用摄动理论和Fourier逆变换公式,研究建立一套快速求解全参数Fréchet导数的 有效算法,并利用规范化处理和奇异值分解技术给出迭代反演过程.最后,给出数值试验证明反演理论的有效性以及反演算法 的抗噪声能力.     

广义Helmholtz方程的Green函数及旋波介质色散关系的另一求解方法

运用并矢法求解广义Helmholtz方程,结合Fourier变换、留数定理等数学处理,给出方程的并矢Green函数解析表达式;并进一步就旋波介质给出其中色散关系的另外一种求解方法．     

利用积分方程法的各向异性地层频率测深三维模拟

应用一种新的积分方程算法对各向异性地层频率域电磁测深三维问题进行正演模拟.利用范数小于或等于1的修正Green算子得到各向异性地层中的新积分方程,由于满足压缩映射条件,该积分方程在任意参数条件下总是迭代收敛的.提出一种应用迭代法求解新积分方程的迭代初值优化选择方法,并通过具体算例对比说明该方法可使频率域电磁测深三维模拟的计算效率得到有效提高.最后应用所述算法对层状各向异性大地轴向频率测深视电阻率响应进行三维数值模拟,考察分析地层的各向异性对视电阻率响应特征的影响,得到一些重要结论.     

任意半径圆形电流在径向成层介质中产生的电磁场

采用递推矩阵方法计算径向成层介质中的二维Green函数.根据层界面处电场和磁场的连续性条件得到确定待定系数的矩阵方程组并通过递推方法快速求解.只需改变方程组中源项元素的位置,就可以方便地得到当源点和场点在任意层时的二维Green函数,并进而得到具有任意半径的圆形电流在介质中产生的电磁场.本文给出的Green函数具有表达方式简洁的优点.     

半无界空间上函数的傅里叶-贝塞尔积分展开

柱坐标系中,本征函数族贝塞尔函数构成完备正交系,因此可作为广义傅里叶级数展开的基.本文从定义在有限区间[0,ρ0]上函数的广义傅里叶级数展开出发,利用贝塞尔函数的渐近展开公式以及贝塞尔函数零点的近似公式,讨论了半无界空间上函数的傅里叶-贝塞尔积分展开问题,得到了本征函数模方的近似表达式.当ρ0趋于无穷时,不连续参量变成连续参量,得到了函数的傅里叶-贝塞尔积分及其展开系数公式.     

跨界面目标电磁散射Sommerfeld积分的双重广义函数束拟合离散复镜像方法

混合势积分方程（mixed-potential electric field integral equation, MPIE）由于其Green函数的低阶奇异性是矩量法（MoM）等数值方法求解散射问题的主要方程.在跨界面目标的散射问题中,MPIE中的矢量势和标量势Green函数包含Sommerfeld类型的谱域积分,利用离散复镜像法（discrete complex image method, DCIM）和Sommerfeld恒等式,将其转化为有限项复镜像Green函数的求和运算,避免了烦杂的谱域积分运算. 关键词： Sommerfeld积分 离散复镜像法 广义函数束     

Scattering and diffraction of a plane wave by a randomly rough half-plane

The scattering and diffraction of a TE (transverse electric) plane wave by a randomly rough half-plane are studied by a combination of three techniques: the Wiener-Hopf technique, the small perturbation method and a probabilistic method based on the shift-invariance of a homogeneous random function. By use of the D^a-Fourier transformation based on the shift-invariance, it is shown that the scattered wave is written by an inverse Fourier transformation of a homogeneous random function with a complex parameter. For a small rough case, such a random function with a complex parameter is expanded in a perturbation series and then the first-order solution is obtained exactly in an integral form. The first-order solution involves two physical processes such that the edge-diffracted wave is scattered by the randomly rough plane and the scattered wave, due to roughness, is diffracted by the half-plane. The solution is transformed into a sum of the Fresnel integrals with complex arguments, an integral along the steepest descent path and a branch-cut integral, which are evaluated numerically. Then, intensities of the coherently scattered wave and incoherent wave are calculated in the region near the edge and illustrated in figures.     

Path integral for spinning particle in the plane wave field: Global and local projections

The Green function related to the problem of a Dirac particle interacting with a plane wave is calculated via the path integral formalism proposed recently by Alexandrou et al. according to the two so-called global and local projections. With the help of the incorporation of two simple identities, it is shown that the contribution to the calculation of the integrals comes essentially from classical solutions projected along the direction of wave propagation.     

有限深水域中三维频域自由面格林函数及其导数的计算方法

三维频域自由面格林函数及其偏导数的数值计算是海洋结构物水动力分析的难点。本文对有限水深格林函数及其偏导数的数值算法作出两点改进：①对原函数的级数表达式进行变形,提出适用于数值计算的形式,解决高频率情况下级数解的计算失真问题;②对Gauss-Laguerre积分法作出改进,解决高频率大水深情况下积分解的数值溢出问题。计算结果表明：本文提出的改进算法可以有效提高有限水深格林函数及其偏导数的数值精度,对于研究较大吃水的浮体或潜体的高频垂荡问题具有重要意义。     

Scattering and diffraction of a plane wave by a randomly rough half-plane*

Abstract

The scattering and diffraction of a TE (transverse electric) plane wave by a randomly rough half-plane are studied by a combination of three techniques: the Wiener-Hopf technique, the small perturbation method and a probabilistic method based on the shift-invariance of a homogeneous random function. By use of the D^a-Fourier transformation based on the shift-invariance, it is shown that the scattered wave is written by an inverse Fourier transformation of a homogeneous random function with a complex parameter. For a small rough case, such a random function with a complex parameter is expanded in a perturbation series and then the first-order solution is obtained exactly in an integral form. The first-order solution involves two physical processes such that the edge-diffracted wave is scattered by the randomly rough plane and the scattered wave, due to roughness, is diffracted by the half-plane. The solution is transformed into a sum of the Fresnel integrals with complex arguments, an integral along the steepest descent path and a branch-cut integral, which are evaluated numerically. Then, intensities of the coherently scattered wave and incoherent wave are calculated in the region near the edge and illustrated in figures.     

High order numerical methods to two dimensional delta function integrals in level set methods

In this paper we design and analyze a class of high order numerical methods to delta function integrals appearing in level set methods in two dimensional case. The methods comprise approximating the mesh cell restrictions of the delta function integral. In each mesh cell the two dimensional delta function integral can be rewritten as a one dimensional ordinary integral with the smooth integrand being a one dimensional delta function integral, and thus is approximated by applying standard one dimensional high order numerical quadratures and high order numerical methods to one dimensional delta function integrals proposed in [X. Wen, High order numerical methods to a type of delta function integrals, J. Comput. Phys. 226 (2007) 1952–1967]. We establish error estimates for the method which show that the method can achieve any desired accuracy by assigning the corresponding accuracy to the sub-algorithms and has better accuracy under an assumption on the zero level set of the level set function which holds generally. Numerical examples are presented showing that the second to fourth order methods implemented in this paper achieve or exceed the expected accuracy and demonstrating the advantage of using our high order numerical methods.     

Diffraction of sound by an elastic or impedance sphere located near an impedance or elastic boundary of a halfspace

An exact solution is obtained to the problem of sound diffraction by an elastic or impedance sphere located near an impedance or elastic boundary of a halfspace. The problem is solved using the Helmholtz integral equation in which the field of a point source in the halfspace with an elastic boundary is used as the Green function. The diffracted field is represented as a series expansion in spherical harmonics. The expansion coefficients are determined from a set of independent algebraic systems of equations. The matrix coefficients of these systems are determined as integrals of the products of the associated Legendre polynomials on the complex plane with respect to the real and complex angles of the sound incidence on the halfspace boundary. To decrease the number of such integrals, expansions using the Klebsh-Gordon coefficients are applied. As a result, algorithms for calculating the scattered field in the halfspace are obtained.     

Wronskian Approach and One-Dimensional Schrodinger Equation with Double-Well Potential

A Wronskian determinant approach is suggested to study the energy and the wave function for onedimensional Schrodinger equation. An integral equation and its corresponding Green function are constructed. As an example, we employed this approach to study the problem of double-well potential with strong coupling. A series of expansion of ground state energy up to the second order approximation of iterative procedure is given.     

基于有限元法的光子并矢格林函数重整化及其在自发辐射率和能级移动研究中的应用

任意微纳结构中量子点的自发辐射率和能级移动均可用并矢格林函数表达.当源点和场点在同一位置时,格林函数的实部是发散的.为解决这一发散问题,可采用重整化格林函数方法.本文提出一种计算重整化格林函数和散射格林函数的方法.该方法利用有限元,计算点电偶极子的辐射场,将其在量子点体积内做平均得到重整化的并矢格林函数,减去均匀空间中解析的重整化格林函数,得到重整化的散射格林函数.在均匀空间情况下,本方法所得数值结果与解析解一致.将该方法应用到银纳米球系统,以解析的散射格林函数作为参考,结果表明该方法能准确处理散射格林函数的重整化问题.将该方法应用到表面等离激元纳米腔中,发现有极大的自发辐射增强和能级移动,且该结果不依赖于量子点的体积.这些研究在光与物质相互作用领域具有积极的意义.