scholarly journals Optimization and parallelization of the discrete ordinate method for radiation transport simulation in OpenFOAM: Hierarchical combination of shared and distributed memory approaches

2021 ◽  
Vol 1 ◽  
pp. 2
Author(s):  
Jose Moreno-SanSegundo ◽  
Cintia Casado ◽  
David Concha ◽  
Antonio S. Montemayor ◽  
Javier Marugán

This paper describes the reduction in memory and computational time for the simulation of complex radiation transport problems with the discrete ordinate method (DOM) model in the open-source computational fluid dynamics platform OpenFOAM. Finite volume models require storage of vector variables in each spatial cell; DOM introduces two additional discretizations, in direction and wavelength, making memory a limiting factor. Using specific classes for radiation sources data, changing the store of fluxes and other minor changes allowed a reduction of 75% in memory requirements. Besides, a hierarchical parallelization was developed, where each node of the standard parallelization uses several computing threads, allowing higher speed and scalability of the problem. This architecture, combined with optimization of some parts of the code, allowed a global speedup of x15. This relevant reduction in time and memory of radiation transport opens a new horizon of applications previously unaffordable.

2017 ◽  
Vol 4 (2) ◽  
pp. 177-181
Author(s):  
A. Mazaheri ◽  
J. Y. Trépanier ◽  
R. Camarero ◽  
P. Robin-Jouan

In this paper, we propose to revisit the method to solve the radiation transport equation in circuit breakers to reduce the computation time. It is based on an explicit approach using a space marching algorithm. The method can further be accelerated using a Cartesian grid and using the axisymmetric assumption. Comparisons performed in terms of accuracy and efficiency between the P1 model, the implicit finite-volume discrete ordinate method and the space-marching finite-volume discrete ordinate method show that the explicit approach is more that an order of magnitude faster than the implicit approach, for the same accuracy.


2007 ◽  
Vol 64 (11) ◽  
pp. 3854-3864 ◽  
Author(s):  
K. Franklin Evans

Abstract The spherical harmonics discrete ordinate method for plane-parallel data assimilation (SHDOMPPDA) model is an unpolarized plane-parallel radiative transfer forward model, with corresponding tangent linear and adjoint models, suitable for use in assimilating cloudy sky visible and infrared radiances. It is derived from the spherical harmonics discrete ordinate method plane-parallel (SHDOMPP, also described in this article) version of the spherical harmonics discrete ordinate method (SHDOM) model for three-dimensional atmospheric radiative transfer. The inputs to the SHDOMPPDA forward model are profiles of pressure, temperature, water vapor, and mass mixing ratio and number concentration for a number of hydrometeor species. Hydrometeor optical properties, including detailed phase functions, are determined from lookup tables as a function of mass mean radius. The SHDOMPP and SHDOMPPDA algorithms and construction of the tangent-linear and adjoint models are described. The SHDOMPPDA forward model is validated against the Discrete Ordinate Radiative Transfer Model (DISORT) by comparing upwelling radiances in multiple directions from 100 cloud model columns at visible and midinfrared wavelengths. For this test in optically thick clouds the computational time for SHDOMPPDA is comparable to DISORT for visible reflection, and roughly 5 times faster for thermal emission. The tangent linear and adjoint models are validated by comparison to finite differencing of the forward model.


Author(s):  
Jian Jiao ◽  
Zhixiong Guo

Thermal analysis of a cylindrical tissue subject to a train of ultrashort pulse irradiations was made by developing a combined time-dependent radiation and conduction bio-heat transfer model. Ultrashort pulsed radiation transport in the cylindrical tissue is simulated using the transient discrete ordinate method. Treatment of focused beam is introduced. The model skin tissue is stratified as three layers with different optical, thermal and physiological properties. Comparisons between the collimated irradiation and focused beam are conducted. The effect of pulse train is investigated.


2021 ◽  
Vol 13 (2) ◽  
pp. 270
Author(s):  
Adrian Doicu ◽  
Dmitry S. Efremenko ◽  
Thomas Trautmann

An algorithm for the retrieval of total column amount of trace gases in a multi-dimensional atmosphere is designed. The algorithm uses (i) certain differential radiance models with internal and external closures as inversion models, (ii) the iteratively regularized Gauss–Newton method as a regularization tool, and (iii) the spherical harmonics discrete ordinate method (SHDOM) as linearized radiative transfer model. For efficiency reasons, SHDOM is equipped with a spectral acceleration approach that combines the correlated k-distribution method with the principal component analysis. The algorithm is used to retrieve the total column amount of nitrogen for two- and three-dimensional cloudy scenes. Although for three-dimensional geometries, the computational time is high, the main concepts of the algorithm are correct and the retrieval results are accurate.


2018 ◽  
Vol 2018 ◽  
pp. 1-16 ◽  
Author(s):  
O. Alvarado-Rolon ◽  
R. Natividad ◽  
R. Romero ◽  
L. Hurtado ◽  
A. Ramírez-Serrano

This work focuses on modeling and simulating the absorption and scattering of radiation in a photocatalytic annular reactor. To achieve so, a model based on four fluxes (FFM) of radiation in cylindrical coordinates to describe the radiant field is assessed. This model allows calculating the local volumetric rate energy absorption (LVREA) profiles when the reaction space of the reactors is not a thin film. The obtained results were compared to radiation experimental data from other authors and with the results obtained by discrete ordinate method (DOM) carried out with the Heat Transfer Module of Comsol Multiphysics® 4.4. The FFM showed a good agreement with the results of Monte Carlo method (MC) and the six-flux model (SFM). Through this model, the LVREA is obtained, which is an important parameter to establish the reaction rate equation. In this study, the photocatalytic oxidation of benzyl alcohol to benzaldehyde was carried out, and the kinetic equation for this process was obtained. To perform the simulation, the commercial software COMSOL Multiphysics v. 4.4 was employed.


Sign in / Sign up

Export Citation Format

Share Document