scholarly journals Accelerated Solution of Discrete Ordinates Approximation to the Boltzmann Transport Equation for a Gray Absorbing–Emitting Medium Via Model Reduction

2017 ◽  
Vol 139 (12) ◽  
Author(s):  
John Tencer ◽  
Kevin Carlberg ◽  
Marvin Larsen ◽  
Roy Hogan

This work applies a projection-based model-reduction approach to make high-order quadrature (HOQ) computationally feasible for the discrete ordinates approximation of the radiative transfer equation (RTE) for purely absorbing applications. In contrast to traditional discrete ordinates variants, the proposed method provides easily evaluated error estimates associated with the angular discretization as well as an efficient approach for reducing this error to an arbitrary level. In particular, the proposed approach constructs a reduced basis from (high-fidelity) solutions of the radiative intensity computed at a relatively small number of ordinate directions. Then, the method computes inexpensive approximations of the radiative intensity at the (remaining) quadrature points of a high-order quadrature using a reduced-order model (ROM) constructed from this reduced basis. This strategy results in a much more accurate solution than might have been achieved using only the ordinate directions used to construct the reduced basis. One- and three-dimensional test problems highlight the efficiency of the proposed method.

Author(s):  
John Tencer ◽  
Kevin Carlberg ◽  
Roy Hogan ◽  
Marvin Larsen

Radiation heat transfer is an important phenomenon in many physical systems of practical interest. When participating media is important, the radiative transfer equation (RTE) must be solved for the radiative intensity as a function of location, time, direction, and wavelength. In many heat transfer applications, a quasi-steady assumption is valid. The dependence on wavelength is often treated through a weighted sum of gray gases type approach. The discrete ordinates method is the most common method for approximating the angular dependence. In the discrete ordinates method, the intensity is solved exactly for a finite number of discrete directions, and integrals over the angular space are accomplished through a quadrature rule. In this work, a projection-based model reduction approach is applied to the discrete ordinates method. A small number or ordinate directions are used to construct the reduced basis. The reduced model is then queried at the quadrature points for a high order quadrature in order to inexpensively approximate this highly accurate solution. This results in a much more accurate solution than can be achieved by the low-order quadrature alone. One-, two-, and three-dimensional test problems are presented.


2021 ◽  
Vol 89 (3) ◽  
Author(s):  
Panagiotis Tsoutsanis ◽  
Ebenezer Mayowa Adebayo ◽  
Adrian Carriba Merino ◽  
Agustin Perez Arjona ◽  
Martin Skote

AbstractIn this paper we extend the application of unstructured high-order finite-volume central-weighted essentially non-oscillatory (CWENO) schemes to multicomponent flows using the interface capturing paradigm. The developed method achieves high-order accurate solution in smooth regions, while providing oscillation free solutions at discontinuous regions. The schemes are inherently compact in the sense that the central stencils employed are as compact as possible, and that the directional stencils are reduced in size, therefore simplifying their implementation. Several parameters that influence the performance of the schemes are investigated, such as reconstruction variables and their reconstruction order. The performance of the schemes is assessed under a series of stringent test problems consisting of various combinations of gases and liquids, and compared against analytical solutions, computational and experimental results available in the literature. The results obtained demonstrate the robustness of the new schemes for several applications, as well as their limitations within the present interface-capturing implementation.


1992 ◽  
Vol 20 (1) ◽  
pp. 33-56 ◽  
Author(s):  
L. O. Faria ◽  
J. T. Oden ◽  
B. Yavari ◽  
W. W. Tworzydlo ◽  
J. M. Bass ◽  
...  

Abstract Recent advances in the development of a general three-dimensional finite element methodology for modeling large deformation steady state behavior of tire structures is presented. The new developments outlined here include the extension of the material modeling capabilities to include viscoelastic materials and a generalization of the formulation of the rolling contact problem to include special nonlinear constraints. These constraints include normal contact load, applied torque, and constant pressure-volume. Several new test problems and examples of tire analysis are presented.


Energies ◽  
2021 ◽  
Vol 14 (3) ◽  
pp. 674
Author(s):  
Paul Caicedo ◽  
David Wood ◽  
Craig Johansen

Solar chimney power plants (SCPPs) collect air heated over a large area on the ground and exhaust it through a turbine or turbines located near the base of a tall chimney to produce renewable electricity. SCPP design in practice is likely to be specific to the site and of variable size, both of which require a purpose-built turbine. If SCPP turbines cannot be mass produced, unlike wind turbines, for example, they should be as cheap as possible to manufacture as their design changes. It is argued that a radial inflow turbine with blades made from metal sheets, or similar material, is likely to achieve this objective. This turbine type has not previously been considered for SCPPs. This article presents the design of a radial turbine to be placed hypothetically at the bottom of the Manzanares SCPP, the only large prototype to be built. Three-dimensional computational fluid dynamics (CFD) simulations were used to assess the turbine’s performance when installed in the SCPP. Multiple reference frames with the renormalization group k-ε turbulence model, and a discrete ordinates non-gray radiation model were used in the CFD simulations. Three radial turbines were designed and simulated. The largest power output was 77.7 kW at a shaft speed of 15 rpm for a solar radiation of 850 W/m2 which exceeds by more than 40 kW the original axial turbine used in Manzanares. Further, the efficiency of this turbine matches the highest efficiency of competing turbine designs in the literature.


Author(s):  
J. R. Beisheim ◽  
G. B. Sinclair ◽  
P. J. Roache

Current computational capabilities facilitate the application of finite element analysis (FEA) to three-dimensional geometries to determine peak stresses. The three-dimensional stress concentrations so quantified are useful in practice provided the discretization error attending their determination with finite elements has been sufficiently controlled. Here, we provide some convergence checks and companion a posteriori error estimates that can be used to verify such three-dimensional FEA, and thus enable engineers to control discretization errors. These checks are designed to promote conservative error estimation. They are applied to twelve three-dimensional test problems that have exact solutions for their peak stresses. Error levels in the FEA of these peak stresses are classified in accordance with: 1–5%, satisfactory; 1/5–1%, good; and <1/5%, excellent. The present convergence checks result in 111 error assessments for the test problems. For these 111, errors are assessed as being at the same level as true exact errors on 99 occasions, one level worse for the other 12. Hence, stress error estimation that is largely reasonably accurate (89%), and otherwise modestly conservative (11%).


Author(s):  
Athanasios Donas ◽  
Ioannis Famelis ◽  
Peter C Chu ◽  
George Galanis

The aim of this paper is to present an application of high-order numerical analysis methods to a simulation system that models the movement of a cylindrical-shaped object (mine, projectile, etc.) in a marine environment and in general in fluids with important applications in Naval operations. More specifically, an alternative methodology is proposed for the dynamics of the Navy’s three-dimensional mine impact burial prediction model, Impact35/vortex, based on the Dormand–Prince Runge–Kutta fifth-order and the singly diagonally implicit Runge–Kutta fifth-order methods. The main aim is to improve the time efficiency of the system, while keeping the deviation levels of the final results, derived from the standard and the proposed methodology, low.


2021 ◽  
Vol 40 (3) ◽  
Author(s):  
Bo Hou ◽  
Yongbin Ge

AbstractIn this paper, by using the local one-dimensional (LOD) method, Taylor series expansion and correction for the third derivatives in the truncation error remainder, two high-order compact LOD schemes are established for solving the two- and three- dimensional advection equations, respectively. They have the fourth-order accuracy in both time and space. By the von Neumann analysis method, it shows that the two schemes are unconditionally stable. Besides, the consistency and convergence of them are also proved. Finally, numerical experiments are given to confirm the accuracy and efficiency of the present schemes.


Sign in / Sign up

Export Citation Format

Share Document