A new approach to solve correlated k-distribution function

The resorting profile of absorption coefficient is fitted with Voigt function by employing the similarity between it and a single mixed broadened Voigt line type. In this way, simple mathematical fitting formulae are obtained to solve k-distribution (k-D) function at any pressure and temperature. Thus, a new correlated k-D method is proposed on the basis of it. Finally, taking mid-latitude summer atmosphere as an example, longwave cooling rates are calculated for three major gases by the new method, and then compared with the corresponding results by line-by-line integration.     

A full spectrum k-distribution based non-gray radiative property model for unburnt char

By using the concept of weighted sum of four gray particles and spectrum k-distribution (WSGP-SK), a non-gray radiative property model for unburnt char particles is developed. Based on the carbon burnout kinetic model for structure during oxidation, and the linear mixed approximation theory for complex index of refraction, spectral radiative properties of unburnt char particles are first calculated as function of the burnout ratio by Mie theory. Referring to the full spectrum k-distribution model, k-distribution is applied to reorder absorption and scattering efficiencies of particles. Then, weighting factors and efficiency factors of the non-gray radiative property model are directly obtained from Gaussian integral points of k-distribution. The model is validated against the benchmark solutions of line-by-line (LBL) model. Maximum relative errors of this model are 3% and 15% for radiative heat fluxes and source terms in non-isothermal inhomogeneous particulate media, respectively. The assumption of linearly varying radiative properties with burnout ratio (Lockwood et al. 1986) will result in a predicted deviation of 53% for radiative source terms. Results also show that this non-gray model is remarkably better than the Planck mean method. Moreover, a satisfactory comparison with LBL solutions is achieved in the gas and particle mixture by combining the non-gray WSGG-SK model (Guo et al. 2015). As a radiation sub-model, this non-gray radiative property model can significantly improve prediction accuracy of radiative heat transfer in oxy-fuel combustion.     

Generalization of the k-moment method using the maximum entropy principle. Application to the NBKM and full spectrum SLMB gas radiation models

The k-moment method is generalized by applying the maximum entropy principle to get several estimates of the k-distribution function on any kind of spectral interval as a function of the first two moments of the absorption coefficient. Corresponding formulations of the blackbody weighted band averaged transmission function of a gaseous uniform path are obtained. Different constraints involving the first and second order positive, first order negative and logarithmic moments are introduced together with a physical meaning whenever it is possible. Different sets of these constraints are considered to get maximum entropy estimates of the distributions functions: the Dirac, exponential, Gamma, inverse Gaussian and reciprocal inverse Gaussian k-distribution functions. Analytical formulas are provided for each of these distributions and for their associated transmission function, as a function of the mean and variance of the absorption coefficient. The methodology can be applied considering any spectral interval: narrow, wide, the full spectrum, continuous or not. Thus the resulting associated transmission and cumulative k-distribution functions can be utilized in the frame of a large variety of gas radiation models. Hence the k-moment method using the maximum entropy principle is assessed in the frame of the NBKM and full spectrum SLMB gas radiation models. A series of test cases implying comparisons with reference Line-by-Line results exhibits which maximum entropy k-distributions are likely to give the best estimations of narrow band or total emitted intensities, curves-of-growth of the total emission function and full spectrum cumulative k-distribution functions. In particular, the inverse Gaussian and Gamma k-distributions seem most of the time to perform very well.     

The k-moment method for the narrow band modeling of radiative properties of nonuniform gaseous media

Moments of the absorption coefficient distribution function are used for the derivation of statistical narrow band (SNB) model parameters of nonuniform optical paths in gases. The approach yields approximations for the path-averaged first- and second-order k-moments from which equivalent SNB parameters are determined in the frame of the Malkmus model. The approach is assessed through comparisons with LBL data. The nonuniform approximation is shown to enable the computation of transmissivities and radiation intensities with accuracy similar to or higher than those achieved by the Curtis-Godson one.     

A new k-interval selection technique for fast atmospheric radiance calculation in remote sensing applications

A new technique is proposed to generate the k-interval parameters, including the number of k-intervals, the equivalent absorption coefficients, and the quadrature weights when using the correlated k-distribution method for the computation of spectrally integrated three-dimensional (3D) atmospheric radiance. The main difference between the proposed technique and the traditional exponential sum fitting technique is that only quadrature weights are involved in the optimization process. To avoid the ill-conditioned problem in the proposed technique, the absorption coefficients with high value are dealt with by the delta log(k) (Δlog(k)) technique instead of involving them in the fitting procedure. The performance of the proposed technique is illustrated by radiance calculation results of cloudless and cloudy atmosphere for three different band settings. Results show that there are less relative errors with the proposed optimization technique than with the Δlog(k) technique under the same number of k-intervals. However, as the absorption becomes stronger, the performance of the proposed technique gradually decreases to the Δlog(k) technique. The relative root-mean-square error (RMSE) of radiance for 3D cloudy atmosphere can be controlled in less than 2% when the number of k-intervals is less than 10 particularly for weak absorption band, the RMSEs are less than 1% with only 6 terms.     

A nonuniform narrow band correlated-k approximation using the k-moment method

Narrow band k-moment (NBKM) formulations of the transmission function for nonuniform gaseous paths are developed. One of them, called the scaled variance (SV) approximation, is based on the assumption of correlated absorption coefficients. Theoretical derivations are properly detailed and several test cases taken from the literature are provided to assess the approximate modeling results against LBL reference calculations. In most cases representative of combustion configurations, it is shown that the scaled variance approximation is more accurate than the Curtis–Godson one for CO₂ but not for H₂O.     

Extension of the weak-line approximation and application to correlated-k methods

Global climate models require accurate and rapid computation of the radiative transfer through the atmosphere. Correlated-k methods are often used. One of the approximations used in correlated-k models is the weak-line approximation. We introduce an approximation T_γ which reduces to the weak-line limit when optical depths are small, and captures the deviation from the weak-line limit as the extinction deviates from the weak-line limit. This approximation is constructed to match the first two moments of the gamma distribution to the k-distribution of the transmission. We compare the errors of the weak-line approximation with T_γ in the context of a water vapor spectrum. The extension T_γ is more accurate and converges more rapidly than the weak-line approximation.     

The effects of the choice of the k-interval number on radiative calculations

The selection of the number of k-interval is a foundation to correlated k-distribution method and the problem of how to do it still remains unsettled. It is pointed out by numerical computation in this work that choosing the number of k-interval is a major factor affecting accuracy and speed in radiative calculation. To increase the number of k-interval is an efficient method to improve the accuracy. However, it is found by this study that there exists a saturation of the accuracy to an increase of the number. The optimal rules on the number of k-interval choosing are proposed in the paper. Then, five versions on atmospheric absorption by gases appropriate for GCMs are given according to them.     

A simplified CFD model for spectral radiative heat transfer in high-pressure hydrocarbon–air combustion systems

Detailed radiation modeling in piston engines has received relatively little attention to date. Recently, it is being revisited in light of current trends towards higher operating pressures and higher levels of exhaust-gas recirculation (EGR), both of which enhance molecular gas radiation. Advanced high-efficiency engines also are expected to function closer to the limits of stable operation, where even small perturbations to the energy balance can have a large influence on system behavior. Detailed radiation modeling using sophisticated tools like photon Monte Carlo/line-by-line (PMC/LBL) is computationally expensive. Here, guided by results from PMC/LBL, a simplified stepwise-gray spectral model in combination with a first-order spherical harmonics (P1 method) radiative transfer equation (RTE) solver is proposed and tested for engine-relevant conditions. Radiative emission, reabsorption and radiation reaching the walls are computed for a heavy-duty compression-ignition engine at part-load and full-load operating conditions with different levels of EGR and soot. The results are compared with those from PMC/LBL, P1/FSK (P1 with a full-spectrum k-distribution spectral model) and P1/Gray radiation models to assess the proposed model’s accuracy and computational cost. The results show that the proposed P1/StepwiseGray model can calculate reabsorption locally and globally with less than 10% error (with respect to PMC/LBL) at a small fraction of the computational cost of PMC/LBL (a factor of 30) and P1/FSK (a factor of 15). In contrast, error in computed reabsorption by the P1/Gray model is as high as 60%. It is expected that the simplified model should be broadly applicable to high-pressure hydrocarbon–air combustion systems, with or without soot.     

Evaluation of solution methods for radiative heat transfer in gaseous oxy-fuel combustion environments

The exact solution to radiative heat transfer in combusting flows is not possible analytically due to the complex nature of the integro-differential radiative transfer equation (RTE). Many different approximate solution methods for the solution of the RTE in multi-dimensional problems are available. In this paper, two of the principal methods, the spherical harmonics (P₁) and the discrete ordinates method (DOM) are used to calculate radiation. The radiative properties of the gases are calculated using a non-gray gas full spectrum k-distribution method and a gray method. Analysis of the effects of numerical quadrature in the DOM and its effect on computation time is performed. Results of different radiative property methods are compared with benchmark statistical narrow band (SNB) data for both cases that simulate air combustion and oxy-fuel combustion. For both cases, results of the non-gray full spectrum k-distribution method are in good agreement with the SNB data. In the case of oxy-fuel simulations with high partial pressures of carbon dioxide, use of gray method for the radiative properties may cause errors and should be avoided.     

Numerical explanation for accurate radiative cooling rates resulting from the correlated k-distribution hypothesis

Numerical simulations are conducted to demonstrate that (a) correlation between spectral absorption at two levels decreases with increasing level distance; and (b) the contribution to the flux gained by one layer from another decreases rapidly with the layer separation. The combination of the two facts explains why the existence of poor correlation between distant layers produces insignificant radiative cooling rate error. Therefore, the overall accuracy of cooling rate calculations under the correlated k-distribution hypothesis is high.     

The relationship between absorption coefficient and temperature and their effect on the atmospheric cooling rate

Several approaches are considered to determine the temperature effect on the absorption coefficient within a correlated k-distribution method. Taking in the 610- region for example, the absorption coefficients and atmospheric cooling rates calculated using these approaches are compared with line-by-line integration. It is emphasized in this paper by numerical calculation that the effect of pressure on absorption coefficient is related to temperature and vise versa; the larger the pressure, the larger the effect of temperature on absorption coefficient. Results show that the temperature effect must be considered in radiative calculations although its effect on the absorption coefficient is much smaller than that of pressure.     

Non-gray chemical composition based radiative property model of fly ash particles

Particle radiation has a spectral dependence and is closely related to the chemical composition of the material. Iron oxide, one of the main components of fly ash, observably affects the complex index of refraction of the particles. In this study, following the theory of the spectrum k-distribution based weighted sum of gray particles model (Guo et al. [4,13]), a non-gray fly ash radiative property model involving the chemical composition was developed. First, four typical fly ash particles with different iron oxide contents were selected, and the corresponding particle radiative parameters were obtained using the Mie theory. Then, the absorption efficiency and weighting factors of the non-gray model were directly obtained from the Gaussian integral points of the k-distribution. The scattering efficiency of the particles was obtained from the Planck mean. The accuracy of the newly developed model was evaluated in a one-dimensional plane-parallel slab system through comparison with the line-by-line (LBL) model and two commonly used gray radiative property models. The results show that the new non-gray model agrees well with the LBL solution and becomes more accurate as the iron oxide content increases. When the iron oxide content of the fly ash increased from 5.47% to 30.50%, the maximum relative error of the radiative heat flux and the radiative source term decreased from 12.50% to 5.68% and from 20.97% to 12.62%, respectively. The new model can improve the prediction accuracy of radiative heat transfer in pulverized coal-fired furnaces.     

A database for the SLMB modeling of the full spectrum radiative properties of CO2

Development and application of a database for the Spectral-Line Moment-Based (SLMB) modeling of the full spectrum radiative properties of mixtures of carbon dioxide and nitrogen is presented. The critical issue of the definition of a reference thermophysical condition is addressed together with the suggestion of a coherent and precise methodology to derive parameters of the model for any other configuration. The database is built accordingly from the CDSD-1000 high temperature spectroscopic databank for gas and blackbody-weighting temperatures in the range [300; 2700 K]. Accuracy of both the modeling and the database is assessed through comparisons with LBL results in terms of full spectrum k-distributions and emission functions. Results obtained from the application of FSK correlations and the Leckner's formula are also provided for extended analyses.     

Application of approximations for joint cumulative k-distributions for mixtures to FSK radiation heat transfer in multi-component high temperature non-LTE plasmas

Approximations for joint cumulative k-distribution for mixtures are efficient for full spectrum k-distribution (FSK) computations. These approximations provide reduction of the database that is necessary to perform FSK computation when compared to the direct approach, which uses cumulative k-distributions computed from the spectrum of the mixture, and also less computational expensive when compared to techniques in which RTE's are required to be solved for each component of the mixture. The aim of the present paper is to extend the approximations for joint cumulative k-distributions for non-LTE media. For doing that, a FSK to non-LTE media formulation well-suited to be applied along with approximations for joint cumulative k-distributions is presented. The application of the proposed methodology is demonstrated by solving the radiation heat transfer in non-LTE high temperature plasmas composed of N, O, N₂, NO, N₂⁺ and mixtures of these species. The two more efficient approximations, that is, the superposition and multiplication are employed and analyzed.     

Multi-scale k-distribution model for gas mixtures in hypersonic nonequilibrium flows

A k-distribution model is presented for gas mixtures in thermodynamic nonequilibrium, containing strongly radiating atomic species N and O together with molecular species of N₂, N₂⁺, NO and O₂. In the VUV range of the spectrum there is strong absorption of atomic radiation by bands of N₂. For this spectral range, a multi-scale model is presented, where RTEs are solved separately for each emitting species and overlap with other species is treated in an approximate way. Methodology for splitting the gas mixture into scales and evaluation of the overlap parameter between different scales is presented. The accuracy of the new model is demonstrated by solving the radiative transfer equation along the stagnation line flow field of the Crew Exploration Vehicle (CEV).     

A study of the radiative forcing and global warming potentials of hydrofluorocarbons

We developed a new radiation parameterization of hydrofluorocarbons (HFCs), using the correlated k-distribution method and the high-resolution transmission molecular absorption (HITRAN) 2004 database. We examined the instantaneous and stratospheric adjusted radiative efficiencies of HFCs for clear-sky and all-sky conditions. We also calculated the radiative forcing of HFCs from preindustrial times to the present and for future scenarios given by the Intergovernmental Panel on Climate Change Special Report on Emission Scenarios (SRES, in short). Global warming potential and global temperature potential were then examined and compared on the basis of the calculated radiative efficiencies. Finally, we discuss surface temperature changes due to various HFC emissions.     

High fidelity radiative heat transfer models for high-pressure laminar hydrogen–air diffusion flames

Radiative heat transfer is studied numerically for high-pressure laminar H₂–air jet diffusion flames, with pressure ranging from 1 to 30 bar. Water vapour is assumed to be the only radiatively participating species. Two different radiation models are employed, the first being the full spectrum k-distribution model together with conventional Radiative Transfer Equation (RTE) solvers. Narrowband k-distributions of water vapour are calculated and databased from the HITEMP 2010 database, which claims to retain accuracy up to 4000 K. The full-spectrum k-distributions are assembled from their narrowband counterparts to yield high accuracy with little additional computational cost. The RTE is solved using various spherical harmonics methods, such as P₁, simplified P₃ (SP₃) and simplified P₅ (SP₅). The resulting partial differential equations as well as other transport equations in the laminar diffusion flames are discretized with the finite-volume method in OpenFOAM^®. The second radiation model is a Photon Monte Carlo (PMC) method coupled with a line-by-line spectral model. The PMC absorption coefficient database is derived from the same spectroscopy database as the k-distribution methods. A time blending scheme is used to reduce PMC calculations at each time step. Differential diffusion effects, which are important in laminar hydrogen flames, are also included in the scalar transport equations. It was found that the optically thin approximation overpredicts radiative heat loss at elevated pressures. Peak flame temperature is less affected by radiation because of faster chemical reactions at high pressures. Significant cooling effects are observed at downstream locations. As pressure increases, the performance of RTE models starts to deviate due to increased optical thickness. SP_N models perform only marginally better than P₁ because P₁ is adequate except at very high pressure.     

Theory of the hot free carrier intraband-absorption coefficient of n-inversion layers

The absorption coefficient K of a quasi two dimensional (2D) hot free electron gas is calculated for the first time as a function of the lattice temperature T, the photon angular frequency w, the carrier density N_s as well as the electron temperature T_e when the carriers are scattered by ionized impurities, acoustic phonons and polar optical phonons. Analytical expressions are derived in the limiting cases of non-degeneracy and degeneracy of the electron system. In the quantum limit ħw/k_BT_e ≳ 1 where the interaction between the electron and the photon is inelastic K sensitively depends on the limiting scattering mechanism showing that the electron motion is completely controlled by the photon field. In the classical limit ħw/k_BT_e ⪡ 1 the absorption decreases proportional to w¹ independent of the limiting scattering mechanism in agreement with the experimental data deduced from far-infrared absorptivity measurements on GaAs heterolayers.     

Numerical study on surface acoustic wave method for determining Young's modulus of low-k films involved in multi-layered structures

The surface acoustic waves (SAWs) technique is becoming an attractive tool for accurately and nondestructively characterizing the mechanical property of the brittle low dielectric constant (low-k) thin film. The theoretical equations for describing SAWs propagating on the multi-layered structure are derived in this study. The dispersion features of SAWs propagating on different structures of low-k/SiO₂/Si substrate, SiO₂/low-k/Si substrate, low-k/Si substrate, and low-k/Cu/Si substrate are investigated to instruct an accurate and facile fitting process for determining Young's modulus of low-k films. The dependence of dispersion relation on the film thickness, elastic modulus of low-k materials as well as frequency are provided and discussed in detail. The study shows an obvious influence of layered structure on the dispersion relation of SAWs. For a fixed structure, the dispersion curvature increases with the decrease of Young's modulus of low-k films.