Ray Tracing Method Research Articles

Semidiurnal Internal Tide Interference in the Northern South China Sea

Multiwave interference plays a crucial role in shaping the spatial variations of internal tides. Based on a combination of in situ mooring and altimeter data, interference of semidiurnal internal tides was investigated in the northern South China Sea. Mooring observations indicate the observed kinetic-to-potential energy ratio and group speed are both relatively lower than the theoretical values of mode-1 semidiurnal internal tides, indicating the presence of partly-standing waves. This is consistent with the altimeter result that the mooring was located at the antinode within the interference pattern formed by the superposition of the westward and southward semidiurnal internal tides from the Luzon Strait and the continental slope of the southern Taiwan Strait. However, the kinetic-to-potential energy ratio and group velocity were notably changed when an anticyclonic eddy passed by the mooring. By employing the ray-tracing method, we identified that mesoscale processes may induce a phase difference in the semidiurnal internal tides between the Luzon Strait and the continental slope of the southern Taiwan Strait. This alteration further leads to changes in the positions of nodes and antinodes within the interference pattern of the semidiurnal internal tides.

Effect of bathymetry data on tsunami wave ray tracing in the western Banten sea

The lack of tsunami wave warning during the Anak Krakatau collapse in 2018 resulted in devastating damage to coastal areas. The tsunami wave arrived in the coastal areas approximately half an hour after the collapse. As a tsunami wave can travel abruptly, an early warning system should warn faster. However, a warning system based on a sophisticated hydrodynamic model in real-time would take time to conduct the numerical tsunami wave simulation. Hence, this study proposes a fast and reliable estimation of tsunami wave propagation through a classic ray tracing analysis. We use two ray-tracing methods to investigate the tsunami wave propagation from Anak Krakatau to the western Banten Sea. The first method follows Snell’s law, considering dispersive waves, and the second assumes non-dispersive waves based on the ray tracing equations, considering the Earth’s sphericity. Both methods are quantitatively evaluated by comparing the travel time measured at Anyer and Marina Jambu. This study finds that non-dispersive wave tracing performs a shorter computational time and slightly better prediction of tsunami wave propagation than dispersive-based wave tracing, with a relative absolute difference of the travel time of 17.9–26.7% in Anyer and 3.6–5.2% in Marina Jambu. This study also confirms the importance of bathymetry validity in wave ray tracing. Two regional bathymetry datasets with a mean difference of less than 5 m result in different wave ray tracing, in which one dataset does not produce the wave ray path towards the Panaitan Island. Based on bathymetric surveys in Anyer and Marina Jambu, the national bathymetry data (BATNAS) has shown its superiority to being used as a bathymetry in the ray-tracing process, with correlations of 81% and 93% in both areas compared to other available bathymetry datasets. We summarize that reliable bathymetry data and the non-dispersive ray tracing method can be used as an initial estimation of tsunami wave propagation efficiently.

Design, analysis, and manufacturing of a glass-plastic hybrid minimalist aspheric panoramic annular lens

We propose a high-performance glass-plastic hybrid minimalist aspheric panoramic annular lens (ASPAL) to solve several major limitations of the traditional panoramic annular lens (PAL), such as large size, high weight, and complex system. The field of view (FoV) of the ASPAL is 360° × (35° ∼ 110°) and the imaging quality is close to the diffraction limit. This large FoV ASPAL is composed of only 4 lenses. Moreover, we establish a physical structure model of PAL using the ray tracing method and study the influence of its physical parameters on compactness ratio. In addition, for the evaluation of local tolerances of annular surfaces, we propose a tolerance analysis method suitable for ASPAL. This analytical method can effectively analyze surface irregularities on annular surfaces and provide clear guidance on manufacturing tolerances for ASPAL. Benefiting from high-precision glass molding and injection molding aspheric lens manufacturing techniques, we finally manufactured 20 ASPALs in small batches. The weight of an ASPAL prototype is only 8.5 g. Our framework provides promising insights for the application of panoramic systems in space and weight-constrained environmental sensing scenarios such as intelligent security, micro-UAVs, and micro-robots.

Overall performance evaluation of a novel optical truncation method for compound parabolic concentrated photovoltaic-thermal system

A variable geometric concentration ratio (VGCR) truncation method was proposed to enhance the performance of compound parabolic concentrated photovoltaic-thermal (CPC-PVT) system in this study. Based on a multi-physics field model of optical-thermal-electrical coupling, numerical investigations using Monte Carlo ray tracing (MCRT) and finite volume method (FVM) were conducted to analyze the effects of the variation rate of truncation ratio (k) and light incident angle (θin) on the overall performance (optical efficiency (ηop), instantaneous thermal efficiency (ηth), instantaneous electrical efficiency (ηele,total), instantaneous electrical power (P), instantaneous exergy efficiency (ηu,total) and economic performance) of CPC-PVT system. Also, the overall performance between the VGCR-CPC-PVT system and the full CPC-PVT (FCPC-PVT) system was compared. The results indicate that, when θin is smaller, there is little disparity in the overall performance between the two systems. However, as k increases, ηop, ηele,total, ηth and ηu,total of VGCR-CPC-PVT system monotonically increase, and the area of the high-temperature region on the PV panel surface decreases. When θin is larger, the VGCR-CPC-PVT system has a smaller shaded area, ηop, ηele,total, P and ηu,total are significantly improved compared with the FCPC-PVT system. Additionally, k has a significant impact on the overall performance of the system, and P, ηu,total and saved reflective surface area (S) are the criteria for selecting k. When using P and ηu,total as the selection criteria, the optimal k value is 0.04. For 0°≤θin≤38°, the average ηop, average ηele,total, average ηth, average P and average ηu,total of the VGCR-CPC-PVT system with k = 0.04 are relatively improved by 1.91 %, 4.59 %, 2.42 %, 3.12 % and 4.04 %, respectively, while the reflector area decreases by 17.5 %, compared with the FCPC-PVT system. The cost-effectiveness of concentrator increases significantly with k.

View-Independent Adjoint Light Tracing for Lighting Design Optimization

Differentiable rendering methods promise the ability to optimize various parameters of 3d scenes to achieve a desired result. However, lighting design has so far received little attention in this field. In this paper, we introduce a method that enables continuous optimization of the arrangement of luminaires in a 3d scene via differentiable light tracing. Our experiments show two major issues when attempting to apply existing methods from differentiable path tracing to this problem: first, many rendering methods produce images, which restricts the ability of a designer to define lighting objectives to image space. Second, most previous methods are designed for scene geometry or material optimization and have not been extensively tested for the case of optimizing light sources. Currently available differentiable ray-tracing methods do not provide satisfactory performance, even on fairly basic test cases in our experience. In this paper, we propose a novel adjoint light tracing method that overcomes these challenges and enables gradient-based lighting design optimization in a view-independent (camera-free) way. Thus, we allow the user to paint illumination targets directly onto the 3d scene or use existing baked illumination data (e.g., light maps). Using modern ray-tracing hardware, we achieve interactive performance. We find light tracing advantageous over path tracing in this setting, as it naturally handles irregular geometry, resulting in less noise and improved optimization convergence. We compare our adjoint gradients to state-of-the-art image-based differentiable rendering methods. We also demonstrate that our gradient data works with various common optimization algorithms, providing good convergence behaviour. Qualitative comparisons with real-world scenes underline the practical applicability of our method.

Design and numerical simulation of a 45 kWel multi-source high-flux solar simulator as a heating source of dielectric materials positioned in a microwave characterization cavity

This study reports the design and modelling of a high-flux solar simulator (HFSS) combined to a resonant cavity for microwave dielectric properties characterization of ceramics at very high temperatures. The HFSS comprises seven xenon arc lamps with ellipsoidal reflectors, delivering a maximal electrical power of 6.5 kWel each. Positioned on a virtual sphere of around 1600 mm radius, the lamps provide concentrated irradiation heating on a sample positioned inside the resonant cavity. The propagation of the lamp’s irradiation to the sample and the cavity is simulated using Monte Carlo ray tracing method. An ANSYS Fluent finite volume method for solving the resulting multimode heat transfer and flow indicate a radiative flux reaching 1.84 MW m−2, and a sample temperature of approximately 1470°C with 10 % non-uniformity. These results suggest that the designed heating system and resonant cavity is suitable for conducting microwave dielectric properties characterization of refracting ceramic materials at high temperatures, addressing a gap in studies focusing on temperatures exceeding 1000°C.

Determination of interparticle radiative heat transfer and radiation absorption for ellipsoidal particle beds using Monte Carlo ray tracing

The calculation of radiative heat transfer between particles in particle-based solar systems presents significant challenges in terms of computational cost and accuracy. This study addresses these challenges by employing the discrete element method to generate random particle packings and the Monte Carlo ray tracing method to investigate the influence of particle shape and volume fraction (ϕ) on interparticle radiation transfer and the absorption of incident radiation. The results reveal that, at low volume fractions, the shape and orientation of ellipsoidal particles have a pronounced impact on direct radiative transfer. Oblate particles exhibit the highest average radiation distribution factor (RDF), and the dispersion of RDF data is more prominent for oblate and elongated prolate particles compared to spheres. As ϕ increases to 0.35, the contribution of indirect transfer from neighboring particles becomes significant for oblate particles, resulting in a reduction in the dispersion of RDF data. Furthermore, when ϕ reaches 0.55, the decisive factor in determining RDF becomes the shielding effect of particles positioned between the emitting and receiving particles. The average RDF tends to converge among particles of different shapes, although oblate and elongated prolate particles still exhibit considerably larger RDF dispersion compared to spherical and near-spherical particles. Additionally, the effective absorptivity of the particle bed for incident radiation decreases as the aspect ratio of particles increases at all volume fractions. Moreover, the transmitted fraction demonstrates an exponential relationship with the depth, and beds composed of spherical particles exhibit better radiation transmittance compared to beds containing non-spherical particles.

ARTPOL: Analytical ray-tracing method for spectro-polarimetric properties of accretion disks around Kerr black holes

Spectro-polarimetric signatures of accretion disks in X-ray binaries and active galactic nuclei contain information on the masses and spins of their central black holes, as well as the geometry of matter in proximity to the compact objects. This information can be extracted by means of X-ray polarimetry. In this work, we present a fast analytical ray-tracing technique for polarized light (ARTPOL) that helps us to obtain the spinning black hole parameters from the observed properties. This technique can replace the otherwise time-consuming numerical ray-tracing calculations for any optically thick or geometrically thin accretion flow. For the purposes of illustration, we considered a standard optically thick, geometrically thin accretion disk in the equatorial plane of the Kerr black hole. We show that ARTPOL proves accurate for dimensionless spin parameter a ≤ 0.94 with a speed that is over four orders of magnitude faster than direct ray-tracing calculations. This approach opens up broader prospects for direct fittings of the spectro-polarimetric data from the Imaging X-ray Polarimetry Explorer.

A conservative first-collision source treatment for ray effect mitigation in discrete-ordinate radiation transport solutions

Deterministic transport codes play a fundamental role in the modeling and simulation of neutron transport. One of the most common deterministic methods is the method of discrete ordinates, also known as the Sn method. While offering significant advantages over other deterministic methods or stochastic methods like Monte Carlo, the method of discrete ordinates suffers from non-physical artifacts in its local solution due to its discretization of angle. These artifacts, referred to as ray effects because of their ray-like appearance, tend to be worse in problems with small sources in areas with little scattering. Significant effort has gone into developing methods to mitigate ray effects, such as the first-collision source treatment, which separates the angular flux into the uncollided and collided fluxes and solves them using non-traditional techniques such as ray tracing. Current ray tracing methods typically trace to a set of points inside a zone to compute an overall flux. However, this approach has significant drawbacks, such as a low potential for optimization, as well as not being conservative. A new method has been developed that traces instead to a set of points on each of a zone's surfaces and computes the currents, before using these to obtain the flux. A comparison between these two ray tracing methods shows significant advantages to the new surface method, including inherent particle conservation, similar convergence rate and increased potential for optimization through ray sharing.

Automatic generation method for long-focal-length unobscured freeform optical systems with small volume

The existing design methods for long-focal-length unobscured freeform systems rarely consider the imaging quality requirements and volume constraints simultaneously, causing most of the final designs to not fulfill the requirement of light weight. This study proposes a method to automatically design a long-focal-length unobscured reflective system that satisfies volume constraints while maintaining high imaging quality. First, a method to adaptively set the structural parameter range is proposed, and multiple parameters for different systemic specifications can be effectively calculated within it. Subsequently, the systemic volume and area functions are constructed using the ray tracing method, where the tilt angles, distances between mirrors, and radii of curvature of the mirrors are chosen as the optimization parameters. Third, a comprehensive objective function is jointly established combining ray obscuration and convergence as performance evaluation factors. Then, the structural parameters of a long-focal-length unobscured system with small volume are easily obtained via the simulated annealing method. Finally, the improved W-W method is used to further enhance the imaging quality of the system, and an unobscured freeform reflective optical system with three mirrors is automatically generated. Experimental results demonstrate that our method can automatically calculate the parameter ranges to facilitate the search for structural parameters, and effectively design the long-focal-length unobscured freeform systems with small volume and high imaging quality.

A RT-FDTD method of analyzing wireless propagation characteristics in underground mine

Efficient communication is crucial in reducing injuries and fatalities in coal mine accidents, necessitating the study of simulation methods for mine communication. When transceiver antennas are positioned close to the same side of the tunnel, the simulation results from the Ray Tracing (RT) method exhibit significant errors. Additionally, the Finite-Difference Time-Domain (FDTD) method demands substantial computational resources. In response to these challenges, we propose a RT-FDTD method, guided by the law of conservation of energy. This approach involves dividing the mine tunnel into a cuboidal region, using the RT method to calculate the electric field strength on the cuboid’s surface, and then employing this as the excitation source for the FDTD method. Subsequently, the FDTD method is used to calculate the electric field strength within the cuboid. Experimental results demonstrate that the RT-FDTD method effectively mitigates the limitations of the RT and FDTD methods, enhancing both the efficiency and accuracy of simulations in underground mine.

Stochastic ray tracing for Fresnel diffraction

We propose stochastic ray tracing for laser beam propagation in Fresnel diffraction to find the duality between wave and ray representations. We transform from the Maxwell equations to the Schrödinger equation for a monochromatic laser beam in the slowly varying envelope approximation. The stochastic ray tracing method interprets this Schrödinger equation as a stochastic process, of an analogy of Nelson’s stochastic mechanics. It can illustrate the stochastic paths and the wavefront of an optical beam. This ray tracing method includes Fresnel diffraction effects naturally. We show its general theoretical construction and numerical tests for a Gaussian laser beam with diffraction, that stochasticity realizes the beam waist around the Rayleigh range.

A New Biomechanical Deformation Response Parameter: Change in Central Corneal Thickness During Air Puff Induced Corneal Deformation

Purpose To investigate the percent change in central corneal thickness (%ΔCCT) during air-puff-induced deformation as an indicator of corneal biomechanical response. Methods Forty ex vivo human eyes from forty donors were imaged using the CorVis ST at experimentally controlled intraocular pressure (IOP) of 10, 20, 30, and 40 mmHg, followed by uniaxial strip testing to calculate tensile modulus. The CorVis ST research software tracked the anterior and posterior cornea edges and determined the dynamic corneal response (DCR) parameters. Eyes were excluded if image quality or posterior tracking issues were present. Custom algorithms were used to calculate CCT during deformation using a ray-tracing method to correct for Scheimpflug and optical distortion within each image. Correlation and stepwise regression analyses between the shape-related DCR parameters and %ΔCCT were conducted. A mixed model analysis was performed to test the effect of IOP and the strongest significant predictors of the stepwise regression on %ΔCCT. The significance threshold was set to p < 0.05. Results Thirty eyes were ultimately analyzed and CCT increased significantly from the pre-deformation state to the highest concavity state at each IOP level (p < 0.001). IOP and multiple shape DCRs were found to be significantly related to %ΔCCT (p < 0.0001). The strongest predictor of %ΔCCT was integrated inverse radius (IIR) (p < 0.0001; partial R2 = 0.4772) with no other parameter having a partial R2 value greater than 0.04. The mixed model analysis showed that IIR was the sole predictor (p = 0.0098) and IOP was no longer significant as a single predictor. However, the interaction of IIR with IOP (p = 0.0023) had a significant effect on %ΔCCT. Conclusion Percent change in CCT is influenced by corneal stiffness as indicated by the significant relationship with IIR. The %ΔCCT may be a potential biomarker for determining differences in corneal deformation response with corneal diseases.

Ray-tracing analytical absorption correction for X-ray crystallography based on tomographic reconstructions

Processing of single-crystal X-ray diffraction data from area detectors can be separated into two steps. First, raw intensities are obtained by integration of the diffraction images, and then data correction and reduction are performed to determine structure-factor amplitudes and their uncertainties. The second step considers the diffraction geometry, sample illumination, decay, absorption and other effects. While absorption is only a minor effect in standard macromolecular crystallography (MX), it can become the largest source of uncertainty for experiments performed at long wavelengths. Current software packages for MX typically employ empirical models to correct for the effects of absorption, with the corrections determined through the procedure of minimizing the differences in intensities between symmetry-equivalent reflections; these models are well suited to capturing smoothly varying experimental effects. However, for very long wavelengths, empirical methods become an unreliable approach to model strong absorption effects with high fidelity. This problem is particularly acute when data multiplicity is low. This paper presents an analytical absorption correction strategy (implemented in new software AnACor) based on a volumetric model of the sample derived from X-ray tomography. Individual path lengths through the different sample materials for all reflections are determined by a ray-tracing method. Several approaches for absorption corrections (spherical harmonics correction, analytical absorption correction and a combination of the two) are compared for two samples, the membrane protein OmpK36 GD, measured at a wavelength of λ = 3.54 Å, and chlorite dismutase, measured at λ = 4.13 Å. Data set statistics, the peak heights in the anomalous difference Fourier maps and the success of experimental phasing are used to compare the results from the different absorption correction approaches. The strategies using the new analytical absorption correction are shown to be superior to the standard spherical harmonics corrections. While the improvements are modest in the 3.54 Å data, the analytical absorption correction outperforms spherical harmonics in the longer-wavelength data (λ = 4.13 Å), which is also reflected in the reduced amount of data being required for successful experimental phasing.

Double‐difference constrained reflection tomography in two‐dimensional elastic media

AbstractThe inverted velocity model obtained from the reflection tomography based on the angle domain common‐image gathers has a certain fuzziness. The inverted velocity model's stratigraphic interface is always not clear enough in areas with complex stratigraphic structure. In order to improve the accuracy and resolution of the inverted velocity model, a double‐difference constraint condition is added on the basis of minimizing the absolute travel‐time residual at the subsurface imaging points. This constraint makes the inverted velocity model local structure information more refined by minimizing the differential travel‐time residual at adjacent imaging points (i.e. closely spaced points within the same layer) and makes the variation of velocity model information within a certain range more accurate. The method in this paper is based on angle domain common‐image gathers, the tomography inversion equation is established by using the ray tracing method, and the conversion relationship between the traveltime residual and the residual curvature of the angle domain common‐image gathers. Then, by adding differential constraint and double‐differential constraint conditions and using the least squares QR decomposition method to solve the set of equations, the inverted velocity model can be obtained through multiple iterations, which provides a high‐precision velocity field for the migration and improves the accuracy of seismic imaging. Numerical experiments on both one typical model and a field data example demonstrate the effectiveness of the proposed double‐difference constrained elastic reflection tomography in generating high‐precision velocity models.

The effect of urban environment on large-scale path loss model’s main parameters for mmWave 5G mobile network in Iraq

Abstract The high speeds resulting from the use of millimeter waves (mmWave) in 5G mobile networks are accompanied by high path loss. The issue of generating a reliable propagation model of radio waves is crucial to the development of cellular networks since it reveals essential information regarding the properties of the wireless channel. The received signal strength, the coverage area, and the outage probability in certain places may all be determined through theoretical or empirical radio frequency propagation models, which offer essential valuable information regarding signal path loss and fading. This work analyzes a comprehensive three-dimensional ray-tracing method at 28 GHz for Najaf city, Iraq. The optimum path loss model for the city of Najaf is evaluated using the close-in (CI) model. On average, the values of the main parameters of CI model n n , X σ CI {X}_{\sigma }^{{\rm{CI}}} accomplished, respectively, 3.461866667 and 11.13958333. The lowest achievable path loss exponent was 3.0619 across all analyzed scenarios, while the highest possible value was 4.1253. The results of this work can serve as a baseline for mmWave measurement campaigns conducted in comparable conditions, and they provide a new avenue for future research into mmWave at 28 GHz in Iraq.

Evaluation of the influence of distortions of the digital object model on the result of synthetic aperture radar image modelling

An essential problem in radar image processing is the synthesis of optimal and quasi-optimal algorithms for the formation of a vector of sufficient statistics and, as a consequence, for obtaining the features used to classify an object of observation. For this purpose it is required to consider the influence of distortions of digital three-dimensional models of objects on the results of modelling images obtained for synthetic aperture radar. Aim of article – to develop and evaluate the possibilities of the method of geometrical optics in the process of application of ray tracing apparatus for the cases of diffuse and specular reflection of waves from the surface of a reflecting object while evaluating the influence of angular distortions of a digital 3D model of an object on the results of modelling. The method of taking into account the theoretical relationship between the magnitude of coordinate distortions of the digital model and the magnitude of distortion of the received radar image is proposed. The paper evaluates the influence of distortions of digital models of three-dimensional (3D) objects on the results of modelling images obtained for synthetic aperture radar. The presented method of taking into account the influence of distortions of the digital 3D model of the object on the results of radar image modelling allows to estimate the conditions of applicability of the ray tracing method.

Influences of tilted thin accretion disks on the observational appearance of hairy black holes in Horndeski gravity

Research on the observational appearance of black holes, both in general relativity and modified gravity, has been in full swing since the Event Horizon Telescope Collaboration announced photos of M87* and Sagittarius A*. Nevertheless, limited attention has been given to the impact of tilted accretion disks on black hole images. This paper investigates the 230 GHz images of non-rotating hairy black holes illuminated by tilted, thin accretion disks in Horndeski gravity with the aid of a ray tracing method. The results indicate that reducing the scalar hair parameter effectively diminishes image luminosity and extends both the critical curve and the inner shadow. This trend facilitates the differentiation between hairy black holes and Schwarzschild black holes, especially in certain parameter spaces where the current Event Horizon Telescope array is capable of capturing such variations. Furthermore, we observe that the inclination of the tilted accretion disk can mimic the observation angle, consequently affecting image brightness and the morphology of the inner shadow. In specific parameter spaces, alterations in the tilt or position of the accretion disk can lead to a drift in the light spot within the images of hairy black holes. This finding may establish a potential correlation between the precession of the tilted accretion disk and image features. Additionally, through an examination of images depicting hairy black holes surrounded by two thin accretion disks, we report the obscuring effect of the accretion environment on the inner shadow of the black hole.

Modeling and control strategy optimizing of solar flux distribution in a four quadrant and adjustable focusing solar furnace

The solar furnace is key equipment for the thermochemical reaction research which uses focused solar energy to drive the reactor, and the solar flux distribution on the reactor surface plays a decisive role in the regulation of the thermochemical reaction process and formation of target products. In this paper, a reflective solar furnace with a thermal power of 30 kWth in Yanqing, Beijing, is taken as the research objective. Firstly, the Monte Carlo Ray Tracing (MCRT) method is used to establish the mathematical model for the optical calculation, and the distribution characteristics of solar flux are deeply analyzed. The accuracy of the model is verified through the comparison with the SolTrace software and experimental measurements respectively. Secondly, an improved four-quadrant and adjustable focusing solar furnace experimental equipment is developed, and a new solar flux distribution regulation strategy is proposed by using the inward crossover and outward expansion offset of the four-quadrant focusing points coupled with the heat-absorbing platform movement. The results show that the non-uniform factor and the peak value of the solar flux distribution in the set target area can be optimized by selecting an appropriate inward crossover offset and movement of the heat-absorbing platform, respectively. Moreover, when the DNI varies between 480 W/m2 and 800 W/m2, through the optimized control strategy, a non-uniformity factor with the value of less than 0.10 can be achieved in the target area of a square reactor surface with a side length of 10 cm, and the peak flux can also be controlled in the range from 3000 kW/m2 to 3300 kW/m2 with a fluctuation range of less than 10.0% as set by the thermochemical reaction. The optimized control strategy can effectively meet the uniformity and stability requirements of solar flux distribution in solar thermochemical reaction and finally achieve the goal of regulating the target products.

Design of heliostat field based on ray tracing

In the tower solar thermal power generation system, the optical efficiency of mirror field is the primary factor that determines the power generation efficiency of the whole system. Given the size and installation height of heliostat, in order to meet the requirements of rated Power 60MW, the parameters of heliostat field are designed and the model of tower solar thermal power generation system is established. Considering the cosine loss, shadow occlusion loss and truncation efficiency, and knowing that these three points are closely related to the position of the Sun, a method of ray tracing is proposed to calculate the sunshine accurately, the optimum parameters of helioscope field were determined by simulation. The simulation results show that the average thermal power output per unit mirror area can be increased by reducing the shadow occlusion between the heliostats.