The heavy quarks present in the quark-gluon plasma(QGP)can act as a probe of relativistic heavy ion collisions as they retain the memory of their interaction history.In a previous study,a stochastic Schrödinger e...The heavy quarks present in the quark-gluon plasma(QGP)can act as a probe of relativistic heavy ion collisions as they retain the memory of their interaction history.In a previous study,a stochastic Schrödinger equation(SSE)has been applied to describe the evolution of heavy quarks,where an external field with random phases is used to simulate the thermal medium.In this work,we study the connection between the SSE and the Boltzmann transport equation(BE)approach in the Keldysh Green’s function formalism.By comparing the Green’s function of the heavy quark from the SSE and the Keldysh Green’s functions leading to the Boltzmann equation,we demonstrate that the SSE is consistent with the Boltzmann equation in the weak coupling limit.We subsequently confirm their consistency through numerical calculations.展开更多
Understanding the effects of point liquid loading on transversely isotropic poroelastic media is crucial for advancing geomechanics and biomechanics, where precise modeling of fluid-structure interactions is essential...Understanding the effects of point liquid loading on transversely isotropic poroelastic media is crucial for advancing geomechanics and biomechanics, where precise modeling of fluid-structure interactions is essential. This paper presents a comprehensive analysis of infinite transversely isotropic poroelasticity under a fluid source, based on Biot's theory, aiming to uncover new and previously unexplored insights in the literature. We begin our study by deriving a general solution for fluid-saturated, transversely isotropic poroelastic materials in terms of harmonic functions that satisfy sixth-order homogeneous partial differential equations, using potential theory and Almansi's theorem. Based on these general solutions and potential functions, we construct a Green's function for a point fluid source, introducing three new harmonic functions with undetermined constants. These constants are determined by enforcing continuity and equilibrium conditions. Substituting these into the general solution yields fundamental solutions for poroelasticity that provide crucial support for a wide range of project problems. Numerical results and comparisons with existing literature are provided to illustrate physical mechanisms through contour plots. Our observations reveal that all components tend to zero in the far field and become singular at the concentrated source. Additionally, the contours exhibit rapid changes near the point fluid source but display gradual variations at a distance from it. These findings highlight the intricate behavior of the system under point liquid loading, offering valuable insights for further research and practical applications.展开更多
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 ...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.展开更多
The singularities and oscillatory performance of translating-pulsating source Green＇s function in Bessho form were analyzed. Relative numerical integration methods such as Gaussian quadrature rule, variable substitut...The singularities and oscillatory performance of translating-pulsating source Green＇s function in Bessho form were analyzed. Relative numerical integration methods such as Gaussian quadrature rule, variable substitution method （VSM）, and steepest descent integration method （SDIM） were used to evaluate this type of Green＇s function. For SDIM, the complex domain was restricted only on the 0-plane. Meanwhile, the integral along the real axis was computed by use of the VSM to avoid the complication of a numerical search of the steepest descent line. Furthermore, the steepest descent line was represented by the B-spline function. Based on this representation, a new self-compatible integration method corresponding to parametric t was established. The numerical method was validated through comparison with other existing results, and was shown to be efficient and reliable in the calculation of the velocity potentials for the 3D seakeeping and hydrodynamic performance of floating struc- tures moving in waves.展开更多
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 co...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.展开更多
In this paper, the Adomian decomposition method with Green’s function (Standard Adomian and Modified Technique) is applied to solve linear and nonlinear tenth-order boundary value problems with boundary conditions de...In this paper, the Adomian decomposition method with Green’s function (Standard Adomian and Modified Technique) is applied to solve linear and nonlinear tenth-order boundary value problems with boundary conditions defined at any order derivatives. The numerical results obtained with a small amount of computation are compared with the exact solutions to show the efficiency of the method. The results show that the decomposition method is of high accuracy, more convenient and efficient for solving high-order boundary value problems.展开更多
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-s...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.展开更多
The accurate analyses for a plate fin heat sink with the ability to control the temperature of the avionics devices within a pre-set controllable temperature range are required both in the process of circuit design an...The accurate analyses for a plate fin heat sink with the ability to control the temperature of the avionics devices within a pre-set controllable temperature range are required both in the process of circuit design and for the real-time temperature monitoring purposes. In order to provide an insight into the behavior of the temperature of a plate fin heat sink subjected non-uniform heat density on the surfaces, it is necessary to obtain accurate analytical solutions yielding explicit formulas relating the dissipated power to the temperature rise at any point of avionics devices. This paper presents a method for thermal simulation of a plate fin heat sink using an analytical solution of the three-dimensional heat equation resulting from an appropriate plate fin heat sink transient thermal model. The entire solution methodology is illustrated in detail on the particular examples of the plate fin heat sink subjected non-uniform heat density on the surfaces. The transient temperature profiles are obtained for different positions at the surface of the plate fin heat sink. The analytical results are compared with measurements made on the surface of the cold plate and it is found that they are in good agreement with an error of less than 3 K.展开更多
The principle objective of the paper is to study the acoustic radiation problem of the 3D space domain with control boundary. By using the conformal transformation theory, the Green's function for acoustic point s...The principle objective of the paper is to study the acoustic radiation problem of the 3D space domain with control boundary. By using the conformal transformation theory, the Green's function for acoustic point source in the control domain space is obtained. With it, the expression of acoustic radiation function of the control domain is formed. Discussion about the acoustic radiation of pulsating sphere in right-angle space is drawn in the end, including the acoustic radiation directivity effect by the boundary characteristics, acoustic radiation frequency and acoustic source location. Numerical results show that: for the lower frequency radiation, the infection of free surface is significant; for the high frequency radiation, the infection of location is significant on the contrary. The research provides a new method for boundary characteristic problem of the structural-acoustic acoustic.展开更多
To reduce computational costs, an improved form of the frequency domain boundary element method（BEM） is proposed for two-dimensional radiation and propagation acoustic problems in a subsonic uniform flow with arbitr...To reduce computational costs, an improved form of the frequency domain boundary element method（BEM） is proposed for two-dimensional radiation and propagation acoustic problems in a subsonic uniform flow with arbitrary orientation. The boundary integral equation（BIE） representation solves the two-dimensional convected Helmholtz equation（CHE） and its fundamental solution, which must satisfy a new Sommerfeld radiation condition（SRC） in the physical space. In order to facilitate conventional formulations, the variables of the advanced form are expressed only in terms of the acoustic pressure as well as its normal and tangential derivatives, and their multiplication operators are based on the convected Green＇s kernel and its modified derivative. The proposed approach significantly reduces the CPU times of classical computational codes for modeling acoustic domains with arbitrary mean flow. It is validated by a comparison with the analytical solutions for the sound radiation problems of monopole,dipole and quadrupole sources in the presence of a subsonic uniform flow with arbitrary orientation.展开更多
Using a polarization method, the scattering problem for a two-dimensional inclusion embedded in infinite piezoelectric/piezomagnetic matrices is investigated. To achieve the purpose, the polarization method for a two-...Using a polarization method, the scattering problem for a two-dimensional inclusion embedded in infinite piezoelectric/piezomagnetic matrices is investigated. To achieve the purpose, the polarization method for a two-dimensional piezoelectric/piezomagnetic ＂comparison body＂ is formulated. For simple harmonic motion, kernel of the polarization method reduces to a 2-D time-harmonic Green＇s function, which is obtained using the Radon transform. The expression is further simplified under conditions of low frequency of the incident wave and small diameter of the inclusion. Some analytical expressions are obtained. The analytical solutions for generalized piezoelectric/piezomagnetic anisotropic composites are given followed by simplified results for piezoelectric composites. Based on the latter results, two numerical results are provided for an elliptical cylindrical inclusion in a PZT-5H-matrix, showing the effect of different factors including size, shape, material properties, and piezoelectricity on the scattering cross-section.展开更多
Abstract： The scattering of plane SH-waves by topographic features in a layered transversely isotropic （TI） half-space is investigated by using an indirect boundary element method （IBEM）. Firstly, the anti-plane ...Abstract： The scattering of plane SH-waves by topographic features in a layered transversely isotropic （TI） half-space is investigated by using an indirect boundary element method （IBEM）. Firstly, the anti-plane dynamic stiffness matrix of the layered TI half-space is established and the free fields are solved by using the direct stiffness method. Then, Green＇s functions are derived for uniformly distributed loads acting on an inclined line in a layered TI half-space and the scattered fields are constructed with the deduced Green＇s functions. Finally, the free fields are added to the scattered ones to obtain the global dynamic responses. The method is verified by comparing results with the published isotropic ones. Both the steady-state and transient dynamic responses are evaluated and discussed. Numerical results in the frequency domain show that surface motions for the TI media can be significantly different from those for the isotropic case, which are strongly dependent on the anisotropy property, incident angle and incident frequency. Results in the time domain show that the material anisotropy has important effects on the maximum duration and maximum amplitudes of the time histories.展开更多
The Menyuan area is an important transportation hub in the Hexi Corridor.The Menyuan M_(S)6.9 earthquake that occurred on January 8,2022 had a major impact on the local infrastructure and transportation of this region...The Menyuan area is an important transportation hub in the Hexi Corridor.The Menyuan M_(S)6.9 earthquake that occurred on January 8,2022 had a major impact on the local infrastructure and transportation of this region.Due to the high possibility of similar strong earthquakes occurring in this area in the future,preliminary assessment of the seismic intensity characteristics of destructive earthquakes in this region is essential for effective disaster control.This paper uses the empirical Green′s function(EGF)method as a numerical simulation tool to predict the ground motion intensity of Datong Autonomous County under the action of the scenario earthquake(M_(S)7.5).Seismic records of aftershocks of the 2016 Menyuan M_(S)6.4 earthquake were used as Green’s functions for this simulation.The uncertainties associated with various source parameters were considered,and 36possible earthquake scenarios were simulated to obtain 72 sets of horizontal ground motions in Datong County.The obtained peak ground acceleration(PGA)vs.time histories of the horizontal ground motion were screened using the attenuation relationships provided by the fifth-edition of China’s Seismic Ground Motion Parameter Zoning Map and the NGA-West2dataset.Ultimately,32 possible acceleration-time histories were selected for further analysis.The screened PGA values ranged from 78.8 to 153 cm/s^(2).The uncertainty associated with the initial rupture point was found to greatly affect the results of the earthquake simulation.The average acceleration spectrum of the selected acceleration-time history exceeded the expected spectrum of a intermediate earthquake,which means that buildings in Datong County might sustain some damage should the scenario earthquake occur.This research can provide reliable ground motion input for urban earthquake damage simulation and seismic design in Datong County.Growing the dataset of small earthquakes recorded in this region will facilitate the large-scale simulation of ground motions under different earthquake scenarios.展开更多
In this paper, the solution of Chebyshev equation with its argument being greater than 1 is obtained. The initial value of the derivative of the solution is the expression of magnetization, which is valid for any spin...In this paper, the solution of Chebyshev equation with its argument being greater than 1 is obtained. The initial value of the derivative of the solution is the expression of magnetization, which is valid for any spin quantum number S. The Chebyshev equation is transformed from an ordinary differential equation obtained when we dealt with Heisenberg model, in order to calculate all three components of magnetization, by many-body Green＇s function under random phase approximation. The Chebyshev functions with argument being greater than 1 are discussed. This paper shows that the Chebyshev polynomials with their argument being greater than 1 have their physical application.展开更多
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 co...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.展开更多
Taking Tianjin as an example, this paper proposed a methodology and process for evaluating near-fault strong ground motions from future earthquakes to mitigate earthquake damage for the metropolitan area and important...Taking Tianjin as an example, this paper proposed a methodology and process for evaluating near-fault strong ground motions from future earthquakes to mitigate earthquake damage for the metropolitan area and important engineering structures. The result of strong ground motion was predicted for Tianjin main faults by the hybrid method which mainly con- sists of 3D finite difference method and stochastic Green＇s function. Simulation is performed for 3D structures of Tianjin re- gion and characterized asperity models. The characterized asperity model describing source heterogeneity is introduced fol- lowing the fault information from the project of Tianjin Active Faults and Seismic Hazard Assessment. We simulated the worst case that two earthquakes separately occur. The results indicate that the fault position, rupture process and the sedi- mentary deposits of the basin significantly affect amplification of the simulated ground motion. Our results also demonstrate the possibility of practical simulating wave propagation including basin induced surface waves in broad frequency-band, for seismic hazard analysis near the fault from future earthquakes in urbanized areas.展开更多
The Cauchy problem of the generalized Kuramoto-Sivashinsky equation in multidimensions(n ≥ 3) is considered. Based on Green's function method, some ingenious energy estimates are given. Then the global existence ...The Cauchy problem of the generalized Kuramoto-Sivashinsky equation in multidimensions(n ≥ 3) is considered. Based on Green's function method, some ingenious energy estimates are given. Then the global existence and pointwise convergence rates of the classical solutions are established. Furthermore, the L^p convergence rate of the solution is obtained.展开更多
In this study, the thermoelectric properties of pyrene molecule doped with boron and nitrogen atom at different sites of molecule are investigated using density functional theory and none-equilibrium Greens function f...In this study, the thermoelectric properties of pyrene molecule doped with boron and nitrogen atom at different sites of molecule are investigated using density functional theory and none-equilibrium Greens function formalism in the linear response regime. Our calculations show that when the impurities are added to the edge of the molecule, the anti-resonant peaks will appear in the transmission diagram in the vicinity of the Fermi energy level. So it increases the thermoelectric figure of merit of the system in comparison with the one that the impurity is located in the center of molecule. Additionally, the seebeck coefficient signs are not the same among the B, N, and N ＆ B doped devices, indicating that the types of the carriers can be changed with different types of doping.展开更多
文摘The heavy quarks present in the quark-gluon plasma(QGP)can act as a probe of relativistic heavy ion collisions as they retain the memory of their interaction history.In a previous study,a stochastic Schrödinger equation(SSE)has been applied to describe the evolution of heavy quarks,where an external field with random phases is used to simulate the thermal medium.In this work,we study the connection between the SSE and the Boltzmann transport equation(BE)approach in the Keldysh Green’s function formalism.By comparing the Green’s function of the heavy quark from the SSE and the Keldysh Green’s functions leading to the Boltzmann equation,we demonstrate that the SSE is consistent with the Boltzmann equation in the weak coupling limit.We subsequently confirm their consistency through numerical calculations.
基金supported by the National Natural Science Foundation of China (Grant Nos. 12272269, 11972257,11832014 and 11472193)the Shanghai Pilot Program for Basic Researchthe Shanghai Gaofeng Project for University Academic Program Development。
文摘Understanding the effects of point liquid loading on transversely isotropic poroelastic media is crucial for advancing geomechanics and biomechanics, where precise modeling of fluid-structure interactions is essential. This paper presents a comprehensive analysis of infinite transversely isotropic poroelasticity under a fluid source, based on Biot's theory, aiming to uncover new and previously unexplored insights in the literature. We begin our study by deriving a general solution for fluid-saturated, transversely isotropic poroelastic materials in terms of harmonic functions that satisfy sixth-order homogeneous partial differential equations, using potential theory and Almansi's theorem. Based on these general solutions and potential functions, we construct a Green's function for a point fluid source, introducing three new harmonic functions with undetermined constants. These constants are determined by enforcing continuity and equilibrium conditions. Substituting these into the general solution yields fundamental solutions for poroelasticity that provide crucial support for a wide range of project problems. Numerical results and comparisons with existing literature are provided to illustrate physical mechanisms through contour plots. Our observations reveal that all components tend to zero in the far field and become singular at the concentrated source. Additionally, the contours exhibit rapid changes near the point fluid source but display gradual variations at a distance from it. These findings highlight the intricate behavior of the system under point liquid loading, offering valuable insights for further research and practical applications.
基金supported by the National Natural Science Foundation of China (Grant No. 50879090)
文摘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.
基金Project supported by the National Natural Science Foundation of China （No. 50879090）, and the Key Research Program of Hydrody- namics of China （No. 9140A 14030712JB 11044）
文摘The singularities and oscillatory performance of translating-pulsating source Green＇s function in Bessho form were analyzed. Relative numerical integration methods such as Gaussian quadrature rule, variable substitution method （VSM）, and steepest descent integration method （SDIM） were used to evaluate this type of Green＇s function. For SDIM, the complex domain was restricted only on the 0-plane. Meanwhile, the integral along the real axis was computed by use of the VSM to avoid the complication of a numerical search of the steepest descent line. Furthermore, the steepest descent line was represented by the B-spline function. Based on this representation, a new self-compatible integration method corresponding to parametric t was established. The numerical method was validated through comparison with other existing results, and was shown to be efficient and reliable in the calculation of the velocity potentials for the 3D seakeeping and hydrodynamic performance of floating struc- tures moving in waves.
基金the National Natural Science Foundation of China(No.U2032141)the Open Project of Guangxi Key Laboratory of Nuclear Physics and Nuclear Technology(No.NLK2022-02)+4 种基金the Central Government Guidance Funds for Local Scientific and Technological Development,China(Guike ZY22096024)the Natural Science Foundation of Henan Province(No.202300410479)the Guizhou Provincial Science and Technology Projects(No.ZK[2022]203)the Foundation of Fundamental Research for Young Teachers of Zhengzhou University(No.JC202041041)the Physics Research and Development Program of Zhengzhou University(No.32410217).
文摘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.
文摘In this paper, the Adomian decomposition method with Green’s function (Standard Adomian and Modified Technique) is applied to solve linear and nonlinear tenth-order boundary value problems with boundary conditions defined at any order derivatives. The numerical results obtained with a small amount of computation are compared with the exact solutions to show the efficiency of the method. The results show that the decomposition method is of high accuracy, more convenient and efficient for solving high-order boundary value problems.
基金National Natural Science Foundation of China under grant No.51578373 and 51578372the Natural Science Foundation of Tianjin Municipality under Grant No.16JCYBJC21600
文摘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.
基金Aeronautical Science Foundation of China （2008ZC52024）
文摘The accurate analyses for a plate fin heat sink with the ability to control the temperature of the avionics devices within a pre-set controllable temperature range are required both in the process of circuit design and for the real-time temperature monitoring purposes. In order to provide an insight into the behavior of the temperature of a plate fin heat sink subjected non-uniform heat density on the surfaces, it is necessary to obtain accurate analytical solutions yielding explicit formulas relating the dissipated power to the temperature rise at any point of avionics devices. This paper presents a method for thermal simulation of a plate fin heat sink using an analytical solution of the three-dimensional heat equation resulting from an appropriate plate fin heat sink transient thermal model. The entire solution methodology is illustrated in detail on the particular examples of the plate fin heat sink subjected non-uniform heat density on the surfaces. The transient temperature profiles are obtained for different positions at the surface of the plate fin heat sink. The analytical results are compared with measurements made on the surface of the cold plate and it is found that they are in good agreement with an error of less than 3 K.
文摘The principle objective of the paper is to study the acoustic radiation problem of the 3D space domain with control boundary. By using the conformal transformation theory, the Green's function for acoustic point source in the control domain space is obtained. With it, the expression of acoustic radiation function of the control domain is formed. Discussion about the acoustic radiation of pulsating sphere in right-angle space is drawn in the end, including the acoustic radiation directivity effect by the boundary characteristics, acoustic radiation frequency and acoustic source location. Numerical results show that: for the lower frequency radiation, the infection of free surface is significant; for the high frequency radiation, the infection of location is significant on the contrary. The research provides a new method for boundary characteristic problem of the structural-acoustic acoustic.
基金supported by National Engineering School of Tunis (No.13039.1)
文摘To reduce computational costs, an improved form of the frequency domain boundary element method（BEM） is proposed for two-dimensional radiation and propagation acoustic problems in a subsonic uniform flow with arbitrary orientation. The boundary integral equation（BIE） representation solves the two-dimensional convected Helmholtz equation（CHE） and its fundamental solution, which must satisfy a new Sommerfeld radiation condition（SRC） in the physical space. In order to facilitate conventional formulations, the variables of the advanced form are expressed only in terms of the acoustic pressure as well as its normal and tangential derivatives, and their multiplication operators are based on the convected Green＇s kernel and its modified derivative. The proposed approach significantly reduces the CPU times of classical computational codes for modeling acoustic domains with arbitrary mean flow. It is validated by a comparison with the analytical solutions for the sound radiation problems of monopole,dipole and quadrupole sources in the presence of a subsonic uniform flow with arbitrary orientation.
基金supported by the National Natural Science Foundation of China (Nos. 10732100, 10572155)the Science and Technology Planning Project of Guangdong Province of China (No. 2006A11001002)the Ph. D. Programs Foundation of Ministry of Education of China (No. 2006300004111179)
文摘Using a polarization method, the scattering problem for a two-dimensional inclusion embedded in infinite piezoelectric/piezomagnetic matrices is investigated. To achieve the purpose, the polarization method for a two-dimensional piezoelectric/piezomagnetic ＂comparison body＂ is formulated. For simple harmonic motion, kernel of the polarization method reduces to a 2-D time-harmonic Green＇s function, which is obtained using the Radon transform. The expression is further simplified under conditions of low frequency of the incident wave and small diameter of the inclusion. Some analytical expressions are obtained. The analytical solutions for generalized piezoelectric/piezomagnetic anisotropic composites are given followed by simplified results for piezoelectric composites. Based on the latter results, two numerical results are provided for an elliptical cylindrical inclusion in a PZT-5H-matrix, showing the effect of different factors including size, shape, material properties, and piezoelectricity on the scattering cross-section.
基金National Natural Science Foundation of China under Grant Nos.51578373 and 51578372
文摘Abstract： The scattering of plane SH-waves by topographic features in a layered transversely isotropic （TI） half-space is investigated by using an indirect boundary element method （IBEM）. Firstly, the anti-plane dynamic stiffness matrix of the layered TI half-space is established and the free fields are solved by using the direct stiffness method. Then, Green＇s functions are derived for uniformly distributed loads acting on an inclined line in a layered TI half-space and the scattered fields are constructed with the deduced Green＇s functions. Finally, the free fields are added to the scattered ones to obtain the global dynamic responses. The method is verified by comparing results with the published isotropic ones. Both the steady-state and transient dynamic responses are evaluated and discussed. Numerical results in the frequency domain show that surface motions for the TI media can be significantly different from those for the isotropic case, which are strongly dependent on the anisotropy property, incident angle and incident frequency. Results in the time domain show that the material anisotropy has important effects on the maximum duration and maximum amplitudes of the time histories.
基金jointly supported by the National Key Research and Development Program(No.2019YFC1511004-02)Youth Fund of the National Natural Science Foundation(No.42104053)the Research Project Fund of the Institute of Geophysics,China Earthquake Administration(No.DQJB22B21)。
文摘The Menyuan area is an important transportation hub in the Hexi Corridor.The Menyuan M_(S)6.9 earthquake that occurred on January 8,2022 had a major impact on the local infrastructure and transportation of this region.Due to the high possibility of similar strong earthquakes occurring in this area in the future,preliminary assessment of the seismic intensity characteristics of destructive earthquakes in this region is essential for effective disaster control.This paper uses the empirical Green′s function(EGF)method as a numerical simulation tool to predict the ground motion intensity of Datong Autonomous County under the action of the scenario earthquake(M_(S)7.5).Seismic records of aftershocks of the 2016 Menyuan M_(S)6.4 earthquake were used as Green’s functions for this simulation.The uncertainties associated with various source parameters were considered,and 36possible earthquake scenarios were simulated to obtain 72 sets of horizontal ground motions in Datong County.The obtained peak ground acceleration(PGA)vs.time histories of the horizontal ground motion were screened using the attenuation relationships provided by the fifth-edition of China’s Seismic Ground Motion Parameter Zoning Map and the NGA-West2dataset.Ultimately,32 possible acceleration-time histories were selected for further analysis.The screened PGA values ranged from 78.8 to 153 cm/s^(2).The uncertainty associated with the initial rupture point was found to greatly affect the results of the earthquake simulation.The average acceleration spectrum of the selected acceleration-time history exceeded the expected spectrum of a intermediate earthquake,which means that buildings in Datong County might sustain some damage should the scenario earthquake occur.This research can provide reliable ground motion input for urban earthquake damage simulation and seismic design in Datong County.Growing the dataset of small earthquakes recorded in this region will facilitate the large-scale simulation of ground motions under different earthquake scenarios.
基金The project supported by the State Key Project of Fundamental Research of China under Grant No. G2000067101
文摘In this paper, the solution of Chebyshev equation with its argument being greater than 1 is obtained. The initial value of the derivative of the solution is the expression of magnetization, which is valid for any spin quantum number S. The Chebyshev equation is transformed from an ordinary differential equation obtained when we dealt with Heisenberg model, in order to calculate all three components of magnetization, by many-body Green＇s function under random phase approximation. The Chebyshev functions with argument being greater than 1 are discussed. This paper shows that the Chebyshev polynomials with their argument being greater than 1 have their physical application.
基金supported by National Natural Science Foundation of China under Grant Nos.10574121,10874160,and 10025420the‘111’Project of the Ministry of Education and the Chinese Academy of Sciences
文摘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.
基金supported by the National Natural Science Foundation of China (40674016)Earthquake Research Fund of Ministry of Science and Technology of China ((2007)203/8-53)
文摘Taking Tianjin as an example, this paper proposed a methodology and process for evaluating near-fault strong ground motions from future earthquakes to mitigate earthquake damage for the metropolitan area and important engineering structures. The result of strong ground motion was predicted for Tianjin main faults by the hybrid method which mainly con- sists of 3D finite difference method and stochastic Green＇s function. Simulation is performed for 3D structures of Tianjin re- gion and characterized asperity models. The characterized asperity model describing source heterogeneity is introduced fol- lowing the fault information from the project of Tianjin Active Faults and Seismic Hazard Assessment. We simulated the worst case that two earthquakes separately occur. The results indicate that the fault position, rupture process and the sedi- mentary deposits of the basin significantly affect amplification of the simulated ground motion. Our results also demonstrate the possibility of practical simulating wave propagation including basin induced surface waves in broad frequency-band, for seismic hazard analysis near the fault from future earthquakes in urbanized areas.
基金supported by the National Natural Science Foundation of China(11271141)Chongqing Science&Technology Commission(cstc2018jcyjAX0787)
文摘The Cauchy problem of the generalized Kuramoto-Sivashinsky equation in multidimensions(n ≥ 3) is considered. Based on Green's function method, some ingenious energy estimates are given. Then the global existence and pointwise convergence rates of the classical solutions are established. Furthermore, the L^p convergence rate of the solution is obtained.
文摘In this study, the thermoelectric properties of pyrene molecule doped with boron and nitrogen atom at different sites of molecule are investigated using density functional theory and none-equilibrium Greens function formalism in the linear response regime. Our calculations show that when the impurities are added to the edge of the molecule, the anti-resonant peaks will appear in the transmission diagram in the vicinity of the Fermi energy level. So it increases the thermoelectric figure of merit of the system in comparison with the one that the impurity is located in the center of molecule. Additionally, the seebeck coefficient signs are not the same among the B, N, and N ＆ B doped devices, indicating that the types of the carriers can be changed with different types of doping.
