Radiative transfer universal 2D - 3D- code raduga 5.1 (P) for multiprocessor computer
Abstract
Problem is presented by the stationary transfer equation with general boundary conditions, describing both the incident radiation and the reflection law. Geometries: (x,y), (r,z), (r,theta), (x,y,z), (r,theta,z). Inhomogeneities of the calculational region are described by elementary geometrical bodies (rectangles, parallelepipeds, prisms, cylinders, spheres and so on). Mesh approximator represents boundaries of these bodies by boundaries of space cells. PL approximation of the scattering indicatrix is used (corresponding errors are unessential at L ~ 20 30 for practically interesting problems). Sources of radiation: inner distributed sources (all geometries), point and ray sources ((r,z) and (x,y,z) geometries); parallel, normal to the boundary surface (all geometries); parallel, inclined to the boundary surface ((x,y,z) geometry). Source Iteration Method is used in Transfer Boundary Problem solving and the representative set of Discrete Ordinates Methods is used for this Problem discretization. Well - known mesh schemes (Step, Diamond Difference, DD/St) and special ones (AWDD, SWDD and others) are realized. Both analytical algorithms to calculate the unscattered radiation in problems with point , ray , parallel sources and semi analytical algorithms to calculate the single scattered radiation in ray source problems are developed. Space meshes are regular. Angle meshes: on the sphere are used triangular (Sn) and rectangular meshes. Parallel algorithms are used MPI procedures and based on the decomposition of the calculation region into subregions. Calculations for each region are realized by the one processor. After each Source Iteration all processors exchange by corresponding boundary values. Code is on Fortan-90 written and run on parallel supercomputer MBC1000M ( http://www/jscc.ru).
- Publication:
-
EGS - AGU - EUG Joint Assembly
- Pub Date:
- April 2003
- Bibcode:
- 2003EAEJA......349B