Defect engineering is a commonly methodology used to enhance the thermoelectric performance of thermoelectric PbTe by improving its electronic transport properties.At the nanoscale,defects can induce quantum tunneling...Defect engineering is a commonly methodology used to enhance the thermoelectric performance of thermoelectric PbTe by improving its electronic transport properties.At the nanoscale,defects can induce quantum tunneling effects that significantly impact the electrical properties of materials.To understand the specific mechanisms underlying the quantum transport properties of PbTe,we employ the non-equilibrium Green's function(NEGF)method to investigate the effects of intrinsic defects(point defects and grain boundaries)on the electronic transport properties of PbTe-based nanodevices from a quantum mechanical perspective.Our results show that the Pb vacancy(VPb)has the highest conduction.The conduction depends on the defect type,chemical potential and bias voltage.The presence of intrinsic point defects introduces impurity levels,facilitating the electron tunneling and leading to an increase in the transmission coefficient,thereby enhancing the electronic transport properties.For PbTe containing grain boundaries,these boundaries suppress the electronic transport properties.The Te occupied twin boundary(Te-TB)exerts a stronger inhibitory effect than the Pb occupied twin boundary(Pb-TB).Nevertheless,the combined effect between twin boundaries and point defects can enhance the electrical properties.Therefore,in order to obtain highly conductive of PbTe materials,a Te-rich synthesis environment should be used to promote the effective formation of Pb vacancy.Our work offers a comprehensive understanding of the impact of defects on electron scattering in thermoelectric materials.展开更多
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 derivation of Green function in a two-layer fluid model has been treated in different ways. In a two-layer fluid with the upper layer having a free surface, there exist two modes of waves propagating due to the fr...The derivation of Green function in a two-layer fluid model has been treated in different ways. In a two-layer fluid with the upper layer having a free surface, there exist two modes of waves propagating due to the free surface and the interface. This paper is concerned with the derivation of Green functions in the three dimensional case of a stationary source oscillating. The source point is located either in the upper or lower part of a two-layer fluid of finite depth. The derivation is carried out by the method of singularities. This method has an advantage in that it involves representing the potential as a sum of singularities or multipoles placed within any structures being present. Furthermore, experience shows that the systems of equations resulted from using a singularity method possess excellent convergence characteristics and only a few equations are needed to obtain accurate numerical results. Validation is done by showing that the derived two-layer Green function can be reduced to that of a single layer of finite depth or that the upper Green function coincides with that of the lower, for each case. The effect of the density on the internal waves is demonstrated. Also, it is shown how the surface and internal wave amplitudes are compared for both the wave modes. The fluid in this case is considered to be inviscid and incompressible and the flow is irrotational.展开更多
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.展开更多
A novel numerical model based on the image Green function and first-order Taylor expansion boundary element method(TEBEM), which can improve the accuracy of the hydrodynamic simulation for the non-smooth body, was dev...A novel numerical model based on the image Green function and first-order Taylor expansion boundary element method(TEBEM), which can improve the accuracy of the hydrodynamic simulation for the non-smooth body, was developed to calculate the side wall effects on first-order motion responses and second-order drift loads upon offshore structures in the wave tank. This model was confirmed by comparing it to the results from experiments on hydrodynamic coefficients, namely the first-order motion response and second-order drift load upon a hemisphere, prolate spheroid, and box-shaped barge in the wave tank. Then,the hydrodynamics of the KVLCC2 model were also calculated in two wave tanks with different widths. It was concluded that this model can predict the hydrodynamics for offshore structures effectively, and the side wall has a significant impact on the firstorder quantities and second-order drift loads, which satisfied the resonant rule.展开更多
To solve the numerical divergence problem of the direct time domain Green function method for the motion simulation of floating bodies with large flare, a time domain hybrid Rankine-Green boundary element method is pr...To solve the numerical divergence problem of the direct time domain Green function method for the motion simulation of floating bodies with large flare, a time domain hybrid Rankine-Green boundary element method is proposed. In this numerical method, the fluid domain is decomposed by an imaginary control surface, at which the continuous condition should be satisfied. Then the Rankine Green function is adopted in the inner domain. The transient free surface Green function is applied in the outer domain, which is used to find the relationship between the velocity potential and its normal derivative for the inner domain. Besides, the velocity potential at the mean free surface between body surface and control surface is directly solved by the integration scheme. The wave exciting force is computed through the convolution integration with wave elevation, by introducing the impulse response function. Additionally, the nonlinear Froude-Krylov force and hydrostatic force, which is computed under the instantaneous incident wave free surface, are taken into account by the direct pressure integration scheme. The corresponding numerical computer code is developed and first used to compute the hydrodynamic coefficients of the hemisphere, as well as the time history of a ship with large flare; good agreement is obtained with the analytical solutions as well as the available numerical results. Then the hydrodynamic properties of a FPSO are studied. The hydrodynamic coefficients agree well with the results computed by the frequency method; the influence of the time interval and the truncated time is investigated in detail.展开更多
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.展开更多
A new type of dual boundary integral equations(DBIE)is presented first,through which,a smaller system of equations needs to be solved in fracture analysis.Then a non-conforming crack tip element in two-dimensional pro...A new type of dual boundary integral equations(DBIE)is presented first,through which,a smaller system of equations needs to be solved in fracture analysis.Then a non-conforming crack tip element in two-dimensional problems is proposed.The exact formula for the hypersingular integral over the non-con- forming crack tip element is given next.By virtue of Green's-function-library strategy,a series of stress in- tensity factors(SIF)of different crack orientations,locations and/or sizes in a complicated structure can be obtained easily and efficiently.Finally,several examples of fracture analysis in two dimensions are given to demonstrate the accuracy and efficiency of the method proposed.展开更多
基金financial support from the National Natural Science Foundation of China(No.12474016)the program of“Distinguished Expert of Taishan Scholar”(No.tstp20221124)+4 种基金the National Natural Science Foundation of China(Nos.52172212,12474017)the Shandong Provincial Science Foundation(ZR2021YQ03)the program for“Young Scientists of Taishan Scholars”(No.tsqn202306184)financial support from the National Natural Science Foundation of China(No.12464034)the Natural Science Foundation of Ningxia,China(No.2024AAC05070)。
文摘Defect engineering is a commonly methodology used to enhance the thermoelectric performance of thermoelectric PbTe by improving its electronic transport properties.At the nanoscale,defects can induce quantum tunneling effects that significantly impact the electrical properties of materials.To understand the specific mechanisms underlying the quantum transport properties of PbTe,we employ the non-equilibrium Green's function(NEGF)method to investigate the effects of intrinsic defects(point defects and grain boundaries)on the electronic transport properties of PbTe-based nanodevices from a quantum mechanical perspective.Our results show that the Pb vacancy(VPb)has the highest conduction.The conduction depends on the defect type,chemical potential and bias voltage.The presence of intrinsic point defects introduces impurity levels,facilitating the electron tunneling and leading to an increase in the transmission coefficient,thereby enhancing the electronic transport properties.For PbTe containing grain boundaries,these boundaries suppress the electronic transport properties.The Te occupied twin boundary(Te-TB)exerts a stronger inhibitory effect than the Pb occupied twin boundary(Pb-TB).Nevertheless,the combined effect between twin boundaries and point defects can enhance the electrical properties.Therefore,in order to obtain highly conductive of PbTe materials,a Te-rich synthesis environment should be used to promote the effective formation of Pb vacancy.Our work offers a comprehensive understanding of the impact of defects on electron scattering in thermoelectric materials.
基金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.
基金supported by the National Natural Science Foundation of China (Grant No. 50779008)
文摘The derivation of Green function in a two-layer fluid model has been treated in different ways. In a two-layer fluid with the upper layer having a free surface, there exist two modes of waves propagating due to the free surface and the interface. This paper is concerned with the derivation of Green functions in the three dimensional case of a stationary source oscillating. The source point is located either in the upper or lower part of a two-layer fluid of finite depth. The derivation is carried out by the method of singularities. This method has an advantage in that it involves representing the potential as a sum of singularities or multipoles placed within any structures being present. Furthermore, experience shows that the systems of equations resulted from using a singularity method possess excellent convergence characteristics and only a few equations are needed to obtain accurate numerical results. Validation is done by showing that the derived two-layer Green function can be reduced to that of a single layer of finite depth or that the upper Green function coincides with that of the lower, for each case. The effect of the density on the internal waves is demonstrated. Also, it is shown how the surface and internal wave amplitudes are compared for both the wave modes. The fluid in this case is considered to be inviscid and incompressible and the flow is irrotational.
基金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 National Natural Science Foundation of China (Grant No.51709064)the Numerical Tank Project sponsored by the Ministry of Industry and Information Technology (MIIT)of P.R.China.
文摘A novel numerical model based on the image Green function and first-order Taylor expansion boundary element method(TEBEM), which can improve the accuracy of the hydrodynamic simulation for the non-smooth body, was developed to calculate the side wall effects on first-order motion responses and second-order drift loads upon offshore structures in the wave tank. This model was confirmed by comparing it to the results from experiments on hydrodynamic coefficients, namely the first-order motion response and second-order drift load upon a hemisphere, prolate spheroid, and box-shaped barge in the wave tank. Then,the hydrodynamics of the KVLCC2 model were also calculated in two wave tanks with different widths. It was concluded that this model can predict the hydrodynamics for offshore structures effectively, and the side wall has a significant impact on the firstorder quantities and second-order drift loads, which satisfied the resonant rule.
基金the financial support provided by the National Basic Research Program of China(No.2011CB3703)the National Natural Science Foundation of China(No.51079034)
文摘To solve the numerical divergence problem of the direct time domain Green function method for the motion simulation of floating bodies with large flare, a time domain hybrid Rankine-Green boundary element method is proposed. In this numerical method, the fluid domain is decomposed by an imaginary control surface, at which the continuous condition should be satisfied. Then the Rankine Green function is adopted in the inner domain. The transient free surface Green function is applied in the outer domain, which is used to find the relationship between the velocity potential and its normal derivative for the inner domain. Besides, the velocity potential at the mean free surface between body surface and control surface is directly solved by the integration scheme. The wave exciting force is computed through the convolution integration with wave elevation, by introducing the impulse response function. Additionally, the nonlinear Froude-Krylov force and hydrostatic force, which is computed under the instantaneous incident wave free surface, are taken into account by the direct pressure integration scheme. The corresponding numerical computer code is developed and first used to compute the hydrodynamic coefficients of the hemisphere, as well as the time history of a ship with large flare; good agreement is obtained with the analytical solutions as well as the available numerical results. Then the hydrodynamic properties of a FPSO are studied. The hydrodynamic coefficients agree well with the results computed by the frequency method; the influence of the time interval and the truncated time is investigated in detail.
文摘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.
基金the Aeronautical Science Foundation of China (No.99C53026).
文摘A new type of dual boundary integral equations(DBIE)is presented first,through which,a smaller system of equations needs to be solved in fracture analysis.Then a non-conforming crack tip element in two-dimensional problems is proposed.The exact formula for the hypersingular integral over the non-con- forming crack tip element is given next.By virtue of Green's-function-library strategy,a series of stress in- tensity factors(SIF)of different crack orientations,locations and/or sizes in a complicated structure can be obtained easily and efficiently.Finally,several examples of fracture analysis in two dimensions are given to demonstrate the accuracy and efficiency of the method proposed.