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.

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 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.

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.

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 Green's function of the harmonic horizontal force applied at the interior of the saturated half-space soil

By using integral transform methods, the Green(s functions of horizontal harmonic force applied at the interior of the saturated half-space soil are obtained in the paper. The general solutions of the Biot dynamic equations in frequency domain are established through the use of Hankel integral transforms technique. Utilizing the above- mentioned general solutions, and the boundary conditions of the surface of the half-space and the continuous conditions at the plane of the horizontal force, the solutions of the boundary value problem can be determined. By the numerical inverse Hankel transforms method, the Green(s functions of the harmonic horizontal force are obtainable. The degenerate case of the results deduced from this paper agrees well with the known results. Two numerical examples are given in the paper.

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.

Green's function of multi-layered poroelastic half-space for models of ground vibration due to railway traffic

This study proposes a Green＇s function, an essential representation of water-saturated ground under moving excitation, to simulate ground borne vibration from trains. First, general solutions to the governing equations of poroelastic medium are derived by means of integral transform. Secondly, the transmission and reflection matrix approach is used to formulate the relationship between displacement and stress of the stratified ground, which results in the matrix of the Green＇s function. Then the Green＇s function is combined into a train-track-ground model, and is verified by typical examples and a field test. Additional simulations show that the computed ground vibration attenuates faster in the immediate vicinity of the track than in the surrounding area. The wavelength of wheel-rail unevenness has a notable effect on computed displacement and pore pressure. The variation of vibration intensity with the depth of ground is significantly influenced by the layering of the strata soil. When the train speed is equal to the velocity of the Rayleigh wave, the Mach cone appears in the simulated wave field. The proposed Green＇s function is an appropriate representation for a layered ground with shallow ground water table, and will be helpful to understand the dynamic responses of the ground to complicated moving excitation.

NONLINEAR EVOLUTION SYSTEMS AND GREEN’S FUNCTION

In this paper, we will introduce how to apply Green＇s function method to get the pointwise estimates for the solutions of Cauchy problem of nonlinear evolution equations with dissipative structure. First of all, we introduce the pointwise estimates of the time-asymptotic shape of the solutions of the isentropic Navier-Stokes equations and show to exhibit the generalized Huygen＇s principle. Then, for other nonlinear dissipative evolution equations, we will only introduce the result and give some brief explanations. Our approach is based on the detailed analysis of the Green＇s function of the linearized system and micro-local analysis, such as frequency decomposition and so on.

Some Kind of Equations Involving Euler＇s Function

For any given positive integer n ≥ 1, the Euler function φ（n） is defined to be the number of positive integers not exceeding n which are relatively prime to n. w（n） is defined to be the number of different prime divisors of n. Some kind of equations involving Euler＇s function is studied in the paper.

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.

Searching for single-particle resonances with the Green’s function method

Single-particle resonances in the continuum are crucial for studies of exotic nuclei.In this study,the Green’s function approach is employed to search for single-particle resonances based on the relativistic-mean-field model.Taking^(120)Sn as an example,we identify singleparticle resonances and determine the energies and widths directly by probing the extrema of the Green’s functions.In contrast to the results found by exploring for the extremum of the density of states proposed in our recent study[Chin.Phys.C,44:084105(2020)],which has proven to be very successful,the same resonances as well as very close energies and widths are obtained.By comparing the Green’s functions plotted in different coordinate space sizes,we also found that the results very slightly depend on the space size.These findings demonstrate that the approach by exploring for the extremum of the Green’s function is also very reliable and effective for identifying resonant states,regardless of whether they are wide or narrow.

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.

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.

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.

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.

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.

ON GREEN'S FUNCTION FOR HYPERBOLIC-PARABOLIC SYSTEMS

We study the Green＇s function for a general hyperbolic-parabolic system, including the Navier-Stokes equations for compressible fluids and the equations for magnetohydrodynamics. More generally, we consider general systems under the basic Kawashima- Shizuta type of conditions. The first result is to make precise the secondary waves with subscale structure, revealing the nature of coupling of waves pertaining to different characteristic families. The second result is on the continuous differentiability of the Green＇s function with respect to a small parameter when the coefficients of the system are smooth functions of that parameter. The results significantly improve previous results obtained by the authors.

A new approach of solving Green's function for wave propagation in an inhomogeneous absorbing medium

A new approach is developed to solve the Green＇s function that satisfies the Hehmholtz equation with complex refractive index. Especially, the Green＇s function for the Helmholtz equation can be expressed in terms of a onedimensional integral, which can convert a Helmholtz equation into a Schrodinger equation with complex potential. And the Schrodinger equation can be solved by Feynman path integral. The result is in excellent agreement with the previous work.