With the development of ultra-wide coverage technology,multibeam echo-sounder(MBES)system has put forward higher requirements for localization accuracy and computational efficiency of ray tracing method.The classical ...With the development of ultra-wide coverage technology,multibeam echo-sounder(MBES)system has put forward higher requirements for localization accuracy and computational efficiency of ray tracing method.The classical equivalent sound speed profile(ESSP)method replaces the measured sound velocity profile(SVP)with a simple constant gradient SVP,reducing the computational workload of beam positioning.However,in deep-sea environment,the depth measurement error of this method rapidly increases from the central beam to the edge beam.By analyzing the positioning error of the ESSP method at edge beam,it is discovered that the positioning error increases monotonically with the incident angle,and the relationship between them could be expressed by polynomial function.Therefore,an error correction algorithm based on polynomial fitting is obtained.The simulation experiment conducted on an inclined seafloor shows that the proposed algorithm exhibits comparable efficiency to the original ESSP method,while significantly improving bathymetry accuracy by nearly eight times in the edge beam.展开更多
With the techniques provided by the complex my method,like co,nplex ray expansion, complex mytracing,complex my paraxial approximation, etc.,the back scatter of dihedral corner reflectors is evaluated.Numerical result...With the techniques provided by the complex my method,like co,nplex ray expansion, complex mytracing,complex my paraxial approximation, etc.,the back scatter of dihedral corner reflectors is evaluated.Numerical results show that this method gives good and useful RCS prediction of the targets.展开更多
The radiated fields from Gaussian Distributed Aperture (GDA) can be approaximately analyzed bythe Complex Ray Method (CRM). In this article, the radiated fields derived in terms of plana wave spectrumintegral are us...The radiated fields from Gaussian Distributed Aperture (GDA) can be approaximately analyzed bythe Complex Ray Method (CRM). In this article, the radiated fields derived in terms of plana wave spectrumintegral are used as standard, and the error source of the replacement by complex ray fields is found. Withdiscussions on the error characteristics of this CIUN approximation, conclusions are drawn which are expected to be helpful in applications of CRM.展开更多
Based on the method of reverberation ray matrix（MRRM）, a reverberation matrix for planar framed structures composed of anisotropic Timoshenko（T） beam members containing completely hinged joints is developed for st...Based on the method of reverberation ray matrix（MRRM）, a reverberation matrix for planar framed structures composed of anisotropic Timoshenko（T） beam members containing completely hinged joints is developed for static analysis of such structures.In the MRRM for dynamic analysis, amplitudes of arriving and departing waves for joints are chosen as unknown quantities. However, for the present case of static analysis, displacements and rotational angles at the ends of each beam member are directly considered as unknown quantities. The expressions for stiffness matrices for anisotropic beam members are developed. A corresponding reverberation matrix is derived analytically for exact and unified determination on the displacements and internal forces at both ends of each member and arbitrary cross sectional locations in the structure. Numerical examples are given and compared with the finite element method（FEM） results to validate the present model. The characteristic parameter analysis is performed to demonstrate accuracy of the present model with the T beam theory in contrast with errors in the usual model based on the Euler-Bernoulli（EB） beam theory. The resulting reverberation matrix can be used for exact calculation of anisotropic framed structures as well as for parameter analysis of geometrical and material properties of the framed structures.展开更多
A procedure of the method of reverberation ray matrix(MRRM)is developed to perform the buckling analysis of thin multi-span rectangular plates having internal line supports or stiffeners.A computation algorithm for th...A procedure of the method of reverberation ray matrix(MRRM)is developed to perform the buckling analysis of thin multi-span rectangular plates having internal line supports or stiffeners.A computation algorithm for the reverberation ray matrix in the MRRM is derived to determine the buckling loading.Specifically,the analytical solutions are presented for the buckling of the structure having two opposite simply-supported or clamped-supported edges with spans,while the constraint condition of two remaining edges may be in any combination of free,simply-supported,and clamped boundary conditions.Furthermore,based on the analysis of matrices relating to the unknown coefficients in the solution form for the deflection in terms of buckling modal functions,some recursive equations(REs)for the MRRM are introduced to generate a reduced reverberation ray matrix with unchanged dimension when the number of spans increases,which promotes the computation efficiency.Several numerical examples are given,and the present results are compared with the known solutions to illustrate the validity and accurateness of the MRRM for the buckling analysis.展开更多
The power-law frequency distributions of the peak flux of solar flare X-ray emission have been studied extensively and attributed to a system having self-organized criticality （SOC）. In this paper, we first show tha...The power-law frequency distributions of the peak flux of solar flare X-ray emission have been studied extensively and attributed to a system having self-organized criticality （SOC）. In this paper, we first show that, so long as the shape of the normalized light curve is not correlated with the peak flux, the flux histogram of solar flares also follows a power-law distribution with the same spectral index as the power- law frequency distribution of the peak flux, which may partially explain why power-law distributions are ubiquitous in the Universe. We then show that the spectral indexes of the histograms of soft X-ray fluxes observed by GOES satellites in two different energy channels are different： the higher energy channel has a harder distribution than the lower energy channel, which challenges the universal power-law distribution predicted by SOC models and implies a very soft distribution of thermal energy content of plasmas probed by the GOES satellites. The temperature （T） distribution, on the other hand, approaches a power-law dis- tribution with an index of 2 for high values of T. Hence the application of SOC models to the statistical properties of solar flares needs to be revisited.展开更多
With the assistance of the identified/associated sources in the second Fermi Large Area Telescope （LAT） catalog, we analyze and resolve the spatial distribution and the distributions of the gamma-ray spectral and va...With the assistance of the identified/associated sources in the second Fermi Large Area Telescope （LAT） catalog, we analyze and resolve the spatial distribution and the distributions of the gamma-ray spectral and variability indices of the remaining 575 unassociated Fermi LAT sources. Consequently, it is suggested that the unassociated sources could statistically consist of Galactic supernova remnants/pulsar wind nebulae, BL Lacertae objects, fiat spectrum radio quasars and other types of active galaxies with fractions of 25%, 29%, 41% and 5%, respectively.展开更多
Ray tracing method is used to study the propagation of collimated beams in a liquid-core cylindrical lens(LCL),which has dual functions of diffusion cell and image formation.The diffusion images on the focal plane of ...Ray tracing method is used to study the propagation of collimated beams in a liquid-core cylindrical lens(LCL),which has dual functions of diffusion cell and image formation.The diffusion images on the focal plane of the used LCL are simulated by establishing and solving both linear and nonlinear ray equations,the calculated results indicate that the complex imaging results of LCL in inhomogeneous media can be treated by the law of ray propagation in homogeneous media under the condition of small refractive index gradient of diffusion solution.Guided by the calculation conditions,the diffusion process of triethylene glycol aqueous solution is experimentally studied at room temperature by using the LCL in this paper.The spatial and temporal concentration profile Ce(z,t)of diffusion solution is obtained by analyzing diffusion image appearing on the focal plane of the LCL;Then,the concentration-dependent diffusion coefficient is assumed to be a polynomial D(C)=D0×(1+α1C+α2C2+α3C3+…).The finite difference method is used to solve the Fick diffusion equation for calculating numerically the concentration profiles Cn(z,t).The D(C)of triethylene glycol aqueous solution is obtained by comparing the Cn(z,t)with Ce(z,t).Finally,the obtained polynomial D(C)is used to calculate the refractive index profiles nn(z,t)s of diffusion solution in the used LCL.Based on the ray propagation law in inhomogeneous media and the calculated n(z,t),the ray tracing method is used again to simulate the dynamic images of the whole experimental diffusion process to varify the correctness of the calculated D(C).The method presented in this work opens up a new way for both measuring and verifying the concentration-dependent liquid diffusion coefficients.展开更多
The study of marine data visualization is of great value. Marine data, due to its large scale, random variation and multiresolution in nature, are hard to be visualized and analyzed. Nowadays, constructing an ocean mo...The study of marine data visualization is of great value. Marine data, due to its large scale, random variation and multiresolution in nature, are hard to be visualized and analyzed. Nowadays, constructing an ocean model and visualizing model results have become some of the most important research topics of ‘Digital Ocean'. In this paper, a spherical ray casting method is developed to improve the traditional ray-casting algorithm and to make efficient use of GPUs. Aiming at the ocean current data, a 3D view-dependent line integral convolution method is used, in which the spatial frequency is adapted according to the distance from a camera. The study is based on a 3D virtual reality and visualization engine, namely the VV-Ocean. Some interactive operations are also provided to highlight the interesting structures and the characteristics of volumetric data. Finally, the marine data gathered in the East China Sea are displayed and analyzed. The results show that the method meets the requirements of real-time and interactive rendering.展开更多
We investigate the relationship between the rise width and the full width of gamma-ray burst pulses. Theoretical analysis shows that either width is proportional to Г^-2△τθ,FWHMRc/c（Г the Lorentz factor of the b...We investigate the relationship between the rise width and the full width of gamma-ray burst pulses. Theoretical analysis shows that either width is proportional to Г^-2△τθ,FWHMRc/c（Г the Lorentz factor of the bulk motion, △τθ,FWHM a local pulse＇s width, Rc the radius of fireballs and c the velocity of light）. We study the relationship for four samples of observed pulses. We find： （1） merely the curvature effect could reproduce the relationship between the rise and full widths with the same slope as derived from the model of Qin et al.; （2） gamma-ray burst pulses, selected from both the short and long GRBs, follow the same sequence in the rise width vs. full width diagram, with the shorter pulses at one end; （3） all GRBs may intrinsically result from local Gaussian pulses. These features place constraints on the physical mechanism（s） for producing long and short GRBs.展开更多
We performed detailed time-resolved spectroscopy of bright tong gamma- ray bursts （GRBs） which show significant GeV emissions （GRB 080916C, GRB 090902B and GRB 090926A）. In addition to the standard Band model, we ...We performed detailed time-resolved spectroscopy of bright tong gamma- ray bursts （GRBs） which show significant GeV emissions （GRB 080916C, GRB 090902B and GRB 090926A）. In addition to the standard Band model, we also use a model consisting of a black body and a power law to fit the spectra. We find that for the latter model there are indications of an additional soft component in the spectra. While previous studies have shown that such models are required for GRB 090902B, here we find that a composite spectral model consisting of two blackbodies and a power law adequately fits the data of all the three bright GRBs. We investigate the evolution of the spectral parameters and find several interesting features that appear in all three GRBs, like （a） temperatures of the blackbodies are strongly correlated with each other, （b） fluxes in the black body components are strongly correlated with each other, （c） the temperatures of the black body trace the profile of the individual pulses of the GRBs, and （d） the characteristics of power law components like the spectral index and the delayed onset bear a close similarity to the emission characteristics in the GeV regions. We discuss the implications of these results and the possibility of identifying the radiation mechanisms during the prompt emission of GRBs.展开更多
The relativistic external shock model of gamma-ray burst （GRB） afterglows has been estab- lished with five free parameters, i.e., the total kinetic energy E, the equipartition parameters for electrons ee and for the...The relativistic external shock model of gamma-ray burst （GRB） afterglows has been estab- lished with five free parameters, i.e., the total kinetic energy E, the equipartition parameters for electrons ee and for the magnetic field eB, the number density of the environment n and the index of the power- law distribution of shocked electrons p. A lot of modified models have been constructed to consider the variety of GRB afterglows, such as： the wind medium environment by letting n change with radius, the energy injection model by letting kinetic energy change with time and so on. In this paper, by as- suming all four parameters （except p） change with time, we obtain a set of formulas for the dynamics and radiation, which can be used as a reference for modeling GRB afterglows. Some interesting results are obtained. For example, in some spectral segments, the radiated flux density does not depend on the number density or the profile of the environment. As an application, through modeling the afterglow of GRB 060607A, we find that it can be interpreted in the framework of the time dependent parameter model within a reasonable range.展开更多
We analyze the emission plateaus in the X-ray afterglow light curves of gamma-ray bursts （GRBs） and those in the optical light curves of type Ⅱ plateau su- pernovae （SNe Ⅱ-P） in order to study whether they have ...We analyze the emission plateaus in the X-ray afterglow light curves of gamma-ray bursts （GRBs） and those in the optical light curves of type Ⅱ plateau su- pernovae （SNe Ⅱ-P） in order to study whether they have similar late energy injection behaviors. We show that correlations of bolometric energies （or luminosities） between the prompt explosions and the plateaus for the two phenomena are similar. The energy emitted by SNe II-P are at the lower end of the range of possible energies for GRBs. The bolometric energies （or luminosities） in the prompt phase Eexpl （or Lexpl） and in the plateau phase E_plateau （or L_plateau） share relations of E_expl ∝E _0.73±0.14_plateau and L_expl ∝ L^-0.70_plateau. These results may indicate a similar late energy injection behavior that produces the observed plateaus in these two phenomena.展开更多
To help optimize the spacecraft design and reduce the risk of spacecraft mission failure, a new approach to assess the survivability of spacecraft in orbit is presented here, including the following three steps: 1)Sen...To help optimize the spacecraft design and reduce the risk of spacecraft mission failure, a new approach to assess the survivability of spacecraft in orbit is presented here, including the following three steps: 1)Sensitivity Analysis of spacecraft. A new sensitivity analysis method, a ray method based on virtual outer wall, is presented here. Using rays to simulate the debris cloud can effectively address the component shadowing issues. 2) Component Vulnerability analysis of spacecraft. A function "Component functional reduction degree-Component physical damage degree" is provided here to clearly describe the component functional reduction. 3) System-level Survivability Assessment of spacecraft. A new method based on expert knowledge reasoning, instead of traditional artificial failure tree method, is presented here to greatly improve the efficiency and accuracy of calculation.展开更多
The generalized ray method(GRM) has been successfully used to study the transient elastic wave transmitting in the beams,planar trusses,space frames and infinite layered media.In this letter,the GRM is extended to inv...The generalized ray method(GRM) has been successfully used to study the transient elastic wave transmitting in the beams,planar trusses,space frames and infinite layered media.In this letter,the GRM is extended to investigate the early short time transient responses of laminated composite cylindrical shells under impact load.By using the Laplace transformation and referring to the boundary conditions,the ray groups transmitting in the finite laminated cylindrical shells under the shock load are obtained and the transient response related to each ray group can be derived via FFT algorithm.From the numerical results,it is shown that the early short time transient accelerations of the laminated composite cylindrical shell under impact loads are very large.But the short time transient shear strain and displacement are very small.展开更多
Using the ray trace method, three-section semiconductor lasers are studied. An analytic expression of output power for the three-section semiconductor lasers is derived for the first time. From this expression, thresh...Using the ray trace method, three-section semiconductor lasers are studied. An analytic expression of output power for the three-section semiconductor lasers is derived for the first time. From this expression, threshold condition is also obtained.展开更多
Multielectrode semiconductor lasers are studied via the ray method.The expression of the output photon number of N -electrode semiconductor lasers has been derived for the first time.When N =1 or 2,the expressio...Multielectrode semiconductor lasers are studied via the ray method.The expression of the output photon number of N -electrode semiconductor lasers has been derived for the first time.When N =1 or 2,the expression of the output photon number fits in that of one-electrode (general) or two-electrode semiconductor lasers perfectly.展开更多
Known as laser trapping,optical tweezers,with nanometer accuracy and pico-newton precision,plays a pivotal role in single bio-molecule measurements and controllable motions of micro-machines.In order to advance the fl...Known as laser trapping,optical tweezers,with nanometer accuracy and pico-newton precision,plays a pivotal role in single bio-molecule measurements and controllable motions of micro-machines.In order to advance the flourishing applications for those achievements,it is necessary to make clear the three-dimensional dynamic process of micro-particles stepping into an optical field.In this paper,we utilize the ray optics method to calculate the optical force and optical torque of a micro-sphere in optical tweezers.With the influence of viscosity force and torque taken into account,we numerically solve and analyze the dynamic process of a dielectric micro-sphere in optical tweezers on the basis of Newton mechanical equations under various conditions of initial positions and velocity vectors of the particle.The particle trajectory over time can demonstrate whether the particle can be successfully trapped into the optical tweezers center and reveal the subtle details of this trapping process.Even in a simple pair of optical tweezers,the dielectric micro-sphere exhibits abundant phases of mechanical motions including acceleration,deceleration,and turning.These studies will be of great help to understand the particle-laser trap interaction in various situations and promote exciting possibilities for exploring novel ways to control the mechanical dynamics of microscale particles.展开更多
Pipelines are an important part of urban infrastructure development.As part of a virtual globe(VG),the high-efficiency and high-quality visualization of 3D large-scale and high-density urban pipelines is of great impo...Pipelines are an important part of urban infrastructure development.As part of a virtual globe(VG),the high-efficiency and high-quality visualization of 3D large-scale and high-density urban pipelines is of great importance.This paper proposes a GPU-based pipeline ray casting method for the visualization of urban-scale pipelines in the framework of a VG.The method involves the initial partitioning of the pipeline data into tiles,based on the relationship between the pipeline layer scale and the discrete global grid system(DGGSs).The pipeline centerline in each tile is then segmented and encoded,and a coarser pipeline bounding volume is subsequently constructed using a geometry shader.Finally,the fine 3D pipeline is rendered using a pixel shader.The results of the experimental implementation of the proposed method show that it satisfies the requirements for the multiscale visualization of pipelines in a VG.Moreover,compared with the traditional polygon-based method,the method facilitates a 20%increase in rendering frame rate for the same pixel level accuracy display effect.It also enables the visualization of the thickness of the 3D pipeline without any obvious effect on the rendering efficiency。展开更多
The space partitioning algorithm based on the rounding and addressing operations has been proved to be an efficient space partitioning algorithm with the potential for real-time calculation.An improvement on this kind...The space partitioning algorithm based on the rounding and addressing operations has been proved to be an efficient space partitioning algorithm with the potential for real-time calculation.An improvement on this kind of space partitioning algorithms for solving complex 3D models is presented.Numerical examples show that the efficiency of the improved algorithm is better than that of the original method.When the size of most target elements is smaller than the size of spatial grids,the efficiency of the improved method can be more than four times of that of the original method.An adaptive method of space partitioning based on the improved algorithm is developed by taking the surface element density or the curvature as the threshold for deep partitioning and conducting the deep partitioning using the octree method.A computer program implementation for applying the method in some typical applications is discussed,and the performance in terms of the efficiency,reliability,and resource use is evaluated.Application testing shows that the results of the adaptive spacing partitioning are more convenient for the follow-up use than that of the basic uniform space partitioning.Furthermore,when it is used to calculate the electromagnetic scattering of complex targets by the ray tracing(RT)method,the adaptive space partitioning algorithm can reduce the calculation time of the RT process by more than 40%compared with the uniform space segmentation algorithm.展开更多
基金The Natural Science Foundation of Shandong Province of China under contract Nos ZR2022MA051 and ZR2020MA090the National Natural Science Foundation of China under contract No.U22A2012+2 种基金China Postdoctoral Science Foundation under contract No.2020M670891the SDUST Research Fund under contract No.2019TDJH103the Talent Introduction Plan for Youth Innovation Team in universities of Shandong Province(innovation team of satellite positioning and navigation)。
文摘With the development of ultra-wide coverage technology,multibeam echo-sounder(MBES)system has put forward higher requirements for localization accuracy and computational efficiency of ray tracing method.The classical equivalent sound speed profile(ESSP)method replaces the measured sound velocity profile(SVP)with a simple constant gradient SVP,reducing the computational workload of beam positioning.However,in deep-sea environment,the depth measurement error of this method rapidly increases from the central beam to the edge beam.By analyzing the positioning error of the ESSP method at edge beam,it is discovered that the positioning error increases monotonically with the incident angle,and the relationship between them could be expressed by polynomial function.Therefore,an error correction algorithm based on polynomial fitting is obtained.The simulation experiment conducted on an inclined seafloor shows that the proposed algorithm exhibits comparable efficiency to the original ESSP method,while significantly improving bathymetry accuracy by nearly eight times in the edge beam.
文摘With the techniques provided by the complex my method,like co,nplex ray expansion, complex mytracing,complex my paraxial approximation, etc.,the back scatter of dihedral corner reflectors is evaluated.Numerical results show that this method gives good and useful RCS prediction of the targets.
文摘The radiated fields from Gaussian Distributed Aperture (GDA) can be approaximately analyzed bythe Complex Ray Method (CRM). In this article, the radiated fields derived in terms of plana wave spectrumintegral are used as standard, and the error source of the replacement by complex ray fields is found. Withdiscussions on the error characteristics of this CIUN approximation, conclusions are drawn which are expected to be helpful in applications of CRM.
基金Project supported by the Program for New Century Excellent Talents in Universities(NCET)by the Ministry of Education of China(No.NCET-04-0373)
文摘Based on the method of reverberation ray matrix（MRRM）, a reverberation matrix for planar framed structures composed of anisotropic Timoshenko（T） beam members containing completely hinged joints is developed for static analysis of such structures.In the MRRM for dynamic analysis, amplitudes of arriving and departing waves for joints are chosen as unknown quantities. However, for the present case of static analysis, displacements and rotational angles at the ends of each beam member are directly considered as unknown quantities. The expressions for stiffness matrices for anisotropic beam members are developed. A corresponding reverberation matrix is derived analytically for exact and unified determination on the displacements and internal forces at both ends of each member and arbitrary cross sectional locations in the structure. Numerical examples are given and compared with the finite element method（FEM） results to validate the present model. The characteristic parameter analysis is performed to demonstrate accuracy of the present model with the T beam theory in contrast with errors in the usual model based on the Euler-Bernoulli（EB） beam theory. The resulting reverberation matrix can be used for exact calculation of anisotropic framed structures as well as for parameter analysis of geometrical and material properties of the framed structures.
文摘A procedure of the method of reverberation ray matrix(MRRM)is developed to perform the buckling analysis of thin multi-span rectangular plates having internal line supports or stiffeners.A computation algorithm for the reverberation ray matrix in the MRRM is derived to determine the buckling loading.Specifically,the analytical solutions are presented for the buckling of the structure having two opposite simply-supported or clamped-supported edges with spans,while the constraint condition of two remaining edges may be in any combination of free,simply-supported,and clamped boundary conditions.Furthermore,based on the analysis of matrices relating to the unknown coefficients in the solution form for the deflection in terms of buckling modal functions,some recursive equations(REs)for the MRRM are introduced to generate a reduced reverberation ray matrix with unchanged dimension when the number of spans increases,which promotes the computation efficiency.Several numerical examples are given,and the present results are compared with the known solutions to illustrate the validity and accurateness of the MRRM for the buckling analysis.
基金supported partially by the Strategic Priority Research Program,the Emergence of Cosmological Structures,of the Chinese Academy of Sciences(Grant No.XDB09000000)MSTC Program2011 CB811402+3 种基金the National Natural Science Foundation of China(NSFC)(Grant Nos.11173063, 11173064,11233008 and 11427803)supported by the NFSC(Grant No.11473070)the Natural Science Foundation of Jiangsu Province (Grant BK2012889)the Youth Innovation Promotion Association,CAS,for financial support
文摘The power-law frequency distributions of the peak flux of solar flare X-ray emission have been studied extensively and attributed to a system having self-organized criticality （SOC）. In this paper, we first show that, so long as the shape of the normalized light curve is not correlated with the peak flux, the flux histogram of solar flares also follows a power-law distribution with the same spectral index as the power- law frequency distribution of the peak flux, which may partially explain why power-law distributions are ubiquitous in the Universe. We then show that the spectral indexes of the histograms of soft X-ray fluxes observed by GOES satellites in two different energy channels are different： the higher energy channel has a harder distribution than the lower energy channel, which challenges the universal power-law distribution predicted by SOC models and implies a very soft distribution of thermal energy content of plasmas probed by the GOES satellites. The temperature （T） distribution, on the other hand, approaches a power-law dis- tribution with an index of 2 for high values of T. Hence the application of SOC models to the statistical properties of solar flares needs to be revisited.
基金supported by the National Natural Science Foundation of China (Grant No. 11103004)the Foundation for the Authors of National Excellent Doctoral Dissertations of China (Grant No. 201225)
文摘With the assistance of the identified/associated sources in the second Fermi Large Area Telescope （LAT） catalog, we analyze and resolve the spatial distribution and the distributions of the gamma-ray spectral and variability indices of the remaining 575 unassociated Fermi LAT sources. Consequently, it is suggested that the unassociated sources could statistically consist of Galactic supernova remnants/pulsar wind nebulae, BL Lacertae objects, fiat spectrum radio quasars and other types of active galaxies with fractions of 25%, 29%, 41% and 5%, respectively.
基金the National Natural Science Foundation of China(Grant No.11804296)the Joint Key Project of Yunnan Province,China(Grant Nos.2018FY001-020 and 2018ZI002)the Fund from the Educational Department of Yunnan Province,China(Grant No.2016CYH05).
文摘Ray tracing method is used to study the propagation of collimated beams in a liquid-core cylindrical lens(LCL),which has dual functions of diffusion cell and image formation.The diffusion images on the focal plane of the used LCL are simulated by establishing and solving both linear and nonlinear ray equations,the calculated results indicate that the complex imaging results of LCL in inhomogeneous media can be treated by the law of ray propagation in homogeneous media under the condition of small refractive index gradient of diffusion solution.Guided by the calculation conditions,the diffusion process of triethylene glycol aqueous solution is experimentally studied at room temperature by using the LCL in this paper.The spatial and temporal concentration profile Ce(z,t)of diffusion solution is obtained by analyzing diffusion image appearing on the focal plane of the LCL;Then,the concentration-dependent diffusion coefficient is assumed to be a polynomial D(C)=D0×(1+α1C+α2C2+α3C3+…).The finite difference method is used to solve the Fick diffusion equation for calculating numerically the concentration profiles Cn(z,t).The D(C)of triethylene glycol aqueous solution is obtained by comparing the Cn(z,t)with Ce(z,t).Finally,the obtained polynomial D(C)is used to calculate the refractive index profiles nn(z,t)s of diffusion solution in the used LCL.Based on the ray propagation law in inhomogeneous media and the calculated n(z,t),the ray tracing method is used again to simulate the dynamic images of the whole experimental diffusion process to varify the correctness of the calculated D(C).The method presented in this work opens up a new way for both measuring and verifying the concentration-dependent liquid diffusion coefficients.
基金supported by the Natural Science Foundation of China under Project 41076115the Global Change Research Program of China under project 2012CB955603the Public Science and Technology Research Funds of the Ocean under project 201005019
文摘The study of marine data visualization is of great value. Marine data, due to its large scale, random variation and multiresolution in nature, are hard to be visualized and analyzed. Nowadays, constructing an ocean model and visualizing model results have become some of the most important research topics of ‘Digital Ocean'. In this paper, a spherical ray casting method is developed to improve the traditional ray-casting algorithm and to make efficient use of GPUs. Aiming at the ocean current data, a 3D view-dependent line integral convolution method is used, in which the spatial frequency is adapted according to the distance from a camera. The study is based on a 3D virtual reality and visualization engine, namely the VV-Ocean. Some interactive operations are also provided to highlight the interesting structures and the characteristics of volumetric data. Finally, the marine data gathered in the East China Sea are displayed and analyzed. The results show that the method meets the requirements of real-time and interactive rendering.
基金Supported by the National Natural Science Foundation of China
文摘We investigate the relationship between the rise width and the full width of gamma-ray burst pulses. Theoretical analysis shows that either width is proportional to Г^-2△τθ,FWHMRc/c（Г the Lorentz factor of the bulk motion, △τθ,FWHM a local pulse＇s width, Rc the radius of fireballs and c the velocity of light）. We study the relationship for four samples of observed pulses. We find： （1） merely the curvature effect could reproduce the relationship between the rise and full widths with the same slope as derived from the model of Qin et al.; （2） gamma-ray burst pulses, selected from both the short and long GRBs, follow the same sequence in the rise width vs. full width diagram, with the shorter pulses at one end; （3） all GRBs may intrinsically result from local Gaussian pulses. These features place constraints on the physical mechanism（s） for producing long and short GRBs.
文摘We performed detailed time-resolved spectroscopy of bright tong gamma- ray bursts （GRBs） which show significant GeV emissions （GRB 080916C, GRB 090902B and GRB 090926A）. In addition to the standard Band model, we also use a model consisting of a black body and a power law to fit the spectra. We find that for the latter model there are indications of an additional soft component in the spectra. While previous studies have shown that such models are required for GRB 090902B, here we find that a composite spectral model consisting of two blackbodies and a power law adequately fits the data of all the three bright GRBs. We investigate the evolution of the spectral parameters and find several interesting features that appear in all three GRBs, like （a） temperatures of the blackbodies are strongly correlated with each other, （b） fluxes in the black body components are strongly correlated with each other, （c） the temperatures of the black body trace the profile of the individual pulses of the GRBs, and （d） the characteristics of power law components like the spectral index and the delayed onset bear a close similarity to the emission characteristics in the GeV regions. We discuss the implications of these results and the possibility of identifying the radiation mechanisms during the prompt emission of GRBs.
基金supported by the National Basic Research Program of China(973 Program,Grant 2014CB845800)the National Natural Science Foundation of China(Grant Nos.11773010,U1431124 and U1738132)
文摘The relativistic external shock model of gamma-ray burst （GRB） afterglows has been estab- lished with five free parameters, i.e., the total kinetic energy E, the equipartition parameters for electrons ee and for the magnetic field eB, the number density of the environment n and the index of the power- law distribution of shocked electrons p. A lot of modified models have been constructed to consider the variety of GRB afterglows, such as： the wind medium environment by letting n change with radius, the energy injection model by letting kinetic energy change with time and so on. In this paper, by as- suming all four parameters （except p） change with time, we obtain a set of formulas for the dynamics and radiation, which can be used as a reference for modeling GRB afterglows. Some interesting results are obtained. For example, in some spectral segments, the radiated flux density does not depend on the number density or the profile of the environment. As an application, through modeling the afterglow of GRB 060607A, we find that it can be interpreted in the framework of the time dependent parameter model within a reasonable range.
基金supported by the National Basic Research Program of China (973 program, Grant Nos. 2012CB821800 and 2009CB824800)the National Natural Science Foundation of China(Grant Nos. 11225314, 11103026 and 10935001)
文摘We analyze the emission plateaus in the X-ray afterglow light curves of gamma-ray bursts （GRBs） and those in the optical light curves of type Ⅱ plateau su- pernovae （SNe Ⅱ-P） in order to study whether they have similar late energy injection behaviors. We show that correlations of bolometric energies （or luminosities） between the prompt explosions and the plateaus for the two phenomena are similar. The energy emitted by SNe II-P are at the lower end of the range of possible energies for GRBs. The bolometric energies （or luminosities） in the prompt phase Eexpl （or Lexpl） and in the plateau phase E_plateau （or L_plateau） share relations of E_expl ∝E _0.73±0.14_plateau and L_expl ∝ L^-0.70_plateau. These results may indicate a similar late energy injection behavior that produces the observed plateaus in these two phenomena.
文摘To help optimize the spacecraft design and reduce the risk of spacecraft mission failure, a new approach to assess the survivability of spacecraft in orbit is presented here, including the following three steps: 1)Sensitivity Analysis of spacecraft. A new sensitivity analysis method, a ray method based on virtual outer wall, is presented here. Using rays to simulate the debris cloud can effectively address the component shadowing issues. 2) Component Vulnerability analysis of spacecraft. A function "Component functional reduction degree-Component physical damage degree" is provided here to clearly describe the component functional reduction. 3) System-level Survivability Assessment of spacecraft. A new method based on expert knowledge reasoning, instead of traditional artificial failure tree method, is presented here to greatly improve the efficiency and accuracy of calculation.
基金supported by the National Basic Research Program of China(2011CB711102)the National Natural Science Foundation of China(10672017)
文摘The generalized ray method(GRM) has been successfully used to study the transient elastic wave transmitting in the beams,planar trusses,space frames and infinite layered media.In this letter,the GRM is extended to investigate the early short time transient responses of laminated composite cylindrical shells under impact load.By using the Laplace transformation and referring to the boundary conditions,the ray groups transmitting in the finite laminated cylindrical shells under the shock load are obtained and the transient response related to each ray group can be derived via FFT algorithm.From the numerical results,it is shown that the early short time transient accelerations of the laminated composite cylindrical shell under impact loads are very large.But the short time transient shear strain and displacement are very small.
文摘Using the ray trace method, three-section semiconductor lasers are studied. An analytic expression of output power for the three-section semiconductor lasers is derived for the first time. From this expression, threshold condition is also obtained.
文摘Multielectrode semiconductor lasers are studied via the ray method.The expression of the output photon number of N -electrode semiconductor lasers has been derived for the first time.When N =1 or 2,the expression of the output photon number fits in that of one-electrode (general) or two-electrode semiconductor lasers perfectly.
基金This work is supported by the National Natural Science Foundation of China(Grant No.11974119 and No.11804399)the Guangdong Innovative and Entrepreneurial Research Team Program(Grant No.2016ZT06C594)+1 种基金the Fundamental Research Funds for the Central Universities,South-Central University for Nationalities(Grant No.CZQ20018)National Key R&D Program of China(No.2018YFA 0306200).
文摘Known as laser trapping,optical tweezers,with nanometer accuracy and pico-newton precision,plays a pivotal role in single bio-molecule measurements and controllable motions of micro-machines.In order to advance the flourishing applications for those achievements,it is necessary to make clear the three-dimensional dynamic process of micro-particles stepping into an optical field.In this paper,we utilize the ray optics method to calculate the optical force and optical torque of a micro-sphere in optical tweezers.With the influence of viscosity force and torque taken into account,we numerically solve and analyze the dynamic process of a dielectric micro-sphere in optical tweezers on the basis of Newton mechanical equations under various conditions of initial positions and velocity vectors of the particle.The particle trajectory over time can demonstrate whether the particle can be successfully trapped into the optical tweezers center and reveal the subtle details of this trapping process.Even in a simple pair of optical tweezers,the dielectric micro-sphere exhibits abundant phases of mechanical motions including acceleration,deceleration,and turning.These studies will be of great help to understand the particle-laser trap interaction in various situations and promote exciting possibilities for exploring novel ways to control the mechanical dynamics of microscale particles.
文摘Pipelines are an important part of urban infrastructure development.As part of a virtual globe(VG),the high-efficiency and high-quality visualization of 3D large-scale and high-density urban pipelines is of great importance.This paper proposes a GPU-based pipeline ray casting method for the visualization of urban-scale pipelines in the framework of a VG.The method involves the initial partitioning of the pipeline data into tiles,based on the relationship between the pipeline layer scale and the discrete global grid system(DGGSs).The pipeline centerline in each tile is then segmented and encoded,and a coarser pipeline bounding volume is subsequently constructed using a geometry shader.Finally,the fine 3D pipeline is rendered using a pixel shader.The results of the experimental implementation of the proposed method show that it satisfies the requirements for the multiscale visualization of pipelines in a VG.Moreover,compared with the traditional polygon-based method,the method facilitates a 20%increase in rendering frame rate for the same pixel level accuracy display effect.It also enables the visualization of the thickness of the 3D pipeline without any obvious effect on the rendering efficiency。
基金This work was supported by the National Natural Science Foundation of China(61601015,91538204).
文摘The space partitioning algorithm based on the rounding and addressing operations has been proved to be an efficient space partitioning algorithm with the potential for real-time calculation.An improvement on this kind of space partitioning algorithms for solving complex 3D models is presented.Numerical examples show that the efficiency of the improved algorithm is better than that of the original method.When the size of most target elements is smaller than the size of spatial grids,the efficiency of the improved method can be more than four times of that of the original method.An adaptive method of space partitioning based on the improved algorithm is developed by taking the surface element density or the curvature as the threshold for deep partitioning and conducting the deep partitioning using the octree method.A computer program implementation for applying the method in some typical applications is discussed,and the performance in terms of the efficiency,reliability,and resource use is evaluated.Application testing shows that the results of the adaptive spacing partitioning are more convenient for the follow-up use than that of the basic uniform space partitioning.Furthermore,when it is used to calculate the electromagnetic scattering of complex targets by the ray tracing(RT)method,the adaptive space partitioning algorithm can reduce the calculation time of the RT process by more than 40%compared with the uniform space segmentation algorithm.
