格林函数倒易性在电磁场理论中的应用

格林函数在电磁场理论中有着非常广泛的应用.如在求解静电场问题时,往往会涉及到求解感应电荷的问题,而一般来说感应电荷的量值是不易求得的,特别是对不规则形状的导体.通过应用格林函数的倒易性来求解某些接地导体上感应电荷,能比较简便地解决这个问题.

零温下在Λ线三角绝缘铁磁体的磁振子软化

在二维绝缘铁磁系统基础上建立了一个磁振子-声子相互作用模型.利用Matsubara格林函数方法研究了磁振子-声子相互作用下的三角绝缘铁磁体的磁振子谱,计算了布里渊区的主要对称点线上的磁振子色散曲线.发现在布里渊区边界区域磁振子谱的软化和磁振子谱线增宽最明显.比较了纵向声子与横向声子对磁振子谱的软化与磁振子谱线增宽的影响,也讨论了各项参数的变化对磁振子谱的软化与磁振子谱线增宽的影响.

Singularity-free Green's function for EM sources embedded in a stratified medium

We present a method to unify the calculation of Green＇s functions for an electromagnetic（EM） transmitting source embedded in a homogeneous stratified medium.A virtual interface parallel to layer interfaces is introduced through the source location.The potentials for Green＇s function are derived by decomposing the partial wave solutions to Helmholtz＇s equations into upward and downward within boundaries.The amplitudes of the potentials in each stratum are obtained recursively from the initial amplitudes at the source level.The initial amplitudes are derived by coupling with the transmitting sources and following the discontinuity of the tangential electric and magnetic fields at the source interface.Only the initial terms are related to the transmitting sources and thus need to be modified for different transmitters,whereas the kernel connected with the stratified media stays unchanged.Hence,the present method can be easily applied to EM transmitting sources with little modification.The application of the proposed method to the marine controlled-source electromagnetic method（MCSEM） demonstrates its simplicity and flexibility.

The research on the highly efficient calculation method of 3-D frequency-domain Green function

The traditional calculation method of frequency-domain Green function mainly utilizes series or asymptotic expansion to carry out numerical approximation, however, this method requires very careful zoning, thus the computing process is complex with many cycles, which has greatly affected the computing efficiency. To improve the computing efficiency, this paper introduces Gaussian integral to the numerical calculation of the frequency-domain Green function and its partial derivatives. It then compares the calculation result with that in existing references. The comparison results demonstrate that, on the basis of its sufficient accuracy, the method has greatly simplified the computing process, reduced the zoning and improved the computing efficiency.

Transport Properties of Two Coupled Quantum Dots Under Optical Pumping

Using the Keldysh-Green function,we present a theoretical study on the electron transport properties of two coupled quantum dots under optical pumping. Plateaus in the I-V curve and resonant peaks in the transmission coefficient occur and can be explained by the local electron density of states in the quantum dots. The effects of the optical pumping frequency and intensity on the transport properties of the system are also discussed. The electron dynamical localization phenomenon occurs when the optical pumping frequency is equal to the discrete hole energy level. This result can be used to realize optical control switches.

Analysis of a Cavity Used for Medical Purposes with the Dyadic Green's Functions and the Moment method

The structure of a microwave radiator used for medical purposes is described. The dyadic Green's function and the method are used to analyze this Kind of multimode rectangular medium-filled cavity. The distribution of electromagnetic field intensity and the power density,as well as the temperature effect in the biological sample load are obtained.OPtimization of the size of cavity and the position of the input aperture have been performed with the computer to optimize the uniformity or microwave effect and the input VSWR.Necessary experiments were performed to compare the data obtained with theoretical analysis.

THE BOUNDARY ELEMENT ANALYSIS PROGRAM FOR 3-D TRANSIENT DYNAMIC FIELD

The matrix expression for the 3 D transient dynamic boundary integral equation in Laplace transform space is obtained and the degenerative element method has been implemented to treat the kernel function over the singular element. In the computer program BEMTDY the Koizumi′s numerical inversion method is used and three examples of the 3 D vibrated foundation under harmonic forces and the influence with both adjacent foundations are studied.

P-FFT and FG-FFT with real coefficients algorithm for the EFIE

In order to reduce the storage amount for the sparse coefficient matrix in pre-corrected fast Fourier transform （P-FFT） or fitting the Green function fast Fourier transform （FG-FFT）, the real coefficients are solved by improving the solution method of the coefficient equations. The novel method in both P-FFT and FG-FFT for the electric field integral equation （EFIE） is employed. With the proposed method, the storage amount for the sparse coefficient matrix can be reduced to the same level as that in the adaptive integral method （AIM） or the integral equation fast Fourier transform （IE-FFT）. Meanwhile, the new algorithms do not increase the number of the FFTs used in a matrix-vector product, and maintain almost the same level of accuracy as the original versions. Besides, in respect of the time cost in each iteration, the new algorithms have also the same level as AIM （or IE- FFF）. The numerical examples demonstrate the advantages of the proposed method.

2.5D scattering of incident plane SH waves by a canyon in layered half-space

This paper presents 2.5D scattering of incident plane SH waves by a canyon in layered half-space by the indirect boundary element method (IBEM). The free field response is carried out to give the displacements and stresses on the line which forms boundary of the canyon. The fictitious uniform moving loads are applied to the same line to calculate the Green's functions for the displacements and stresses. The amplitudes of the loads are determined by the boundary conditions. The displacements due to the free field and from the fictitious uniform moving loads have to be added to obtain the whole motion. The numerical results are carried out for the cases of a canyon in homogenous and in one layer over bedrock. The results show that the 2.5D wave scattering problem is essentially different from the 2D case, and there exist distinct differences between the wave amplification by a canyon in layered half-space and that in homogeneous half-space. The reasons for the distinct difference are explored, and the effects of the thickness and stiffness of the layer on the amplification are discussed.

Decomposition of the Unsteady Wave Patterns for Bessho form Translating-Pulsating Source Green Function

In order to interpret the physical feature of Bessho form translating-pulsating source Green function, the phase function is extracted from the integral representation and stationary-phase analysis is carried out in this paper. The complex characteristics of the integral variable and segmentation of the integral intervals are discussed in m complex plane. In θ space, the interval [-π/2+φ,-π/2+φ-iε] is dominant in the near-field flow, and there is a one-to-one correspondence between the real intervals in m space and the unsteady wave patterns in far field. If 4τ>1(τ is the Brard number), there are three kinds of propagation wave patterns such as ring-fan wave pattern, fan wave pattern and inner V wave pattern, and if 0<4τ<1, a ring wave pattern, an outer V and inner V wave pattern are presented in far field. The ring-fan or ring wave pattern corresponds to the interval [-π+α,-π/2+φ] for integral terms about k2, and the fan or outer V wave pattern and inner V wave pattern correspond to [-π+α,-π/2) and(-π/2,-π/2+φ] respectively for terms about k1. Numerical result shows that it is beneficial to decompose the unsteady wave patterns under the condition of τ≠0 by converting the integral variable θ to m. In addition, the constant-phase curve equations are derived when the source is performing only pulsating or translating.

A Free Surface Frequency Domain Green Function with Viscous Dissipation and Partial Reflections from Side Walls

The free-surface Green function method is widely used in solving the radiation or diffraction problems caused by a ship or ocean structure oscillating on the waves. In the context of inviscid potential flow, hydrodynamic problems such as multi-body interaction and tank side wall effect cannot be properly dealt with based on the traditional free-surface frequency domain Green function method, in which the water viscosity is omitted and the energy dissipation effect is absent. In this paper, an open-sea Green function with viscous dissipation was presented within the theory ofvisco-potential flow. Then the tank Green function with a partial reflection from the side walls in wave tanks was formulated as a formal sum of open-sea Green functions representing the infinite images between two parallel side walls of the source in the tank. The new far-field characteristics of the tank Green function is vitally important fur improving the validity of side-wall effects evaluation, which can be used in supervising the tank model tests.

Study of algorithms of new slender ship theory of wave resistance

In this paper, Noblesse's New Slender-Ship Wave-Making Theory was investigated numerically. Detailed expressions of zeroth and lst order wave resistance have been derived and calculation programs have also been compiled. In the single and double integral terms of Green function, the kernel function of wave resistance expression, special function expansion method and Chebyshev polynomials approach have been adopted respectively, which greatly simplify the calculation and increase the convergence speed.

The Interaction of Oblique Flexural Gravity Waves With a Small Bottom Deformation on a Porous Ocean-Bed： Green＇s Function Approach

The interaction of oblique incident water waves with a small bottom deformation on a porous ocean-bed is examined analytically here within the framework of linear water wave theory. The upper surface of the ocean is assumed to be covered by an infinitely extended thin uniform elastic plate, while the lower surface is bounded by a porous bottom surface having a small deformation. By employing a simplified perturbation analysis, involving a small parameter c^（〈〈l ）, which measures the smallness of the deformation, the governing Boundary Value Problem （BVP） is reduced to a simpler BVP for the first-order correction of the potential function. This BVP is solved using a method based on Green＇s integral theorem with the introduction of suitable Green＇s function to obtain the first-order potential, and this potential function is then utilized to calculate the first-order reflection and transmission coefficients in terms of integrals involving the shape function c（x） representing the bottom deformation. Consideration of a patch of sinusoidal ripples shows that when the quotient of twice the component of the incident field wave number propagating just below the elastic plate and the ripple wave number approaches one, the theory predicts a resonant interaction between the bed and the surface below the elastic plate. Again, for small angles of incidence, the reflected wave energy is more as compared to the other angles of incidence. It is also observed that the reflected wave energy is somewhat sensitive to the changes in the flexural rigidity of the elastic plate, the porosity of the bed and the ripple wave numbers. The main advantage of the present study is that the results for the values of reflection and transmission coefficients obtained are found to satisfy the energy-balance relation almost accurately.

Fast Prediction of Acoustic Radiation from a Hemi-capped Cylindrical Shell in Waveguide

In order to predict acoustic radiation from a structure in waveguide, a method based on wave superposition is proposed, in which the free-space Green＇s function is used to match the strength of equivalent sources. In addition, in order to neglect the effect of sound reflection from boundaries, necessary treatment is conducted, which makes the method more efficient. Moreover, this method is combined with the sound propagation algorithms to predict the sound radiated from a cylindrical shell in waveguide. Numerical simulations show the effect of how reflections can be neglected if the distance between the structure and the boundary exceeds the maximum linear dimension of the structure. It also shows that the reflection from the bottom of the waveguide can be approximated by plane wave conditionally. The proposed method is more robust and efficient in computation, which can be used to predict the acoustic radiation in waveguide.

Heat Generation by Electric Current in Normal-Metal-Molecular Quantum Dot-Superconductor System

We investigate the heat generation induced by electrical current in a normal-metal-molecular quantum dot-superconductor （NDS） system. By using nonequilibrium Green＇s function method, the heat generation Q is derived and studied in detail. The superconducting lead influences the heat generation significantly. An obvious step appears in Q - eV characteristics and the iocation of this step is related with the phonon frequency ωo. The heat generations exhibit very different behaviour in the condition eV 〈 △ and eV 〉 △ due to different tunneling mechanism. From the study of Q - eVg curves, there is an extra peak as eV 〉 △. The difference in this two cases is also shown in Q - ωo curve, an extra peak emerges as eV 〉 △.

Solution of pavement temperature field in“Environment-Surface”system through Green's function

In order to simplify the boundary conditions of pavement temperature field,the "Environment-Surface" system which considered the natural environment and pavement surface was established.Based on this system,the partial differential equations of the one-dimensional heat conduction in the pavement were established on the basis of the heat transfer theory.Furthermore,the function forms of the initial and boundary conditions of the equations were created through the field experiments.The general solution of the pavement one-dimensional heat conduction partial differential equations was acquired by using Green's function,and the explicit expression of pavement temperature field under specific constraint conditions was derived.For the purpose of analysis,the pavement temperatures in different seasons were calculated using the explicit expression of pavement temperature field,and the calculation accuracy was analyzed through the comparison between measured and calculated values.Then,the relationship between fitting accuracy and calculation accuracy of pavement temperatures was analyzed.The analysis results show that: the usage of "Environment-Surface" system simplifies the calculation of pavement temperature field; the relative error between calculated and measured values is generally less than 7% and is seldom influenced by seasons; there is a positive correlation between the calculation accuracy and the fitting accuracy of pavement surface temperature; high fitting accuracy would result in less error of pavement temperature prediction.

Characteristics of unsteady flow in porous media while considering threshold pressure gradient with Green＇s function

The flow behavior in porous media with threshold pressure gradient(TPG) is more complex than Darcy flow and the equations of motion, and outer boundary and inner boundary with TPG are also different from Darcy flow for unsteady flow of a producing well in a reservoir. An analytic method to solve this kind of problem is in a need of reestablishment. The classical method of Green's function and Newman product principle in a new way are used to solve the unsteady state flow problems of various shapes of well and reservoir while considering the TPG. Four Green's functions of point, line, band and circle while considering the TPG are achieved. Then, two well models of vertical well and horizontal well are built and simultaneously the function to calculate the moving boundary of each well model is provided. The results show that when considering TPG the pressure field is much different, which has a sudden pressure change, with a moving boundary in it. And the moving boundary of each well model increases with time but slows down rapidly, especially when the TGP is large.

Theoretical Investigation of Exchange Bias in Compensated Cases

The exchange bias （EB） of the ferromagnetic （FM）/antiferromagnetic （AFM） bilayers in a compensated case is studied by use of the many-body Green＇s function method of quantum statistical theory. The so-called compensated case is that there is no net magnetization on the AFM side of the interface. Our conclusion is that the EB in this case is primarily from the asymmetry of the interracial exchange coupling strengths between the FM and the two sublattices of the AFM. The effects of the layer thickness, temperature and the interracial coupling strength oi2 the exchange bias HE are investigated. The dependence of HE on the FM layer thickness and temperature is qualitatively in agreement with experimental results. HE is nearly inversely proportional to FM thickness. When temperature varies, both HE and He decrease with temperature increasing. The anisotropy of the FM layer only slightly influence He, but does not influence HE.

Perturbative Correction to Transverse Ward-Takahashi Relation for the Vector Vertex

We re-derive exactly the transverse Ward–Takahashi relation for the vector vertex in momentum space. The result shows that this transverse Ward–Takahashi relation in momentum space involves a perturbative correction term. We demonstrate explicitly that this transverse Ward–Takahashi relation is satisfied indeed at one-loop order.

Solving Electromagnetism Differential Equations Based on Random Walk

An approach of solving the finite difference equations with Monte Carlo method was presented. The accuracy of algorithm is guaranteed by the Central Limit Theorem. The computation of the value on each single node is independent of each other, which makes it easy to realize the algorithm in parallel processing and greatly improves the efficiency while dealing with local area computation. As long as different kinds of boundary conditions are statistics modeled,wide applications can the be made where the finite difference method is of competence.