Numerical method for problems of diffraction on semi-infinite structures

We propose a method which transforms homogeneous integral equations into inhomogeneous ones for problems of diffraction by semi-infinite structures. New integral equations and the corresponding stationary functionals dependent on the desired scattering parameters are obtained. The consideration is performed for the open-end diffraction problem of a parallel-plate waveguide which has a rigorous solution, but the method has sufficient generality to use it for two-dimensional surface integral equations describing planar and nonplanar structures as well as for an arbitrary structure of waveguide transformer type, the solution for which can be sought in the finite domain. The method is based on field representation at the infinity as incident and scattered waves. __________ Translated from Izvestiya Vysshikh Uchebnykh Zavedenii, Radiofizika, Vol. 49, No. 3, pp. 235–245, March 2006.     

Emptiness formation probability in the domain-wall six-vertex model

The emptiness formation probability in the six-vertex model with domain wall boundary conditions is considered. This correlation function allows one to address the problem of limit shapes in the model. We apply the quantum inverse scattering method to calculate the emptiness formation probability for the inhomogeneous model. For the homogeneous model, the result is given both in terms of certain determinant and as a multiple integral representation.     

Solution of the kinetic equation for thin two-and three-dimensional metallic slabs with diffuse boundary scattering and arbitrary impurity scattering in the bulk

The integrodifferential equation in the phase space, describing the response to a homogeneous electric field of two-or three-dimensional isotropic metallic slab with diffuse scattering on the border and arbitrary scattering on impurities, is transformed into the system of inhomogeneous integral equations in the configuration space. It is shown that one of the functions fulfilling this system is simply connected with the current response of the slab. The general solution of this system is found in the form of asymptotical expansion provided that the mean-free path is much greater than the width of the slab. In the opposite limit, the asymptotic solution is found only at almost isotropic impurity scattering. In other words, the Fuchs problem (1938) has been solved for thin slabs at a real low-temperature collision integral.     

Exact wave field simulation for finite-volume scattering problems

An exact boundary condition is presented for scattering problems involving spatially limited perturbations of arbitrary magnitude to a background model in generally inhomogeneous acoustic media. The boundary condition decouples the wave propagation on a perturbed domain while maintaining all interactions with the background model, thus eliminating the need to regenerate the wave field response on the full model. The method, which is explicit, relies on a Kirchhoff-type integral extrapolation to update the boundary condition at every time step of the simulation. The Green's functions required for extrapolation through the background model are computed efficiently using wave field interferometry.     

On multidimensional ultrasonic scattering in an inhomogeneous elastic background.

This work is concerned with the modeling of elastic wave scattering by solid or fluid-filled objects embedded in an inhomogeneous elastic background. The medium is probed by a monochromatic force and the scattered field is computed (forward problem) or observed (inverse problem) at some known receiver locations. Based on vector integral equations for elastic scattering, a general framework is developed, independent of both the problem geometry and the transmitter-receiver characteristics. This framework encompasses both forward and inverse modeling. In the forward model, a Born approximation for an inhomogeneous background is applied to obtain a closed form expression for the scattered field. In the inverse model, this approximation is also invoked to linearize for the multiparameter characteristic of the object. Finally, an iterative inversion scheme alternating forward and inverse modeling is proposed to improve the resolution and accuracy of the reconstruction algorithm.     

A unified optical theorem for scalar and vectorial wave fields

The generalized optical theorem is an integral relation for the angle-dependent scattering amplitude of an inhomogeneous scattering object embedded in a homogeneous background. It has been derived separately for several scalar and vectorial wave phenomena. Here a unified optical theorem is derived that encompasses the separate versions for scalar and vectorial waves. Moreover, this unified theorem also holds for scattering by anisotropic elastic and piezoelectric scatterers as well as bianisotropic (non-reciprocal) EM scatterers.     

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.     

Particle transfer problem in a two-region inhomogeneous slab with generalized boundary condition

The problem of particle transfer in a two-region inhomogeneous slab of equal or unequal thickness has been solved with general boundary conditions. A Galerkin-iterative technique is used to solve the coupled integral equations for two-regions. Numerical results are obtained for the reflection and the transmission coefficients for the homogeneous and inhomogeneous cases with diffusely and specularly reflecting boundaries condition. Results obtained for homogeneous and inhomogeneous medium are compared well with the published calculations.     

A Riemann–Hilbert Problem for Propagation of Electromagnetic Waves in an Inhomogeneous,Dispersive Ω Waveguide

We consider the inverse scattering problem for a model of electromagnetic wave propagation in a rectangular waveguide filled with dispersive material. The waveguide is inhomogeneous in the longitudinal direction but homogeneous in the transverse directions. Dispersive properties of the material are described by a single-resonance Lorentz model. By reformulating the scattering problem in the frequency domain as a Riemann–Hilbert problem, we prove that the constitutive parameters of the inhomogeneous waveguide are reconstructed uniquely from the scattering data.     

Pseudopulse near-field subsurface tomography

Decisive success has been achieved in developing the subsurface near-field scanning tomography that overcomes the Rayleigh diffraction limit of a resolution. It is related to the transformation of the multifrequency inverse scattering problem to that for a complex-valued synthesized pulse (pseudopulse). It leads to the integral equation that has maxima in the depth dependence of its kernel and, hence, to the much better depth resolution of tomography. Moreover, the noise related to surface scattering is mainly suppressed in such an approach. This idea is realized here in the microwave subsurface tomography of 3D inhomogeneous dielectric structures. For homogeneous dielectric targets, this approach is applied to obtain holography images of their shape.     

A NUMERICAL METHOD FOR SCATTERING FROM ACOUSTICALLY SOFT AND HARD THIN BODIES IN TWO DIMENSIONS

This paper presents a numerical method for predicting the acoustic scattering from two-dimensional (2-D) thin bodies. Both the Dirichlet and Neumann problems are considered. Applying the thin-body formulation leads to the boundary integral equations involving weakly singular and hypersingular kernels. Completely regularizing these kinds of singular kernels is thus the main concern of this paper. The basic subtraction-addition technique is adopted. The purpose of incorporating a parametric representation of the boundary surface with the integral equations is two-fold. The first is to facilitate the numerical implementation for arbitrarily shaped bodies. The second one is to facilitate the expansion of the unknown function into a series of Chebyshev polynomials. Some of the resultant integrals are evaluated by using the Gauss-Chebyshev integration rules after moving the series coefficients to the outside of the integral sign; others are evaluated exactly, including the modified hypersingular integral. The numerical implementation basically includes only two parts, one for evaluating the ordinary integrals and the other for solving a system of algebraic equations. Thus, the current method is highly efficient and accurate because these two solution procedures are easy and straightforward. Numerical calculations consist of the acoustic scattering by flat and curved plates. Comparisons with analytical solutions for flat plates are made.     

Anisotropic scattering in inhomogeneous media—1. A new representation formula for the solution of the auxiliary integral equation for the source function

A new representation formula for the solution of the auxiliary integral equation for the source function in inhomogeneous, anisotropically scattering media is presented. It involves two new functions Φ and ψ of two variables instead of the original five variables. This generalizes earlier results of Kagiwada et al. (1969) and Sobolev (1972) applicable to homogeneous atmospheres. The corresponding Bellman-Krein formula for the resolvent kernel is also derived. The present representation for the solution of Fredholm integral equations of second kind with unsymmetric kernels provides a new approach to radiative transfer in anisotropic inhomogeneous media.     

A method of cauchy integral equation for non-coherent transfer in half-space

A technique is presented which allows easy construction of solutions for various half-space problems arising in non-coherent radiative transfer with complete redistribution. By use of an inverse Laplace transform method, Wiener-Hopf integral equations are reduced to Cauchy-type singular integral equations. The factorization technique used by Case and Zweifel for coherent scattering can then be carried over to non-coherent transfer. The method is applied to the inhomogeneous integral equation for the source function of a two-level atom, previously solved by Ivanov. It is also applied to the conservative, homogeneous case and to singular Wiener-Hopf equations arising from asymptotic expansions in the limit of vanishing probability of collisional destruction ?. Consequences for the scaling laws in a finite slab are examined in a companion paper.     

Method for solving the problems of multiple scattering by several bodies in a homogeneous unbounded medium

A method is developed for solving problems of multiple scattering by an aggregate of bodies in a homogeneous unbounded medium. For this purpose, the problem on the multiple scattering produced by two bodies in the field of a plane wave is first considered under the assumption that the initial unperturbed scattering amplitudes of both scatterers are known. The solution is constructed by considering plane waves multiply rescattered by the scatterers. Integral equations are obtained that allow one to calculate the resulting scattering amplitude of each scatterer and the combined scattering amplitude of the system of two scatterers. It is shown that knowledge of the solution to this problem is sufficient to solve the problem on the scattering field of a system consisting of an arbitrary number of scatterers. Expressions for the scattering amplitude in the case of an arbitrary primary field are presented. The relationship between the integral equations describing the multiple scattering in a homogeneous space and the multiple scattering by a single scatterer located near an interface is demonstrated. Approximate expressions are given for calculating the scattering amplitude in the case of multiple scattering.     

Some comments on modeling the linearized Boltzmann equation

Some exact solutions of the homogeneous and the inhomogeneous linearized Boltzmann equation (LBE) for rigid-sphere collisions are used to define two model equations in the general area of rarefied-gas dynamics. These equations are obtained from a systematic development of two synthetic scattering kernels that yield model equations that have as exact solutions certain known exact solutions of the homogeneous and of the inhomogeneous LBE. The first model established is defined in terms of the collisional invariants and the Chapman-Enskog integral equations for viscosity and for heat conduction. An extended model is defined also in terms of the collisional invariants and the Chapman-Enskog functions for viscosity and heat conduction, but the first and second Burnett functions are also included in the model. The variable collision frequency or generalized BGK model is also obtained as a special case. In addition, the exact mean-free paths defined, for rigid-sphere collisions and the LBE, in terms of viscosity or heat conduction are employed to define approximations of these quantities that are consistent with the use of the variable collision frequency model.     

金属圆柱腔体中使用非均一背景增强微波断层成像

针对基于圆柱金属腔体的微波断层成像系统,提出了一种利用非均一背景增强系统获取目标信息能力的方法.该方法通过在腔体内放置已知物体构成非均一背景,这样不但能利用背景的先验信息,而且可以增加等效辐射源对目标进行探测.首先,利用矩量法计算圆柱金属腔体内非均一背景的格林函数和离散积分算子,并对离散积分算子进行奇异值谱和条件数分析,在理论上证明该方法的可行性;然后,利用基于有限元的对比源逆成像法对均一背景、有耗非均一背景和无耗非均一背景三种情况进行仿真研究;最后对仿真结果进行了误差分析和比较.仿真结果表明,该方法可以提高反演收敛速度和结果准确度,有耗非均一背景略优于无耗非均一背景.该方法可以在不改变硬件系统和算法的情况下得到更准确的反演结果,可应用于医学成像与工业无损探测.     

Unitarity and the inhomogeneous equations method for identical particle scattering

A K-matrix solution to the coupled, inhomogeneous equations describing the scattering of a particle by a system of identical particles is developed. It is shown that K is a sum of two terms, one arising from the homogeneous solution and one from the particular integral. The former is a direct contribution, i.e., with no exchange, while the latter is a pure exchange contribution. Thus, as in the previously studied case of the T matrix arising from this system of equations, the direct and exchange portions of K are additive, and can be computed separately. A unitary S matrix is obtained from K in the usual way: S = (1 + iK)(1 − iK)⁻¹. The problem of how to calculate K when an apparent two-channel problem is actually a two-particle problem with the channels referring to the identical particle labels is also solved.     

一种金属腔体中微波断层成像的最优分层非均一背景

针对基于金属腔体的微波断层成像系统,提出了一种最优分层非均一背景的设计方法.该方法使用一种新的微波断层成像积分算子评价方法和模拟退火法等最优化方法.首先,介绍了一种基于有限元法的微波断层成像积分算子计算方法.然后,提出一种新的微波断层成像积分算子度量,该度量可以综合评价整个积分算子奇异值谱,并通过一组仿真研究证明该度量与反演结果的误差具有相关性;该度量用一个数值综合评价一个积分算子,可以方便地应用于最优化算法中;利用模拟退火法选择圆形金属腔体中分层非均一背景的每一层介质的相对介电常数,从而获得一个最优分层非均一背景.最后,对尺寸小于半波长的圆柱目标和"凹"字形复杂目标进行仿真研究,仿真结果证明该最优分层非均一背景可以提高微波断层成像算法的收敛速度,提高反演结果的准确性.     

A simple reduction procedure for the three- and N-body scattering equations

A generalized form of the two-body Kowalski-Noyes method is shown to provide a both simple and powerful unitary reduction of the three- and N-body scattering equations. Employing generalized half-off-shell functions that satisfy of-sshell but real and non-singular integral equations, the reduction directly leads to on-shell integral equations for the scattering amplitudes. Physically, it is simple example of how the scattering problem can be split into an internal and an external part.