Finite-volume ray tracing using Computational Fluid Dynamics in linear focus CSP applications

The modelling of solar irradiation in concentrated solar power (CSP) applications is traditionally done with ray-tracing methods, e.g. the Monte Carlo method. For the evaluation of CSP receivers, the results from ray-tracing codes are typically used to provide boundary conditions to Computational Fluid Dynamics (CFD) codes for the solution of conjugate heat transfer in the receivers. There are both advantages and disadvantages to using separate software for the irradiation and heat transfer modelling. For traditional ray-tracing methods, advantages are the cost-effectiveness of the Monte Carlo method in modelling reflections from specular surfaces; the ability to statistically assign a sun shape to the rays; the statistical treatment of reflectivity and optical errors (e.g. surface slope errors), to name a few. When considering a complex mirror field and a complex receiver with secondary reflective surfaces, especially with selective coatings to enhance absorption and limit re-radiation losses, standard ray tracers may be limited in specifying emissivity and absorptivity, which are both specular and temperature dependent, and are hence not suitable as radiation analysis tool. This type of scenario can be modelled accurately using CFD, through the finite volume (FV) treatment of the radiative transfer equation (RTE) and a banded spectrum approach at an increased computational cost. This paper evaluates the use of CFD in the form of the commercial CFD code ANSYS Fluent v15 and v16 to model the reflection, transmission and absorption of solar irradiation from diffuse and specular surfaces found in linear CSP applications. 2-D CFD solutions were considered, i.e. line concentration. To illustrate and validate the method, two sources were used. The first source was test cases from literature with published solutions and the second a combined modelling approach where solutions were obtained using both FV and ray tracing (with SolTrace). For all the test cases, good agreement was found when suitable modelling settings were used to limit both ray-effect and false scattering errors.

[1]  Y. Çengel.  Heat and Mass Transfer: Fundamentals and Applications , 2000 .

[2]  Ahmed Amine Hachicha,et al.  Numerical modelling of a parabolic trough solar collector , 2013 .

[3]  J. Muñoz,et al.  Analysis of internal helically finned tubes for parabolic trough design by CFD tools , 2011 .

[4]  Josua P. Meyer,et al.  3-D CFD Modeling of a Slanted Receiver in a Compact Linear Fresnel Plant with Etendue-Matched Mirror Field , 2015 .

[5]  Josua P. Meyer,et al.  A novel computational approach to combine the optical and thermal modelling of Linear Fresnel Collectors using the finite volume method , 2015 .

[6]  Pierre Garcia,et al.  Codes for solar flux calculation dedicated to central receiver system applications : A comparative review , 2008 .

[7]  J. L. Bhagoria,et al.  Heat transfer and fluid flow analysis of solar air heater: A review of CFD approach , 2013 .

[8]  Tae-Kuk Kim,et al.  Radiative transfer in two-dimensional anisotropic scattering media with collimated incidence , 1989 .

[9]  Assensi Oliva,et al.  Numerical simulation of wind flow around a parabolic trough solar collector , 2013 .

[10]  Mohamed B. Gadi.  Design and simulation of a new energy-conscious system (CFD and solar simulation) , 2000 .

[11]  Suhas V. Patankar,et al.  RAY EFFECT AND FALSE SCATTERING IN THE DISCRETE ORDINATES METHOD , 1993 .

[12]  Markus Eck,et al.  Linear Fresnel collector demonstration on the PSA, Part I – Design, construction and quality control , 2008 .

[13]  R. Kneer,et al.  CFD-simulation of a new receiver design for a molten salt solar power tower , 2013 .

[14]  S. Camporeale,et al.  CFD analysis of melting process in a shell-and-tube latent heat storage for concentrated solar power plants , 2016 .

[15]  Ya-Ling He,et al.  Study on optical and thermal performance of a linear Fresnel solar reflector using molten salt as HTF with MCRT and FVM methods , 2015 .

[16]  Josua P. Meyer,et al.  Introducing CFD in the optical simulation of linear Fresnel collectors , 2016 .

[17]  S. Report.  Forms of Approximate Radiation Transport , 2002 .

[18]  H. E. Bennett,et al.  Relation between Surface Roughness and Specular Reflectance at Normal Incidence , 1961 .

[19]  Kenneth E. Torrance,et al.  Radiative heat exchange modeling inside an oven , 2009 .

[20]  M. Moghimi,et al.  Optimization of a trapezoidal cavity absorber for the Linear Fresnel Reflector , 2015 .

[21]  Gilles Flamant,et al.  SOLFAST, a Ray-Tracing Monte-Carlo software for solar concentrating facilities , 2012 .

[22]  Luisa F. Cabeza,et al.  Computational efficiency in numerical modeling of high temperature latent heat storage: Comparison of selected software tools based on experimental data , 2016 .

[23]  P. Gauché,et al.  CFD analysis of solar tower Hybrid Pressurized Air Receiver (HPAR) using a dual-banded radiation model , 2014 .

[24]  M. Modest.  Radiative heat transfer , 1993 .

[25]  Alan W. Weimer,et al.  Rapid High Temperature Solar Thermal Biomass Gasification in a Prototype Cavity Reactor , 2010 .

[26]  Armando C. Oliveira,et al.  Numerical simulation of a trapezoidal cavity receiver for a linear Fresnel solar collector concentrator , 2011 .

[27]  Chiheb Bouden,et al.  A CFD analysis of the flow structure inside a steam ejector to identify the suitable experimental operating conditions for a solar-driven refrigeration system , 2014 .

[28]  J. P. Jessee,et al.  Bounded, High-Resolution Differencing Schemes Applied to the Discrete Ordinates Method , 1997 .

[29]  Pedro J. Coelho,et al.  The role of ray effects and false scattering on the accuracy of the standard and modified discrete ordinates methods , 2002 .

[30]  Ya-Ling He,et al.  A MCRT and FVM coupled simulation method for energy conversion process in parabolic trough solar collector , 2011 .

[31]  Emilio Baglietto,et al.  Modeling direct steam generation in solar collectors with multiphase CFD , 2014 .

[32]  Juan Xiao,et al.  Three-dimensional numerical study of heat transfer characteristics in the receiver tube of parabolic trough solar collector , 2010 .

[33]  Hong-Shun Li.  Reduction of false scattering in arbitrarily specified discrete directions of the discrete ordinates method , 2004 .

[34]  W. H. Reed,et al.  Ray-Effect Mitigation Methods for Two-Dimensional Neutron Transport Theory , 1977 .

[35]  H. Tan,et al.  Radiation performance of dish solar concentrator/cavity receiver systems , 2008 .

[36]  Julio Chaves,et al.  Etendue-matched two-stage concentrators with multiple receivers , 2010 .

[37]  Björn Laumert,et al.  Conjugate heat transfer analysis of an impinging receiver design for a dish-Brayton system , 2015 .

[38]  Aron Dobos,et al.  SolTrace: A Ray-Tracing Code for Complex Solar Optical Systems , 2013 .

[39]  John C. Chai,et al.  Discrete‐Ordinates and Finite‐Volume Methods for Radiative Heat Transfer , 2009 .

[40]  A. Gross,et al.  CFD analysis for solar chimney power plants , 2013 .

[41]  Susan L. Brantley,et al.  Evaluation of surface preparation methods for glass , 2001 .

[42]  Peiwen Li,et al.  Application of phase change materials for thermal energy storage in concentrated solar thermal power plants: A review to recent developments , 2015 .

[43]  Alan W. Weimer,et al.  Evaluation of finite volume solutions for radiative heat transfer in a closed cavity solar receiver for high temperature solar thermal processes , 2013 .