Radiative transfer in cylindrical media

A numerical method for the solution of the radiative transfer equation in a circularly symmetric, cylindrical region is developed. The transfer equation is formulated as a second-order differential equation resulting in a set of tridiagonal difference equations. This form is particularly well suited to line formation and energy balance calculations using the complete linearization method. Several numerical examples are presented.     

Theoretical study of combined radiative conductive and convective heat transfer in semi-transparent porous medium in a spherical enclosure

A new method for the solution of the radiative transfer equation in spherical media based on a modified discrete ordinates method is extended to study radiative, conductive and convective heat transfer in a semi-transparent scattering porous medium. The set of differential equations is solved using the fourth-order Runge-Kutta method. Various results are obtained for the case of combined radiative and conductive heat transfer, as well as for the interaction of those modes with convection. The effects of some radiative properties of the medium on the heat transfer rate are examined.     

Fundamentals of surface thermodynamics

The nonequilibrium behavior for mixtures of fluids in interfaces is discussed. In particular, a thermodynamic field theory is given for media in thin, curved regions (interfaces with finite thickness), which separates two media with different physical properties. The moving interface is considered as semipermeable and a generalized transport equation and specific balance equations are derived. A systematic investigation of constitutive equations is made and in the limit as the thickness of the interface goes to zero it is shown that all relevant interfacial relations can be found.     

Meshless analysis of three-dimensional steady-state heat conduction problems

Steady-state heat conduction problems arisen in connection with various physical and engineering problems where the functions satisfy a given partial differential equation and particular boundary conditions, have attracted much attention and research recently. These problems are independent of time and involve only space coordinates, as in Poisson’s equation or the Laplace equation with Dirichlet, Neuman, or mixed conditions. When the problems are too complex, it is difficult to find an analytical solution, the only choice left is an approximate numerical solution. This paper deals with the numerical solution of three-dimensional steady-state heat conduction problems using the meshless reproducing kernel particle method (RKPM). A variational method is used to obtain the discrete equations. The essential boundary conditions are enforced by the penalty method. The effectiveness of RKPM for three-dimensional steady-state heat conduction problems is investigated by two numerical examples.     

Effect of reflecting modes on combined heat transfer within an anisotropic scattering slab

Under various interface reflecting modes, different transient thermal responses will occur in the media. Combined radiative-conductive heat transfer is investigated within a participating, anisotropic scattering gray planar slab. The two interfaces of the slab are considered to be diffuse and semitransparent. Using the ray tracing method, an anisotropic scattering radiative transfer model for diffuse reflection at boundaries is set up, and with the help of direct radiative transfer coefficients, corresponding radiative transfer coefficients (RTCs) are deduced. RTCs are used to calculate the radiative source term in energy equation. Transient energy equation is solved by the full implicit control-volume method under the external radiative-convective boundary conditions. The influences of two reflecting modes including both specular reflection and diffuse reflection on transient temperature fields and steady heat flux are examined. According to numerical results obtained in this paper, it is found that there exits great difference in thermal behavior between slabs with diffuse interfaces and that with specular interfaces for slabs with big refractive index.     

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.     

The energy equation of a sphere in an unsteady and nonuniform temperature field

The equation for the unsteady heat transfer from a sphere in a viscous/conducting fluid at finite Biot numbers is developed. This process has two characteristic times, one for the diffusion of heat inside the sphere and the other for the diffusion of heat in the external fluid. The solution of the governing equations proves that a general analytical solution may be obtained in the Laplace domain, but not in the time domain. This is due to the complexity of the heat transfer process and the two time scales involved. Asymptotic analytical solutions may be obtained at short and long times. It is observed that the complete form of the energy equation of a small sphere, even at zero Peclet numbers (creeping flow condition) is very complex. Several history terms appear in the energy equation. These emanate from temperature gradients diffused in the two media since the commencement of the heat transfer process.     

Effect of Hall current on the velocity and temperature distributions of Couette flow with variable properties

The unsteady hydromagnetic Couette flow and heat transfer between two parallel porous plates is studied with the Hall effect and temperature dependent properties. The fluid is acted upon by a constant pressure gradient and an external uniform magnetic field as well as uniform suction and injection applied perpendicular to the parallel plates. A numerical solution for the governing non-linear equations of motion and the energy equation are obtained. The effect of the Hall term and the temperature-dependent viscosity and thermal conductivity on both velocity and temperature distributions is examined.     

Thermally stratified flow of hybrid nanofluids with radiative heat transport and slip mechanism: multiple solutions

Research on flow and heat transfer of hybrid nanofluids has gained great significance due to their efficient heat transfer capabilities.In fact,hybrid nanofluids are a novel type of fluid designed to enhance heat transfer rate and have a wide range of engineering and industrial applications.Motivated by this evolution,a theoretical analysis is performed to explore the flow and heat transport characteristics of Cu/Al₂O₃ hybrid nanofluids driven by a stretching/shrinking geometry.Further,this work focuses on the physical impacts of thermal stratification as well as thermal radiation during hybrid nanofluid flow in the presence of a velocity slip mechanism.The mathematical modelling incorporates the basic conservation laws and Boussinesq approximations.This formulation gives a system of governing partial differential equations which are later reduced into ordinary differential equations via dimensionless variables.An efficient numerical solver,known as bvp4c in MATLAB,is utilized to acquire multiple(upper and lower)numerical solutions in the case of shrinking flow.The computed results are presented in the form of flow and temperature fields.The most significant findings acquired from the current study suggest that multiple solutions exist only in the case of a shrinking surface until a critical/turning point.Moreover,solutions are unavailable beyond this turning point,indicating flow separation.It is found that the fluid temperature has been impressively enhanced by a higher nanoparticle volume fraction for both solutions.On the other hand,the outcomes disclose that the wall shear stress is reduced with higher magnetic field in the case of the second solution.The simulation outcomes are in excellent agreement with earlier research,with a relative error of less than 1%.     

Boundary conditions for differential approximations

Low order spherical harmonic (P-N) approximations are applied to a radiative transfer Marshak wave problem. A modified Milne boundary condition is developed for the P-2 approximation, similar to one suggested earlier for the P-1 approximation. Comparison with exact Monte Carlo results suggests that this modified P-2 method may be an accurate and generally applicable differential approximation to the equation of transfer. The Monte Carlo results presented should be useful for testing other approximate formulations of radiative transfer and validating time dependent numerical solution methods for the equation of transfer.     

Higher-order implicit strong numerical schemes for stochastic differential equations

Higher-order implicit numerical methods which are suitable for stiff stochastic differential equations are proposed. These are based on a stochastic Taylor expansion and converge strongly to the corresponding solution of the stochastic differential equation as the time step size converges to zero. The regions of absolute stability of these implicit and related explicit methods are also examined.     

Diffusion front capturing schemes for a class of Fokker–Planck equations: Application to the relativistic heat equation

In this research work we introduce and analyze an explicit conservative finite difference scheme to approximate the solution of initial-boundary value problems for a class of limited diffusion Fokker–Planck equations under homogeneous Neumann boundary conditions. We show stability and positivity preserving property under a Courant–Friedrichs–Lewy parabolic time step restriction. We focus on the relativistic heat equation as a model problem of the mentioned limited diffusion Fokker–Planck equations. We analyze its dynamics and observe the presence of a singular flux and an implicit combination of nonlinear effects that include anisotropic diffusion and hyperbolic transport. We present numerical approximations of the solution of the relativistic heat equation for a set of examples in one and two dimensions including continuous initial data that develops jump discontinuities in finite time. We perform the numerical experiments through a class of explicit high order accurate conservative and stable numerical schemes and a semi-implicit nonlinear Crank–Nicolson type scheme.     

A numerical scheme for the solution of axisymmetric radiative transfer problems in irregular domains filled by media with opaque and transparent diffuse and specular boundaries

This paper presents a new numerical scheme of the discrete ordinates method for the solution of axisymmetric radiative transfer problems in irregular domains filled by media with opaque and transparent diffuse and specular (Fresnel) boundaries and interfaces. New test problems of radiative transfer, which describe radiative transfer in domains with Fresnel interfaces, are proposed in this paper. These problems admit analytic solutions and can be used as benchmark ones. The proposed scheme is applied to the solution of the problems. Numerical results show that the presence of Fresnel interfaces leads to an appreciably larger error in numerical solution. This is connected with the “discontinuity” of the Fresnel reflectivity, which, through numerical diffusion, leads to the distortion of numerical solution. Modification of the scheme allows to reduce the numerical error.     

A plasma resistive diffusion model

A self-consistent study of the slow resistive evolution of an axisymmetric toroidal plasma gives rise to a set of transport equations involving one space variable which require input from the solution of a generalized differential equation obtained from the time-differentiated Grad-Shafranov equation. An iterative scheme is presented for the numerical solution of this generalized differential equation which overcomes the problems of the non-standard boundary conditions. As an illustration this method is used to compute the instantaneous diffusion velocity of a class of model toroidal equilibria. A more detailed study is presented of the time evolution of this model in the cylindrical limit in order to illustrate techniques which can be used in a more complete toroidal simulation.     

Similarity solutions for boundary layer flow and heat transfer of a FENE-P fluid with thermal radiation

This Letter presents a numerical study of the flow and heat transfer of an incompressible FENE-P fluid over a non-isothermal surface. The governing partial differential equations are converted into ordinary differential equations by a similarity transformation. The effects of the thermal radiation are considered in the energy equation, and the variations of dimensionless surface temperature and dimensionless surface temperature gradient, as well as the heat transfer characteristics with various physical parameters are graphed and tabulated. Two cases are studied, namely, (i) the sheet with prescribed surface temperature (PST case) and (ii) the sheet with prescribed heat flux (PHF case). Moreover, the mechanical characteristics of the corresponding flow are also presented.     

变系数瞬态热传导问题边界元格式的特征正交分解降阶方法

提出了一种基于边界元法求解变系数瞬态热传导问题的特征正交分解(POD)降阶方法,重组并推导出变系数瞬态热传导问题适合降阶的边界元离散积分方程,建立了变系数瞬态热传导问题边界元格式的POD降阶模型,并用常数边界条件下建立的瞬态热传导问题的POD降阶模态,对光滑时变边界条件瞬态热传导问题进行降阶分析.首先,对一个变系数瞬态热传导问题,建立其边界域积分方程,并将域积分转换成边界积分;其次,离散并重组积分方程,获得可用于降阶分析的矩阵形式的时间微分方程组;最后,用POD模态矩阵对该时间微分方程组进行降阶处理,建立降阶模型并对其求解.数值算例验证了本文方法的正确性和有效性.研究表明:1)常数边界条件下建立的低阶POD模态矩阵,能够用来准确预测复杂光滑时变边界条件下的温度场结果;2)低阶模型的建立,解决了边界元法中采用时间差分推进技术求解大型时间微分方程组时求解速度慢、算法稳定性差的问题.     

Numerical analysis of thermogravitational turbulent convection in a closed rectangular region with radiation source of energy

The mathematical modeling of the conjugate heat transfer in a closed rectangular region has been carried out under the conditions of the radiation supply of energy. The temperature and stream function fields obtained by the modeling illustrate a substantially unsteady nature of the conjugate heat exchange process under study. An analysis of temperature distributions in typical cross sections of the solution domain has shown a considerable inhomogeneity of the temperature field. It is found that an increase in the Rayleigh number leads to substantial modifications of the temperature and stream function fields. The influence of the distribution of radiation fluxes over the internal interfaces on the temperature fields and the airflow character is shown. The influence of the turbulization on the heat transfer intensity near the interfaces between media has been estimated. Comparisons of the obtained numerical results with experimental data have shown their good agreement.     

The development and investigation of a strongly non-equilibrium model of heat transfer in fluid with allowance for the spatial and temporal non-locality and energy dissipation

The differential equation of heat transfer with allowance for energy dissipation and spatial and temporal nonlocality has been derived by the relaxation of heat flux and temperature gradient in the Fourier law formula for the heat flux at the use of the heat balance equation. An investigation of the numerical solution of the heat-transfer problem at a laminar fluid flow in a plane duct has shown the impossibility of an instantaneous acceptance of the boundary condition of the first kind — the process of its settling at small values of relaxation coefficients takes a finite time interval the duration of which is determined by the thermophysical and relaxation properties of the fluid. At large values of relaxation coefficients, the use of the boundary condition of the first kind is possible only at Fo → ∞. The friction heat consideration leads to the alteration of temperature profiles, which is due to the rise of the intervals of elevated temperatures in the zone of the maximal velocity gradients. With increasing relaxation coefficients, the smoothing of temperature profiles occurs, and at their certain high values, the fluid cooling occurs at a gradientless temperature variation along the transverse spatial variable and, consequently, the temperature proves to be dependent only on time and on longitudinal coordinate.     

Unsteady MHD Non-Darcian Flow over a Vertical Stretching Plate Embedded in a Porous Medium with Thermal Non-Equilibrium Model

An analysis is performed to study the influence of local thermal non-equilibrium (LTNE) on unsteady MHD laminar boundary layer flow of viscous, incompressible fluid over a vertical stretching plate embedded in a sparsely packed porous medium in the presence of heat generation/absorption. The flow in the porous medium is governed by Brinkman-Forchheimer extended Darcy model. A uniform heat source or sink is presented in the solid phase. By applying similarity analysis, the governing partial differential equations are transformed into a set of time dependent non-linear coupled ordinary differential equations and they are solved numerically by Runge-Kutta Fehlberg method along with shooting technique. The obtained results are displayed graphically to illustrate the influence of different physical parameters on the velocity, temperature profile and heat transfer rate for both fluid and solid phases. Moreover, the numerical results obtained in this study are compared with the existing literature in the case of LTE and found that they are in good agreement.     

Two Unconditional Stable Schemes for Simulation of Heat Equation on Manifold Using DEC

To predict the heat diffusion in a given region over time, it is often necessary to find the numerical solution for heat equation. However, the computational domain of classical numerical methods are limited to flat spacetime. With the techniques of discrete differential calculus, we propose two unconditional stable numerical schemes for simulation heat equation on space manifoldand time. The analysis of their stability and error is accomplished by the use of maximum principle.