iat-cener / tonatiuh

A Monte Carlo ray tracer for the optical simulation of solar concentrating systems
http://iat-cener.github.io/tonatiuh/
GNU General Public License v3.0
55 stars 15 forks source link

how to choose grid divisons in run flux analysis option? #116

Closed AditiiiGarg closed 1 year ago

AditiiiGarg commented 2 years ago

In Tonatiuh, there is an option of run flux analysis, which gives total power absorbed by the receiver along with maximum, minimum, and average flux values. Besides, it also provides a flux profile along the circumference and length of the cylindrical absorber with coordinates. My query is that there is an option for grid divisions in run flux analysis and changing the grid divisions impacts the maximum flux value. So how should I decide which grid divisions to choose? I have attached the two figures for grid division of 20 x 20 where max. flux is 5676.8 W/sq.m and 200 x 20 where the max flux is 19905.1 W/sq.m. Also, can anyone suggest any technique to develop an expression for flux as a function of the angle around the absorber circumference? ( I have tried using excel that creates a flux equation best fit as a polynomial of 6th order for angles ranging from 0-360 degree)

20 x 20 200 x 20

amutuberria commented 1 year ago

The discretization that you should define is completely related to the application and, above all, to the needed precision of your results. It makes sense by changing the number of cells that you are defining to the grid the maximum flux value changes since this value is the highest value between the cells. If we suppose only a single cell, it would cover the entire surface and therefore the maximum flux would be equal to the average flux. Make sure that when you choose the appropriate number of cells, to run all the required number of rays to obtain fully converged solution.

AditiiiGarg commented 1 year ago

Thank you amutuberria