Spectral spherical harmonics discrete ordinate method

Author(s):  
Adrian Doicu ◽  
Michael I. Mishchenko ◽  
Dmitry S. Efremenko ◽  
Thomas Trautmann
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.


2020 ◽  
Vol 12 (22) ◽  
pp. 3703
Author(s):  
Adrian Doicu ◽  
Dmitry S. Efremenko ◽  
Thomas Trautmann

A spectral acceleration approach for the spherical harmonics discrete ordinate method (SHDOM) is designed. This approach combines the correlated k-distribution method and some dimensionality reduction techniques applied on the optical parameters of an atmospheric system. The dimensionality reduction techniques used in this study are the linear embedding methods: principal component analysis, locality pursuit embedding, locality preserving projection, and locally embedded analysis. Through a numerical analysis, it is shown that relative to the correlated k-distribution method, PCA in conjunction with a second-order of scattering approximation yields an acceleration factor of 12. This implies that SHDOM equipped with this acceleration approach is efficient enough to perform spectral integration of radiance fields in inhomogeneous multi-dimensional media.


Atmosphere ◽  
2019 ◽  
Vol 10 (6) ◽  
pp. 292 ◽  
Author(s):  
Adrian Doicu ◽  
Dmitry S. Efremenko

Linearizations of the spherical harmonic discrete ordinate method (SHDOM) by means of a forward and a forward-adjoint approach are presented. Essentially, SHDOM is specialized for derivative calculations and radiative transfer problems involving the delta-M approximation, the TMS correction, and the adaptive grid splitting, while practical formulas for computing the derivatives in the spherical harmonics space are derived. The accuracies and efficiencies of the proposed methods are analyzed for several test problems.


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.


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.


Sign in / Sign up

Export Citation Format

Share Document