skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: Modeling Thermal Radiation in Combustion Environments: Progress and Challenges
Modeling thermal radiation in combustion environments can be extremely challenging for two main reasons. First, the radiative transfer equation (RTE), which is the cornerstone of modeling radiation in such environments, is a five-dimensional integro-differential equation. Second, the absorption and scattering coefficients of molecular gases and particulates prevalent in combustion environments oscillate strongly with the wavenumber (or wavelength), i.e., the medium is strongly nongray, requiring the solution of the RTE for a large number of wavenumbers. This article reviews the progress that has been made in this area to date with an emphasis on the work performed over the past three decades. Progress in both deterministic and stochastic (Monte Carlo) solutions of the RTE is reviewed, in addition to the review of the treatment of the spectral properties of gases, soot, and fuel droplets that dominate combustion environments, i.e., spectral or nongray models. The application of the various state-of-the-art nongray models and RTE solution methods to flames (particularly turbulent), fires, combustors, and other combustion systems are summarized along with a critical discussion of the pros and cons of the models and methods. Finally, the challenges that remain in modeling thermal radiation in combustion systems are highlighted and future outlooks are shared.  more » « less
Award ID(s):
2144290
PAR ID:
10579922
Author(s) / Creator(s):
;
Publisher / Repository:
MDPI
Date Published:
Journal Name:
Energies
Volume:
16
Issue:
10
ISSN:
1996-1073
Page Range / eLocation ID:
4250
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract Numerical modeling of radiative transfer in nongray reacting media is a challenging problem in computational science and engineering. The choice of radiation models is important for accurate and efficient high-fidelity combustion simulations. Different applications usually involve different degrees of complexity, so there is yet no consensus in the community. In this paper, the performance of different radiative transfer equation (RTE) solvers and spectral models for a turbulent piloted methane/air jet flame are studied. The flame is scaled from the Sandia Flame D with a Reynolds number of 22,400. Three classes of RTE solvers, namely the discrete ordinates method, spherical harmonics method, and Monte Carlo method, are examined. The spectral models include the Planck-mean model, the full-spectrum k-distribution (FSK) method, and the line-by-line (LBL) calculation. The performances of different radiation models in terms of accuracy and computational cost are benchmarked. The results have shown that both RTE solvers and spectral models are critical in the prediction of radiative heat source terms for this jet flame. The trade-offs between the accuracy, the computational cost, and the implementation difficulty are discussed in detail. The results can be used as a reference for radiation model selection in combustor simulations. 
    more » « less
  2. Abstract We prove existence, uniqueness and non-negativity of solutions of certain integral equations describing the density of states u ( z ) in the spectral theory of soliton gases for the one dimensional integrable focusing nonlinear Schrödinger equation (fNLS) and for the Korteweg–de Vries (KdV) equation. Our proofs are based on ideas and methods of potential theory. In particular, we show that the minimising (positive) measure for a certain energy functional is absolutely continuous and its density u ( z ) ⩾ 0 solves the required integral equation. In a similar fashion we show that v ( z ), the temporal analog of u ( z ), is the difference of densities of two absolutely continuous measures. Together, the integral equations for u , v represent nonlinear dispersion relation for the fNLS soliton gas. We also discuss smoothness and other properties of the obtained solutions. Finally, we obtain exact solutions of the above integral equations in the case of a KdV condensate and a bound state fNLS condensate. Our results is a step towards a mathematical foundation for the spectral theory of soliton and breather gases, which appeared in work of El and Tovbis (2020 Phys. Rev. E 101 052207). It is expected that the presented ideas and methods will be useful for studying similar classes of integral equation describing, for example, breather gases for the fNLS, as well as soliton gases of various integrable systems. 
    more » « less
  3. The rare-earth tritellurides (RTe 3 ) are a distinct class of 2D layered materials that recently gained significant attention due to hosting such quantum collective phenomena as superconductivity or charge density waves (CDWs). Many members of this van der Waals (vdW) family crystals exhibit CDW behavior at room temperature, i.e. , RTe 3 compound where R = La, Ce, Pr, Nd, Sm, Gd, and Tb. Here, our systematic studies establish the CDW properties of RTe 3 when the vdW spacing/interaction strength between adjacent RTe 3 layers is engineered under extreme hydrostatic pressures. Using a non-destructive spectroscopy technique, pressure-dependent Raman studies first establish the pressure coefficients of phonon and CDW amplitude modes for a variety of RTe 3 materials, including LaTe 3 , CeTe 3 , PrTe 3 , NdTe 3 , SmTe 3 , GdTe 3 , and TbTe 3 . Results further show that the CDW phase is eventually suppressed at high pressures when the interlayer spacing is reduced and interaction strength is increased. Comparison between different RTe 3 materials shows that LaTe 3 with the largest thermodynamic equilibrium interlayer spacing (smallest chemical pressure) exhibits the most stable CDW phases at high pressures. In contrast, CDW phases in late RTe 3 systems with the largest internal chemical pressures are suppressed easily with applied pressure. Overall results provide comprehensive insights into the CDW response of the entire RTe 3 series under extreme pressures, offering an understanding of CDW formation/engineering in a unique class of vdW RTe 3 material systems. 
    more » « less
  4. Abstract Numerical hydrodynamics simulations of gases dominated by ideal, nondegenerate matter pressure and thermal radiation pressure in equilibrium entail finding the temperature as part of the evolution. Since the temperature is not typically a variable that is evolved independently, it must be extracted from the evolved variables (e.g., the rest-mass density and specific internal energy). This extraction requires solving a quartic equation, which, in many applications, is done numerically using an iterative root-finding method. Here we show instead how the equation can be solved analytically and provide explicit expressions for the solution. We also derive Taylor expansions in limiting regimes and discuss the respective advantages and disadvantages of the iterative versus analytic approaches to solving the quartic. 
    more » « less
  5. For the inverse problem in physical models, one measures the solution and infers the model parameters using information from the collected data. Oftentimes, these data are inadequate and render the inverse problem ill-posed. We study the ill-posedness in the context of optical imaging, which is a medical imaging technique that uses light to probe (bio-)tissue structure. Depending on the intensity of the light, the forward problem can be described by different types of equations. High-energy light scatters very little, and one uses the radiative transfer equation (RTE) as the model; low-energy light scatters frequently, so the diffusion equation (DE) suffices to be a good approximation. A multiscale approximation links the hyperbolic-type RTE with the parabolic-type DE. The inverse problems for the two equations have a multiscale passage as well, so one expects that as the energy of the photons diminishes, the inverse problem changes from well- to ill-posed. We study this stability deterioration using the Bayesian inference. In particular, we use the Kullback–Leibler divergence between the prior distribution and the posterior distribution based on the RTE to prove that the information gain from the measurement vanishes as the energy of the photons decreases, so that the inverse problem is ill-posed in the diffusive regime. In the linearized setting, we also show that the mean square error of the posterior distribution increases as we approach the diffusive regime. 
    more » « less