Ray Tracing Algorithm Research Articles

A high-performance ray tracing particle tracking model for the simulation of microplastics in inland and coastal aquatic environments

In this study, we present a high-performance Particle Tracking Model (PTM) designed for simulating any type of particles, with a focus on microplastics. The PTM is efficient compared to existing models, parallelized, and utilizes a ray tracing algorithm incorporating both ray reflection and ray refraction in order to traverse particles as well as find the location of each particle over three-dimensional unstructured grids. Various numerical corrections are implemented in the model to address computational round-off errors and discontinuities in the water surface level of the input hydrodynamic models. To increase the accuracy of the model, partially reflective boundary conditions are imposed as well as the capability to simulate microplastics beaching and washout in very shallow areas or dry computational cells. Several tests are conducted to study the performance, scalability, and accuracy of the model. The proposed model is tested with over 3.88 billion double-precision particles on three-dimensional computational grids with up to approximately one million cells. The tests show that the ray tracing approach is efficient, achieves over 17× faster runtime, and offers greater accuracy compared to using an auxiliary grid for particle location finding. For larger timesteps, the ray tracing PTM with refraction shows improved accuracy compared to the ray tracing PTM without refraction. The model's capabilities are tested in a real-world case study over the Saguenay Fjord, Quebec, Canada. The model is utilized to reproduce the paths of five surface drifters. A second numerical test is conducted in the Fjord and high particle concentration areas are identified.

A New Framework for Understanding Systematic Errors in Cluster Lens Modeling. III. Deflection from Large-Scale Structure

Interpreting and reconstructing distant sources that are gravitationally lensed by galaxy clusters requires accurate and precise lens models. While high-quality data sets have reduced statistical errors in such models, systematic errors remain important. We examine systematic lensing effects caused by density fluctuations due to large-scale structure along the line of sight. We use a multiplane ray-tracing algorithm with the IllustrisTNG 100-3 cosmological simulation of matter distribution and compute the statistical distributions of shear, convergence, and higher-order deflections using two Hubble Frontier Field clusters as examples (A2744 and MACS J0416.1−2403). The cosmic shear distribution is Gaussian in each component, while the cosmic convergence distribution is skewed such that 1 + κ is consistent with a log-normal distribution; the standard deviations for these quantities are at the level of a few to 10%, depending on the redshift of the source. The deflection from higher-order terms beyond convergence and shear has significant scatter: the rms deflection is ∼15″, considerably larger than the image position residuals for current lens models. These results indicate that line-of-sight deflection effects due to large-scale structure can significantly impact lens models and should not be neglected. We present results in forms that can be incorporated into future cluster lens models.

Research on wireless channel model based on improved ray tracing algorithm that considers multiple diffuse scattering and employs pyramid-shaped ray tubes as the carrier

This paper extensively utilizes fine three dimensional environmental data obtained from laser point clouds. Based on theories such as geometrical optics and effective roughness theory, a deterministic wireless channel model is established, which integrates higher-order diffuse scattering. This model is referred to as the ray tracing fusion with higher-order diffuse scattering model. To expedite the collision calculation between rays and the scene, this paper introduces a combined approach of voxelization and signed distance field, resulting in a remarkable 16-fold improvement in computational speed. Moreover, aiming to balance accuracy and efficiency, the paper systematically analyzes the optimization computation parameters of the model. Finally, the proposed model is validated using measurement data in the frequency range of 1 GHz to 6 GHz in mountainous terrain. The results indicate that the predicted outcomes of the proposed model have an accuracy within 6 dB compared to the measurement results, and are superior to ITU-R P.1546, which is an international standard recommended by the International Telecommunication Union for modeling electromagnetic wave propagation in undulating terrain. This provides necessary technical support for network planning and optimization.

Secondary Gravity Wave Propagation in Tropical Thermospheric Region: Role of Varying Kinematic Viscosity

AbstractThe current study has investigated the thunderstorm induced atmospheric gravity waves (AGWs) over Indian region based on the perturbation signatures in ionospheric total electron content (TEC) measurement. Robust traveling ionospheric disturbance (TID) signature has been identified along the east side of the thunderstorm affected area. Neutral wind was found to have a favorable impact in this aspect for a certain time duration of the day by modulating the vertical wavelength. The role of temperature was analyzed in terms of kinematic viscosity which is a crucial component, especially over tropical region, for wave dissipation and reflection along its propagation path. Ray tracing algorithm is also applied with varying kinematic viscosity and thermal diffusivity for retrieval of possible ray paths and source location of observed waves. A statistical investigation has been carried out to identify the dissipation altitude of observed waves along the ray paths. It has been found that all waves dissipated at almost a constant altitude for a specific kinematic viscosity and above this altitude vertical wavelength was found to decrease. The ray paths interacted at a common point which was located at about 125 km altitude and was very close to the region of maximum lightning activity. It can also be noted that the observed phase velocities can't be achieved by a wave below the turbopause. It indicates that the observed waves were excited from a secondary source and not directly connected to convective system. The study provides an in‐depth analysis of mesoscale system induced gravity wave propagation and dissipation over tropical region.

Simulating the Universe from the cosmological horizon to halo scales

Ultra-large scales close to the cosmological horizon will be probed by the upcoming observational campaigns. They hold the promise to constrain single-field inflation as well as general relativity, but in order to include them in the forthcoming analyses, their modelling has to be robust. In particular, general relativistic effects may be mistaken for primordial signals, and no consensus has emerged either from analytical modelling nor from the numerical route, obstructed by the large dynamical range to be simulated. In this work, we present a numerical technique to overcome the latter limitation: we compute the general relativistic displacement field with the N-body relativistic code gevolution and combine it with the accurate Newtonian simulation Gadget-4. This combination leads to an effective simulation reproducing the desired behaviour at the level of the matter power spectrum and bispectrum. We then measure, for the first time in a simulation, the relativistic scale-dependent bias in Poisson gauge; at redshift z = 0, we find b 1 GR = -8.1 ± 2.8. Our results at the field level are only valid in the Poisson gauge and need to be complemented with a relativistic ray tracing algorithm to compute the number count observable.

Development of an algorithm for proton dose calculation in magnetic fields.

The advantages of proton therapy can be further enhanced with online magnetic resonance imaging (MRI) guidance. One of the challenges in the realization of MRI-guided proton therapy (MRPT) is accurately calculating the radiation dose in the presence of magnetic fields. This study aims to develop an efficient and accurate proton dose calculation algorithm adapted to the presence of magnetic fields. An analytical-numerical radiation dose calculation algorithm, Proton and Ion Dose Engine (PRIDE), was developed. The algorithm combines the pencil beam algorithm (PBA) with a novel iterative voxel-based ray-tracing algorithm. The new ray-tracing method uses fewer assumptions and ensures broader applicability for proton beam trajectory prediction in magnetic fields, and has been compared to Wolf's method and Schellhammer's method. The accuracy of PRIDE algorithm was validated on three phantoms and two practical plans (one single-field water plan and one prostate tumor plan) in different magnetic field strengths up to 3.0 T. The validation was performed by comparing the results against the Monte Carlo (MC) simulations, using the global gamma index criteria of 2%/2mm and 3%/3mm with a 10% threshold. PRIDE showed good agreement with MC in homogeneous and slab heterogeneous phantom, achieving gamma passing rates (%GPs) above 99% for 2%/2mm criteria when magnetic field strength is not greater than 1.5 T. Although the agreement decreased for scenarios involving high proton energy (240 MeV) and strong magnetic field (3.0 T), the 2%/2mm %GPs still remained above 98%. In lateral heterogeneous phantom, the accuracy of PRIDE decreased due to the PBA's limitation. For the two practical plans in different magnetic fields, %GPs exceeded 98% and 99% for 2%/2mm and 3%/3mm criteria, respectively. PRIDE can perform efficient and accurate proton dose calculation in magnetic fields up to 3.0 T, and is expected to work as a useful tool for proton dose calculation in MRPT.

A GPU based 3D raytracing algorithm for DUED laser fusion code

Abstract These days, graphical processing units (GPUs) deliver performance
comparable to that of hundreds of CPU cores. This level of performance allows
certain classes of simulations to be run in-house on a standard consumer workstation,
eliminating the need for a cluster. In this paper, it is shown that medium-resolution,
2D radiation hydrodynamics simulations for laser-driven inertial confinement fusion
with realistic 3D laser raytracing can now be conducted on a single consumer device.
A novel raytracing module has indeed been developed for the 2D Lagrangian radiation-
hydro-nuclear code DUED to leverage the computational power of GPUs. By
employing 3D raytracing, more realistic investigations of laser-driven plasmas become
feasible, with a particular focus on perturbations resulting from non-uniform laser
irradiation.

Alignment of a 360° image with posed color images for locally accurate texturing of 3D mesh

With the popularity of 3D content like virtual tours, the challenges of 3D data registration have become increasingly significant. The registration of heterogeneous data obtained from 2D and 3D sensors is required to create photo-realistic 3D models. However, the alignment of 2D images with 3D models introduces a significant challenge due to their inherent differences. This article introduces a rigorous mathematical approach to align a 360° image with its corresponding 3D model generated from images with known camera poses followed by texture projection on the model. We use Scale-Invariant Feature Transform (SIFT) feature descriptors enhanced with a homography-based metric to establish correspondences between the faces of a cubemap and the posed images. To achieve optimal alignment, we use a non-linear least squares optimization technique with a custom objective function. Subsequently, the outcomes of the alignment process are evaluated through texturing using a customized raytracing algorithm. The resulting projections are compared against the original textures, with a comprehensive assessment of the alignment's fidelity and precision.

Ray-Tracing-Assisted SAR Image Simulation under Range Doppler Imaging Geometry

In order to achieve an effective balance between SAR image simulation fidelity and efficiency, we proposed a ray-tracing-assisted SAR image simulation method under range doppler (RD) imaging geometry. This method utilizes the spatial traversal mode of RD imaging geometry to transmit discrete electromagnetic (EM) waves into the SAR radiation area and follows the Nyquist sampling law to set the density of transmitted EM waves to effectively identify the beam radiation area. The ray-tracing algorithm is used to obtain the backscatter amplitude and real-time slant range of the transmitted EM wave, which can effectively record the multiple backscattering among the components of the distributed target so that the backscattering subfields of each component can be correlated. According to the RD condition equation, the backscattering amplitude is assigned to the corresponding range gate, and the three-dimensional (3D) target is mapped into the two-dimensional (2D) SAR slant-range coordinate system, and the SAR target simulated image is directly obtained. Finally, the simulation images of the proposed method are compared qualitatively and quantitatively with those obtained by commercial simulation software, and the effectiveness of the proposed method is verified.

Ray tracing routing using packet reception timing in dense nanonetworks

This paper presents a novel routing method called ray tracing. The method targets wireless nanonetworks, which are networks of nanometer-sized nodes with applications in various domains such as future medicine and metamaterial. Due to their tiny size, nanonodes have quite limited resources, in particular energy, and it is crucial to use as few resources as possible. In multi-hop nanonetworks, resource consumption for routing depends on several factors, such as the number of nodes that forward packets, and the number of packets sent and received. To reduce the energy used, in the proposed ray tracing routing, only the nodes that receive duplicate packets at the same time forward them. This leads to the selection of the forwarding nodes over a quasistraight line between the source and the destination, which can bend at the edges of the network. The proposed method is a major improvement over its previous version, aiming to fix several crucial issues such as double propagation, large forwarding angle, and die-out. Moreover, the applicability of ray tracing in heterogeneous and 3D nanonetworks is also analyzed and a practical application for the protocol is proposed. Extensive simulation results demonstrate the fixing of these issues, and the superiority of the proposed improved ray tracing algorithm over other routing methods in the long and narrow nanonetworks, such as blood vessels and branches.

Terahertz channel modeling based on surface sensing characteristics

The dielectric properties of environmental surfaces, including walls, floors and the ground, etc., play a crucial role in shaping the accuracy of terahertz (THz) channel modeling, thereby directly impacting the effectiveness of communication systems. Traditionally, acquiring these properties has relied on methods such as terahertz time-domain spectroscopy (THz-TDS) or vector network analyzers (VNA), demanding rigorous sample preparation and entailing a significant expenditure of time. However, such measurements are not always feasible, particularly in novel and uncharacterized scenarios. In this work, we propose a new approach for channel modeling that leverages the inherent sensing capabilities of THz channels, specifically by obtaining channel measurement data through the analysis of refractive indices. By comparing the results obtained through channel sensing with that derived from THz-TDS measurements, we demonstrate the its ability to yield dependable surface property information. Integrating it into a ray-tracing algorithm for channel modeling in both a miniaturized cityscape scenario and an indoor environment, the results show consistency with experimental measurements, thereby validating its effectiveness in real-world settings.

An adaptive stratification algorithm based on gradient fitting deviation and its application to acoustic ray-tracing algorithm

Accurate underwater positioning is a challenging task for complex marine environments where the study of sound velocity profiles (SVPs) is essential. To overcome the problem that the original SVP reduces the computational efficiency of the ray-tracing algorithm due to great capacity of data, this paper simplifies the SVP. The nature of the SVP and its influence on accurate underwater positioning are discussed. An adaptive stratification algorithm applicable to the ray-tracing algorithm is proposed, based on the gradient fitting deviation and characterized by the fitting treatment of the sound velocity gradient by the Fourier method. Measured SVP data and acoustic positioning data were used to validate this algorithm in the context of acoustic ray-tracing algorithm, to analyze and compare the simplification rates before and after SVP simplification, and the effect on the positioning accuracy. The experimental results in the shallow sea (water depth of 60 m) show that the simplification rate of the SVP is more than 95.56%, and the RMS of the deviation of the ray-tracing algorithm is better than 0.0112 m when determining the optimal number of unfolding levels and threshold; the experimental results in the deep sea (water depth of 1900 m) show that the simplification rate of the SVP is more than 99%, and the RMS of the deviation of the ray-tracing algorithm is better than 0.5 m. This indicates that, the algorithm can reconstruct the original SVP well and implement the automatic hierarchical simplification of the SVP while ensuring the positioning accuracy.

Numerical investigation in an indirect concentrated solar thermal dryer incorporated with discrete mirrors

The objective of this study is to promote the production for small scale drying industries by reducing the drying time of agricultural products. Discrete mirrors are arranged in a parabolic fashion and used as reflectors to concentrate the incoming solar radiation. The available solar flux is doubled by using this discrete arrangement of mirrors. Optical simulations are performed using Monte Carlo Ray tracing algorithm and various effects were studied to select the best arrangement of mirrors. The concentrated type dryer has been numerically simulated and no-load testing is performed using ANSYS Fluent. Simulation is carried out for various inlet velocities from 0.1 to 1 m s−1. The concentrated setup operates between 0.4 and 1 m s−1 of velocity to obtain temperature in the range of 40–70 °C. The maximum and minimum increase in outlet temperature is 23.52% and 2.87%, and the Nusselt number enhancement is 18.94% and 4.15%, respectively.

Clinical application of a GPU-accelerated monte carlo dose verification for cyberknife M6 with Iris collimator

PurposeTo apply an independent GPU-accelerated Monte Carlo (MC) dose verification for CyberKnife M6 with Iris collimator and evaluate the dose calculation accuracy of RayTracing (TPS-RT) algorithm and Monte Carlo (TPS-MC) algorithm in the Precision treatment planning system (TPS).MethodsGPU-accelerated MC algorithm (ArcherQA-CK) was integrated into a commercial dose verification system, ArcherQA, to implement the patient-specific quality assurance in the CyberKnife M6 system. 30 clinical cases (10 cases in head, and 10 cases in chest, and 10 cases in abdomen) were collected in this study. For each case, three different dose calculation methods (TPS-MC, TPS-RT and ArcherQA-CK) were implemented based on the same treatment plan and compared with each other. For evaluation, the 3D global gamma analysis and dose parameters of the target volume and organs at risk (OARs) were analyzed comparatively.ResultsFor gamma pass rates at the criterion of 2%/2 mm, the results were over 98.0% for TPS-MC vs.TPS-RT, TPS-MC vs. ArcherQA-CK and TPS-RT vs. ArcherQA-CK in head cases, 84.9% for TPS-MC vs.TPS-RT, 98.0% for TPS-MC vs. ArcherQA-CK and 83.3% for TPS-RT vs. ArcherQA-CK in chest cases, 98.2% for TPS-MC vs.TPS-RT, 99.4% for TPS-MC vs. ArcherQA-CK and 94.5% for TPS-RT vs. ArcherQA-CK in abdomen cases. For dose parameters of planning target volume (PTV) in chest cases, the deviations of TPS-RT vs. TPS-MC and ArcherQA-CK vs. TPS-MC had significant difference (P < 0.01), and the deviations of TPS-RT vs. TPS-MC and TPS-RT vs. ArcherQA-CK were similar (P > 0.05). ArcherQA-CK had less calculation time compared with TPS-MC (1.66 min vs. 65.11 min).ConclusionsOur proposed MC dose engine (ArcherQA-CK) has a high degree of consistency with the Precision TPS-MC algorithm, which can quickly identify the calculation errors of TPS-RT algorithm for some chest cases. ArcherQA-CK can provide accurate patient-specific quality assurance in clinical practice.

Graphical Heatmap-Based Approach to Indoor Radio Signal Propagation: Adapting Advanced Ray Tracing and Global Illumination Algorithms

Graphical Heatmap-Based Approach to Indoor Radio Signal Propagation: Adapting Advanced Ray Tracing and Global Illumination Algorithms

PyC[formula omitted]Ray: A flexible and GPU-accelerated radiative transfer framework for simulating the cosmic epoch of reionization

Detailed modeling of the evolution of neutral hydrogen in the intergalactic medium during the Epoch of Reionization, 5≤z≤20, is critical in interpreting the cosmological signals from current and upcoming 21-cm experiments such as the Low-Frequency Array (LOFAR) and the Square Kilometre Array (SKA). Numerical radiative transfer codes provide the most physically accurate models of the reionization process. However, they are computationally expensive as they must encompass enormous cosmological volumes while accurately capturing astrophysical processes occurring at small scales (≲Mpc). Here, we present pyC2Ray, an updated version of the massively parallel ray-tracing and chemistry code, C2-Ray, which has been extensively employed in reionization simulations. The most time-consuming part of the code is calculating the hydrogen column density along the path of the ionizing photons. Here, we present the Accelerated Short-characteristics Octahedral ray-tracing (ASORA) method, a ray-tracing algorithm specifically designed to run on graphical processing units (GPUs). We include a modern Python interface, allowing easy and customized use of the code without compromising computational efficiency. We test pyC2Ray on a series of standard ray-tracing tests and a complete cosmological simulation with volume size (349Mpc)3, mesh size of 2503 and approximately 106 sources. Compared to the original code, pyC2Ray achieves the same results with negligible fractional differences, ∼10−5, and a speedup factor of two orders of magnitude. Benchmark analysis shows that ASORA takes a few nanoseconds per source per voxel and scales linearly for an increasing number of sources and voxels within the ray-tracing radii.

Elemental image array generation based on BVH structure combined with spatial partition and display optimization

Integral imaging display has been widely used because of its features such as full parallax viewing, high comfort without glasses, simple structure and easy implementation. This paper enhances the speed of EIA generation by optimizing the acceleration structure of the ray tracing algorithm. By considering the characteristics of segmental objects captured by the camera during integral image rendering, a novel accelerating structure is constructed by combining BVH with the camera space. The BVH traversal is expanded into a 4-tree based on depth priority order to reduce hierarchy and expedite hit point detection. Additionally, the parameters of the camera array are constrained according to the reconstructed three-dimensional (3D) image range, ensuring optimal object coverage on screen. Experimental results demonstrate that this algorithm reduces ray tracing time for hitting triangle grid of collision objects while automatically determining display range for stereo images and adjusting camera parameters accordingly, thereby maximizing utilization of integrated imaging display resources.

Universal simulation of absorption effects for X-ray diffraction in reflection geometry.

Analytical calculations of absorption corrections for X-ray powder diffraction experiments on non-ideal samples with surface roughness, porosity or absorption contrasts from multiple phases require complex mathematical models to represent their material distribution. In a computational approach to this problem, a practicable ray-tracing algorithm is formulated which is capable of simulating angle-dependent absorption corrections in reflection geometry for any given rasterized sample model. Single or multiphase systems with arbitrary surface roughness, porosity and spatial distribution of the phases in any combination can be modeled on a voxel grid by assigning respective values to each voxel. The absorption corrections are calculated by tracing the attenuation of X-rays along their individual paths via a modified shear-warp algorithm. The algorithm is presented in detail and the results of simulated absorption corrections on samples with various surface modulations are discussed in the context of published experimental results.

Comprehensive and strict validation of neutral atmosphere ray tracing algorithms and some improvements for space geodetic radio signals

This study presents a more thorough and rigorous method for the validation of atmospheric ray tracing algorithms used for computing the path delays of space geodetic microwave signals. Romberg integration through a hypothetical atmospheric structure, e.g. the US standard atmosphere, is recommended to derive the reference values of slant total delay (STD) and geometric delay (GD) for validations. Enlighten by the validation results, simplification for the ray tracing algorithm APM-3D and revision for the piecewise linear (PWL) method are proposed and both are proved to be valid. Comparison results show that in a spherically symmetric atmosphere (one-dimensional (1D) media), if the atmosphere is vertically stratified with sufficient height layers, all 1D, two-dimensional (2D), and three-dimensional (3D) ray tracing algorithms tested yield identical STDs and GDs with reference values at an elevation angle of 5°; while with limited height layers, the revised PWL method exhibits better performance than other 1D algorithms considered. The 2D and 3D algorithms are further compared using realistic 3D atmosphere from ERA5 reanalysis data, and it indicates mm to cm level of differences in GDs; while the STD differences range between <0.1 mm to cm level, depending on weather condition and whether shooting is implemented during the tracing. Both STD and GD discrepancies become larger under severe weather, and the magnitude of 2D-3D difference is found to be consistent with the azimuthal variation of STDs. In comparison of results obtained without shooting, shooting procedure tends to decrease the STD differences but has little impact on the differences in GD. The simplified APM-3D can replace the original version with 10−5 mm level of agreements in STD.

Temporal vectorized visibility for direct illumination of animated models

Direct illumination rendering is an important technique in computer graphics. Precomputed radiance transfer algorithms can provide high quality rendering results in real time, but they can only support rigid models. On the other hand, ray tracing algorithms are flexible and can gracefully handle animated models. With NVIDIA RTX and the AI denoiser, we can use ray tracing algorithms to render visually appealing results in real time. Visually appealing though, they can deviate from the actual one considerably. We propose a visibility-boundary edge oriented infinite triangle bounding volume hierarchy (BVH) traversal algorithm to dynamically generate visibility in vector form. Our algorithm utilizes the properties of visibility-boundary edges and infinite triangle BVH traversal to maximize the efficiency of the vector form visibility generation. A novel data structure, temporal vectorized visibility, is proposed, which allows visibility in vector form to be shared across time and further increases the generation efficiency. Our algorithm can efficiently render close-to-reference direct illumination results. With the similar processing time, it provides a visual quality improvement around 10 dB in terms of peak signal-to-noise ratio (PSNR) w.r.t. the ray tracing algorithm reservoir-based spatiotemporal importance resampling (ReSTIR).