This thesis deals with the development of optical models for solar power tower technology. Specifically, this work is focused on modeling flux mapping and aiming strategies for central receiver systems (CRS). The resulting codes are applicable to CRS design and operation. This dissertation essentially presents four computational models. The first model, on which the rest of the models are built up, computes the flux density distribution incident on any kind of central receiver which is caused by a single heliostat. The procedure relies on the oblique projection of the receiver mesh onto the image plane, where an accurate analytic function, e.g. UNIZAR, is evaluated. Oblique projection is accomplished by transformation of coordinate systems. The 4-step projection method remarkably reproduces the distorted spot found for large incidence angles on the heliostat and the receiver. This basic model was validated against flux measurements on a flat receiver and Monte Carlo Ray Tracing simulations on a cylindrical receiver. Compared to SolTrace, the model takes 50 times less computation time and higher level of resolution. The second model was developed to determine canting errors in the facets of real heliostats. Based on a deterministic optimization algorithm, a procedure was set up to minimize the difference between computed flux maps and captured images on a lambertian target. Experimental images from THEMIS plant were employed to find out canting errors in selected CETHEL heliostats. From results of the model, one of the heliostats was successfully readjusted, significantly improving its optical quality, and validating the proposed methodology. The third model extends the basic model to superpose single heliostat flux maps in a whole field of heliostats. Shading and blocking losses are computed by parallel projection of neighbor heliostats. An aiming strategy, symmetric about the receiver equator, was developed on the basis of only one parameter: k, aiming factor. Nearly single equatorial aiming is achieved with k>3, while k=0 results in pointing to either upper or lower receiver edges. For the Gemasolar case study, an aiming factor equal to 2 yielded the most uniform flux maps, i.e. flat profile in the central region, and negligible increase in spillage losses compared to equatorial aiming. An optimal aiming strategy for molten salt receivers was implemented in the fourth model. An algorithm was developed to maximize receiver thermal output, while meeting at the same time corrosion and thermal stress limits; which were translated into allowable flux densities, AFD. Compared to unreliable single aiming, the optimized aiming strategy ensures receiver integrity and spillage losses only increase up to 4 percentage points. It was found that optimal aim points are, on average, slightly shifted towards the panel entrance. Despite the conflicting demand between adjacent panels in multi-panel receivers with serpentine flow pattern, the fit algorithm performs noticeable matching to the AFD profile. The resulting code takes around 2 minutes in a standard PC to compute the optimal aim points for a field made up of 2650 heliostats.
Read full abstract