An improved algorithm based on equivalent sound velocity profile method at large incident angle

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.

Method of reverberation ray matrix for static analysis of planar framed structures composed of anisotropic Timoshenko beam members

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 for the buckling analysis of a thin multi-span plate

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.

Indirect-approach method for specified-end points seismic ray tracing in three dimen-sional inhomogeneous media

Generally speaking, the factors of both medium and tectonic that give rise to heterogeneity of the earth crust and mantle structures should be taken into account simultaneously in three dimensional seismic ray tracing. In this paper, the three dimensional structure models are constructed with the model similar to generation system in computer aid design and manufacturing (CAD/CAM). Based on the algorithm proposed by Cerveny et al . for complete seismic ray tracing in complex three dimensional structures, a new technique called the indirect approach method for two point seismic ray tracing in three dimensional laterally heterogeneous media has been put forward, and its numerical computing examples were given.

On the power-law distributions of X-ray fluxes from solar flares observed with GOES

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.

Measurement and verification of concentration-dependent diffusion coefficient:Ray tracing imagery of diffusion process

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 Relationship between the Rise Width and the Full Width of γ-ray Burst Pulses and Its Implications

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.

Time-resolved spectral analysis of prompt emission from long gamma-ray bursts with GeV emission

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.

A statistical classification of the unassociated gamma-ray sources in the second Fermi Large Area Telescope Catalog

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.

Revisiting gamma-ray burst afterglows with time-dependent parameters

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.

X-ray and optical plateaus following the main bursts in GRBs and SNe Ⅱ-P: a hint about similar late injection behaviors?

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.

Survivability assessment of spacecraft impacted by orbit debris

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.

Transient impact responses of laminated composite cylindrical shells

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.

3D dynamic motion of a dielectric micro-sphere within optical tweezers

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.

GPU ray casting method for visualizing 3D pipelines in a virtual globe

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。

Basement interface structural characteristics beneath Jiashi strong earthquake swarm area in Xinjiang, China

The seismic data obtained from high resolution seismic refraction profile in Jiashi strong earthquake swarm area in Xinjiang, China were further processed with ray hit analysis method and more complete basement interface structural characteristics beneath Jiashi strong earthquake swarm area were determined. The results show that there are two clear basement interfaces at the upper crust in Jiashi strong earthquake swarm area. The first one with buried depth ranging from 2.6 km to 3.3 km presents integral and continuous structure, and it appears an inclined plane interface and smoothly rises up toward Tianshan Mountain. The second basement interface with buried depth from 8.5 km to 11.8 km, is the antiquated crystalline basement of Tarim basin. Near the post number of 37 km, the buried depth of the crystalline basement changed abruptly by 2.5 km, which maybe result from an ultra crystalline basement fault. If taking this fault as a boundary, the crystalline basement could be divided into two parts, i.e. the southwestern segment with buried depth about 11.5 km, and the northeastern segment with buried depth approxi-mately from 8.5 km to 9.0 km. That is to say, in each segment, the buried depth changes not too much. The northeast segment rises up as a whole and upheaves slightly from southwest to northeast, which reflects the upper crustal deformation characteristics under the special tectonic background at the northwestern edge of Tarim basin.

An efficient adaptive space partitioning algorithm for electromagnetic scattering calculation of complex 3D models

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.

GRB jet beaming angle statistics

Existing theory and models suggest that a Type I （merger） GRB should have a larger jet beaming angle than a Type II （collapsar） GRB, but so far no statistical evidence is available to support this suggestion. In this paper, we obtain a sample of 37 beaming angles and calculate the probability that this is true. A correction is also devised to account for the scarcity of Type I GRBs in our sample. The probability is calculated to be 83% without the correction and 71% with it.

ANALYSIS AND EXPERIMENT ON TRANSIENT DYNAMIC RESPONSE IN FINITE MINDLIN PLATE

The transient wave propagation in the finite rectangular Mindlin plate is investi- gated by the analytical and experimental methods. The generalized ray method （GRM） which has been successfully applied to study the transient responses of beams, planar trusses, space frames and infinite layered media is extended to investigate the transient wave propagation and early short time transient response in finite Mindlin plate. Combining the wave solution, the shock source and the boundary conditions, the ray groups transmitted in the finite rectangular plate can be determined. Numerical simulations and experiments are performed and compared with each other. The results show that the transient wave propagation and early short time transient responses in the finite plate can be studied using the GRM. The early short time transient acceler- ations are very large for the finite plate subjected to the unit impulse, while the early short time transient displacements are very small. The early short time transient accelerations under the unit impulse are much larger than those under the unit step impulse. The thickness and material characteristics have remarkable effects on the early short time transient responses.