RADIATION HEAT TRANSFER- SHAPE FACTOR CALCULATION (MECHANICAL ENGINEERING)
sfactor.m ; (example) xxx = sfactor(a).
calculates shape factor between two flat surfaces with arbitrary orientation and location.
(sfactor12.m and sfactor16.m are similar to sfactor.m but with 12 and 16 gauss-point integration. They are added newly for precission improvement. Check.dat is a sample data file. )
( Note: open(?radiation.fig?) command will give you the details of matrix a- radiation.fig file is newly added. It is a description with a figure)
It gives a1*f12, a2*f21, a1, a2, f12,f21. (a1-area of emitter, a2-area of receiver, f12- emitter factor, f21- receiver factor)
The gemetry is homotopically parametrized. The integration is achieved through numerical methods. (Gauss ?8-poits).
Two quadrilateral plates are considered for computations- an emitter and a receiver.
The coordinates of the nodes of the emitter and receiver form the matrix A. Vector methods are used to calculate necessary sub calculations. They are in the subroutines
Dot.m, cross.m , area.m, fun.m gass.dat.
For triangular plates, two nodes of the quadrilateral must be merged, so that two colums will have the same coordinates. The colums of A must be along a particular orientation. A is 3X8 matrix. All nodes are given 3 dimensional coordinates.
(Note : Numerical integration has its own limits, the user can check for the known geometry and decide)
It is a more generalized version, in the sense that any geometry can be factored into triangular or quadrilateral simplexes and then the data can be fed to this program..
Mathan Chandran (2021). Radiation Heat Transfer Shape Factor (Mechnical Engineering) (https://www.mathworks.com/matlabcentral/fileexchange/8853-radiation-heat-transfer-shape-factor-mechnical-engineering), MATLAB Central File Exchange. Retrieved .
Very useful routine
Find the treasures in MATLAB Central and discover how the community can help you!Start Hunting!