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

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

非均质各向异性地层中方位随钻电磁测井响应三维有限体积法数值模拟算法

方位随钻电磁测井是一种能够实时探测地层边界、实现地质导向与井眼成像的新型测井技术. 本文根据方位随钻电磁测井仪器的典型线圈系结构, 首先引入柱坐标系下非均质完全各向异性地层中电流源并矢Green函数, 并利用电磁场叠加原理给出倾斜发射线圈激发的电场以及倾斜接收线圈上感应电动势的计算公式; 然后应用电流源电场并矢Green函数的混合势克服非均质地层中电磁数值模拟的低感应数问题, 通过ρ 和z方向上Lebedev网格设法降低网格节点个数, 并且利用标准化算法确定柱坐标系下非均质单元上的等效电导率. 在此基础上, 用三维有限体积法建立柱坐标系电场混合势的离散方法, 得到一个交错网格上电场矢势和标势大型代数方程, 并用不完全LU分解以及稳定双共轭梯度法确定数值解. 最后, 通过数据模拟结果对算法的有效性进行检验, 并考察钻铤、线圈倾斜角度以及地层各向异性等参数对仪器响应的影响. 数值结果表明: 在柱坐标系下用三维有限体积法的数值模拟算法处理非均质各向异性层中方位随钻电磁测井响应可以得到很好的结果. 钻铤、电导率各向异性、层边界均对方位随钻电磁波测井响应产生较大的影响; 在电阻率较大的地层, 幅度比和相位差响应越小; 发射线圈和接收线圈同时倾斜时, 幅度比和相位差响应受地层的影响更灵敏.     

Equivalence of integral and integrofunctional equations in electrodynamic problems of diffraction by inhomogeneous bodies

Relations for fields with specified distributions of volume electric and magnetic currents are investigated. The dependence of the term outside the integral on the limit method for calculation of improper integrals is shown; an evaluation method is proposed. The relations are used to establish equivalence between integral equations with volume polarization currents and integrofunctional equations of the standard boundary-value problem of electrodynamic for an inhomogeneous ferrite body. It is shown that the solutions of both types of equations satisfy the boundary conditions. The practical use of integrofunctional equations in the form of systems of integral equations of the second kind and functional equations is discussed.Kazan' Aviation Institute. Translated from Izvestiya Vysshikh Uchebnykh Zavedenii, Radiofizika, Vol. 34, No. 8, pp. 936–946, August, 1991.     

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

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

Radiative transfer in absorbing, emitting and linearly anisotropic-scattering inhomogeneous cylindrical medium

Radiative integral transfer equations for one-dimensional solid cylinder with absorbing, emitting and linearly anisotropic-scattering inhomogeneous medium were derived by Abulwafa et al. (JQSRT 62 (1999) 755). The anisotropic terms in the integral equations and their results for anisotropic benchmark problems (JQSRT 66 (2000) 487) are inaccurate. In this study, the integral equations for absorbing, emitting and linearly anisotropic-scattering medium are rederived, and the integral equations for one-dimensional solid cylindrical medium are solved. The results are compared with selected cases using the discrete ordinates S₁₆ and the exact solutions available in the literature.     

Efficient omnidirectional couplers achieved by anisotropic photonic crystal waveguides

Photonic crystals (PCs) have many potential applications because of their ability to control light-wave propagation. We have investigated omnidirectional couplers in two-dimensional anisotropic PC structures. The anisotropic PC coupler composed of an anisotropic-dielectric cylinder in air is studied by solving Maxwell's equations using the plane wave expansion method and finite-difference time-domain method. The photonic band structures are found to exhibit absolute bandgaps for the square lattices. Numerical simulations show that the incident light-waves at both transverse electric and transverse magnetic modes have efficient coupling in anisotropic PC coupler with square lattices. The guided modes and coupling length are analyzed by considering various line defect anisotropic PC waveguides and interaction regions of couplers. Such a mechanism of omnidirectional coupling should open up a new application for designing components in photonic integrated circuits.     

Analytical methods for modeling of ultrasonic nondestructive testing of anisotropic media

Many modern structural materials exhibit anisotropic elastic behavior leading to complicated wave propagation phenomena. To ensure the reliability of ultrasonic nondestructive testing techniques, these material properties as well as the influence of microstructural inhomogeneities and the effects of interfaces on ultrasonic wave propagation have to be taken into account. In this respect, mathematical modeling provides an efficient method of assisting analysis. Two computationally efficient analytical approaches--a Gaussian beam and a point source superposition technique--are presented, which are well-suited for performing ultrasonic wave propagation and scattering simulations for anisotropic media. Results for homogeneous as well as inhomogeneous anisotropic media like composites and weld material are presented.     

Split-step finite-difference time-domain method with perfectly matched layers for efficient analysis of two-dimensional photonic crystals with anisotropic media

This Letter presents a split-step (SS) finite-difference time-domain (FDTD) method for the efficient analysis of two-dimensional (2-D) photonic crystals (PhCs) with anisotropic media. The proposed SS FDTD method is formulated with perfectly matched layer boundary conditions and caters for inhomogeneous anisotropic media. Furthermore, the proposed method is derived using the efficient SS1 splitting formulas with simpler right-hand sides that are more efficient and easier to implement. A 2-D PhC cavity with anisotropic media is used as an example to validate the efficiency of the proposed method.     

A kernel-free boundary integral method for elliptic boundary value problems

This paper presents a class of kernel-free boundary integral (KFBI) methods for general elliptic boundary value problems (BVPs). The boundary integral equations reformulated from the BVPs are solved iteratively with the GMRES method. During the iteration, the boundary and volume integrals involving Green’s functions are approximated by structured grid-based numerical solutions, which avoids the need to know the analytical expressions of Green’s functions. The KFBI method assumes that the larger regular domain, which embeds the original complex domain, can be easily partitioned into a hierarchy of structured grids so that fast elliptic solvers such as the fast Fourier transform (FFT) based Poisson/Helmholtz solvers or those based on geometric multigrid iterations are applicable. The structured grid-based solutions are obtained with standard finite difference method (FDM) or finite element method (FEM), where the right hand side of the resulting linear system is appropriately modified at irregular grid nodes to recover the formal accuracy of the underlying numerical scheme. Numerical results demonstrating the efficiency and accuracy of the KFBI methods are presented. It is observed that the number of GMRES iterations used by the method for solving isotropic and moderately anisotropic BVPs is independent of the sizes of the grids that are employed to approximate the boundary and volume integrals. With the standard second-order FEMs and FDMs, the KFBI method shows a second-order convergence rate in accuracy for all of the tested Dirichlet/Neumann BVPs when the anisotropy of the diffusion tensor is not too strong.     

A kernel-free boundary integral method for elliptic boundary value problems

This paper presents a class of kernel-free boundary integral (KFBI) methods for general elliptic boundary value problems (BVPs). The boundary integral equations reformulated from the BVPs are solved iteratively with the GMRES method. During the iteration, the boundary and volume integrals involving Green’s functions are approximated by structured grid-based numerical solutions, which avoids the need to know the analytical expressions of Green’s functions. The KFBI method assumes that the larger regular domain, which embeds the original complex domain, can be easily partitioned into a hierarchy of structured grids so that fast elliptic solvers such as the fast Fourier transform (FFT) based Poisson/Helmholtz solvers or those based on geometric multigrid iterations are applicable. The structured grid-based solutions are obtained with standard finite difference method (FDM) or finite element method (FEM), where the right hand side of the resulting linear system is appropriately modified at irregular grid nodes to recover the formal accuracy of the underlying numerical scheme. Numerical results demonstrating the efficiency and accuracy of the KFBI methods are presented. It is observed that the number of GMRES iterations used by the method for solving isotropic and moderately anisotropic BVPs is independent of the sizes of the grids that are employed to approximate the boundary and volume integrals. With the standard second-order FEMs and FDMs, the KFBI method shows a second-order convergence rate in accuracy for all of the tested Dirichlet/Neumann BVPs when the anisotropy of the diffusion tensor is not too strong.     

高频电磁干扰对传输线耦合全波建模方法

对于场线耦合问题,经典传输线理论不适用于求解高频电磁干扰辐照下传输线负载上的电压和电流响应。针对这一问题,首先介绍了一种基于天线理论和模拟行为建模（ABM）的时域全波建模方法。该方法利用Harrington矩量法将电流积分方程离散并推导得到宏模型时域表达式,然后利用ABM频域功能实现频变参数的傅里叶逆变换和时域卷积计算。利用电路求解器,该建模方法可直接求解任意结构传输线耦合的负载处瞬态响应;与传统全波算法相比,模型一旦建立便可应用于任意入射场和线性/非线性负载的情况,无需重复耗时地求解电流积分方程。该方法可简化全波算法求解过程,提高仿真计算效率,尤其便于在入射场和负载存在不确定参数时进行高效重复抽样计算以获得统计特性。然后以高频电磁干扰耦合有损大地上的双导体传输线为例,通过与数值电磁代码和传统传输线理论方法的求解结果对比,验证了所提宏模型的有效性以及传输线理论在解决场线耦合问题时的局限性。结果表明,基于全波方法构建的宏模型可在时域内高效准确地求解高频电磁干扰辐照下任意形状传输线负载上的瞬态响应。     

Decoherence Effect in An Anisotropic Two-Qubit Heisenberg XYZ Model with Inhomogeneous Magnetic Field

Taking the decoherence effect into account, the entanglement evolution of a two-qubit anisotropic Heisenberg XYZ chain in the presence of inhomogeneous magnetic field is investigated. The time evolution of concurrence is studied for the initial state cos θ｜01） ＋ sin θ｜10） at zero temperature. The influences of inhomogeneous magnetic field, anisotropic parameter and decoherence on entanglement dynamic are addressed in detail, and a concurrence formula of the steady state is found. It is shown that the entanglement sudden death （ESD） and entanglement sudden birth （ESB） appear with the decoherence effect, and the stable concurrence depends on the uniform magnetic field B, anisotropic parameter △ and environment coupling strength γ, which is independent of different initial states and nonuniform magnetic field b.     

High-performance modeling acoustic and elastic waves using the parallel Dichotomy Algorithm

A high-performance parallel algorithm is proposed for modeling the propagation of acoustic and elastic waves in inhomogeneous media. An initial boundary-value problem is replaced by a series of boundary-value problems for a constant elliptic operator and different right-hand sides via the integral Laguerre transform. It is proposed to solve difference equations by the conjugate gradient method for acoustic equations and by the GMRES(k) method for modeling elastic waves. A preconditioning operator was the Laplace operator that is inverted using the variable separation method. The novelty of the proposed algorithm is using the Dichotomy Algorithm [26], which was designed for solving a series of tridiagonal systems of linear equations, in the context of the preconditioning operator inversion. Via considering analytical solutions, it is shown that modeling wave processes for long instants of time requires high-resolution meshes. The proposed parallel fine-mesh algorithm enabled to solve real application seismic problems in acceptable time and with high accuracy. By solving model problems, it is demonstrated that the considered parallel algorithm possesses high performance and efficiency over a wide range of the number of processors (from 2 to 8192).     

Scattering Parameters of FSS on Anisotropic Layers at Millimeter Wave Frequencies

The analysis of periodic structures on anisotropic substrates is presented for incident millimeter waves. Frequency selective surface structures composed by a periodic array of conducting crossed dipoles on anisotropic layers are analyzed. A full-wave analysis for the problem of scattering is performed. The formulation uses the moment method in combination with the immitance approach, in the Fourier domain, to determine the scattering parameters of the considered structures. A set of entire-domain type basis functions is used in the expansion for the unknown induced current. This analysis is computationally efficient, and its accuracy is verified by the agreement between results obtained in this work and those available in the literature, for particular cases.     

Modeling magnetic photonic crystals with lossy ferrites using an efficient complex envelope alternating-direction-implicit finite-difference time-domain method

In this Letter, we present an efficient complex-envelope alternating-direction-implicit finite-difference time-domain (CE-ADI-FDTD) method for the transient analysis of magnetic photonic crystals with lossy ferrites. The proposed CE-ADI-FDTD method is generally formulated for a saturated ferrite with anisotropic permittivity tensor and ferrite loss. Auxiliary differential equations for modeling saturated ferrite and Maxwell's curl equations are first cast into a first-order differential system in a CE form. Then, by using an efficient ADI splitting formulas, the proposed CE-ADI-FDTD method is attained in a very concise form with few and simple right-hand side terms. The performance of the proposed method is validated and compared with the explicit FDTD method.     

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

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

Numerical calculation of the magnetic field and magnetic moment of ferromagnetic bodies with a complex spatial configuration

A method for numerical calculation of the magnetic field and magnetic moment of 3D ferromagnetic objects is described using the method of finite volume and absorbing boundary conditions. The magnetic moments and magnetic field of some ferromagnetic systems imitating various objects are calculated. It is shown that the configuration of ferromagnetic structures consisting of various combinations of rods determines their magnetic moments.     

Analysis of stress field in a single anisotropic diamond/square shaped inclusion using the volume integral equation method and the numerical equivalent inclusion method

A volume integral equation method (VIEM) is used to study elastostatic problems in an unbounded elastic solid containing a single diamond/square shaped inclusion subject to uniform tensile stress at infinity. The inclusion is assumed to be a long parallel diamond/square cylinder composed of isotropic or anisotropic elastic materials and perfectly bonded to the isotropic matrix. The solid is assumed to be under plane strain on the plane normal to the cylinder. A detailed analysis of the stress field at the interface between the isotropic matrix and the single isotropic/orthotropic diamond/square shaped inclusion is carried out. The effects of a single isotropic/orthotropic diamond/square shaped inclusion on the stress field at the interface between the matrix and the inclusion are investigated in detail. The accuracy of the volume integral equation method for the interfacial stress field is validated and compared by the numerical equivalent inclusion method (NEIM) and the finite element method (FEM) using ADINA. Through detailed analysis of plane elastostatic problems using the parallel volume integral equation method (PVIEM) in an unbounded isotropic matrix with multiple isotropic diamond shaped inclusions under uniform remote tensile loading, it is demonstrated that the volume integral equation method can also be applied to solve general two- and three-dimensional elastostatic problems involving multiple isotropic/anisotropic inclusions whose shape and number are arbitrary.     

Spatially modulated antiferromagnetic structures in an easy-plane multiferroic

Possible types of spatially modulated periodic antiferromagnetic structures in a uniaxial rhombohedral multiferroic with BiFeO₃ crystal symmetry have been studied depending on the ratio of the uniaxial anisotropy and magnetoelectric interaction parameters. It has been shown that, along with symmetric cycloid antiferromagnetic structures with zero transverse component of the antiferromagnetism vector, there are changes in the antiferromagnetism vector direction with both right and left nonzero components of the antiferromagnetic moment, which are branched from the high-symmetry spatially modulated distribution. These solutions degenerate into a homogeneous state at a critical value of the normalized easy-plane anisotropy parameter. The existence of the found spatially inhomogeneous antiferromagnetic states with an incommensurate period can lead to additional features in magnetoelectric properties in multiferroics of the type under consideration near magnetic phase transitions in electric and magnetic fields.     

Inverse problem of reflection of electromagnetic waves from layered inhomogeneous plasma-like structures with refractive-index discontinuity

To reconstruct the distribution of free carriers in a layered inhomogeneous plasma-like half-space during the electromagnetic-wave reflection, we derive the Gelfrand-Levitan-Marchenko integral equation for structures with refractive-index discontinuity by the inverse problem method. Exact model solutions of the generalized equations are analyzed. The applicability limits and scope of these exact solutions for reconstruction problems are shown. __________ Translated from Izvestiya Vysshikh Uchebnykh Zavedenii, Radiofizika, Vol. 50, No. 1, pp. 31–39, January 2007