Performance improvement of three-body radiative diodes driven by graphene surface plasmon polaritons.

As an analogue to an electrical diode, a radiative thermal diode allows radiation to transfer more efficiently in one direction than in the opposite direction by operating in a contactless mode. In this study, we demonstrated that within the framework of three-body photon thermal tunneling, the rectification performance of a three-body radiative diode can be greatly improved by bringing graphene into the system. The system is composed of three parallel slabs, with the hot and cold terminals of the diode coated with graphene films and the intermediate body made of vanadium dioxide (VO2). The rectification factor of the proposed radiative thermal diode reaches 300% with a 350 nm separation distance between the hot and cold terminals of the diode. With the help of graphene, the rectification performance of the radiative thermal diode can be improved by over 11 times. By analyzing the spectral heat flux and energy transmission coefficients, it was found that the improved performance is primarily attributed to the surface plasmon polaritons (SPPs) of graphene. They excite the modes of insulating VO2 in the forward-biased scenario by forming strongly coupled modes between graphene and VO2 and thus dramatically enhance the heat flux. However, for the reverse-biased scenario, the VO2 is at its metallic state, and thus, graphene SPPs cannot work by three-body photon thermal tunneling. Furthermore, the improvement was also investigated for different chemical potentials of graphene and geometric parameters of the three-body system. Our findings demonstrate the feasibility of using thermal-photon-based logical circuits, creating radiation-based communication technology and implementing thermal management approaches at the nanoscale.

Lattice Boltzmann modeling of two-phase electrohydrodynamic flows under unipolar charge injection.

In this work, a two-dimensional droplet confined between two parallel electrodes under the combined effects of a nonuniform electric field and unipolar charge injection is numerically investigated using the lattice Boltzmann method (LBM). Under the non-Ohmic regime, the interfacial tension and electric forces at the droplet surface cooperate with the volumetric Coulomb force, leading to complex deformation and motion of the droplet while at the same time inducing a bulk electroconvective flow. After we validate the model by comparing with analytical solutions at the hydrostatic state, we perform a quantitative analysis on the droplet deformation factor D and bulk flow stability criteria T_{c} under different parameters, including the electric capillary number Ca, the electric Rayleigh number T, the permittivity ratio ɛ_{r}, and the mobility ratio K_{r}. It is found that the bulk flow significantly modifies the magnitude of D, which in turn decreases T_{c} of the electroconvective flow. For a droplet repelled by the anode, ɛ_{r}>1, an interesting linear relationship can be observed in the D-Ca curves. However, for a droplet attracted to the anode, ɛ_{r}<1, the system is potentially unstable. After first evolving into a quasisteady state, the droplet successively experiences steady flow, periodic flow, second steady flow, and oscillatory flow with increasing T. Moreover, discontinuities can be observed in the D-T curves due to the transitions of bulk flow.

Efficient lattice Boltzmann method for electrohydrodynamic solid-liquid phase change.

Melting in the presence of electrohydrodynamic (EHD) flow driven by the Coulomb force in dielectric phase change material is numerically studied. A model is developed for the EHD flow in the solid-liquid phase change process. The fully coupled equations including mechanical equations, electrical equations, energy equations, and the continuity equations in the solid-liquid interface are solved using a unified lattice Boltzmann model (LBM). Firstly, the numerical model is validated by several cases in the hydrostatic state, and all LBM results are found to be highly consistent with analytical solutions. Besides, our LBM code is able to reproduce the step changes in the distribution of charge density and electric field due to the discontinuous distribution of physical properties at the interface. Then, a systematical investigation is conducted on various nondimensional parameters, including electric Rayleigh number T, Prandtl number Pr, and Stefan number St. Results are presented for the transient evolutions of temperature, fluid flow, charge density fields, and liquid fraction. Four flow stages in the melting process together with three kinds of flow instabilities are observed. It is found that the electric field has significant influence on the melting, especially at high T and Pr and low St. Over the tested cases, a maximum melting time saving of around 50% is found.

Simple and fast approach to exploit the spectral reflection properties of liquid media.

Bidirectional reflectance distribution functions (BRDFs) are of importance for their wide applications. In this study, we presented a simple and fast approach to measure the spectral BRDF of both solid and liquid samples. Based on this approach, we fabricated a prototype and measured the BRDF value of some liquid samples such as water and NaCl solution at different wavelengths. According to the experimental data, we discussed the trend of the BRDF value of the NaCl solution of different concentrations. Then, the experimental data of the different NaCl solution at 637 nm were used to invert the parameters of a five-parameter model. Additionally, we fitted the parameters as a polynomial.

Multiscale solutions of radiative heat transfer by the discrete unified gas kinetic scheme.

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.

Responses transition in a monolayer Al-Al2O3 nanoparticle-crystal due to oxidation.

Nanoparticle is a promising candidate for large scale fabrication of metamaterial. However, optical responses for metamaterial made of abound metal like Al can be thoroughly changed due to oxidization. Especially for nanoparticle whose aspect ratio is extremely high, oxidation usually occurs. So to understand how the responses shift in a nanoparticle system due to oxidization is essential for large scale application of metamaterial. In this paper, we have concluded and quantified two general principles describing this transition in a monolayer Al-Al2O3 nanoparticle-crystal, which can be used in a thermophotovoltaic system. Square pattern, in which the unit of changing crystal is a square cell made up of Al and Al2O3 particles, is firstly demonstrated. A double oscillators model has been proposed to understand the interference between different absorption modes and their coupling. Using near-field distribution, equivalent inductor-capacitor model and dispersion relationship of surface Plasmon polariton, we have distinguished the resonance modes, concluded the transition principles in a simple case. Then the two principles are applied in a larger cell to verify its university. After detailed demonstration of symmetric square pattern, models and principles are extrapolated to more complex non-symmetric systems. The basic understanding gained here will help the design of robust large-scale metamaterial.

Transient/time-dependent radiative transfer in a two-dimensional scattering medium considering the polarization effect.

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.

Nanoparticle-crystal towards an absorbing meta-coating.

In this paper, a double layer nanoparticle-crystal has been proposed, which shown incident and polarization angle, substrate independences for spectral absorptivity. Such phenomenon originates from the near-field light redistribution and excitation of internal collective oscillating. This kind of nanoparticle-crystal can be made of various types of metal with similar optical responses. A three oscillators mode has been proposed in this paper to understand the shift between global and internal collective oscillating, and verify the physical picture demonstrated. That kind of near-field redistribution result in a prototype of novel meta-coating, and facilitates the large scale application of metamaterial.

Transient polarized radiative transfer analysis in a scattering medium by a discontinuous finite element method.

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.

Multi-focused microlens array optimization and light field imaging study based on Monte Carlo method.

Plenoptic cameras are used for capturing flames in studies of high-temperature phenomena. However, simulations of plenoptic camera models can be used prior to the experiment improve experimental efficiency and reduce cost. In this work, microlens arrays, which are based on the established light field camera model, are optimized into a hexagonal structure with three types of microlenses. With this improved plenoptic camera model, light field imaging of static objects and flame are simulated using the calibrated parameters of the Raytrix camera (R29). The optimized models improve the image resolution, imaging screen utilization, and shooting range of depth of field.

Direct wavefront manipulating for a transverse electric wave microlens.

Due to the polarization nature of the transverse electric electromagnetic wave, manipulating it has been a difficult task and can be even more challenging for integrated on-chip optics. In this Letter, a transverse electric wave manipulating method based on direct wavefront bending and its physical picture have been proposed. Even with only five cells, the microlens can exhibit a focusing pattern and retrieve sub-wavelength spatial features. An analytical mode has been proposed to help understand the physical picture and verify the result. This Letter facilitates the basic understanding for transverse electric wave manipulating and the design of integrated optical elements.

Correction model for microlens array assembly error in light field camera.

In a light field camera, a microlens array (MLA) assembly error can affect the quality of the image. In this study, aiming to ensure corrective imaging using a light field camera, we accurately evaluate and eliminate the assembly error. We used an error image and a standard image to confirm the MLA assembly error, and we developed an assembly error correction model combined with an image quality evaluation index to correct the error. The proposed error correction model can be employed for various assembly errors and different error ranges. Quantitative analyses are performed for these different scenarios. The proposed model can be applied in accurate imaging using a light field camera, four-dimensional optical radiation field information reconstitution, MLA manufacturing and assembly processes, etc.

Lattice Boltzmann model for a steady radiative transfer equation.

A complete lattice Boltzmann model (LBM) is proposed for the steady radiative transfer equation (RTE). The RTE can be regarded as a pure convection equation with a source term. To derive the expressions for the equilibrium distribution function and the relaxation time, an artificial isotropic diffusion term is introduced to form a convection-diffusion equation. When the dimensionless relaxation time has a value of 0.5, the lattice Boltzmann equation (LBE) is exactly applicable to the original steady RTE. We also perform a multiscale analysis based on the Chapman-Enskog expansion to recover the macroscopic RTE from the mesoscopic LBE. The D2Q9 model is used to solve the LBE, and the numerical results obtained by the LBM are comparable to the results obtained by other methods or analytical solutions, which demonstrates that the proposed model is highly accurate and stable in simulating multidimensional radiative transfer. In addition, we find that the convergence rate of the LBM depends on the transport properties of RTE: for diffusion-dominated RTE with a large optical thickness, the LBM shows a second-order convergence rate in space, while for convection-dominated RTE with a small optical thickness, a lower convergence rate is observed.

Vector radiative transfer in a multilayer medium by natural element method.

The vector radiative transfer problem in a vertically multilayer scattering medium with spatial changes in the index of refraction is solved by the natural element method (NEM). The top boundary of the multilayer medium is irradiated by a collimated beam. In our model, the angular space is discretized by the discrete ordinates approach, and the spatial discretization is conducted by the Galerkin weighted residuals approach. In the solution procedure, the collimated component for the Stokes parameters is first solved by NEM, and then it is embedded into the vector radiative transfer equation for the diffuse component as a source term. To keep the consistency of the directions in all the layers, angular interpolation of the Stokes parameters at the interfaces is adopted. The NEM approach for the collimated component is first validated. Then, the classical coupled atmosphere-water system irradiated by different states of collimated beam is examined to verify the numerical performance of the method. Numerical results show that the NEM is accurate, flexible, and effective in solving polarized radiative transfer in a multilayer medium. Finally, polarized radiative transfer in a four-layer system is investigated and analyzed.

Lattice Boltzmann model for Coulomb-driven flows in dielectric liquids.

In this paper, we developed a unified lattice Boltzmann model (LBM) to simulate electroconvection in a dielectric liquid induced by unipolar charge injection. Instead of solving the complex set of coupled Navier-Stokes equations, the charge conservation equation, and the Poisson equation of electric potential, three consistent lattice Boltzmann equations are formulated. Numerical results are presented for both strong and weak injection regimes, and different scenarios for the onset and evolution of instability, bifurcation, and chaos are tracked. All LBM results are found to be highly consistent with the analytical solutions and other numerical work.

Local radial basis function meshless scheme for vector radiative transfer in participating media with randomly oriented axisymmetric particles.

A local radial basis function meshless scheme (LRBFM) is developed to solve polarized radiative transfer in participating media containing randomly oriented axisymmetric particles in which radial basis functions augmented with polynomial basis are employed to construct the trial functions, and the vector radiative-transfer equation based on the discrete-ordinates approach is discretized directly by collocation method. The LRBFM belongs to a class of truly meshless methods that do not need any mesh or any numerical integration scheme. Performances of the LRBFM are verified with analytical solutions and other numerical results reported earlier in the literature via five various test cases. The predicted angular distribution of brightness temperature and Stokes vector by the LRBFM agree very well with the benchmark. It is demonstrated that the LRBFM is accurate to solve vector radiative transfer in participating media with randomly oriented axisymmetric particles.

Optical coherent thermal emission by excitation of magnetic polariton in multilayer nanoshell trimer.

A theoretical demonstration is given of coherent thermal emission via the visible region by exciting magnetic polaritons in isolated metal-dielectric-metal multilayer nanoshells and the collective behavior in a trimer comprising multilayer nanoshells. The dipolar metallic core induces magnetic polaritons in the dielectric shell creating a large enhancement of the emissivity, whose mechanism is different from that of film-coupled metamaterials. The coupling effect of the magnetic polaritons and the electric/magnetic modes of symmetric nanoparticle trimers is discussed to understand the collective behavior in self-assembled nanoparticle clusters with potential solar energy utilizations. The concept of hybridization is employed to understand the collective magnetic polaritons of a multilayer nanoshell trimer. The fundamental understanding gained herein opens up new ways to explore, control, and tailor spectral absorptance, thus facilitating rational design of novel self-assembled nanoclusters for energy harvesting.

Short-pulsed laser transport in two-dimensional scattering media by natural element method.

The natural element method (NEM) is extended to solve transient radiative transfer (TRT) in two-dimensional semitransparent media subjected to a collimated short laser irradiation. The least-squares (LS) weighted residuals approach is employed to spatially discretize the transient radiative heat transfer equation. First, for the case of the refractive index matched boundary, LSNEM solutions to TRT are validated by comparison with results reported in the literature. Effects of the incident angle on time-resolved signals of transmittance and reflectance are investigated. Afterward, the accuracy of this algorithm for the case of the refractive index mismatched boundary is studied. Finally, the LSNEM is extended to study the TRT in a two-dimensional semitransparent medium with refractive index discontinuity irradiated by the short pulse laser. The effects of scattering albedo, optical thickness, scattering phase function, and refractive index on transmittance and reflectance signals are investigated. Several interesting trends on the time-resolved signals are observed and analyzed.

Polarized radiative transfer in an arbitrary multilayer semitransparent medium.

Polarized radiative transfer in a multilayer system is an important problem and has wide applications in various fields. In this work, a Monte Carlo (MC) model is developed to simulate polarized radiative transfer in a semitransparent arbitrary multilayer medium with different refractive indices in each layer. Two kinds of polarization mechanisms are considered: scattering by particles and reflection and refraction at the Fresnel surfaces or interfaces. The MC method has an obvious superiority in that complex mathematical derivations can be avoided in solving the polarization by Fresnel reflection and refraction in an arbitrary multilayer system. We define the vector radiative transfer matrix (VRTM), which describes the polarization characteristics of radiative transfer, and obtain four elements of Stokes vector using the VRTM. The results for the two-layer model by MC method are compared against those for coupled atmosphere-ocean model by the discrete-ordinate method available in the literature, which validates the correctness of the MC multilayer model of polarized radiative transfer. Finally, the results for three-layer, five-layer, and ten-layer models are presented in graphical form. Results show that in the multilayer system, total reflections occurring at the surfaces/interfaces have significant effects on the polarized radiative transfer, which causes abrupt changes or fluctuations like waves in the curves of the Stokes vector.

Approximate solution to vector radiative transfer in gradient-index medium.

Within a gradient-index medium, the radiative rays propagate in curved paths, which makes polarized states change continuously and the solution to the radiative transfer be thus more complex and difficult. In this paper, an arbitrary multilayer model is developed to approximately simulate vector (polarized) radiative transfer in a gradient-index plane-parallel medium. The gradient-index medium is divided into an arbitrary number of sublayers, and each sublayer has a uniform refractive index and two virtual Fresnel's interfaces where only transmission (refraction) is considered. Thus the polarization caused by the curved propagation of lights is approximated by that resulting from refraction on the interfaces. Radiative transfer with consideration of polarization caused by particle scattering and refraction (reflection) on the interfaces (surfaces) in the multilayer model is solved by the MC method. The grid independence of results obtained by the multilayer model for vector radiative transfer in gradient-index medium shows that the convergent solution of Stokes vector will be achieved provided that the sublayer number is large enough. The results for apparent emissivity of gradient-index medium and Stokes vector for two-layer medium are compared well with those in published literatures. Finally, we investigate polarized behaviors of radiative transfer in Rayleigh scattering slabs with linear and sinusoidal gradient indexes and present angular distributions of Stokes vector. Results show that total reflection inside the gradient-index medium resulting from the curved paths traveled by the photons affects greatly the angular distribution characteristics of Stokes vector.