Implementation of the Integrated Green’s Function Method for 3D Poisson’s Equation in a Large Aspect Ratio Computational Domain

The solution of Poisson’s Equation plays an important role in many areas, including modeling high-intensity and high-brightness beams in particle accelerators. For the computational domain with a large aspect ratio, the integrated Green’s function method has been adopted to solve the 3D Poisson equation subject to open boundary conditions. In this paper, we report on the efficient implementation of this method, which can save more than a factor of 50 computing time compared with the direct brute force implementation and its improvement under certain extreme conditions.

Existence of Monotone Positive Solution for a Fourth-Order Three-Point BVP with Sign-Changing Green’s Function

This paper is concerned with the following fourth-order three-point boundary value problem , where , we discuss the existence of positive solutions to the above problem by applying to the fixed point theory in cones and iterative technique.

Green's functions of two-dimensional piezoelectric quasicrystal in half-space and bimaterials

In this paper,we obtain Green’s functions of two-dimensional(2D)piezoelectric quasicrystal(PQC)in half-space and bimaterials.Based on the elastic theory of QCs,the Stroh formalism is used to derive the general solutions of displacements and stresses.Then,we obtain the analytical solutions of half-space and bimaterial Green’s functions.Besides,the interfacial Green’s function for bimaterials is also obtained in the analytical form.Before numerical studies,a comparative study is carried out to validate the present solutions.Typical numerical examples are performed to investigate the effects of multi-physics loadings such as the line force,the line dislocation,the line charge,and the phason line force.As a result,the coupling effect among the phonon field,the phason field,and the electric field is prominent,and the butterfly-shaped contours are characteristic in 2D PQCs.In addition,the changes of material parameters cause variations in physical quantities to a certain degree.

Spectral-domain dyadic Green’s functions in chiral media

A novel method of deriving the electromagnetic dyadic Green＇s functions in an unbounded, lossless, reciprocal and homogeneous chiral media described by the constitutive relations D = εE ＋ jγB and H = jγE ＋ μ^-1B - （ωε）^-1γJ is given. The divergenceless and irrotational splitting of dyadic Dirac 8 function and Fourier transformation are used to directly obtain the divergenceless and irrotational component of spectral-domain dyadic Green＇s functions in chiral media. This method avoids using the dyadic Green＇s function eigenfunction expansion technique. The method given here can be generalized to a source-free region and an achiral case.

Continuum Skyrme Hartree-Fock-Bogoliubov theory with Green’s function method for neutron-rich Ca,Ni,Zr,and Sn isotopes

The possible exotic nuclear properties in the neutron-rich Ca,Ni,Zr,and Sn isotopes are examined with the continuum Skyrme Hartree-Fock-Bogoliubov theory in the framework of the Green’s function method.The pairing correlation,the couplings with the continuum,and the blocking effects for the unpaired nucleon in odd-A nuclei are properly treated.The Skyrme interaction SLy4 is adopted for the ph channel and the density-dependentinteraction is adopted for the pp chan-nel,which well reproduce the experimental two-neutron separation energies S_(2n)and one-neutron separation energies Sn.It is found that the criterion S_(n)>0 predicts a neutron drip line with neutron numbers much smaller than those for S_(2n)>0.Owing to the unpaired odd neutron,the neutron pairing energies−E_(pair)in odd-A nuclei are much lower than those in the neighbor-ing even-even nuclei.By investigating the single-particle structures,the possible halo structures in the neutron-rich Ca,Ni,and Sn isotopes are predicted,where sharp increases in the root-mean-square(rms)radii with significant deviations from the traditional rA^(1∕3)rule and diffuse spatial density distributions are observed.Analyzing the contributions of various partial waves to the total neutron densityρlj(r)∕ρ(r)reveals that the orbitals located around the Fermi surface-particularly those with small angular momenta-significantly affect the extended nuclear density and large rms radii.The number of neutrons Nλ(N_(0))occupying above the Fermi surfacen(continuum threshold)is discussed,whose evolution as a function of the mass number A in each isotope is consistent with that of the pairing energy,supporting the key role of the pairing correlation in halo phenomena.

GREEN'S FUNCTION AND THE POINTWISE BEHAVIORS OF THE VLASOV-POISSON-FOKKER-PLANCK SYSTEM

The pointwise space-time behaviors of the Green’s function and the global solution to the Vlasov-Poisson-Fokker-Planck(VPFP)system in three dimensional space are studied in this paper.It is shown that the Green’s function consists of the diffusion waves decaying exponentially in time but algebraically in space,and the singular kinetic waves which become smooth for all(t,x,v)when t>0.Furthermore,we establish the pointwise space-time behaviors of the global solution to the nonlinear VPFP system when the initial data is not necessarily smooth in terms of the Green’s function.

Green's functions for uniformly distributed loads acting on an inclined line in a poroelastic layered site

Based on one type of practical Biot＇s equation and the dynamic-stiffness matrices ofa poroelastic soil layer and half-space, Green＇s functions were derived for unitformly distributed loads acting on an inclined line in a poroelastie layered site. This analysis overcomes significant problems in wave scattering due to local soil conditions and dynamic soil-structure interaction. The Green＇s functions can be reduced to the case of an elastic layered site developed by Wolf in 1985. Parametric studies are then carried out through two example problems.

In-plane dynamic Green's functions for inclined and uniformly distributed loads in a multi-layered transversely isotropic half-space

The dynamic stiffness method combined with the Fourier transform is utilized to derive the in-plane Green’s functions for inclined and uniformly distributed loads in a multi-layered transversely isotropic（TI）half-space.The loaded layer is fixed to obtain solutions restricted in it and the corresponding reactions forces,which are then applied to the total system with the opposite sign.By adding solutions restricted in the loaded layer to solutions from the reaction forces,the global solutions in the wavenumber domain are obtained,and the dynamic Green’s functions in the space domain are recovered by the inverse Fourier transform.The presented formulations can be reduced to the isotropic case developed by Wolf（1985）,and are further verified by comparisons with existing solutions in a uniform isotropic as well as a layered TI halfspace subjected to horizontally distributed loads which are special cases of the more general problem addressed.The deduced Green’s functions,in conjunction with boundary element methods,will lead to significant advances in the investigation of a variety of wave scattering,wave radiation and soil-structure interaction problems in a layered TI site.Selected numerical results are given to investigate the influence of material anisotropy,frequency of excitation,inclination angle and layered on the responses of displacement and stress,and some conclusions are drawn.

TIME-HARMONIC DYNAMIC GREEN'S FUNCTIONS FOR ONE-DIMENSIONAL HEXAGONAL QUASICRYSTALS

Quasicrystals have additional phason degrees of freedom not found in conventional crystals. In this paper, we present an exact solution for time-harmonic dynamic Green＇s function of one-dimensional hexagonal quasicrystals with the Laue classes 6/mh and 6/mhmm. Through the introduction of two new functions φ and ψ, the original problem is reduced to the determination of Green＇s functions for two independent Helmholtz equations. The explicit expressions of displacement and stress fields are presented and their asymptotic behaviors are discussed. The static Green＇s function can be obtained by letting the circular frequency approach zero.

Decoupling coefficients of dilatational wave for Biot's dynamic equation and its Green's functions in frequency domain

Green＇s functions for Blot＇s dynamic equation in the frequency domain can be a highly useful tool for the investigation of dynamic responses of a saturated porous medium. Its applications are found in soil dynamics, seismology, earthquake engineering, rock mechanics, geophysics, and acoustics. However, the mathematical work for deriving it can be daunting. Green＇s functions have been presented utilizing an analogy between the dynamic thermoelasticity and the dynamic poroelasticity in the frequency domain using the u-p formulation. In this work, a special term ＂decoupling coefficient＂ for the decomposition of the fast and slow dilatational waves is proposed and expressed to present a new methodology for deriving the poroelastodynamic Green＇s functions. The correct- ness of the solution is demonstrated by numerically comparing the current solution with Cheng＇s previous solution. The separation of the two waves in the present methodology allows the more accurate evaluation of Green＇s functions, particularly the solution of the slow dilatational wave. This can be advantageous for the numerical implementation of the boundary element method （BEM） and other applications.

Comprehensive form of solution for Lamb's dynamic problem expressed by Green's functions

By the analysis for the vectors of a wave field in the cylindrical coordinate and Sommerfeld＇s identity as well as Green＇s functions of Stokes＇ solution pertaining the conventional elastic dynamic equation, the results of Green＇s function in an infinite space of an axisymmetric coordinate are shown in this paper. After employing a supplementary influence field and the boundary conditions in the free surface of a senti-space, the authors obtain the solutions of Green＇s function for Lamb＇s dynamic problem. Besides, the vertical displacement uzz and the radial displacement urz can match Lamb＇s previous results, and the solutions of the linear expansion source u^r and the linear torsional source uee are also given in the paper. The authors reveal that Green＇s function of Stokes＇ solution in the semi-space is a comprehensive form of solution expressing the dynamic Lamb＇s problem for various situations. It may benefit the investigation of deepening and development of Lamb＇s problems and solution for pertinent dynamic problems conveniently.

Study on the 3D Green＇s functions of the fluid and piezoelectric bimaterials

Because most piezoelectric devices have interfaces with fluid in engineering, it is valuable to study the coupled field between fluid and piezoelectric media. As the fundamental problem, the 3D Green＇s functions for point forces and point charge loaded in the fluid and piezoelectric bimaterials are studied in this paper. Based on the 3D general solutions expressed by harmonic functions, we constructed the suitable harmonic functions with undetermined constants at first. Then, the couple field in the fluid and piezoelectric bimaterials can be derived by substitution of harmonic functions into general solutions. These constants can be obtained by virtue of the compatibility, boundary, and equilibrium conditions. At last, the characteristics of the electromechanical coupled fields are shown by numerical results.

PERIODIC DYADIC GREEN’S FUNCTION FOR FIELD ANALYSIS IN EMC CHAMBER

Herein a novel Dyadic Green’s Function (DGF) is presented to calculate the field in ElectroMag-netic Compatibility (EMC) chamber. Due to the difficulty of simulating the whole chamber environment, the analysis combines the DGF formulation and the FEM method, with the latter deals with the reflection from absorbers. With DGF formulation for infinite periodic array structures, this paper investigates electromagnetic field in chamber with truncated arrays. The reflection from the absorber serves as the virtual source contribut-ing to the total field. Hence the whole chamber field calculation can be separated from the work of absorber model set-up. Practically the field homogeneity test and Normal Site Attenuation (NSA) test are carried out to evaluate the chamber performance. Based on the method in this paper, the simulation results agree well with the test, and predict successfully the victim frequency points of the chamber.

A NEW METHOD FOR SOLVING DYADIC GREEN'S FUNCTION OF ELECTROMAGNETIC FIELD

A new method for solving electromagnetic field boundary value problem is given.Byusing this method,the boundary value problem of the vector wave equation can be transformedinto the independent boundary value problem of scalar wave equations and the two additionalvector differential operations.All the dyadic Green’s functions got by eigenfunction expansionof the dyadic Green’s function can be got by this method easily and some of the dyadic Green’sfunctions for complex systems which are very difficult to get by the ordinary method have beengot by this new method.The dyadic Green’s function for a dielectric loaded cavity is one of thegiven examples.

THE DYADIC GREEN'S FUNCTION IN A LOADED RECTANGULAR WAVEGUIDE AND ITS APPLICATION

This paper presents a method to derive the Dyadic Green’s Function(DGF)ofa loaded rectangular waveguide by using the image method.In the calculation of the DGF,we use the integral transformation and replace the multi-infinite summation by a single one;thus it greatly simplifies the calculation and saves computer time.As an example of the DGF’sapplication,we give the moment method’s scattering field calculation of a metal sphere resting onthe broad wall of the loaded rectangular waveguide.Results of our calculations well agree withboth data of experiments performed in our laboratory and those are published.It is easy to seethat the method used in this paper can be expanded to other related waveguide problems.

Spectral Theorem of Many-Body Green's Functions When Complex Eigenvalues Appear

In this paper, an extended spectral theorem is given, which enables one to calculate the correlation functions when complex eigenvalues appear. To do so, a Fourier transformation with a complex argument is utilized. We treat all the Matsbara frequencies, including Fermionic and Bosonic frequencies, on an equal footing. It is pointed out that when complex eigenvalues appear, the dissipation of a system cannot simply be ascribed to the pure imaginary part of the Green function. Therefore, the use of the name fluctuation-dissipation theorem should be careful.

Asymmetric Green's functions for exponentially graded transversely isotropic substrate–coating system

By virtue of a complete set of two displacement potentials,an analytical derivation of the elastostatic Green’s functions of an exponentially graded transversely isotropic substrate–coating system is presented.Three-dimensional point–load and patch–load Green’s functions for stresses and displacements are given in line-integral representations.The formulation includes a complete set of transformed stress–potential and displacement–potential relations,with utilizing Fourier series and Hankel transforms.As illustrations,the present Green’s functions are degenerated to the special cases such as an exponentially graded half-space and a homogeneous two-layered half-space Green’s functions.Because of complicated integrand functions,the integrals are evaluated numerically and for numerical computation of the integrals,a robust and effective methodology is laid out which gives the necessary account of the presence of singularities of integration.Comparisons of the existing numerical solutions for homogeneous two-layered isotropic and transversely isotropic half-spaces are made to confirm the accuracy of the present solutions.Some typical numerical examples are also given to show the general features of the exponentially graded two-layered half-space Green’s functions that the effect of degree of variation of material properties will be recognized.

Controllable precision of the projective truncation approximation for Green's functions

Recently, we developed the projective truncation approximation for the equation of motion of two-time Green's functions(Fan et al., Phys. Rev. B 97, 165140(2018)). In that approximation, the precision of results depends on the selection of operator basis. Here, for three successively larger operator bases, we calculate the local static averages and the impurity density of states of the single-band Anderson impurity model. The results converge systematically towards those of numerical renormalization group as the basis size is enlarged. We also propose a quantitative gauge of the truncation error within this method and demonstrate its usefulness using the Hubbard-I basis. We thus confirm that the projective truncation approximation is a method of controllable precision for quantum many-body systems.

Some Universal Properties of the Green’s Functions Associated with the Wave Equation in Bounded Partially-Homogeneous Domains and Their Use in Acoustic Tomography

Direct and inverse scattering problems connected with the wave equation in non-homogeneous bounded domains constitute challenging actual subjects for both mathematicians and engineers. Among them one can mention, for example, inverse source problems in seismology, nondestructive archeological probing, mine prospecting, inverse initial-value problems in acoustic tomography, etc. In spite of its crucial importance, almost all of the available rigorous investigations concern the case of unbounded simple domains such as layered planar or cylindrical or spherical structures. The main reason for the lack of the works related to non-homogeneous bounded structures is the extreme complexity of the explicit expressions of the Green’s functions. The aim of the present work consists in discovering some universal properties of the Green’s functions in question, which reduce enormously the difficulties arising in various applications. The universality mentioned here means that the properties are not depend on the geometrical and physical properties of the configuration. To this end one considers first the case when the domain is partially-homogeneous. Then the results are generalized to the most general case. To show the importance of the universal properties in question, they are applied to an inverse initial-value problem connected with photo-acoustic tomography.

Numerical Solution of Green’s Function for Solving Inhomogeneous Boundary Value Problems with Trigonometric Functions by New Technique

A numerical technique is presented for solving integration operator of Green’s function. The approach is based on Hermite trigonometric scaling function on [0,2π], which is constructed for Hermite interpolation. The operational matrices of derivative for trigonometric scaling function are presented and utilized to reduce the solution of the problem. One test problem is presented and errors plots show the efficiency of the proposed technique for the studied problem.