Improved Discrete Ordinates Method for Ray Effects Mitigation

2011 ◽  
Vol 133 (4) ◽  
Author(s):  
Zhi-Feng Huang ◽  
Huai-Chun Zhou ◽  
Pei-feng Hsu

A new and improved method based on the concept of discrete ordinates scheme with infinitely small weights (DOS+ISW) is developed for modeling radiative heat transfer in three-dimensional participating media. To demonstrate the effectiveness of the method in mitigating ray effects, the ray effects caused by (1) abrupt step changes in the boundary conditions and (2) the stepwise variation of the medium emissive power are considered. In this work, angular quadrature sets with large number of discrete ordinate directions are chosen to mitigate ray effects while at the same time keeping the computational time increase to a minimum. Comparing with the conventional discrete ordinates method, the difference is that intensities in these directions are calculated by DOS+ISW method. Intensity with fine directional resolution calculated by this method is validated by comparing with that of reverse Monte Carlo method. The large number of discrete ordinate directions used in the new method becomes computationally prohibitive in the conventional discrete ordinates method due to the increased computer memory and computation time requirements.

Author(s):  
Zhi-Feng Huang ◽  
Huai-Chun Zhou ◽  
Pei-feng Hsu

A new and improved method based on the discrete ordinates scheme with infinitely small weights (DOS+ISW) is developed for radiative heat transfer in three-dimensional participating media. To demonstrate the effectiveness of the method, ray effects caused by 1) abrupt step changes in the boundary conditions and 2) the stepwise variation of the medium emissive power are discussed. In this work, angular quadrature sets with large number of discrete ordinate directions are chosen to mitigate ray effects, while at the same time keeping the computational time increase to a minimum. Comparing with the conventional discrete ordinates method, the difference is that intensities in these directions are calculated by DOS+ISW. Intensity with fine directional resolution calculated by this method is validated by comparing with that of Reverse Monte Carlo method. The large number of discrete ordinates directions used in the new method becomes computationally prohibitive in discrete ordinates method due to the increased computer memory and computation time requirements.


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.


2021 ◽  
Author(s):  
Brett W. Larsen ◽  
Shaul Druckmann

AbstractLateral and recurrent connections are ubiquitous in biological neural circuits. The strong computational abilities of feedforward networks have been extensively studied; on the other hand, while certain roles for lateral and recurrent connections in specific computations have been described, a more complete understanding of the role and advantages of recurrent computations that might explain their prevalence remains an important open challenge. Previous key studies by Minsky and later by Roelfsema argued that the sequential, parallel computations for which recurrent networks are well suited can be highly effective approaches to complex computational problems. Such “tag propagation” algorithms perform repeated, local propagation of information and were introduced in the context of detecting connectedness, a task that is challenging for feedforward networks. Here, we advance the understanding of the utility of lateral and recurrent computation by first performing a large-scale empirical study of neural architectures for the computation of connectedness to explore feedforward solutions more fully and establish robustly the importance of recurrent architectures. In addition, we highlight a tradeoff between computation time and performance and demonstrate hybrid feedforward/recurrent models that perform well even in the presence of varying computational time limitations. We then generalize tag propagation architectures to multiple, interacting propagating tags and demonstrate that these are efficient computational substrates for more general computations by introducing and solving an abstracted biologically inspired decision-making task. More generally, our work clarifies and expands the set of computational tasks that can be solved efficiently by recurrent computation, yielding hypotheses for structure in population activity that may be present in such tasks.Author SummaryLateral and recurrent connections are ubiquitous in biological neural circuits; intriguingly, this stands in contrast to the majority of current-day artificial neural network research which primarily uses feedforward architectures except in the context of temporal sequences. This raises the possibility that part of the difference in computational capabilities between real neural circuits and artificial neural networks is accounted for by the role of recurrent connections, and as a result a more detailed understanding of the computational role played by such connections is of great importance. Making effective comparisons between architectures is a subtle challenge, however, and in this paper we leverage the computational capabilities of large-scale machine learning to robustly explore how differences in architectures affect a network’s ability to learn a task. We first focus on the task of determining whether two pixels are connected in an image which has an elegant and efficient recurrent solution: propagate a connected label or tag along paths. Inspired by this solution, we show that it can be generalized in many ways, including propagating multiple tags at once and changing the computation performed on the result of the propagation. To illustrate these generalizations, we introduce an abstracted decision-making task related to foraging in which an animal must determine whether it can avoid predators in a random environment. Our results shed light on the set of computational tasks that can be solved efficiently by recurrent computation and how these solutions may appear in neural activity.


Author(s):  
Andrew M. Feldick ◽  
Gopalendu Pal

Abstract The introduction of higher fidelity spectral models into a Discrete Ordinates Method (DOM) RTE solver introduces the challenge of solving the N(N+2) coupled equations in intensity over many spectral points. The inability to store intensity fields leads to a nonlinear increase in computational cost as compared to basic gray models, as the solution in an evolving field must be recalculated at each radiation time step. In this paper an approximate initialization approach is used to a reconstructed values of the intensities. This approach is particularly well suited to spectrally reordered methods, as the boundary conditions and scattering coefficients are gray. This approach leads to more tractable computational time, and is demonstrated using on two industrial scale flames.


2019 ◽  
Vol 20 (2) ◽  
pp. 89
Author(s):  
Gede A Widyadana ◽  
Audrey Tedja Widjaja ◽  
Kun Jen Wang

A single container loading problem is a problem to effectively load boxes in a three-dimensional container. There are many researchers in this problem try to find the best solution to solve the problem with feasible computation time and to develop some models to solve real case problem. Heuristics are the most method used to solve this problem since the problem is an NP-hard. In this paper, we introduce a hybrid greedy algorithm and simulate annealing algorithm to solve a real container loading problem in one flexible packaging company in Indonesia. Validation is used to show that the method can be applied practically. We use seven real cases to check the validity and performance of the model. The proposed method outperformed the solution developed by the company in all seven cases with feasible computational time.


Author(s):  
Mahesh Ravishankar ◽  
Sandip Mazumder

The first-order spherical harmonics method (or P1 approximation) has found prolific usage for approximate solution of the radiative transfer equation (RTE) in participating media. However, the accuracy of the P1 approximation deteriorates as the optical thickness of the medium is decreased. The Modified Differential Approximation (MDA) was originally proposed to remove the shortcomings of the P1 approximation in optically thin situations. This article presents algorithms to apply the MDA to arbitrary geometry—in particular, three-dimensional (3D) geometry with obstructions, and inhomogeneous media. The wall-emitted component of the intensity was computed using a combined view-factor and ray-tracing approach. The Helmholtz equation, arising out of the medium-emitted component, was solved using an unstructured finite-volume procedure. The general procedure was validated against benchmark Monte Carlo results. The accuracy of MDA was found to be far superior to the standard P1 approximation in optically thin situations, and comparable to the P1 approximation in optically thick situations. Calculation and storage of the view-factor matrix was found to be a major shortcoming of the MDA, and several strategies to reduce both memory and computational time are discussed and demonstrated. In addition, for inhomogeneous media, calculation of optical distances requires a ray-tracing procedure, which was found to be a bottleneck from a computational efficiency standpoint.


Sign in / Sign up

Export Citation Format

Share Document