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.

An improved empirical green function method to extract microseismic signals

In mining, reservoir impoundment and hydraulic fracturing exploitation, the induced microseismic clusters may present complex waveforms at receivers since the individual events generally arrive very close. And the routine methods for arrival picking are insufficiently efficient due to these strong influences. Here, we modified the empirical green function method and applied it to extract the green functions and the radiation coefficients of microseismic events with the largest energy in clusters.Multiple-channel records were used to estimate an "average" source spectrum and then the "average" source wavelet was removed from the records by deconvolution. We applied this method to the real data,and the result indicated clear improvement in extracting the dominant event of the clusters.

Source time functions of the Gonghe，China earthquake retrieved from long－period digital waveform data using empirical Green＇s function technique

An earthquake of Ms= 6, 9 occurred at the Gonghe, Qinghai Province, China on April 26, 1990. Three larger aftershocks took place at the same region, Ms= 5. 0 on May 7, 1990, Ms= 6. 0 on Jan. 3, 1994 and Ms= 5. 7on Feb. 16, 1994. The long-period recordings of the main shock from China Digital Seismograph Network (CDSN) are deconvolved for the source time functions by the correspondent0 recordings of the three aftershocks asempirical Green's functions (EGFs). No matter which aftershock is taken as EGF, the relative source time functions (RSTFs) Obtained are nearly identical. The RSTFs suggest the Ms= 6. 9 event consists of at least two subevents with approximately equal size whose occurrence times are about 30 s apart, the first one has a duration of 12 s and a rise time of about 5 s, and the second one has a duration of 17 s and a rise time of about & s. COmParing the RSTFs obtained from P- and SH-phases respectively, we notice that those from SH-phases are a slightly more complex than those from p-phases, implying other finer subevents exist during the process of the main shock. It is interesting that the results from the EGF deconvolution of long-Period way form data are in good agreement with the results from the moment tensor inversion and from the EGF deconvolution of broadband waveform data. Additionally, the two larger aftershocks are deconvolved for their RSTFs. The deconvolution results show that the processes of the Ms= 6. 0 event on Jan. 3, 1994 and the Ms= 5. 7 event on Feb. 16,1994 are quite simple, both RSTFs are single impulses.The RSTFs of the Ms= 6. 9 main shock obtained from different stations are noticed to be azimuthally dependent, whose shapes are a slightly different with different stations. However, the RSTFs of the two smaller aftershocks are not azimuthally dependent. The integrations of RSTFs over the processes are quite close to each other, i. e., the scalar seismic moments estimated from different stations are in good agreement. Finally the scalar seismic moments of the three aftershocks are compared. The relative scalar seismic moment Of the three aftershocks deduced from the relative scalar seismic moments of the Ms=6. 9 main shock are very close to those inverted directly from the EGF deconvolution. The relative scalar seismic moment of the Ms =6. 9 main shock calculated using the three aftershocks as EGF are 22 (the Ms= 6. 0 aftershock being EGF), 26 (the Ms= 5. 7 aftershock being EGF) and 66 (the Ms= 5. 5 aftershock being EGF), respectively. Deducingfrom those results, the relative scalar sesimic moments of the Ms= 6. 0 to the Ms= 5. 7 events, the Ms= 6. 0 tothe Ms= 5. 5 events and the Ms= 5. 7 to the Ms= 5. 5 events are 1. 18, 3. 00 and 2. 54, respectively. The correspondent relative scalar seismic moments calculated directly from the waveform recordings are 1. 15, 3. 43, and 3. 05.

Rupture process of November 6,1988,Lan-cang-Gengma,Yunnan,China,earthquake of M_s=7.6 using empirical Green’s function de-convolution method

The empirical Green's function deconvolution technique is applied to the November 6, 1988, Ms=7.6, Yunnan, China earthquake to study the rupture process of this event. An aftershock of Ms=6.3 is taken as the empirical Green's function. Recordings of these two events at six stations of China Digital Seismographic Network (CDSN) are employed. Deconvolution results show that this event is a relatively simple event. Directivity analysis indicates that the rupture was initiated at hypocenter and propagated bilaterally and nearly symmetrically towards the northwest and southeast directions with a total length of 70 km and a time duration of 19 s. The rupture velocity is estimated to be about 2.0 km/s.

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.

A NEW ONE-STEP BAND-LIMITED EXTRAPOLATION PROCEDURE USING EMPIRICAL ORTHOGONAL FUNCTIONS

A one-step band-limited extrapolation procedure is systematically developed under an a priori assumption of bandwidth. The rationale of the proposed scheme is to expand the known signal segment based on a band-limited basis function set and then to generate a set of Empirical Orthogonal Functions (EOF’s) adaptively from the sample values of the band-limited function set. Simulation results indicate that, in addi- tion to the attractive adaptive feature, this scheme also appears to guarantee a smooth result for inexact data, thus suggesting the robustness of the proposed procedure.

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 Function模型

将纳米分子器件的电极和分子本体作为一个系统,考虑到入射电子在原子电极和分子本体中的传导,在电极与分子之间界面上的透射的量子状态,在满足Fisher-Lee关系式的基础上,一个新型的Green’sFunction被设计出来。这一理论模型把三个区域上的不同量子状态完整地表达在一个方程式里,符合物理学基本原理,可以系统地用来研究纳米分子器件的量子传导特性。

Synthesis of accelerations of the 1976 Tangshan earthquake(M_S 7.8)in near-and far-field by using semi-empirical method

In this paper a heterogeneous fault model of the Tangshan earthquake is suggested, which consists of two southern sub faults striking N30°E and two northern sub faults striking N50°E. Total length of the main shock fault is 114 km and seismic moment is about 1.4×10 20 N·m. The epicentre of the main shock is located at the southern part, near the intersection of the two bands. Accelerations of two aftershocks ( M L 5.5, M S 6.9 ) were used as empirical Green's functions to synthesize the accelerations of the main shock in near and far field. A method that small events and main shock are considered not satisfying the similarity relationship in the improved empirical Green's function is also applied in this paper. Peak values, duration and response spectra of synthesized accelerations in far field are in agreement with the observed records. The synthesized results in near field are also in agreement with the epicentral intensity distribution of the main shock. The results show that the peak acceleration of Tangshan earthquake in epicentral region exceeds 1.1 g. It is consistent with the peak accelerations recorded in some large earthquakes occurred in recent years.

基于小震经验格林函数研究中等地震破裂方向性:以2008年攀枝花M_(S)6.1地震为例

中小地震破裂方向性研究在确定实际发震面以及地震灾害评估中具有非常重要的意义.本文以2008年8月30日攀枝花M_(S)6.1地震为研究实例,以较小余震作为经验格林函数地震,利用表观震源时间函数随台站方位角变化特性,确定了主震的破裂方向、破裂尺度以及破裂速度,验证了该方法的适用性和可行性.研究结果表明,主震的破裂方向为几乎正南的184°,破裂长度8.05 km,破裂速度3.35 km·s^(-1),进一步确认了主震震源机制解其中沿南北向的节面是地震实际发震面,主震在水平方向沿着南北向元谋断裂从北往南破裂.攀枝花地震震源破裂速度接近地壳剪切波速度,这可能是造成此次地震虽然震级不大但破坏性严重的重要原因.本文综合分析主震的矩心深度、震源机制解、破裂方向性以及构造应力背景等信息,探讨了地震发生的破裂过程以及灾害严重的原因.

Study on Characteristics of 3-D Translating-Pulsating Source Green Function of Deep-Water Havelock Form and Its Fast Integration Method

The singularities, oscillatory performances and the contributing factors to the 3-＇D translating-pulsating source Green function of deep-water Havelock form which consists of a local disturbance part and a far-field wave-like part, are analyzed systematically. Relative numerical integral methods about the two parts are presented in this paper. An improved method based on LOBATTO rule is used to eliminate singularities caused respectively by infinite discontinuity and jump discontinuous node from the local disturbance part function, which makes the improvement of calculation efficiency and accuracy possible. And variable substitution is applied to remove the singularity existing at the end of the integral interval of the far-field wave-like part function. Two auxiliary techniques such as valid interval calculation and local refinement of integral steps technique in narrow zones near false singularities are applied so as to avoid unnecessary integration of invalid interval and improve integral accordance. Numerical test results have proved the efficiency and accuracy in these integral methods that thus can be applied to calculate hydrodynamic performance of floating structures moving in waves.

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.

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.