发动机热喷流红外辐射计算与仿真

利用求每个小视场视线方向辐射亮度的方法计算喷流红外辐射的光谱分布。以辐射传递方程数值和形式为基础,采用Malkmus统计窄谱带模型和Curtis-Godson(CG)近似求视线方向的辐射强度。采用CFD分析软件FLUENT模拟流场和组分摩尔分数分布。建立喷流红外成像仿真模型,仿真生成了液体火箭发动机热喷流红外图像。结果表明,该方法可以很好地分辨出流场的细微结构。该模型也适用于航空发动机喷流红外辐射计算与仿真。     

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.     

碳氢扩散火焰实验研究：燃烧气氛对火焰光谱与温度的影响

辐射是各种燃烧过程中热传递的主要方式。在不同的火焰中,辐射光谱分布十分复杂。在这项工作中,利用光谱仪测量了可见光(200～900 nm),近红外(900～1 700 nm)和中红外(2 500～5 000 nm)波段火焰的光谱强度,分析了空气和富氧气氛下扩散火焰的光谱特征。并基于光谱分析,定量得到了火焰中碳烟以及气体发射的辐射力,计算了火焰的温度分布。结果表明,空气燃烧中的火焰温度低于富氧燃烧中的火焰温度。在空气气氛下,火焰中的碳烟和气体均对中的热辐射起着重要作用。而在富氧气氛下,气体对于火焰热辐射更为重要。在可见光和近红外波段,由于在空气气氛下火焰中碳烟的大量形成,光谱曲线显示出了良好连续性。而富氧气氛下火焰的辐射光谱降低。在中红外波段,空气气氛下火焰的气体辐射明显弱于富氧气氛下火焰的气体辐射。     

Coupled radiation and turbulent multiphase flow in an aluminised solid propellant rocket engine

We present in this paper numerical simulations of coupled radiative transfer and turbulent flows at high temperature and pressure, typical of multiphase flows encountered in aluminised solid propellant rocket engines. The radiating medium is constituted of gases and of liquid or solid particles of oxidised aluminum. The turbulent flow of the gaseous phase is treated by using a four equation, low Reynolds number, boundary-layer-type turbulence model. The distributions of concentrations, temperatures, and temperature fluctuation variances of particles are calculated from a Lagrangian approach and a turbulence dispersion model. Thermal and mechanical non-equilibrium between the gas and different classes of particles is allowed. A locally one dimensional, iteratively based, radiative transfer solver is developed to compute wall fluxes and radiative source terms. It is shown that the thermal boundary layer attenuates significantly the radiative fluxes coming from the outer regions. Particle radiation is found to be much more important than gas radiation. Turbulent dispersion of particles in the boundary layer induces a decrease of particle concentration in the region of maximum turbulent kinetic energy, and then, decreases the attenuation effect of wall fluxes due to the boundary layer. The effects of turbulent temperature fluctuations are found to be small in the problem under consideration.     

Matrix formulations of radiative transfer including the polarization effect in a coupled atmosphere-ocean system

A vector radiative transfer model has been developed for a coupled atmosphere-ocean system. The radiative transfer scheme is based on the discrete ordinate and matrix operator methods. The reflection/transmission matrices and source vectors are obtained for each atmospheric or oceanic layer through the discrete ordinate solution. The vertically inhomogeneous system is constructed using the matrix operator method, which combines the radiative interaction between the layers. This radiative transfer scheme is flexible for a vertically inhomogeneous system including the oceanic layers as well as the ocean surface. Compared with the benchmark results, the computational error attributable to the radiative transfer scheme has been less than 0.1% in the case of eight discrete ordinate directions. Furthermore, increasing the number of discrete ordinate directions has produced computations with higher accuracy. Based on our radiative transfer scheme, simulations of sun glint radiation have been presented for wavelengths of 670 nm and 1.6 μm. Results of simulations have shown reasonable characteristics of the sun glint radiation such as the strongly peaked, but slightly smoothed radiation by the rough ocean surface and depolarization through multiple scattering by the aerosol-loaded atmosphere. The radiative transfer scheme of this paper has been implemented to the numerical model named Pstar as one of the OpenCLASTR/STAR radiative transfer code systems, which are widely applied to many radiative transfer problems, including the polarization effect.     

Simultaneous solution of Kompaneets equation and radiative transfer equation in the photon energy range 1-125 keV

Radiative transfer equation in plane parallel geometry and Kompaneets equation is solved simultaneously to obtain theoretical spectrum of 1-125 keV photon energy range. Diffuse radiation field are calculated using time-independent radiative transfer equation in plane parallel geometry, which is developed using discrete space theory (DST) of radiative transfer in a homogeneous medium for different optical depths. We assumed free-free emission and absorption and emission due to electron gas to be operating in the medium. The three terms n, n² and (∂n/∂x_k) where n is photon phase density and x_k=(hν/kT_e), in Kompaneets equation and those due to free-free emission are utilized to calculate the change in the photon phase density in a hot electron gas. Two types of incident radiation are considered: (1) isotropic radiation with the modified black body radiation I^MB[1] and (2) anisotropic radiation which is angle dependent. The emergent radiation at τ=0 and reflected radiation τ=τ_max are calculated by using the diffuse radiation from the medium. The emergent and reflected radiation contain the free-free emission and emission from the hot electron gas. Kompaneets equation gives the changes in photon phase densities in different types of media. Although the initial spectrum is angle dependent, the Kompaneets equation gives a spectrum which is angle independent after several Compton scattering times.     

A numerical model for multigroup radiation hydrodynamics

We present in this paper a multigroup model for radiation hydrodynamics to account for variations of the gas opacity as a function of frequency. The entropy closure model (M₁) is applied to multigroup radiation transfer in a radiation hydrodynamics code. In difference from the previous grey model, we are able to reproduce the crucial effects of frequency-variable gas opacities, a situation omnipresent in physics and astrophysics. We also account for the energy exchange between neighbouring groups which is important in flows with strong velocity divergence. These terms were computed using a finite volume method in the frequency domain. The radiative transfer aspect of the method was first tested separately for global consistency (reversion to grey model) and against a well-established kinetic model through Marshak wave tests with frequency-dependent opacities. Very good agreement between the multigroup M₁ and kinetic models was observed in all tests. The successful coupling of the multigroup radiative transfer to the hydrodynamics was then confirmed through a second series of tests. Finally, the model was linked to a database of opacities for a Xe gas in order to simulate realistic multigroup radiative shocks in Xe. The differences with the previous grey models are discussed.     

气粒混合物非灰辐射特性合并宽窄谱带K分布模型

气粒混合物辐射问题具有全场性、非灰性、耦合性等特点,准确预估高温燃气/粒子非灰辐射特性是非常重要的。本文将合并宽窄谱带K分布模础(CWNBCK)与离散坐标法(DOM)结合,开展了非灰气粒混合物辐射换热问题的模拟工作,分别验证了一维和三维情况下应用该模型的准确性,给出不同工况下的热流源项、壁面热流或辐射热流等。结果表明:该模型能够给出与SNB模型精度基本相同的结果,考虑其计算效率的提高,可以在工程实际中应用该模型计算非灰气粒混合物辐射换热。     

灰气体加权和辐射模型综合评估及分析

在O_2/CO_2气氛下,参与性介质的非灰气体辐射特性表现出不同于空气气氛下的特性,因此,非灰气体辐射模型的选择和应用在换热过程中将变得十分重要.基于统计窄谱带模型,本文综合评估近年发展应用较广的灰气体加权和(WSGG)模型.结果表明,几种WSGG模型的预测值总体趋势正确,但仍存在着一定的差别.对于发射率,Dorigon等(2013 Int.J.Heat Mass Transfer 64 863)和Bordbar等(2014 Combust.Flame 161 2435)的WSGG模型与基准模型符合较好,相对误差小于20%.与离散坐标法结合,本文求解了PH_2O/PCO_2=1,2时的一维平行平板间辐射换热问题.结果显示,由Dorigon等和Bordbar等的WSGG模型得到的辐射热源和热流密度分布的相对误差均较小(10%左右).Johansson等(2011 Combust.Flame 158893)和Bordbar等的WSGG模型具有更广的适用范围.     

Application of adaptive multi-frequency grids to three-dimensional astrophysical radiative transfer

We present a method to solve the three-dimensional (3D) radiative transfer equation for astrophysical applications using adaptive photon transport grids. Contrary to earlier treatments, they are calculated for each frequency separately. Generated minimizing the first-order discretization error in the scattered radiation intensity, they provide global error control for solutions of radiative transfer problems on the grid. We discuss minimization of the grid point number in regions where the optical depth becomes large and show that the method allows for treating applications with optical depth of any value using the concept of penetration depth. The proposed grid generation algorithm is easy to implement, allows pre-calculation of the grids and storage in integer arrays, making a fast solution of the 3D radiative transfer equation possible. The grid generation algorithm is suitable for optimization in cases where simple radiation source distributions are given. Besides discussing application to simple density distribution commonly occurring in astrophysical objects, we illustrate the capabilities of the method by generating grids for an accretion disk around a young star.     

A concept of equivalent absorption and emission coefficients for gray analysis of radiative heat transfer

Due to the non-gray of gas radiation, the total emissivity differs from the total absorptivity. Therefore, in gray analysis, the equivalent absorption coefficient is not equal to the equivalent emission coefficient. Based on this idea, in this paper, a concept of equivalent absorption coefficient and equivalent emission coefficients is presented for gray analysis of gas radiative heat transfer. The equivalent emission coefficients are calculated by Leckner's formula and the equivalent absorption coefficients are estimated by inverse analysis. A one-dimensional gas radiation is taken as an example to show the efficiency of this concept. The results show that the concept of equivalent absorption and emission coefficients is feasible. It is necessary to use both the equivalent absorption coefficient and the equivalent emission coefficient for gray analysis.     

Radiative transfer in finite participating atmospheric aerosol media

The properties of radiation transfer through a plane-parallel atmospheric aerosol medium has been studied. It has been done by employing Mie theory to calculate the radiation transfer scattering parameters of the medium in the form of extinction, scattering, and absorption efficiencies. Then, the equation of radiative transfer through a plane-parallel atmosphere of aerosol has been solved for partial heat fluxes using two different analytical techniques, namely, the Variational Pomraning -Eddington approximation and Galerkin technique. Average efficiencies over log-normal and modified gamma size distributions are calculated. Therefore, the radiative properties of Carbon, Anthracite, Bituminous, Lignite, and Fly ash have been calculated. The obtained numerical results show very good agreement with each other in addition to the previous published work.     

Coupled radiation and natural convection: Different approaches of the slw model for a non-gray gas mixture

The coupling between non-gray radiation heat transfer and convection-conduction heat transfer is studied. The spectral line weighted sum of gray gases model (slw) is used to account for non-gray radiation properties. The aim of this work is to analyze the influence of the different approaches used when calculating the parameters of the slw model. Such strategies include the use of optimized model coefficients to reduce the number of operations, and the interpolation of the distribution function instead of the use of mathematical correlations. Non-gray calculations are also compared to gray solutions using the Planck mean absorption coefficient, which can be also calculated with the slw model. The radiative transfer equation (rte) is solved by means of the discrete ordinates method (dom). A natural convection driven cavity is chosen to couple radiation and conduction-convection energy transfer. Several cases, with a significant variation of the ratio between radiation to convection heat transfer, as well as the ratio between radiation to conduction heat transfer, are discussed.     

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.     

Transient coupled radiation-conduction in infinite semitransparent cylinders

A method is developed to analyze the transient coupled radiation-conduction in infinite semitransparent cylinders surrounded by isothermal black walls. The radiative heat source term is calculated by the radiative transfer coefficients and the transient energy equation is solved by an implicit finite difference method. The radiative transfer coefficients are deduced by use of the ray tracing method in combination with the Hottel and Sarofim zonal method. The effects of the related parameters on the transient radiative heat source and temperature distribution are analyzed. It is found that the peak of the dimensionless radial radiative heat source can be located at the interior shell of the cylinder with small optical thickness when heated by the surrounding irradiation. Treating the volume radiation as a surface radiation will result in large errors of transient temperature distribution for the cylinder with small optical thickness.     

Radiative heat transfer between two concentric spheres separated by a two-phase mixture of non-gray gas and particles using the modified discrete-ordinates method

The radiative heat transfer between two concentric spheres separated by a two-phase mixture of non-gray gas and a cloud of particles is investigated by using the combined finite-volume and discrete-ordinates method, named modified discrete-ordinates method (MDOM), which integrates the radiative transfer equation (RTE) over a control volume and a control angle simultaneously like in the finite-volume method (FVM) and treats the angular derivative terms due to spherical geometry as the conventional discrete-ordinates method (DOM). The radiative properties involving non-gray gas and particle behavior are modeled by using the extended weighted sum of gray gases model (WSGGM) with particles. Mathematical formulation and final discretization equations for the RTE are introduced by considering the behavior of a two-phase mixture of non-gray gas and particles in a spherically symmetric concentric enclosure. The present approach is validated by comparing with the results of previous works including gray and non-gray radiative heat transfer. Finally, a detailed investigation of the radiative heat transfer with non-gray gases and/or a two-phase mixture is conducted to examine the dependence of the radiative heat transfer upon temperature ratio between inner and outer spherical enclosure, particle concentration, and particle temperature.     

Coupling heat conduction and radiative transfer

A method to compute combined heat conduction and radiative transfer in a semitransparent medium is proposed. The heat equation is solved by an implicit finite-difference scheme and the radiation is simulated using a ray-tracing method. Stability and convergence for the coupling of both methods is shown. The crucial point for ray-tracing as well as for discrete-ordinate methods is the selection of the transfer directions and their associated weights. Using a model problem different direction sets proposed in the literature are compared with respect to their accuracy and the computational effort.     

Heat transfer by radiation and conduction in fibrous media without axial symmetry

Coupled radiative and conductive heat transfer in a fibrous medium formed by silica fibres is investigated in this paper by not taking account of the axial symmetry for the distribution of fibres or the boundary conditions. Radiative properties of the medium are calculated by using the Mie theory. The model obtained depends only on optical parameters (indices of silica) and on morphological parameters (diameter and orientation of the fibres, density of the medium). Simulations make it possible to study the strongly anisotropic behaviour of the scattering of the radiation by a fibre and to study the influence of various parameters on the radiative properties of the medium. The results of the Mie theory make possible the simulation of the heat transfer coupled by radiation and conduction. To do this, we introduce a new numerical scheme able to simulate heat transfer in the lack of axial symmetry. With this model, we can show the effects of distribution of fibres and temperature on the thermal behaviour of the medium as well as showing the importance of the phenomenon of scattering in fibrous media.     

Comparison of the standard weighted-sum-of-gray-gases with the absorption-line blackbody distribution function for the computation of radiative heat transfer in H2O/CO2 mixtures

This paper presents the computation of radiation heat transfer in a cylindrical enclosure in which the dimensions, the chemical species concentrations and the temperature fields make a realistic representation of an actual combustion chamber. Two gas models are applied and compared: the absorption-line blackbody distribution function (ALBDF), and the standard weighted-sum-of-gray-gases (WSGG) based on coefficients and correlations that are widely used in engineering. While the standard WSGG is restricted to the assumption of homogeneous gas mixture, the ALBDF can be applied to both homogeneous and non-homogeneous media. For the two gas models, the radiative exchanges are computed with the aid of the Monte Carlo method. The results show considerable discrepancies between the WSGG and the ALBDF models for the homogeneous medium. In addition, the importance of considering the non-homogeneity of the medium for an accurate computation of the radiative heat transfer is shown.     

A k-distribution-based radiation code and its computational optimization for an atmospheric general circulation model

The gas absorption process scheme in the broadband radiative transfer code “mstrn8”, which is used to calculate atmospheric radiative transfer efficiently in a general circulation model, is improved. Three major improvements are made. The first is an update of the database of line absorption parameters and the continuum absorption model. The second is a change to the definition of the selection rule for gas absorption used to choose which absorption bands to include. The last is an upgrade of the optimization method used to decrease the number of quadrature points used for numerical integration in the correlated k-distribution approach, thereby realizing higher computational efficiency without losing accuracy. The new radiation package termed “mstrnX” computes radiation fluxes and heating rates with errors less than 0.6 W/m² and 0.3 K/day, respectively, through the troposphere and the lower stratosphere for any standard AFGL atmospheres. A serious cold bias problem of an atmospheric general circulation model using the ancestor code “mstrn8” is almost solved by the upgrade to “mstrnX”.