Key Takeaways & Executive Findings
- •• • Single irregular heliostats achieve an optical efficiency of 0.8022, outperforming square (0.7372), pentagonal (0.7453), hexagonal (0.7485), heptagonal (0.7491), octagonal (0.7502), and circular (0.7513) geometries by 6.8–8.8% relative improvement, directly reducing the mirror area required per unit of thermal power and lowering capital expenditure on heliostat fields. • • The reverse projection method eliminates ray-intersection checks with the receiver by transforming the problem into ray-heliostat intersection tests, reducing computational overhead for optical efficiency evaluation; this enables rapid iteration over thousands of heliostat shapes and layouts without sacrificing accuracy, as validated by the consistency of the north-high/south-low and center-high/edge-low energy density patterns with established models. • • The no-blocking dense layout, combining Campo and EB arrangements, avoids the need to compute shadowing and blocking arrays, simplifying the efficiency calculation to a product of effective flux density, shadowing array, and shape array; this reduces simulation complexity while maintaining a stable optical efficiency of 0.8022 for irregular heliostats, demonstrating robustness across varying field positions. • • The effective energy flux density distribution on heliostats, derived by reordering the total power integral to integrate over the receiver first, provides a direct spatial map of irradiance that can be thresholded to define optimal heliostat shapes; this data-driven shaping yields a 0.8022 efficiency versus 0.7372 for squares, translating to a 8.8% reduction in mirror area for equivalent power and a corresponding decrease in land use and structural costs.
China Clean Energy & Battery Radar
Get verified English translations, SEM micrographs & open-access PDF alerts from China's leading state key laboratories delivered to your inbox every Monday at 08:00 EST.
Abstract
This study addresses the computational bottleneck in heliostat field optical efficiency assessment by proposing a reverse projection method for calculating effective energy flux density on heliostat surfaces. The method replaces conventional cosine efficiency and truncation efficiency calculations with an irradiance function, determines shadowing and blocking occurrences, computes single heliostat power, and accumulates total field power. By reordering the summation in the total power integral, the method derives the energy flux density distribution on the receiver surface and, through an alternative reordering, the effective energy flux density distribution on each heliostat. This distribution enables determination of the heliostat shape that maximizes power under given area or other constraints, and rapid evaluation of optical efficiency across different shapes and layouts. Numerical simulations demonstrate that heliostats shaped according to the proposed algorithm require smaller mirror areas to achieve equivalent power output compared to conventional rectangular and polygonal geometries, while maintaining superior stability. The study employs a no-blocking dense layout combining Campo and EB arrangements. Results show that a single irregular heliostat achieves an optical efficiency of 0.8022, significantly exceeding square (0.7372), pentagonal (0.7453), hexagonal (0.7485), heptagonal (0.7491), octagonal (0.7502), and circular (0.7513) configurations. The method also reveals that heliostats closer to the receiver exhibit higher energy flux density, with the field energy density in the northern hemisphere displaying a north-high/south-low and center-high/edge-low pattern, consistent with other modeling approaches.
1. Introduction
Concentrated solar power tower plants rely on heliostat fields to concentrate direct normal irradiance onto a central receiver, yet the optical efficiency of these fields is governed by a complex interplay of cosine losses, shadowing and blocking, truncation, and atmospheric attenuation. Conventional modeling approaches—such as Monte Carlo ray tracing, two-dimensional normal distribution fitting, and microintegration—impose substantial computational burdens when evaluating thousands of heliostats, particularly for non-rectangular geometries. Commercial software like SolarPilot mitigates this by partitioning the field into zones and applying a single heliostat's truncation efficiency to each zone, but this sacrifices spatial fidelity and cannot resolve the effective flux density distribution on individual heliostat surfaces. Consequently, the design of heliostat shapes that minimize mirror area while maximizing intercepted power remains an unresolved bottleneck, with most commercial fields defaulting to rectangular or polygonal mirrors that are suboptimal for the local flux profile.
This study introduces a reverse projection method that inverts the optical path: instead of tracing rays from the sun to the receiver, it projects the receiver onto each heliostat and computes the effective energy flux density distribution directly. By replacing cosine and truncation efficiency calculations with an irradiance function and incorporating shadowing and blocking as a binary array, the method reduces the total power computation to a sum of Hadamard products over heliostat surface elements. Reordering the resulting integral yields both the receiver surface flux distribution and the heliostat surface effective flux density, enabling threshold-based shape optimization. The approach is validated on a no-blocking dense layout combining Campo and EB arrangements, and the resulting irregular heliostat shape achieves an optical efficiency of 0.8022—significantly higher than conventional geometries—while maintaining stability and reducing required mirror area for equivalent power output.
Loading authentic research manuscript (Pages 1–5)...
LIU Jinzi, LI Chentao, GUO Han (2026). A Method for Heliostat Field Modeling and Effective Energy Flux Density Calculation. Acta Energiae Solaris Sinica. https://doi.org/10.19912/j.0254-0096.tynxb.202608_9684
Research & Educational Purpose Only: The translations, structured abstracts, analytical annotations, and data reports provided by SinoGreenTechare intended exclusively for academic research, internal corporate R&D, and educational benchmarking. They do not constitute formal engineering, chemical safety, legal, or professional advice.
Copyright & Intellectual Property Notice: Original copyright of the underlying source articles and experimental data remains with the respective authors, institutions, and original publishing journals. SinoGreenTech claims intellectual property only over its proprietary translations, analytical syntheses, and AEO structured enhancements in accordance with international fair use and academic citation principles.
Frequently Asked Questions
What is the quantified optical efficiency gain of the proposed irregular heliostat shape compared to conventional geometries, and what is the industrial implication for heliostat field capital cost?
The irregular heliostat achieves an optical efficiency of 0.8022, versus 0.7372 for square, 0.7453 for pentagonal, 0.7485 for hexagonal, 0.7491 for heptagonal, 0.7502 for octagonal, and 0.7513 for circular heliostats. This represents a 6.8–8.8% relative improvement over the best conventional shape (circular). Industrially, a higher optical efficiency directly reduces the mirror area required to achieve a given thermal power at the receiver, lowering heliostat field capital expenditure by a proportional amount, assuming mirror cost dominates field cost. The method also enables rapid evaluation of shape and layout trade-offs, potentially reducing design cycle time.
How does the reverse projection method reduce computational cost compared to conventional ray tracing, and what is the trade-off in accuracy?
Conventional ray tracing requires checking each ray for intersection with the receiver, which scales with the number of rays and receiver complexity. The reverse projection method transforms this into a ray-heliostat intersection problem by projecting the receiver onto the heliostat surface, eliminating the need for receiver intersection tests. This reduces the computational burden from O(N_rays × N_receiver_elements) to O(N_rays × N_heliostat_elements), where N_heliostat_elements is typically much smaller. Accuracy is maintained because the method computes the effective flux density distribution directly, and the resulting energy density patterns (north-high/south-low, center-high/edge-low) match those from established models, as validated in the study.
What are the failure mechanisms or limitations of the no-blocking dense layout used in this study, and how does it affect the generalizability of the optical efficiency results?
The no-blocking dense layout, combining Campo and EB arrangements, is designed to avoid shadowing and blocking by spacing heliostats such that no heliostat blocks the reflected beam of another. This eliminates the need to compute shadowing and blocking arrays, simplifying the efficiency calculation. However, this layout may require larger land area per heliostat, increasing land cost. The optical efficiency of 0.8022 for irregular heliostats was obtained under this layout; in denser layouts where blocking occurs, the efficiency would be lower and the optimal shape might differ. The method itself can incorporate shadowing and blocking via the binary array B_{t,i,j,k}, so the limitation is not inherent to the algorithm but to the specific layout chosen for validation.
How does the effective energy flux density distribution vary spatially across the heliostat field, and what operational thresholds should be used for shape optimization?
The computed effective energy flux density is higher for heliostats closer to the receiver, and the field energy density exhibits a north-high/south-low and center-high/edge-low pattern in the northern hemisphere, consistent with other modeling methods. For shape optimization, the study demonstrates three different threshold cuts on the flux density map to define heliostat shapes. The optimal threshold depends on the trade-off between mirror area and intercepted power; a higher threshold yields a smaller, more efficient shape but may increase manufacturing complexity. The irregular shape achieving 0.8022 efficiency was derived from one such threshold, but the method allows rapid evaluation of multiple thresholds to identify the shape that maximizes power under given area constraints.
What is the scalability bottleneck for implementing the reverse projection method in a commercial-scale heliostat field with tens of thousands of heliostats, and how can it be mitigated?
The primary scalability bottleneck is the computation of the four-dimensional array J_{t,i,j,k} for each time point t, heliostat i, and surface element (j,k). For a field with N heliostats and M surface elements per heliostat, the memory and computational requirements scale as O(T × N × M), where T is the number of time points. For a commercial field with 10,000 heliostats and 100 elements each, this is 10^6 elements per time point, which is manageable with modern GPUs. The method can be further optimized by exploiting the spatial smoothness of J across neighboring heliostats, as noted in the paper, allowing interpolation from a subset of heliostats. Additionally, the no-blocking layout reduces the need for shadowing calculations, further lowering computational load. Parallelization across heliostats and time points is straightforward, enabling real-time or near-real-time optimization.
Related Chinese Research & Cross-Citations
Wind Turbine Gearbox Fault Diagnosis Method Based on Improved CNN-XGBoost Fusion Model Under Gramian Angular Difference Field
Gearbox failures account for 20–30% of total wind turbine faults and incur maintenance costs equivalent to 10–15% of overall turbine value. Conventional vibration diagnostic pipelines—complementary ensemble empirical mode decomposition with singular value energy spectrum, time-varying filtering empirical mode decomposition, and Teager energy spectrum analysis—remain bounded below 90% accuracy and depend on expert-driven feature engineering that is sensitive to non-stationary operating conditions and noise. This study proposes an intelligent diagnostic architecture that converts one-dimensional gearbox vibration signals into two-dimensional images via Gramian angular difference field (GADF) transformation, preserving intrinsic temporal correlation and time-frequency structure while exploiting matrix sparsity to suppress interference. An improved convolutional neural network (CNN) extracts multi-dimensional features: a convolutional block attention module (CBAM) is embedded in the convolutional layers to weight critical channels and focus on fault-sensitive spatial regions, and a modified βc-ACONC activation function replaces ReLU to mitigate neuron necrosis and enable selective activation. The extracted composite features are then fed into an XGBoost network whose hyperparameters are optimized by an improved sparrow search algorithm (ISSA). Validation on a laboratory wind turbine gearbox dataset yields diagnostic accuracy exceeding 99%, demonstrating robust fault identification capability under complex operating conditions.
Joint Forecasting of Wind and Photovoltaic Power Considering Complementarity
The inherent spatiotemporal complementarity between wind and solar resources offers a theoretical basis for improving renewable power forecasting accuracy. This study proposes a joint wind-photovoltaic (PV) power forecasting strategy that explicitly exploits this complementarity. A bidirectional long short-term memory (BiLSTM) neural network serves as the baseline forecasting model, and a novel sorting and comparative optimization (SCO) algorithm is developed to optimize the model's hyperparameters. The SCO algorithm ranks individuals in ascending order and compares adjacent fitness values to escape local optima, a known deficiency in conventional metaheuristics such as genetic algorithms and particle swarm optimization. For wind farms and PV plants exhibiting significant complementarity, the joint forecasting strategy first aggregates their power outputs, normalizes the combined signal, and then feeds it into the optimized BiLSTM model. Experimental results demonstrate that the proposed SCO-BiLSTM model reduces the eMAPE by 10.313% compared with PSO-BiLSTM. Furthermore, joint forecasting under SCO-BiLSTM lowers the eRMSE by 27.443% relative to standalone PV power forecasting. The study also establishes that forecasting accuracy improves with stronger wind-solar complementarity but degrades as the forecasting horizon extends. These findings confirm that exploiting complementarity in joint forecasting substantially enhances predictive performance for renewable energy integration.
Ultra-Short-Term Wind Power Forecasting Based on Fluctuation Continuation Scenario Identification
Existing ultra-short-term wind power forecasting methods exhibit limited performance due to insufficient extraction of fluctuation information and inadequate analysis of evolution patterns. This paper proposes an ultra-short-term wind power forecasting method based on fluctuation continuation scenario identification. First, the coupling mechanism of wind power fluctuations under multiple turbulence processes is investigated, and historical wind power dynamics are decoupled into a combination of nonlinear and linear fluctuation components. A fluctuation continuation concept is introduced, and the future continuation scale of wind power fluctuations is derived from nonlinear and linear decoupling parameters, thereby classifying fluctuation continuation scenarios. A sparse neural network (SNN) oriented to high-dimensional sparse features is constructed to identify historical fluctuation continuation scenarios, and ultra-short-term power forecasting is conducted separately for each scenario. Validation using measured wind speed and power data from three wind farms shows that, compared with baseline models, the proposed model improves root mean square error (RMSE), mean absolute error (MAE), and mean absolute percentage error (MAPE) by at least 1.46%, 2.44%, and 14.67%, respectively, demonstrating superior accuracy and stability. The method addresses the limitations of signal decomposition techniques that lack physical interpretability and are sensitive to hyperparameters, and overcomes the high-dimensional sparsity challenges faced by traditional scenario identification models.
Multi-Classifier Open Adversarial Network for Rolling Bearing Fault Diagnosis in Wind Turbine Generator Systems
Rolling bearings in wind turbine generator systems operate under variable speed and load conditions that induce significant data distribution shifts between training and field data, while unknown fault modes absent from the source domain are frequently misclassified as known classes. This study proposes a multi-classifier open adversarial network (MCOAN) for open-set fault diagnosis. Within an adversarial domain adaptation framework, a K-way classifier and an additional K+1-way one-vs-all classifier independently estimate target-sample similarity to the source domain. These similarity scores drive a dynamic weighting mechanism that adaptively reweights target samples during open-set adversarial training and supplies per-sample dynamic thresholds for known/unknown discrimination, thereby promoting cross-domain alignment of shared-class features while suppressing negative transfer from unknown samples. A non-adversarial domain classifier is introduced to stabilize dynamic weight estimation. Validation on two datasets demonstrates high-precision shared-class distribution alignment and unknown-class recognition with favorable robustness. The method removes reliance on empirically preset thresholds that plague conventional open-set back-propagation approaches, where a fixed threshold of 0.5 in the binary cross-entropy adversarial loss provides no per-sample adaptivity. By coupling K-way and K+1-way similarity estimates, MCOAN achieves simultaneous known-class alignment and unknown-class separation without prior knowledge of the unknown-class cardinality, addressing a persistent bottleneck in wind turbine drivetrain condition monitoring where unanticipated bearing failure modes emerge under field conditions.
Unsupervised Automated Identification Method for Abnormal States of Wind Turbine Gearboxes
Addressing the scarcity of labeled data for training classification models in wind turbine planetary gearbox anomaly identification, this study proposes an unsupervised automated detection method. Log Mel-band energy features are extracted from raw vibration signals and fed into an unsupervised anomaly recognition model centered on a U-net autoencoder. A health-state threshold is established based on reconstruction error between model input and output, enabling anomaly identification. The method is validated using factory gearbox test data and operational data from a wind farm in Yangtouya, Shanxi. For factory gearboxes, dual validation is performed using a spectrum amplitude modulation-based signal processing method. Results demonstrate that the proposed method achieves 93.34% recognition accuracy on both factory and wind farm test sets, confirming its capability to automatically and correctly separate abnormal wind turbine gearboxes. The approach eliminates reliance on labeled fault data, offering a scalable solution for full-lifecycle health monitoring, from factory acceptance testing to in-service early anomaly detection, adaptable across different operating conditions and turbine models.
Improved Adaptive Super-Twisting Sliding Mode Control for Permanent Magnet Synchronous Motors
Adaptive super-twisting sliding mode control (ASTSMC) for permanent magnet synchronous motors (PMSM) suffers from prolonged convergence and insufficient disturbance rejection under complex operating conditions. This paper proposes an improved ASTSMC incorporating a fixed-time disturbance observer (FTDO). A power term is introduced into the adaptive super-twisting controller to accelerate convergence far from the origin, while the discontinuous sign function is replaced by a continuous h(s) function to mitigate chattering. The FTDO ensures disturbance estimation converges within a fixed time independent of initial states, overcoming the limitations of traditional and finite-time observers. The estimated disturbance is fed forward to the sliding mode controller for compensation, enhancing robustness. The FTDO design is based on an auxiliary state variable z and its derivative, with error dynamics analyzed via Lyapunov stability. Comparative simulations against conventional disturbance observers and sliding mode controllers validate the proposed strategy. The results demonstrate shorter convergence time, improved dynamic performance, and superior disturbance rejection, making the approach suitable for high-performance servo systems. The method addresses the critical need for robust, fast-response control in electric vehicles, rail transit, and aerospace applications where PMSM drives face significant uncertainties and external disturbances.