Temperature field of radiative equilibrium in a two-dimensional graded index medium with gray boundaries

In this paper, a numerical method is presented for the study of the radiative transfer in a two-dimensional graded index semitransparent medium with diffuse gray boundaries. The numerical method is a combination of the linear refractive index bar model, the discrete curved ray-tracing technique and the pseudo source adding method (LRIB-CRTP). In the traditional ray-tracing technique, it is difficult to deal with the diffuse gray boundary while solving the radiative transfer. Using the pseudo source adding method, the diffuse gray boundary of the medium can be treated as a black boundary. We have also studied the radiative equilibrium temperature field of the medium and analyzed the influence of some parameters involved. The results show that the directional discrete number is important for the medium having a large absorption coefficient. The results also show that the refractive index distribution greatly influences the temperature field, whereas the linear absorption coefficient distribution has little influence on the temperature field.     

Calculation of radiative properties of nonequilibrium hydrogen plasma

A computer program called NEQRAP is described that calculates the radiative properties of nonequilibrium ionized hydrogen. From the given electron temperature, electron density, and atom density values (which do not necessarily satisfy the equilibrium relationship) and intensities of incident radiation, the non-Boltzmann populations of electronic states are computed by solving the equation of quasi-steady- state population distribution. Emission and absorption coefficients are determined as functions of wavelength by invoking the principle of detailed balance between the upper and lower states of each radiative transition. Radiative transport through the medium is computed assuming a one-dimensional uniform slab. The rate of ionic reaction is also computed. When used on a sample case, the program shows that there is a large difference between the calculated intensities of radiation emitted by a bulk of equilibrium and nonequilibrium hydrogen. Accuracy of the program is estimated to be better than 10%.     

The effects of gradients on the diagnostic use of spectral features from laser compressed plasmas

The use of spectral features as diagnostics of the plasma temperature and density is studied using a non-LTE radiative transfer model, which solves the coupled statistical equilibrium and radiative transfer equation consistently. The effect of gradients on the hydrodynamic data, i.e. temperature and mass density, is studied to establish the accuracy of the diagnostic procedures.     

Spectral emission from MHD combustion gas mixed with pulverized coal

The spectral emission from MHD combustion gas mixed with pulverized coal was measured between 0.38 and 2.7 μm. Twenty lines of K, eight lines of Li, Na, Rb, and Ca, and two molecular spectra of calcium compounds were observed. The experimental line shape for the 4P-4S transition of K agreed well with the Lorentzian shape combined with the measured farwing absorption cross sections. The population temperature of the K electronic state in the optically thin limit agreed well with the wing-reversal temperature. Continuous emission from coal-ash particles was not observed.     

Transient coupled radiative and conductive heat transfer in a semitransparent layer of ceramic

This work considers transient radiative and conductive heat transfer in a semitransparent layer of ceramic, submitted to several thermal and radiative boundary conditions. Each side of the layer is exposed to hot or cold radiative surroundings, while each boundary is heated or cooled by convection. The solution procedure must provide accurate temperature distribution in the layer, so a nodal analysis based on Hottel's zonal method extended by ray tracing method is carried out. A finite difference method with non-uniform space and time increments is used to solve the transient energy equation, including a radiative heat source, coupled to a equation of radiative transfer. Variable spacing was used to concentrate grid points in regions with large temperature gradients. The influence of refractive index, optical thicknesses and conduction-radiation parameters is investigated.     

Monte Carlo ray-tracing simulation for radiative heat transfer in turbulent fluctuating media under the optically thin fluctuation approximation

The Monte Carlo ray-tracing method (MCRT) based on the concept of radiation distribution factor is extended to solve radiative heat transfer problem in turbulent fluctuating media under the optically thin fluctuation approximation. A one-dimensional non-scattering turbulent fluctuating media is considered, in which the mean temperature and absorption coefficient distribution are assumed and the shape of probability density function is given. The distribution of the time-averaged volume radiation heat source is solved by MCRT and direct integration method. It is shown that the results of MCRT based on the concept of radiation distribution factor agree with these of integration solution very well, but results of MCRT based on the concept of radiative transfer coefficient do not agree with these of integration solution. The solution of time-averaged radiative transfer equation by the concept of radiative transfer coefficient should be treated with caution.     

Influence of the boundary layer of a shock-heated plasma on the stark profiles of the A1(I) resonance lines

The influence of temperature and density gradients in the boundary layer of an Al(CH₃)₃-seeded shock-heated plasma on the Stark profiles of the A1(I) resonance lines λ = 3961.5 Å and λ = 3944 Å is investigated by solving the radiative transfer equation. The thermal boundary layer thickness δ is varied from 0 to 1 cm, assuming complete local thermodynamic equilibrium (LTE) throughout the inhomogeneous region. The calculations show that the line profiles are distorted. For the plasma conditions of an experiment made recently by one of us, it is shown that reabsorption diminishes the intensity and displaces the maxima of the profiles of both resonance lines in the same relative proportions. Thus, the ratios of their half-widths and of their shifts are practically independent of the boundary-layer effect. We conclude, therefore, that the measured difference of the half-widths of the A1(I) resonance doublet must be an intrinsic plasma effect.     

管内高温介质层流入口段中的热辐射作用

数值研究了高温介质密度随温度变化时,管内层流入口段耦合换热中的热辐射作用。采用离散坐标法、控制容积法耦合求解辐射传递方程、能量方程及N-S方程。考察了中等大小光学厚度下,热辐射作用对介质内速度分布、温度分布以及换热的影响。结果表明,即使在不大的光学厚度下,热辐射作用对管内高温介质层流入口段耦合换热的速度场与换热强度都有明显影响。     

漫反射各向异性散射介质内的温度热流场

本文采用射线踪迹结合节点分析法和谱带模型,研究了漫反射不透明边界下吸收、发射、各向异性散射介质内的热辐射传递过程。考虑介质辐射能的入射和散射方向,导出漫反射、不透明边界、各向异性散射介质的辐射传递系数。在辐射平衡的情况下,考察了表面发射率和散射反照率对介质内辐射热流和温度场的影响。研究表明,介质不透明边界处存在温度跃迁现象,而且,内界面发射率越大,相应界面温度跃迁越小。     

EAST位形磁流体平衡的模拟

对托卡马克平衡反演数值计算代码EFIT进行了改进,使之适用于EAST装置的磁流体平衡研究.用改进后的平衡反演代码EFIT模拟了EAST装置稳态阶段的偏滤器平衡位形.提出极向场线圈电流作为线性方程组主约束条件,期望位形作为从约束条件,求解非对称电流条件下自由边界平衡问题的方法.结果表明,该方法能较为容易地得到非对称电流条件产生的平衡位形,在非对称电流条件下自由边界平衡计算是收敛的.     

An immersed boundary-thermal lattice Boltzmann method using an equilibrium internal energy density approach for the simulation of flows with heat transfer

In present paper, a novel immersed boundary-thermal lattice Boltzmann method by the name of “an equilibrium internal energy density approach” is proposed to simulate the flows around bluff bodies with the heat transfer. The main idea is to combine the immersed boundary method (IBM) with the thermal lattice Boltzmann method (TLBM) based on the double population approach. The equilibrium internal energy density approach based on the equilibrium velocity approach [X. Shan, H. Chen, Lattice Boltzmann model for simulating flows with multiple phases and components, Phys. Rev. E 47 (1993) 1815] is used to combine IBM with TLBM. The idea of the equilibrium internal energy density approach is that the satisfaction of the energy balance between heat source on the immersed boundary point and the amount of change of the internal energy density according to time ensures the temperature boundary condition on the immersed boundary. The advantages of this approach are the simple concept, easy implementation and the utilization of original governing equation without modification. The simulation of natural convection in a square cavity with various body shapes for different Rayleigh numbers has been conducted to validate the capability and the accuracy of present method on solving heat transfer problems. Consequently, the present results are found to be in good agreement with those of previous studies.     

Theoretical and Experimental Investigations of the Suitability of Low‐Current Z‐Pinch Plasma as an Absorption Medium for Laser Radiation

In this study the electron density of z‐pinch plasmas driven at relatively low currents (ca. 2‐5 kA) was determined using only emission spectroscopy. The suitability of a hollow‐cathode‐triggered z‐pinch plasma as an absorption medium for laser radiation was investigated. The temporal and spatial behaviors of electron temperature and density profiles were estimated using magnetohydrodynamic (MHD) simulations to evaluate the experimental results. Temperature measurements were performed according to the Boltzmann plot method in the visible spectrum range, using the fact that, in low‐current z‐pinch plasma, a local thermodynamic equilibrium prevails for states at high principal quantum numbers (partial local thermodynamic equilibrium). In this case, the Saha equation can be used to determine the electron density. The results demonstrate that this method of determining the electron temperature and density of z‐pinch plasmas is only applicable during the pinching phase. However, in this case the experimentally determined values are in fairly good agreement with the values determined using the MHD model. A user‐oriented 1‐D radiation MHD code was used to simulate the dynamic evolution of the plasma. The experimentally determined maximum electron temperature of approximately 12 eV is in fairly good agreement with the simulated value. (© 2008 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Two analytic methods applied to radiative transfer in the linear-anisotropic scattering medium with graded index and diffuse gray boundaries

The curved ray-tracing method is extended to radiative transfer in the graded index medium with diffuse gray boundary conditions instead of black boundary conditions and the pseudo-source adding method is extended to the case of the linear-anisotropic scattering medium with graded index from non-scattering medium. Furthermore, the equivalence of the two methods is verified by formulation derivation. As exact analytical solutions, both the methods have high accuracy and fast computational speed. The predicted temperature distributions and dimensionless radiative heat flux at radiative equilibrium are determined by the proposed methods, and the numerical results are compared with the data in references. The results show that the present methods have a good accuracy. Influences of various combinations of refractive index and boundary emissivities on the temperature distributions and dimensionless radiative heat flux are also investigated.     

A current density conservative scheme for incompressible MHD flows at a low magnetic Reynolds number. Part I: On a rectangular collocated grid system

A consistent, conservative and accurate scheme has been designed to calculate the current density and the Lorentz force by solving the electrical potential equation for magnetohydrodynamics (MHD) at low magnetic Reynolds numbers and high Hartmann numbers on a finite-volume structured collocated grid. In this collocated grid, velocity (u), pressure (p), and electrical potential (φ) are located in the grid center, while current fluxes are located on the cell faces. The calculation of current fluxes on the cell faces is conducted using a conservative scheme, which is consistent with the discretization scheme for the solution of electrical potential Poisson equation. A conservative interpolation is used to get the current density at the cell center, which is used to conduct the calculation of Lorentz force at the cell center for momentum equations. We will show that both “conservative” and “consistent” are important properties of the scheme to get an accurate result for high Hartmann number MHD flows with a strongly non-uniform mesh employed to resolve the Hartmann layers and side layers of Hunt’s conductive walls and Shercliff’s insulated walls. A general second-order projection method has been developed for the incompressible Navier–Stokes equations with the Lorentz force included. This projection method can accurately balance the pressure term and the Lorentz force for a fully developed core flow. This method can also simplify the pressure boundary conditions for MHD flows.     

An inverse radiative transfer problem of simultaneously estimating profiles of temperature and radiative parameters from boundary intensity and temperature measurements

A parallel-plane space filled with absorbing, emitting, isotropically scattering, gray medium is studied in this paper. The boundary intensity and boundary temperature profiles are calculated for the inverse analysis. For the simultaneous estimation of temperature, absorption and scattering coefficient profiles in the medium, the sum of residuals of boundary intensity and temperature after being weighted by a balance factor is minimized through using a Newton-type iteration algorithm and the least-squares method. To avoid over-updating for the parameters, the relative updating magnitude during the iteration process is constrained not to be >0.5. It is shown that the boundary intensity measurement alone is not enough to estimate simultaneously the temperature (source) and the radiative properties (both absorption and scattering coefficients) when the measurement data contain sensitive random errors. The boundary temperature measurement can serve as a necessary supplementation to the boundary intensity to make this kind of inverse radiative transfer problem resolvable. It was shown that a compensation relationship between absorption and scattering coefficients makes it difficult to fix them accurately. Parabolic profiles for the three parameters are used to validate the estimation method. When the optical thickness approaches 4.0, the results for the radiative properties are not acceptable, although the result for temperature profile is reasonable. This means the method needs further improvements.     

A comparison of differential and integral equations of radiative transfer

The equations of radiative transfer and of statistical equilibrium of a two-level atom are solved by means of differential and integral equations for a one-dimensional medium. The numerical solutions are compared to the analytic solution. It is found that the integral equation for piecewise quadratic source functions gives more accurate results than does the differential equation.     

Variable viscosity effects on the flow of MHD hybrid nanofluid containing dust particles over a needle with Hall current—a Xue model exploration

This study scrutinizes the flow of engine oil-based suspended carbon nanotubes magneto-hydrodynamics (MHD) hybrid nanofluid with dust particles over a thin moving needle following the Xue model. The analysis also incorporates the effects of variable viscosity with Hall current. For heat transfer analysis, the effects of the Cattaneo-Christov theory and heat generation/absorption with thermal slip are integrated into the temperature equation. The Tiwari-Das nanofluid model is used to develop the envisioned mathematical model. Using similarity transformation, the governing equations for the flow are translated into ordinary differential equations. The bvp4c method based on Runge-Kutta is used, along with a shooting approach. Graphs are used to examine and depict the consequences of significant parameters on involved profiles. The results revealed that the temperature of the fluid and boundary layer thickness is diminished as the solid volume fraction is raised. Also, with an enhancement in the variable viscosity parameter, the velocity distribution becomes more pronounced. The results are substantiated by assessing them with an available study.     

多维梯度折射率介质内辐射传递的扩散近似

推导多维梯度折射率介质内稳态辐射传递的扩散近似方程.使用有限元法对扩散近似进行离散和求解,利用两个二维半透明介质的稳态辐射传递问题验证该扩散近似的精度及适用性.算例考虑介质为均匀折射率及梯度折射率两种情况.利用扩散近似分别求解辐射平衡时的边界热流、介质内温度场分布,并与辐射传递方程的求解结果进行对比分析.结果表明:介质折射率变化、散射特性、光学厚度及散射反照率均直接影响扩散近似的精度;在光学厚及强散射条件下,该扩散近似可以作为一种快速算法应用于梯度折射率介质稳态辐射传递的求解.     

Transient coupled heat transfer in a rectangular medium with black surfaces

The main objective of this paper is to extend to two-dimensional (2-D) medium the ray tracing-node analyzing method, which has already been successfully used to solve one-dimensional (1-D) problem of coupled heat transfer in a semitransparent medium. For simplicity, an infinitely long rectangular semitransparent medium with four black opaque surfaces is chosen as our studying object. A control volume method in the implicit scheme is adopted for discretizing the partial transient energy equation. In combination with spectral band model, the radiative heat source term is calculated using the radiative transfer coefficients (RTCs), which are deduced by the ray tracing method. The Partankar's linearization method is used to linearize the radiative source term and the opaque boundary condition, and the linearized equations are solved by the ADI method. Effects of absorption coefficient, refractive index and conductivity on transient cooling process in the 2-D gray rectangular medium are investigated under the condition that the radiation and convection processes cool one side of the rectangular medium while heat the remaining three sides.     

A variational treatment of the mixed boundary condition for the equilibrium diffusion equation

We formulate the asymptotic boundary layer analysis which leads to a mixed (Robbin) boundary condition for the equilibrium diffusion equation of radiative transfer. The information required from the nonlinear boundary layer equations to obtain the boundary condition is obtained by deriving and applying an appropriate variational principle. Based upon an examination of certain limiting cases and the good accuracy previously reported in the literature for similar variational treatments, the accuracy of this variational procedure is expected to be quite good over a wide range of conditions.