A boundary integral method to compute Green's functions for the radiative transport equation

The Green's function for the time-independent radiative transport equation in the whole space can be computed as an expansion in plane wave solutions. Plane wave solutions are a general class of solutions for the radiative transport equation. Because plane wave solutions are not known analytically in general, we calculate them numerically using the discrete ordinate method. We use the whole space Green's function to derive boundary integral equations. Through the solution of the boundary integral equations, we compute the Green's function for bounded domains. In particular we compute the Green's function for the half space, the slab, and the two-layered half space. The boundary conditions used here are in their most general form. Hence, this theory can be applied to boundaries with any kind of reflection and transmission law.     

Solutions of the Dirac equation and Green's electron function in an external electromagnetic field of special form

A complete orthonormal set of c-numerical solutions of the Dirac equation is constructed and integral representations are obtained for the Green's electron function in an external electromagnetic field representing the combination of a longitudinal electrical wave and a plane wave being propagated in one direction (along the x₃ axis).     

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

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

Scattering of a scalar wave from a two-dimensional randomly rough Neumann surface

Abstract

We present a reciprocity and unitarity preserving formulation of the scattering of a scalar plane wave from a two-dimensional, randomly rough surface on which the Neumann boundary condition is satisfied. The theory is formulated on the basis of the Rayleigh hypothesis in terms of a single-particle Green's function G(q_‖|k_‖) for the surface electromagnetic waves that exist at the surface due to its roughness, where k_‖ and q_‖ are the projections on the mean scattering plane of the wave vectors of the incident and scattered waves, respectively. The specular scattering is expressed in terms of the average of this Green's function over the ensemble of realizations of the surface profile function (G(q_‖|k_‖)). The Dyson equation satisfied by (G(q_‖|k_‖)) is presented, and the properties of the solution are discussed, with particular attention to the proper self-energy in terms of which the averaged Green's function is expressed. The diffuse scattering is expressed in terms of the ensemble average of a two-particle Green's function, which is the product of two single-particle Green's functions. The Bethe-Salpeter equation satisfied by the averaged two-particle Green's function is presented, and properties of its solution are discussed. In the small roughness limit, and with the irreducible vertex function approximated by the sum of the contribution from the maximally-crossed diagrams, which represent the coherent interference between all time-reversed scattering sequences, the solution of the Bethe-Salpeter equation predicts the presence of enhanced backscattering in the angular dependence of the intensity of the waves scattered diffusely.     

Small-slope expansion for electromagnetic-wave diffraction on a rough surface

Abstract

The diffraction and absorption of the plane electromagnetic wave on a rough surface is considered to find the scattering and emissivity of the surface. For this purpose a system of integral equations for unknown surface fields is derived from Green's formula for the Helmholtz equation. The small-slope approach is used to find a solution, i.e. the solution is determined from an expansion over the roughness spectrum that, in the limit of the large-scale roughness, turns out to be the expansion over the slope spectrum.     

Green's electron function in a quantized plane wave field

It is shown that the Green's function of an electron that interacts with a quantized plane wave can be expressed in terms of the corresponding Green's function of a scalar particle. By using the known expression for the Green's function of a scalar particle, an integral representation is found with respect to the intrinsic time for the Green's electron function in a quantized plane wave of arbitrary form.     

The Hadamard Construction of Green's Functions on a Curved Space-Time with Symmetries

The Hadamard constituents of Green's functions for a ζ-parametrized generalization of the massless scalar d'Alembert equation to a curved space-time including the conformally invariant wave equation: the world function of space-time, the transport scalar, and the tail-term coefficients, being simultaneously coefficients in the Schwinger-DeWitt expansion of the Feynman propagator for the corresponding invariant Klein-Gordon equation, are considered on a general static spherically symmetric and (2,2)-decomposable metric. The construction equations determining the Hadamard building elements are cast into a symmetry-adapted form and used to obtain, on a specific model metric, exact explicit solutions.     

Use of computer simulations to determine steady-state distributions in ion bombardment

Abstract

Computer simulations typically determine the particle-distribution function corresponding to a single-event point source. This distribution function contains all the information needed to determine the corresonding steady-state distribution. From the single-event point-source distribution one can also deduce the distributions resulting from plane and volume-homogeneous sources, both for single event and steady state. A formalism for making these deductions is developed, based on the observation that the distribution determined by computer simulation obeys the Boltzmann transport equation, and on the properties of the Green's functions solutions of that equation, particularly the consequences of space- and time-translation invariance. The treatment deals explicitly with a homogeneous target infinite in every direction, but is easily generalized to cases with boundary conditions, such as a beam bombarding a layered target in slab geometry.     

Simulation of electromagnetic scattering by random discrete particles using a hybrid FE-BI-CBFM technique

A hybrid finite element–boundary integral–characteristic basis function method (FE-BI-CBFM) is proposed for an efficient simulation of electromagnetic scattering by random discrete particles. Specifically, the finite element method (FEM) is used to obtain the solution of the vector wave equation inside each particle and the boundary integral equation (BIE) using Green's functions is applied on the surfaces of all the particles as a global boundary condition. The coupling system of equations is solved by employing the characteristic basis function method (CBFM) based on the use of macro-basis functions constructed according to the Foldy–Lax multiple scattering equations. Due to the flexibility of FEM, the proposed hybrid technique can easily deal with the problems of multiple scattering by randomly distributed inhomogeneous particles that are often beyond the scope of traditional numerical methods. Some numerical examples are presented to demonstrate the validity and capability of the proposed method.     

有界空间中的非线性声学和收敛的积累解

在微扰近似下,拉格朗日体系下的一阶、二阶波动方程解是具势运动,应用拉格朗日变动参数法来寻求积累解。在一般情况下,二阶波的波动方程在半空间会出现各式各样的积累解,它们沿着3个坐标变量的方向都有积累,在理想介质中它们不满足辐射条件。本文的结果表明,在考虑到介质的非理想性之后,也只有沿着平面边界面法线方向有积累的积累解才满足辐射条件,因而是收敛的。      

Tempered wave functions: Schrödinger's equation within the light cone

The standard derivation of Schrödinger's equation from a Lorentz-invariant Feynman path integral consists in taking first the limit of infinite speed of light and then the limit of short time slice. In this order of limits the light cone of the path integral disappears, giving rise to an instantaneous spread of the wave function to the entire space. We ascribe the failure of the propagation in time according to Schrödinger's equation to retain the light cone of the path integral to the very nature of the limiting process: it is a regular expansion of a singular approximation problem, because the time-dependent boundary conditions of the path integral on the light cone are lost in this limit. We propose a distinguished limit of the time-sliced relativistic path integral, which produces Schrödinger's equation and preserves the zero boundary conditions on and outside the original light cone of the path integral. This produces an intermediate model between non-relativistic and relativistic mechanics of a single particle quantum particle. These boundary conditions relieve the solutions of Schrödinger's equation in the entire space of several annoying, seemingly unrelated unphysical artifacts, including non-analytic wave functions, spontaneous appearance of discontinuities, non-existence of moments when the initial wave function has a jump discontinuity (e.g., a collapsed wave function after a measurement), and so on. The practical implications of the present formulation are yet to be seen.     

Specular propagation over rough surfaces: numerical assessment of Uscinski and Stanek's mean Green's function technique

The purpose of this paper is to numerically evaluate the effectiveness and accuracy of Uscinski and Stanek's mean Green's function technique for computing the mean field of a wave scattered by a rough surface. We present here a direct comparison of this technique with a rigorous numerical method, the forward scattering integral equation method, and another analytical method, the first-order smoothing approximation. Furthermore, we compare the roughness generated equivalent admittance using the three methods. Numerical computations reveal that the scattered field calculated by this technique is not accurate particularly for the equivalent admittance at low grazing angles, even though the mean surface current density is recovered when the wave has traversed several correlation lengths on the surface.     

Analytical approach for solving the radiative transfer equation in two-dimensional layered media

This study presents an analytical approach for obtaining Green's function of the two-dimensional radiative transfer equation to the boundary-value problem of a layered medium. A conventional Fourier transform and a modified Fourier series which is defined in a rotated reference frame are applied to derive an analytical solution of the radiance in the transformed space. The Monte Carlo method was used for a successful validation of the derived solutions.     

基于特征基函数的球面共形微带天线阵列分析

采用全波分析法对球面共形微带天线阵列进行了分析.相比体-面积分方程,采用球并矢格林函数的面积分方程法可以大幅减少未知量的数目,进而缓解计算机内存压力.微带天线阵列表面采用曲面三角形剖分,可较精确地模拟球面特性.首先,引入边界电荷以及半Rao-Wilton-Glisson基函数,成功实现了探针馈电,并采用镜像法解决了馈电边处线积分奇异问题.然后,采用特征基函数法降低了阻抗矩阵的阶数,并采取有效措施进一步节省内存和计算时间.最后,分析计算了不同尺寸的球面共形微带天线阵列的输入阻抗及远区场特性.与文献和仿真软件结果进行比较,证明了所提出的处理方法的正确性和有效性.     

The prediction of jet noise ground effects using an acoustic analogy and a tailored Green's function

An assessment of an acoustic analogy for the mixing noise component of jet noise in the presence of an infinite surface is presented. The reflection of jet noise by the ground changes the distribution of acoustic energy and is characterized by constructive and destructive interference patterns. The equivalent sources are modeled based on the two-point cross-correlation of the turbulent velocity fluctuations and a steady Reynolds-Averaged Navier–Stokes (RANS) solution. Propagation effects, due to reflection by the surface and refraction by the jet shear layer, are taken into account by calculating the vector Green's function of the linearized Euler equations (LEE). The vector Green's function of the LEE is written in relation to that of Lilley's equation; that is, it is approximated with matched asymptotic solutions and Green's function of the convective Helmholtz equation. The Green's function of the convective Helmholtz equation in the presence of an infinite flat plane with impedance is the Weyl–van der Pol equation. Predictions are compared with measurements from an unheated Mach 0.95 jet. Microphones are placed at various heights and distances from the nozzle exit in the peak jet noise direction above an acoustically hard and an asphalt surface. The predictions are shown to accurately capture jet noise ground effects that are characterized by constructive and destructive interference patterns in the mid- and far-field and capture overall trends in the near-field.     

Introducing a Green–Volterra series formalism to solve weakly nonlinear boundary problems: Application to Kirchhoff's string

This paper introduces a formalism which extends that of “Green's function” and that of “the Volterra series”. These formalisms are typically used to solve, respectively, linear inhomogeneous space–time differential equations in physics and weakly nonlinear time-differential input-to-output systems in automatic control. While Green's function is a space–time integral kernel which fully characterizes a linear problem, the Volterra series expansions involve a sequence of multi-variate time integral kernels (of convolution type for time-invariant systems). The extension proposed here consists in combining the two approaches, by introducing a series expansion based on multi-variate space–time integral kernels. This series allows the representation of the space–time solution of weakly nonlinear boundary problems excited by an “input” which depends on space and time.     

Field of a point source in a semi-infinite elastic medium

An exact solution for the tensor Green's function of a harmonic field for a semi-infinite elastic medium is presented in an easy-to-use form in the theory of wave scattering. The solution is derived in the form of a sum of the Green's functions for an infinite medium and the term satisfying the homogeneous wave equation for a semi-infinite elastic medium. The results reproduce the known far-field asymptotics containing longitudinal, transversal and surface Rayleigh-type wave modes. The near-field asymptotic is essentially different for the regions far and near the boundary.     

Transformation of Gorkov's equation for type II superconductors into transport-like equations

The stationary behavior of type II superconductors is completely described by Gorkov's equations for a set of four Green's functions, supplemented by two self-consistency equations for gap parameterΔ(r) and vector potentialA(r). Expanding all quantities as usual at the Fermi surface and averaging over impurity positions, this set of equations is transformed into a simpler set for integrated Green's functions (which still contain much more information than is needed in most cases). The resulting equations, when linearized, yield essentially Lüders' transport equation for de Gennes' correlation function. The full equations contain all the known results and provide a promising starting point for numerical calculations. The thermodynamic potential is constructed as a functional of the integrated Green's functions and the mean fieldsΔ andA and a variational principle is formulated which uses this functional. Finally, paramagnetic scatterers are included (in Born approximation) as an example for possible generalizations of the new equations.     

Homogenization of multicoated inclusion-reinforced linear elastic composites with eigenstrains: Application to thermoelastic behavior

A new micromechanical approach for arbitrary multicoated ellipsoidal elastic inclusions with general eigenstrains is developed. We start from the integral equation of the linear elastic medium with eigenstrains adopting the Green's function technique and we apply an ‘(n+1)-phase’ model with a self-consistent condition to determine the homogenized behavior of multicoated inclusion-reinforced composites. The effective elastic moduli and eigenstrains are obtained as well as the residual stresses through the local stress concentration equations. The effective eigenstrains are determined either with the concentration tensors obtained here by the present model, or, more classically, with Levin's formula. In order to assess our micromechanical model, some applications to the isotropic thermoelastic behavior of composites with and without interphase are given. In particular, ‘four-phase’ and ‘three-phase’ models are derived for isotropic homothetic spherical inclusion-reinforced materials, and the results are successfully compared to exact analytical solutions regarding the effective elastic moduli and the effective thermal expansion.