RESUMEN
Research on conjugated radiation-conduction (CRC) heat transfer in participating media is of vital scientific and engineering significance due to its extensive applications. Appropriate and practical numerical methods are essential to forecast the temperature distributions during the CRC heat-transfer processes. Here, we established a unified discontinuous Galerkin finite-element (DGFE) framework for solving transient CRC heat-transfer problems in participating media. To overcome the mismatch between the second-order derivative in the energy balance equation (EBE) and the DGFE solution domain, we rewrite the second-order EBE as two first-order equations and then solve both the radiative transfer equation (RTE) and the EBE in the same solution domain, resulting in the unified framework. Comparisons between the DGFE solutions with published data confirm the accuracy of the present framework for transient CRC heat transfer in one- and two-dimensional media. The proposed framework is further extended to CRC heat transfer in two-dimensional anisotropic scattering media. Results indicate that the present DGFE can precisely capture the temperature distribution at high computational efficiency, paving the way for a benchmark numerical tool for CRC heat-transfer problems.
RESUMEN
The lattice Boltzmann method (LBM) has been developed as a powerful solution method in computational fluid dynamics and heat transfer. However, the development of the LBM for solving radiative transfer problems has been far from perfect. This paper proposes a generalized form of the lattice Boltzmann model for the multidimensional radiative transfer equation (RTE) in irregular geometry with a graded index based on body-fitted coordinates. The macroscopic RTE is recovered from Chapman-Enskog analysis, which provides two possible procedures to formulate the Boltzmann equation in graded-index media and irregular geometries. These proposed models have been tested by considering one-and two-dimensional problems of the RTE, and the benchmark solutions reported in the literature were used for comparisons. Afterwards, the LBM is used to analyze the radiation transport in graded-index media for various forms of scattering law, refractive index, boundary reflection, laser and optical properties, and temperatures. The graded-index function and the geometry type have a significant effect on radiative transport in cases in which the refractive index matches or mismatches the boundary. It is also apparent that the developed LBM is an efficient, powerful, robust, and accurate solver for radiative transport in inhomogeneous media with a graded-index function and irregular geometries.
RESUMEN
Based on our previously proposed modified Monte Carlo method, which is efficient to simulate the time-dependent polarized radiative transfer problem in an atmosphere-ocean model with a reflective/refractive interface, we further investigate the square pulse effect on the polarized radiative transfer in an atmosphere-ocean model. A short square pulse, with a duration of nanoseconds, is assumed to be incident at the top of the atmosphere. The polarized signals varying with time and directions are presented for the locations just above and below the atmosphere-water interface and at the bottom of the ocean, and effects of the incidence and disappearance of the external pulse on the Stokes vector components are analyzed. Results in this paper present the general distribution of square-pulse induced polarized signals and they are important for signal analysis in the field of remote sensing using nanosecond pulses.
RESUMEN
Time-dependent polarized radiative transfer in an atmosphere-ocean system exposed to external illumination is numerically investigated. The specular reflection and transmission effects based on the relative refractive index between the atmosphere and water are considered. A modified Monte Carlo (MMC) algorithm combined with time shift and superposition principle, which significantly improves the computational efficiency of the traditional Monte Carlo (TMC) method, is developed to simulate the time-dependent polarized radiative transfer process. The accuracy and computational superiority of the MMC for polarized radiative transfer in the atmosphere-ocean system are validated, and the time-resolved polarized radiative signals are discussed.
RESUMEN
The radiative transfer equation (RTE) has two asymptotic regimes characterized by the optical thickness, namely, optically thin and optically thick regimes. In the optically thin regime, a ballistic or kinetic transport is dominant. In the optically thick regime, energy transport is totally dominated by multiple collisions between photons; that is, the photons propagate by means of diffusion. To obtain convergent solutions to the RTE, conventional numerical schemes have a strong dependence on the number of spatial grids, which leads to a serious computational inefficiency in the regime where the diffusion is predominant. In this work, a discrete unified gas kinetic scheme (DUGKS) is developed to predict radiative heat transfer in participating media. Numerical performances of the DUGKS are compared in detail with conventional methods through three cases including one-dimensional transient radiative heat transfer, two-dimensional steady radiative heat transfer, and three-dimensional multiscale radiative heat transfer. Due to the asymptotic preserving property, the present method with relatively coarse grids gives accurate and reliable numerical solutions for large, small, and in-between values of optical thickness, and, especially in the optically thick regime, the DUGKS demonstrates a pronounced computational efficiency advantage over the conventional numerical models. In addition, the DUGKS has a promising potential in the study of multiscale radiative heat transfer inside the participating medium with a transition from optically thin to optically thick regimes.
RESUMEN
Transient/time-dependent radiative transfer in a two-dimensional scattering medium is numerically solved by the discontinuous finite element method (DFEM). The time-dependent term of the transient vector radiative transfer equation is discretized by the second-order central difference scheme and the space domain is discretized into non-overlapping quadrilateral elements by using the discontinuous finite element approach. The accuracy of the transient DFEM model for the radiative transfer equation considering the polarization effect is verified by comparing the time-resolved Stokes vector component distributions against the steady solutions for a polarized radiative transfer problem in a two-dimensional rectangular enclosure filled with a scattering medium. The transient polarized radiative transfer problems in a scattering medium exposed to an external beam and in an irregular emitting medium are solved. The distributions of the time-resolved Stokes vector components are presented and discussed.
RESUMEN
Transient (time-dependent) polarized radiative transfer in a scattering medium exposed to an external collimated beam illumination is conducted based on the time-dependent polarized radiative transfer theory. The transient term, which persists the nanosecond order time and cannot be ignored for the time-dependent radiative transfer problems induced by a short-pulsed beam, is considered as well as the polarization effect of the radiation. A discontinuous finite element method (DFEM) is developed for the transient vector radiative transfer problem and the derivation of the discrete form of the governing equation is presented. The correctness of the developed DFEM is first verified by comparing the DFEM solutions with the results from the literature. The DFEM is then applied to study the transient polarized radiative transfer induced by a pulsed beam. The time-dependent Stokes vector components are calculated, plotted and analyzed as functions of the axis coordinate and discrete direction. Effects of the diffuse/specular boundary and the incident beam polarization state with respect to the Stokes vector components are further analyzed for cases of different boundary reflection modes and incident beam illuminations.
RESUMEN
Transient radiative transfer induced by a short-pulsed laser in a one-dimensional graded-index medium is investigated by the discontinuous finite element method (DFEM). The boundaries of the medium are Fresnel reflectors, and the incident pulse is considered as the combination of the collimated and the diffuse parts after its first interaction with the medium. The correctness and accuracy of the DFEM solutions for time-resolved reflectance and transmittance are first validated by comparisons with the results obtained by the Monte Carlo method, and the DFEM is then employed to investigate the transient radiative transfer in a graded-index medium with Fresnel boundaries. Effects of the refractive index distributions, the pulse width, the optical thickness, and the scattering phase functions on the transient radiative signals are examined. Several meaningful trends on the time-resolved reflectance and transmittance are observed and analyzed.
