Although ray tracing produces high-fidelity, realistic images, it is considered computationally burdensome when implemented on a high rendering rate system. Perception-driven rendering techniques generate images with ...Although ray tracing produces high-fidelity, realistic images, it is considered computationally burdensome when implemented on a high rendering rate system. Perception-driven rendering techniques generate images with minimal noise and distortion that are generally acceptable to the human visual system, thereby reducing rendering costs. In this paper, we introduce a perception-entropy-driven temporal reusing method to accelerate real-time ray tracing. We first build a just noticeable difference(JND) model to represent the uncertainty of ray samples and image space masking effects. Then, we expand the shading gradient through gradient max-pooling and gradient filtering to enlarge the visual receipt field. Finally, we dynamically optimize reusable time segments to improve the accuracy of temporal reusing. Compared with Monte Carlo ray tracing, our algorithm enhances frames per second(fps) by 1.93× to 2.96× at 8 to 16 samples per pixel, significantly accelerating the Monte Carlo ray tracing process while maintaining visual quality.展开更多
To address the problem of subdividing inflexible rectangular grid models and their poor definition of velocity interfaces,we propose a complex structure triangular net for a minimum traveltime ray tracing global algor...To address the problem of subdividing inflexible rectangular grid models and their poor definition of velocity interfaces,we propose a complex structure triangular net for a minimum traveltime ray tracing global algorithm.Our procedure is：（1） Subdivide a triangle grid based on the Delaunay triangular subdivision criterion and the relationships of the points,lines,and the surfaces in the subdividing area.（2） Define the topology relationships and related concepts of triangular unit ray tracing.（3） The source point and wave arrival points at any time compose the propagating plane wave and the minimum traveltime and secondary source positions are calculated during the plane wave propagation.We adopt the hyperbolic approximation global algorithm for secondary source retrieving.（4） By minimum traveltime ray tracing,collect the path from receiver to source points with the neighborhood point＇s traveltime and the direction of the secondary source.Numerical simulation examples are given to test the algorithm.The results show that the triangular net ray tracing method demonstrates model subdivision flexibility,precise velocity discontinuity interfaces,and accurate computations.展开更多
For large-scale 3D seismic data,target-oriented reservoir imaging is more attractive than conventional full-volume migration,in terms of computation efficiency.Gaussian beam migration(GBM)is one of the most robust dep...For large-scale 3D seismic data,target-oriented reservoir imaging is more attractive than conventional full-volume migration,in terms of computation efficiency.Gaussian beam migration(GBM)is one of the most robust depth imaging method,which not only keeps the advantages of ray methods,such as high efficiency and flexibility,but also allows us to solve caustics and multipathing problems.But conventional Gaussian beam migration requires slant stack for prestack data,and ray tracing from beam center location to subsurface,which is not easy to be directly applied for target-oriented imaging.In this paper,we modify the conventional Gaussian beam migration scheme,by shooting rays from subsurface image points to receivers to implement wavefield back-propagation.This modification helps us to achieve a better subsurface illumination in complex structure and allows simple implementation for target reservoir imaging.Significantly,compared with the wavefi eld-based GBM,our method does not reconstruct the subsurface snapshots,which has higher efficiency.But the proposed method is not as efficient as the conventional Gaussian beam migration.Synthetic and field data examples demonstrate the validity and the target-oriented imaging capability of our method.展开更多
A real-time ray-based hardware emulator for Unmanned Aerial Vehicle(UAV) communication channels which suits for the Three-Dimensional(3D) dynamic scenarios and considers the movements of both terminals is developed in...A real-time ray-based hardware emulator for Unmanned Aerial Vehicle(UAV) communication channels which suits for the Three-Dimensional(3D) dynamic scenarios and considers the movements of both terminals is developed in this paper. The time-variant channel parameters, i.e.,ray delay, ray gain, and ray Doppler frequency are precalculated in the host by using the Ray Tracing(RT) method. Meanwhile, RT simulation dramatically increases the number of valid rays. To address the problem of resource limitation and huge computational burden in the implementation,an efficient ray coefficients generation method based on iteration is proposed and implemented.With the advantages of low cost and high flexibility, a Software Defined Radio(SDR) hardware platform is used to emulate the ray-based UAV channels by utilizing the compact architecture including the Time-Division(TD) scheme and Tapped-Delay Line(TDL) for channel convolution.Finally, hardware measurement results demonstrate that the properties of emulated channel, i.e.,Power Delay Profile(PDP) and Doppler Power Spectrum Density(DPSD) consist with the simulated ones, which verifies the correctness of hardware implementation. The proposed channel emulator provides an efficient way for optimization, verification, and evaluation of UAV communication systems.展开更多
Incorporation of the Monte Carlo(MC)algorithm in optimizing CyberKnife(CK)plans is cumbersome,and early models unconfigured MC calculations,therefore,this study investigated algorithm-based dose calculation discrepanc...Incorporation of the Monte Carlo(MC)algorithm in optimizing CyberKnife(CK)plans is cumbersome,and early models unconfigured MC calculations,therefore,this study investigated algorithm-based dose calculation discrepancies by selecting different prescription isodose lines(PIDLs)in head and lung CK plans.CK plans were based on anthropomorphic phantoms.Four shells were set at 2-60 mm from the target,and the constraint doses were adjusted according to the design stratcgy.After optimization,30%-90%PIDL plans were generated by ray tracing(RT).In the evaluation module,CK plans were recalculated using the MC algorithm.Therefore,the dosimetric parameters of different PIDL plans based on the RT and MC algorithms were obtained and analyzed.The discrepancies(mean+SD)were 3.72%+0.31%,3.40%+0.11%,3.47%+0.32%,0.17%+0.11%,0.64%+3.60%,7.73%+1.60%,14.62%+3.21%and 10.10%+1.57%for Djs,Dmeam),Dys,and coverage of the PTV,DGI,V,,V;and V,in the head plans and-6.32%+1.15%,-13.46%+0.98%,-20.63%+2.25%,-34.78%+25.03%,12248%+175.60%,-12.92%+5.41%,3.19%+4.67%and 7.13%+1.56%in the lung plans,respectively.The following parameters were significantly correlated with PIDL:dp98%at the 0.05 level and dpal,dys and dv3 at the 0.01 level for the head plans;dp98e%at the 0.05 level and do1e%,dpmeam,Ccoweange,dool,dvs and dv;at the 0.01 level for the lung plans.RT may be used to calculate the dose in CK head plans,but when the dose of organs at risk is close to the limit,it is necessary to refer to the MC results or to further optimize the CK plan to reduce the dose.For lung plans,the MC algorithm is recommended.For early models without the MC algorithm,a lower PIDL plan is recommended;otherwise,a large PIDL plan risks serious underdosage in the target area.展开更多
Based on Fermat’s principle, two-point ray tracing method was studied in three-dimensional structure. By means of first order Taylor’s incomplete series expansion (i.e. no expansion to the length of the ray), a symm...Based on Fermat’s principle, two-point ray tracing method was studied in three-dimensional structure. By means of first order Taylor’s incomplete series expansion (i.e. no expansion to the length of the ray), a symmetry block tridiagonal matrix equation set was deduced. Further, the positive definiteness of coefficient matrix was discussed, and the positive definiteness was accurately proved in a mathematical way. It assured that the algorithm was well-posed. Associated with iterative method, the solution to ray tracing can be got through step-by-step linearized iteration of the nonlinear problem. An algorithm of the whole path iterative ray tracing method in three-dimensional velocity structure was obtained. This method shows a clear and simple as well as explicit computation formula, which makes ray tracing computation easily applicable in practice. The correction vector is obtained through finding the solution to the positive definite block tridiagonal equation set, which ensures the method is robust convergence. This study offers a new kind of feasible and efficient ray tracing method for three dimensional seismic migration and tomography. Meanwhile, it also provides the prerequisite guarantee to design a fast algorithm.展开更多
A new algorithm for ray tracing bicubic Bezier surface intersection is presented. In order to find the intersected patches more effectively, a quadtree representation for surface patches is utilized. The introduction ...A new algorithm for ray tracing bicubic Bezier surface intersection is presented. In order to find the intersected patches more effectively, a quadtree representation for surface patches is utilized. The introduction of the alternative binary tree subdivision of Bezier surface effectively improves the performance of the ray surface intersection, with the benefit of avoiding the problem that the binary tree subdivision would produce long narrow patches and speeding the intersection finding process. The algorithm has been programmed in FORTRAN-77 and on prime-550 II computer. The result shows that its structure is simple and it is easy to implement with good performance.展开更多
This paper describes a technique to estimate surface-based duct parameters by using a simple ray tracing/correlation method. The approach is novel in that it incorporates the Spearman rank-order correlation scheme bet...This paper describes a technique to estimate surface-based duct parameters by using a simple ray tracing/correlation method. The approach is novel in that it incorporates the Spearman rank-order correlation scheme between the observed surface clutter and the surface ray density for a given propagation path. The simulation results and the real data results both demonstrate the ability of this method to estimate surface-based duct parameters. Compared with the results obtained by a modified genetic algorithm combined with the parabolic wave equation, the results retrieved from the ray tracing/correlation scheme show a minor reduction in accuracy but a great improvement on computation time. Therefore the ray tracing/correlation method might be used as a precursor to more sophisticated and slower techniques, such as genetic algorithm and particle filters, by narrowing the parameter search space and providing a comprehensive and more efficient estimation algorithm.展开更多
Seismic ray tracing in anisotropic media with irregular surface is crucial for the exploration of the fine crustal structure. Elliptically anisotropic medium is the type of anisotropic media with only four independent...Seismic ray tracing in anisotropic media with irregular surface is crucial for the exploration of the fine crustal structure. Elliptically anisotropic medium is the type of anisotropic media with only four independent elastic parameters. Usually, this medium can be described by only the vertical phase velocity and the horizontal phase velocity for seismic wave propagation. Model parameteri- zation in this study is described by flexible triangular grids, which is beneficial for the description of irregular surface with high degree of approximation. Both the vertical and horizontal phase velocities are defined in the triangular grids, respectively, which are used for the description of phase velocity distribution everywhere in the model by linear interpolation. We develop a shooting ray tracing method of turning wave in the elliptically anisotropic media with irregular surface. Runge-Kutta method is applied to solve the partial differential equation of seismic ray in elliptically anisotropic media. Linearly modified method is used for adjusting emergent phase angles in the shooting scheme. Numerical tests demonstrate that ray paths coincide well with analytical trajectories in trans- versely homogeneous elliptically anisotropic media. Seis- mic ray tracing results in transversely inhomogeneous elliptically anisotropic media demonstrate that our method is effective for further first-arrival tomography in ellipti- cally anisotropic media with an irregular surface.展开更多
Ultrasonic inspections on nickel-based alloy weldments meets some difficulties due to the curvilinear propagating paths will appear in both anisotropic and inhomogeneons welded joints. Thus, it is difficult to determi...Ultrasonic inspections on nickel-based alloy weldments meets some difficulties due to the curvilinear propagating paths will appear in both anisotropic and inhomogeneons welded joints. Thus, it is difficult to determine the exact location of flaws using the traditional ultrasonic testing method which is based on the assumption that the sound beam will propagate in straight lines. In order to overcome this problem, we provide a new model-based inspection approach to locate the flaws in the nickel-based alloy weldments. Furthermore, some experimental examinations are carried out to compare the locating accuracy between the traditional and the modeled-based approach. It shows that the provided model-based inspection method is more accurate than the traditional method in inspection of the nickel-based alloy weldments.展开更多
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 tracin...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.展开更多
The Doppler backscattering system has been widely used for turbulence measurements,and the microwave beam will be backscattered near the cut-off layer when the Brag condition is fulfilled.In tokamak,the ray-tracing co...The Doppler backscattering system has been widely used for turbulence measurements,and the microwave beam will be backscattered near the cut-off layer when the Brag condition is fulfilled.In tokamak,the ray-tracing code is used to obtain the radial position and perpendicular wave number of the scattering layer for turbulence velocity measurement and the WKB（Wentzel-Kramers-Brillouin） approximation should be satisfied for optical propagation.To calculate the backscattering location and wave number at the cut-off layer only,a single ray tracing in the cross section is enough,while for spatial and wave number resolution calculation,multiple rays reflecting the microwave beam size should be used.Considering the angle between the wave vector and the magnetic field,a three-dimension quasi-optical Gaussian ray tracing is sometimes needed.展开更多
This paper presents a gradient-descent travel time tomography method for solving the acoustictype velocity model inversion problem.Similarly to the adjoint-state method,the proposed method is based on the Eikonal equa...This paper presents a gradient-descent travel time tomography method for solving the acoustictype velocity model inversion problem.Similarly to the adjoint-state method,the proposed method is based on the Eikonal equation,enabling simultaneous calculation of contributions from all common-source receivers to the gradient.This overcomes the inefficiency inherent in conventional travel time tomography methods,which rely on a two-point ray tracing process.By directly calculating Fréchet derivatives,our method avoids the complex derivation processes associated with the adjoint-state method.The key to calculating the Fréchet derivatives is to calculate a so-called ray-path term.Consequently,compared to the adjoint-state method,the proposed method can explicitly obtain the ray paths,resulting in a more concise and intuitive derivation process.Furthermore,our method retains the benefits of the adjoint-state method,such as speed,low memory usage,and robustness.This paper focuses on elucidating the principles and algorithms for calculating the raypath term based on the fast sweeping method.The algorithms could be further speeded up by using parallel computational techniques.Synthetic tests demonstrate that our proposed travel time tomographic method accurately calculates ray paths,regardless of the complexity of the model and recording geometry.展开更多
Kinetic simulation is a powerful tool to study the excitation and propagation of whistler-mode waves in the Earth’s inner magnetosphere.This method typically applies a scaled-down dipole magnetic field to save comput...Kinetic simulation is a powerful tool to study the excitation and propagation of whistler-mode waves in the Earth’s inner magnetosphere.This method typically applies a scaled-down dipole magnetic field to save computational time.However,it remains unknown whether whistler wave propagation in the scaled-down dipole field is consistent with that in the realistic dipole field.In this work,we develop a ray-tracing code with a scalable dipole magnetic field to address this concern.The simulation results show that parallel whistler waves at different frequencies gradually become oblique after leaving the equator and propagate in different raypaths in a dipole magnetic field.During their propagation,the higher frequency waves tend to have larger wave normal angles at the same latitude.Compared with the wave propagation in a realistic dipole field,the wave raypath and wave normal remain the same,whereas the wave amplification or attenuation is smaller because of the shorter propagation time in a scaled-down dipole field.Our study provides significant guidance for kinetic simulations of whistler-mode waves.展开更多
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.展开更多
Power deposition profiles generated by Ion Cyclotron Resonance Heating （ICRH） in non-circular tokamaks are studied using a ray-tracing technique. The simulation results for the Experimental Advanced Superconductor T...Power deposition profiles generated by Ion Cyclotron Resonance Heating （ICRH） in non-circular tokamaks are studied using a ray-tracing technique. The simulation results for the Experimental Advanced Superconductor Tokamak （EAST） D-shaped plasma are presented. It is indicated that the spatial distributions of plasma parameters （plasma density, species temperature, minority ion concentration, etc.） have an significant influence on the power deposition profiles. The findings may be highly useful to the planned plasma heating and experiments in EAST.展开更多
A fast algorithm for ray tracing is presented, with which the specular reflection term of global illumination model is improved. A hybrid technique combining hierarchical bounding volumes and constant size box partiti...A fast algorithm for ray tracing is presented, with which the specular reflection term of global illumination model is improved. A hybrid technique combining hierarchical bounding volumes and constant size box partitioning is presented and a fast box traversal algorithm is used. By this technique multiple ray intersections with objects that are in more than one box can be avoided. As a result, the speed of ray tracing is considerably increased.展开更多
Ray tracing study of electromagnetic ion cyclotron （EMIC） waves is conducted based on a realistic plasma density model. The simulation result shows that EMIC waves propagate away from the equatorial source region to...Ray tracing study of electromagnetic ion cyclotron （EMIC） waves is conducted based on a realistic plasma density model. The simulation result shows that EMIC waves propagate away from the equatorial source region to higher latitudes basically along geomagnetic field lines, and are reflected at the region where their frequency matches the local bi-ion frequency. H＋ band suffers H＋-He＋ bi-ion frequency reflection at lower latitudes, whereas He＋ band suffers He＋-O＋ bi-ion frequency reflection at higher latitudes. Moreover, the concentration of heavy ions slightly affects the bi-ion frequencies and then slightly determines the reflection location of ray paths of EMIC waves. The current results present the first detailed study on the propagation characteristics of EMIC waves associated with bi-ion frequencies.展开更多
A three-dimensional ray tracing study of a whistler-mode chorus is conducted for different geomagnetic activities by using a global core plasma density model. For the upperband chorus, the initial azimuthal wave angle...A three-dimensional ray tracing study of a whistler-mode chorus is conducted for different geomagnetic activities by using a global core plasma density model. For the upperband chorus, the initial azimuthal wave angle affects slightly the projection of ray trajectories onto the plane （Z, √（x^2 ＋ y^2））, but controls the longitudinal propagation. The trajectory of the upper-band chorus is strongly associated with the plasmapause and the magnetic local time （MLT） of chorus source region. For the high geomagnetic activity, the chorus trajectory moves inward together with the plasmapause. In the bulge region, the plasmapause extends outward, while the chorus trajectory moves outward together with the plasmapause. For moderately or high geomagnetic activity, the lower-band chorus suffers low hybrid resonance （LHR） reflection before it reaches the plasmapause, leading to a weak correlation with the geomagnetic activity and magnetic local time of the chorus source region. For low geomagnetic activity, the lower-band chorus may be reflected firstly at the plasmapause instead of suffering LHR reflection, exhibiting a propagation characteristic similar to that of the upper-band chorus. The results provide a new insight into the propagation characteristics of the chorus for different geomagnetic activities and contribute to further understanding of the acceleration of energetic electron by a chorus wave.展开更多
An advanced method for rapidly computing holograms of large three-dimensional(3D)objects combines backward ray tracing with adaptive resolution wavefront recording plane(WRP)and adaptive angular spectrum propagation.I...An advanced method for rapidly computing holograms of large three-dimensional(3D)objects combines backward ray tracing with adaptive resolution wavefront recording plane(WRP)and adaptive angular spectrum propagation.In the initial phase,a WRP with adjustable resolution and sampling interval based on the object’s size is defined to capture detailed information from large 3D objects.The second phase employs an adaptive angular spectrum method(ASM)to efficiently compute the propagation from the large-sized WRP to the small-sized computer-generated hologram(CGH).The computation process is accelerated using CUDA and OptiX.Optical experiments confirm that the algorithm can generate high-quality holograms with shadow and occlusion effects at a resolution of 1024×1024 in 29 ms.展开更多
基金supported by the National Natural Science Foundation of China (No.U19A2063)the Jilin Provincial Science&Technology Development Program of China (No.20230201080GX)。
文摘Although ray tracing produces high-fidelity, realistic images, it is considered computationally burdensome when implemented on a high rendering rate system. Perception-driven rendering techniques generate images with minimal noise and distortion that are generally acceptable to the human visual system, thereby reducing rendering costs. In this paper, we introduce a perception-entropy-driven temporal reusing method to accelerate real-time ray tracing. We first build a just noticeable difference(JND) model to represent the uncertainty of ray samples and image space masking effects. Then, we expand the shading gradient through gradient max-pooling and gradient filtering to enlarge the visual receipt field. Finally, we dynamically optimize reusable time segments to improve the accuracy of temporal reusing. Compared with Monte Carlo ray tracing, our algorithm enhances frames per second(fps) by 1.93× to 2.96× at 8 to 16 samples per pixel, significantly accelerating the Monte Carlo ray tracing process while maintaining visual quality.
基金the National Natural Science Foundation of China(Grant Nos.50974081,50774051, 51034003)the Foundation for the Author of National Excellent Doctoral Dissertation of PR China(Grant No.200958)+1 种基金the Program for Changjiang Scholars and Innovative Research Team in University(Grant No.IRT0843)the Open Research Fund Program of the State Key Laboratory of Coal Resources and Safe Mining(China University of Mining and Technology, Beijing)(No.2007-07)
文摘To address the problem of subdividing inflexible rectangular grid models and their poor definition of velocity interfaces,we propose a complex structure triangular net for a minimum traveltime ray tracing global algorithm.Our procedure is：（1） Subdivide a triangle grid based on the Delaunay triangular subdivision criterion and the relationships of the points,lines,and the surfaces in the subdividing area.（2） Define the topology relationships and related concepts of triangular unit ray tracing.（3） The source point and wave arrival points at any time compose the propagating plane wave and the minimum traveltime and secondary source positions are calculated during the plane wave propagation.We adopt the hyperbolic approximation global algorithm for secondary source retrieving.（4） By minimum traveltime ray tracing,collect the path from receiver to source points with the neighborhood point＇s traveltime and the direction of the secondary source.Numerical simulation examples are given to test the algorithm.The results show that the triangular net ray tracing method demonstrates model subdivision flexibility,precise velocity discontinuity interfaces,and accurate computations.
文摘For large-scale 3D seismic data,target-oriented reservoir imaging is more attractive than conventional full-volume migration,in terms of computation efficiency.Gaussian beam migration(GBM)is one of the most robust depth imaging method,which not only keeps the advantages of ray methods,such as high efficiency and flexibility,but also allows us to solve caustics and multipathing problems.But conventional Gaussian beam migration requires slant stack for prestack data,and ray tracing from beam center location to subsurface,which is not easy to be directly applied for target-oriented imaging.In this paper,we modify the conventional Gaussian beam migration scheme,by shooting rays from subsurface image points to receivers to implement wavefield back-propagation.This modification helps us to achieve a better subsurface illumination in complex structure and allows simple implementation for target reservoir imaging.Significantly,compared with the wavefi eld-based GBM,our method does not reconstruct the subsurface snapshots,which has higher efficiency.But the proposed method is not as efficient as the conventional Gaussian beam migration.Synthetic and field data examples demonstrate the validity and the target-oriented imaging capability of our method.
基金supported in part by CEMEE State Key Laboratory Fund(No.2020Z0207B)in part by the National Key Scientific Instrument and Equipment Development Project,China(No.61827801)+1 种基金in part by the State Key Laboratory of Integrated Services Networks Fund,China(No.ISN22-11)in part by the Fundamental Research Funds for the Central Universities,China(Nos.NS2020026 and NS2020063)。
文摘A real-time ray-based hardware emulator for Unmanned Aerial Vehicle(UAV) communication channels which suits for the Three-Dimensional(3D) dynamic scenarios and considers the movements of both terminals is developed in this paper. The time-variant channel parameters, i.e.,ray delay, ray gain, and ray Doppler frequency are precalculated in the host by using the Ray Tracing(RT) method. Meanwhile, RT simulation dramatically increases the number of valid rays. To address the problem of resource limitation and huge computational burden in the implementation,an efficient ray coefficients generation method based on iteration is proposed and implemented.With the advantages of low cost and high flexibility, a Software Defined Radio(SDR) hardware platform is used to emulate the ray-based UAV channels by utilizing the compact architecture including the Time-Division(TD) scheme and Tapped-Delay Line(TDL) for channel convolution.Finally, hardware measurement results demonstrate that the properties of emulated channel, i.e.,Power Delay Profile(PDP) and Doppler Power Spectrum Density(DPSD) consist with the simulated ones, which verifies the correctness of hardware implementation. The proposed channel emulator provides an efficient way for optimization, verification, and evaluation of UAV communication systems.
基金This study was supported by grants from National Key Research and Development Plan for Digital Diagnostic Equipment Research and Development(No.2016YFC0106700)the Natural Science Foundation of Union Hospital,Tongji Medical College,Huazhong University of Science and Technology(No.02.03.2018-131).
文摘Incorporation of the Monte Carlo(MC)algorithm in optimizing CyberKnife(CK)plans is cumbersome,and early models unconfigured MC calculations,therefore,this study investigated algorithm-based dose calculation discrepancies by selecting different prescription isodose lines(PIDLs)in head and lung CK plans.CK plans were based on anthropomorphic phantoms.Four shells were set at 2-60 mm from the target,and the constraint doses were adjusted according to the design stratcgy.After optimization,30%-90%PIDL plans were generated by ray tracing(RT).In the evaluation module,CK plans were recalculated using the MC algorithm.Therefore,the dosimetric parameters of different PIDL plans based on the RT and MC algorithms were obtained and analyzed.The discrepancies(mean+SD)were 3.72%+0.31%,3.40%+0.11%,3.47%+0.32%,0.17%+0.11%,0.64%+3.60%,7.73%+1.60%,14.62%+3.21%and 10.10%+1.57%for Djs,Dmeam),Dys,and coverage of the PTV,DGI,V,,V;and V,in the head plans and-6.32%+1.15%,-13.46%+0.98%,-20.63%+2.25%,-34.78%+25.03%,12248%+175.60%,-12.92%+5.41%,3.19%+4.67%and 7.13%+1.56%in the lung plans,respectively.The following parameters were significantly correlated with PIDL:dp98%at the 0.05 level and dpal,dys and dv3 at the 0.01 level for the head plans;dp98e%at the 0.05 level and do1e%,dpmeam,Ccoweange,dool,dvs and dv;at the 0.01 level for the lung plans.RT may be used to calculate the dose in CK head plans,but when the dose of organs at risk is close to the limit,it is necessary to refer to the MC results or to further optimize the CK plan to reduce the dose.For lung plans,the MC algorithm is recommended.For early models without the MC algorithm,a lower PIDL plan is recommended;otherwise,a large PIDL plan risks serious underdosage in the target area.
基金Project(40674071) supported by the National Natural Science Foundation of Chinaproject(KFAS2002-2003) supported by Korea Foundation for Advanced Studies
文摘Based on Fermat’s principle, two-point ray tracing method was studied in three-dimensional structure. By means of first order Taylor’s incomplete series expansion (i.e. no expansion to the length of the ray), a symmetry block tridiagonal matrix equation set was deduced. Further, the positive definiteness of coefficient matrix was discussed, and the positive definiteness was accurately proved in a mathematical way. It assured that the algorithm was well-posed. Associated with iterative method, the solution to ray tracing can be got through step-by-step linearized iteration of the nonlinear problem. An algorithm of the whole path iterative ray tracing method in three-dimensional velocity structure was obtained. This method shows a clear and simple as well as explicit computation formula, which makes ray tracing computation easily applicable in practice. The correction vector is obtained through finding the solution to the positive definite block tridiagonal equation set, which ensures the method is robust convergence. This study offers a new kind of feasible and efficient ray tracing method for three dimensional seismic migration and tomography. Meanwhile, it also provides the prerequisite guarantee to design a fast algorithm.
文摘A new algorithm for ray tracing bicubic Bezier surface intersection is presented. In order to find the intersected patches more effectively, a quadtree representation for surface patches is utilized. The introduction of the alternative binary tree subdivision of Bezier surface effectively improves the performance of the ray surface intersection, with the benefit of avoiding the problem that the binary tree subdivision would produce long narrow patches and speeding the intersection finding process. The algorithm has been programmed in FORTRAN-77 and on prime-550 II computer. The result shows that its structure is simple and it is easy to implement with good performance.
基金Project supported by the National Natural Science Foundation of China (Grant No. 40775025)
文摘This paper describes a technique to estimate surface-based duct parameters by using a simple ray tracing/correlation method. The approach is novel in that it incorporates the Spearman rank-order correlation scheme between the observed surface clutter and the surface ray density for a given propagation path. The simulation results and the real data results both demonstrate the ability of this method to estimate surface-based duct parameters. Compared with the results obtained by a modified genetic algorithm combined with the parabolic wave equation, the results retrieved from the ray tracing/correlation scheme show a minor reduction in accuracy but a great improvement on computation time. Therefore the ray tracing/correlation method might be used as a precursor to more sophisticated and slower techniques, such as genetic algorithm and particle filters, by narrowing the parameter search space and providing a comprehensive and more efficient estimation algorithm.
基金financial support for this work contributed by the National Key Research and Development Program of China(Grants Nos.2016YFC0600101,2016YFC0600201 and 2016YFC0600302)the National Natural Science Foundation of China(Grants Nos.41522401 and 41474068)
文摘Seismic ray tracing in anisotropic media with irregular surface is crucial for the exploration of the fine crustal structure. Elliptically anisotropic medium is the type of anisotropic media with only four independent elastic parameters. Usually, this medium can be described by only the vertical phase velocity and the horizontal phase velocity for seismic wave propagation. Model parameteri- zation in this study is described by flexible triangular grids, which is beneficial for the description of irregular surface with high degree of approximation. Both the vertical and horizontal phase velocities are defined in the triangular grids, respectively, which are used for the description of phase velocity distribution everywhere in the model by linear interpolation. We develop a shooting ray tracing method of turning wave in the elliptically anisotropic media with irregular surface. Runge-Kutta method is applied to solve the partial differential equation of seismic ray in elliptically anisotropic media. Linearly modified method is used for adjusting emergent phase angles in the shooting scheme. Numerical tests demonstrate that ray paths coincide well with analytical trajectories in trans- versely homogeneous elliptically anisotropic media. Seis- mic ray tracing results in transversely inhomogeneous elliptically anisotropic media demonstrate that our method is effective for further first-arrival tomography in ellipti- cally anisotropic media with an irregular surface.
基金This study was supported by State Key Laboratory of Advanced Welding Production Technology and National Natural Science Foundation of China （Grant No. 50975028, 50775054）.
文摘Ultrasonic inspections on nickel-based alloy weldments meets some difficulties due to the curvilinear propagating paths will appear in both anisotropic and inhomogeneons welded joints. Thus, it is difficult to determine the exact location of flaws using the traditional ultrasonic testing method which is based on the assumption that the sound beam will propagate in straight lines. In order to overcome this problem, we provide a new model-based inspection approach to locate the flaws in the nickel-based alloy weldments. Furthermore, some experimental examinations are carried out to compare the locating accuracy between the traditional and the modeled-based approach. It shows that the provided model-based inspection method is more accurate than the traditional method in inspection of the nickel-based alloy weldments.
文摘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.
基金supported by National Natural Science Foundation of China(Nos.10990211 and 11105146)the ITER-CN Project973 Program of China(No.2013GB106002)
文摘The Doppler backscattering system has been widely used for turbulence measurements,and the microwave beam will be backscattered near the cut-off layer when the Brag condition is fulfilled.In tokamak,the ray-tracing code is used to obtain the radial position and perpendicular wave number of the scattering layer for turbulence velocity measurement and the WKB（Wentzel-Kramers-Brillouin） approximation should be satisfied for optical propagation.To calculate the backscattering location and wave number at the cut-off layer only,a single ray tracing in the cross section is enough,while for spatial and wave number resolution calculation,multiple rays reflecting the microwave beam size should be used.Considering the angle between the wave vector and the magnetic field,a three-dimension quasi-optical Gaussian ray tracing is sometimes needed.
基金supported by 14th Five-Year Plan major science and technology projects(no.KJGG2022-0201)。
文摘This paper presents a gradient-descent travel time tomography method for solving the acoustictype velocity model inversion problem.Similarly to the adjoint-state method,the proposed method is based on the Eikonal equation,enabling simultaneous calculation of contributions from all common-source receivers to the gradient.This overcomes the inefficiency inherent in conventional travel time tomography methods,which rely on a two-point ray tracing process.By directly calculating Fréchet derivatives,our method avoids the complex derivation processes associated with the adjoint-state method.The key to calculating the Fréchet derivatives is to calculate a so-called ray-path term.Consequently,compared to the adjoint-state method,the proposed method can explicitly obtain the ray paths,resulting in a more concise and intuitive derivation process.Furthermore,our method retains the benefits of the adjoint-state method,such as speed,low memory usage,and robustness.This paper focuses on elucidating the principles and algorithms for calculating the raypath term based on the fast sweeping method.The algorithms could be further speeded up by using parallel computational techniques.Synthetic tests demonstrate that our proposed travel time tomographic method accurately calculates ray paths,regardless of the complexity of the model and recording geometry.
基金supported by the National Natural Science Foundation of China (Grant No. 42104155)the China Postdoctoral Science Foundation (Grant No. 2021M693049)+1 种基金the Fundamental Research Funds for the Central Universities (Grant Nos. WK2080000150 and WK3420000013)the USTC (University of Science and Technology of China) Tang Scholar Program
文摘Kinetic simulation is a powerful tool to study the excitation and propagation of whistler-mode waves in the Earth’s inner magnetosphere.This method typically applies a scaled-down dipole magnetic field to save computational time.However,it remains unknown whether whistler wave propagation in the scaled-down dipole field is consistent with that in the realistic dipole field.In this work,we develop a ray-tracing code with a scalable dipole magnetic field to address this concern.The simulation results show that parallel whistler waves at different frequencies gradually become oblique after leaving the equator and propagate in different raypaths in a dipole magnetic field.During their propagation,the higher frequency waves tend to have larger wave normal angles at the same latitude.Compared with the wave propagation in a realistic dipole field,the wave raypath and wave normal remain the same,whereas the wave amplification or attenuation is smaller because of the shorter propagation time in a scaled-down dipole field.Our study provides significant guidance for kinetic simulations of whistler-mode waves.
基金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.
文摘Power deposition profiles generated by Ion Cyclotron Resonance Heating （ICRH） in non-circular tokamaks are studied using a ray-tracing technique. The simulation results for the Experimental Advanced Superconductor Tokamak （EAST） D-shaped plasma are presented. It is indicated that the spatial distributions of plasma parameters （plasma density, species temperature, minority ion concentration, etc.） have an significant influence on the power deposition profiles. The findings may be highly useful to the planned plasma heating and experiments in EAST.
文摘A fast algorithm for ray tracing is presented, with which the specular reflection term of global illumination model is improved. A hybrid technique combining hierarchical bounding volumes and constant size box partitioning is presented and a fast box traversal algorithm is used. By this technique multiple ray intersections with objects that are in more than one box can be avoided. As a result, the speed of ray tracing is considerably increased.
基金supported by National Natural Science Foundation of China(Nos.40925014,41204114,and 41274165)the Aid Program for Scienceand Technology Innovative Research Team in Higher Educational Institutions of Hunan Province,Chinathe Construct Program of theKey Discipline in Hunan Province,China
文摘Ray tracing study of electromagnetic ion cyclotron （EMIC） waves is conducted based on a realistic plasma density model. The simulation result shows that EMIC waves propagate away from the equatorial source region to higher latitudes basically along geomagnetic field lines, and are reflected at the region where their frequency matches the local bi-ion frequency. H＋ band suffers H＋-He＋ bi-ion frequency reflection at lower latitudes, whereas He＋ band suffers He＋-O＋ bi-ion frequency reflection at higher latitudes. Moreover, the concentration of heavy ions slightly affects the bi-ion frequencies and then slightly determines the reflection location of ray paths of EMIC waves. The current results present the first detailed study on the propagation characteristics of EMIC waves associated with bi-ion frequencies.
基金supported by National Natural Science Foundation of China(Nos.40874076,40925014)
文摘A three-dimensional ray tracing study of a whistler-mode chorus is conducted for different geomagnetic activities by using a global core plasma density model. For the upperband chorus, the initial azimuthal wave angle affects slightly the projection of ray trajectories onto the plane （Z, √（x^2 ＋ y^2））, but controls the longitudinal propagation. The trajectory of the upper-band chorus is strongly associated with the plasmapause and the magnetic local time （MLT） of chorus source region. For the high geomagnetic activity, the chorus trajectory moves inward together with the plasmapause. In the bulge region, the plasmapause extends outward, while the chorus trajectory moves outward together with the plasmapause. For moderately or high geomagnetic activity, the lower-band chorus suffers low hybrid resonance （LHR） reflection before it reaches the plasmapause, leading to a weak correlation with the geomagnetic activity and magnetic local time of the chorus source region. For low geomagnetic activity, the lower-band chorus may be reflected firstly at the plasmapause instead of suffering LHR reflection, exhibiting a propagation characteristic similar to that of the upper-band chorus. The results provide a new insight into the propagation characteristics of the chorus for different geomagnetic activities and contribute to further understanding of the acceleration of energetic electron by a chorus wave.
基金Project supported by the Special Project of Central Government Guiding Local Science and Technology Development in Beijing 2020(Grant No.Z201100004320006).
文摘An advanced method for rapidly computing holograms of large three-dimensional(3D)objects combines backward ray tracing with adaptive resolution wavefront recording plane(WRP)and adaptive angular spectrum propagation.In the initial phase,a WRP with adjustable resolution and sampling interval based on the object’s size is defined to capture detailed information from large 3D objects.The second phase employs an adaptive angular spectrum method(ASM)to efficiently compute the propagation from the large-sized WRP to the small-sized computer-generated hologram(CGH).The computation process is accelerated using CUDA and OptiX.Optical experiments confirm that the algorithm can generate high-quality holograms with shadow and occlusion effects at a resolution of 1024×1024 in 29 ms.
