Integral diffraction coefficients of the crystal are the essential data of a crystal spectrometer which is extensively used to measure quantitative x-ray spectra of high temperature plasmas in kilo-electron-volt regio...Integral diffraction coefficients of the crystal are the essential data of a crystal spectrometer which is extensively used to measure quantitative x-ray spectra of high temperature plasmas in kilo-electron-volt region. An experimental method has been developed to measure the integral diffraction coefficients of crystals on beamline 4B7 of Beijing Synchrotron Radiation Facility. The integral diffraction coefficients of several crystals including polyethylene terephthalate （PET）, thallium acid phthalate （T1AP） and rubidium acid phthalate （RAP） crystals have been measured in the x-ray energy range 2100-5600 eV and compared with the calculations of the ＇Darwin Prins＇ and the ＇Mosaic＇ models. It is shown that the integral diffraction coefficients of these crystals are between the calculations of the ＇Darwin Prins＇ and the ＇Mosaic＇ models, but more close to the ＇Darwin Prins＇ model calculations.展开更多
In this paper, Kirchhoff formula has been transformed from surface integral form into a line integral form. The new form of the formula can be applied to separate geometrical optical field from diffraction field, and ...In this paper, Kirchhoff formula has been transformed from surface integral form into a line integral form. The new form of the formula can be applied to separate geometrical optical field from diffraction field, and reduce the time of numerical computation greatly. Based on the new form, an analytical formula of diffraction field in the far zone has been presented for the polygonal aperture illuminated by a uniform plane wave.展开更多
A Fast Multipole Method (FMM) is developed as a numerical approach to the reduction of the computational cost and requirement memory capacity for a large in solving large-scale problems. In this paper it is applied to...A Fast Multipole Method (FMM) is developed as a numerical approach to the reduction of the computational cost and requirement memory capacity for a large in solving large-scale problems. In this paper it is applied to the boundary integral equation method (BIEM) for current diffraction from arbitrary 3D bodies. The boundary integral equation is discretized by higher order elements, the FMM is applied to avoid the matrix/vector product, and the resulting algebraic equation is solved by the Generalized Conjugate Residual method (GCR). Numerical examination shows that the FMM is more efficient than the direct evaluation method in computational cost and storage of computers.展开更多
Diffraction of plane P waves around an alluvial valley of arbitrary shape in poroelastic half-space is investigated by using an indirect boundary integral equation method. Based on the Green＇s fimctions of line sourc...Diffraction of plane P waves around an alluvial valley of arbitrary shape in poroelastic half-space is investigated by using an indirect boundary integral equation method. Based on the Green＇s fimctions of line source in poroelastic half-space, the scattered waves are constructed using the fictitious wave sources close to the interface of the valley and the density of ficti- tious wave sources are determined by boundary conditions. The precision of the method is verified by the satisfaction extent of boundary conditions, and the comparison between the degenerated solutions and available results in single-phase case. Finally, the nature of diffraction of plane P waves around an alluvial valley in poroelastic half-space is investigated in detail through nu- merical examples.展开更多
This review analyzes following numerical methods of a solution of problems of a sound diffraction on ideal and elastic scatterers of a non-analytical form: a method of integral equations, a method of Green’s function...This review analyzes following numerical methods of a solution of problems of a sound diffraction on ideal and elastic scatterers of a non-analytical form: a method of integral equations, a method of Green’s functions, a method of finite elements, a boundary elements method, a method of Kupradze, a T-matrix method and a method of a geometrical theory of a diffraction.展开更多
Silicon(Si)diffraction microlens arrays are usually used to integrating with infrared focal plane arrays(IRFPAs)to improve their performance.The errors of lithography are unavoidable in the process of the Si diffrac-t...Silicon(Si)diffraction microlens arrays are usually used to integrating with infrared focal plane arrays(IRFPAs)to improve their performance.The errors of lithography are unavoidable in the process of the Si diffrac-tion microlens arrays preparation in the conventional engraving method.It has a serious impact on its performance and subsequent applications.In response to the problem of errors of Si diffraction microlens arrays in the conven-tional method,a novel self-alignment method for high precision Si diffraction microlens arrays preparation is pro-posed.The accuracy of the Si diffractive microlens arrays preparation is determined by the accuracy of the first li-thography mask in the novel self-alignment method.In the subsequent etching,the etched area will be protected by the mask layer and the sacrifice layer or the protective layer.The unprotection area is carved to effectively block the non-etching areas,accurately etch the etching area required,and solve the problem of errors.The high precision Si diffraction microlens arrays are obtained by the novel self-alignment method and the diffraction effi-ciency could reach 92.6%.After integrating with IRFPAs,the average blackbody responsity increased by 8.3%,and the average blackbody detectivity increased by 10.3%.It indicates that the Si diffraction microlens arrays can improve the filling factor and reduce crosstalk of IRFPAs through convergence,thereby improving the perfor-mance of the IRFPAs.The results are of great reference significance for improving their performance through opti-mizing the preparation level of micro nano devices.展开更多
The development of an in-house computer program for determining the motions and loads of advancing ships through sea waves in the frequency domain,is described in this paper.The code is based on the potential flow for...The development of an in-house computer program for determining the motions and loads of advancing ships through sea waves in the frequency domain,is described in this paper.The code is based on the potential flow formulation and originates from a double-body code enhanced with the regular part of the velocity potential computed using the pulsing source Green function.The code is fully developed in C++language with extensive use of the object-oriented paradigm.The code is capable of estimating the excitation and inertial radiation loads or arbitrary incoming wave frequencies and incidence angles.The hydrodynamic responses such as hydrodynamic coefficients,ship motions,the vertical shear force and the vertical bending moment are estimated.A benchmark container ship and an LNG carrier are selected for testing and validating the computer code.The obtained results are compared with the available experimental data which demonstrate the acceptable compliance for the zero speed whereas there are some discrepancies over the range of frequencies for the advancing ship in different heading angles.展开更多
Graphene oxide （GO） ultrathin flat lenses have provided a new and viable solution to achieve high resolution, high efficiency, ultra-light weight, integratable and flexible optical systems. Current GO lenses are des...Graphene oxide （GO） ultrathin flat lenses have provided a new and viable solution to achieve high resolution, high efficiency, ultra-light weight, integratable and flexible optical systems. Current GO lenses are designed based on the Fresnel diffraction model, which uses a paraxial approximation for low numerical aperture （NA） focusing process. Herein we develop a lens design method based on the Rayleigh-Sommerfeld （RS） diffraction theory that is able to unambiguously determine the radii of each ring without the optimization process for the first time. More importantly, the RS design method is able to accurately design GO lenses with arbitrary NA and focal length. Our design is experimentally confirmed by fabricating high NA GO lenses with both short and long focal lengths. Compared with the conventional Fresnel design methods, the differences in ring positions and the resulted focal length are up to 13.9% and 9.1%, respectively. Our method can be further applied to design high performance flat lenses of arbitrary materials given the NA and focal length requirements, including metasurfaces or other two-dimensional materials.展开更多
This paper is concerned with the problem of scattering of time-harmonic electromag- netic waves from penetrable diffraction gratings in the 2D polarization case. We propose a new, weakly singular, integral equation fo...This paper is concerned with the problem of scattering of time-harmonic electromag- netic waves from penetrable diffraction gratings in the 2D polarization case. We propose a new, weakly singular, integral equation formulation for the scattering problem which is proved to be uniquely solvable. A main feature of the new integral equation formula- tion is that it avoids the computation of the normal derivative of double-layer potentials which is difficult and time consuming. A fast numerical algorithm is also developed for the scattering problem, based on the NystrSm method for the new integral equation. Nu- merical examples are also shown to illustrate the applicability of the new integral equation formulation.展开更多
The integral equation method for the simulation of the diffraction by optical gratings is an efficient numerical tool if profile gratings determined by simple crosssection curves are considered.This method in its rece...The integral equation method for the simulation of the diffraction by optical gratings is an efficient numerical tool if profile gratings determined by simple crosssection curves are considered.This method in its recent version is capable to tackle profile curves with corners,gratings with thin coated layers,and diffraction scenarios with unfavorably large ratio period over wavelength.We discuss special implementational issues including the efficient evaluation of the quasi-periodic Green kernels,the quadrature algorithm,and the iterative solution of the arising systems of linear equations.Finally,as an example we present the simulation of echelle gratings which demonstrates the efficency of our approach.展开更多
A novel method to generate a collimated hollow-laser-beam (HLB) by only a single axicon is proposed. With some reasonable assumptions, the radial light intensity distribution is calculated in detail by diffraction int...A novel method to generate a collimated hollow-laser-beam (HLB) by only a single axicon is proposed. With some reasonable assumptions, the radial light intensity distribution is calculated in detail by diffraction integral theory.The result of numerical simulation shows that this method is valid.Compared with other methods of generating HLB,this scheme is extraordinarily simple in principle and can be utilized experimentally to construct a light trap in atomic fountain for convenience.展开更多
In this letter, we propose a method for the numerical calculations of the femtosecond laser pulse passed through a subwavelength aperture. The time-dependent laser pulse is decomposed into a series of monochromatic si...In this letter, we propose a method for the numerical calculations of the femtosecond laser pulse passed through a subwavelength aperture. The time-dependent laser pulse is decomposed into a series of monochromatic simple harmonic waves. For the light field of the harmonic wave with a single frequency, the numerical calculation is made based on the solution of the Green's integral equation set of the electromagnetic waves. Such numerical solution is iterated for all the waves with different frequencies, and all the numerical solutions are transformed into the light fields in the time domain by inverse Fourier transform. The light intensity distributions transmitted the subwavelength aperture are calculated and the results show the propagation of the light field is along the direction of the medium interface.展开更多
基金Project supported by the Foundation of National Key Laboratory of High Temperature and Density Plasma Physics(Grant No.9140C6804020704)the National Natural Science Foundation of China(Grant No.10874156)one of the authors,GanXin-Shi,was also partially supported by the School of Physical Science and Technology in Sichuan University,China
文摘Integral diffraction coefficients of the crystal are the essential data of a crystal spectrometer which is extensively used to measure quantitative x-ray spectra of high temperature plasmas in kilo-electron-volt region. An experimental method has been developed to measure the integral diffraction coefficients of crystals on beamline 4B7 of Beijing Synchrotron Radiation Facility. The integral diffraction coefficients of several crystals including polyethylene terephthalate （PET）, thallium acid phthalate （T1AP） and rubidium acid phthalate （RAP） crystals have been measured in the x-ray energy range 2100-5600 eV and compared with the calculations of the ＇Darwin Prins＇ and the ＇Mosaic＇ models. It is shown that the integral diffraction coefficients of these crystals are between the calculations of the ＇Darwin Prins＇ and the ＇Mosaic＇ models, but more close to the ＇Darwin Prins＇ model calculations.
文摘In this paper, Kirchhoff formula has been transformed from surface integral form into a line integral form. The new form of the formula can be applied to separate geometrical optical field from diffraction field, and reduce the time of numerical computation greatly. Based on the new form, an analytical formula of diffraction field in the far zone has been presented for the polygonal aperture illuminated by a uniform plane wave.
文摘A Fast Multipole Method (FMM) is developed as a numerical approach to the reduction of the computational cost and requirement memory capacity for a large in solving large-scale problems. In this paper it is applied to the boundary integral equation method (BIEM) for current diffraction from arbitrary 3D bodies. The boundary integral equation is discretized by higher order elements, the FMM is applied to avoid the matrix/vector product, and the resulting algebraic equation is solved by the Generalized Conjugate Residual method (GCR). Numerical examination shows that the FMM is more efficient than the direct evaluation method in computational cost and storage of computers.
基金supported by National Natural Science Foundation of China (50978183)
文摘Diffraction of plane P waves around an alluvial valley of arbitrary shape in poroelastic half-space is investigated by using an indirect boundary integral equation method. Based on the Green＇s fimctions of line source in poroelastic half-space, the scattered waves are constructed using the fictitious wave sources close to the interface of the valley and the density of ficti- tious wave sources are determined by boundary conditions. The precision of the method is verified by the satisfaction extent of boundary conditions, and the comparison between the degenerated solutions and available results in single-phase case. Finally, the nature of diffraction of plane P waves around an alluvial valley in poroelastic half-space is investigated in detail through nu- merical examples.
文摘This review analyzes following numerical methods of a solution of problems of a sound diffraction on ideal and elastic scatterers of a non-analytical form: a method of integral equations, a method of Green’s functions, a method of finite elements, a boundary elements method, a method of Kupradze, a T-matrix method and a method of a geometrical theory of a diffraction.
基金Supported by the National Natural Science Foundation of China(NSFC 62105100)the National Key research and development program in the 14th five year plan(2021YFA1200700)。
文摘Silicon(Si)diffraction microlens arrays are usually used to integrating with infrared focal plane arrays(IRFPAs)to improve their performance.The errors of lithography are unavoidable in the process of the Si diffrac-tion microlens arrays preparation in the conventional engraving method.It has a serious impact on its performance and subsequent applications.In response to the problem of errors of Si diffraction microlens arrays in the conven-tional method,a novel self-alignment method for high precision Si diffraction microlens arrays preparation is pro-posed.The accuracy of the Si diffractive microlens arrays preparation is determined by the accuracy of the first li-thography mask in the novel self-alignment method.In the subsequent etching,the etched area will be protected by the mask layer and the sacrifice layer or the protective layer.The unprotection area is carved to effectively block the non-etching areas,accurately etch the etching area required,and solve the problem of errors.The high precision Si diffraction microlens arrays are obtained by the novel self-alignment method and the diffraction effi-ciency could reach 92.6%.After integrating with IRFPAs,the average blackbody responsity increased by 8.3%,and the average blackbody detectivity increased by 10.3%.It indicates that the Si diffraction microlens arrays can improve the filling factor and reduce crosstalk of IRFPAs through convergence,thereby improving the perfor-mance of the IRFPAs.The results are of great reference significance for improving their performance through opti-mizing the preparation level of micro nano devices.
文摘The development of an in-house computer program for determining the motions and loads of advancing ships through sea waves in the frequency domain,is described in this paper.The code is based on the potential flow formulation and originates from a double-body code enhanced with the regular part of the velocity potential computed using the pulsing source Green function.The code is fully developed in C++language with extensive use of the object-oriented paradigm.The code is capable of estimating the excitation and inertial radiation loads or arbitrary incoming wave frequencies and incidence angles.The hydrodynamic responses such as hydrodynamic coefficients,ship motions,the vertical shear force and the vertical bending moment are estimated.A benchmark container ship and an LNG carrier are selected for testing and validating the computer code.The obtained results are compared with the available experimental data which demonstrate the acceptable compliance for the zero speed whereas there are some discrepancies over the range of frequencies for the advancing ship in different heading angles.
文摘Graphene oxide （GO） ultrathin flat lenses have provided a new and viable solution to achieve high resolution, high efficiency, ultra-light weight, integratable and flexible optical systems. Current GO lenses are designed based on the Fresnel diffraction model, which uses a paraxial approximation for low numerical aperture （NA） focusing process. Herein we develop a lens design method based on the Rayleigh-Sommerfeld （RS） diffraction theory that is able to unambiguously determine the radii of each ring without the optimization process for the first time. More importantly, the RS design method is able to accurately design GO lenses with arbitrary NA and focal length. Our design is experimentally confirmed by fabricating high NA GO lenses with both short and long focal lengths. Compared with the conventional Fresnel design methods, the differences in ring positions and the resulted focal length are up to 13.9% and 9.1%, respectively. Our method can be further applied to design high performance flat lenses of arbitrary materials given the NA and focal length requirements, including metasurfaces or other two-dimensional materials.
文摘This paper is concerned with the problem of scattering of time-harmonic electromag- netic waves from penetrable diffraction gratings in the 2D polarization case. We propose a new, weakly singular, integral equation formulation for the scattering problem which is proved to be uniquely solvable. A main feature of the new integral equation formula- tion is that it avoids the computation of the normal derivative of double-layer potentials which is difficult and time consuming. A fast numerical algorithm is also developed for the scattering problem, based on the NystrSm method for the new integral equation. Nu- merical examples are also shown to illustrate the applicability of the new integral equation formulation.
基金the German Federal Ministry of Education and Research under Grant number 13N8478.
文摘The integral equation method for the simulation of the diffraction by optical gratings is an efficient numerical tool if profile gratings determined by simple crosssection curves are considered.This method in its recent version is capable to tackle profile curves with corners,gratings with thin coated layers,and diffraction scenarios with unfavorably large ratio period over wavelength.We discuss special implementational issues including the efficient evaluation of the quasi-periodic Green kernels,the quadrature algorithm,and the iterative solution of the arising systems of linear equations.Finally,as an example we present the simulation of echelle gratings which demonstrates the efficency of our approach.
基金This work was supported by the National Natural Science Foundation of China (No. 60008002)the Key Oriental Project of Chinese Academy of Sciences (KGCX_2-SW-110)the Shanghai Optical-Tech Special Project (01DJGK015).
文摘A novel method to generate a collimated hollow-laser-beam (HLB) by only a single axicon is proposed. With some reasonable assumptions, the radial light intensity distribution is calculated in detail by diffraction integral theory.The result of numerical simulation shows that this method is valid.Compared with other methods of generating HLB,this scheme is extraordinarily simple in principle and can be utilized experimentally to construct a light trap in atomic fountain for convenience.
文摘In this letter, we propose a method for the numerical calculations of the femtosecond laser pulse passed through a subwavelength aperture. The time-dependent laser pulse is decomposed into a series of monochromatic simple harmonic waves. For the light field of the harmonic wave with a single frequency, the numerical calculation is made based on the solution of the Green's integral equation set of the electromagnetic waves. Such numerical solution is iterated for all the waves with different frequencies, and all the numerical solutions are transformed into the light fields in the time domain by inverse Fourier transform. The light intensity distributions transmitted the subwavelength aperture are calculated and the results show the propagation of the light field is along the direction of the medium interface.
